[R] nlme: "Deficient rank in gls_loglik" when creating corAR1()

From: David Hugh-Jones <davidhughjones_at_gmail.com>
Date: Wed 04 May 2005 - 01:53:40 EST


I have a bunch of data which is structured by year and US state, so I have created a nlme groupedData object for it:

formula(gd2)
DEPVAR ~ YEAR | ABREV Now I am trying to run a gls regression on it. I want the error correlation structure to be AR1 with a different rho for each state, so I do

> mdyn.1.1 = gls(model = DEPVAR ~ BLAH + BLAH, data=gd2, corr=corAR1(form= ~ YEAR | ABREV),na.action=na.omit)

YEAR and ABREV are always present; DEPVAR is absent for one state.

I get the following error message:

Error in logLik.glsStruct(glsSt, glsPars) :

        Deficient rank in gls_loglik

Can anyone enlighten me? The error message goes away if I just do corAR1(form = ~1), but this is not meaningful for my data.

Cheers
David



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed May 04 02:32:00 2005

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