与随机效应项统计模型分析相关的编程问题,也有不同的:重复测量,分层,多层次模型
如何在 lmer 混合模型中获得每个随机效应水平的 95% CI
我在 lme4 中运行了几个混合模型,其中包含一些主效应、一些交互项和两个随机效应(animalID 和年份)。我的模型看起来像这样: 型号 <- lmer(response ~ ns(area_sc...
考虑这个玩具数据集,使用交叉随机效应进行模拟,然后使用 lme4::lmer 拟合一个混合模型: 设置.种子(123) # 科目和项目数量 科目数 <- 30 num_items &...
我有以下概念模型: IV:人格特质(以李克特量表衡量) 主持人:低信任度与高信任度 DV:谈判中的初始报价(值可以在 5-15 之间) 我...
sommer 包相当于 lmer(反应 ~ 天数 + (天数 | 主题), 数据 = sleepstudy)
在 lme4 中,lmer(Reaction ~ Days + (Days | subject), data = sleepstudy) 指定一个线性混合模型,该模型具有主题的随机截距和天数的随机斜率,具有相关的随机截距和 r...
在 GPBoost 的 GPModel 中创建 group_rand_coef_data 的适当方法
下面是我从 GPBoost 站点复制的代码。 https://github.com/fabsig/GPBoost/blob/master/examples/python-guide/generalized_linear_Gaussian_process_mixed_effects_models.py 我认为输入
我有以下数据: str(增长数据) 小标题 [92 × 4] (S3: tbl_df/tbl/data.frame) $ 人:数字 [1:92] 1 1 1 1 2 2 2 2 3 3 ... $ 性别:chr [1:92]“F”“F”“F”...
我一直在缺失值的数据集上运行lmer。我已经成功地估算了缺失值并成功运行了我的模型并获得了汇总拟合的摘要。但是,我需要方差分析
R sommer 包:用于基因组预测的带有 Predict.mmer 的 Dtable 错误
我有一个多环境试验数据集,在四个环境中测试基因型。我想预测所有基因型的环境特异性 BLUP。我使用 sommer-p 中的 mmer 函数...
glmmTMB - AR1协方差结构 - 条件公式和零膨胀公式中参数数量不同
我试图理解为什么在条件公式中指定 AR-1 协方差结构估计 2 个参数(就像它应该的那样),但在零通货膨胀中指定时估计许多参数......
使用 R 中的模拟数据使用混合效应模型计算多中心 RCT 的样本量
我目前正在建立一项治疗与安慰剂(1:1)的多中心(九个中心)随机对照试验,主要结果是延长妊娠天数。最小的临床相关差异...
使用 Stata 建立反映年龄和性别混合效应线性回归模型作为固定效应和体重作为响应的图表的步骤 cpde 是什么? 混合体重年龄_月数....
R 中的逻辑混合模型收敛警告,但 Stata 中没有 [已关闭]
我使用带有 lme4 的单变量逻辑混合模型来研究血浆参数 (M) 对协变量疾病结果的影响: 公式<- Disease ~ Age + Sex + BMI + M + (1 | Family) The sample ...
为什么我在运行零膨胀负二项式混合模型时收到“二元运算符的非数字参数”错误?
我正在尝试运行零膨胀负二项式混合模型。我进行了一项实验,给成年蜗牛提供富磷饮食和贫磷饮食(P)。这些蜗牛都是
以“时间”作为连续固定效应,但分类随机效应来拟合广义线性混合效应模型是否有效?
我想使用 R 包 glmmTMB 中的 AR1 相关结构来拟合广义线性混合效应模型。检查文档后 ……看来“时间”一定是……
我正在使用 lme4 包进行线性混合效应建模 混合效应模型为: FM01 <- lmer(sublat <- goal + (1|userid)) The above command returns an S4 object called fm01. This obj...
这个问题是关于如何从S4 lme模型摘要对象中删除部分字符串。 假设这是我的模型。 # ----预备工作-------------------------------------------- ------------ 图书馆...
nlminb 问题,收敛错误代码 = 1 消息 = 达到迭代极限但未收敛 (10)
我正在尝试使用 R 中的库(nlme)和 lme 函数找到适合我的数据的最佳模型。这是当斜率固定时我的模型: 固定绳长 <- lme(EnergyCost~ RopeLength, data = ...
我正在尝试在Python中实现混合效应逻辑回归。作为比较,我使用 R 中 lme4 包中的 glmer 函数。 我发现 statsmodels 模块 h...
使用非线性混合效应模型 (nlme) 的预测变量进行预测的置信区间
我正在尝试使用 Ben Bolker 描述的 bootstrap 方法为非线性混合效应模型生成 95% 的置信区间。我可能误解了所使用的一些功能。我...