从 `zlogit$sim()` (Zelig) 中提取输出
Extracting output from `zlogit$sim()` (Zelig)
我仍在努力适应 Zelig
的开发人员正在努力的全新语法(在 Zelig5 中,安装当前开发版本 here 的说明)。感觉很 Pythonic,除了,不是...
不管怎样,我只是想存储一个sim
练习的结果,但只能弄清楚如何将结果打印到控制台。
让我们使用引用的示例 in the documentation (well, sort of--updated to reflect the Zelig5 syntax seen, e.g., 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]]
有什么吗?
在这种情况下,str
是你的朋友。
您可以获得所有的值:
x<-unlist(ztob[["sim.out"]][["x"]][["ev"]])
平均数:
mean(x)
我仍在努力适应 Zelig
的开发人员正在努力的全新语法(在 Zelig5 中,安装当前开发版本 here 的说明)。感觉很 Pythonic,除了,不是...
不管怎样,我只是想存储一个sim
练习的结果,但只能弄清楚如何将结果打印到控制台。
让我们使用引用的示例 in the documentation (well, sort of--updated to reflect the Zelig5 syntax seen, e.g., 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]]
有什么吗?
在这种情况下,str
是你的朋友。
您可以获得所有的值:
x<-unlist(ztob[["sim.out"]][["x"]][["ev"]])
平均数:
mean(x)