从lmer输出中外推单位降低和百分比降低的95% CI

k4ymrczo  于 2023-05-04  发布在  其他
关注(0)|答案(1)|浏览(98)

我正在测试不透水表面积(伊萨)的百分比和每次繁殖事件的鸡蛋数量之间的关联。这是我的模型结构:

fit=lmer(CS~ISA+Mass_Fem*LayDate+Year+(1|Site),data=DataGT_noNA)

其中CS是窝卵数,伊萨是不透水性(以%计),雌性质量是以克为单位的雌性体重,产蛋日期是序数,年份是5个水平因子。随机效应“站点”有8个水平。
使用ggplot和ggeffects获得某个%伊萨的损失估计值后,如下所示:

effects_ISA <- effects::effect(term= "ISA", mod= fit, xlevels=list(ISA=c(0,10,20,30,40,50,60,70)))

我创建了一个包含所有这些值的数据框,如下所示
x_伊萨〈- as.data.frame(effects_伊萨),并获得此表作为结果:

我想获得给定伊萨百分比的鸡蛋数量(包括95% CI)的单位减少(和相应的置信区间)。参见以下示例:
伊萨的30%的增加对应于0。54(CI 95??)卵数减少。我想获得相同的信息作为百分比(%损失和% CI95%)

tv6aics1

tv6aics11#

我使用R包marginaleffects解决了这个问题:
具体来说,我写了这样一段代码:

avg_comparisons(mymodel, variables = list(ISA = 30))

得到了这些结果:

Term Contrast Estimate Std. Error     z Pr(>|z|) 2.5 %  97.5 %
  ISA      +30   -0.537      0.252 -2.13    0.033 -1.03 -0.0432

相关问题