#全患者の背景 比較群なし

tbl_1a <-
  CreateTableOne(vars=c("age_adm_d","code","sex","weight","los_JIPAD","adm_reason","adm_from","adm_transport","adm_urgent",
                        "NPPV_bu","HFNC_bu","ECMO_bu","CRRT_bu","Hyogo","after_hour","outcome_icu_dc","outcome_hp_dc", "PIM3_death_rate"),
                 factorVars =
                   c("code","sex","adm_reason","adm_from","adm_transport","NPPV_bu","HFNC_bu","ECMO_bu","CRRT_bu","Hyogo", 
                     "after_hour","outcome_icu_dc","outcome_hp_dc"),
                 data = merged_df)
print(tbl_1a,nonnormal =
        c("age_adm_d","weight","los_JIPAD"),showAllLevels = TRUE,includeNA = TRUE)
##                              
##                               level         Overall            
##   n                                           174              
##   age_adm_d (median [IQR])                   6.00 [3.00, 15.00]
##   code (%)                    no              174 (100.0)      
##   sex (%)                     male             96 ( 55.2)      
##                               female           78 ( 44.8)      
##   weight (median [IQR])                      2.80 [2.40, 3.10] 
##   los_JIPAD (median [IQR])                  10.00 [7.00, 14.75]
##   adm_reason (%)              post_op         127 ( 73.0)      
##                               Resp              6 (  3.4)      
##                               shock            28 ( 16.1)      
##                               monitoring       13 (  7.5)      
##   adm_from (%)                ward              4 (  2.3)      
##                               ER                8 (  4.6)      
##                               OR              129 ( 74.1)      
##                               referral         18 ( 10.3)      
##                               NICU             15 (  8.6)      
##   adm_transport (%)           in_hosp         148 ( 85.1)      
##                               self              3 (  1.7)      
##                               amb_referral     19 ( 10.9)      
##                               amb_kch           2 (  1.1)      
##                               amb_scene         1 (  0.6)      
##                               heli_referral     1 (  0.6)      
##   adm_urgent (%)              planned         119 ( 68.4)      
##                               urgent           55 ( 31.6)      
##                               procedure         0 (  0.0)      
##   NPPV_bu (%)                 No              168 ( 96.6)      
##                               Yes               6 (  3.4)      
##   HFNC_bu (%)                 No               64 ( 36.8)      
##                               Yes             110 ( 63.2)      
##   ECMO_bu (%)                 No              170 ( 97.7)      
##                               Yes               4 (  2.3)      
##   CRRT_bu (%)                 No              169 ( 97.1)      
##                               Yes               5 (  2.9)      
##   Hyogo (%)                   non Hyogo         9 ( 32.1)      
##                               Hyogo            19 ( 67.9)      
##   after_hour (%)              No              110 ( 63.2)      
##                               Yes              64 ( 36.8)      
##   outcome_icu_dc (%)          ward            117 ( 67.2)      
##                               transfer          2 (  1.1)      
##                               home              3 (  1.7)      
##                               death             8 (  4.6)      
##                               other_ICU         1 (  0.6)      
##                               CCU               8 (  4.6)      
##                               PICU             10 (  5.7)      
##                               NICU             25 ( 14.4)      
##   outcome_hp_dc (%)           alive           155 ( 89.1)      
##                               transfer          8 (  4.6)      
##                               death            11 (  6.3)      
##   PIM3_death_rate (mean (SD))                3.93 (5.83)

#胎児診断あり群となし群での比較

tbl_2 <-
  CreateTableOne(vars=c("age_adm_d","code","sex","weight","los_JIPAD","adm_reason","adm_from","adm_transport","adm_urgent",
                        "NPPV_bu","HFNC_bu","ECMO_bu","CRRT_bu","Hyogo","after_hour","outcome_icu_dc","outcome_hp_dc","PIM3_death_rate"), 
                 strata = "prenatal_diagnosis",
                 factorVars =
                   c("code","sex","adm_reason","adm_from","adm_transport","NPPV_bu","HFNC_bu","ECMO_bu","CRRT_bu","Hyogo", 
                     "after_hour","outcome_icu_dc","outcome_hp_dc"),
                 data = merged_df)
print(tbl_2,nonnormal =
        c("age_adm_d","weight","los_JIPAD"),showAllLevels = TRUE,includeNA = TRUE)
