我仍然在尝试习惯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]]
什么吗?
答案 0 :(得分:2)
在这种情况下,str
是您的朋友。
您可以获得所有值:
x<-unlist(ztob[["sim.out"]][["x"]][["ev"]])
平均值:
mean(x)