/* PSPP - a program for statistical analysis.
- Copyright (C) 2005, 2009 Free Software Foundation, Inc.
+ Copyright (C) 2005, 2009, 2011 Free Software Foundation, Inc.
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
int
reg_sweep (gsl_matrix * A, int last_col)
{
- if (A == NULL)
- return GSL_EFAULT;
-
- if (A->size1 != A->size2)
- return GSL_ENOTSQR;
-
int i;
int j;
int k;
gsl_matrix *B;
+ if (A == NULL)
+ return GSL_EFAULT;
+
+ if (A->size1 != A->size2)
+ return GSL_ENOTSQR;
assert (last_col < A->size1);
gsl_matrix_swap_rows (A, A->size1 - 1, last_col);
const struct fmt_spec *wfmt =
cmd->wv ? var_get_print_format (cmd->wv) : &F_8_0;
+ double intercept_ssq;
+ double ssq_effects;
double n_total, mean;
double df_corr = 1.0;
double mse = 0;
mse = gsl_vector_get (ws->ssq, 0) / (n_total - df_corr);
- const double intercept_ssq = pow2 (mean * n_total) / n_total;
+ intercept_ssq = pow2 (mean * n_total) / n_total;
- double ssq_effects = 0.0;
+ ssq_effects = 0.0;
if (cmd->intercept)
{
const double df = 1.0;
double df = categoricals_df (ws->cats, f);
double ssq = gsl_vector_get (ws->ssq, f + 1);
+ double F;
+
ssq_effects += ssq;
if (! cmd->intercept)
ssq += intercept_ssq;
}
- const double F = ssq / df / mse;
+ F = ssq / df / mse;
interaction_to_string (cmd->interactions[f], &str);
tab_text (t, 0, r, TAB_LEFT | TAT_TITLE, ds_cstr (&str));
ds_destroy (&str);
/* Model / Corrected Model */
double df = df_corr;
double ssq = ws->total_ssq - gsl_vector_get (ws->ssq, 0);
+ double F;
+
if ( cmd->intercept )
df --;
else
ssq += intercept_ssq;
- const double F = ssq / df / mse;
+ F = ssq / df / mse;
tab_double (t, 1, heading_rows, 0, ssq, NULL);
tab_double (t, 2, heading_rows, 0, df, wfmt);
tab_double (t, 3, heading_rows, 0, ssq / df, NULL);
covariance_dump_enc_header (const struct covariance *cov, int length)
{
struct tab_table *t = tab_create (cov->dim, length);
+ int n;
int i;
+
tab_title (t, "Covariance Encoding");
tab_box (t,
tab_vline (t, TAL_1, i + 1, 0, tab_nr (t) - 1);
}
- int n = 0;
+ n = 0;
while (i < cov->dim)
{
struct string str;
int idx = i - cov->n_vars;
const struct interaction *iact =
categoricals_get_interaction_by_subscript (cov->categoricals, idx);
+ int df;
ds_init_empty (&str);
interaction_to_string (iact, &str);
- int df = categoricals_df (cov->categoricals, n);
+ df = categoricals_df (cov->categoricals, n);
tab_joint_text (t,
i, 0,