##                              Stratified by prenatal_diagnosis
##                               level         0                  
##   n                                            72              
##   age_adm_d (median [IQR])                   7.00 [3.00, 19.00]
##   code (%)                    no               72 (100.0)      
##   sex (%)                     male             46 ( 63.9)      
##                               female           26 ( 36.1)      
##   weight (median [IQR])                      2.80 [2.53, 3.02] 
##   los_JIPAD (median [IQR])                  12.00 [9.00, 14.50]
##   adm_reason (%)              post_op          44 ( 61.1)      
##                               Resp              5 (  6.9)      
##                               shock            14 ( 19.4)      
##                               monitoring        9 ( 12.5)      
##   adm_from (%)                ward              1 (  1.4)      
##                               ER                8 ( 11.1)      
##                               OR               45 ( 62.5)      
##                               referral         13 ( 18.1)      
##                               NICU              5 (  6.9)      
##   adm_transport (%)           in_hosp          51 ( 70.8)      
##                               self              3 (  4.2)      
##                               amb_referral     15 ( 20.8)      
##                               amb_kch           1 (  1.4)      
##                               amb_scene         1 (  1.4)      
##                               heli_referral     1 (  1.4)      
##   adm_urgent (%)              planned          43 ( 59.7)      
##                               urgent           29 ( 40.3)      
##                               procedure         0 (  0.0)      
##   NPPV_bu (%)                 No               71 ( 98.6)      
##                               Yes               1 (  1.4)      
##   HFNC_bu (%)                 No               25 ( 34.7)      
##                               Yes              47 ( 65.3)      
##   ECMO_bu (%)                 No               68 ( 94.4)      
##                               Yes               4 (  5.6)      
##   CRRT_bu (%)                 No               67 ( 93.1)      
##                               Yes               5 (  6.9)      
##   Hyogo (%)                   non Hyogo         7 ( 30.4)      
##                               Hyogo            16 ( 69.6)      
##   after_hour (%)              No               46 ( 63.9)      
##                               Yes              26 ( 36.1)      
##   outcome_icu_dc (%)          ward             53 ( 73.6)      
##                               transfer          1 (  1.4)      
##                               home              3 (  4.2)      
##                               death             4 (  5.6)      
##                               other_ICU         0 (  0.0)      
##                               CCU               3 (  4.2)      
##                               PICU              4 (  5.6)      
##                               NICU              4 (  5.6)      
##   outcome_hp_dc (%)           alive            62 ( 86.1)      
##                               transfer          6 (  8.3)      
##                               death             4 (  5.6)      
##   PIM3_death_rate (mean (SD))                4.32 (6.89)       
##                              Stratified by prenatal_diagnosis
##                               1                  p      test   
##   n                            102                             
##   age_adm_d (median [IQR])    6.00 [3.00, 13.00]  0.353 nonnorm
##   code (%)                     102 (100.0)           NA        
##   sex (%)                       50 ( 49.0)        0.074        
##                                 52 ( 51.0)                     
##   weight (median [IQR])       2.70 [2.30, 3.10]   0.217 nonnorm
##   los_JIPAD (median [IQR])    9.00 [6.00, 14.50]  0.026 nonnorm
##   adm_reason (%)                83 ( 81.4)        0.008        
##                                  1 (  1.0)                     
##                                 14 ( 13.7)                     
##                                  4 (  3.9)                     
##   adm_from (%)                   3 (  2.9)       <0.001        
##                                  0 (  0.0)                     
##                                 84 ( 82.4)                     
##                                  5 (  4.9)                     
##                                 10 (  9.8)                     
##   adm_transport (%)             97 ( 95.1)        0.001        
##                                  0 (  0.0)                     
##                                  4 (  3.9)                     
##                                  1 (  1.0)                     
##                                  0 (  0.0)                     
##                                  0 (  0.0)                     
##   adm_urgent (%)                76 ( 74.5)          NaN        
##                                 26 ( 25.5)                     
##                                  0 (  0.0)                     
##   NPPV_bu (%)                   97 ( 95.1)        0.407        
##                                  5 (  4.9)                     
##   HFNC_bu (%)                   39 ( 38.2)        0.754        
##                                 63 ( 61.8)                     
##   ECMO_bu (%)                  102 (100.0)        0.058        
##                                  0 (  0.0)                     
##   CRRT_bu (%)                  102 (100.0)        0.025        
##                                  0 (  0.0)                     
##   Hyogo (%)                      2 ( 40.0)        1.000        
##                                  3 ( 60.0)                     
##   after_hour (%)                64 ( 62.7)        1.000        
##                                 38 ( 37.3)                     
##   outcome_icu_dc (%)            64 ( 62.7)        0.080        
##                                  1 (  1.0)                     
##                                  0 (  0.0)                     
##                                  4 (  3.9)                     
##                                  1 (  1.0)                     
##                                  5 (  4.9)                     
##                                  6 (  5.9)                     
##                                 21 ( 20.6)                     
##   outcome_hp_dc (%)             93 ( 91.2)        0.138        
##                                  2 (  2.0)                     
##                                  7 (  6.9)                     
##   PIM3_death_rate (mean (SD)) 3.66 (4.97)         0.465