123456789101112131415161718192021222324252627282930313233 |
- function [Pvalue, result, gnames] = Stat_ANOVA2(Y,GROUP,NAME)
-
- % Example 1
- % a=[1;2;3;4;5]
- % b=[2;2;3;4;4]
- % c=[1;3;3;2;2]
- % d=[2;3;4;4;1]
- % A=[1;2;3;4;5]
- % B=[2;2;3;4;4]
- % C=[1;3;3;2;2]
- % D=[2;3;4;4;1]
- % Y = [a; b; c; d;...
- % A; B; C; D];
- % Species = [GetGroupingVar(a,1); GetGroupingVar(b,1); GetGroupingVar(c,1); GetGroupingVar(d,1); ...
- % GetGroupingVar(A,2); GetGroupingVar(B,2); GetGroupingVar(C,2); GetGroupingVar(D,2)];
- % Block = [GetGroupingVar(a,1); GetGroupingVar(b,2); GetGroupingVar(c,3); GetGroupingVar(d,4); ...
- % GetGroupingVar(A,1); GetGroupingVar(B,2); GetGroupingVar(C,3); GetGroupingVar(D,4)];
- % NAME = {'Species', 'Block'};
- % [Pvalue, result, gnames] = Stat_ANOVA2(Y,[Species Block],NAME)
- % Example 2
- % Y = [dHP_Odd; dHP_Even; iHP_Odd; iHP_Even];
- % Region = [GetGroupingVar(dHP_Odd,1); GetGroupingVar(dHP_Even,1); ...
- % GetGroupingVar(iHP_Odd,2); GetGroupingVar(iHP_Even,2)];
- % Trajectory = [GetGroupingVar(dHP_Odd,1); GetGroupingVar(dHP_Even,2);...
- % GetGroupingVar(iHP_Odd,1); GetGroupingVar(iHP_Even,2)];
- % NAME = {'Region', 'Trajectory'};
- % [Pvalue, result, gnames] = Stat_ANOVA2(Y,[Region Trajectory],NAME)
- [Pvalue.Main_TwoWay,result,stats] = anovan(Y,GROUP,'model','interaction','varnames',NAME,'display','off');
- [Pvalue.MultipleComparison, ~, ~, gnames] = multcompare(stats,'Dimension',[1,2],'display','off');
- end
|