問題タブ [r-zelig]

私がこれを行うとき：

次のエラーが表示されます。

何か案は？

私は R パッケージ stargazer() を使用しています。スターゲイザーに Latex コードを作成してもらうことができない zelig モデルがあります。

model="normal" に変更すると、同じエラー メッセージが表示されます。model="gaussian" は、zelig() でサポートされているモデルではありません。stargazer のドキュメントには、zelig オブジェクトで動作すると書かれています。

私は初心者のRユーザーです。Zeligバージョン4.1-3とAmeliaIIバージョン1.7をインストールしました。Rパッケージと関数を使用して、結合された多重代入データの自由度、t統計量、およびf値を取得する方法に戸惑っています。

まず、AmeliaとZeligをロードしました。

次に、Ameliaに付属のサンプルデータをロードしました。

amelia関数を使用して、このデータセットに対して5つの代入を作成しました。

次に、代入を組み合わせるために、zelig関数を使用しました。

ただし、このコードを使用したときに、結合されたセットの係数ではなく、個々の代入の係数のように見える係数を取得しました。

それらは次のとおりでした：

アメリアのマニュアルには、Rを使用してすべての係数を取得する方法についての説明はありませんが、結合された多重代入データセットの係数の一部がどのようになるかが示されています（http://cran.r-project.org/web/の46ページを参照）。 packages / Amelia / vignettes / amelia.pdf）

アメリア関数はモンテカルロシミュレーションを使用するため、実行間のわずかな違いが予想されます。ただし、切片の大きな違いは、zelig関数が結合されたセット以外のものの回帰統計を返すという手がかりでした。

アメリアのマニュアルには、次のコードが記載されています。

使ってみました。返される結果は、46ページに示されている値と標準誤差に非常に近いものです。

ただし、t統計、自由度、またはf値は含まれていません。

手動計算を行わなくても自由度、t統計量、およびf値を取得できるように、Rで利用可能なオープンソースのパッケージまたは関数はありますか？

ありがとう。

乗算代入データ (Amelia で作成) でマルチレベル モデルを実行しようとしています。サンプルは、グループ = 24、N = 150 のクラスター化されたサンプルに基づいています。

このコードは、次のエラー コードを生成します。

OLS 回帰を実行すると、次のように機能します。

実例を提供させていただきます。

問題は、Zelig が Amelia の mi クラスとどのようにやり取りするかにあると思います。したがって、代替の R パッケージである lme4 に目を向けました。

結果は次のとおりです。

などで置き換えdiff[[5]]てもdiff[[4]]、結果は同じままdiff[[3]]です。それでも、これが実際に結合されたデータセットの結果なのか、1 つの帰属データセットの結果なのか疑問に思っています。何かご意見は？ありがとう！

以前の投稿のフォローアップとして、 (Amelia で作成された) 多重代入データを使用して、回帰のために zelig で統計モデルの相対的な品質の通常の尺度を取得する方法を学習することに興味があります。

質問

(1) AIC、F統計量、およびMIデータの自由度の値を取得する方法を知っている人はいますか?

(2)個々のデータセットの要約によって生成された測定値の単純な平均を取ることができると今井幸助が答えた同様の質問を見つけました。「おそらく」は、私を少し疑わしくさせます。これについて何か考えはありますか？

どうもありがとう！！

パッケージと同じようにMCMCglmm多かれ少なかれ推定されたモデルから、関心のある量をシミュレートしたいと思います。Zelig独立した値にZelig必要な値を設定すると、ソフトウェアが結果変数 (期待値、確率など) の結果を計算します。例：

ご覧のとおり、z = 10 の場合、y は約 17 です。

によって推定された事後分布で z = 10 をシミュレートしMCMCglmm、y の期待値を取得する方法はありますか? どうもありがとうございました！

Zelig Packageでの後続のモデリングのために、一致したデータセット ( MatchIt Packageを使用して、観測が傾向スコアで一致している場合) のリストを取得したいと考えています。

この例では、(t1とt2)、2 つの独立変数 (x1とx2)、および結果 ( y1) で照合する 2 つの処理があります。

まず、一致するデータ フレームのリストを作成します。

の各要素にlist.dfsは適切なクラスがあります。

[1] "matchit"

次に、このリストから要素が一致するデータ フレームを取得し、Zelig モデル オブジェクトのリストを作成します。

次のエラーが発生します。

Error in match.data(i) : object 'i' not found

しかし、これは明らかにリストと関係があります。同じ関数を の単一の要素に実行するとすべてが機能するためですlist.dfs。

class(zelig(y1 ~ x1 + x2, model = "logit", data = match.data(list.dfs[[1]])))

[1] "zelig" "logit"

私は何が欠けていますか？このリストの個別の項目で Zelig に作業させるにはどうすればよいですか?

背景線形回帰のためのAmelia出力でのZeligパッケージの使用については、 Amelia II のドキュメントで詳しく説明されています。

問題残念ながら、この出力を他の分析に使用する方法に関するドキュメントを見つけることができません。問題は、他のパッケージが Amelia の出力をどのように見るかにあるようです。

Zeligパッケージはすべての代入を 1 つのデータ セットとして扱うことができるように思われるため、その単一のデータ セットを他のパッケージ (たとえばdplyr )での操作に使用できるようにしたいと考えていますが、単一のデータ セットを取得することはできません。パッケージからのdata.frame。