Re: [R] an unknown error message when using gamm function

From: Lung-Chang Chien <chienorz_at_gmail.com>
Date: Wed, 21 May 2008 11:14:20 -0400


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_at_bath.ac.uk>
> To: r-help_at_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]]

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Wed 21 May 2008 - 17:35:49 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 23 May 2008 - 11:30:40 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive