问题描述
我想知道我的相关矩阵的 p 值四舍五入到三位小数和使用 Benjamini Hochberg 校正进行调整。
为了这个目标,我尝试使用
psych::corr.p(as.matrix(cormatrix),n = 999,adjust = "BH")
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)
有效