我想知道我的相关矩阵的p值四舍五入到三位小数并使用Benjamini Hochberg校正进行调整。
为此,我尝试使用
psych::corr.p(as.matrix(cormatrix), n = 999, adjust = "BH")
然而,在该函数中,p值四舍五入为两位小数,
round(sych::corr.p(as.matrix(cormatrix), n = 99, adjust = "BH"), 3)
给我错误
数学函数的非数字参数
我想我可以使用p.adjust(pvector, method = "BH"), 3)
。但是corr.test
只给了我以某种形式调整的p值,所以另一次调整可能会导致高估p值。
我不能简单地使用round(corr.test(), 3)
,它会给我调整后的p值和三个小数,因为在我的相关矩阵中是二元相关性和偏相关性,我想同时调整p。
round(corr.p(cormatrix, n = 999)$p, 3)
工作