2

Rの反復測定AnovaでHolm-Sidak検定を実行するにはどうすればよいですか?

MeanSEP は、私がクラスで統計テストを再現しなければならない論文の結果からのボルタグ電位です。条件間の大きな違いを探しています。ベースライン、ABCD、DCBA、および 300。

これが私が得た限りです:

control1 = data.frame(
  c("m1","m2","m3","m4"),
  c(129.43,152.48,118.79,156.02),
  c(180.85,228.72,257.09,278.37),
  c(198.58,230.5,235.82,304.96),
  c(134.75,159.57,223.4,255.32))
colnames(control1) = c("subject","Control","ABCD","DCBA","300")
control=stack(control1) 

subject = rep(control1$subject,4)        # create the "subject" variable
control[3] = subject                   # add it to the new data frame
rm(subject)                               # clean up your workspace
colnames(control) = c("MeanSEP", "Condition", "subject") 

testc= aov(MeanSEP ~ Condition + Error(subject/Condition), data=control)
summary(testc)
4

0 に答える 0