我使用R中mgcv包的gam函数将gam拟合到区间(0,1)上的数据。我的模型代码如下所示：
mod <- gam(y ~ x1 + x2 + s(latitude, longitude), faimly=betar(link='logit'), data = data)模型拟合得很好，但当我绘制<em
浏览 26提问于2016-09-01得票数 1
回答已采纳
1回答
在绘制GAM模型的预测拟合时，有没有办法设置x轴限制？更具体地说，我使用'by =‘为因子的每个级别拟合一个平滑器，但是，每个因子级别都有不同的值范围。在ggplot中绘制变量会导致x轴自动适应‘x’的不同范围；然而，在拟合GAM (mgcv::gam())之后，plot.gam()的默认行为似乎是预
浏览 38提问于2021-06-02得票数 1
回答已采纳
1回答
在GAM (和GLM )中，我们正在拟合一个条件似然模型。因此，在对模型进行拟合后，对于一个新的输入x和响应y，我应该能够计算给定x的y特定值的预测概率或密度。例如，为了比较各种模型在验证数据上的适合性，我可能会这样做。是否有一种方便的方法在mgcv中安装了GAM？否则，我如何计算出所使用的密
浏览 2提问于2017-06-21得票数 1
回答已采纳
我用mgcv::gam拟合负二项式模型，我注意到零偏差从fit到fit变化。如果我使用negbin家族函数而不是nb，那么问题就没有了。set.seed(3)dat <- gamSim(1, n=n)
## negative binomial data(paste("y ~ x0 +", pred
浏览 2提问于2017-10-09得票数 1
回答已采纳
我刚刚完成了安装使用mgcv包的GAM (我将称之为模型gam1.5)。我一直在玩vis.gam函数，我有一个问题是我无法解决的。我想标准化我的模型的拟合值，所以当我使用vis.gam时，z轴有限值0，1。我的想法是将规范化公式应用于GAM模型的$fitte
浏览 1提问于2017-06-22得票数 2
回答已采纳
建立了具有固定效应和随机截距效应的广义加性混合效应模型(即范畴变量)。在运行模型之后，我能够使用ranef(m1$lme)$x[[1]]提取每个类别的随机截取。但是，当我试图使用se.ranef(m1$lme)提取随机效应的标准错误时，该函数不起作用。其他使用se.ranef(m1)和se.ranef(m1$gam)的尝试也不起作用。我不知道这是不是因为这些函数只适用于类lmer的
浏览 2提问于2018-01-19得票数 7
我正在用GAM来模拟时间趋势的logistic回归。然而，我想从它中提取拟合样条，将它添加到另一种模型中，这在GAM或GAMM中是无法安装的。因此，我有两个问题：
我如何从合适的GAM中提取矩阵，以便我可以将它作为不同模型的</
浏览 2提问于2013-03-23得票数 18
回答已采纳
我试图同时对两个约束条件下的数据拟合一个GAM模型：(1)拟合是单调的(递增)，(2)拟合通过一个不动点，比如(x0,y0)。
到目前为止，我设法使这两个约束分别工作。对于(1)，在的基础上，利用mgcv::mono.con()得到了满足单调性的线性约束，并利用mgcv::pcls()对模型协方差进行了估计。对
浏览 1提问于2021-03-16得票数 10
1回答
我使用来自gam包的mgcv来拟合负二项式族中的广义加法模型。我有一个包含因变量y、自变量x、因子fac和随机变量ran的数据框架。我符合以下模型我在负二项式回归书中读到，模型</e
浏览 0提问于2019-12-15得票数 3
回答已采纳
我用geom_line和geom_ribbon (图1)创建了一个情节，结果是可以的，但为了美观起见，我希望线条和丝带更加流畅。我知道我可以用geom_smooth来表示线条(图2)，但我不确定是否可以平滑核糖核酸，我可以为丝带的顶部和底部创建一条geom_smooth线(图3)，但这两条线之间是否存在填充空间呢？
浏览 6提问于2020-04-03得票数 1
回答已采纳
2回答
根据通气量的大小和通气量的大小，收集通气量数据。采集了几个样本，并将其拟合为三次方程。用Excel编写，得到了第三个回归方程。 ? 然而，从图中可以看到，90-95%的通风量高于100%。数据不应高于100%，但自回归的最大顶点是凸的，因此它以曲线的形式超过100%。 有没有办法减少最大顶点并对其进行拟合？按原样使用测量数据，但不要超过100%。 也欢迎使用R或其他统计程序。R值可以稍微低一点。 谢谢。
浏览 49提问于2019-09-27得票数 0
回答已采纳
3回答
几年前的描述了如何提取用于绘制拟合gam模型的平滑分量的数据。它是有效的，但只有当有一个平滑的变量时。我有多个平滑变量，不幸的是，我只能从系列的最后一个变量中提取平滑。下面是一个示例：a = rnorm(100)y = a*b/(a+b)
summary(mod)
浏览 3提问于2013-04-06得票数 16
回答已采纳
由于我的数据具有二元响应，但很少发生事件，因此我希望通过拟合bgeva模型而不是gam模型来改进其预测。为了证明和比较它的预测准确性，并将其与我尝试过的其他模型进行比较，我需要计算AUC并绘制ROC曲线。
问题是，我的代码可以与glm和gam一起工作，但不能与bgeva对象一起工作。下面是bgeva()包中的一个简单示例，以及我用来计算glm和gam对象的
浏览 3提问于2015-08-30得票数 5
1回答
我和一位统计学家谈过，问他，我怎样才能识别和删除那些对gam (广义加法模型)有太多影响的点。library(mgcv)
set.seed(11)
浏览 0提问于2017-05-04得票数 1
2回答
我正在尝试对比例数据(介于0和1之间的数字)运行GAM。m5 <- gam(prop ~ s(SST_mean) + s(NAO, bs="re") + WarmCold + Cycle6 + Region,
family=binomial, data=DAT_WC, method = "REML") prop =发生的比例，WarmCold =物种是热适应的还是冷适应的，Cycle6
浏览 158提问于2020-04-17得票数 0
回答已采纳
我使用mgcv包对一些数据进行多项式样条拟合，方法是： family =poisson(link = "log"))
我在试着提取合适的功能形式。x.gam是一个gamObject，我一直在阅读文档，但还没有找到足够的信息来手动重构合适的函数。x.gam$smo
浏览 2提问于2016-05-22得票数 13
回答已采纳
1回答
mgcv后部覆盖
、、、
我正在用R中的mgcv包拟合一个空间二项式模型，并想要模拟预测点的后验分布(代码如下)。我一直在使用模拟数据来测试后验的覆盖特性。我发现，当总流行率在0.5左右(50%)时，覆盖率很低(大约35%的真实值在95%的后间隔内)，但当你离开0.5时，这种情况会有所改善。例如，当平均患病率为1%时，~97%在95%的后部。我想我的问题是：
，这是在这种方法中使用GAMs
浏览 1提问于2019-12-05得票数 0
2回答
任务：在移动窗口中找到最佳线性拟合的斜率(例如，最小化误差方差)。X值是等距的，例如随着时间的推移自动测量。朴素实现:y值的循环。#some dataset.seed(42)
plot(y~x,pch=optWinLinFit0
浏览 11提问于2013-06-25得票数 5
1回答
R软件中GAM模型的故障
、、、、
我尝试在R上运行以下代码：我的数据集是这样的：
"Error in if (abs(old.score - score) > score.scale * conv.tolIn addition: There were 50
浏览 34提问于2021-11-01得票数 0
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储
腾讯云开发者
