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@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.
