問題タブ [gnm]

を使用してロジスティック回帰がglmあり、次の形式の用語を追加したいと思います

c ₁ (k+ac ₂ )/(t+c ₂ )

ここで、k と t はデータ フレームの列、a は定数です。_{R に c 1}と c ₂の最適な値を見つけてもらいたいと思います。これは可能ですか？

固定値のみが必要な場合は、c ₂ = 2 とします。

c ₁ (k+2a)/(t+2)

私はちょうど書くことができました

これは私が今していることと似ています。しかし、2が最適であるとは思わず、「手動で」反復するのは非常に時間がかかります.

私はgnmR の関数を使用しています。列が各時間間隔でその共変量のバイナリ インジケーターである 6 つの行列を使用して、離散化した整数部分でモデルを適合させようとしています。これらの各列にはa1, a2、マトリックス 1 の場合b1, b2は など、マトリックス 2 の場合は などのラベルが付けられます。これらの行列ごとに約 250 の列があります。

これらの行列を文字列として gnm コマンドに渡します。

ここexprで、 はすべての行列列変数名のリストを含む文字列です (つまりa1, a2, ..., b1, b2, ....)。

私のコードは、5 つの行列を使用すると適切な結果を生成しますが、6 つの行列を使用しようとすると、gnmコマンドを実行しようとすると次のエラーが発生します。

これをバイパスする方法はありますか？

一般化された Lee-Carter モデルを R パッケージ gnm に当てはめようとしています。データフレームを構築しました

さて、これが機能している間

これに移動すると、エラーが発生します

具体的には、

数式で「D/E」を「base$D/base$E」に置き換えようとしましたが、うまくいきませんでした...