我正在构建一个 GAM，其数据集的分布类似于泊松分布数据。 但是，我的数据是连续的，即它包含以立方米为单位的树木体积信息。 那么，在做R中的GAM代码（带mgcv库）时，我可以使用泊松作为族吗？ 或者我应该选择别的东西，因为数据不是计数数据？ 我确实发现了一些讨论类似问题的线程，但他们没有提供答案 ...
我有一个存在/不存在数据的大型数据集（~100k 观察值），我试图用具有共享惩罚的单个效应（例如 Pedersen 等人 2019 中的“S”）来拟合分层 GAM。 数据由作为数字的温度、作为因素的区域（5 组）组成。 这是我尝试安装的 model 的简单版本。 在第一种情况下，我收到以下错误：我假 ...
这个问题是从 Stack Overflow 迁移过来的，因为它可以在 Cross Validated 上回答。 5 天前迁移。 我可以获得带有预测变量的部分依赖图，也可以使用预测 function 进行 model 预测。 但是我想在研究期刊上发表 model，因此我需要将方程式格式的 mode ...
我想根据 x0 为 241 个概率值（'prob'）建模第 25、50 和 75 分位数回归曲线（q25、q50、q75）。 为此，我使用了 qgamV 包如下。 然而，这种方法导致一些 q25、q50、q75 值 <0 和 >1，这对于概率来说是不期望的。 从图形上看，人们会期望 q ...
我正在尝试使用 GAM {mgcv} 了解 NDVI 和海拔之间的关系。 首先，我强制将 NDVI 限制在 0 和 1 之间，以便我可以使用 beta 分布 然后我使用 mgcv 运行一个 gam 我试图绘制一个可解释的图，也就是说，y 轴具有 NDVI 值，这是我们根据模型预测“在地面上”所期望的 ...
我不确定我的查询是否适合这里。 我正在自学广义加性模型 (GAM)。 现在我的目标是模拟广义加性模型的时间序列数据。 但是我没有成功。 我试图从 GAM 书和 Simon N Wood 的软件包mgcv中找到 R ode，但找不到任何代码来模拟时间序列。 如果有人给我任何关于如何做到这一点的想法，我 ...
我有一个关于在带有 2 个线性预测变量的mgcv package 中使用交互项的问题。 我编写了代码以适应x1和x2之间的交互， mgcv::gam(y ~ x1 + x2 + ti(x1, x2, k = 3), method = "REML", data = DATA) 但我真的不确定这是在mg ...
这个问题是这篇文章的后续问题： Using gratia::draw() in R to display partial effect plots within an HGAM that are not relative to global smooth 我有一个看起来像这样的数据集： 我正在运行一 ...
我有一个简单的 GAM，其中一个很流畅。 当我在 model 上运行summary()或anova.gam()时，它为我提供了平滑显着性的 F 检验。 但是，尚不清楚使用什么 dfs 来计算 f 分布和 p 值。 我试图了解mgcv中的 mgcv package 如何确定 F 检验中平滑项重要性的 ...
我正在尝试从具有 GitHub 版本的 gratia 上的一些功能的 gam 中获取我的拟合值，并且在使用具有偏移量的 model 的data_slice() function 时遇到问题。 当我的 model 有一个偏移量时，我不确定我需要做什么来获取fitted_values() functi ...
我有一个如下所示的数据集： 我正在运行一个 HGAM（我认为它是 Pederson 等人的 GI model 2019 https://peerj.com/articles/6876/ ），如下所示： 我的部分效果图如下所示： 我对这些部分效果图的理解是，单个平滑 LakeTBend 与全局平滑 ...
我想通过使用gam结果制作带有ggplot2的轮廓 plot 。 以下是我想要的详细说明： 为了根据gam_fit预测 z 值，我从https://drmowinckels.io/blog/2019-11-16-plotting-gamm-interactions-with-ggplot2/中找到 ...
首先，我对 R 非常陌生，非常基本的统计知识，因此在我的分析中一直在努力。 这意味着搜索结果所需的编码，并且由于某些样本太小，我必须稍后检查它们是否具有任何统计相关性。 不过，就目前而言，我只是想达到在屏幕上显示图表的目标。 我有两个要运行游戏的数据集 - 一个有 9 个 obs。 22 个变量 ...
将by=参数添加到平滑和添加随机效果平滑有什么区别？ 我都试过了，得到了不同的结果。 例如： 产生： 另一方面： 产生： 我知道by=显示每个因素级别的摘要，但 model 的整体结果（例如 R^2）不应该相同吗？ ...
我经常在我的线性回归建模中使用指标函数，以便仅在次要协变量为 TRUE 时才可以估计系数。 这总是在协变量为 FALSE 时考虑影响没有生物学意义的情况下。 我最近想在 gam 框架中实现类似的东西，但是我通常使用的语法并没有扩展到 gam- 或特别是平滑器。 换句话说，是否可以拟合 model 以 ...
我正在使用vis.gam() ) 为名为 a1b 的广义添加剂 model (GAM) 制作透视图 plot，后者又利用persp中的透视图 function。 代码如下：library(mgcv) vis.gam(x = a1b, view = c("wbgt_max", "l ...
我正在尝试将平滑样条曲线拟合到看起来像具有两个峰值的数据。 首先，我将平滑样条曲线拟合到我的数据中，以识别结的潜在位置。 然后，我想使用 gam 函数进行回归，我可以在其中指定结的位置。 我收到一个错误， Error in smooth.construct.bs.smooth.spec(ob ...
我无法解决从mgcv运行bam()时遇到的错误。 我注意到 14 个月前在这里报告了一个类似的错误，似乎没有就解决方案达成一致 - 建议向 Simon Wood 发送电子邮件。 我的数据在这里。 数据集太大，无法粘贴dput()的输出。 如果我使用整个数据集运行以下模型，则会出现以下错误 ...
我有兴趣使用广义加法模型（游戏）估计在几个不同地点监测的计数随着时间的推移共享的全球趋势。 我已经阅读了Pederson 等人对分层游戏 (hgams) 的精彩介绍。 （2019） ，我相信我可以按如下方式设置模型（Pederson et al. (2019) GS 模型）， 我可以绘制局部效 ...
我使用下面的代码将平滑参数仅应用于交互项。 我的 model 的k值非常小，因为我只想获得非常平滑的变量样条曲线。 在此代码中，我尝试将sp = 1应用于交互项。 这是使用sp选项控制平滑度参数的正确方法吗？ 而且，我可以在ti()交互项中分别设置 X1、X2 的 sp 参数值吗？ （如下图） ...