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

Question

我仍在努力适应 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.Zelig 或summary.zelig 这样的东西；我试过了：

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

有什么事吗？

Answer 1

在这种情况下，str 是你的朋友。

你可以得到所有的值：

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

平均数：

mean(x)