-
Notifications
You must be signed in to change notification settings - Fork 0
/
stats.m
111 lines (84 loc) · 3.06 KB
/
stats.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
% Azadeh example for using multcompare for stats
load('AvgStartCount.mat')
% labeling the groups
Label = {'ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl',...
'Sham','Sham','Sham','Sham','Sham','Sham','Sham','Sham','Sham','Sham',...
'2W','2W','2W','2W','2W','2W','2W','2W','2W','2W','2W',...
'1M','1M','1M','1M','1M','1M','1M','1M','1M',...
'EE','EE','EE','EE','EE','EE','EE','EE','EE','EE',...
};
Label2 = {'ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl',...
'ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl',...
'Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk',...
'Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk',...
'Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk',...
};
Label3 = {'ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl','ctrl',...
'sham','sham','sham','sham','sham','sham','sham','sham','sham','sham',...
'Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk',...
'Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk',...
'Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk','Stk',...
};
%% all groups: label
% Left
% values for each group
Left = FullStartCountL;
Left(find(isnan(Left))) = [];
% one-way anova for gruops with different sizes
[p,a,s] = anova1(Left,Label);
% multiple comparison between groups
[c,m,h,nms] = multcompare(s);
% table of comparisons
[nms(c(:,1)), nms(c(:,2)), num2cell(c(:,3:6))]
% Right
% values for each group
Right = FullStartCountR;
Right(find(isnan(Right))) = [];
% one-way anova for gruops with different sizes
[p,a,s] = anova1(Right,Label);
% multiple comparison between groups
[c,m,h,nms] = multcompare(s);
% table of comparisons
[nms(c(:,1)), nms(c(:,2)), num2cell(c(:,3:6))]
%% sham and control versus all stroke: label 2
% Left
% values for each group
Left = FullStartCountL;
Left(find(isnan(Left))) = [];
% one-way anova for gruops with different sizes
[p,a,s] = anova1(Left,Label2);
% multiple comparison between groups
[c,m,h,nms] = multcompare(s);
% table of comparisons
[nms(c(:,1)), nms(c(:,2)), num2cell(c(:,3:6))]
% Right
% values for each group
Right = FullStartCountR;
Right(find(isnan(Right))) = [];
% one-way anova for gruops with different sizes
[p,a,s] = anova1(Right,Label2);
% multiple comparison between groups
[c,m,h,nms] = multcompare(s);
% table of comparisons
[nms(c(:,1)), nms(c(:,2)), num2cell(c(:,3:6))]
%% control versus sham and all stroke: label 3
% Left
% values for each group
Left = FullStartCountL;
Left(find(isnan(Left))) = [];
% one-way anova for gruops with different sizes
[p,a,s] = anova1(Left,Label3);
% multiple comparison between groups
[c,m,h,nms] = multcompare(s);
% table of comparisons
[nms(c(:,1)), nms(c(:,2)), num2cell(c(:,3:6))]
% Right
% values for each group
Right = FullStartCountR;
Right(find(isnan(Right))) = [];
% one-way anova for gruops with different sizes
[p,a,s] = anova1(Right,Label2);
% multiple comparison between groups
[c,m,h,nms] = multcompare(s);
% table of comparisons
[nms(c(:,1)), nms(c(:,2)), num2cell(c(:,3:6))]