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