R mgcv程序包中的广义加性混合模型(GAMM)分组变量的公式实现

文章/答案/技术大牛

发布

1回答

r、gam、mgcv、random-effects、multilevel-analysis

我正在尝试使用mgcv对一个以Time为主要自变量+2级协变量(Xc)的非线性多水平模型进行建模。因为这些数据是人内的，所以我希望该模型能够反映级别1 (Person_ID)的人内响应。问题是我无法理解mgcv包文档，因为我不知道如何为组编码(请参阅我代码中的random=~1|Person_ID，这是错误的)。以下是详细信息，数据和代码如下:因变量= DV，时间变量</e

浏览 204提问于2020-09-03得票数 1

回答已采纳

1回答

从r中类GAMM的随机效应中提取标准误差

r、random

建立了具有固定效应和随机截距效应的广义加性混合效应模型(即范畴变量)。在运行模型之后，我能够使用ranef(m1$lme)$x[[1]]提取每个类别的随机截取。我不知道这是不是因为这些函数只适用于类lmer的模型？有人能帮我从“游戏”类中提取出随机拦截的标准错误吗？我想使用随机截取和标准错误来绘制我的gamm模型的最佳线性无偏

浏览 2提问于2018-01-19得票数 7

1回答

在gamm4 R包中进行GAM-GEE？

r、statistics、spatial、gam

由于这些“跟随”之间的自相关性，我希望使用类似于Pirotta等人的GEE方法。2011年，使用包'yags‘和’样条‘()。他们的R脚本显示在这里()。重要的是要注意，模型有一个二元响应(生物体在横断面上的存在)，因此我认为gamm4比gamm更合适。在gamm中，它提供了下列因素内的自相关示例： ## more complicated autocor

浏览 7提问于2012-05-03得票数 2

1回答

型号选择(gamm4)疏导功能(MuMIn R包)出错:无法识别系列，已跳过型号

r、github、mumin

我正在尝试为广义的可加性混合模型(使用R中的包使用制作)进行模型选择。我本质上是在尝试遵循的一篇文献，用于使用MuMIn和gamm4进行模型选择。我正在创建一个有9个变量的模型，以及一个随机的个体效应-它看起来像这样：library(MuMIn) SouthFu

浏览 10提问于2015-03-17得票数 3

1回答

高斯GAM (MGCV包)因变量的方差估计

r、variance、gam

考虑下面的守则：X = runif(300, 0, 1)Y = X^3 + 2*X^2 + 1 + rnorm(300) 因此，model是一个高斯广义加性模型。如何求出因变量(Y)在model中的估计方差更新：在广义加性模

浏览 2提问于2020-04-12得票数 0

回答已采纳

1回答

如何在R中分析不规则时间序列

r、time-series、zoo

我在R中有一个zoo时间序列：50914, 50914, 50915892548000, 892591200, 892612800, 892634400我想执行ts包允许我执行的一些分析，例如将时间序列分解为趋势和季节性，并查看自相关函数。

浏览 1提问于2012-09-27得票数 14

回答已采纳

1回答

如何将gam()合并到来自格包的xyplot()中？

r、lattice、gam、mgcv

我试图将广义加性模型gam()从mgcv包合并到xyplot()函数或R中的lattice包中的coplot()函数。这是我的内核平滑代码。(x, y, span = .8, ...), xlab="Solar radiation (langleys)", ylab="Ozone (cube root p

浏览 2提问于2014-06-01得票数 3

回答已采纳

1回答

零充气模型(ziP)中的bam误差

r、poisson、gam、mgcv

我试图使用bam运行以下广义加性模型：但是得到以下错误： bam的文档提到了以

浏览 3提问于2017-11-07得票数 3

回答已采纳

1回答

扩展for-循环以适应R中的列表？

r、list、for-loop、matrix、nested-for-loop

我最近对开发一个for循环感兴趣，它可以运行多个广义加性模型，然后在一个表中根据模型中每一个光滑的p值、总体模型的偏差等对它们进行排序。我在堆栈溢出中发现了与相关的问题--这基本上是我想要的--并且能够很好地运行gam()而不是gamm()，但是我想将它扩展到包括模型中的多个独立变量，而不仅仅

浏览 1提问于2022-09-29得票数 0

回答已采纳

1回答

Lme错误：“reStruct中的错误”

r、mixed-models

这些是响应变量。我很难运行混合效果模型。所有蜂箱的响应变量都有一个数据框架。wifi“是一个与暴露时间相对应的二分法变量(1=on 0=off)，而治疗蜂箱是另一个暴露于wifi的蜂箱(1=exposure，0=control)的二分法变量

