計算R中各組之間平方和的公式

Question

誰能告訴我如何編碼R之間的SS以便手工計算，它是∑ ni（meanXi-均值）2

謝謝，lp

Answer 1

如果您有一個向量x和x_mean的平均值， x_mean可以像下面這樣手動計算SS誤差：

> x=c(1,2,3,4,5)
> x_mean = mean(x)
> x-x_mean
[1] -2 -1  0  1  2
> (x-x_mean)^2
[1] 4 1 0 1 4
> sum((x-x_mean)^2)
[1] 10

Answer 2

不確定這是您想要的，但是

#  create sample dataset: 5 groups, 10 values per group
set.seed(1)
df   <- data.frame(group=rep(LETTERS[1:5],each=10),value=rnorm(50))
# calculate between-group sum of squares (SSB)
sum((aggregate(value~group,df,mean)$value-mean(df$value))^2)
# [1] 0.07938908

這將使用aggregate(...)按組計算平均值，然后求和該平均值與總體平均值（ mean(df$value) ）的平方差之和。