library(mgcv);### get package
dat<-read.csv("c:\\users\\Administrator\\Desktop\\heatwave.csv",head=T);### get data
dat$stratadow<-paste(dat$strata28,dat$dow,sep="");### create strata and day of week #Unadjusted library(tsModel); ttt2=mutate(dat,tm1=Lag(hw,1)); model=gam(death ~ tm1+as.factor(stratadow), family=poisson, matchdow=true,data=ttt2); summary(model);这是我的程序