mgcv_1.8-24:“fREML"或“REML"bam() 的方法给出了错误的解释偏差
使用“fREML"和“REML"方法将相同的模型与 bam 拟合给了我接近的结果，但是解释的偏差与 summary.gam 返回的不同. “fREML"的数量是~3.5%(不好)，而“REML"的数量是~50%(不是那么糟糕).怎么可能?哪一个是正确的? 很遗憾，我无法提供一个简单的可重现示例. #######################################方法=“f
..
mgcv相关内容
mgcv:如何在 predict.gam 中使用“排除"参数?
我有一个结构如下的模型，我想在忽略随机效应的同时提取预测值.如?predict.gam 和此处所述，我正在使用 exclude 参数，但出现错误.我的错误在哪里? dt 解决方案 exclude 与您假设的方式不同.您仍然需要在 newd 中为 predict.gam 提供所有变量.请参阅我的此答案，了解 predict.gam 背后的内容. 这是您需要做的: ## 用变量 c2
..
将多边形分配给 R 数据框中的数据点
我有两个数据框: points 包含一系列具有 x, y 坐标的点. poly 包含两个多边形的坐标(我实际上有 100 多个，但在这里保持简单). 我希望能够向数据帧 points 添加一个名为 Area 的额外列，其中包含点所在多边形的名称. poly 所以本质上在这个例子中，前 2 行应该有 Area="P1" 而最后一个点应该是空白的，因为该点不包含在任何多边形中.
..
stat_smooth gam与gam {mgcv}不同
我在ggplot2中使用stat_smooth函数，决定我想要“拟合优度"，并为此使用了mgvc gam.我想到我应该检查以确保它们是相同的模型(stat_smooth与mgvc gam)，因此我使用下面的代码进行检查.似乎它们有不同的结果，如该图所证明(图:stat_smoother gam(红色)，mgcv gam(黑色)).但是，我不知道为什么它们会有不同的结果.两者之间的某些默认参数是否不
..
Caret软件包-使用平滑和线性预测变量对GAM进行交叉验证
我想使用插入符号对GAM模型进行交叉验证。我的GAM模型具有一个二进制结果变量，一个纬度和经度坐标对的各向同性平滑，然后是线性预测变量。使用mgcv时的典型语法为： gam1
..
mgcv:如何返回估计的平滑参数?
考虑简单的GAM拟合如下: library(mgcv) my.gam
..
R中带有mgcv的粗薄板样条拟合(薄板样条插值)
背景 我正在尝试在 一个粗糙的薄板样条曲线拟合图2.3中的收入数据.这 拟合使训练数据错误为零. 到目前为止，我尝试了什么? 我试图复制先前的图2.5，它是平滑的薄板样条曲线拟合，不确定是否成功. income_2
..
具有"gp"的GAM更平滑:如何检索方差图参数?
我正在使用以下地理叠加模型 library(gamair) library(mgcv) data(mack) mack$log.net.area
..
尝试使用mgcv :: gam评估Markov随机字段时发生错误"nb/polys提供的区域名称与数据区域名称之间不匹配".
我尝试实施此很棒的博客加文·辛普森(Gavin Simpson)发布的文章使用通过cancensus软件包下载的数据，但在尝试评估gam时出现以下错误: Error in smooth.construct.mrf.smooth.spec(object, dk$data, dk$knots) : mismatch between nb/polys supplied area names
..
自定义链接功能适用于GLM，但不适用于mgcv GAM
很抱歉，答案很明显，但是我花了很多时间尝试在mgcv.gam中使用自定义链接功能 简而言之， 我想使用来自软件包 psyphy (我想使用 psyphy.probit_2asym ，我称为custom_link) 我可以使用此链接创建一个{stats} family对象，并将其用于glm的"family"参数中. m
..
为什么使用"mgcv :: s"?在"gam(y〜mgcv :: s ...)"中导致错误?
我想清楚一点，并在各行中使用::表示法来适合mgcv::gam.在模型调用mgcv::s中使用表示法时，我偶然发现了一件事情.带有可复制示例/错误的代码如下所示. 原因可能是因为我在模型公式中使用了这种表示法，但是我无法弄清楚为什么它不起作用/不允许这样做.这可能是有关语法的非常具体的内容(我想可能不是特定于mgcv的内容)，但是也许有人可以帮助我理解这一点以及我对R的理解. libr
..
R软件包mgcv是否存在已知的兼容性问题?是否有通用的兼容性规则?
我使用R版本2.15.1(2012-06-22)和mgcv版本1.7-22 我在R中加载了以下软件包: library(sqldf) library(timeDate) library(forecast) library(xts) library(tseries) library(MASS) library(mgcv) 碰巧我不能运行一个简单的模型(我省略了代码).甚至包括从帮
..
使用扫帚和tidyverse总结r平方的gams
我发布了一个问题 library(tidyverse) library(broom) iris %>% nest(-Species) %>% mutate(fit = map(data, ~mgcv::gam(Sepal.Width ~ s(Sepal.Length, bs = "cs"), data = .)), results = map(fit, glance),
..
如何获得smooth.spline的置信区间?
我已经使用smooth.spline来估计我的数据的三次样条.但是，当我使用公式计算90％的逐点置信区间时，结果似乎有些偏离.有人可以告诉我我做错了吗?我只是想知道是否有一个函数可以自动计算与smooth.spline函数相关的逐点间隔带. boneMaleSmooth = smooth.spline( bone[males,"age"], bone[males,"spnbmd"], cv=
..
从GAM光滑对象确定导数
我有一个非常简单的时间序列数据集，其中包含一个单变量("AVERAGE")的年平均值.我希望研究时间序列“趋势"分量的变化率(一阶导数)和加速度(二阶导数)以及相关的标准误差.我已经使用MGCV的GAM和PREDICT功能获得了“趋势"，如下所示: A
..
提取GAM估算值
我对R并不陌生，目前正在阅读一本书"Generalized Additive Models"，这是Wood by R的简介(2006年)，并进行了一些练习，特别是有关空气污染和死亡的部分，这是我的专长.兴趣.使用mgcv软件包，我运行以下模型. library(gamair) library(mgcv) data(chicago) ap1
..
如何仅从gam.check中获取地块
在mgcv包中应用gam.check时，R会生成一些残差图和基本尺寸输出.有没有办法仅生成图而不是打印输出? library(mgcv) set.seed(0) dat
..
为什么来自mgcv的bam对于某些数据比较慢?
我正在使用mgcv中的bam函数将相同的通用加性模型拟合到多个数据集.虽然对于我的大多数数据集，拟合过程都在10到20分钟之间的合理时间内完成.对于一些数据集，运行需要10多个小时才能完成.我发现慢速情况之间没有任何相似之处，最终的拟合度既非优劣，也未包含任何明显的离群值. 我如何弄清楚为什么在这些情况下拟合如此之慢?以及我如何能够加快这些速度? 我的模型包含两个平滑项(使用循环三次样
..
mgcv gam()错误:模型的系数比数据多
我正在为数据集使用GAM(广义加性模型).该数据集具有32个观测值，其中包含6个预测变量和一个响应变量(即幂). 我正在使用mgcv包的gam()函数来适应模型.每当我尝试拟合模型时，都会收到以下错误消息: Error in gam(formula.hh, data = data, na.action = na.exclude, : Model has more coefficien
..
如何为R中的每个变量提取GAM {mgcv}的拟合值?
我正在寻找一种添加模型中每个变量的预测值(实际的，不是标准化的)的方法 > model summary(m13) Family: gaussian Link function:
..