Я следую замечательному учебнику Питера ДеВитта по qwraps2 и summary_table, но не могу продвинуться дальше.
Вот мои данные и код:
data(mtcars)
mtcars2 <- dplyr::mutate(mtcars,
cyl_factor = factor(cyl,
levels = c(6, 4, 8),
labels = paste(c(6, 4, 8), "cylinders")),
cyl_character = paste(cyl, "cylinders"),
gear_factor = factor(gear,
levels = c(3, 4, 5),
labels = paste(c(3, 4, 5), "gears")))
new_summary <- mtcars2 %>%
dplyr::select(.data$mpg, .data$wt, .data$gear_factor) %>%
qsummary(.)
by_cyl <- mtcars2 %>%
dplyr::group_by(.data$cyl_factor) %>%
summary_table(., new_summary)
В учебнике он вычисляет значения p для групповых сравнений и добавляет значения p в новый столбец в таблицах. Я хотел бы остановиться на этом, добавив больше результатов из сравнений (d Коэна и 95% -й ДИ, а также значение p). Затем я хотел бы добавить эти результаты в новую строку под каждой переменной, а не как новый столбец рядом с каждой переменной). Поэтому я бы хотел, чтобы результат выглядел примерно так (я составил цифры для теста сравнения групп):
6 cylinders (N = 7) 4 cylinders (N = 11) 8 cylinders (N = 14)
mpg
minimum 17.80 21.40 10.40
median (IQR) 19.70 (18.65, 21.00) 26.00 (22.80, 30.40) 15.20 (14.40, 16.25)
mean (sd) 19.74 ± 1.45 26.66 ± 4.51 15.10 ± 2.56
maximum 21.40 33.90 19.20
comparison d = 0.87, 95% CI [0.80, 0.94], p = 0.001
wt
minimum 2.62 1.51 3.17
median (IQR) 3.21 (2.82, 3.44) 2.20 (1.88, 2.62) 3.75 (3.53, 4.01)
mean (sd) 3.12 ± 0.36 2.29 ± 0.57 4.00 ± 0.76
maximum 3.46 3.19 5.42
comparison d = 0.87, 95% CI [0.80, 0.94], p = 0.001
Итак, у меня есть два вопроса:
-
Как добавить строку в таблицы и заполнить ее некоторым содержанием
Как запустить тесты сравнения групп, поместить их в правильный формат и ввести в Таблицы?
Моя главная проблема - это вопрос 1, это очень я застрял сейчас. Если мне удастся получить помощь в ее решении, я смогу самостоятельно разобраться с вопросом 2, поиграв с примерами Деввита в mpvals. Хотя я был бы очень рад получить помощь и по вопросу 2.
До сих пор я пытался добавить пустую строку в qsummary (), но не смог сделать это. Я пытался манипулировать матрицей символов, созданной summary_table, но не мог понять, как с ней работать. Любая помощь приветствуется!