問題タブ [marginal-effects]

私はRを使用して研究を再現し、著者が報告したのとほぼ同じ結果を得ています。しかし、ある時点で、非現実的に小さいように見える限界効果を計算します。私の推論と以下のコードを見て、私がどこかで間違っているかどうかを確認していただければ幸いです。

私のサンプルには24535個の観測値が含まれており、従属変数「x028bin」は値0と1をとるバイナリ変数であり、さらに10個の説明変数があります。これらの独立変数のうち9つには数値レベルがあり、独立変数「f025grouped」はさまざまな宗教宗派からなる要素です。

宗教宗派のダミーを含むプロビット回帰を実行してから、限界効果を計算したいと思います。そのためには、最初に欠落している値を削除し、従属変数と独立変数の間のクロスタブを使用して、小さいセルまたは0個のセルがないことを確認します。次に、正常に動作するプロビットモデルを実行すると、妥当な結果も得られます。

ただし、プロビット係数とスケールファクターからすべての変数を平均して限界効果を計算する場合、得られる限界効果は小さすぎます（例：2.6042e-78）。コードは次のようになります。

データセットが大きすぎるため、実際の例を提供できないことをお詫び申し上げます。コメントをいただければ幸いです。どうもありがとう。

一番、

トビアス

R の効果パッケージを使用して、切片を含まない線形モデルの (限界) 効果をプロットできますか?

このようなモデルで試した R コードを次に示します。

それが生成するエラーは次のとおりです。

このエラーを克服する方法について誰かがアイデアを持っている場合は、私に知らせてください。

ありがとう、

イザベラ

私はRが初めてで、その計算の論理をまだ完全には理解していません...

以前の投稿の助けを借りても問題を解決できません。

11 の変数に対して約 600 の観測値のデータ セットがあります。mean()多項モデルを正常に実行しましたが、コマンドが NA を取得しているため、限界効果を達成できません。

データセット:

データを変換するためのコマンド:

平均値の行列のコマンド:

私が得る出力：

このコマンドの NA の理由を教えてください。

ほとんどがNAであるDEBTの正しい出力が得られない場合は理解できますが、RSVS、CHINNなどの平均を計算しないのはなぜですか?

RGDP の平均を計算するのはなぜですか? 両方の変数は次のとおりです。

この問題を克服するにはどうすればよいですか？

アップデート：

David Arenburg のコメントのおかげで、平均を計算することができました。

ただし、限界効果の計算中に発生する別のエラーがあります。

手段を 5 回コピーして、6 つの選択肢に対して 6 行の行列 zz を取得します。

Effects() コマンドは次の出力を提供します。

ただし、この出力について結論を出すことはできません。これらの限界効果をどの変数に割り当てることができるかわかりません。

1つずつ実行しようとすると、別のエラーが発生します：

モデル コマンド:

私が使用しているパッケージ：

事前にありがとう、Zyta

R で mlogit モデルを正常に実行した後、限界効果を取得しようとすると、次のようなエラーが表示されます。

別の投稿で説明されているように、ソース コードの行番号 16 を変更しようとしても、同じエラーが発生します。どんな助けでも大歓迎です。私の変数はすべて代替固有のものではありません。私には4つの選択肢があります。

以下は、Web に投稿された同様の質問へのリンクですが、そこに投稿された解決策に苦労しています。 R における mlogit の限界効果 よろしくお願いし ます

mlogit() パッケージを使用して、R でネストされた Logit モデルを実行しました。現在、限界効果/弾力性を測定しようとしていますが、引き続きエラーが発生します。ここでは、パッケージの作成者がビネットを変更してエラーを再現しました。

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

これは、ネストされたモデル (通常の多項ロジットなど) がない場合に正常に機能し、以前のスタックオーバーフローの質問でカバーされていますが、変更されたデータ フレーム (でこの場合平均フレーム z)。

ここでの解決策は次のとおりです。Rでのmlogitの限界効果は役に立ちませんでした。