Displaying 1 result from an estimated 1 matches for "g3mean".
Did you mean:
g1mean
2007 Nov 22
3
anova planned comparisons/contrasts
...= sum(Group %in% G2)
GrandMean = mean(Var[Group %in% c(G1, G2, G3)])
G1Mean = mean(Var[Group %in% G1])
G2Mean = mean(Var[Group %in% G2])
if(is.null(G3))
MScontr = nG1 * ((G1Mean - GrandMean)^2) +
nG2 * ((G2Mean - GrandMean)^2)
else {
nG3 = sum(Group %in% G3)
G3Mean = mean(Var[Group %in% G3])
MScontr = (nG1 * ((G1Mean - GrandMean)^2) +
nG2 * ((G2Mean - GrandMean)^2) +
nG3 * ((G3Mean - GrandMean)^2))/2
}
An <- anova(lm(Var ~ Group))
MSwithin = An[3]['Residuals',]
DegF = An$Df[length(An$Df)]
Fval = MScontr...