R如何估計/計算nls.lm模型的summary（）中的std.error？

Question

我正在使用nls.lm（）函數來獲取MODEL。 后來，我寫了“摘要（模型）”，我得到了參數列表，std.error ......以及其他收斂和模型細節。

問題是，

有人知道R如何計算這些用summary（）顯示的std.error（對於每個參數）？

謝謝！

Answer 1

如果鍵入函數名稱， summary.nls.lm以包名稱minpack.lm並通過:::函數連接，則會看到代碼。

minpack.lm:::summary.nls.lm

這是計算標准誤差的部分

ibb <- chol(object$hessian)
ih <- chol2inv(ibb)
p <- length(param)
rdf <- length(object$fvec) - p
resvar <- deviance(object)/rdf
se <- sqrt(diag(ih) * resvar)