X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=src%2Flanguage%2Fstats%2Fregression.q;h=4ef12c98479ec2aac3b0711aaedef9a7410899c5;hb=fc01177a9fccd839df87c5784abec33b69b4019c;hp=b879cae90871658313ac47b76b270474ec6119a3;hpb=8b42eb7dfd30afe72c307f91c298abe811befc60;p=pspp diff --git a/src/language/stats/regression.q b/src/language/stats/regression.q index b879cae908..4ef12c9847 100644 --- a/src/language/stats/regression.q +++ b/src/language/stats/regression.q @@ -21,26 +21,27 @@ #include #include #include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include -#include + +#include "data/case.h" +#include "data/casegrouper.h" +#include "data/casereader.h" +#include "data/dataset.h" +#include "data/dictionary.h" +#include "data/missing-values.h" +#include "data/transformations.h" +#include "data/value-labels.h" +#include "data/variable.h" +#include "language/command.h" +#include "language/data-io/file-handle.h" +#include "language/dictionary/split-file.h" +#include "language/lexer/lexer.h" +#include "libpspp/compiler.h" +#include "libpspp/message.h" +#include "libpspp/taint.h" +#include "math/covariance.h" +#include "math/linreg.h" +#include "math/moments.h" +#include "output/tab.h" #include "gl/intprops.h" #include "gl/xalloc.h" @@ -791,6 +792,7 @@ identify_indep_vars (const struct variable **indep_vars, } return n_indep_vars; } + static double fill_covariance (gsl_matrix *cov, struct covariance *all_cov, const struct variable **vars, @@ -803,12 +805,15 @@ fill_covariance (gsl_matrix *cov, struct covariance *all_cov, size_t dep_subscript; size_t *rows; const gsl_matrix *ssizes; - gsl_matrix *cm; const gsl_matrix *mean_matrix; const gsl_matrix *ssize_matrix; double result = 0.0; - cm = covariance_calculate_unnormalized (all_cov); + gsl_matrix *cm = covariance_calculate_unnormalized (all_cov); + + if ( cm == NULL) + return 0; + rows = xnmalloc (cov->size1 - 1, sizeof (*rows)); for (i = 0; i < n_all_vars; i++)