浏览 2提问于2014-11-29得票数 3

回答已采纳

1回答

有办法在mgcv的gam函数中包含自相关结构吗？

r、modeling、binary-data、gam、mgcv

我正在使用r中的mgcv包构建一个模型。这些数据有一系列的度量(数据是在15分钟的时间间隔扫描期间收集的，但不连续，例如，一天可能有5次连续扫描，然后直到第二天才进行任何扫描，等等)。当我拟合模型时，我的剩余自相关时间为2(即30分钟)，用ACF和PACF图进行评估。所以，我们想在模型中包含一个自相关结构，但我的理解是，只有gamm函数，而不是

浏览 17提问于2022-10-28得票数 2

回答已采纳

3回答

biglm预测无法分配大小为xx.xmb的向量

r、regression、linear-regression、lm

我有这样的代码：library(ff) testData <- read.csv- dependent ~ .predictedData <- predict(model, newdata=testData)无法分配大小为xx.xMB

浏览 4提问于2016-07-01得票数 2

1回答

如何从`mgcv::logLik.gam`输出中获取df

r、output、numeric、sentence

在R中，package mgcv具有logLik.gam，它给出了广义加性模型的对数似然值以及df。但是，我不能只取df值。但是要使用df的值(即3.612741)，有没有办法从输出中提取这个值呢？这里的输出是数字，实际上是一个句子。那么，有没有办法削减这个输出线，只取最后一个值呢？

浏览 7提问于2022-11-27得票数 0

1回答

从gamm、gam和lme模型中获取AIC或BIC柠檬酸:如何在mgcv中？我怎么能相信结果呢？

r、gam、mgcv

ML", na.action= na.omit, data= Pilot_fitbit2)警告消息: 1:在smooth.construct.tp.smooth.spec中(2:在smooth.construct.tp.smooth.spec中(object，dk$data，dk$knots)：，，我的主要问题是，我怎样才能从这里得到一个问题是，我想将这个值与gam模型的BIC值进行比较，并且我不确定哪个值(BIC(M.<em

浏览 5提问于2020-01-20得票数 2

回答已采纳

1回答

从包含矩阵"by“变量的mgcv::gam拟合中获取预测

r、mgcv

我刚刚发现，mgcv::s()允许为其by参数提供一个矩阵，允许对连续变量进行平滑处理，对变量组合中的每个变量(如果需要，还可以对它们的交互进行单独平滑)。然而，我很难从这样的模型中获得合理的预测，例如：library(ggplot2) #for plotting

浏览 2提问于2014-05-29得票数 0

1回答

我能用GAMM估计R中的时变季节效应吗？

r、time-series、gam、cyclic、stl-decomposition

我想用一个广义的加性模型来研究R中的时间序列数据，我的数据是每月的，我想估计一个季节效应和一个长期趋势效应。我关注了Gavin Simpson 和的一些有用的文章。我的数据如下所示：我试图用平滑的季节和趋势条件来指定一个广义的<em

浏览 8提问于2017-10-17得票数 3

回答已采纳

2回答

R图中轴标题和勾标之间的空格(使用mgcv包中的vis.gam )

r、plot、mgcv

我正在用a1b绘制广义加性模型(GAM)的透视图，使用vis.gam()，这反过来又利用了R中的persp函数。守则如下：vis.gam(x = a1b, plot.type = "perspylab = "lag (days)",

浏览 22提问于2022-07-06得票数 2

1回答

GAM (mgcv封装)预测的重现性

r、logistic-regression、prediction、gam、mgcv

我使用在bam() R包中使用mgcv函数建立的广义加性模型来预测二进制响应的概率。model <- bam(response ~ categorical_predictor1 + s(con

浏览 6提问于2021-10-04得票数 0

回答已采纳

1回答

更改默认plot.gam图的Y轴

r、mgcv

我已经使用mgcv包在R中运行了一个GAM，其形式如下： s(Mean.Front.density, bs = "cr"), data=r,family=gaussian) 在运行该模型并计算每个变量</

浏览 5提问于2013-04-02得票数 4

1回答

R中具有空间自相关的GAMM

r、spatial、gam、autocorrelation

据我所知，GAMM有一个随机误差或空间自相关误差结构。但我被GAM的模特儿弄糊涂了。下面是我的GAM (通用加性模型)代码。当我尝试使用GAMM (如下面提到的)时，我没有看到太多的变化： model <- gamm(Chave~s(Band_3) + s(Band_7)

浏览 0提问于2019-09-06得票数 1

回答已采纳

点击加载更多