Bonjour,
J’ai quelques problèmes avec la fonction lme sous R et sous Splus.
J’étudie la relation allométrique entre le diamètre des branches et la biomasse soutenue par ces branches chez le manguier.
Je voudrais regarder l’influence de la variété sur la relation entre le diamètre et la biomasse.
J’ai des données pour 6 variétés. Environ 30 données par variétés pour 4 variétés et 10 par variétés pour 2 variétés.
Sous R le modèle ne converge pas
lme1<-lme(lmstoter~cldbm,random=~cldbm|variete)
Error in lme.formula(lmstoter ~ cldbm, random = ~cldbm | variete) :
iteration limit reached without convergence (9)
avec lmstoter=log(biomasse totale de l’ensemble raméal )
cldbm=log(diamètre à la base de l’ensemble raméal ), centré par la moyenne des log(dbm)
Sous S il converge mais donne une corrélation entre l’intercept et cldbm de –1.
Je ne trouve pas ce qui dans les données peut poser problème.
Est ce quelqu’un pourrait m’aider à comprendre ce qui ne va pas ?
Merci par avance
Claire