我已经能够使用以下相关矩阵生成:
attach(iris)
library(corrplot)
library(Hmisc)
library(Formula)
library(survival)
#FOR SETOSA:
m<-levels(Species)
setosaCor=cor(iris[Species==m[1],1:4],method = "pearson")
corrplot(setosaCor,method="number",mar=c(0,0,1,0),tl.col="black")
但是,努力获取该矩阵的 p 值。我需要矩阵形式的 p 值。这是我尝试过但没有任何运气的方法
p.value<-rcorr(as.matrix(iris[c(1,2,3,4)]), type=c("pearson"))
cor_mat(iris, vars = NULL, method = "pearson", alternative = "two.sided",conf.level = 0.95)