search for: ssq_rsd

...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), ssq_smp-ssq_rsd), "") # ...and it works like...