我正在运行2个连续VAR之间的部分相关性,我想使用从中提取的估计来可视化关联:
library(ppcor)
pcor.test(Y~X, Z, data, method="spearman")
estimate p.value statistic n gp Method
0.3997551 0.0007318872 3.543037 69 1 spearman
使用偏回归图显示
Y_resid<-resid(lm(Y~Z,data))
X_resid<-resid(lm(X~Z,data))
library(ggplot2)
m<-ggplot(data, aes(x=X_resid, y=Y_resid)) +
geom_point() +
labs(x="X | Z", y = "Y | Z")+
scale_size_manual(values=c(15))+
theme_classic()
m + geom_smooth(method=lm)
geom_smooth(method=lm, formula= y~x)
然而,选项显然不符合pcor函数的估计值,该函数基于Spearman。
有没有办法做XY(渣油)图并使用pcor.test的估算拟合一条线?
提前致谢