1. Import library

suppressMessages(library(data.table))
suppressMessages(library(tidyverse))
suppressMessages(library(mice))
suppressMessages(library(naniar))

2. Import dataset

## [1] 7672  169
## [1] 51558   151

3. Create new dataset consists of baseline and follow-up data

dat1 <- mutate_all(dat1, as.character)
dat2 <- mutate_all(dat2, as.character)

# Select common variables
common_vars <- intersect(names(dat1), names(dat2))

# Select different variables from dat1
diff_vars_dat1 <- setdiff(names(dat1), common_vars)

# Select different variables from dat2
diff_vars_dat2 <- setdiff(names(dat2), common_vars)


# Create a new dataset with common and different variables
dat3 <- merge(dat1, dat2, by = common_vars, all = TRUE)

# Save dat3 with the code below
# write.csv(dat3, file = "~/Desktop/Dat3.csv", row.names = FALSE)

4. Checking missing

4.1 For baseline

# Calculate missing data
dat1 %>%
  summarise_all(~sum(is.na(.))) %>%
  gather(key = "Variable", value = "MissingCount") %>%
  mutate(MissingPercentage = MissingCount / nrow(dat1) * 100) %>%
  arrange(desc(MissingPercentage))
##      Variable MissingCount MissingPercentage
## 1        uwbc         7672      100.00000000
## 2       d_stk         7672      100.00000000
## 3        d_mi         7672      100.00000000
## 4      f_lost         7672      100.00000000
## 5         ph8         7671       99.98696559
## 6     dth_dur         7666       99.92179353
## 7         ph7         7666       99.92179353
## 8         ph6         7651       99.72627737
## 9         paf         7617       99.28310740
## 10     dth_mi         7606       99.13972888
## 11        ph5         7574       98.72262774
## 12    mi_area         7572       98.69655892
## 13      mine5         7569       98.65745568
## 14     af_day         7512       97.91449426
## 15    int_knd         7512       97.91449426
## 16   int_oped         7508       97.86235662
## 17      mine4         7449       97.09332638
## 18   ci_scale         7405       96.51981230
## 19    ci_type         7396       96.40250261
## 20    ci_area         7396       96.40250261
## 21        ph4         7368       96.03753910
## 22   mi_onset         7367       96.02450469
## 23         ct         7340       95.67257560
## 24        mri         7339       95.65954119
## 25   st_onset         7234       94.29092805
## 26      mine3         7116       92.75286757
## 27     tobsty         6991       91.12356621
## 28     cp_dur         6935       90.39363921
## 29        ph3         6858       89.38998957
## 30      tobyr         6449       84.05891554
## 31   ido_date         6138       80.00521376
## 32     inslin         5949       77.54171011
## 33      mine2         5908       77.00729927
## 34      seiri         5824       75.91240876
## 35        ph2         5819       75.84723670
## 36    death_d         5801       75.61261731
## 37        ldl         5543       72.24973931
## 38       vldl         5532       72.10636079
## 39     tobhon         5336       69.55161627
## 40    heikeiy         5126       66.81438999
## 41      hba1c         4930       64.25964546
## 42      memai         3846       50.13034411
## 43   nocturia         3846       50.13034411
## 44   cp_sleep         3844       50.10427529
## 45   alc_freq         3809       49.64807091
## 46     alcgou         3787       49.36131387
## 47        ph1         3750       48.87904067
## 48     heikei         3640       47.44525547
## 49   deth_moy         3490       45.49009385
## 50   deth_fay         2118       27.60688217
## 51       wt20         2038       26.56412930
## 52    icd_stk         1704       22.21063608
## 53    icd_ihd         1622       21.14181439
## 54    icd_cvd         1394       18.16996872
## 55       hako         1378       17.96141814
## 56      elbow         1376       17.93534932
## 57         bf         1122       14.62460897
## 58    icd_can         1115       14.53336809
## 59    unhappy          841       10.96193952
## 60     stress          805       10.49270073
## 61     satisf          793       10.33628780
## 62       drug          758        9.88008342
## 63      sleep          731        9.52815433
## 64      vital          721        9.39781022
## 65     social          721        9.39781022
## 66     arith2          719        9.37174140
## 67      dx_sd          711        9.26746611
## 68     dx_pmi          685        8.92857143
## 69     dx_sah          680        8.86339937
## 70     arith1          643        8.38112617
## 71     dx_ich          638        8.31595412
## 72      dx_mi          631        8.22471324
## 73     dx_int          605        7.88581856
## 74      scott          590        7.69030240
## 75    sheieht          579        7.54692388
## 76    sheieas          578        7.53388947
## 77       frct          479        6.24348279
## 78      bruit          477        6.21741397
## 79      dx_ci          441        5.74817518
## 80        alb          420        5.47445255
## 81        glb          420        5.47445255
## 82     dx_ihd          412        5.37017727
## 83     kaidan          388        5.05735141
## 84       exer          312        4.06673618
## 85      dx_st          279        3.63660063
## 86        tob          130        1.69447341
## 87        alc          121        1.57716371
## 88  dysphagia           78        1.01668405
## 89        uph           58        0.75599583
## 90       mahi           58        0.75599583
## 91     murmur           57        0.74296142
## 92   sft_kata           53        0.69082377
## 93    deth_mo           52        0.67778936
## 94    deth_fa           43        0.56047967
## 95  chestpain           41        0.53441084
## 96        hdl           35        0.45620438
## 97        rbc           27        0.35192909
## 98        wbc           26        0.33889468
## 99         hg           26        0.33889468
## 100       hct           26        0.33889468
## 101    pulse2           23        0.29979145
## 102   sft_ude           23        0.29979145
## 103       ztt           20        0.26068822
## 104       cpk           20        0.26068822
## 105      fast           19        0.24765381
## 106     waist           18        0.23461940
## 107       hip           18        0.23461940
## 108   agratio           17        0.22158498
## 109       gpt            5        0.06517205
## 110        ca            5        0.06517205
## 111    pulse1            4        0.05213764
## 112    height            3        0.03910323
## 113    weight            3        0.03910323
## 114        ua            3        0.03910323
## 115      sbp2            2        0.02606882
## 116      dbp1            2        0.02606882
## 117      dbp2            2        0.02606882
## 118        tp            2        0.02606882
## 119        tc            2        0.02606882
## 120        tg            2        0.02606882
## 121       got            2        0.02606882
## 122      rgtp            2        0.02606882
## 123      crea            2        0.02606882
## 124        bs            2        0.02606882
## 125     mine1            2        0.02606882
## 126      sbp1            1        0.01303441
## 127       fh1            1        0.01303441
## 128       fh2            1        0.01303441
## 129       fh3            1        0.01303441
## 130       fh4            1        0.01303441
## 131       fh5            1        0.01303441
## 132       fh6            1        0.01303441
## 133       fh7            1        0.01303441
## 134       fh8            1        0.01303441
## 135       age            0        0.00000000
## 136       sex            0        0.00000000
## 137       usg            0        0.00000000
## 138       upr            0        0.00000000
## 139       uob            0        0.00000000
## 140     dx_ht            0        0.00000000
## 141    dx_tia            0        0.00000000
## 142     dx_bd            0        0.00000000
## 143    dx_chf            0        0.00000000
## 144     dx_af            0        0.00000000
## 145     dx_ar            0        0.00000000
## 146     dx_hl            0        0.00000000
## 147     dx_dm            0        0.00000000
## 148    dx_ckd            0        0.00000000
## 149     dx_li            0        0.00000000
## 150    dx_ane            0        0.00000000
## 151    dx_hua            0        0.00000000
## 152    dx_ulc            0        0.00000000
## 153   dx_etc1            0        0.00000000
## 154       icd            0        0.00000000
## 155   icd_all            0        0.00000000
## 156    cohort            0        0.00000000
## 157  last_mod            0        0.00000000
## 158   baseday            0        0.00000000
## 159       uac            0        0.00000000
## 160      ubil            0        0.00000000
## 161    dx_cvd            0        0.00000000
## 162  last_std            0        0.00000000
## 163  last_mid            0        0.00000000
## 164  last_cvd            0        0.00000000
## 165      pycv            0        0.00000000
## 166      pymi            0        0.00000000
## 167      pyst            0        0.00000000
## 168      pymo            0        0.00000000
## 169        id            0        0.00000000

