我使用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
中变量的图中预测值范围的任何建议?
非常感谢,
安东尼奥佩德罗。答案 0 :(得分:0)
看看代码看起来图形范例是基本图形(特别是persp
函数)并且有尾点参数,所以为什么不添加:
..., `xlim=range(columbia1(maternal_age), ylim=range(columbia1$birth_year)`