2012-06-19 16 views
20

Wir haben ein lm Objekt aus und wollen die StandardfehlerR: Standardfehlerausgabe aus lm Objekt

lm_aaa<- lm(aaa~x+y+z) 

ich die Funktion Zusammenfassung, Namen und Koeffizienten kennen zu extrahieren. Zusammenfassung scheint jedoch die einzige Möglichkeit zu sein, manuell auf den Standardfehler zuzugreifen. Habt ihr eine Idee, wie ich das einfach ausgeben kann?

danke!

Antwort

16

Der Ausgang von der summary Funktion ist nur ein R Liste. Sie können also alle Standardlistenoperationen verwenden. Zum Beispiel:

#some data (taken from Roland's example) 
x = c(1,2,3,4) 
y = c(2.1,3.9,6.3,7.8) 

#fitting a linear model 
fit = lm(y~x) 
m = summary(fit) 

Das m Objekt oder Liste hat eine Reihe von Attributen. Sie können sie Zugriff auf die Konsole oder den Namen Ansatz:

m$sigma 
m[[6]] 

Eine praktische Funktion kennen heißt, str. Diese Funktion bietet eine Zusammenfassung der Objekte Attribute, dh

str(m) 
+0

Allerdings bezieht sich @csgillespie auf die ** Rest ** Standardabweichung des Modells, nicht auf die Standardabweichung der einzelnen Koeffizienten. Die Funktion 'm $ sigma' entspricht' sigma (fit) ', siehe [hier] (https://stat.ethz.ch/R-manual/R-devel/library/stats/html/sigma.html). Ich glaube, die Frage war wirklich die Standardabweichung der ** individuellen ** Koeffizienten. –

8
#some data 
x<-c(1,2,3,4) 
y<-c(2.1,3.9,6.3,7.8) 

#fitting a linear model 
fit<-lm(y~x) 

#look at the statistics summary 
summary(fit) 

#get the standard error of the slope 
se_slope<-summary(fit)$coef[[4]] 
#the index depends on the model and which se you want to extract 

#get the residual standard error 
rse<-summary(fit)$sigma 
11

Um eine Liste der Standardfehler für alle Parameter zu erhalten, haben Sie

summary(lm_aaa)$coefficients[, 2] 

Wie andere darauf hingewiesen, verwenden können, str(lm_aaa) erzählen Sie sind so ziemlich alle Informationen, die aus Ihrem Modell extrahiert werden können.

0

Wenn Sie nicht den Standardfehler/Abweichung desModell erhalten wollen, sondern die Standardfehler/Abweichung der einzelnen Koeffizienten verwenden

# some data (taken from Roland's example) 
x = c(1, 2, 3, 4) 
y = c(2.1, 3.9, 6.3, 7.8) 

# fitting a linear model 
fit = lm(y ~ x) 

# get vector of all standard errors of the coefficients 
coef(summary(fit))[, "Std. Error"] 

Für weitere Informationen über die Standardfehler/Abweichung des Modells, siehe here. Weitere Informationen zum Standardfehler/Abweichung der Koeffizienten finden Sie unter here.