Как сгенерировать t-значение, F-значение или хи-квадрат в сводной таблице с помощью пакета R "gtsummary"?

Question

Я работаю над созданием сводной таблицы с использованием отличного пакета R "gtsummary", он действительно очень помогает мне в эффективном и точном создании сводных таблиц. Но мне интересно, могут ли некоторые статистические данные, такие как t-значение, F-значение и хи-квадрат, автоматически генерироваться так же, как p-значение?

library(gtsummary)

add_p_ex1 <-
  trial[c("age", "grade", "response", "trt")] %>%
  tbl_summary(by = trt) %>%
  add_p()

Вот сводная таблица, созданная с помощью gtsummary.

Answer 1

К сожалению, нет простого способа добавить статистику теста в таблицу. Ниже я описал обходной путь, но это непросто. В gtsummary мы сосредотачиваемся на таблицах, которые готовы к публикации в журнале, а медицинские журналы, в которых я публикую свои собственные работы, не сообщают статистику тестов вместе с p-значениями.

Удачного кодирования!

1. Создайте сводную таблицу с tbl_summary() а также add_p().
2. Создайте вторую сводную таблицу, используя tbl_summary() а также add_p(), но на этот раз вам нужно будет предоставить настраиваемую функцию pvalue для add_p(): см. второй пример в add_p()файл. http://www.danieldsjoberg.com/gtsummary/dev/reference/add_p.tbl_summary.html
3. Пользовательская функция p-значения НЕ вернет p-значение, а вернет статистику теста.
4. Теперь у вас есть две сводные таблицы: одна с p-значением, а другая с тестовой статистикой. Объедините эти две таблицы, используяtbl_merge().
5. Теперь у вас есть таблица, которая включает все столбцы, которые вы ищете. Используйте функциюmodify_header()для обновления заголовка столбца тестовой статистики. Вы можете использоватьas_gt() чтобы преобразовать объект gtsummary в объект gt, затем gt::cols_hide() чтобы скрыть повторяющиеся столбцы.

Это аналогичный подход для создания таблицы ниже. Код для создания этой таблицы находится в галерее таблиц gtsummary. http://www.danieldsjoberg.com/gtsummary/dev/articles/gallery.html

Answer 2

mod_fmt_fun(statistic ~ style_sigfig) эта часть кода не работает. однако остальные дают результат с 4 d