lme4 相关问题

lme4是一个R包,用于拟合和分析线性,非线性和广义线性混合模型。

BACI 设计:如何解释前后控制的差异?

我进行了一项 BACI 实验,研究夜间人造光对鱼类繁殖的影响,我们的响应变量之一是鱼卵中的蛋黄面积。我的原始数据如附件所示

回答 1 投票 0

手动添加回归

嗨! 我刚刚将回归手动添加到图中。但是如何去掉间隔之间的红色部分呢?

回答 1 投票 0

R 到 SAS:从 glmer 到 PROC NLMIXED

我知道我的问题可能不完全是关于通常的话题,但我真的很绝望。 几个月来我一直在试图找出我的代码有什么问题。 我需要与 PROC NLMIXED 等效的 SAS ...

回答 1 投票 0

在 Pymer4 中将残差分解为组间和组内

我有一个关于 Pymer4 中混合效果的问题(与 lme4 相同,但在 python 中)。我的模型在机器学习下拟合以获得预测值(Observed_values = Predicted_values + Total_residu...

回答 1 投票 0

在 R 中使用 lme4 混合效应模型时 !is.na(v.e) && v.e > 0 中出现错误

我没有找到与我的问题类似的内容。 我有来自 20 个不同实验的单细胞表达数据(29733 个基因 x 24489 个细胞)。我们想看看基因表达是否存在差异

回答 1 投票 0

在 Gamma GLMM 中处理空间自相关

问题概述 我正在使用一个数据集来检查植物生长与 各种因素,包括 14 个不同的海冰范围 (SeaIce) 采样点。为了分析数据,我...

回答 1 投票 0

使用 ggplot2 进行分层模型 (lmerMod) 的交互图,并突出显示统计显着性

我正在使用一个分层模型(lmerMod),它有一个结果和两个分类变量,cat1 和 cat2。我想使用 ggplot2:: 制作一个交互图来可视化模型

回答 1 投票 0

找出某个协变量相对于 R 中所有其他协变量和交互作用的平均效果

我正在 R 中运行一个线性混合模型,使用一个项目的重复测量数据,作为我分析的一部分,我需要确定时间的影响是什么(结果与一个相比有多少变化

回答 1 投票 0

将 aov 翻译为 lmer

我被困在这里了。如何将以下命令行转换为 lmer 语法。 型号 <- aov( response ~ treatment*time + Error(id/time), data = selfesteem2) summary(model) I've tried this among

回答 1 投票 0

指定嵌套在固定效应相互作用下的随机效应

可能很简单。 我有具有固定和随机效应的数据,我想将混合效应模型拟合到: 设置.种子(1) df <- data.frame(group = c(rep("A",40),rep("B",40)), tre...

回答 1 投票 0

如何在函数内安全地将惰性对象导出到函数中的工作人员?

我在函数 g 中有一个函数 f,它首先操纵来自 lm、lme4::lmer 或 lmerTest::lmer 的调用并重新评估它。 f 将使用 vapply(类型 1)或 mclapply(类型...

回答 1 投票 0

如何用 2x2 实验设计写出 GLM 公式

我进行了一项研究,实验设计是2x2。 我有一个名为 strum 的自变量,它是一个两级因子,并且位于参与者内部(因此两个

回答 1 投票 0

嵌套后无法访问因变量以拟合 lmer 模型

我一直在尝试拟合这样的模型 mtcars %>% 变异(vs = 因子(vs), gear = 因子(gear), carb = 因子(carb), am = 因子(am), cyl = 因子(cyl)) %>% Nest_by(vs) %>% ...

回答 1 投票 0

使用 lmer 绘制与 sjPlot 交互的多级建模

我正在运行一个多级模型,其中包含两个连续的 2 级预测变量(A 和 B)和一个二分的 1 级预测变量(条件,编码为 -1 和 1)。 Y 是因变量。 M2<-lmer(Y ~ 1 + (1|sub...

回答 1 投票 0

在使用混合效应模型时如何解释样本随时间的重复测量

我有一个实验,在三年的时间里,我在春季和夏季反复从同一块地的内部和外部采集沉积物样本。例如,sample_IA 表示来自 plo 内部的样本...

回答 1 投票 0

使用 R 的 lme4 包形成尽可能多的混合效应模型

我想使用 R 中的 lme4 包和 lmer() 函数形成尽可能多的模型。 图书馆(tidyverse);图书馆(lme4); 门 = read_csv('gate_data.csv') m1=lmer(CHG ~ 年*BLAREA + 年:LOC ...

回答 1 投票 0

如何提取相当于lme4::VarCorr()的metafor::rma.mv(),用于具有随机斜率的多级混合效应元回归

我正在尝试使用与...类似的方法计算多级混合效应元回归的伪 R 平方,其中包括 metafor 包(即 rma.mv() 对象)中的随机斜率

回答 1 投票 0

如果我重新调整了回归量,如何使用 visreg 函数绘制预测响应和 95%CI?

当我缩放回归量时,有什么方法可以使用 visreg 函数制作条件图吗? 我想绘制重新调整后的响应以及 95% CI、预测曲线和部分残差。 我...

回答 1 投票 0

平均分类数据的系列/链接

我有一些分类数据,在观察级别编码为 0/1。我对试验中的这些值进行了平均,这为我提供了 0 到 1 之间的概率值。 我正在运行 GLME...

回答 1 投票 0

lmer 输出的方差分量的标准误差

我需要从 lmer 的输出中提取方差分量的标准误差。 图书馆(lme4) 型号 <- lmer(Reaction ~ Days + (1|Subject), sleepstudy) The following produces estimates of

回答 3 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.