myscratchpad

grp <- factor(rep(c(51,55,57,59,61), each=3))
val <- c(12.2,18.8,18.2,22.2,20.5,14.6, 20.8,19.5,26.3, 26.4,32.6,31.3, 24.5,21.2,22.4)

## Leveneの等分散性の検定。こっちよりもBartlett検定のほうがいいらしい。
## car, lawstat, s20xに同じ関数がある
library(car)
leveneTest(val, grp)
leveneTest(lm(val~grp)) # lm関数のオブジェクトも受け取れる。便利
levene.test(val, grp) # もうすぐなくなるらしい

library(lawstat)
levene.test(val, grp) # デフォルト
levene.test(val, grp, "median") # デフォルトと同じ
levene.test(val, grp, "mean") # spssと同じ出力

# s20xパッケージにも同じ関数がある

## Bartlettの等分散性の検定。Rの組み込み関数
bartlett.test(val, grp)

## 2要因参加者間も同じ。以下のデータを参考。感謝
http://academic.udayton.edu/gregelvers/psy216/SPSS/2wayanovabs.htm

cls <- factor(rep(c("dst", "lct"), each=10))
gpa <- rep(factor(rep(c("high", "low"), each=5)), 2)
pnt <- c(332,380,371,366,354,259.5,302.5,296,349,309,354.67,353.5,304,365,339,306,339,353,351,333)

library(psych)
describe.by(pnt, list(cls, gpa))

fct <- factor(paste(cls, gpa, sep=""))
library(lawstat)
levene.test(pnt, fct)
levene.test(pnt, fct, "mean")
bartlett.test(pnt, fct)