Lung-Chang Chien
2008-May-21 05:16 UTC
[R] an unknown error message when using gamm function
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]]
What version of mgcv are you using, please? On Wednesday 21 May 2008 06:16, Lung-Chang Chien wrote:> 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),da >ta=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--> Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK > +44 1225 386603 www.maths.bath.ac.uk/~sw283
Lung-Chang Chien
2008-May-21 15:14 UTC
[R] an unknown error message when using gamm function
The version is mgcv 1.3-27. Last night, I dropped out many arguments from my program to test which part might be the source of causing this error message, and found that the program works when just using corAR1(). This program has no this kind of problem in the beginning, but doesn't work after I impute some missing data in PM10 and tmean by multiple imputation approach. Hope this information can also help identify the element of this problem. LC> From: Simon Wood <s.wood@bath.ac.uk> > To: r-help@r-project.org > Date: Wed, 21 May 2008 08:55:02 +0100 > Subject: Re: [R] an unknown error message when using gamm function > What version of mgcv are you using, please?> On Wednesday 21 May 2008 06:16, Lung-Chang Chien wrote: > 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),da>ta=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]]