Sara Mouro wrote: > > I have this 2 models that fit to my data: > > M3varI <- update (M3, weights=varIdent(form= ~ 1|SITE)) > > M3AR1<-update(M3,correlation=corAR1()) > > How can I put them toghether in one final model? > >
This looks like lme in package nlme, but it is pure guesswork, since your example is incomplete. Probably you should not use update, but something like lme(xxxxx~xxxx, random =....., weights=varIdent..., correlation=corAR1()) but I am not sure if this somewhat overloaded fit will converge. Dieter -- View this message in context: http://www.nabble.com/how-to-join-these-two-models--tp20516771p20516980.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.