4.2 For follow-up

# Calculate missing data
dat2 %>%
  summarise_all(~sum(is.na(.))) %>%
  gather(key = "Variable", value = "MissingCount") %>%
  mutate(MissingPercentage = MissingCount / nrow(dat2) * 100) %>%
  arrange(desc(MissingPercentage))
##     Variable MissingCount MissingPercentage
## 1       lf11        46945      91.052794911
## 2        med        46187      89.582605997
## 3       exer        46185      89.578726871
## 4       ls02        46184      89.576787307
## 5       ls03        46184      89.576787307
## 6       ls01        46183      89.574847744
## 7     tobhon        41304      80.111718841
## 8      alc_q        40156      77.885100275
## 9       lf01        40123      77.821094689
## 10      lf02        40123      77.821094689
## 11      lf03        40123      77.821094689
## 12      lf04        40123      77.821094689
## 13      lf05        40123      77.821094689
## 14      lf06        40123      77.821094689
## 15      lf07        40123      77.821094689
## 16      lf08        40123      77.821094689
## 17      lf09        40123      77.821094689
## 18      lf10        40123      77.821094689
## 19      lh01        40123      77.821094689
## 20      lh02        40123      77.821094689
## 21      lh03        40123      77.821094689
## 22      lh04        40123      77.821094689
## 23      lh05        40123      77.821094689
## 24      hg01        40123      77.821094689
## 25      hg02        40123      77.821094689
## 26  hypliver        34798      67.492920594
## 27   sprenom        34797      67.490981031
## 28   etcphex        34772      67.442491951
## 29   brsound        34761      67.421156755
## 30     d_ecg        34761      67.421156755
## 31      ps08        34755      67.409519376
## 32      pa02        34753      67.405640250
## 33    stress        34752      67.403700687
## 34      pa03        34751      67.401761123
## 35     d_liv        34750      67.399821560
## 36     d_etc        34750      67.399821560
## 37      pa04        34750      67.399821560
## 38      ls04        34750      67.399821560
## 39      d_hl        34749      67.397881997
## 40      pa01        34749      67.397881997
## 41      ps01        34748      67.395942434
## 42      ps02        34748      67.395942434
## 43      ps03        34748      67.395942434
## 44      ps04        34748      67.395942434
## 45      ps05        34748      67.395942434
## 46      ps06        34748      67.395942434
## 47      ps07        34748      67.395942434
## 48      fd01        34747      67.394002871
## 49      fd02        34747      67.394002871
## 50      fd03        34747      67.394002871
## 51      fd04        34747      67.394002871
## 52      fd05        34747      67.394002871
## 53      fd06        34747      67.394002871
## 54      fd07        34747      67.394002871
## 55      fd08        34747      67.394002871
## 56      fd09        34747      67.394002871
## 57      fd10        34747      67.394002871
## 58      fd11        34747      67.394002871
## 59      fd12        34747      67.394002871
## 60      fd13        34747      67.394002871
## 61      fd14        34747      67.394002871
## 62      fd16        34747      67.394002871
## 63      fd15        34747      67.394002871
## 64     d_cdv        34746      67.392063307
## 65    arith2        34618      67.143799216
## 66    arith1        34592      67.093370573
## 67   alcfreq        32940      63.889212149
## 68     edema        28861      55.977733814
## 69    anemia        28856      55.968035998
## 70    hba1c2        14436      27.999534505
## 71      t_ua        11442      22.192482253
## 72     dx_st         7871      15.266302029
## 73     waist         2077       4.028472788
## 74   arhythm          882       1.710694752
## 75   fasting          277       0.537259009
## 76       tob          214       0.415066527
## 77       alc          190       0.368517010
## 78   casound           97       0.188137631
## 79    pulse2           52       0.100857287
## 80    height           23       0.044609954
## 81    pulse1           22       0.042670391
## 82      dbp2           20       0.038791264
## 83    weight           18       0.034912138
## 84      sbp2           18       0.034912138
## 85      dbp1           12       0.023274759
## 86      sbp1           10       0.019395632
## 87      t_st            7       0.013576942
## 88     d_ane            6       0.011637379
## 89     dx_ci            6       0.011637379
## 90      d_dm            5       0.009697816
## 91      t_ht            5       0.009697816
## 92      t_hl            5       0.009697816
## 93      t_dm            5       0.009697816
## 94      t_mi            5       0.009697816
## 95     dx_ht            5       0.009697816
## 96    dx_ich            5       0.009697816
## 97    dx_sah            5       0.009697816
## 98    dx_tia            5       0.009697816
## 99     dx_ap            5       0.009697816
## 100    dx_mi            5       0.009697816
## 101    dx_bd            5       0.009697816
## 102    dx_af            5       0.009697816
## 103    dx_ar            5       0.009697816
## 104    dx_hl            5       0.009697816
## 105    dx_dm            5       0.009697816
## 106    dx_th            5       0.009697816
## 107   dx_ost            5       0.009697816
## 108   dx_ckd            5       0.009697816
## 109   dx_cat            5       0.009697816
## 110    dx_li            5       0.009697816
## 111   dx_ane            5       0.009697816
## 112   dx_hua            5       0.009697816
## 113   dx_ulc            5       0.009697816
## 114   dx_can            5       0.009697816
## 115  dx_etc1            5       0.009697816
## 116  dx_etc2            5       0.009697816
## 117      age            4       0.007758253
## 118      sex            4       0.007758253
## 119   dx_chf            1       0.001939563
## 120       id            0       0.000000000
## 121  examday            0       0.000000000
## 122    u_sug            0       0.000000000
## 123    u_pro            0       0.000000000
## 124    u_occ            0       0.000000000
## 125   u_look            0       0.000000000
## 126     u_ph            0       0.000000000
## 127    u_ket            0       0.000000000
## 128    u_uro            0       0.000000000
## 129    u_bil            0       0.000000000
## 130      u_g            0       0.000000000
## 131    u_no2            0       0.000000000
## 132    u_wbc            0       0.000000000
## 133       tp            0       0.000000000
## 134  agratio            0       0.000000000
## 135      ztt            0       0.000000000
## 136       tc            0       0.000000000
## 137      hdl            0       0.000000000
## 138      ldl            0       0.000000000
## 139       tg            0       0.000000000
## 140      got            0       0.000000000
## 141      gpt            0       0.000000000
## 142     rgtp            0       0.000000000
## 143       ua            0       0.000000000
## 144    hba1c            0       0.000000000
## 145     crea            0       0.000000000
## 146       bs            0       0.000000000
## 147      wbc            0       0.000000000
## 148      rbc            0       0.000000000
## 149       hg            0       0.000000000
## 150      hct            0       0.000000000
## 151      plt            0       0.000000000

