寻找与缺失值的部分相关性

【问题标题】：Finding partial correlation with missing value寻找与缺失值的部分相关性
【发布时间】：2017-10-08 19:38:46
【问题描述】：

我想找到两个变量之间的偏相关，使用 pearson 方法修复第三个变量。

这是我遇到的错误：

if (det(cvx)

事实上，我的数据中有一个缺失值，并且 (na.rm=T) 对我不起作用。

【问题讨论】：

你能发布一个可重现的例子吗？
例如，如果我有 x=c(2,NA,4,5) , y=c(3,2,6,6) 和 z=c(12,22,14,11 ) 我如何找到使用 pcor.test 函数固定 z 的 x 和 y 之间的相关性？ @jruf003

标签： r

【解决方案1】：

这是 OP 问题的可重现示例：

library("ppcor")

varX <- seq(from=1, to=10, length=10)
varY <- seq(from=20, to=50, length=10)
varZ <- rnorm(10)

varZ[c(1,5,7)] <- NA

pcor.test(x=varX, y=varY, z=varZ)

if (det(cvx)

这在注释部分的 ppcor 包 (https://www.rdocumentation.org/packages/ppcor/versions/1.1/topics/pcor.test) 的文档中得到了解决：Note: Missing values are not allowed

您可以查看 psych 包的 partial.r() 函数。

library("psych")
mydata <- cbind(varX, varY, varZ)
partial.r(data=mydata, x=c("varX","varY"), y="varZ")

【讨论】：