Calcul intervalle de confiance de paramètres estimés

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Julie Vercelloni
Messages : 27
Enregistré le : 01 Avr 2010, 23:29

Calcul intervalle de confiance de paramètres estimés

Messagepar Julie Vercelloni » 17 Juin 2011, 13:43

Bonjour,

A partir de régressions non linéaires, j'essaye de calculer les intervalles de confiance de mes paramètres estimés.

A partir de cette commande:

Code : Tout sélectionner

model_pocmat2 <- nls(SPMat ~ SSgompertz(anmat, Asymmat, b2mat, b3mat),
                 data = wmat,trace=TRUE)
summary(model_pocmat2)


J'otiens:

Code : Tout sélectionner

Formula: SPMat ~ SSgompertz(anmat, Asymmat, b2mat, b3mat)

Parameters:
        Estimate Std. Error t value Pr(>|t|)   
Asymmat  44.7411    16.2394   2.755  0.11036   
b2mat     6.5982     1.6884   3.908  0.05968 .
b3mat     0.7891     0.0621  12.706  0.00614 **
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 1.019 on 2 degrees of freedom

Number of iterations to convergence: 1
Achieved convergence tolerance: 9.89e-07


Pour calculer les intervalles de confiance je fais:

Code : Tout sélectionner

confint(model_pocmat2)


Et la, j'otiens le message d'erreur suivant:

Error in prof$getProfile() : number of iterations exceeded maximum of 50.

Je n'arrive pas à trouver sur le net ce que ça signifie. Est ce que ça pourrait être du au fait que j'ai un paramètre de mon modèle qui n'est pas significativement différent de 0 ? Pourquoi même si je spécifie un paramètre significativement différent de 0 comme par exemple:

Code : Tout sélectionner

confint(model_pocmat2, b3mat)


J'obtiens le même message d'erreur !

Avez vous des solutions, idées, suggestions à mon problème ?

En vous remerciant
Cordialement

Julie Vercelloni

Renaud Lancelot
Messages : 2484
Enregistré le : 16 Déc 2004, 08:01
Contact :

Messagepar Renaud Lancelot » 17 Juin 2011, 15:50

Voir http://finzi.psych.upenn.edu/R/Rhelp02/ ... 38295.html et autres messages de ce genre, mais ça risque d'être décevant. Vous fournissez - éventuellement sans le savoir via profile() qui est appelé par confint, des valeurs des paramètres pour lesquelles l'algorithme ne converge pas.

Essayer peut-être d'augmenter le nb d'itérations max dans nls (voir ?nls et ?nls.control). Je doute que ça marche :-(
Renaud

Matthieu Lesnoff
Messages : 120
Enregistré le : 29 Nov 2004, 12:41

Re: Calcul intervalle de confiance de paramètres estimés

Messagepar Matthieu Lesnoff » 17 Juin 2011, 17:52

Julie Vercelloni a écrit :A partir de régressions non linéaires, j'essaye de calculer les intervalles de confiance de mes paramètres estimés. Pour calculer les intervalles de confiance je fais:

Code : Tout sélectionner

confint(model_pocmat2)



La fonction confint calcule des IC de type "profile likelihood" ; si besoin et en alternative, vous pouvez éventuellement vous retrancher vers des IC supposant la normalité asymptotique des estimateurs des paramètres. Par exemple :

Code : Tout sélectionner

z <- coef(summary(model_pocmat2))
# z est une matrice
# l'IC95% des paramètres est donné par :
cbind(
  z[, 1] - 1.96 * z[, 2],
  z[, 1] + 1.96 * z[, 2]
  )
ML

Julie Vercelloni
Messages : 27
Enregistré le : 01 Avr 2010, 23:29

Messagepar Julie Vercelloni » 20 Juin 2011, 12:56

Ok, encore un grand merci pour vos conseils !

Julie Vercelloni

Renaud Lancelot
Messages : 2484
Enregistré le : 16 Déc 2004, 08:01
Contact :

Messagepar Renaud Lancelot » 21 Juin 2011, 10:57

Est-ce que ça marche en augmentant le nb max d'itérations avec nls.control ?
Renaud

Julie Vercelloni
Messages : 27
Enregistré le : 01 Avr 2010, 23:29

Messagepar Julie Vercelloni » 22 Juin 2011, 11:46

Bonjour Renaud,

Oui j'ai essayé le nls.control en augmentant le nombre d'itérations mais j'obtiens le même message d'erreur...

Merci pour vos conseils

Julie


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité