[NLS] NaN ?!

[Coralie Sellos]

Bonjour,

Décidément la fonction nls me donne des difficultés...(je commence sur R)

nls (Yobs ~ model1(x1, x2, a, b, sigma), start=c(a=1, b=3, sigma=0.6), trace=T)
70121.53 : 1.0 3.0 0.6
Erreur dans numericDeriv(form[[3]], names(ind), env) :
Valeur manquante ou infinie obtenue au cours du calcul du modèle
De plus : Warning message:
In log(sigma) : production de NaN

Je ne comprends pas la 'production de NaN' ??

D'avance merci de votre aide (encore)

[Pierre Bady]

bonjour,

sigma semble négatif.

HTH

pierre

[Coralie Sellos]

Oui mais sigma n'est pas négatif... sigma =0.4
Donc c'est bien pour ca que je ne comprends pas le problème

[Pierre Bady]

re

> log(0.4)
[1] -0.9162907
> log(-1)
[1] NaN
Warning message:
In log(-1) : production de NaN
>

il faut peut-être revoir vos conditions initiales.

Sans information sur model1, c'est relativement difficile de comprendre l'apparition des NaN.


HTH

pierre

[Coralie Sellos]

Oui ca je le sais ...
Mais moi je fais pas log(-1)
log(0.4) c'est un nombre je vois pas pourquoi il calcule pas ...

[Pierre Bady]

re...

lors de l'estimation de vos paramètres, il se peut que sigma devienne négatif
(par exemple si les conditions initiales sont mal choisies).

D'après ce que l'on peut voir de votre premier post, sigma n'est pas un nombre "défini", mais un paramètre à estimer.


HTH

pierre

[Coralie Sellos]

C'est ce que je me suis dit ... Du coup j'ai fait le calcul sur un autre logiciel pour vérifier que ma fonction était ok et ca marche
Bon c'est pas grave je vais essayer autre chose

Désolée du dérangement et merci ...