Bonjour,
Je suis actuellement occupé de travailler avec la fonction lme() dans R (package nlme).
Lors de l'appel de la fonction, j'ai des problèmes de convergence.
J'utilise la commance control = list(tolerance = 0.001, MaxIter = 500) afin d'essayer de pallier à ce problème.
Cependant, le rajout de cette comande ne fait aucun effet. (Pourtant, je l'ai trouvé sur des tutoriels...)
Quelqu'un pourrait-il m'éclairer?
Merci.
quand ça ne converge pas, dans la plupart des cas c'est que le modèle est mauvais. Il faut s'interroger sur les hypothèses sous-jacentes, il faut essayer d'enlever des termes, peut-être aussi d'en ajouter, il faut essayer diverses transformations des variables...
Oui, tu as raison.
Je teste en fait la propriété d'un nouveau type de modèle (fractional polynomial) pour modéliser certaine étude.
J'ai voulu appliquer l'effet aléatoire sur tout les paramètres, ce qui était fort présomptueux.
Je me contente, dans un premier temps, de considérer l'intercept comme effet aléatoire. Je verrai bien si ca me donne des résultats pertinents.
Dans le cas contraire, j'avancerais petit à petit dans les hypothèses.
Merci pour te réponse!
Vince.
ça me semble une démarche saine. De toutes façons les modèles trop complexes (par exemple avec des interactions d'ordre 3 ou 4) on ne sait pas les interpréter.
