modèle non linéaire a effet mixte.

Bonjour ,

J'ai le modèle non linéaire suivant :

Code:

Modèlenl1<-nlmer(y~a*(Z)*exp(-b*Z),start=c(a=0.1,b=0.1),data=MergeTestP11)


je voudrais appliquer un effet mixte(par groupe g) sur les paramètres(a et b) avec la fonction "nlmer" mais j'arrive pas a trouver la bonne écriture ,quelqu'un a une idée ?
Merci beaucoup d'avance.

autrement comment introduire effet mixte sur ce type de modèle non linéaire :

Code:

Modèlenl1<-nls(y~a*(Z)*exp(-b*Z),start=c(a=0.1,b=0.1),data=MergeTestP11)

?

Salut,

Tout est dit dans l'aide de la fonction...
Qqchose du type :

Code:

Modèlenl1<-nlmer(y~a*(Z)*exp(-b*Z)~(1|groupe) ,start=c(a=0.1,b=0.1),data=MergeTestP11)

Là ce serait pour un modèle de type "random intercept model"
Il te faut un facteur de regroupement des observations par groupe.

Nik

ça ne marche pas en plus j'ai pas compris l'erreur affiché :

Code:

Modèlenl1<-nlmer(DV_max~a*(Z)*exp(-b*Z)~(a|FELD_NR) ,start=c(a=0.1,b=0.1),data=MergeTestP11)
Erreur dans nlmer(DV_max ~ a * (Z) * exp(-b * Z) ~ (a | FELD_NR), start = c(a = 0.1,  :
  gradient attribute of evaluated model must be a numeric matrix


On ne met pas un paramètre à estimer dans la part aléatoire du modèle. Cette denière permet de spécifier une structure de la matrice de variance-covariance. Il y a donc déjà suffisamment de choses à estimer comme ça.

Es tu sûr de ce que tu fais ? Les modèles mixtes ne sont pas une chose facile à appréhender donc si tu ne connais pas les bases tu va au devant de sérieuses erreurs.

Nik