在nlme和lme4中拟合相同的模型
数据来自这里。 library(nlme) dat0 <- read.table("aids.dat2",head=T) dat1 <- dat0[dat0$day<=90, ] # use only first 90-day data dat2 <- dat1[...
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::lmer 在以下格式的 6000 万个观测值上运行混合模型；所有预测/依赖变量都是分类变量（因子），除了连续依赖变量 tc；patient 是随机拦截项的分组变量。我有 64 位 R 和 16Gb RAM，并且从中央服务器工作。RStudio 是最新的服务器版本。 mo...
在R中（使用lme4），将时间序列纳入混合效应模型
我已经搜索了类似的问题，但没有找到相关的问题，如果我错过了相关的问题，请谅解。 我正在研究在不同条件下饲料器使用时间（因变量），每个受试者访问30次饲料器。 受试者会接触一种类型的饲料器，该饲料器具有不同的组合，包括有香味/无香味、有视觉图案/无视觉图案，以及这些视觉或香味图案在两种空间安...
在LME模型中，第0级，第1块背推中存在奇异性。
dput用于数据传输，从https://pastebin.com/1f7VuBkx复制（太大无法在此处包含） data.frame': 972 obs. of 7 variables: $ data_mTBS : num 20.3 22.7 0 47.8 58.7 ... $ da...
从元分析的混合效应多层模型中获取R平方
我正在使用R进行特定治疗对森林的元分析。为此模型，我需要拟合随机效应以解释方法和地点年龄差异，因为这两个是混淆变量，而我不想明确调查它们引起的变异。但据我所知，包[metfor]不允许在多层模型中计算类似于R平方的统计量。无论如何，为了更清楚地描述我的问题，这里有一个模拟数据集。 Log&l...
在lme模型拟合中提取每个单位的系数及其标准误差
在类似于以下的线性混合模型中，我如何提取每个实验单位（plot）的系数（b0和b1）及其相应的标准误差： Better fits for a linear model 使用与此数据集（df）相同的数据，以及拟合的模型（fitL1）：我如何获得像这样的数据帧... plot b...
如何在nlme和lme4中指定不同的随机效应？
我希望使用nlme::lme在模型中指定不同的随机效应（底部有数据）。这些随机效应是：1）intercept和position在subject上变化；2）intercept在comparison上变化。使用lme4::lmer很容易实现： lmer(rating ~ 1 + position...
lme4::glmer.nb函数根据我运行模型的顺序，会出现“Error in family$family : $ operator not defined for this S4 class”的错误提示。
library(lme4) dummy <- as.data.frame(cbind(speed = rpois(100, 10), pop = rep(1:4, each = 25), season = rep(1:2, each = 50), id = seq(1, 100, by...