Displaying 1 result from an estimated 1 matches for "gnls4a".
Did you mean:
gnls
2005 Jun 10
1
Problems with corARMA
Dear all
I am tryiing to fit the following lme with an ARMA correlation structure:
test <- lme(fixed=fev1f~year, random=~1|id2, data=pheno2,
correlation=corARMA(value=0.2, form=~year|id2), na.action=na.omit)
But I get the following error message:
Error in getGroupsFormula.default(correlation, asList = TRUE) :
"Form" argument must be a formula
I have used this same form