eta平方是如何计算的？

Question

我想知道这张幻灯片中的 eta 平方是如何计算的。 我认为这是组间变异与总平方和的比率。

我尝试重现他们的数据，然后计算 eta 平方。但我最终得到了不同的结果。

G1<-c(6.3,2.8,7.8,7.9,4.9)
G2<-c(9.9,4.1,3.9,6.3,6.9)
G3<-c(5.1,2.9,3.6,5.7,4.5)
G4<-c(1.0,2.8,4.8,3.9,1.6)

n1<-length(G1);n2<-length(G2);n3<-length(G3); n4<-length(G4)

m1<-sum(G1)/n1;m2<-sum(G2)/n2;m3<-sum(G3)/n3;m4<-sum(G4)/n4
Gmean<-sum(c(G1,G2,G3,G4))/(n1+n2+n3+n4)

SSb<-n1*((m1-Gmean)^2)+n2*((m2-Gmean)^2)+n3*((m3-Gmean)^2)+n4*((m4-Gmean)^2)

S1<-var(G1);S2<-var(G2);S3<-var(G3);S4<-var(G4)
SSw<-(n1-1)*S1+(n2-1)*S2+(n3-1)*S3+(n4-1)*S4

SSt<-SSb+SSw

我的公式可能是错误的。

Answer 1

幻灯片似乎是错误的。我使用两种方法计算 eta2 只是为了确保我没有遗漏某些内容，并且它们会显示与您的幻灯片不同的相似结果。因此，如果我犯了错误，我会用两种方法都犯错误。

G1 <- c(6.3, 2.8, 7.8, 7.9, 4.9)
G2 <- c(9.9, 4.1, 3.9, 6.3, 6.9)
G3 <- c(5.1, 2.9, 3.6, 5.7, 4.5)
G4 <- c(1.0, 2.8, 4.8, 3.9, 1.6)

#method 1 (sum of squares as sum of suared deviations):
n1 <- length(G1)
n2 <- length(G2)
n3 <- length(G3)
n4 <- length(G4)

m1 <- sum(G1) / n1
m2 <- sum(G2) / n2
m3 <- sum(G3) / n3
m4 <- sum(G4) / n4

Gmean <- sum(c(G1, G2, G3, G4)) / (n1 + n2 + n3 + n4)

SSb <- n1 * ((m1 - Gmean)^2) + 
  n2 * ((m2 - Gmean)^2) + 
  n3 * ((m3 - Gmean)^2) + 
  n4 * ((m4 - Gmean)^2)

SSw <- sum((G1 - m1)^2) + sum((G2 - m2)^2) + sum((G3 - m3)^2) + sum((G4 - m4)^2)

SSt <- SSb + SSw

eta2 <- SSb / SSt

#method 2 (sum of squares as differences between sum of squares:

Y <- sum(c(G1, G2, G3, G4)^2)
A <- m1^2*n1 + m2^2*n2 + m3^2*n1 + m4^2*n2
Te <- Gmean^2*sum(n1, n2, n3, n4)

SSW <- Y-A
SSB <- A-Te
SST <- Y-Te

eta2_n <- SSb/SSt