如何实现分类变量引入GAM模型

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我做PM10对心脑血管疾病门诊人数的影响，将PM10根据浓度大小分成四类，0-0.1为第一类，0.1-0.15为第二类，0.15=0.25第三类，大于0.25为第四类，我现在要求分别求出第一类，第二类，第三类，第四类分别对应得RR值，请问我如何实现呢
脚本如下
rm(dow)
rm(time)
rm(wind)
rm(press)
rm(vi)
rm(atemp)
rm(htemp)
rm(dtemp)
read.table("C:/Users/14916/Desktop/jieguo/XZshuju/XINZONGLF.txt",header=T,fill=T)->fitty
fix(fitty)
attach(fitty)
downum=factor(dow)
downum
library(mgcv)
fitty=gam(low~downum,family=poisson,na.action=na.omit)
summary(fitty)
AIC(fitty)
stat <-NULL
mod<-list()
for (i in 1:50){
mod [[i]]=fit=gam(low~downum+s(time,k=11,fx=T)+s(wind,k=3,fx=T)+s(press,k=3,fx=T)+s(dew,k=3,fx=T)+s(vi,k=11,fx=T)+s(htemp,k=i,fx=T),family=poisson,na.action=na.omit)
tt=AIC(mod[[i]])
stat<-append(stat,tt)
}
plot(x=(1:50),y=stat)
stat
summary(stat)
fitty=gam(low~downum+s(time,k=11,fx=T)+s(wind,k=3,fx=T)+s(press,k=3,fx=T)+s(dew,k=3,fx=T)+s(vi,k=11,fx=T)+s(atemp,k=7,fx=T)+PM10,family=poisson,na.action=na.omit)
summary(fitty)
AIC(fitty)

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：心脑血管疾病 脑血管疾病 如何实现 心脑血管 血管疾病