Models mixtes en R

Porto una bona tanda de models mixtes.

Una conclusió que he arribat: Son molt complicats.

Una altra el paquet nlme de R té tot alló com hom pot desitjar:

Si es vol crear un model de mesures repetides:
summary(lme(Y~X*temps,random=list(Subject=pdDiag(form=~1+temps)),correlation=corAR1(form=~temps|Subjecte),data=dades))


Continuarà,...
Ferran Carrascosa