为了测试疾病组(categorical_variable)和疾病(结果;计数)之间是否存在关联,我运行了3个负二项回归模型。
为了显示OR和CI,我使用了gtsummary
包中的tbl_regression
函数。然而,这个函数显示用逗号分隔的配置项,我希望它们放在括号中,并显示p值,我希望将它们全部抑制在一起。
我的代码:
library(gtsummary)
model <- glm(data=data, formula=outcome ~ categorical_variable)
tbl_regression(model,
exponentiate = TRUE,
include="categorical_variable")
在格式化配置项和抑制p值列方面有什么帮助吗?
我尝试添加%>% as_gt() %>% cols_hide("p-value")
和%>% as_gt() %>% cols_hide(columns=vars("p-value"))
,但都无济于事。它说它不将p值识别为列(没有括号也不起作用)。
发布于 2021-04-07 11:26:40
如果res
是上述问题中的对象,即gtsummary
对象,则可以修改查看器窗格中显示的表体,如下所示。在这里,我使用dplyr
的select
函数(modify_table_body
需要一个函数)来排除p值列。
gtsummary::modify_table_body(res, dplyr::select, -p.value)
我只想使用base
p_val_col<-which(names(res$table_body)=="p.value")
gtsummary::modify_table_body(res, `[`, -p_val_col)
要修改ci
的显示(您可能可以编写一个比这个简单得多的正则表达式),可以按如下方式运行(同时)。你可以一次写一个函数来完成这个任务,而不是两次调用modify_table_body
:
gtsummary::modify_table_body(res, `[`, -23) %>%
gtsummary::modify_table_body(., dplyr::mutate,
ci=gsub("(\\d\\.\\d{,4})(, )(\\d\\.\\d{,4})"
,"\\[\\1 \\3\\]",ci))
发布于 2021-11-30 14:23:06
正如Elin在单表情况下所提到的:
tbl_regression(model,
exponentiate = TRUE,
include="categorical_variable") %>%
modify_column_hide(columns = p.value)
对于具有两个(或更多)回归模型的情况:
tbl_merge(
tbls = list(tbl_regression(model1),
tbl_regression(model2))
) %>%
modify_column_hide(columns = c(p.value_1, p.value_2)
https://stackoverflow.com/questions/66984675
复制相似问题