[R] please help for mgcv package

Simon Wood s.wood at bath.ac.uk
Thu Jun 23 12:58:14 CEST 2011


In the 'inst/scripts' directory of package 'gamair' there are files 
containing the R code for each chapter of the book (questionable layout 
uncorrected, but typo free, at least!).

Simon


On 21/06/11 16:27, pigpigmeow wrote:
> i read a book from WOOD, there's an example which is talking about the
> pollutant.
>
> library(gamair)
> library(mgcv)
> y<-gam(death~s(time,bs="cr",k=200)+s(pm10median,bs="cr")+s(so2median,bs="cr")+s(o3median,bs="cr")+s(tmpd,bs="cr"),data=chicago,family=Possion)
>
> lag.sum<-function(a,10,11)
> {n<-length(a)
> b<-rep(0,n-11)
> for(i in 0:(11-10))
> b<-b+a[(i+1):(n-11+i)]
>
> b}
>
> death<-chicago$death[4:5114]
> time<-chicago$time[4:5114]
> o3<-lag.sum(chicago$o3median,0,3)
> tmp<-lag.sum(chicago$tmpd,0,3)
> pm10<-lag.sum(log(chicago$pm10median+40),0,3)
> so2<-lag.sum(log(chicago$so2median+10),0,)
>
> I don't know what is the script (Bold font ) used for......
>
> and it shows "Error: unexpected numeric constant in "lag.sum<-function(a,10"
> ", why?
>
> anyone can answer me?
>
> --
> View this message in context: http://r.789695.n4.nabble.com/please-help-for-mgcv-package-tp3614485p3614485.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> 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.
>


-- 
Simon Wood, Mathematical Science, University of Bath BA2 7AY UK
+44 (0)1225 386603               http://people.bath.ac.uk/sw283



More information about the R-help mailing list