## Not run:
# library(modTempEff)
# data(dataDeathTemp)
# o1<-tempeff(dec1~day+factor(dweek)+factor(year)+factor(month)+
# csdl(mtemp,L=c(60,60),psi=20),
# data=dataDeathTemp, fcontrol = fit.control(display=TRUE))
#
# #add a ridge penalty: note how you *can* specify ridge!
# #you do NOT need to use csdl(..,ridge=..)
# o2<-update(o1, ridge=list(cold="l^2", heat="l^2"))
#
# #a model without temperature effects (the first drop.L obs are dropped)
# o3<-tempeff(dec1~day+factor(dweek)+factor(year)+factor(month),
# data=dataset,drop.L=60)
#
# #see ?anova.modTempEff for model comparisons
# ## End(Not run)
Run the code above in your browser using DataLab