关于DLNM数据包的ACF代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

最近在学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()等也不会，求大神指教。。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：dlnm ACF 数据包 Library BRARY

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)))
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",]

就没人回复吗

ACluo 发表于 2017-12-10 10:24
最近在学R语言中的DLNM数据包，其中关于模型代码如下：
library(dlnm)
basis.o3

能请教一下您关于DLNM模型的相关问题吗，非常感谢， DLNM模型中的fun是怎么选择的呢