如何在 R 粘合語法中設置“na.rm = TRUE”

Question

我正在使用gtsummary::tbl_summary創建一個統計匯總表。

stats_summary_table <-
  dat_wide%>%
  tbl_summary(by = mss_plate_id, 
              missing = "no",
              statistic = list(all_continuous() ~ " {mean} ± {sd}   {cv}"  ))

cv function 來自goeveg庫。 報告表顯示 cv 值為 0。我想這是因為我的數據中缺少值。 所以，我的問題是如何用膠水語法編寫cv(x, na.rm = TRUE) 。

編輯：事實證明，挑戰是由於精度問題而不是設置na.rm=TRUE 。 因此解決方案是將 cv 的位數設置為 2 或 3。 （請參閱標記的答案及其評論）

Answer 1

這是一個很好的問題。

tbl_summary()嘗試做的一件困難的事情是猜測要將統計信息四舍五入到多少位數。 這個例子很棘手，因為您通常希望以不同的精度級別顯示均值和 CV，因為它們不在同一比例上。 在下面的示例中， tbl_summary()猜測年齡應該顯示為四舍五入到最接近的 integer。 這是均值和 SD 的合理假設，但我們需要告訴它為 CV 顯示更多數字。 為此，我們使用tbl_summary(digits=)參數。 顯示了年齡的三個統計數據（平均值、標准差和 cv），我們將傳遞一個長度為 3 的向量，指示將每個統計數據四舍五入的位數。

library(goeveg)
library(gtsummary)

trial %>%
  select(age, marker) %>%
  tbl_summary(
    statistic = list(all_continuous() ~ " {mean} ± {sd} [{cv}]"),
    digits = list(age ~ c(0, 0, 2))
  )

希望這能回答你的問題！ 快樂編碼！