1

私は学生と一緒に作業し、GLMER でいくつかのモデルを実行していましたが、同じコードを使用すると、モデルは彼ではなく私のために収束することがわかりました。つまり、次のようなエラー メッセージが表示されます。 = 0.001、コンポーネント 1)

確認したところ、私は古いバージョンの R を実行していたため、パッケージを更新したときに、彼が 1.1-11 を実行していたときに LME4 が 1.1-7 までしか上がらなかったことに気付きました。そこで、R と LME4 を更新したところ、同じ結果が得られました。私たちが実行しているモデルに固有のものかもしれないと思っていましたが、現在、LME4 1.1-7 を実行している別の学生と一緒に作業しており、彼女のモデルは再び彼女に収束しますが、私には収束しません。1.1-11 で、これらの警告が表示される可能性が高くなるような変更はありますか? もしそうなら、変更の性質から、1.1-11 で警告が表示される理由についてのヒントを得ることができますか? 最後に、どの結果を信じますか?収束の警告に対処するためにできることはありますか?

4

1 に答える 1

2

LMEのCRANに関するニュース ページは、おそらく変更点を探すのに最適な場所です。

例えばこんなことに気付きました。

収束チェックの 1.1-8 勾配スケーリングでは、ヘッセ行列のコレスキー係数が使用されるようになりました。より正確ですが、これにより、追加の (おそらく誤検知の) 収束警告が発生します。

収束の警告を回避するために、できることはいくらでもあります。

  • 連続予測子の再スケーリング
  • さまざまなオプティマイザーを試す
  • 予測子の組み合わせでまばらなセルを探す
  • 入れ子になった変量効果が実験計画に必要かどうかを検討する

収束の失敗に関するこの相互検証済みの議論も参照してください。

于 2016-02-19T16:49:20.607 に答える