ltx5151 发表于 2013-12-29 10:58 
您好,您这样没有任何细节和具体内容,别的会员很难帮助您。建议将问题描述清楚,也更容易让大家了解情况, ...
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);这是我的程序