v0.0.63 Now fstats.apply() has three methods: regular, Rle, and Matrix (default). The newer methods are either faster or less memory intensive, although which one of the new methods to use depends on the size of the data set. The old method (regular) is there for scenarios where the projection matrix does not have row sums equal to 0 (for example, when the model matrices do not include an intercept) and for comparison purposes. Several functions and manual pages were updated accordingly.