+ model = pspp_linreg_cache_alloc (n_data, n_indep);
+ model->depvar = v_dependent;
+ /*
+ For large data sets, use QR decomposition.
+ */
+ if (n_data > sqrt (n_indep) && n_data > GLM_LARGE_DATA)
+ {
+ model->method = PSPP_LINREG_QR;
+ }
+ coeff_init (model, X);
+ pspp_linreg_with_cov (X, model);