Compare nested multiply imputed models

gistfile1.txt

      ## Likelihood test p values for categorical variables
      model_uni <- outcomes %>% map(function(.) model_outcome_uni(., variable = 'ageatvac', LRT = TRUE))
      model_adj <- outcomes %>% map(function(.) model_outcome_adj(., variable = 'ageatvac', LRT = TRUE))
   
      ageatvac_lrt_uni <- map2(model_ageatvac_uni, model_uni, function(x, y) pool.compare(x, y, data = imp_df, method = "likelihood")) %>% 
    map(function(x) ifelse(x < 0.001, "0.001", pretty_round(x, digits = 3))) 
      ageatvac_lrt_adj <- map2(model_ageatvac_adj, model_adj, function(x, y) pool.compare(x, y, data = imp_df, method = "likelihood")$pvalue) %>% 
        map(function(x) ifelse(x < 0.001, "0.001", pretty_round(x, digits = 3)))