Displaying 1 result from an estimated 1 matches for "rsd_grp".
2012 May 11
1
ANOVA question
...{
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 mean, variance, and difference between SSQs from the two
# methods.
write(sprintf("%s: mean=%5.1f var=%5.1f (%.2g)", l,
mean(smp), var(smp)...