我正在测试不透水表面积(伊萨)的百分比和每次繁殖事件的鸡蛋数量之间的关联。这是我的模型结构:
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%)
1条答案
按热度按时间tv6aics11#
我使用R包marginaleffects解决了这个问题:
具体来说,我写了这样一段代码:
得到了这些结果: