Fix logLik and summary to allow for unknown number of parameters. Since it can't be read from REMLF90 output.
Removed sparse computing of B-splines incidence matrix for compatibility with older versions of R.