update I_ANOVATest and put in error for unequal sample sizes for sp ANOVA

Author: Eric Fields

perm_spANOVA.m

@@ -38,6 +38,10 @@
 %This code is free and open source software made available under the 3-clause BSD license.
 
 function [F_obs, F_dist, df_effect, df_res, exact_test] = perm_spANOVA(data, cond_subs, dims, n_perm, reduce)
+    
+    if ~all(cond_subs == cond_subs(1))
+        error('Split plot ANOVA is currently only available for designs with equal sample sizes');
+    end
 
     %Eliminate factors not involved in this effect and reduce interactions
     %via subtraction

testing/I_ANOVATest.m

@@ -15,7 +15,8 @@
 %% RB DESIGNS PARAMETRIC ANOVA
 
 %Designs to test
-anova_designs = {5, ...
+anova_designs = {2, ...
+                 5, ...
                  [2, 2], ...
                  [4, 2], ...
                  [3, 4], ...
@@ -67,6 +68,7 @@
     oneway_data = reshape(data, n_electrodes, n_time_pts, [], n_subs);
     rm_data = squeeze(oneway_data(e,t,:,:))';
     [rm, ranovatbl] = matlab_ANOVA(rm_data, wg_design, cond_subs, var_names);
+    %disp(ranovatbl);
 
     %Calculate all effects in model and compare to MATLAB ANOVA
     for i = 1:length(effects)
@@ -86,7 +88,8 @@
         end
         assert(abs(test_results.p(e,t) - ranovatbl{rm_table_row, ['pValue' suff]}) < 1e-9);
 
-        %%% Between subjects effects
+        %%% Group interactions %%%
+        
         if ~isempty(cond_subs)
 
             %FMUT calculations
@@ -106,6 +109,11 @@
 
     end
 
+    %%% Group main effect %%%
+        
+    test_results = calc_param_ANOVA(data, cond_subs, ndims(data), 0.05, 'none', sphericity_corr);
+    assert(abs(test_results.p(e,t) - ranovatbl{'group', ['pValue' suff]}) < 1e-9);
+    
 end