Ich benutze die Funktion nls.lm(), um ein MODELL zu erhalten. Später schreibe ich "Zusammenfassung (MODEL)" und ich bekomme die Liste der Parameter, std.error ... und andere Konvergenz und Modelldetails.Wie schätzt/berechnet R den std.error in summary() eines nls.lm Modells?
Die Frage ist,
jemand weiß, wie R diese std.error (für jeden Parameter) mit Zusammenfassung gezeigt berechnet() ???
Danke!
Wenn Sie den Funktionsnamen eingeben, summary.nls.lm mit dem vorangestellten Paketnamen, minpack.lm und mit der :::-Funktion verbunden, wird der Code angezeigt.
minpack.lm:::summary.nls.lm
Dies ist der Abschnitt, der die Standardfehler
ibb <- chol(object$hessian)
ih <- chol2inv(ibb)
p <- length(param)
rdf <- length(object$fvec) - p
resvar <- deviance(object)/rdf
se <- sqrt(diag(ih) * resvar)
Danke !! Ich denke, dass ich mir helfen werde, weil diese Funktion (nls.lm) im minpack.lm Paket ist, können Sie es in R-Dokumentation sehen (https://cran.r-project.org/web/packages/minpack.lm /index.html). Wirklich, ich denke, dass nlsLM nls.lm aufrufen, um zu arbeiten. Grüße!! – Sergio
@Roland berechnet: Ich glaube, Sie haben durch die Lockerheit von R Namenskonventionen die Irre geführt worden, die nicht die Verbindung zwischen dem Punkt-Formalismus erzwingt in Funktion Versand. Es gibt keine "nls.lm" -Funktion in einem der Standardpakete, und die "nls" -Funktion würde kein Objekt der Klasse "lm" akzeptieren. Das "LM" steht für den Levenberg-Marquardt-Algorithmus und nicht für das "lineare Modell". –
@ 42- Nein, ich habe es nur mit 'nlsLM' verwechselt, was die Funktion dieser Pakete ist, die ich normalerweise verwende. – Roland