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

Question

我尝试通过 $teststat 从我已运行 10 000 次的测试中提取测试统计信息，但我遇到了问题:

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))