問題タブ [logistic-regression]

解決しようとしている問題がありますが、成功していません。2日以上検索しましたが、手がかりは1つもありませんでした。答えがそこにあり、私が見つけられなかった場合は申し訳ありません。

数年前に推定した古いモデルからのロジスティック方程式回帰 (バイナリ モデル) があるとします。したがって、パラメータ βk (k = 1, 2, ..., p) は過去に推定されているため、既知です。しかし、モデルを適合させるために使用されたデータがありません。

私の質問は、R でこの古い推定ロジスティック モデルをオブジェクト (ロジスティック回帰モデルに対応) として導入できますか?

「予測」機能を使用して、このロジスティック回帰を新しいデータ セット (現在のデータ) で証明し、この古いモデルの有効性を時の試練に耐えられるかどうかを確認したいと思います。この関数を使用するには、ロジスティック回帰モデルのオブジェクトが必要です。

事前にどうもありがとうございました。

次の回帰を実行するときに、変数の下のすべてのカテゴリの N が必要です。

logit RIC i.BMF2 Kerosene Male ib1.agegroup3 ib5.Education2 ib3.Occupation2 i.Household_monthly Smoker, or

各カテゴリに数人しかいないのか、それとも多くの人がいるのかを知ることができるように. 各変数に対してこれらの行に沿って何かを行うことができることは理解していtab RIC variable, by(Sex)ますが、それが回帰が使用している実際の人数であることを保証するものではありません回帰、正しいですか?)

私は正則化されたLRで手を試しています.matlabのこの式で簡単です:

コスト関数:

勾配:

これは matlab コードではなく、単なる数式です。

これまでのところ、私はこれを行ってきました：

そして、コスト関数は正しい結果を与えていますが、勾配（1ステップ）がそうでない理由がわかりません。コストは正しいJ = 0.6931を与え、勾配grad = 0.3603 -0.1476 0.0320を与えます。パラメータ theta(1) を正則化する必要がないため、2 から コードに何か問題があると思いますが、4日経っても表示されません。ありがとう

7 日間の広告期間が 5 日間の広告期間よりも優れているかどうかをテストしようとしているデータ セットがあります。これをテストするには、ロジスティック回帰が最適な方法だと思います。テストをそれぞれ 2 週間実行しました。トラフィック、サインアップ、減少などのデータがあります。

データは次のようになります。

1 ははい、0 はいいえを意味します。

Rでこのコードを実行しました：

それから

私の出力は次のとおりです。

また

1 または 0 を伝えるために 1 つの出力のみを取得するように実行するにはどうすればよいですか (IE では、7 日間で全体の増加率が 5% を超えていますか?)

ありがとう

実行例として、多項ロジットに関する UCLA の例を使用してみましょう---

どうすれば 95% 信頼区間を取得できますか?

混合効果のロジスティック回帰モデルを実行しようとしていますが、各クラスター/グループの可変サンプル サイズと、一部のモデルの「成功」数が非常に少ないことを懸念しています。

2004 年から 2011 年まで毎年訪れた 163 のフィールド プロット (つまり、クラスター/グループ) に約 700 本の木が分散しています。私は、研究の各年に個別の混合効果ロジスティック回帰モデル (以下、GLMM) を適合させて、この出力を共有虚弱モデル (つまり、ランダム効果による生存分析) からの推論と比較しています。

プロットあたりの木の数は 1 ～ 22 です。また、年によっては、「成功」(つまり、病気の木) の数が非常に少ない年もあります。たとえば、2011 年には、694 回の「失敗」(つまり、健全な木) のうち、成功したのは 4 回だけでした。

私の質問は次のとおりです: (1) 推論の焦点が GLMM の固定効果の推定のみにある場合のサンプル|グループの理想的な数の一般的な規則はありますか? (2) GLMM にそのような極端な違いがある場合、GLMM は安定していますか?成功と失敗の比率。

ソースのアドバイスや提案をありがとう。

-サラ

Incanter データ分析ライブラリを使用して、単純なロジスティック回帰の例を Clojure に実装しようとしています。シグモイド関数とコスト関数のコーディングは成功しましたが、Incanter の BFGS 最小化関数がかなりの問題を引き起こしているようです。

(minimize (cost-func X y) (matrix [0 0]))関数を指定してパラメーターを開始しようとminimizeすると、REPL がエラーをスローします。

最小化機能が正確に何を期待しているのかについて、私は非常に混乱しています。

参考までに、すべて Python で書き直したところ、同じ最小化アルゴリズムを使用して、すべてのコードが期待どおりに実行されました。

出力中

Python コードが正常に実行できる理由がわかりませんが、私の Clojure 実装は失敗します。助言がありますか？

アップデート

新しいエラーをスローするminimize導関数を計算しようとしていたため、ドキュメント文字列を読み直しました。cost-func

を使用transして 1xn col マトリックスを nx1 行マトリックスに変換すると、反対のエラーで同じエラーが発生します。

:asize [1 2], :bsize [2 1]}

私はここでかなり迷っています。