1

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

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

但命令

df.test$teststats

给出输出:NULL

有人可以帮忙吗?

4

2 回答 2

3

这是一个 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))
于 2013-04-21T20:52:24.237 回答
0

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

sapply(df.test, function(z) z@teststat)
于 2013-04-21T20:57:36.950 回答