Displaying 1 result from an estimated 1 matches for "mod_aov".
2012 May 11
1
ANOVA question
...t;- data.frame(Label=rep(LABELS, each=N))
df$Value <- NA
# fill in random data for each factor level
for (i in 1:length(MEAN)) {
df$Value[(1+N*(i-1)):(N*i)] <- rnorm(N, MEAN[i], sqrt(VAR[i]))
}
par(mfrow=c(2,2))
plot(df) # Box plot of the data
plot(df$Value) # scatter plot
mod_aov <- aov(Value ~ Label, data=df)
print(summary(mod_aov))
print(mod_aov$coefficients)
rsd <- mod_aov$residuals
plot(rsd)
# find and print mean() and var() for each level
for (l in levels(df$Label)) {
index <- df$Label == l
# Method 1: directly from data
smp <- df$Value[ind...