Как присвоить значения N, логарифмического правдоподобия, AIC и BIC каждой модели многомерной регрессии в выходных данных объединенной таблицы `gtsummary`?

Пожалуйста, как мне назначить значения N, логарифма правдоподобия, AIC и BIC для каждой модели многомерной регрессии в выходных данных объединенной gtsummary таблицы?

Я Rmerteur пытаюсь подогнать аддитивные многомерные модели, добавляя переменные более высокого уровня к последующим моделям (см. Пример ниже). Между тем, я не мог сохранить значения N, логарифма правдоподобия, AIC и BIC для каждой подходящей модели после того, как я объединю отдельные модели. Мне нужны оценки, чтобы сравнить соответствие между моделями.

Есть ли простой способ сделать это?

Спасибо за помощь!

m1<-
  tbl_regression(mod1,
    exponentiate = T, hide_n = T, 
    add_estimate_to_reference_rows = T,
    estimate_fun = ~style_ratio(.x, digits = 2), 
    pvalue_fun = ~style_pvalue(.x, digits = 2), 
    label = list(age_yr_grp ~ "Age",
                 educ ~ "Education")) %>%
  bold_p(t = 0.05) %>% 
  add_glance_source_note()  

m2<-
  tbl_regression(mod2,
    exponentiate = T, hide_n = T, 
    add_estimate_to_reference_rows = T,
    estimate_fun = ~style_ratio(.x, digits = 2), 
    pvalue_fun = ~style_pvalue(.x, digits = 2), 
    label = list(age_yr_grp ~ "Age",
                 educ ~ "Education",
                 hhwealth ~ "Household wealth",
                 hhsize ~ "Household size")) %>%
  bold_p(t = 0.05) %>% 
  add_glance_source_note()  

m3<-
  tbl_regression(mod3,
    exponentiate = T, hide_n = T, 
    add_estimate_to_reference_rows = T,
    estimate_fun = ~style_ratio(.x, digits = 2), 
    pvalue_fun = ~style_pvalue(.x, digits = 2), 
    label = list(age_yr_grp ~ "Age",
                 educ ~ "Education",
                 hhwealth ~ "Household wealth",
                 hhsize ~ "Household size",
                 placeres~ "Area of residence")) %>%
  bold_p(t = 0.05) %>% 
  add_glance_source_note()     


tab3_multivars_models= 
  tbl_merge(list(m1,m2,m3),
  tab_spanner= c("**Mod 1**","**Mod 2**","**Mod 3**"))%>%
  as_gt() %>%
  tab_header(title = "Table title here") %>%
  tab_source_note(gt::md("*Notes here"))

gt::gtsave(tab3_multivars_models, "dir\\tab3_multivars_models.html", inline_css = TRUE)

r gtsummary
person user14896383    schedule 17.04.2021    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Вы захотите использовать add_glance_table() вместо версии исходной заметки.

library(gtsummary)
#> #Uighur
packageVersion("gtsummary")
#> [1] '1.4.0'

tbl1 <-
  lm(age ~ marker + grade, trial) %>% 
  tbl_regression() %>%
  add_glance_table(
    label = list(sigma ~ "\U03C3"),
    include = c(r.squared, AIC, sigma)
  )

tbl2 <- tbl_merge(list(tbl1, tbl1))

введите здесь описание изображения Создано 17 апреля 2021 года пакетом репекс < / а> (v2.0.0)

person Daniel D. Sjoberg    schedule 17.04.2021
comment
Сработало как по волшебству !!!! Большое тебе спасибо. @ Дэниел Д. Шоберг - person user14896383; 17.04.2021