在R中加速lmer函数
我想分享一些关于如何使用lme4软件包在R中改进线性混合效应模型拟合时间的想法。 数据集大小：该数据集大约有400,000行和32列。不幸的是，无法透露数据的性质。 假设和检查：假定响应变量来自正态分布。在模型拟合过程之前，使用相关表和alias函数检查变量是否存在共线性和多重共线性。 ...
lmer模型拟合的R平方
我有一个混合效应模型，想要查看R²值和p值。我以为可以通过summary()访问，但实际上不行。或者至少我没有找到。 > summary(fit1.lme <- lmer(log(log(Amplification)) ~ poly(Voltage, 3) + (1 | Seri...
尝试在R中运行glmer时出现警告信息。
我正在尝试在最新版本的R和lme4上重新运行一份旧的数据分析报告，其中包括二项式glmer模型（来自2013年初），因为我已经不再有旧版本的R和lme4。然而，我遇到了与dmartin和carine以前的帖子中类似的警告信息（第一个警告信息），以及其他stack overflow之外的帖子中的...
lme4中的多元线性混合模型
我想知道如何使用lme4拟合多元线性混合模型。我已经使用以下代码来拟合单变量线性混合模型:library(lme4) lmer.m1 <- lmer(Y1~A*B+(1|Block)+(1|Block:A), data=Data) summary(lmer.m1) anova(lmer....
增加新版本lmer的迭代次数？
我刚刚更新了lme4到1.0-4版本，当我运行我的混合效应模型lmer()时，它以前收敛但现在打印出这个警告：Warning message: In (function (fn, par, lower = rep.int(-Inf, n), upper = rep.int(Inf, : ...
PLM或LME4适用于面板数据的随机效应和固定效应模型。
我能否使用lme4在面板数据上指定随机效应和固定效应模型？ 我正在使用r重新执行Wooldridge（2013年，第494-5页）的示例14.4。感谢this site和this blog post，我已经成功在plm包中完成了它，但我想知道是否可以在lme4包中完成相同的操作？ 这是我在...
lme4错误：边界（奇异）拟合：参见？isSingular
我正尝试在R中运行lme4包。 我总共有10条线，每条线上都有两个重复的实验，每个实验中有四棵植物。 但是一些植物死了，有一些缺失值。Weight是响应变量。以下是数据的一些行：Line Rep Weight PLANT Line 1 1 NA 1 Line 1 1 NA...
使用lme4进行新水平的预测
我正在尝试拟合一个混合效应模型，然后使用该模型在可能具有不同水平的新数据集上生成估计值。我原以为在新数据集上的估计值将使用估计参数的平均值，但事实并非如此。以下是最小工作示例：library(lme4) d = data.frame(x = rep(1:10, times = 3), ...
glmer logit - 用概率刻度测量的交互作用效应（使用`predict`复制“effects”）
我正在使用lme4包运行glmer logit模型。我对各种二元和三元交互作用效应及其解释感兴趣。为简化问题，我只关心固定效应系数。 我设法编写了一段代码来计算并绘制这些效应的logit比例，但我在将它们转换为预测概率比例时遇到了困难。最终，我希望复制“effects”包的输出结果。 该示...
在Ubuntu上从CRAN安装lme4
我在Ubuntu 18.04上安装CRAN的lme4包时遇到了问题。我在R 3.4.4上运行install.packages("lme4")并收到以下错误：ERROR: configuration failed for package ‘nloptr’ * removing...