我正在尝试拟合一个混合效应模型，然后使用该模型在可能具有不同级别的新数据集上生成估计值。我预计新数据集的估计值将使用估计参数的平均值，但事实似乎并非如此。下面是一个最低限度的工作示例： librar...

为什么我不能在线性混合效应模型上运行 Breusch-Pagan 检验来检验异方差性？该函数适用于使用 和 构建的模型，但不能。我应该使用其他功能吗？bptest()lmer()bptestlmglm...

我正在以下问题上寻求一些帮助。 我必须分析SNPs与抑郁症状之间的关联。SNP 的遗传易感性将被归类为 0,1 或 2 个次要等位基因（暴露变量）。抑郁症状将随着时间的推移（T0、T1、T2、T3）进...

我一直在尝试拟合这样的模型 mtcars %>% mutate(vs = factor(vs), gear = factor(gear), carb = factor(carb), am = f...

我正在使用一个分层模型 （lmerMod），该模型具有一个结果和两个分类变量，以及 .我想制作一个交互作用图，用它来可视化模型系数并考虑两个特定的复杂性：cat1cat2ggplot2:: 成对比较...

我收到 PIRLS 步进减半错误，使用 .在我看来，我生成的数据可能在下面模型过于复杂，这可能是完整的答案，但我对错误注释的理解不够好，无法自信。我看到了另一个答案，它引用了不太规范的链接函数的行为，...

不幸的是，没有一个可重现的例子：我有一个包含 480,000 个观测值的数据帧。由于数据的嵌套结构，我使用依赖于 lme4 包的分层模型。如果我运行线性模型，则没有错误。但是，由于我有一个二分因变量，...

我在 R 中运行。glmer() 这是我的空模型 empty <- glmer(DV ~ +(1|year), weights = weights, data= df, famil...