使用gtsummary :: tbl_svy_summary时如何报告未加权的缺失值数量？答案

【问题标题】：How to report the unweighted number of missing values when using gtsummary :: tbl_svy_summary?使用gtsummary :: tbl_svy_summary时如何报告未加权的缺失值数量？
【发布时间】：2021-06-03 18:10:42
【问题描述】：

我们分析加权调查数据，并希望报告分类变量的未加权频率（除了加权百分比）。但 statistic 参数中的 {n_unweighted} 选项不适用于缺失值。

请在下面找到将代码应用到包问题 r 中的 hdv2003 数据时显示的问题

library(questionr)
library(survey)
library(gtsummary)

data(hdv2003)
dp <- svydesign(ids=~1, weight= ~poids, data=hdv2003)
dp %>%
tbl_svysummary(
include=c("sexe", "nivetud"),
by="sexe",
statistic=all_categorical()~"{n_unweighted} ({p}%)",
missing = "ifany"
)

table of result

【问题讨论】：

您能否更新您的帖子以包含可重现的代码（即我们都可以在我们的机器上运行的代码）？另外，您指的是哪个未加权统计数据？您的示例中引用了三个。为此，请尽量减少您的示例，这意味着它仅包含说明您遇到的问题的最低限度。
抱歉，我已通过将代码应用于包问题器中包含的 hdv2003 数据来更新我的帖子

标签： r gtsummary

【解决方案1】：

无法通过统计参数修改未知行 N（好吧，根本无法修改，只报告加权 N！）。但是，您可以通过 add_n() 函数添加未加权的缺失 N。下面的例子！

library(questionr)
library(survey)
library(gtsummary)

data(hdv2003)
dp <- svydesign(ids=~1, weight= ~poids, data=hdv2003)
tbl <-
  dp %>%
  tbl_svysummary(
    include=c("sexe", "nivetud"),
    by="sexe",
    statistic=all_categorical()~"{n_unweighted} ({p}%)",
    missing = "no"
  ) %>%
  add_n("{N_miss_unweighted}") %>%
  modify_header(n ~ "**Unweighted N Missing**")

^{由reprex package (v2.0.0) 于 2021-06-04 创建}

【讨论】：

感谢您的回复！我们将添加未加权缺失 N 的列