R:从摘要中提取值()

Question

我正在使用urca包对单位根进行一些计算。例如，让我们使用urca包中的以下数据集：

library(urca)
data(Raotbl3)
attach(Raotbl3)

应用增强的Dickey Fuller测试如下所示：

lc.df <- ur.df(y=lc, lags=3, type='trend')
summary(lc.df)

作为回报，我们得到所有值的摘要。但是，我只想提取特定的值。可能的字段如下：

slotNames(summary(lc.df))
#  [1] "classname" "test.name" "testreg"   "teststat"  "cval"      "bpoint"        "signif"    "model"     "type"      "auxstat"  
# [11] "lag"       "H"         "A"         "lambda"    "pval"      "V"             "W"         "P"       

有趣的是，

ur.df(y=lc, lags=3, type='trend')@teststat 
#             tau3  phi2  phi3
# statistic -2.239 3.738 2.597

或

ur.df(y=lc, lags=3, type='trend')@cval
#       1pct  5pct 10pct
# tau3 -4.04 -3.45 -3.15
# phi2  6.50  4.88  4.16
# phi3  8.73  6.49  5.47

工作得很好。但是，从"pval“或"lag”中提取值不起作用(返回"NULL")。

当我查看here时，我发现每个"=NULL“的变量都会出现这个问题。不确定是不是很重要。

setMethod("summary", "ur.df", function(object){
  return(new("sumurca", classname="ur.df", test.name=object@test.name,     
             testreg=object@testreg, teststat=object@teststat, 
             cval=object@cval, bpoint=NULL, signif=NULL, 
             model=object@model, type=NULL, auxstat=NULL, lag=NULL, 
             H=NULL, A=NULL, lambda=NULL, pval=NULL, V=NULL, W=NULL, P=NULL))
})

你知道怎么解决这个问题吗？

Answer 1

像往常一样用@获取插槽：

summary(lc.df)@testreg$coefficients
summary(lc.df)@testreg$coefficients[-c(1,3),4]