トップ
新規
一覧
単語検索
最終更新
ヘルプ
ログイン
アールメカブ
R_二元配置の分散分析計算式
をテンプレートにして作成
開始行:
[[Rの備忘録]]
# 手計算で2要因の分散分析を行う．手計算といっても，掛け算...
# 山上他 心理統計法 p.155 より
dat <- c(8,10,11, 16,14,15, 2,3,4, 6,6,3, 11,12,8, 11...
dat155 <- data.frame(labB = c(rep("B1",9), rep("B2",9)),
labA = rep(c(rep("A1",3) ,rep("A2", 3),rep("A3",3)),2...
data = dat)
# 本当はこれで済むのだが
summary(aov(data ~ labA*labB, data = dat155))
# 上述書の記載どおりに計算してみる
fA <- 3
fB <- 2
N <- 3
(CT <- (sum(dat))^2 / (fA * fB * N))
(SST <- sum(dat^2) - CT)
(SSA <- (sum(dat155[dat155$labA == "A1",3])^2 +
sum(dat155[dat155$labA == "A2",3])^2 +
sum(dat155[dat155$labA == "A3",3])^2) / (fB*N) - CT)
(SSB <- (sum(dat155[dat155$labB == "B1",3])^2 +
sum(dat155[dat155$labB == "B2",3])^2) / (fA*N) - CT)
(SSC <- (
sum(dat155[dat155$labA == "A1" & dat155$labB == "B1",3])...
sum(dat155[dat155$labA == "A1" & dat155$labB == "B2",3])...
sum(dat155[dat155$labA == "A2" & dat155$labB == "B1",3])...
sum(dat155[dat155$labA == "A2" & dat155$labB == "B2",3])...
sum(dat155[dat155$labA == "A3" & dat155$labB == "B1",3])...
sum(dat155[dat155$labA == "A3" & dat155$labB == "B2",3])...
/ N - CT)
(SSAB <- SSC - (SSA+SSB))
(SSE <- SST - SSC)
終了行:
[[Rの備忘録]]
# 手計算で2要因の分散分析を行う．手計算といっても，掛け算...
# 山上他 心理統計法 p.155 より
dat <- c(8,10,11, 16,14,15, 2,3,4, 6,6,3, 11,12,8, 11...
dat155 <- data.frame(labB = c(rep("B1",9), rep("B2",9)),
labA = rep(c(rep("A1",3) ,rep("A2", 3),rep("A3",3)),2...
data = dat)
# 本当はこれで済むのだが
summary(aov(data ~ labA*labB, data = dat155))
# 上述書の記載どおりに計算してみる
fA <- 3
fB <- 2
N <- 3
(CT <- (sum(dat))^2 / (fA * fB * N))
(SST <- sum(dat^2) - CT)
(SSA <- (sum(dat155[dat155$labA == "A1",3])^2 +
sum(dat155[dat155$labA == "A2",3])^2 +
sum(dat155[dat155$labA == "A3",3])^2) / (fB*N) - CT)
(SSB <- (sum(dat155[dat155$labB == "B1",3])^2 +
sum(dat155[dat155$labB == "B2",3])^2) / (fA*N) - CT)
(SSC <- (
sum(dat155[dat155$labA == "A1" & dat155$labB == "B1",3])...
sum(dat155[dat155$labA == "A1" & dat155$labB == "B2",3])...
sum(dat155[dat155$labA == "A2" & dat155$labB == "B1",3])...
sum(dat155[dat155$labA == "A2" & dat155$labB == "B2",3])...
sum(dat155[dat155$labA == "A3" & dat155$labB == "B1",3])...
sum(dat155[dat155$labA == "A3" & dat155$labB == "B2",3])...
/ N - CT)
(SSAB <- SSC - (SSA+SSB))
(SSE <- SST - SSC)
ページ名:
