文章目錄

• 偏相關(guān)系數(shù)的計(jì)算
• • corpcor包
  • ggm包
  • psych包

偏相關(guān)系數(shù)的計(jì)算

corpcor包

例子：

M=cbind(y,x1,x2) xcor=cor(M)#相關(guān)系數(shù)矩陣 xpcor=cor2pcor(xcor)#偏相關(guān)矩陣

ggm包

#函數(shù)調(diào)用格式：pcor(u,s) #其中u是一個(gè)數(shù)值向量，前兩個(gè)數(shù)值表示要計(jì)算相關(guān)系數(shù)的變量下標(biāo)，其余數(shù)值為條件變量的下標(biāo)jsbl <- c(1,5) #要計(jì)算的相關(guān)系數(shù)的變量下標(biāo) tjbl <- c(2,3,6) #條件(控制)變量的下標(biāo)，即要排除影響的變量的下標(biāo) u <- c(jsbl,tjbl) s <- cov(pcordata) #變量的協(xié)方差 r <- pcor(u,s) #偏相關(guān)系數(shù)

psych包

##偏相關(guān)系數(shù)顯著性檢驗(yàn) #函數(shù)調(diào)用格式：pcor.test(r,q,n) #其中r是由pcor()函數(shù)計(jì)算得到的偏相關(guān)系數(shù)，q為要控制的變量數(shù)(以數(shù)值表示位置)，n為樣本大小 q <- length(tjbl) #計(jì)算要控制的變量數(shù) n <- dim(pcordata)[1] #計(jì)算樣本量 pcor_test <- pcor.test(r,q,n) #偏相關(guān)系數(shù)顯著性檢驗(yàn)結(jié)果

總結(jié)

以上是生活随笔為你收集整理的R语言中用于计算相关系数和偏相关系数的包的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。