我们能否从R中的混合模型计算BLUP的标准误差(SE)？我正在使用一个名为ASReml的自定义包，它可以根据预测值计算SE。然而，我不确定在哪里可以找到BLUP值的SE。到目前为止，我有这样的东西...
df<-data.frame(inbred=c("x1","x2","x3","x4",&qu
浏览 6提问于2020-02-27得票数 2
在元分析包元中的R中，rma.uni()有一个blup()函数，但没有rma.mv()。我有一个带有随机截距项的rma.mv模型。如果我用：然后，我从固定效果版主那里得到了对原始数据集中的每个数据点的估计值，并且：给我预测我的随机效应的
浏览 6提问于2021-03-14得票数 1
回答已采纳
1回答
因此，当我认为某一特定特性组存在依赖时，就会使用混合效应回归模型。我附加了Wiki链接，因为它比我解释得更好。()在机器学习设置中，如何处理需要考虑混合效果的情况？有没有其他模型<
浏览 1提问于2017-09-08得票数 2
我用R中的来分析的社会调查。调查的指定(在第45页)权重已被缩放为平均值为1。当使用svydesign()函数时，我将权重变量传递给weight参数。在中，在surveysummary()函数下，它声明：
请注意，设计效果将是不正确的，如果权重已被重新标度，使它们不是抽样概率的倒数。因此，在使用诸如svyglm()、等函数时，会得到不正确的估计和/或标准错误吗？这引起了我的注意，因为当使
浏览 7提问于2021-08-14得票数 0
回答已采纳
我正在使用R包mgcv拟合GAM，以估计均值和SD： b <- gam(list(Y ~ s(X1)+s(X2)+s(X3), ~ s(X1)+s(X2)+s(X3)),relative_i
浏览 29提问于2020-12-04得票数 2
我的混合模型公式在SAS proc mixed语法中的随机效果如下所示：/type = vc subject = group solution; ((1|group) + (0 + color|group)) + ((1|group) + (0 + size|group))
这是正确的吗？我可以将sas
浏览 0提问于2020-11-22得票数 0
我正在使用lme4包进行线性混合效果建模fm01 <- lmer(sublat <- goal + (1|userid))此对象包括固定效果及其OLS标准误差(如下所示)goaltypeF4 -20.219 8.196 -2.467
goaltypeF5 -1
浏览 3提问于2014-10-17得票数 6
1回答
我正在进行一项元分析，其中我感兴趣的是X对年龄对生境使用的影响(原始平均值和方差)。我的一个模型的例子是： rma.mv( V = Used_variance, modsSlope, degrees = 2), data = vel.focal, )
我不使用引
浏览 6提问于2022-08-01得票数 0
2回答
我已经使用lmer在R中运行了一个线性混合模型。我正在尝试可视化随机效果结构。为了生成图形，我使用了print(dotplot(ranef(RT.model.4, condVar=T))[['part_no']])，其中part_no是来自混合模型的随机效果。它创建了类似这样的东西：这太棒了。然而，我希望能够在图形中直观地区分我的两
浏览 0提问于2020-06-05得票数 2
我希望将summary()的固定效果部分提取为data.frame。我使用lme4运行以下模型：summary(SleepStudy)SleepStudy_RE <- as.data.frame(VarCorr(SleepStudy))
对于固定效果是否有类似的代码行，包括估计
浏览 3提问于2016-05-20得票数 5
回答已采纳
我试着用随机效应考虑生理学上的个体差异(让我们暂时满足于拦截)，也许用另一种随机效应考虑子弹的疲劳。因此，我想在R中运行的模型是，在SPSS中： /FIXED=CONDhier CONDabund CONDhier*CONDabund, data=phys2)
但一个更复杂的模型超出了我的能力范围。我读
浏览 0提问于2013-11-21得票数 41
这可能是一个显而易见的问题，我还没有找到(R newbie)的答案，但是当使用lmer函数生成混合效果模型，然后使用以下命令显示结果时：我们得到了以恒星形式表示的beta估计值、标准误差和显著性水平的列表。使用什么测试来确定这些p值来标记恒星(重要的是，我认识到关于如何使用这些模型适当地确定显着影响存在一些争议
浏览 1提问于2017-05-12得票数 1
1回答
我想要拟合一个混合效果模型，并获得AIC排名的前十个候选模型的列表。对于固定效果模型，我使用glmulti:glmulti函数来拟合模型，并使用glmulti:weightable函数来获得AIC排序的候选模型列表，这将给出以下输出：## 8 outcome ~ 1 201.4210 0
浏览 26提问于2021-11-17得票数 1
1回答
我正在努力研究如何使用REML方法来适应复制交叉设计的模型。FDA建议的模型如上所述，有人可以帮助如何将其编码为R编码吗？这是我的代码，我想知道它是正确的还是错误的？
浏览 2提问于2020-12-21得票数 1
我已经成功地拟合了一个线性混合效果模型，我正在寻找提取单个组的随机效果分量。我知道可以使用以下命令提取随机效果的完整列表然后print(random.effects(模型))给出一个包含两列的组名称和随机效果列表，即使数据本身似乎只有一列我的问题是，是否可以通过组名“查找”与特定组相关联的</
浏览 2提问于2010-09-10得票数 2
回答已采纳
nlme包为我提供了一种使用resid(fitted对象，type=“规范化”)对归一化残差进行合成的方法，但是lme4没有选择这样做。如果没有lme4中的这个特性，我就无法诊断自相关。我不认为R stats包resid(lme4object,type="normalized")工作，lme4-object$residuals是不正确的语法。拟合混合效应模型描述的类"merMod“--混合效应模型
浏览 8提问于2022-08-23得票数 1
回答已采纳
因此，我在R中拟合了两个随机拦截的线性混合模型：其中b ~ MVN (0, Sigma)；X和Z分别是固定和随机效应模型矩阵，beta和b是固定效应参数和随机效应BLUP/条件模式。我想了解一下b的基本协方差矩阵，这在lme4包中似乎不是一件小事。您只能通过VarCorr获得方差，而不是实际的相关矩阵
浏览 2提问于2016-06-13得票数 6
回答已采纳
我将GLMM与spaMM package的fitme()-function一起使用。模型运行良好，但是没有发出p值(见下文)：p_v(h) (marginal L): -577.0237有没有人可以帮助我，如何获得这个模型的此外，您是否有一个想法，如何计算fitme(
浏览 0提问于2017-03-19得票数 2
1回答
我想知道有没有人能帮上忙 我需要建立一个基于以下数据结构的lmer模型： 3个小区，每个小区2个小区(给定2个不同的气候)，每个小区3个小区(给定3个不同的施肥量)，每个小区2个小区(有2个不同的品种我不明白“正常变量”和“混合变量”之间的区别。有人能给出缺失的代码吗： fit.cucumber.1 <- lmer(rate ~ ........., data = cucumber) 以下内容作为cucumber数据帧的一部分
浏览 22提问于2020-11-11得票数 1
1回答
R中采样纬向的纵向分析
、、、、
我有两个调查的纵向数据，我想做一个事后分析。通常，我会使用survey::svyglm()或svyVGAM::svy_vglm (对于多项式族)来包含抽样权值，但是这些函数并不考虑随机效应。另一方面，lme4::lmer占重复测量的比重，而不是抽样权重。那么，继续:你能告诉我如何在测试前分析2级或2级以上的分类结果吗？任何关于R中的包/函数以及如何使用/编写它们的提示都将不胜感激。，所以我在本文和其他文章中再看一看
浏览 3提问于2021-07-11得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
腾讯会议
云直播
对象存储活动推荐
腾讯云开发者
