Question

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

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

应用Augmented 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]

R：从摘要（）

1 个答案: