我正在使用 vis.gam() 绘制以下模型的结果:
gam.5 <- gam(mortality.under.2~s(maternal_age)+ s(birth_year) +
te(birth_year,maternal_age) + wealth + sex +
residence+ maternal_educ
,data=colombia1,family="binomial")
制作情节的代码如下:
vis.gam(gam.5,view=c("maternal_age","birth_year"),
plot.type="contour",type="response",
cond=list(wealth="Lowest quintile",sex="Female",
residence="Rural",maternal_educ="Primary"))
不幸的是,它产生的预测值超出了变量“maternal_age”和“birth_year”的数据范围。关于如何限制图中变量的预测值范围的任何建议view
?
非常感谢,
安东尼奥·佩德罗。