在情况a中，mgcv R包中的gam代码工作正常。library(mgcv)fit <- gam(y~s(x0,bs = "cr", k = num_knots, m=2),data=dat)
summary(fit) 但是，当我在gam()中添加
浏览 43提问于2019-07-04得票数 0
回答已采纳
我有一个线性回归问题的预测。在这个问题上，一周中的几天很重要。通过这种方式，不同的拦截器(虚拟人前面的系数)给出了周相关性。但是，我想计算一周中每一天的x的不同系数。当我在样条函数中使用GAM (库: mgcv)时，我可以执行此操作，其中" day“是包含一周中的一天名称的一个分类变量。gam.mod
浏览 0提问于2019-02-02得票数 0
回答已采纳
我正在R中查看R包，我想知道如何根据新的数据更新模型。例如，假设我有以下数据，并且我有兴趣拟合一个三次回归样条。, lwd = 2, col = "red")我不清楚的是，当新的数据(即y)出现时，如何更新模型。mgcv::s()创建的任何内容，并在mgcv::gam的后续调用中重用它？然而，我仍
浏览 6提问于2022-03-06得票数 2
回答已采纳
1回答
当我试图为相关的错误结构拟合带有corARMA函数的模型时，我正在使用mgcv::gamm在GAMM上得到一个错误收敛。响应数据是降水量，并且是伽玛分布的.预测因素是年复一年(朱利安)。我正在将模型拟合到降水发生的观测的子集(即，我只对非0降水日进行建模)。你会发现数据中有很多NAs，这是因为我只对夏季的降水模式感兴趣。相关<e
浏览 5提问于2022-01-01得票数 1
回答已采纳
1回答
我使用来自gam包的mgcv来拟合负二项式族中的广义加法模型。我有一个包含因变量Y、自变量X、其他自变量Oth和因子Fac的数据框架。我想要符合以下模式每个因素级别都有不同的theta。换句话说，我用但这只为整个数据集提供了一个色散参数theta。相反，我相信不同因素<e
浏览 2提问于2015-12-07得票数 0
回答已采纳
我正在用R中的mgcv软件包来拟合GAM模型，我的一些预测器是圆形的，所以我使用的是一个周期平滑器。在交叉验证中遇到了一个问题，在交叉验证中，我的保留数据集可以包含超出培训数据范围的值。由于gam包自动选择打结作为平滑，这将导致一个错误(请参阅我的相关问题 -感谢@no葡萄和
浏览 1提问于2012-08-01得票数 0
回答已采纳
我进行了一项分析，需要绘制两个变量之间的非线性关系。我读过样条回归，其中一个挑战是找到节点的数目和位置。因此，我很高兴在中读到广义加性模型适合“具有自动选择节点的样条模型”。因此，我开始阅读如何在R中进行GAM分析，我惊讶地看到gam()函数有一个knots参数。
现在我很困惑。R中的gam()函数是否运行一个原子式查找最佳结的GAM？如果是的话
浏览 3提问于2020-02-05得票数 1
回答已采纳
1回答
我正在创建一个函数，它将未评估的对回归函数的调用作为输入，创建一些数据，然后对调用进行评估。函数是内部包的一部分，我想提供模型列表，然后在包中进行评估。
包含模型响应变量和公式所需协变量的数据
浏览 2提问于2018-12-12得票数 3
回答已采纳
最近，我在运行GAM模型时遇到了一个问题，该模型来自以前工作过的代码。我相信这与更新的R版本和mgcv包的更新版本有关。因此，如果有人有同样的问题或有解决方案，那就太好了。我目前运行的是:R version3.2.2 (2015-08-14) -- Windows上的"Fire Safety“。我使用的是mgcv包1.8-7。下面是一个示例代码，当在我<
浏览 2提问于2015-08-31得票数 0
我正在用GAM来模拟时间趋势的logistic回归。然而，我想从它中提取拟合样条，将它添加到另一种模型中，这在GAM或GAMM中是无法安装的。因此，我有两个问题：
我如何从合适的GAM中提取矩阵，以便我可以将它作为不同模型的估算？我正在
浏览 2提问于2013-03-23得票数 18
回答已采纳
让我列出我们的应用场景，以说明为什么只要我们想要使用全局匹配的参数协变量，这就很容易成为一个问题。下面的玩具数据代码(具有一维x中的空间趋势)适用于两个或三个偏移组。不幸的是，对
浏览 6提问于2011-06-16得票数 3
回答已采纳
我想比较一下spgwr和mgcv之间的GWR配件，但是我发现mgcv的gam函数有一个错误。下面是一个例子：require(mgcv)col.bw <- gwr.sel(crime ~ income<-gam(crime ~s(x,y)+s(x,y)*income+s(x,y)*housing, data=columbus)#mgc
浏览 1提问于2013-12-13得票数 0
回答已采纳
我尝试使用gam包来安装GAM (我知道mgcv更灵活，但这里需要使用gam )。我现在有一个问题，这个模型看起来很好，但是与原始数据相比，它似乎被一个常数沿y轴偏移，我不知道它是从哪里来的。col=1:2, legend=c("Original", "Minus intercept"))钱伯斯，J.M.和Hastie，T.J. (1993)在S (Chapman & Hall)中的<
浏览 2提问于2016-10-21得票数 2
回答已采纳
1回答
我正在研究一个模型，其中包括几个REs和一个变量的样条，所以我尝试使用gam()。但是，我达到了内存耗尽限制错误(即使在使用128 on的集群上运行时也是如此)。即使我只使用一个RE运行最简单的模型，也会发生这种情况。当我使用lmer()时，相同的模型(减去样条)运行平稳，只需几秒钟(或整个模型只需几分钟)。我想知道是否有人知道为什么gam()和lmer()之间的差异和任何潜在<e
浏览 2提问于2021-08-26得票数 5
一开始，我想要一个简单的函数来拟合曲线数据(我对数据一无所知)，并希望了解如何使用nls()或optim()来实现这一点。在我发现的类似问题中，似乎每个人都在暗示这一点。更新:stackoverflow上关于同一主题的一些有用的页面：
smooth.spline“开箱即用”在我的第一个和第三个例子中给出了很好的结果，但在第二个例子中很糟糕(它只是连接了几个点)。更新#2: gam() (来自<e
浏览 2提问于2011-09-26得票数 11
回答已采纳
我正在尝试将GAM模型拟合到由两对(x，y)值组成的数据集，即(x1，y1)和(x2，y2)，方法是首先拟合第一对，然后移动到第二对。当我在‘for’循环中调用gam函数时，它给出了一个错误“没有足够的(非NA)数据来做任何有意义的事情”。我怀疑这与我构造列的x1、y1、x2和y2标签的方式有关，因为在‘for’循环之外，gam函数可以工作。
谢谢!library
浏览 0提问于2020-12-09得票数 2
我想使用GAM函数在mgcv包中运行一个分级gam。我在brms中使用了相同的模型形式，没有问题，我最终会在brms中重新运行相同的模型，但是在周日是抽象提交的最后期限，所以我想尝试mgcv中的模型，以获得更快的结果。")) + s(INTERTIDAL_TRANSECT, bs = "re",
m
浏览 9提问于2022-01-28得票数 0
回答已采纳
3回答
在R中，可以使用包含转换(如mgcv或sqrt )的公式从GAM包中拟合GAM模型，默认情况下返回model.frame (只返回公式中指定的变量以及应用的转换)。示例：返回
> head(reg$model,3) loguntransforme
浏览 4提问于2017-03-18得票数 3
回答已采纳
1回答
在GAM (和GLM )中，我们正在拟合一个条件似然模型。因此，在对模型进行拟合后，对于一个新的输入x和响应y，我应该能够计算给定x的y特定值的预测概率或密度。例如，为了比较各种模型在验证数据上的适合性，我可能会这样做。是否有一种方便的方法在mgcv中安装了GAM？否则，我如何计算出所使用的密度的确切形式，以便适当地插入
浏览 2提问于2017-06-21得票数 1
回答已采纳
我的问题是关于不必要的预测器，即不提供任何新的线性信息的变量，或者是其他预测器的线性组合的变量。如您所见，swiss数据集有六个变量。它是Examination和Education的线性组合。ec <- swiss$Examination + swiss$Catholic
当我们运行一个带有不必要变量的线性回归时，R会删除其他项的线性组合，并返回NA作为它们的系数。Catho
浏览 1提问于2017-06-23得票数 5
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
腾讯会议
云直播
对象存储活动推荐
腾讯云开发者
