偏相关系数的计算
corpcor包
例子:
M=cbind(y,x1,x2)
xcor=cor(M)#相关系数矩阵
xpcor=cor2pcor(xcor)#偏相关矩阵
ggm包
#函数调用格式:pcor(u,s)
#其中u是一个数值向量,前两个数值表示要计算相关系数的变量下标,其余数值为条件变量的下标
jsbl <- c(1,5) #要计算的相关系数的变量下标
tjbl <- c(2,3,6) #条件(控制)变量的下标,即要排除影响的变量的下标
u <- c(jsbl,tjbl)
s <- cov(pcordata) #变量的协方差
r <- pcor(u,s) #偏相关系数
psych包
##偏相关系数显著性检验
#函数调用格式:pcor.test(r,q,n)
#其中r是由pcor()函数计算得到的偏相关系数,q为要控制的变量数(以数值表示位置),n为样本大小
q <- length(tjbl) #计算要控制的变量数
n <- dim(pcordata)[1] #计算样本量
pcor_test <- pcor.test(r,q,n) #偏相关系数显著性检验结果