GLMER警告：方差-协方差矩阵[...]不是正定的或包含NA值。
有时候我发现使用lme4包中的glmer函数进行建模后，当调用其summary时，会出现以下警告信息： Warning messages: 1: In vcov.merMod(object, use.hessian = use.hessian) : variance-covariance...
在R中，lmer（或lme）中的随机效应变量是否自动被视为因子？
我了解在混合效应模型中，将连续或数值变量作为随机效应并不太合理（例如，请参见此处）。 但我想知道的是，R中的lme4::lmer或nlme::lme是否有意防止你这样做... 具体来说，我的问题是：如果我将任何非因子（非分类）变量作为随机效应提供给lmer（或lme），该函数会自动将其视为...
使用lme4的遗传算法时，glmuli运行时间不确定。
我正在使用glmuiti在R中进行模型平均。我的模型中有约10个变量，因此无法进行全面的筛选 - 因此我需要使用遗传算法（GA）（调用：method="g"）。 我需要包括随机效应，因此我使用glmuiti作为lme4的包装器。这里提供了执行此操作的方法：http://www.inside-...
如何从使用缩放响应拟合的lmer()模型中取消缩放系数
我使用lme4包中的lmer()函数在R中拟合了一个模型。我对因变量进行了缩放： mod <- lmer(scale(Y) ~ X + (X | Z), data = df, ...
拟合线性混合模型到一个非常大的数据集中
我想使用 lme4::lmer 在以下格式的 6000 万个观测值上运行混合模型；所有预测/依赖变量都是分类变量（因子），除了连续依赖变量 tc；patient 是随机拦截项的分组变量。我有 64 位 R 和 16Gb RAM，并且从中央服务器工作。RStudio 是最新的服务器版本。 mo...
在R中（使用lme4），将时间序列纳入混合效应模型
我已经搜索了类似的问题，但没有找到相关的问题，如果我错过了相关的问题，请谅解。 我正在研究在不同条件下饲料器使用时间（因变量），每个受试者访问30次饲料器。 受试者会接触一种类型的饲料器，该饲料器具有不同的组合，包括有香味/无香味、有视觉图案/无视觉图案，以及这些视觉或香味图案在两种空间安...
glmer的事后检验
我正在使用广义线性混合模型（glmer，lme4包）在R中分析我的二项式数据集。我想使用Tukey的事后检验（glht，multcomp包）对某个固定效应（“声音”）进行成对比较。 大部分都很好用，但我的其中一个固定效应变量（“SoundC”）根本没有方差（96次为“1”，0次为“0”），而...
一个混合效应模型的“显著”随机效应的毛毛虫图
我之前在这里得到了很好的帮助经验，现在又希望能得到一些帮助。 我正在估计一个相当大的混合效应模型，其中一个随机效应有超过150个不同的水平。这将使得标准的毛虫图变得难以阅读。 如果可能的话，我想要一个毛虫图，只显示随机效应中“显著”的水平（缺乏更好的术语）。也就是说：我需要一个毛虫图，其中...
如何预测merMod对象的项（lme4）？
对于简单的glm对象，我可以使用predict(fit, type = "terms")来检索每个术语的拟合值矩阵。 那么对于lmer或者glmer拟合模型，有相应的等效物吗？据我所见，predict.merMod函数不支持type=terms。
广义线性混合模型误差（二元响应）
我正在R中运行一个广义线性混合模型，用于二元响应变量，但是出现了错误信息。 我的代码如下： library('lme4') m1<-glmer(data=mydata, REPRODUCE~F1TREAT*SO+(1|LINE/MATERNAL_ID), family=binomia...