跪求大神指导R软件DLNM结果分析，3D图分析

model1<-gam(case~cb1.Tem+ns(time,7*8)+as.factor(dow)+as.factor(doy),
family=quasipoisson,data=zz)
pre1.Tem<-crosspred(cb1.Tem,model1,cen=19,by=1)
plot(pre1.Tem,xlab="Tem",zlab="RR",theta=200,phi=40,lphi=30,
main="3D Graph of Tem effect")
R软件中的DLNM模型，使用GAM，最后得到的3D图如下，跪求大神指导怎么查看RR值最高时所对应的滞后天数
以及Tem的值。谢谢

您好，我的程序和您写得差不多，但是我总是报错，请问您知道为什么吗？
model1 <- glm(num ~  cb.tavg + cb.rh  +factor(dow) , family=quasipoisson(),data1)
summary(model1)
pred1.tavg <- crosspred(cb.tavg,model1)
plot(pred1.tavg, xlab="Temperature", zlab="RR", theta=200, phi=40, lphi=30, main="3D graph of temperature effect")