无法通过$ teststat从R中的列表中提取测试统计信息

Question

我试图通过$teststat从我已经运行了10000次的测试中提取测试统计信息，但是我遇到了问题：

library(urca)
df.test<-replicate(10000,ur.df(cumsum(rnorm(n=100, mean=0))))

但是命令

df.test$teststats

给出输出：NULL

有人可以帮忙吗？

Answer 1

这是一个s4对象，测试数据位于插槽中。我想你想要这样的东西：

df.test<-replicate(10,ur.df(cumsum(rnorm(n=100, mean=0))))
df.test[[1]]@teststat

#               tau1
#statistic 0.8491917

获得所有结果然后你可能想要像

这样的东西

unlist( lapply(df.test, function(x) x@teststat))

Answer 2

要在一个向量中一次获取所有teststat值，请使用

sapply(df.test, function(z) z@teststat)