最近在学R语言中的DLNM数据包，其中关于模型代码如下：
library(dlnm)
basis.o3 <- crossbasis(chicagoNMMAPS$o3, lag=10, argvar=list(fun="thr",
            side="h",thr=40.3),
            arglag=list(fun="strata",breaks=c(2,6)))##cen=?为默认
varknots<-equalknots(chicagoNMMAPS$temp,fun="bs",df=6,degree=3)
lagknots<-logknots(30,3)
basis.temp<-crossbasis(chicagoNMMAPS$temp,argvar=list(fun="bs",knots=varknots),
        cen=25,arglag=list(knots=lagknots),lag=30)
summary(basis.temp)
library(splines)
model<-glm(death~basis.temp+basis.o3+ns(time,7*14)+dow,family=quasipoisson(),
           chicagoNMMAPS)
pred.o3<-crosspred(basis.o3,model,cumul=TRUE,at=c(0:65,40.3,50.3))
pred.temp<-crosspred(basis.temp,model,cen=25,by=2)
pred.o3$allRRfit["50.3"]
cbind(pred.o3$allRRlow,pred.o3$allRRhigh)["50.3",]
但是在做模型的敏感性分析时，我想得到模型"model"中不同参数的df值的ACF值，那么他的代码怎么编写呢，用library(MASS)的stepAIC()等也不会，求大神指教。。