我有一个Logistic回归模型预测兴趣(0=无兴趣,1=兴趣),随机实验条件(0=对照,1=实验),性别,性别:条件交互作用,以及其他几个协变量作为预测因素.

在R中使用marginaleffects包,我需要解决两件事:

  1. 性别特定风险比(RR),反映实验条件对利息的影响(相对于其他协变量而言处于边缘地位).
  2. 交互作用项(即风险比:RR男性/RR女性),表示男性的RR是女性的多少倍.

以下是Logistic回归模型:

model <- glm(data=dt, formula=interest ~ condition*sex + other_covariates,
              family=binomial(link="logit"))

我已经成功地完成了目标(1)--计算按性别划分的边际风险比率(RR男性和RR女性)--使用以下方法:

marginaleffects::avg_comparisons(model, comparison="ratio",
                variables="condition",by="sex")

然而,我似乎想不出如何计算(2)--风险比率(RR男性/RR女性)--及其相关的可信区间.

任何帮助都将不胜感激.如果使用另一个像emmeans这样的包更简单,我也很乐意使用它.谢谢!

推荐答案

您可以使用hypothesis参数对marginaleffects包生成的任何量进行(非线性)线性测试.请参阅此处非常详细的教程:

https://marginaleffects.com/vignettes/hypothesis.html

您没有提供一个最小的工作示例,所以我不能提供一个完全工作的解决方案,但是您的代码最终可能如下所示:

avg_comparisons(model, 
   comparison="ratio",
   variables="condition",
   by="sex",
   hypothesis = "b1 / b2 = 1")

请注意,我只添加了最后一个参数,告诉marginaleffects我们要判断第一个估计值与第二个估计值的比率是否不同于1.当我说"第一个"估计值时,我指的是同一命令输出中的第一行(如果您只删除hypothesis个参数).

R相关问答推荐

从具有随机模式的字符串中提取值

查找满足SpatRaster中条件的单元格位置

是否可以 Select 安装不带文档的R包以更有效地存储?

按R中的组查找相邻列的行累积和的最大值

如果列中存在相同的字符串,则对行值进行总和

在R中创建一个包含转换和转换之间的时间的列

pickerInput用于显示一条或多条geom_hline,这些线在图中具有不同 colored颜色

将嵌套列表子集化为嵌套列表

在R中为马赛克图中的每个字段着色

在R中,如何将变量(A,B和C)拟合在同一列中,如A和B,以及A和C在同一面板中?

TreeNode打印 twig 并为其上色

有没有办法一次粘贴所有列

为什么我对圆周率图的蒙特卡罗估计是空的?

删除数据帧中特定行号之间的每第三行和第四行

如何显示准确的p值而不是<;0.001*?

如何移动点以使它们的打印不重叠

R-找出存在其他变量的各种大小的所有组合

将R中对象的CSV数组转换为JSON数组

列间序列生成器的功能

R中的交叉表