比如在上圖中,我想求出向下比較與主觀幸福感之間的關系
#object: 在元分析中,通過方差分析得到相關系數 #writer: mike #time: 2020,11,20 #向量的拼接 a <- c(12,32,45,23,554) b <- c(65,32,567,87,34) total <- c(a,b) #合並成一組不能用 rbind(), cbind(), 這兩個函數無論如何,會增加維度 total1 <- rbind(a,b) total2 <- cbind(a,b) mean(a) sd(a) #計算平方和 var(a)*4 mean(b) sd(b) var(b)*4 mean(total) sd(total) #計算平方和 var(total)*9 #為什么結果對不上? 我之所以對不上是因為var() 的分子其實是 n-1 #下面是書上的例子,單因素完全隨機實驗設計 a1 <- c(3,6,4,3,5,7,5,2) a2 <- c(4,6,4,2,4,5,3,3) a3 <- c(8,9,8,7,5,6,7,6) a4 <- c(9,8,8,7,12,13,12,11) atatal <- c(a1,a2,a3,a4) #試試看到底是 n 還是 n-1 var(atatal)*32 var(atatal)*31 #計算組內SS平方和 (var(a1)+var(a2)+var(a3)+var(a4))*7
總方差分解成組間方差與組內方差是沒問題的,只不過var() 使用的是 n-1
然后計算出 t 值,進而計算出 r 值