Dear everyone,

I'm encountering an unknown error message when using gamm function:

>         fitoutput <-
gamm(cvd~as.factor(dow)+pm10+s(time,bs="cr",k=15,fx=TRUE)+s(tmean,bs="cr",k=7,fx=TRUE)
+
,correlation=corAR1(form=~1|city),family=poisson,random=list(city=~pm10),data=mimp)


 Maximum number of PQL iterations:  20
iteration 1
iteration 2
iteration 3
iteration 4
iteration 5
iteration 6
Error in `*tmp*`[[k]] : attempt to select less than one element

I never saw this error message before, and there is also no related
information about this error from the Internet. Hope some people who once
encountered this problem can tell me what's wrong.

Many thanks,

LC

-- 
Lung-Chang Chien
Doctoral candidate, Department of Biostatistics
Statistical assistant, School of Nursing
University of North Carolina at Chapel Hill
Office: Carrington Hall, RSC, Room 2014
TEL: 919-843-2480

	[[alternative HTML version deleted]]

