[R] Problems with corARMA

From: Pamela McCaskie <pmccask_at_cyllene.uwa.edu.au>
Date: Fri 10 Jun 2005 - 15:43:18 EST

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 argument with differerent correlation structures and it has worked fine. Does anyone know why it won't recognise ~year | id2 (or even ~ 1 | id2) as a formula?

Any help would be great

Pamela A McCaskie
BSc(Mathematical Sciences)(Hons)

Western Australian Institute for Medical Research
University of Western Australia
SCGH Campus
Ground Floor, B Block
QE-II Medical Centre
Hospital Avenue, Nedlands
Western Australia  6009
Email:        pmccask@cyllene.uwa.edu.au
Phone:        +61-8-9346 1612
Mob:          0417 926 607

R-help@stat.math.ethz.ch mailing list
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Fri Jun 10 15:51:24 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:28 EST