問題タブ [hessian-matrix]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
r - R: LTM: 標準化が失敗した場合、どうすれば奇妙な動作をするヘッセ行列を収束させることができますか?
R パッケージ ltm で段階的応答モデルを適合させようとしています。問題は、ヘッセ行列が収束しないことであり、その理由がわかりません。私が使用するコードは次のとおりです。
これにより、次のエラー メッセージが表示されます。「収束時のヘッセ行列には無限値または欠損値が含まれています。解が不安定です。
grmのヘルプページで提案されているように追加しました。ここでstart.val = "random"提案されているように、を使用して変数をスケーリングしようとしましたが、成功しませんでした。standardize()
コードはその変数がなくても正常に動作するため、犯人は変数 Apathy8 です。ヘッセ行列の出力も、最後の 5 列が壊れていることを示しています。
私がテストしたデータの多くのサンプルで問題が発生することは興味深いことです。私は1000個の観測を行っていますが、たとえば次のような60個のランダムな観測でも問題は解決しません( csvへのリンク):
これを解決する方法を知っている人はいますか？