我在R中使用mgcv包运行GAM,表格如下:
shark.gamFINAL <- gam(ln.raw.CPUE...0.1 ~ Year + Month +
s(Mean.Temp, bs = "cr") + s(Mean.Chl.a, bs = "cr") +
s(Mean.Front.density, bs = "cr"), data=r, family=gaussian)
在运行该模型并计算每个变量解释方差百分比后,我想绘制每个变量对响应的影响图。然而,在使用R中的plot.gam函数时,我的图表出现了一个y轴,其标签为“s(预测变量,edf)”。我不确定y轴刻度代表什么意思?是否有办法将y轴范围更改为代表响应的范围,就像在Walsh和Kleiber(2001)的论文中所做的那样:'Generalized additive model and regression tree analyses of blue shark (Prionace glauca) catch rates by the Hawaii-based commercial longline fishery' 的第三幅图中一样。我本想上传一些我描述的图表示例,但由于这是我的第一篇文章,我没有至少10个声望,因此无法这样做!我已经搜索了许多网站和论坛,试图找到答案,但没有成功,因此任何帮助都将非常感激!
termplot
方法来绘制gam
对象的图形。Year
和Month
在模型中可能是数字或因子 - 从您提供的信息中不清楚,但它们都是参数项。您需要根据每个月份和年份组合创建预测数据,如 PDF 中所述,然后绘制图形。我认为现在您需要一个新问题,并提供一个可重现的示例,这样我们才能更好地帮助您 - 评论真的只适合这样的扩展跟进。 - Gavin Simpson