Я использую модель 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)