r - 从具有多个统计测试结果列表的列表中提取统计值

Question

我在 r 中进行了 Ljung-Box 独立性测试，有 36 个滞后，并将结果存储在列表中。

for (lag in c(1:36)){
box.test.list[[lag]] <- (Box.test(btcr, type = "Ljung", lag))
}

我想提取 p 值以及检验统计量(X 平方)并将它们打印出来，如下所示:

X 平方 = 100，p 值 = 0.0001

我也想单独提取 p 值，但不仅仅是吐出数字，我想要类似的东西:

[1] p 值 = 0.001

[2] p 值 = 0.0001

等等。这可以吗？

Answer 1

附测试数据

set.seed(7)
btcr <- rnorm(100)

您可以使用以下命令执行所有测试

box.test.list <- lapply(1:36, function(i) Box.test(btcr, type = "Ljung", i))

然后将所有结果放入 data.frame 中

results <- data.frame(
  lag = 1:36, 
  xsquared = sapply(box.test.list, "[[", "statistic"), 
  pvalue = sapply(box.test.list, "[[", "p.value")
)

然后你可以对结果做你喜欢的事情

head(results)
#   lag  xsquared     pvalue
# 1   1  3.659102 0.05576369
# 2   2  7.868083 0.01956444
# 3   3  8.822760 0.03174261
# 4   4  9.654935 0.04665920
# 5   5 11.190969 0.04772238
# 6   6 12.607454 0.04971085