1

我正在使用 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

非常感谢,

安东尼奥·佩德罗。

4

1 回答 1

0

查看代码,图形范例似乎是基本图形(特别是persp函数),并且有尾随点参数,所以为什么不添加:

..., `xlim=range(columbia1(maternal_age), ylim=range(columbia1$birth_year)`
于 2013-02-11T04:47:29.150 回答