R语言中用于计算相关系数和偏相关系数的包

偏相关系数的计算

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) #偏相关系数显著性检验结果