1

我用rcorr24 来关联 54 个项目。

别担心，我知道这是很多相关性，可能需要执行额外的统计来解决这个问题。这主要是为了了解我的数据。

使用rcorr我想将生成的表减少到仅低于 0.05 的 p 值。我对与此相关的项目不感兴趣。该命令会是什么样子？现在这是我正在使用的命令

library(Hmisc)
x=matrix(sample.int(7, size = 100*79, replace = TRUE), nrow = 100, ncol = 79)
rcorr(as.matrix(x[,2:25]),as.matrix(x[,26:79]))

我还运行了所有项目的 PCA 以查看它们如何一起加载，但我也想查看相关性。

4

1 回答 1

1

我相信 Frank Harrell 会对此感到震惊，因为他以直指直觉而闻名，这种多重测试策略将无法成为推理的有效基础，但作为编程练习，这可能很有用：

pres <- rcorr(as.matrix(x[,2:25]),as.matrix(x[,26:79]))$P
is.na(pres) <- pres < 0.05
pres

如果您不希望以矩阵格式显示这些值的行号和列号，则应为：

 which( pres < 0.05, arr.ind=TRUE )

（弗兰克，请原谅他们，因为他们不知道自己在做什么。）

于 2014-08-04T19:03:44.103 回答