2

我仍在努力适应开发人员正在努力的全新语法Zelig(在 Zelig5 中,在此处安装当前开发版本的说明。感觉很 Pythonic,除了,不是……

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

让我们使用文档中引用的示例(好吧,有点更新以反映看到的 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.Zeligor summary.zelig; 我试过了:

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

任何事物?

4

1 回答 1

2

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

您可以获得所有值:

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

意思是:

mean(x)
于 2015-06-24T03:42:22.577 回答