4.3 For new dataset

# Calculate missing data
dat3 %>%
  summarise_all(~sum(is.na(.))) %>%
  gather(key = "Variable", value = "MissingCount") %>%
  mutate(MissingPercentage = MissingCount / nrow(dat3) * 100) %>%
  arrange(desc(MissingPercentage))
##      Variable MissingCount MissingPercentage
## 1        uwbc        59230      1.000000e+02
## 2       d_stk        59230      1.000000e+02
## 3        d_mi        59230      1.000000e+02
## 4      f_lost        59230      1.000000e+02
## 5         ph8        59229      9.999831e+01
## 6     dth_dur        59224      9.998987e+01
## 7         ph7        59224      9.998987e+01
## 8         ph6        59209      9.996454e+01
## 9         paf        59175      9.990714e+01
## 10     dth_mi        59164      9.988857e+01
## 11        ph5        59132      9.983454e+01
## 12    mi_area        59130      9.983117e+01
## 13      mine5        59127      9.982610e+01
## 14     af_day        59070      9.972987e+01
## 15    int_knd        59070      9.972987e+01
## 16   int_oped        59066      9.972311e+01
## 17      mine4        59007      9.962350e+01
## 18   ci_scale        58963      9.954921e+01
## 19    ci_type        58954      9.953402e+01
## 20    ci_area        58954      9.953402e+01
## 21        ph4        58926      9.948675e+01
## 22   mi_onset        58925      9.948506e+01
## 23         ct        58898      9.943947e+01
## 24        mri        58897      9.943778e+01
## 25   st_onset        58792      9.926051e+01
## 26      mine3        58674      9.906129e+01
## 27     tobsty        58549      9.885024e+01
## 28     cp_dur        58493      9.875570e+01
## 29        ph3        58416      9.862570e+01
## 30      tobyr        58007      9.793517e+01
## 31   ido_date        57696      9.741010e+01
## 32     inslin        57507      9.709100e+01
## 33      mine2        57466      9.702178e+01
## 34      seiri        57382      9.687996e+01
## 35        ph2        57377      9.687152e+01
## 36    death_d        57359      9.684113e+01
## 37       vldl        57090      9.638697e+01
## 38    heikeiy        56684      9.570150e+01
## 39      memai        55404      9.354044e+01
## 40   nocturia        55404      9.354044e+01
## 41   cp_sleep        55402      9.353706e+01
## 42   alc_freq        55367      9.347797e+01
## 43     alcgou        55345      9.344082e+01
## 44        ph1        55308      9.337836e+01
## 45     heikei        55198      9.319264e+01
## 46   deth_moy        55048      9.293939e+01
## 47       lf11        54617      9.221172e+01
## 48        med        53859      9.093196e+01
## 49       ls02        53856      9.092690e+01
## 50       ls03        53856      9.092690e+01
## 51       ls01        53855      9.092521e+01
## 52   deth_fay        53676      9.062300e+01
## 53       wt20        53596      9.048793e+01
## 54    icd_stk        53262      8.992402e+01
## 55    icd_ihd        53180      8.978558e+01
## 56    icd_cvd        52952      8.940064e+01
## 57       hako        52936      8.937363e+01
## 58      elbow        52934      8.937025e+01
## 59         bf        52680      8.894141e+01
## 60    icd_can        52673      8.892960e+01
## 61    unhappy        52399      8.846699e+01
## 62     satisf        52351      8.838595e+01
## 63       drug        52316      8.832686e+01
## 64      sleep        52289      8.828128e+01
## 65      vital        52279      8.826439e+01
## 66     social        52279      8.826439e+01
## 67      dx_sd        52269      8.824751e+01
## 68     dx_pmi        52243      8.820361e+01
## 69     dx_int        52163      8.806855e+01
## 70      scott        52148      8.804322e+01
## 71    sheieht        52137      8.802465e+01
## 72    sheieas        52136      8.802296e+01
## 73       frct        52037      8.785582e+01
## 74      bruit        52035      8.785244e+01
## 75        alb        51978      8.775620e+01
## 76        glb        51978      8.775620e+01
## 77     dx_ihd        51970      8.774270e+01
## 78     kaidan        51946      8.770218e+01
## 79  dysphagia        51636      8.717879e+01
## 80        uph        51616      8.714503e+01
## 81       mahi        51616      8.714503e+01
## 82     murmur        51615      8.714334e+01
## 83   sft_kata        51611      8.713659e+01
## 84    deth_mo        51610      8.713490e+01
## 85    deth_fa        51601      8.711970e+01
## 86  chestpain        51599      8.711633e+01
## 87    sft_ude        51581      8.708594e+01
## 88        cpk        51578      8.708087e+01
## 89       fast        51577      8.707918e+01
## 90        hip        51576      8.707749e+01
## 91         ca        51563      8.705555e+01
## 92      mine1        51560      8.705048e+01
## 93        fh1        51559      8.704879e+01
## 94        fh2        51559      8.704879e+01
## 95        fh3        51559      8.704879e+01
## 96        fh4        51559      8.704879e+01
## 97        fh5        51559      8.704879e+01
## 98        fh6        51559      8.704879e+01
## 99        fh7        51559      8.704879e+01
## 100       fh8        51559      8.704879e+01
## 101       usg        51558      8.704710e+01
## 102       upr        51558      8.704710e+01
## 103       uob        51558      8.704710e+01
## 104       icd        51558      8.704710e+01
## 105   icd_all        51558      8.704710e+01
## 106    cohort        51558      8.704710e+01
## 107  last_mod        51558      8.704710e+01
## 108   baseday        51558      8.704710e+01
## 109       uac        51558      8.704710e+01
## 110      ubil        51558      8.704710e+01
## 111    dx_cvd        51558      8.704710e+01
## 112  last_std        51558      8.704710e+01
## 113  last_mid        51558      8.704710e+01
## 114  last_cvd        51558      8.704710e+01
## 115      pycv        51558      8.704710e+01
## 116      pymi        51558      8.704710e+01
## 117      pyst        51558      8.704710e+01
## 118      pymo        51558      8.704710e+01
## 119     alc_q        47828      8.074962e+01
## 120      lf01        47795      8.069391e+01
## 121      lf02        47795      8.069391e+01
## 122      lf03        47795      8.069391e+01
## 123      lf04        47795      8.069391e+01
## 124      lf05        47795      8.069391e+01
## 125      lf06        47795      8.069391e+01
## 126      lf07        47795      8.069391e+01
## 127      lf08        47795      8.069391e+01
## 128      lf09        47795      8.069391e+01
## 129      lf10        47795      8.069391e+01
## 130      lh01        47795      8.069391e+01
## 131      lh02        47795      8.069391e+01
## 132      lh03        47795      8.069391e+01
## 133      lh04        47795      8.069391e+01
## 134      lh05        47795      8.069391e+01
## 135      hg01        47795      8.069391e+01
## 136      hg02        47795      8.069391e+01
## 137    tobhon        46640      7.874388e+01
## 138      exer        46497      7.850245e+01
## 139  hypliver        42470      7.170353e+01
## 140   sprenom        42469      7.170184e+01
## 141   etcphex        42444      7.165963e+01
## 142   brsound        42433      7.164106e+01
## 143     d_ecg        42433      7.164106e+01
## 144      ps08        42427      7.163093e+01
## 145      pa02        42425      7.162755e+01
## 146      pa03        42423      7.162418e+01
## 147     d_liv        42422      7.162249e+01
## 148     d_etc        42422      7.162249e+01
## 149      pa04        42422      7.162249e+01
## 150      ls04        42422      7.162249e+01
## 151      d_hl        42421      7.162080e+01
## 152      pa01        42421      7.162080e+01
## 153      ps01        42420      7.161911e+01
## 154      ps02        42420      7.161911e+01
## 155      ps03        42420      7.161911e+01
## 156      ps04        42420      7.161911e+01
## 157      ps05        42420      7.161911e+01
## 158      ps06        42420      7.161911e+01
## 159      ps07        42420      7.161911e+01
## 160      fd01        42419      7.161742e+01
## 161      fd02        42419      7.161742e+01
## 162      fd03        42419      7.161742e+01
## 163      fd04        42419      7.161742e+01
## 164      fd05        42419      7.161742e+01
## 165      fd06        42419      7.161742e+01
## 166      fd07        42419      7.161742e+01
## 167      fd08        42419      7.161742e+01
## 168      fd09        42419      7.161742e+01
## 169      fd10        42419      7.161742e+01
## 170      fd11        42419      7.161742e+01
## 171      fd12        42419      7.161742e+01
## 172      fd13        42419      7.161742e+01
## 173      fd14        42419      7.161742e+01
## 174      fd16        42419      7.161742e+01
## 175      fd15        42419      7.161742e+01
## 176     d_cdv        42418      7.161574e+01
## 177   alcfreq        40612      6.856660e+01
## 178     edema        36533      6.167989e+01
## 179    anemia        36528      6.167145e+01
## 180    stress        35557      6.003208e+01
## 181    arith2        35337      5.966064e+01
## 182    arith1        35235      5.948843e+01
## 183    hba1c2        22108      3.732568e+01
## 184      t_ua        19114      3.227081e+01
## 185   arhythm         8554      1.444201e+01
## 186     dx_st         8150      1.375992e+01
## 187   fasting         7949      1.342056e+01
## 188   casound         7769      1.311666e+01
## 189      t_st         7679      1.296471e+01
## 190     d_ane         7678      1.296303e+01
## 191      d_dm         7677      1.296134e+01
## 192      t_ht         7677      1.296134e+01
## 193      t_hl         7677      1.296134e+01
## 194      t_dm         7677      1.296134e+01
## 195      t_mi         7677      1.296134e+01
## 196     dx_ap         7677      1.296134e+01
## 197     dx_th         7677      1.296134e+01
## 198    dx_ost         7677      1.296134e+01
## 199    dx_cat         7677      1.296134e+01
## 200    dx_can         7677      1.296134e+01
## 201   dx_etc2         7677      1.296134e+01
## 202   examday         7672      1.295290e+01
## 203     u_sug         7672      1.295290e+01
## 204     u_pro         7672      1.295290e+01
## 205     u_occ         7672      1.295290e+01
## 206    u_look         7672      1.295290e+01
## 207      u_ph         7672      1.295290e+01
## 208     u_ket         7672      1.295290e+01
## 209     u_uro         7672      1.295290e+01
## 210     u_bil         7672      1.295290e+01
## 211       u_g         7672      1.295290e+01
## 212     u_no2         7672      1.295290e+01
## 213     u_wbc         7672      1.295290e+01
## 214       plt         7672      1.295290e+01
## 215       ldl         5543      9.358433e+00
## 216     hba1c         4930      8.323485e+00
## 217     waist         2095      3.537059e+00
## 218    dx_sah          685      1.156509e+00
## 219    dx_ich          643      1.085599e+00
## 220     dx_mi          636      1.073780e+00
## 221     dx_ci          447      7.546851e-01
## 222       tob          344      5.807868e-01
## 223       alc          311      5.250718e-01
## 224    pulse2           75      1.266250e-01
## 225       hdl           35      5.909168e-02
## 226       rbc           27      4.558501e-02
## 227    height           26      4.389667e-02
## 228    pulse1           26      4.389667e-02
## 229       wbc           26      4.389667e-02
## 230        hg           26      4.389667e-02
## 231       hct           26      4.389667e-02
## 232      dbp2           22      3.714334e-02
## 233    weight           21      3.545501e-02
## 234      sbp2           20      3.376667e-02
## 235       ztt           20      3.376667e-02
## 236   agratio           17      2.870167e-02
## 237      dbp1           14      2.363667e-02
## 238      sbp1           11      1.857167e-02
## 239       gpt            5      8.441668e-03
## 240     dx_ht            5      8.441668e-03
## 241    dx_tia            5      8.441668e-03
## 242     dx_bd            5      8.441668e-03
## 243     dx_af            5      8.441668e-03
## 244     dx_ar            5      8.441668e-03
## 245     dx_hl            5      8.441668e-03
## 246     dx_dm            5      8.441668e-03
## 247    dx_ckd            5      8.441668e-03
## 248     dx_li            5      8.441668e-03
## 249    dx_ane            5      8.441668e-03
## 250    dx_hua            5      8.441668e-03
## 251    dx_ulc            5      8.441668e-03
## 252   dx_etc1            5      8.441668e-03
## 253       age            4      6.753334e-03
## 254       sex            4      6.753334e-03
## 255        ua            3      5.065001e-03
## 256        tp            2      3.376667e-03
## 257        tc            2      3.376667e-03
## 258        tg            2      3.376667e-03
## 259       got            2      3.376667e-03
## 260      rgtp            2      3.376667e-03
## 261      crea            2      3.376667e-03
## 262        bs            2      3.376667e-03
## 263    dx_chf            1      1.688334e-03
## 264        id            0      0.000000e+00