keep rank of X when fitting model

This is required to calculate the same p-values as those returned by R. The only change is to use "gsl_multifit_wlinear_svd" instead of "gsl_multifit_wlinear".