Comment imprimer la variance d'un lm dans R sans calculer à la main l'erreur standard?
Simple question vraiment! Je suis en cours d'exécution beaucoup de régressions linéaires de y~x
et veulent obtenir la variance pour chaque régression sans calcul de la part de la sortie d'Erreur Standard donné dans le summary.lm
de commande. Juste pour économiser un peu de temps :-). Toutes les idées de la commande pour le faire? Ou devrai-je écrire une fonction pour le faire moi-même?
m<-lm(Alopecurus.geniculatus~Year)
> summary(m)
Call:
lm(formula = Alopecurus.geniculatus ~ Year)
Residuals:
Min 1Q Median 3Q Max
-19.374 -8.667 -2.094 9.601 21.832
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 700.3921 302.2936 2.317 0.0275 *
Year -0.2757 0.1530 -1.802 0.0817 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 11.45 on 30 degrees of freedom
(15 observations deleted due to missingness)
Multiple R-squared: 0.09762, Adjusted R-squared: 0.06754
F-statistic: 3.246 on 1 and 30 DF, p-value: 0.08168
Si je reçois une sortie d'Erreur Standard et j'ai été l'espoir d'obtenir une Variance de sortie sans calculer à la main...
source d'informationauteur Sarah
Vous devez vous connecter pour publier un commentaire.
Je ne suis pas sûr de ce que vous voulez de la variance de.
Si vous voulez le résiduelle de la variance, c'est:
(summary(m)$sigma)**2
.Si vous voulez la variance de votre pentec'est:
(summary(m)$coefficients[2,2])**2
ouvcov(m)[2,2]
.vcov(m)
donne la matrice de covariance des coefficients – écarts sur la diagonale.
si vous faites référence à la norme des erreurs pour les estimations des coefficients, la réponse est
et si vous faites allusion à l'estimation de la variance résiduelle, c'est
Type
names( summary(m) )
etnames(m)
pour d'autres informations que vous pouvez accéder.