Проблема построения кривой реагирования xposure с использованием кросспред или перекрестного - PullRequest
0 голосов
/ 21 января 2020

Я использую модель GLM с 2 кросс-базисными матрицами. Я хочу построить кривую изменения экспозиции для конкретных c лагов (т.е. лаг = 5, лаг = 20, лаг = 25). Моя проблема в том, что я получаю разные графики для функций кросспред и кроссредуцирования. Я не уверен, что я понимаю их различия и что я делаю неправильно ..

Мой код:

model<-glm(AllCauseMort ~ Temp.basis + PM.basis + ns(AvgHum1, df=4) + dow + ns(time, df=98), 
          family=quasipoisson, data=subdata, na.action="na.exclude")
cen<-median(subdata$AvgTemp,na.rm=T)

pred <- crosspred(Temp.basis,model,by=0.5,cen=cen, lag=30, cumul=TRUE)   
plot(pred, lag=20, ylab="RR", xlab="Temp") 

rd<-crossreduce(Temp.basis,model,type="lag",lag=30,value=20,bylag=0.2,cen=cen)
plot(rd,xlab="Temp",ylab="RR",col=2,lwd=2)
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...