我尝试$teststat
从运行了 10 000 次的测试中提取测试统计信息,但遇到了问题:
library(urca)
df.test<-replicate(10000,ur.df(cumsum(rnorm(n=100, mean=0))))
但命令
df.test$teststats
给出输出:NULL
有人可以帮忙吗?
这是一个 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))
teststat
要一次在一个向量中获取所有值,请使用
sapply(df.test, function(z) z@teststat)