如何从 cor.test() 中提取 p.value 和估计值？

我在 for 循环中对数据集执行 cor.test，但我不知道如何从测试中提取估计和 tau 等信息。

在数据集中执行 for 循环之前，cor.test() 函数返回如下：

cor.test(armpit$Corynebacterium.1, armpit$Staphylococcus.1, alterantive="two-sided", method="kendall", exact=FALSE, continuity=TRUE)

返回结果

这是我执行 for 循环的代码。现在我想从我的测试中提取估计值和 tau 值。

for (i in 1:8) {
  for (j in 1:8) {
    if (j != i)
    cor.test( as.numeric(unlist(armpit[i])),
        as.numeric(unlist(armpit[j])), alterantive="two-sided",
        method="kendall", exact=FALSE, continuity=TRUE)
   }
}

我已经检查了类似的问题

类似的问题

然后我将代码更改为：

estimates = numeric(50)
pvalues = numeric(50)
for (i in 1:8) {
  for (j in 1:8) {
    if (j != i)
    cor.test( as.numeric(unlist(armpit[i])),
        as.numeric(unlist(armpit[j])), alterantive="two-sided",
         method="kendall", exact=FALSE, continuity=TRUE)
    estimates[i] = cor.test$estimate
    pvalues[i]= cor.test$p-value
   }
   }

但它返回：

Error in cor.test$estimate : object of type 'closure' is not subsettable

有人可以为我提供一些关于如何在 for 循环中从 cor.test() 函数中提取估计值和 tau 值的帮助吗？提前致谢。

cor.test返回一个列表。您可以创建一个对象来捕获此列表：

cor_test <- cor.test( as.numeric(unlist(armpit[i])), as.numeric(unlist(armpit[j])), alterantive="two-sided", method="kendall", exact=FALSE, continuity=TRUE)

然后使用cor_test之后与$访问列表中的每个元素：

estimates[i] = cor_test$estimate  
pvalues[i]= cor_test$p.value  # note the ., not the -

最初的错误非常神秘，所以你对此感到困惑是可以理解的。你写了cor.test$estimate，它要求 R 访问estimate的组成部分cor.test function，不是测试的结果。

estimates = numeric(50)
pvalues = numeric(50)
for (i in 1:8) {
  for (j in 1:8) {
    if (j != i)
    cor_test <- 
      cor.test( as.numeric(unlist(armpit[i])),
        as.numeric(unlist(armpit[j])), alterantive="two-sided",
         method="kendall", exact=FALSE, continuity=TRUE)
    estimates[i] = cor_test$estimate
    pvalues[i]= cor_test$p.value
   }
 }