correct AIC in scr_predval, clean up manual versions

git-svn-id: https://svn.code.sf.net/p/pspm/svn/trunk@92 23ec9670-3a1e-4a12-9187-a0ba8bce340d
bachlab · Jun 5, 2015 · 16431e1 · 16431e1
1 parent 2435f93
commit 16431e1
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/Manual/PsPM 3 Manual.pdf → Manual/PsPM 3.0.1 Manual.pdf b/Manual/PsPM 3 Manual.pdf → Manual/PsPM 3.0.1 Manual.pdf
diff --git a/scr_predval.m b/scr_predval.m
@@ -7,7 +7,8 @@
 %                   
 %           with X: a 2-column vector for paired tests, a 2-element cell
 %                   array for unpaired tests
-%                AIC: computed from RSS of the predictive model
+%                AIC: computed from RSS of the predictive model, up to a
+%                     constant (i. e. only AIC comparisons are meaningful)
 %                T: t-value of the associated t-test
 %                df: df of the associated t-test
 %__________________________________________________________________________
@@ -64,7 +65,8 @@
 [b, dev, stat] = glmfit(X, Y);
 T = stat.t(2); % first column is intercept
 RSS = sum(stat.resid(:).^2);
-AIC = n * log(RSS/n);
+% number of parameters: number of columns in X, plus intercept
+AIC = n * log(RSS/n) + 2 * (1 + size(X, 2));