Glm $ Residuen und Resid (Glm) verstehen

8

Können Sie mir sagen, was von glm $ residues und resid (glm) zurückgegeben wird, wobei glm ein Quasipoisson-Objekt ist. z.B. Wie würde ich sie mit glm $ y und glm $ linear.predictors erstellen.

Glm $ Residuen

%Vor%

37715 10042 2174 -0,2574 -2,7538 -2,2661 -1,4480 -0,4381 0,7542 1,9845 2,7749

niedrigste: -4.243 -3.552 -3.509 -3.481 -3.464 Höchster: 8,195 8,319 8,592 9,089 9,416

resid (glm)

%Vor%

-0,2080 0,4106 1,1766 1,7333

niedrigste: -1,0000 -0,8415 -0,8350 -0,8333 -0,8288 Höchste: 7.2491 7.6110 7.6486 7.9574 10.1932

    
Michael Bishop 28.03.2010, 00:40
quelle

2 Antworten

18

Das Aufrufen von resid (Modell) wird standardmäßig auf die Devianz-Residuen angewendet, während das Modell $ resid die Arbeitsresiduen liefert. Aufgrund der Verknüpfungsfunktion gibt es keine eindeutige Definition dessen, was ein Modellrest ist. Es gibt die Devianz-, Arbeits-, Teil-, Pearson- und Response-Residuen. Da diese nur auf der mittleren Struktur (nicht der Varianz) beruhen, haben die Residuen für das Quasipoisson und das Poisson die gleiche Form. Sie können die Funktion residuals.glm für Details betrachten, aber hier ist ein Beispiel:

%Vor%     
Ian Fellows 28.03.2010, 05:39
quelle
4

Ich weiß nicht genug über Poisson- und Quasi-Poisson-Verteilungen, um Ihre Frage in der geforderten Tiefe zu beantworten (dh eine exakte Gleichung, die die Variablen unter Verwendung des Modells in die Residuen umwandelt), aber wenn eine der Verwirrungen ist aufgrund dessen, welche Resttypen verwendet werden und warum die beiden Befehle eine andere Antwort geben, könnte dies helfen:

resid () verwendet standardmäßig einen "deviance" -Typ in R. Allerdings weist glm () dem $ residals-Vektor verschiedene Residuen zu.

Wenn Sie die Quasi-Poisson-Familie verwenden, weist glm () Residuen des Arbeitstyps zu, während resid () den Abweichungstyp als Standard angibt.

Um dies auszuprobieren, können Sie Folgendes verwenden:

  

resid (glm, type="arbeiten")

und

  

glm $ residuen

und das sollte Ihnen die gleiche Antwort geben (zumindest bei einem Beispieldatensatz, den ich verwendet habe).

Nach R sind Arbeitsresiduen: "Die Residuen in der letzten Iteration des IWLS passen"

Wenn Sie das Buch "Verallgemeinerte lineare Modelle und Erweiterungen" (von Hardin und Hilbe) auf googlebooks nachschlagen, können Sie Abschnitt 4.5 aufrufen, in dem die verschiedenen Arten von Residuen erläutert werden.

    
Adam 28.03.2010 02:36
quelle

Tags und Links