Fix a bug with incorrect degrees of freedom used in standard error
calculation. ncol(model) is dangerous!
Add imidas_r_fast function, but do not export it yet.