adjusted degrees of freedom

Author: MSenden

code/PEA.m

@@ -162,8 +162,8 @@ function set_direction(self,direction)
             results.F_stat = zeros(self.n_total,1);
             results.P_value = ones(self.n_total,1);
 
-            df1 = 2;
-            df2 = self.n_samples-1;
+            df1 = 1;
+            df2 = self.n_samples-2;
             for v=1:self.n_total
                 if std_signal>0
                     b = XX * data(:,v);

code/RRT.m

@@ -308,8 +308,9 @@ function optimize_lambda(self,data,range,varargin)
             results.F_stat = zeros(self.n_total,1);
             results.P_value = ones(self.n_total,1);
 
-            df1 = 2;
-            df2 = self.n_samples-1;
+            [~,p] = size(self.X);
+            df1 = p-1;
+            df2 = self.n_samples-p;
             for v=1:self.n_total
                 if mask(v)
                     b = XX * data(:,v);