问题标签 [mgcv]

在使用 mgcv 包运行 gam 模型时，我遇到了一个我无法理解的奇怪错误消息：

“model.frame.default 中的错误（公式 = 死亡 ~ pm10 + Lag(resid1, 1) + ：可变长度不同（为 'Lag(resid1, 1)' 找到）”。

model1 中使用的观察次数与偏差残差的长度完全相同，因此我认为这个误差与数据大小或长度的差异无关。

我在网上找到了一个相当相关的错误消息here，但该帖子没有得到足够的答案，所以它对我的问题没有帮助。

可重现的示例和数据如下：

运行模型

获得偏差残差

为模型 1 添加一天滞后的偏差

这两种模型都产生了相同的错误消息。

我正在使用 GAMM 并希望使用方差-协方差矩阵提取模型的区间，但这样做时遇到问题

编码：

但最后一行代码返回，

Error in intervals.lme(gammdl$lme) : cannot get confidence intervals on var-cov components: Non-positive definite approximate variance-covariance

我不明白为什么会出现此错误消息。

我正在尝试复制 Simon Wood 在Generalized Additive Models: An Introduction with R using my data 第 316 页所做的事情。

我想比较在 spgwr 和 mgcv 之间产生的 GWR 配件，但我在 mgcv 的 gam 函数中遇到了错误。这是一个例子：

问题：

1. 如何使用 gam 和 bayesx 函数拟合相同的 gwr（位置的平滑函数）

2. 如何控制参数尽可能相似，包括最佳带宽

我有一个模型，y=f(x,z,a)。我想优化该模型（最终受到约束）。当具有梯度函数时，R 中的数值优化器要快得多。但是我已经非参数地拟合了我的模型，因此我不能轻易地分析得到梯度。是否有某种方法可以从拟合模型中获取梯度函数——类似于我想的拟合模型对象，并为其定义了预测方法？

这是一些虚拟代码：

定义变量：

拟合模型：

最小化以获得最小的 y：

我想要的是一个渐变对象，这样我就可以定义

然后运行

...考虑到大量数据以及复杂的模型和目标函数，这会快得多。

我想做的事可能吗？

这个

VGAM加载包时发生故障：

两者都实现s()了功能，但这不应该发生对吗？这是错误mgcv还是VGAM包装错误？

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

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

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

谢谢您的帮助。

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

我的模型是：

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

但它不起作用。

谢谢！！

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

这是我感兴趣的模型：

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

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

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

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

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

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

gam.check在mgcv包中应用时，R会产生一些残差图和基础维度输出。有没有办法只生成图而不是打印输出？