nls; erreur dans nlsmodel

[ali bouali]

bonjour,
je suis un débutant de R, je travail avec la dernière version 2.5, j'essaye de faire des régressions. dans un première temps j'ai fait une régression linéaire ca marche,par contre au moment de régression non linéaire j'ai un message qui apparaitre :Erreur dans nlsModel(formula, mf, start, wts) :
singular gradient matrix at initial parameter estimates


mon code est:
Test<-read.table("R:/users/Linear_328.txt", header=TRUE)

library(quantreg)

library(stats)

DML2<-exp(DML)
head(DML2)
# Regression moindres carrés
Des<-lm(Age~DML2, data=Test)
init<-coef(Des)
names(init)<-c("b","a")
print(init)
Reg.nls<-nls(Age~a*log(DML2)+b, data=Test, start = init)
Reg.nls

quand je modifie les valeurs de start pour nls, le message n'apparu plus
le fichier de donnes est jont dans la liste
merci a l'avance de votre aide

DML Age
31 9.38709677419355
32 7.93548387096774
36 6.48387096774194
23 8.61290322580645
25 6.06451612903226
29 8.54838709677419
27 9.09677419354839
20 6.64516129032258
20 6.29032258064516
22 9.29032258064516
16 6.48387096774194
18 7.7741935483871
17 8.90322580645161
16 7.45161290322581
15 7.51612903225806
19 8.67741935483871
11 9.45161290322581
9 6.06451612903226
11 6.12903225806452
13 6.83870967741935
14 6.06451612903226
35 11.7096774193548
22 7.80645161290323
21 10.5806451612903
23 5.96774193548387
23 7.45161290322581
21 7.54838709677419
18 8.29032258064516
17 9.64516129032258
13 6.7741935483871
14 6.64516129032258
15 6.16129032258065
26 7.90322580645161
25 7.2258064516129
22 8.32258064516129
22 8.93548387096774
22 8.00
18 7.64516129032258
12 6.64516129032258
34 9.93548387096774
40 10.5483870967742
26 11.4516129032258
25 10.5806451612903
28 10.8709677419355
23 10.1612903225806
17 6.96774193548387
38 12.5806451612903
28 8.96774193548387
25 7.67741935483871

[Pierre Bady]

bonjour,

quand je modifie les valeurs de start pour nls, le message n'apparu plus
le fichier de donnes est jont dans la liste
merci a l'avance de votre aide

normal


@++

Pierre

PS: consulté la doc, svp:

Code : Tout sélectionner

?nls

[ali bouali]

merci pierre
en faite ce que je vous c'est que nls prend les valeurs de départ trouve par lm pour me faire le calcul de nls, mais il y tjs ce message d'erreur suivant:
Erreur dans nlsModel(formula, mf, start, wts) :
singular gradient matrix at initial parameter estimates

mon code est:
Test<-read.table("R:/users/Linear_328.txt", header=TRUE)
attach(Test)
library(stats)
DML2<-exp(DML)
head(DML2)
library(quantreg)
library(stats)
DML2<-exp(DML)
head(DML2)
# Regression moindres carrés
Des<-lm(Age~DML2, data=Test)
init<-coef(Des)
names(init)<-c("b","a")
print(init)
Reg.nls<-nls(Age~a*log(DML2)+b, data=Test, start = init)
Reg.nls

quand je modifie les valeurs de start pour nls, le message n'apparu plus
le fichier de donnes est jont dans la liste
merci a l'avance de votre aide

[ali bouali]

désole pour la répétition de même message

[Pierre Bady]

Rebonjour,

- Si vous avez le message " singular gradient matrix", c'est que les valeurs initiales sont mal choisies, donc il n'y a pas de convergence possible.

- Il me semble qu'il y a aussi un soucis avec les variables "DML" et "DML2" qui sont passé en "exp", puis retransformées en "log" .... à voir
En tout cas, ça me semble un peu confus.

- votre modèle non-linéaire ressemble furieusement à un modèle linéaire, non ? (y= a*X+b)


en espérant avoir aidé un peu

@+

Pierre


PS : désolé d’insister


une jolie doc sur le modèle non-linéaire sous S-plus :
http://www.stat.ucl.ac.be/cours/stat243 ... on_lin.pdf

[ali bouali]

merci de votre reponse