不能由$ teststat從列表中提取檢驗統計量中的R

Question

我嘗試$teststat從我已經運行10萬次測試，以提取檢驗統計量，但比我遇到的問題：

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

但命令

df.test$teststats 

給輸出：NULL

任何人可以幫助？

Answer 1

這是一個s4對象，testsats位於插槽中。我想你想是這樣的：

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

#    tau1 
#statistic 0.8491917 

把所有的結果，那麼你prabably要像

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

Answer 2

爲了得到一個向量一次性的值，使用

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