問題タブ [covariance-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.
covariance-matrix - rmvnorm 警告メッセージ: sigma は数値的に半正定値ではありません
コード 1 で収束したモデルに適合した平均と共分散行列を使用して、パラメトリック ブートストラップを設定しようとしてnlm()
います。
rmvnorm(1,mean=mod$estimate,sigma=solve(mod$hessian))
次の警告メッセージが表示されます。
Warning message:
In rmvnorm(1, mod$estimate, solve(mod$hessian)) :
sigma is numerically not positive semidefinite
ブートストラップはエラーなしで完了し、妥当な結果が得られます。私の質問は、結果がどれほど信頼できるかです。