-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcjmP.R
168 lines (167 loc) · 6.23 KB
/
cjmP.R
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
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
rm(list=ls())
library(readxl)
Cuisse_jambe_metrics_CMT1B <- read_excel("//139.124.150.240/data_project/SUPERSEG/CLIN_SEG/For_Aynaz/Excels/Cuisse_jambe_metrics_CMT1B.xlsx",
col_types = c("text", "text", "text",
"text", "numeric", "numeric", "numeric"))
View(Cuisse_jambe_metrics_CMT1B)
cjmP<-Cuisse_jambe_metrics_CMT1B
head(cjmP)
attach(cjmP)
###########################################################################################################################
MT2<-subset(cjmP,Modality=="T2")
MT2
length(MT2$Modality)
MDP<-subset(cjmP,Modality=="DP")
MDP
length(MDP$Modality)
MFF<-subset(cjmP,Modality=="FF")
MFF
length(MFF$Modality)
MMTR<-subset(cjmP,Modality=="MTR")
MMTR
length(MMTR$Modality)
###########################################################################################################################
names(cjmP)
class(ID)
str(cjmP)
ID<-as.factor(ID)
class(ID)
Part<-as.factor(Part)
class(Part)
Labels<-as.factor(Labels)
class(Labels)
Modality<-as.factor(Modality)
class(Modality)
str(cjmP)
levels(ID)
levels(Part)
levels(Labels)
table(Labels)
levels(Modality)
table(Modality)
summary(Labels)
#########################################################################################################################
#test normality
qqplot(mean[Modality=="DP"],Labels[Modality=="DP"])
qqplot(mean[Modality=="FF"],Labels[Modality=="FF"])
qqplot(mean[Modality=="T2"],Labels[Modality=="T2"])
qqplot(mean[Modality=="MTR"],Labels[Modality=="MTR"])
#########################################################################################################################
boxplot(mean[Modality=="DP"]~Labels[Modality=="DP"],las=2)
anova3<-aov(mean[Modality=="DP"]~Labels[Modality=="DP"])
summary(anova3)
TukeyHSD(anova3)
model3<-lm(mean[Modality=="DP"]~Labels[Modality=="DP"])
summary(model3)
pairwise.t.test(mean[Modality=="DP"],Labels[Modality=="DP"])
plot(TukeyHSD(anova3))
plot(anova3,1)
library(car)
leveneTest(mean[Modality=="DP"],Labels[Modality=="DP"])
#res3<-oneway.test(mean[Modality=="DP"],Labels[Modality=="DP"])
#res3
plot(anova3,2)
aov3_residuals<-residuals(object = anova3)
shapiro.test(x=aov3_residuals)
kruskal.test(mean[Modality=="DP"]~Labels[Modality=="DP"])
pairwise.wilcox.test(mean[Modality=="DP"],Labels[Modality=="DP"],p.adjust.method = "BH")
##########################################################################################################################
boxplot(mean[Modality=="FF"]~Labels[Modality=="FF"],las=2)
anova4<-aov(mean[Modality=="FF"]~Labels[Modality=="FF"])
summary(anova4)
TukeyHSD(anova4)
plot(TukeyHSD(anova4))
model4<-lm(mean[Modality=="FF"]~Labels[Modality=="FF"])
summary(model4)
pairwise.t.test(mean[Modality=="FF"],Labels[Modality=="FF"])
library(car)
plot(anova4,1)
leveneTest(mean[Modality=="FF"],Labels[Modality=="FF"])
res4<-oneway.test(mean[Modality=="FF"]~Labels[Modality=="FF"],data = cjmP)
res4
pairwise.t.test(mean[Modality=="FF"],Labels[Modality=="FF"],p.adjust.method = "BH",pool.sd = FALSE)
plot(anova4,2)
aov4_residuals<-residuals(object = anova4)
shapiro.test(x=aov4_residuals)
kruskalwallisFF<-kruskal.test(mean[Modality=="FF"]~Labels[Modality=="FF"])
kruskalwallisFF
pairwise.wilcox.test(mean[Modality=="FF"],Labels[Modality=="FF"],p.adjust.method = "BH")
##########################################################################################################################
boxplot(mean[Modality=="MTR"]~Labels[Modality=="MTR"],las=2)
anova5<-aov(mean[Modality=="MTR"]~Labels[Modality=="MTR"])
summary(anova5)
pairwise.t.test(mean[Modality=="MTR"],Labels[Modality=="MTR"])
TukeyHSD(anova5)
plot(TukeyHSD(anova5))
model5<-lm(mean[Modality=="MTR"]~Labels[Modality=="MTR"])
summary(model5)
library(car)
plot(anova5,1)
leveneTest(mean[Modality=="MTR"],Labels[Modality=="MTR"])
res5<-oneway.test(mean[Modality=="MTR"]~Labels[Modality=="MTR"])
res5
pairwise.t.test(mean[Modality=="MTR"],Labels[Modality=="MTR"],p.adjust.method = "BH",pool.sd = FALSE)
plot(anova5,2)
aov5_residuals<-residuals(object = anova5)
shapiro.test(x=aov5_residuals)
kruskalwallisMTR<-kruskal.test(mean[Modality=="MTR"]~Labels[Modality=="MTR"])
kruskalwallisMTR
pairwise.wilcox.test(mean[Modality=="MTR"],Labels[Modality=="MTR"],p.adjust.method = "BH")
##########################################################################################################################
boxplot(mean[Modality=="T2"]~Labels[Modality=="T2"],las=2)
anova6<-aov(mean[Modality=="T2"]~Labels[Modality=="T2"])
summary(anova6)
pairwise.t.test(mean[Modality=="T2"],Labels[Modality=="T2"])
TukeyHSD(anova6)
plot(TukeyHSD(anova6))
model6<-lm(mean[Modality=="T2"]~Labels[Modality=="T2"])
summary(model6)
library(car)
plot(anova6,1)
leveneTest(mean[Modality=="T2"],Labels[Modality=="T2"])
#res6<-oneway.test(mean[Modality=="T2"]~Labels[Modality=="T2"])
#res6
plot(anova6,2)
aov6_residuals<-residuals(object = anova6)
shapiro.test(x=aov6_residuals)
kruskalwallisT2<-kruskal.test(mean[Modality=="T2"]~Labels[Modality=="T2"])
kruskalwallisT2
pairwise.wilcox.test(mean[Modality=="T2"],Labels[Modality=="T2"],p.adjust.method = "BH")
#########################################################################################################################
plot(mean~Modality)
anova1<-aov(mean~Modality)
summary(anova1)
model1<-lm(mean~Modality-1)
summary(model1)
pairwise.t.test(mean,Modality)
help(levenetest)
TukeyHSD(anova1)
plot(TukeyHSD(anova1))
plot(anova1,1)
library(car)
leveneTest(mean~Modality)
res1<-oneway.test(mean~Modality)
res1
plot(anova1,2)
aov1_residuals<-residuals(object = anova1)
shapiro.test(x=aov1_residuals)
kruskal.test(mean~Modality)
#########################################################################################################################
plot(mean~Labels)
anova2<-aov(mean~Labels)
summary(anova1)
model2<-lm(mean~Labels-1)
summary(model2)
pairwise.t.test(mean,Labels)
TukeyHSD(anova2)
plot(TukeyHSD(anova2))
plot(anova2,1)
install.packages("car")
library(car)
leveneTest(mean~Labels)
#res2<-oneway.test(mean~Labels)
#res2
plot(anova2,2)
aov2_residuals<-residuals(object = anova2)
shapiro.test(x=aov2_residuals)
#kruskal.test(mean~Labels)