fixed predict arg list; undid bad indent
authorJason Stover <jhs@math.gcsu.edu>
Sun, 9 Apr 2006 14:00:04 +0000 (14:00 +0000)
committerJason Stover <jhs@math.gcsu.edu>
Sun, 9 Apr 2006 14:00:04 +0000 (14:00 +0000)
src/math/linreg/linreg.c
src/math/linreg/linreg.h

index b2c3c0d76de9fdf1670fcdf55e2cdc709997b465..0406cad2183778097fc1783d016b1ad0c1c3eaa8 100644 (file)
@@ -1,21 +1,21 @@
 /*
* lib/linreg/linreg.c
- * 
* Copyright (C) 2005 Free Software Foundation, Inc. Written by Jason H. Stover.
- * 
* This program is free software; you can redistribute it and/or modify it under
* the terms of the GNU General Public License as published by the Free
* Software Foundation; either version 2 of the License, or (at your option)
* any later version.
- * 
* This program is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License for
* more details.
- * 
* You should have received a copy of the GNU General Public License along with
* this program; if not, write to the Free Software Foundation, Inc., 51
* Franklin Street, Fifth Floor, Boston, MA 02111-1307, USA.
+  lib/linreg/linreg.c
+  
+  Copyright (C) 2005 Free Software Foundation, Inc. Written by Jason H. Stover.
+  
+  This program is free software; you can redistribute it and/or modify it under
+  the terms of the GNU General Public License as published by the Free
+  Software Foundation; either version 2 of the License, or (at your option)
+  any later version.
+  
+  This program is distributed in the hope that it will be useful, but WITHOUT
+  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License for
+  more details.
+  
+  You should have received a copy of the GNU General Public License along with
+  this program; if not, write to the Free Software Foundation, Inc., 51
+  Franklin Street, Fifth Floor, Boston, MA 02111-1307, USA.
  */
 
 #include <gsl/gsl_fit.h>
@@ -102,9 +102,11 @@ pspp_linreg_cache_alloc (size_t n, size_t p)
   c->indep_means = gsl_vector_alloc (p);
   c->indep_std = gsl_vector_alloc (p);
   c->ssx = gsl_vector_alloc (p);       /* Sums of squares for the
-                                        * independent variables. */
+                                          independent variables. 
+                                       */
   c->ss_indeps = gsl_vector_alloc (p); /* Sums of squares for the
-                                        * model parameters. */
+                                          model parameters. 
+                                       */
   c->cov = gsl_matrix_alloc (p + 1, p + 1);    /* Covariance matrix. */
   c->n_obs = n;
   c->n_indeps = p;
@@ -180,7 +182,8 @@ pspp_linreg (const gsl_vector * Y, const gsl_matrix * X,
   cache->dfm = cache->n_indeps;
   cache->dfe = cache->dft - cache->dfm;
   cache->n_coeffs = X->size2 + 1;      /* Adjust this later to allow for
-                                        * regression through the origin. */
+                                          regression through the origin. 
+                                       */
   if (cache->method == PSPP_LINREG_SWEEP)
     {
       gsl_matrix *sw;
index b84c3393ba26d5ef6bc2a81d27398f3010bcd66c..940476424cbdb15d6f0b48f890c08fd41eba5c63 100644 (file)
@@ -1,22 +1,22 @@
 /*
* lib/linreg/linreg.h
- * 
* Copyright (C) 2005 Free Software Foundation, Inc. Written by Jason H. Stover.
- * 
* This program is free software; you can redistribute it and/or modify it under
* the terms of the GNU General Public License as published by the Free
* Software Foundation; either version 2 of the License, or (at your option)
* any later version.
- * 
* This program is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License for
* more details.
- * 
* You should have received a copy of the GNU General Public License along with
* this program; if not, write to the Free Software Foundation, Inc., 51
* Franklin Street, Fifth Floor, Boston, MA 02111-1307, USA.
- */
+  lib/linreg/linreg.h
+  
+  Copyright (C) 2005 Free Software Foundation, Inc. Written by Jason H. Stover.
+  
+  This program is free software; you can redistribute it and/or modify it under
+  the terms of the GNU General Public License as published by the Free
+  Software Foundation; either version 2 of the License, or (at your option)
+  any later version.
+  
+  This program is distributed in the hope that it will be useful, but WITHOUT
+  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License for
+  more details.
+  
+  You should have received a copy of the GNU General Public License along with
+  this program; if not, write to the Free Software Foundation, Inc., 51
+  Franklin Street, Fifth Floor, Boston, MA 02111-1307, USA.
+*/
 
 #ifndef LINREG_H
 #define LINREG_H
@@ -46,25 +46,25 @@ struct pspp_linreg_opts_struct
 
   int get_depvar_mean_std;
   int *get_indep_mean_std;     /* Array of booleans
-                                * dictating which
-                                * independent variables need
-                                * their means and standard
-                                * deviations computed within
-                                * pspp_linreg. This array
-                                * MUST be of length
-                                * n_indeps. If element i is
-                                * 1, pspp_linreg will
-                                * compute the mean and
-                                * variance of indpendent
-                                * variable i. If element i
-                                * is 0, it will not compute
-                                * the mean and standard
-                                * deviation, and assume the
-                                * values are stored.
-                                * cache->indep_mean[i] is
-                                * the mean and
-                                * cache->indep_std[i] is the
-                                * sample standard deviation. */
+                                  dictating which
+                                  independent variables need
+                                  their means and standard
+                                  deviations computed within
+                                  pspp_linreg. This array
+                                  MUST be of length
+                                  n_indeps. If element i is
+                                  1, pspp_linreg will
+                                  compute the mean and
+                                  variance of indpendent
+                                  variable i. If element i
+                                  is 0, it will not compute
+                                  the mean and standard
+                                  deviation, and assume the
+                                  values are stored.
+                                  cache->indep_mean[i] is
+                                  the mean and
+                                  cache->indep_std[i] is the
+                                  sample standard deviation. */
 };
 typedef struct pspp_linreg_opts_struct pspp_linreg_opts;
 
@@ -103,8 +103,8 @@ struct pspp_linreg_cache_struct
   int n_coeffs;
 
   /*
-   * The variable struct is ignored during estimation. It is here so
-   * the calling procedure can find the variable used in the model.
+    The variable struct is ignored during estimation. It is here so
+    the calling procedure can find the variable used in the model.
    */
   const struct variable *depvar;
 
@@ -130,17 +130,18 @@ struct pspp_linreg_cache_struct
    */
   double ssm;                  /* Sums of squares for the overall model. */
   gsl_vector *ss_indeps;       /* Sums of squares from each
-                                * independent variable. */
+                                  independent variable. */
   double sst;                  /* Sum of squares total. */
   double sse;                  /* Sum of squares error. */
   double mse;                  /* Mean squared error. This is just sse /
-                                * dfe, but since it is the best unbiased
-                                * estimate of the population variance, it
-                                * has its own entry here. */
+                                  dfe, but since it is the best unbiased
+                                  estimate of the population variance, it
+                                  has its own entry here. */
   gsl_vector *ssx;             /* Centered sums of squares for independent
-                                * variables, i.e. \sum (x[i] - mean(x))^2. */
+                                  variables, i.e. \sum (x[i] - mean(x))^2. */
   double ssy;                  /* Centered sums of squares for dependent
-                                * variable. */
+                                  variable. 
+                               */
   /*
      Covariance matrix of the parameter estimates.
    */
@@ -158,7 +159,7 @@ struct pspp_linreg_cache_struct
    */
   gsl_matrix *hat;
 
-  double (*predict) (const struct variable *, const union value *,
+  double (*predict) (const struct variable **, const union value **,
                     const void *, int);
 };