-
1 votesanswersviews
随机效应中的交互术语和生长模型中的白天编码lme4
我对R中的模型设置有疑问，经过长时间彻底的搜索，我没有找到回答我的两个问题的线程： 我将首先描述设置： 它是一个重复测量数据集，有两种不同的干预措施（食物和培训），每种干预都有两个层次 . 所有参与者都经历了各种条件的组合 . 在一天中收集激素样品 . 年龄和BMI作为协变量包括在内 . 这是一个可重复的数据集： ID <- rep(rep(c("A","B&q... -
0 votesanswersviews
R：混合（多级）模型（lme4）：固定效应和随机效应之间的相互作用
我正在分析一个纵向数据集，其中时间t（级别1）的观测值嵌套在单位i（级别2）内 . 具有随机截距和斜率的基线模型看起来像这样： y_ti = beta_0i + beta_1i * x1_ti + e_ti 同 beta_0i = gamma_00 + gamma_10 * a_i + u_0i beta_1i = gamma_10 + gamma_11 * a_i + u_1i 和二阶误差项... -
0 votesanswersviews
模拟“zelig风格”用于r中的GLMER多级
我用r运行逻辑混合效应回归 . 回归在某种程度上是这样的： glmer（Y~X1 X2 X1：X2（1 | country），data = hdp，family = binomial） 现在，通过固定效果，我想绘制Y的预测概率 . 我尝试使用Zelig，因为这是我学习的最简单的模拟方法和预测概率，但我看到新版本不包括多级模型和前Zelig Multilevel非常“不稳定” . 有没有简单的选择？... -
0 votesanswersviews
多级数据的多重插补
我在R中使用Mice包来多次插入多级数据，其中重复测量嵌套在个体内 . 但是在小鼠中存在一个错误，我们需要将组变量（此处为个体）转换为整数类型，以便使用“mice.2l.norm / mice.2l.pan”方法对小鼠进行插补以考虑聚类 . （我无法发布我的数据的快照）我的分组变量（STIDnum）是一个因子变量，因为我想检查插补前后的ICC值 对于Imputation我已将STIDnum转换为整...