コード 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
ブートストラップはエラーなしで完了し、妥当な結果が得られます。私の質問は、結果がどれほど信頼できるかです。