Forum de Statistiques
Vous souhaitez réagir à ce message ? Créez un compte en quelques clics ou connectez-vous pour continuer.
-17%
Le deal à ne pas rater :
SSD interne Crucial SSD P3 1To NVME à 49,99€
49.99 € 59.99 €
Voir le deal

Erreur standard du paramètre estimé

2 participants

Aller en bas

Erreur standard du paramètre estimé Empty Erreur standard du paramètre estimé

Message par coralie Mer 19 Mar 2014 - 8:56

Bonjour,
Tout le monde!
J'ai deux questions par rapport à ce sujet.
1. Dans un modèle de la régression linéaire avec variables qualitatives, j'ai calculé les erreur standards de mes trois paramètres estimés, mais je trouve que l'erreur standard des paramètres estimés  avec summary() et celle calculée par la matrice d'information de fisher est différent. Je voudrais savoir pourquoi
avec summary(glm)
Call:
glm(formula = yl ~ x1 + x2, family = binomial(logit), data = dfl)

Deviance Residuals:
    Min        1Q    Median        3Q       Max  
-2.02909  -0.02690   0.00002   0.06501   1.85812  

Coefficients:
              Estimate    Std. Error  z value Pr(>|z|)    
(Intercept) -79.2404   23.2708  -3.405 0.000661 ***
x1             2.1419     0.6588     3.251 0.001150 **
x2             5.1187     1.5124     3.385 0.000713 ***
calcul des erreurs standards avec la matrice d'information de fisher
> sqrt(diag(solve(V)))
[1] 20.9344641  0.5770681  1.3979011

2.
Pourquoi les paramètres eux peuvent avoir une erreur standard, ça veut dire que les paramètres peuvent aussi osciller? tongue 
s

Merci pour m'aider!

coralie

Nombre de messages : 2
Date d'inscription : 19/03/2014

Revenir en haut Aller en bas

Erreur standard du paramètre estimé Empty Re: Erreur standard du paramètre estimé

Message par droopy Mer 19 Mar 2014 - 9:25

bonjour,

il faudrait nous dire comment tu as calculée cette matrice V. Dans un glm il est important que tu prennes en compte les poids dans le calcul de V, est-ce que c'est le cas ?
Code:
# avec un jeu de données issu du net :
mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv")
mylogit <- glm(admit ~ gre + gpa + rank, data = mydata, family = "binomial")
summary(mylogit)
...
Coefficients:
            Estimate Std. Error z value Pr(>|z|)   
(Intercept) -3.449548  1.132846  -3.045  0.00233 **
gre          0.002294  0.001092  2.101  0.03564 * 
gpa          0.777014  0.327484  2.373  0.01766 * 
rank        -0.560031  0.127137  -4.405 1.06e-05 ***
...

W <- diag(mylogit$weights) # la matrice des poids
X <- model.matrix(mylogit) # la matrice des variables explicatives
sqrt(diag(solve(t(X)%*%W%*%X)))
(Intercept)        gre        gpa        rank
1.132846009 0.001091839 0.327483878 0.127136989

Ca correspond bien. Par contre si tu ne prends pas en compte ces poids tu risque d'avoir une différence.

2. tu fais l'hypothèse que les paramètres sont issus d'une certaine distribution (normale)
Cdlt
droopy
droopy

Nombre de messages : 1156
Date d'inscription : 04/09/2009

Revenir en haut Aller en bas

Erreur standard du paramètre estimé Empty Re: Erreur standard du paramètre estimé

Message par coralie Lun 24 Mar 2014 - 9:41

Smile Merci beaucoup pour votre aide!!!!!!
J'ai refait mon exercice, ça donne le même résultat maintenant
Cordialement!Bonne semaine!

coralie

Nombre de messages : 2
Date d'inscription : 19/03/2014

Revenir en haut Aller en bas

Erreur standard du paramètre estimé Empty Re: Erreur standard du paramètre estimé

Message par Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum