2

私はまだ、開発者Zeligが取り組んでいるまったく新しい構文に慣れようとしています (Zelig5 では、現在の開発バージョンのインストール手順はこちら)。非常にPythonicに感じますが、そうではありません...

とにかく、演習の結果を保存したいだけですがsim、結果をコンソールに出力する方法しかわかりません。

ドキュメントで引用されている例を使用してみましょう(まあ、一種の -- 見られる Zelig5 構文を反映するように更新されます。たとえば、こちら):

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]]

なんでも？

4

1 に答える 1

2

このような場合、strあなたの味方です。

すべての値を取得できます。

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

そして平均：

mean(x)
于 2015-06-24T03:42:22.577 に答える