問題タブ [mlogit]

少なくとも 1 つの変数がクラスから独立して推定されるモデルを作成する際に問題があるため、すべてのクラスで 1 つの同じ係数を使用します。どうすればこれを行うことができますか？R パッケージ gmnl を使用しています。

クラスから独立して変数 pf、cl、loc、wk、tod、または seas の 1 つをどのようにモデル化しますか? ありがとうございました！

ネストされたロジット モデルでは、ツリーの各レベルでリグレッサーを定義できます。マニュアルと追加の例で読んだすべての例では、リグレッサーは最後のレベルに対してのみ定義されています。よく議論される釣りモードの例を使用します。

巣:

回帰子Price、CatchRate、およびがあるとしIncomeます。最後のレベルを説明し、最初のレベルを説明するにはPrice、どのように使用すればよいですか。CatchRateIncome

RI でできること:

しかし、変数をどこに貼り付けるかわかりませんincome。

3 つの連続変数と 5 つの名義変数 (各変数に 5 つのカテゴリ) を持つ 1 つの結果変数 (はい/いいえ) を持つデータ ファイル (100 万行) があります。モデルを構築するのにどのタイプの分析が適しているかを知りたかったのです。ロジット、プロビット、ロジスティック回帰を見てきました。分析に役立つ可能性が高い変数を何から始めて分析するかについて、私は混乱しています。

データ ファイル: 性別、地域、年齢、会社、専門分野、職務、診断、ラボ、注文、ステータス

M、西、41、PA、FPC、アシスタント、code18、27、3、はい

M、サウスウェスト、65、CV、FPC、ワーカー、コード18、69、11、なし

M、南、27、DV、IMC、アシスタント、無効、62、13、いいえ

M、サウスウェスト、18、CV、IMC、ワーカー、コード8、6、1、はい

PS: R 言語を使用します。どんな助けでも大歓迎ですありがとう！

以前は機能していましたが、パッケージのオンとオフで一日中問題が発生していましたが、それとは関係がないかもしれません. 以前に以下のコードを実行しましたが、うまくいきました。しかし、過去数時間、同じエラーが発生し続けます

「if (abs(x - oldx) < ftol) { のエラー: TRUE/FALSE が必要な場所に値がありません」</p>

代替の特定の定数なしでロジット モデルを実行しようとしています (ただし、生成しようとしている係数に関係なく、同じエラーが発生し続けます)。

ポインタをいただければ幸いです、ありがとう！

R でオンライン データセットを使用して単純な多項ロジ​​ット モデルを実行したいとします。

次の出力が得られます

summary(ml) 呼び出し: multinom(式 = educg ~ 性別 + prestigeg + reasong, data = データ)

Stata で同様のルーチンを実行すると、次のようになります。

次の出力が得られます。

これらの値が完全に異なるのはなぜですか? Rで多項ロジットモデルのStataのような出力を得るにはどうすればよいですか?

私はサッカーの試合のブックメーカーのオッズにおける市場効率の仮説を検証しようとしています. mlogit パッケージを使用して多項ロジット モデルを推定しました。

モデル: 結果=ログ(P1/Px)+ログ(P2/Px)

ここで、P1 はホームでの勝利の暗黙のブッキー確率、Px は引き分けの暗黙のブッキー確率などです。ドロー (x) は参照カテゴリです。

ここで、次の仮説に対して尤度ベースのテスト (LR、Wald、または LM) を使用したいと考えています。

H0: β1=(0,1,0)、β2=(0,0,1)

つまり、帰無仮説の下では、切片係数は両方の回帰で 0 です。ホームウィンのロジットの係数は、ｙ＝ホームウィンのとき１、ｙ＝アウェイウィンのとき０である。アウェイ勝利のロジットの係数は、y=ホーム勝利のとき0、y=アウェイ勝利のとき1である。

制約付きモデル (H0 モデル) を適合させる方法を理解するのに問題があります。そこから対数尤度を抽出して、LR テストで ML 推定モデルから受け取った同上と比較します。

ここの57ページの指示に従ってみました： https://cran.r-project.org/web/packages/mlogit/vignettes/mlogit.pdf

しかし、update() 関数を使用して H0 モデルを指定する方法がわかりません。出来ますか？

おそらく「オフセット」を使用して、nnet (multinom) パッケージを使用して同等のテストを行う方法を知っている場合は、その方法の説明も非常に高く評価されます。

助けてくれてありがとう！

Python または R でサンプルの重みが異なる多項ロジット推定を実行したいと思います (これらは調査の母集団の重みです)。statsmodels'MNLogitはこれを提供していないようです。Rでもそうではありませんmlogit。それを行う別の方法はありますか?