×
[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。
a=gl(2, 10, labels=c("a1", "a2"))
b=rep(gl(2, 5, labels=c("b1", "b2")), 2)
value=c(1:5, 2:6, 6:10, 4:8)
aovres <- aov(value~a*b)
summary(aovres)
model.tables(aovres, "means")
# 1要因のプロット
library(gplots)
plotmeans(value~a)
# 2要因のプロット
interaction.plot(x.factor=a, trace.factor=b, response=value)
b=rep(gl(2, 5, labels=c("b1", "b2")), 2)
value=c(1:5, 2:6, 6:10, 4:8)
aovres <- aov(value~a*b)
summary(aovres)
model.tables(aovres, "means")
# 1要因のプロット
library(gplots)
plotmeans(value~a)
# 2要因のプロット
interaction.plot(x.factor=a, trace.factor=b, response=value)
PR
Comment
Trackback
Trackback URL
Comment form
