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

Question

我使用nls.lm（）函数来获取MODEL。后来，我写了＃34;摘要（MODEL）＆＃34;我得到参数列表，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)