我使用心理软件包中的statsBy函数来计算个体之间和个体内部的相关性，如下所示：
statsBy(df[,c(1,10,12,13,14,17,20,24:25)], "group", na.rm = TRUE)
然后我通过以下方式提取p值：
print(stats$pwg, short=FALSE) #P-values within group cor
这显示为矩阵。但是，如statsBy函数的帮助中所述，这些p值未针对多个比较进行校正。我希望通过bonferonni方法这样做。
我尝试通过使用标准p.adjust命令来执行此操作。但是，这会将我的矩阵转换为数值的类（我认为至少是一个类？）。这会导致2个问题：
我不知道哪个p值属于哪个两个变量。
我不能将p值用作corrplot函数的矩阵（来自corrplot程序包）。
有没有一种方法可以在保留p值矩阵的结构的同时使用p.adjust？最好保留完整的原始标头（尽管本身并不是必需的）。
或者，是否有另一个函数可以在不更改矩阵设置的情况下计算经bonferonni调整的p值？