単純主効果の検定など

SPF.tsmeとRBF.tsmeを更新しました。
使用方法は:
SPF.tsme(分散分析の結果、データ名、従属変数名)
RBF.tsme(分散分析の結果、データ名、従属変数名)
具体例は以下の通りです。
なお、これらを使用する前に
source(“http://peach.l.chiba-u.ac.jp/course_folder/tsme.txt”)
を必ず実行してください。

source("http://peach.l.chiba-u.ac.jp/course_folder/tsme.txt")

# SPFの例
dat<-read.csv("http://www.matsuka.info/data_folder/dktb3211.txt")
dat.aov<-aov(result~method*duration+Error(s+s:duration),dat)
tsme <- SPF.tsme(dat.aov, dat, "result")
simple main effect test for BETWEEEN subject factor 
           ss df   ms       f         p
1hr       2.5  1  2.5  1.0417 0.3150887
2hr       0.4  1  0.4  0.1667 0.6858104
3hr      12.1  1 12.1  5.0417 0.0317831
4hr      32.4  1 32.4 13.5000 0.0008662
residual 76.8 32  2.4                  

 Tukey HSD test - between subject factor @ duration = 3hr 
      D     E
D FALSE  TRUE
E  TRUE FALSE

 Tukey HSD test - between subject factor @ duration = 4hr 
      D     E
D FALSE  TRUE
E  TRUE FALSE

 simple main effect test for WITHIN subject factor 
           ss df      ms      f         p
D        49.2  3 16.4000 12.990 3.047e-05
E         1.0  3  0.3333  0.264 8.506e-01
residual 30.3 24  1.2625                 

 Tukey HSD test - within subject factor @ method = D 
      1hr   2hr   3hr   4hr
1hr FALSE FALSE  TRUE  TRUE
2hr FALSE FALSE FALSE  TRUE
3hr  TRUE FALSE FALSE FALSE
4hr  TRUE  TRUE FALSE FALSE

# RBFの例
dat<-read.csv("http://www.matsuka.info/data_folder/dktb3218.txt")
dat.aov<-aov(result~method*duration+Error(s+s:duration+s:method),dat)
tsme <- RBF.tsme(dat.aov, dat, "result")
simple main effect test for the 1st within-subject factor 
           ss df   ms      f         p
1hr       2.5  1  2.5  1.563 0.2292743
2hr       0.4  1  0.4  0.250 0.6238816
3hr      12.1  1 12.1  7.562 0.0142341
4hr      32.4  1 32.4 20.250 0.0003635
residual 25.6 16  1.6                 

 Tukey HSD test - 1st within-subject factor @ duration = 3hr 
      D     E
D FALSE  TRUE
E  TRUE FALSE

 Tukey HSD test - 1st within-subject factor @ duration = 4hr 
      D     E
D FALSE  TRUE
E  TRUE FALSE

 simple main effect test for the 2nd within-subject factor 
           ss df      ms      f         p
D        49.2  3 16.4000 12.990 3.047e-05
E         1.0  3  0.3333  0.264 8.506e-01
residual 30.3 24  1.2625                 

 Tukey HSD test - 2nd within-subject factor @ method = D 
      1hr   2hr   3hr   4hr
1hr FALSE FALSE  TRUE  TRUE
2hr FALSE FALSE FALSE  TRUE
3hr  TRUE FALSE FALSE FALSE
4hr  TRUE  TRUE FALSE FALSE