从`zlogit $ sim()`(Zelig)中提取输出

时间:2015-06-24 03:22:14

标签: r r-zelig

我仍然在尝试习惯Zelig的开发人员正在努力的全新语法(在Zelig5中,有关安装当前开发版本here的说明)。感觉非常Pythonic,除了,不是......

无论如何,我只想存储sim练习的结果,但只能弄清楚如何将结果打印到控制台。

让我们使用引用的示例in the documentation(好吧,有点 - 更新以反映所见的Zelig5语法,例如here):

set.seed(1234)
library(Zelig) #Zelig_5.0-5
ztob<-ztobit$new() 
ztob$zelig(durable~age+quant,data=tobin)
ztob$setx(ztob)
ztob$sim()
summary(ztob)

 sim x :
 -----
ev
      mean        sd      50%      2.5%    97.5%
1 1.534273 0.6350075 1.451001 0.5103966 3.042459
pv
         mean       sd      50% 2.5%    97.5%
[1,] 3.002031 4.027547 1.310886    0 13.19713

我真的不知道pv的意思(没有真正记录),但我很确定我想要的预期值是1.53(在ev下,mean)。

有人能弄明白如何提取这个价值吗?我找不到summary.Zeligsummary.zelig之类的内容;我试过了:

  • summary(ztob)$ev / ztob$ev
  • print(summary(ztob))
  • summary(ztob)[1] / summary(ztob)[[1]]

什么吗?

1 个答案:

答案 0 :(得分:2)

在这种情况下,str是您的朋友。

您可以获得所有值:

x<-unlist(ztob[["sim.out"]][["x"]][["ev"]])

平均值:

mean(x)