1回答
我想在R中“可交换”类型的边际模型中估计一些协变量N的平滑效果，其中聚类变量是S。根据我所能找到的，这应该是可能的： geeglm(..., id = S, corstr = "exchangeable") 以及： gamm(..., correlation = corCompSymm然而，如果我使用我们的项目所涉及的真实数据，估计的平滑效果往往会有很大不同。我不能在这里发布，但
浏览 37提问于2019-09-23得票数 0
我用各种函数来比较我的离散计数数据。我在R中用geepack拟合了GEE模型，在log(count)上拟合了线性混合效应模型lme (nlme)，用glmer (lme4)拟合了GLMM，用gamm4 (<code>E 110</code>gamm4</code>E 211</code>)拟合</
浏览 3提问于2012-02-08得票数 4
回答已采纳
1回答
该模型具有二项响应、日随机效应、固定效应和三种平滑效应。我的理解是，REML是二项式模型的最佳拟合方法，但是这种方法不能使用二项式模型的gamm函数来指定。因此，我正在使用gam函数，以允许使用REML拟合。当我拟合模型时，我的剩余自相关时间为2(即30分钟)，用ACF和PACF图进行评估。所以，我们想在模型中包含一个自相关结构，但我的<
浏览 17提问于2022-10-28得票数 2
回答已采纳
我正在一个包含13,500个观察值(这里是学生)的数据集上拟合一个gee模型。学生被分成52所不同的学校。我知道有证据表明学生在学校内嵌套(低ICC)，因此我应该在方差协方差矩阵中调整这种嵌套效应。我计划做的是首先拟合一个具有可交换var-cov结构的gee模型。然后，在此基础上，我将运行Huber-White三明治估计器，也称为鲁棒方差估计器。我为鲁棒方差估计器编写了自己的代码，
浏览 0提问于2012-07-13得票数 4
我有一个类似于下面的数据，我想用对数链接拟合一个具有正态方差分布的可交换相关矩阵。replace = T), gee(events ~ sport + sex, id=id, data = x,
corstr = "exch
浏览 2提问于2020-03-05得票数 0
建立了具有固定效应和随机截距效应的广义加性混合效应模型(即范畴变量)。在运行模型之后，我能够使用ranef(m1$lme)$x[[1]]提取每个类别的随机截取。但是，当我试图使用se.ranef(m1$lme)提取随机效应的标准错误时，该函数不起作用。其他使用se.ranef(m1)和se.ranef(m1$gam)的尝试也不起作用。我不知道这是不是因为这些函数只适用于类lmer的模型？
有人能帮我从“游戏”类中提取出随机拦截的标准错误吗？我想使用随机
浏览 2提问于2018-01-19得票数 7
我希望能够用可交换的var-cov矩阵拟合gee模型，然后在结果模型上运行Huber-White三明治估计器，以防止偏差结果。我的GEE模型的代码如下：class SSID SCHIID0809 Ethnicity(ref = "500") ELLbaseline有没有什么办法，我可以通过一个宏来传递结果，这个宏基于从上面的GENM
浏览 8提问于2012-07-13得票数 0
回答已采纳
我安装了两种型号，一种是gam，另一种是gamm。因此，唯一的区别是分布假设。我想比较这两种模型，但我猜AIC不会起作用，因为这两种模型基于不同的方法？那么有没有其他合适的估计可以用来进行比较呢？
我知道我可以用gam来拟合第二个问题，但是为了这个问题，让我们忽略它。
浏览 2提问于2017-06-04得票数 0
1回答
我正在研究一个模型，其中包括几个REs和一个变量的样条，所以我尝试使用gam()。但是，我达到了内存耗尽限制错误(即使在使用128 on的集群上运行时也是如此)。即使我只使用一个RE运行最简单的模型，也会发生这种情况。当我使用lmer()时，相同的模型(减去样条)运行平稳，只需几秒钟(或整个模型只需几分钟)。我想知道是否有人知道为什么gam()和lmer()之间的差异和任何潜在的解决方案。下面是一些具有模拟数据和最简
浏览 2提问于2021-08-26得票数 5
我试图在一个具有固定效果的GAMM模型中指定随机截距和随机斜率项。我已经使用mgcv库中的以下代码成功地拟合了一个具有随机截距的模型，但现在无法确定gamm()函数中随机斜率的语法是什么：如果我在线性混合效果模型中同时使用随机截
浏览 28提问于2014-09-16得票数 5
我正在尝试对广义估计方程(GEE)执行对数回归。Politician_ID,repeated=Country_ID)The response variable should have more than 2 categories10 3945 1 0 0.549 -1.30
由于我是回归方法论的新手
浏览 1提问于2018-08-07得票数 0
回答已采纳
1回答
当我试图为相关的错误结构拟合带有corARMA函数的模型时，我正在使用mgcv::gamm在GAMM上得到一个错误收敛。响应数据是降水量，并且是伽玛分布的.预测因素是年复一年(朱利安)。我正在将模型拟合到降水发生的观测的子集(即，我只对非0降水日进行建模)。你会发现数据中有很多NAs，这是因为我只对夏季的降水模式感兴趣。一位同事建议我用NAs代替其他日期(有人告诉
浏览 5提问于2022-01-01得票数 1
回答已采纳
我使用mgcv运行了一个包含大型数据集(超过20,000个案例)的GAMM模型。由于数据点的数量很多，很难在图中的残差点之间看到平滑的线条。是否可以为点和平滑拟合线指定不同的颜色？下面是一个取自mgcv文档的示例：set.seed(0)
dat<- gamSim(1,n=20
浏览 4提问于2015-08-07得票数 1
我正在尝试利用来自模型的拟合值，并需要它们与我正在处理的数据帧中的正确行相匹配。下面是我运行的模型： random = ~ (1|school),对于lmer对象，&q
浏览 1提问于2012-11-19得票数 1
回答已采纳
在基本环境中使用stats::glm估计的glm对象具有一个方便的属性，用于因子的编码级别，甚至是根据效果的字母顺序编码为因子的字符变量。该属性是使用不可见函数.getXlevels计算的xlevels槽，该函数使用模型术语和模型框架。
模型框架是复杂的大型对象，回想起来很难计算。我正在使用gee软件包来获得使用可交换相关结构估计的模型系数<em
浏览 0提问于2014-08-21得票数 0
1 26.5 40.29999 28.500000 4我正在尝试拟合一个GEE模型来预测数据帧每一行的Weight。当我尝试使用mutate和predict函数执行此操作时，我得到一个错误，指出模型拟合中使用的观测值数量(789)与原始数据帧中的观测值数量(861)不同。: 1.如何提取模型fit中使用
浏览 8提问于2019-08-17得票数 0
回答已采纳
9001695 28.6 5 NA 4 6.5 9002316 25.1 5 89.3 6 library (gee)
model1<- gee(A~time+BMI, id=ID, corstr= "independence", data = long)Number of Iterations: 1
浏览 5提问于2014-11-14得票数 0
使用lme4软件包提取随机效果和随机效果的方差-协方差矩阵如下：fm1 <- lmer(Reaction ~ Days + (1|Subject), sleepstudy)gam.vcomp‘函数确实提取了估计的方差分量，但不是每个随机效应级别的方差分量。
浏览 5提问于2018-01-10得票数 1
回答已采纳
N个蜂箱= 4，n个传感器在每个蜂箱中收集的数据=~百万。所有蜂箱的响应变量都有一个数据框架。wifi“是一个与暴露时间相对应的二分法变量(1=on 0=off)，而治疗蜂箱是另一个暴露于wifi的蜂箱(1=exposure，0=control)的二分法变量。Object must be a
浏览 2提问于2014-11-29得票数 3
回答已采纳
这可能是一个愚蠢的问题，但经过几个小时的寻找，我仍然没有找到解决方案。我已经运行了一个gamm模型，该模型似乎工作得很好。现在，我想使用在模型中计算的对数似然来计算模型的AIC (mod1)。模型运行如下所示
mod1<-gamm(Biomass..g. ~ Atoll + SizeCm + Depth + OuterOrLagoon,random=list(SiteID=~1),family
浏览 1提问于2015-01-05得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
腾讯会议
云直播
对象存储活动推荐
腾讯云开发者
