0

我想在使用“AER”包中的 ivreg 或“sem”包中的 tsls 后提取 Pr(>|t|) 列。两者都给出了一个相似的术语列表，似乎并没有提供我正在寻找的内容。

ivregest <- ivreg(mdetect~bednet | treat1+treat2, data=simdata)
> names(ivregest)
 [1] "coefficients"  "residuals"     "fitted.values" "weights"      
 [5] "offset"        "n"             "nobs"          "rank"         
 [9] "df.residual"   "cov.unscaled"  "sigma"         "call"         
 [13] "formula"       "terms"         "levels"        "contrasts"    
 [17] "model"         "y"

tslsest <-  tsls(mdetect~bednet , ~ treat1+treat2, data=simdata)
> names(tslsest)
[1] "n"             "p"             "coefficients"  "V"            
[5] "s"             "residuals"     "response"      "model.matrix" 
[9] "instruments"   "weights"       "response.name" "formula"

p 虽然看起来很有希望，但仅提供了在第二阶段回归中估计的参数数量的计数。然而，如果我对这些对象中的任何一个使用 summary 命令，它将返回 ap 值。

所以我真的很想回答两个问题： 1. 我在哪里可以找到 p 值？2. 我怎样才能找到对象的所有隐藏属性，以便下次查找 F-stat 或任何我知道在哪里查找的东西？names() 似乎还不够。

非常感谢您提供的任何帮助！

4

1 回答 1

3

首先，一些模型取自 function 的帮助文件ivreg。

library(AER)
data("CigarettesSW")
CigarettesSW$rprice <- with(CigarettesSW, price/cpi)
CigarettesSW$rincome <- with(CigarettesSW, income/population/cpi)
CigarettesSW$tdiff <- with(CigarettesSW, (taxs - tax)/cpi)

## model 
fm <- ivreg(log(packs) ~ log(rprice) + log(rincome) | log(rincome) + 
            tdiff + I(tax/cpi),
            data = CigarettesSW, subset = year == "1995")

summary()仅当您在模型上调用函数时才会计算 p 值、t 值等。这时候实际上函数summary.ivreg()被调用了。因此，如果您只需要获取 p 值，则应将结果保存summary()为某个对象。该对象（列表）包含几个部分，系数存储在名为 的矩阵中coefficients。

sum.res<-summary(fm)
names(sum.res)
 [1] "call"          "terms"         "residuals"     ""              "coefficients"  "sigma"        
 [7] "df"            "r.squared"     "adj.r.squared" "waldtest"      "vcov"

要获得所有系数：

sum.res$coefficients
               Estimate Std. Error   t value     Pr(>|t|)
(Intercept)   9.8949555  1.0585599  9.347563 4.120910e-12
log(rprice)  -1.2774241  0.2631986 -4.853461 1.496034e-05
log(rincome)  0.2804048  0.2385654  1.175379 2.460247e-01

p 值存储在该矩阵的 4. 列中：

sum.res$coefficients[,4]
 (Intercept)  log(rprice) log(rincome) 
4.120910e-12 1.496034e-05 2.460247e-01
于 2013-07-24T05:56:47.073 回答