X-Git-Url: https://pintos-os.org/cgi-bin/gitweb.cgi?a=blobdiff_plain;f=lib%2Flinreg%2Fsweep.c;h=61326d78dc63ef70d20149b4c98abec97c659fd1;hb=f5c108becd49d78f4898cab11352291f5689d24e;hp=8805840e2c4f0261af1cc3df92267919a20d626d;hpb=7eee0554f378481faf447e2d2e940f389d6b05ec;p=pspp-builds.git diff --git a/lib/linreg/sweep.c b/lib/linreg/sweep.c index 8805840e..61326d78 100644 --- a/lib/linreg/sweep.c +++ b/lib/linreg/sweep.c @@ -24,7 +24,7 @@ Y = Xb + Z - where Y is an n-by-1 column vector, X is an n-by-p matrix of + where Y is an n-by-1 column vector, X is an n-by-p matrix of independent variables, b is a p-by-1 vector of regression coefficients, and Z is an n-by-1 normally-distributed random vector with independent identically distributed components with mean 0. @@ -58,7 +58,7 @@ -- -- X refers to the design matrix and Y to the vector of dependent - observations. reg_sweep sweeps on the diagonal elements of + observations. reg_sweep sweeps on the diagonal elements of X'X. The matrix A is assumed to be symmetric, so the sweep operation is @@ -95,7 +95,7 @@ reg_sweep (gsl_matrix * A) for (j = i; j < A->size2; j++) { /* - Use only the upper triangle of A. + Use only the upper triangle of A. */ if (j < k) {