Displaying 1 result from an estimated 1 matches for "ssq_smp".
Did you mean:
set_smp
2012 May 11
1
ANOVA question
...(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[index] # extract sample for this label
ssq_smp <- var(smp)*(length(smp)-1) # sum of squares is variance
# times d.f.
# Method 2: from ANOVA residuals
rsd_grp <- rsd[index] # extract residuals
ssq_rsd <- sum(rsd_grp **2) # compute sum of squares
# print m...