@appelmar Could you check this part? Achim implemented two optional parameters to efp[.formula], could this be ported to the matrix function as well? Or maybe we should merge the two functions so that changes like that are easier to integrate?
Originally posted by @GreatEmerald in #29 (comment)