Bonjour,
Comment faire un graphique qui représente simultanément les betas et leurs IC issus d'un GLM family=poisson?
Merci d'avance
Modérateur : Groupe des modérateurs
Code : Tout sélectionner
library(sandwich)
cov.m1<-vcovHC (mod, type="HC0")
std.err<-sqrt(diag(cov.m1))
r.est<-cbind(estimate=mod$coefficients, std.err,
pvalue=round(2*(1-pnorm(abs(mod$coefficients/std.err))),5),
lower=mod$coefficients-1.96*std.err,
upper=mod$coefficients+1.96*std.err)
r.est
estimate std.err pvalue lower upper
(Intercept) 3.20128338 0.09936486 0.00000 3.00652825 3.39603851
lage30-40 ans 0.15844727 0.04669330 0.00069 0.06692841 0.24996613
lage40-50 ans 0.38685465 0.04606461 0.00000 0.29656802 0.47714128
lage50-60 ans 0.48844478 0.04655663 0.00000 0.39719379 0.57969578
lage60 ans + 0.55131438 0.10499819 0.00000 0.34551793 0.75711082
Retourner vers « Questions en cours »
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité