问题标签 [gam]

在r中的mgcv包中使用gam时，有什么方法可以识别结的数量和位置吗？所以我有这样的事情：

我想找到打结的地方，这样我就可以在这里使用它们

并节省大量时间试错，确定多少节以及它们应该在哪里提供最佳拟合。

谢谢您的帮助。

我已经安装了一个COZIGAM模型，它类似于 GAM 模型，但用于零膨胀数据。

我的模型是：

如何s(lon1,lat1,sst1)用 3d 图表示第二条样条 ( )？我试过这个：

但它不起作用。

谢谢！！

我vis.gam()用来绘制模型的结果。代码如下：

我想从我自己的调色板中设置新颜色，如以下代码中所述：

不幸的是，我有这个错误：

vis.gam(model, view = c("CHLA", "SST"), xlim = c(0, 15), ylim = c(10, : 无法识别配色方案) 中的错误

我不知道如何使用该jet.colors功能的调色板vis.gam...

谢谢你的帮助 ！

我正在尝试运行 GAM 模型并根据此线程对包 mgcv 的建议以给定值获得精确的估计值。使用 predict 函数时，结果与 GAM 的估计值不匹配。

这是我感兴趣的模型：

当我使用以下方法获得结果时：

结果与我尝试使用以下方法估计变量 a 上的单个绘图点时不同：

有人对我可能做错的事情有任何建议吗？

分离包R不是很好的做法（请参阅参考资料?detach），但由于某些原因，我必须在包gam和mgcv. 一旦mgcv被附加和分离（并且命名空间中的所有依赖项都被卸载！），函数会gam产生一些奇怪的错误（请原谅术语）。似乎 - 即使在前一步卸载 -mgcv并且朋友又回到了命名空间并且函数调度出错了。以前有人遇到过同样的问题吗？

我的会话信息是（新会话）：

我使用最新版本的 gam (1.09) 和 mgcv (1.7-28)。任何提示表示赞赏！

我正在做物种分布建模（生态位模型），我将模型投影到当前或未来的气候栅格（Bioclim 变量）。

当我预测当前栅格（object=bio_stack）时，一切正常，使用以下代码：

但是，当我预测未来的栅格 ( object=future_bio_stack) 时，此代码不会产生结果：

相反，我收到以下警告消息：

我投影到的未来生物堆栈层在图像中看起来不错，具有正常的最小值和最大值（如果没有信誉分数，我无法发布图像）。

然而，future_bio_stack 的最小值和最大值在摘要中看起来很奇怪，尤其是与 bio_stack 相比。以下是每个堆栈的摘要，并提供了最小值和最大值：

任何想法为什么bio_stack正在使用该predict功能但future_bio_stack不是？

我仍在 R 的学习曲线中（并且在这个论坛中是新的），并且一直在尝试循环一个 gam 模型，在该模型中我试图评估对环境预测变量的不同响应（因变量 - 在这种情况下是鱼类）。所以基本上我想在一个循环中重新创建以下内容：

g1<-gam(var1~s(x1)+s(x2))

g2<-gam(var2~s(x1)+s(x2))

g3<-gam(var2~s(x1)+s(x2))

我的数据框由单独列中每个物种的丰度组成，然后是单独列中的环境预测因子。我遵循了这个论坛中的一些建议，使用 melt 和 dlply 按物种迭代 gam 模型。我使用了以下代码：

.上面的代码创建了一个data.frame，每个物种观察一行

代码运行没有错误，然而，当我检查结果时，我注意到不同物种的结果是相似的（即，相同的调整 R、变量 p 值、GCV 分数等）。有人知道如何运行多个游戏吗？或找出错误的来源？

欢迎和赞赏所有想法！希望我没有含糊其辞。如果您需要任何进一步的信息，请告诉我。

最好的，罗洛

我能够在 r 中绘制一个广义线性模型。

这是我所做的：

这是错误：

尝试过 - 重新启动 r - 使用 mgcv 包而不是 gam

不久前我能够绘制它，但我一定做了一些阻止我绘制的事情。

我正在gam使用以下方法进行绘图：

我想根据变量的值在同一个图上有不同的线TargetGroup2012。

我可以这样做：

有没有更快的方法来做到这一点？也许使用类似 groupby 的东西？

编辑：正确的解决方案

当我尝试使用 ggplot2 时收到错误消息，但我不明白为什么。我可能在某个地方犯了一些小错误......

我的代码有什么问题？