如何绘制 corr.test 输出？

Question

我有两个行数相同但列数不同的数据框。我从 psych 包执行 corr.test 使用：

Correlations <- corr.test(final_x, final_y, use = "pairwise", alpha=0.05, method = "spearman")

现在我在运行 corr.test 之后有了输出（例如相关性）。我可以单独保存 Correlation$p 或 Correlation$r 并制作热图。

但我的问题是我想要 Correlation 本身的对角线图，Correlation$p 在图的右上角，Correlations$r 在图的左下角.

我试过了，但没用：

corPlot(Correlations,numbers=TRUE,colors=TRUE,n=51,main=NULL,zlim=c(-1,1),
        show.legend=TRUE, labels=NULL,n.legend=10,keep.par=TRUE,select=NULL, pval=NULL,  
        cuts=c(.001,.01),scale=TRUE,cex,MAR,upper=TRUE,diag=TRUE, symmetric=TRUE,stars=FALSE,
        adjust="holm",xaxis=1, xlas=0,ylas=2,gr=NULL,alpha=.75,min.length=NULL)

如果你能帮我一些忙就好了。

谢谢

Answer 1

这个例子可能有帮助：

mcor <- cor(mtcars) # correlation matrix
mcor
mydata=mtcars
library(corrplot)
corrplot(mcor, type="upper", order="hclust", tl.col="black", tl.srt=45) # print correlation

要打印散点图：

library("ggpubr")
data=iris
str(data)
ggscatter(data, x = "Sepal.Length", y = "Petal.Length", 
          add = "reg.line", conf.int = TRUE, 
          cor.coef = TRUE, cor.method = "pearson",
          xlab = "xlab", ylab = "ylab")

用于绘图；

library("ggpubr")
library("Hmisc")
mydata=mtcars
mydata.rcorr = rcorr(as.matrix(mydata))
mydata.coeff = mydata.rcorr$r
mydata.p = mydata.rcorr$P
plot(mydata.coeff,mydata.p)

如何绘制 corr.test 输出？

How to plot corr.test output?

plot

r

ggplot2

correlation