如何通過方差分析計算出相關系數,總方差是不是能夠分解成組間方差與組內方差(R)


 

 比如在上圖中,我想求出向下比較與主觀幸福感之間的關系

 

#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  值

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM