data<-read.csv("C:/R_work/ECMO903.csv")
head(data)
  Unit.number   Name Sex Age     Date Year Holiday Insertion.time
1   `08243022 이영순   1  49 20090511 2009       1              1
2    13034861 윤기정   1  77 20090518 2009       1             NA
3    14186610 임지현   1  33 20090722 2009       1              1
4    14340394 이수맹   0  67 20091007 2009       1              1
5    10870744 김창길   0  56 20100212 2010       1              0
6    11071670 이영익   0  67 20100217 2010       1              0
  ECMO.duration Death Height Weight   BMI   Diagnosis Smoking CAD CVA PAD
1             8     1     NA     NA    NA Myocarditis       2   1   1   1
2            24     0     NA     NA    NA         ACS       2   1   1   1
3             8     2     NA     NA    NA        ARDS       2   1   1   1
4            15     0     NA     NA    NA         ACS       0   0   1   1
5             5     0   1.65     68 29.76         ACS       2   1   1   1
6             2     0   1.70     65 26.37        ARDS       2   1   1   1
  DM HBP Indication Mode CA OHCA CPR.duration CA.ECMO.Time EKG Weaning
1  1   1          0    0  1   NA           NA           NA   2       0
2  1   0          0    0  1   NA           NA           NA   2       0
3  1   1          1    1  1   NA           NA           NA   2       0
4  0   0          0    0  0    1           37           37   0       0
5  1   0          0    0  0    1           25           25   0       1
6  1   1          1    1  0    1          210          210   2       1
  Vent Intu.day HD ICU.stay CRRT CRRT.duration CAG PCI GCS Bleeding
1    0        8 30       11    0             8   1   1  15        1
2    0       27 31       31    1            NA   0   1   5        1
3    0       27 28       28    1            NA   1   1  15        0
4    0       17 17       17    0            13   0   0   3        1
5    0        5  5        5    0             3   0   0   3        1
6    0        3 43        3    1            NA   1   1  12        1
  Tamponade Pulmo AKI Liver.failure MOF Sepsis Stroke Leg.ischemia   WBC
1         1     1   0             1   1      1      1            1 16400
2         1     1   0             1   1      1      1            1  5230
3         1     0   1             1   0      0      1            1 18030
4         1     0   0             1   0      1      1            1  9670
5         1     1   0             1   1      1      1            1 10930
6         1     1   0             1   1      0      1            1 17950
   Hct PLT Sodium   K HCO3  BUN  Cr Glucose Albumin Total..bilirubin   CRP
1 43.0 248  143.2 3.9 15.9 20.1 2.3     184     4.7              1.1  0.31
2 30.1 154  143.2 3.1 18.0  6.3 1.0     153     4.2              0.5  1.03
3 37.8  55  148.9 3.4 23.3 22.3 1.3     150     3.1              1.1  8.22
4 36.0  54  144.8 3.4  9.0 20.1 1.5      56     2.2              0.5 13.34
5 44.8 293  136.8 3.2 13.4 15.9 1.5     200     2.5              0.6  0.07
6 22.3  70  143.2 6.5 10.9 48.9 2.5     488     1.6              0.3 23.27
  D.dimer  CK.MB Troponin PRC FFP  PC Cryo
1    2.29  30.95    3.050   5   3  72    0
2      NA   3.97    0.298  22   7  61    0
3   11.28 219.50    1.300   4  19  40    0
4    1.03   8.91    0.088  35   9 136    0
5      NA   3.47    0.020  10   0  18    0
6    0.44   4.38    0.036  24   5  42    0
str(data)
'data.frame':   78 obs. of  65 variables:
 $ Unit.number     : Factor w/ 78 levels "`08243022","10279578",..: 1 9 14 15 3 4 68 5 11 17 ...
 $ Name            : Factor w/ 76 levels "강신윤","강신훈",..: 44 37 56 42 13 45 38 30 17 62 ...
 $ Sex             : int  1 1 1 0 0 0 0 0 0 1 ...
 $ Age             : int  49 77 33 67 56 67 76 48 67 69 ...
 $ Date            : int  20090511 20090518 20090722 20091007 20100212 20100217 20100223 20100317 20100323 20100408 ...
 $ Year            : int  2009 2009 2009 2009 2010 2010 2010 2010 2010 2010 ...
 $ Holiday         : int  1 1 1 1 1 1 1 1 1 1 ...
 $ Insertion.time  : int  1 NA 1 1 0 0 0 1 1 0 ...
 $ ECMO.duration   : num  8 24 8 15 5 2 5 1 14 6 ...
 $ Death           : int  1 0 2 0 0 0 0 0 0 0 ...
 $ Height          : num  NA NA NA NA 1.65 1.7 1.64 1.75 1.65 1.55 ...
 $ Weight          : num  NA NA NA NA 68 65 68 75 60 53 ...
 $ BMI             : num  NA NA NA NA 29.8 ...
 $ Diagnosis       : Factor w/ 9 levels "ACS","ARDS","Drug Intoxication",..: 5 1 2 1 1 2 1 1 1 1 ...
 $ Smoking         : int  2 2 2 0 2 2 1 0 2 2 ...
 $ CAD             : int  1 1 1 0 1 1 1 1 1 1 ...
 $ CVA             : int  1 1 1 1 1 1 1 1 1 1 ...
 $ PAD             : int  1 1 1 1 1 1 1 1 1 1 ...
 $ DM              : int  1 1 1 0 1 1 0 0 0 0 ...
 $ HBP             : int  1 0 1 0 0 1 0 0 0 1 ...
 $ Indication      : int  0 0 1 0 0 1 0 0 0 0 ...
 $ Mode            : int  0 0 1 0 0 1 0 0 0 0 ...
 $ CA              : int  1 1 1 0 0 0 1 0 1 0 ...
 $ OHCA            : int  NA NA NA 1 1 1 NA 0 NA 0 ...
 $ CPR.duration    : int  NA NA NA 37 25 210 NA 30 NA 20 ...
 $ CA.ECMO.Time    : int  NA NA NA 37 25 210 NA 240 NA 540 ...
 $ EKG             : int  2 2 2 0 0 2 2 2 2 2 ...
 $ Weaning         : int  0 0 0 0 1 1 0 1 0 1 ...
 $ Vent            : int  0 0 0 0 0 0 0 0 0 0 ...
 $ Intu.day        : Factor w/ 29 levels "","0","1","10",..: 27 17 17 10 24 19 28 12 11 25 ...
 $ HD              : int  30 31 28 17 5 43 12 2 20 6 ...
 $ ICU.stay        : int  11 31 28 17 5 3 12 2 20 6 ...
 $ CRRT            : int  0 1 1 0 0 1 1 1 1 0 ...
 $ CRRT.duration   : int  8 NA NA 13 3 NA NA NA NA 5 ...
 $ CAG             : int  1 0 1 0 0 1 0 1 0 1 ...
 $ PCI             : int  1 1 1 0 0 1 0 1 1 1 ...
 $ GCS             : int  15 5 15 3 3 12 4 3 3 3 ...
 $ Bleeding        : int  1 1 0 1 1 1 1 0 0 1 ...
 $ Tamponade       : int  1 1 1 1 1 1 1 1 1 1 ...
 $ Pulmo           : int  1 1 0 0 1 1 1 1 1 1 ...
 $ AKI             : int  0 0 1 0 0 0 1 1 1 0 ...
 $ Liver.failure   : int  1 1 1 1 1 1 1 1 1 1 ...
 $ MOF             : int  1 1 0 0 1 1 0 0 0 1 ...
 $ Sepsis          : int  1 1 0 1 1 0 1 1 0 1 ...
 $ Stroke          : int  1 1 1 1 1 1 1 1 1 1 ...
 $ Leg.ischemia    : int  1 1 1 1 1 1 1 1 1 1 ...
 $ WBC             : int  16400 5230 18030 9670 10930 17950 6790 20410 18230 21500 ...
 $ Hct             : num  43 30.1 37.8 36 44.8 22.3 40.1 32.1 26.5 30.2 ...
 $ PLT             : int  248 154 55 54 293 70 173 203 334 235 ...
 $ Sodium          : num  143 143 149 145 137 ...
 $ K               : num  3.9 3.1 3.4 3.4 3.2 6.5 3.8 6.6 4.1 5.4 ...
 $ HCO3            : num  15.9 18 23.3 9 13.4 10.9 14.9 10.5 12.7 12.5 ...
 $ BUN             : num  20.1 6.3 22.3 20.1 15.9 48.9 15.8 54.6 49.9 40.7 ...
 $ Cr              : num  2.3 1 1.3 1.5 1.5 2.5 0.9 3.9 1.3 1.6 ...
 $ Glucose         : int  184 153 150 56 200 488 178 256 230 208 ...
 $ Albumin         : num  4.7 4.2 3.1 2.2 2.5 1.6 3.8 3 4.1 2.8 ...
 $ Total..bilirubin: num  1.1 0.5 1.1 0.5 0.6 0.3 1.1 0.6 0.5 1 ...
 $ CRP             : num  0.31 1.03 8.22 13.34 0.07 ...
 $ D.dimer         : num  2.29 NA 11.28 1.03 NA ...
 $ CK.MB           : num  30.95 3.97 219.5 8.91 3.47 ...
 $ Troponin        : num  3.05 0.298 1.3 0.088 0.02 ...
 $ PRC             : int  5 22 4 35 10 24 14 9 42 10 ...
 $ FFP             : int  3 7 19 9 0 5 0 0 5 0 ...
 $ PC              : int  72 61 40 136 18 42 46 0 170 54 ...
 $ Cryo            : int  0 0 0 0 0 0 0 0 0 0 ...
summary(data)
    Unit.number      Name         Sex            Age      
 `08243022: 1   이영순 : 2   Min.   :0.00   Min.   : 2.0  
 10279578 : 1   이영화 : 2   1st Qu.:0.00   1st Qu.:49.2  
 10870744 : 1   강신윤 : 1   Median :0.00   Median :63.0  
 11071670 : 1   강신훈 : 1   Mean   :0.41   Mean   :59.6  
 11144868 : 1   강정애 : 1   3rd Qu.:1.00   3rd Qu.:72.8  
 11620651 : 1   김기수 : 1   Max.   :1.00   Max.   :86.0  
 (Other)  :72   (Other):70                                
      Date               Year         Holiday      Insertion.time 
 Min.   :20090511   Min.   :2009   Min.   :0.000   Min.   :0.000  
 1st Qu.:20110644   1st Qu.:2011   1st Qu.:1.000   1st Qu.:0.000  
 Median :20121022   Median :2012   Median :1.000   Median :0.000  
 Mean   :20120930   Mean   :2012   Mean   :0.859   Mean   :0.416  
 3rd Qu.:20131028   3rd Qu.:2013   3rd Qu.:1.000   3rd Qu.:1.000  
 Max.   :20140804   Max.   :2014   Max.   :1.000   Max.   :1.000  
                                                   NA's   :1      
 ECMO.duration      Death           Height         Weight    
 Min.   : 0.2   Min.   :0.000   Min.   :0.91   Min.   :11.5  
 1st Qu.: 2.0   1st Qu.:0.000   1st Qu.:1.60   1st Qu.:55.0  
 Median : 6.0   Median :0.000   Median :1.65   Median :60.0  
 Mean   : 7.4   Mean   :0.372   Mean   :1.64   Mean   :60.8  
 3rd Qu.:10.8   3rd Qu.:1.000   3rd Qu.:1.70   3rd Qu.:68.0  
 Max.   :24.0   Max.   :2.000   Max.   :1.83   Max.   :80.0  
                                NA's   :11     NA's   :9     
      BMI              Diagnosis     Smoking          CAD       
 Min.   :12.5   ACS         :34   Min.   :0.00   Min.   :0.000  
 1st Qu.:24.6   ARDS        :14   1st Qu.:0.00   1st Qu.:1.000  
 Median :27.1   PTE         : 8   Median :2.00   Median :1.000  
 Mean   :26.8   Myocarditis : 7   Mean   :1.32   Mean   :0.795  
 3rd Qu.:29.1   Pump failure: 7   3rd Qu.:2.00   3rd Qu.:1.000  
 Max.   :33.5   VF, VT      : 4   Max.   :2.00   Max.   :1.000  
 NA's   :11     (Other)     : 4                                 
      CVA             PAD          DM             HBP      
 Min.   :0.000   Min.   :1   Min.   :0.000   Min.   :0.00  
 1st Qu.:1.000   1st Qu.:1   1st Qu.:1.000   1st Qu.:0.00  
 Median :1.000   Median :1   Median :1.000   Median :1.00  
 Mean   :0.949   Mean   :1   Mean   :0.756   Mean   :0.59  
 3rd Qu.:1.000   3rd Qu.:1   3rd Qu.:1.000   3rd Qu.:1.00  
 Max.   :1.000   Max.   :1   Max.   :1.000   Max.   :1.00  
                                                           
   Indication         Mode             CA             OHCA     
 Min.   :0.000   Min.   :0.000   Min.   :0.000   Min.   :0.00  
 1st Qu.:0.000   1st Qu.:0.000   1st Qu.:0.000   1st Qu.:1.00  
 Median :0.000   Median :0.000   Median :0.000   Median :1.00  
 Mean   :0.192   Mean   :0.192   Mean   :0.487   Mean   :0.88  
 3rd Qu.:0.000   3rd Qu.:0.000   3rd Qu.:1.000   3rd Qu.:1.00  
 Max.   :1.000   Max.   :1.000   Max.   :1.000   Max.   :1.00  
                                                 NA's   :38    
  CPR.duration    CA.ECMO.Time         EKG          Weaning     
 Min.   :  2.0   Min.   :   0.0   Min.   :0.00   Min.   :0.000  
 1st Qu.: 10.8   1st Qu.:  19.5   1st Qu.:1.00   1st Qu.:0.000  
 Median : 25.0   Median :  39.0   Median :2.00   Median :1.000  
 Mean   : 30.6   Mean   :  92.3   Mean   :1.53   Mean   :0.603  
 3rd Qu.: 40.2   3rd Qu.:  48.5   3rd Qu.:2.00   3rd Qu.:1.000  
 Max.   :210.0   Max.   :1080.0   Max.   :2.00   Max.   :2.000  
 NA's   :38      NA's   :38                                     
      Vent           Intu.day        HD           ICU.stay   
 Min.   :0.0000   2      : 7   Min.   : 0.00   Min.   : 0.0  
 1st Qu.:0.0000   1      : 6   1st Qu.: 4.25   1st Qu.: 4.0  
 Median :0.0000   3      : 6   Median :18.50   Median :10.0  
 Mean   :0.0385   4      : 6   Mean   :22.59   Mean   :12.9  
 3rd Qu.:0.0000   6      : 6   3rd Qu.:34.00   3rd Qu.:19.0  
 Max.   :1.0000   10     : 4   Max.   :86.00   Max.   :44.0  
                  (Other):43                                 
      CRRT       CRRT.duration        CAG           PCI       
 Min.   :0.000   Min.   : 1.00   Min.   :0.0   Min.   :0.000  
 1st Qu.:0.000   1st Qu.: 3.00   1st Qu.:0.0   1st Qu.:0.000  
 Median :1.000   Median : 5.00   Median :0.5   Median :1.000  
 Mean   :0.603   Mean   : 6.39   Mean   :0.5   Mean   :0.718  
 3rd Qu.:1.000   3rd Qu.: 8.00   3rd Qu.:1.0   3rd Qu.:1.000  
 Max.   :1.000   Max.   :20.00   Max.   :1.0   Max.   :1.000  
                 NA's   :47                                   
      GCS           Bleeding       Tamponade         Pulmo      
 Min.   : 3.00   Min.   :0.000   Min.   :0.000   Min.   :0.000  
 1st Qu.: 4.00   1st Qu.:0.000   1st Qu.:1.000   1st Qu.:1.000  
 Median : 7.00   Median :1.000   Median :1.000   Median :1.000  
 Mean   : 8.26   Mean   :0.667   Mean   :0.949   Mean   :0.872  
 3rd Qu.:13.00   3rd Qu.:1.000   3rd Qu.:1.000   3rd Qu.:1.000  
 Max.   :15.00   Max.   :1.000   Max.   :1.000   Max.   :1.000  
                                                                
      AKI        Liver.failure        MOF            Sepsis     
 Min.   :0.000   Min.   :0.000   Min.   :0.000   Min.   :0.000  
 1st Qu.:0.000   1st Qu.:1.000   1st Qu.:1.000   1st Qu.:1.000  
 Median :1.000   Median :1.000   Median :1.000   Median :1.000  
 Mean   :0.577   Mean   :0.949   Mean   :0.769   Mean   :0.859  
 3rd Qu.:1.000   3rd Qu.:1.000   3rd Qu.:1.000   3rd Qu.:1.000  
 Max.   :1.000   Max.   :1.000   Max.   :1.000   Max.   :1.000  
                                                                
     Stroke       Leg.ischemia        WBC             Hct      
 Min.   :0.000   Min.   :0.000   Min.   :  610   Min.   :18.1  
 1st Qu.:1.000   1st Qu.:1.000   1st Qu.: 8808   1st Qu.:30.7  
 Median :1.000   Median :1.000   Median :12330   Median :35.1  
 Mean   :0.962   Mean   :0.949   Mean   :13493   Mean   :34.8  
 3rd Qu.:1.000   3rd Qu.:1.000   3rd Qu.:16602   3rd Qu.:40.0  
 Max.   :1.000   Max.   :1.000   Max.   :48290   Max.   :53.8  
                                                               
      PLT            Sodium          K             HCO3     
 Min.   : 17.0   Min.   :128   Min.   :2.40   Min.   : 5.9  
 1st Qu.: 97.5   1st Qu.:139   1st Qu.:3.62   1st Qu.:10.9  
 Median :144.5   Median :141   Median :4.00   Median :13.9  
 Mean   :162.1   Mean   :141   Mean   :4.20   Mean   :14.8  
 3rd Qu.:215.5   3rd Qu.:144   3rd Qu.:4.50   3rd Qu.:18.1  
 Max.   :414.0   Max.   :157   Max.   :7.60   Max.   :32.3  
                                              NA's   :2     
      BUN             Cr          Glucose       Albumin    
 Min.   : 6.3   Min.   :0.40   Min.   : 50   Min.   :1.40  
 1st Qu.:15.6   1st Qu.:0.90   1st Qu.:136   1st Qu.:2.40  
 Median :20.4   Median :1.20   Median :182   Median :2.95  
 Mean   :25.7   Mean   :1.44   Mean   :203   Mean   :3.03  
 3rd Qu.:32.2   3rd Qu.:1.50   3rd Qu.:247   3rd Qu.:3.77  
 Max.   :88.9   Max.   :8.30   Max.   :488   Max.   :4.70  
                               NA's   :2                   
 Total..bilirubin      CRP           D.dimer          CK.MB       
 Min.   :0.30     Min.   : 0.02   Min.   : 0.05   Min.   :  0.70  
 1st Qu.:0.50     1st Qu.: 0.43   1st Qu.: 1.21   1st Qu.:  3.01  
 Median :0.90     Median : 2.52   Median : 4.84   Median : 16.52  
 Mean   :1.24     Mean   : 7.05   Mean   : 9.84   Mean   : 42.42  
 3rd Qu.:1.40     3rd Qu.:13.34   3rd Qu.:15.34   3rd Qu.: 40.55  
 Max.   :7.10     Max.   :30.25   Max.   :45.00   Max.   :300.00  
 NA's   :1        NA's   :1       NA's   :26      NA's   :9       
    Troponin           PRC             FFP           PC       
 Min.   : 0.010   Min.   :  0.0   Min.   : 0   Min.   :  0.0  
 1st Qu.: 0.034   1st Qu.:  8.0   1st Qu.: 0   1st Qu.: 10.0  
 Median : 0.232   Median : 18.0   Median : 7   Median : 32.0  
 Mean   : 1.430   Mean   : 23.9   Mean   :11   Mean   : 47.1  
 3rd Qu.: 1.225   3rd Qu.: 31.0   3rd Qu.:17   3rd Qu.: 60.0  
 Max.   :17.190   Max.   :150.0   Max.   :72   Max.   :240.0  
 NA's   :12       NA's   :1       NA's   :1    NA's   :1      
      Cryo     
 Min.   : 0.0  
 1st Qu.: 0.0  
 Median : 0.0  
 Mean   : 1.7  
 3rd Qu.: 0.0  
 Max.   :20.0  
 NA's   :1     
data1=data[data$Death<2,-c(1,2,5)]
data1$Survive=ifelse(data1$Death==0,0,1)
table(data1$Death)

 0  1 
52 23 
library(compareGroups)
data1$Sex=ifelse(data1$Sex==1,"Female","Male")
data1$Holiday=ifelse(data1$Holiday==1,"Not Holiday","Holiday")
data1$Insertion.time=ifelse(data1$Insertion.time==1,"Emergency","Elective")
data1$Death=ifelse(data1$Death==1,"Survivor","Non-survivor")
data1$Smoking=ifelse(data1$Smoking==1,"No","Yes")
data1$CAD=ifelse(data1$CAD==1,"No","Yes")
data1$CVA=ifelse(data1$CVA==1,"No","Yes")
data1$PAD=ifelse(data1$PAD==1,"No","Yes")
data1$DM=ifelse(data1$DM==1,"No","Yes")
data1$HBP=ifelse(data1$HBP==1,"No","Yes")
data1$Indication=ifelse(data1$Indication==1,"Respiratory","Cardiology")
data1$Mode=ifelse(data1$Mode==1,"VV","VA")
data1$CA=ifelse(data1$CA==1,"No","Yes")
data1$OHCA=ifelse(data1$OHCA==1,"IHCA","OHCA")
data1$EKG=ifelse(data1$EKG==0,"Asystole",
                 ifelse(data1$EKG==1,"PEA","ETC"))
data1$Weaning=ifelse(data1$Weaning==1,"Fail","Success")
data1$Vent=ifelse(data1$Vent==1,"No","Yes")
data1$CRRT=ifelse(data1$CRRT==1,"No","Yes")
data1$CAG=ifelse(data1$CAG==1,"No","Yes")
data1$PCI=ifelse(data1$PCI==1,"No","Yes")
data1$Bleeding=ifelse(data1$Bleeding==1,"No","Yes")
data1$Tamponade=ifelse(data1$Tamponade==1,"No","Yes")
data1$Pulmo=ifelse(data1$Pulmo==1,"No","Yes")
data1$AKI=ifelse(data1$AKI==1,"No","Yes")
data1$Liver.failure=ifelse(data1$Liver.failure==1,"No","Yes")
data1$MOF=ifelse(data1$MOF==1,"No","Yes")
data1$Sepsis=ifelse(data1$Sepsis==1,"No","Yes")
data1$Stroke=ifelse(data1$Stroke==1,"No","Yes")
data1$Leg.ischemia=ifelse(data1$Leg.ischemia==1,"No","Yes")
mosaicplot(~Year+Death,data=data1,color=TRUE,xlab="Year",ylab="Survival",main="Distribution of ECMO")

plot of chunk unnamed-chunk-1

res=compareGroups(Survive~.,data=data1)
res


-------- Summary of results by groups of 'Survive'---------


   var              N  p.value  method            selection
1  Sex              75 0.878    categorical       ALL      
2  Age              75 0.018**  continuous normal ALL      
3  Year             75 0.015**  continuous normal ALL      
4  Holiday          75 1.000    categorical       ALL      
5  Insertion.time   74 0.673    categorical       ALL      
6  ECMO.duration    75 0.241    continuous normal ALL      
7  Death            75 <0.001** categorical       ALL      
8  Height           65 0.708    continuous normal ALL      
9  Weight           67 0.980    continuous normal ALL      
10 BMI              65 0.694    continuous normal ALL      
11 Diagnosis        75 0.237    categorical       ALL      
12 Smoking          75 0.728    categorical       ALL      
13 CAD              75 0.369    categorical       ALL      
14 CVA              75 1.000    categorical       ALL      
15 PAD              75 .        categorical       ALL      
16 DM               75 0.180    categorical       ALL      
17 HBP              75 0.126    categorical       ALL      
18 Indication       75 0.750    categorical       ALL      
19 Mode             75 0.750    categorical       ALL      
20 CA               75 0.083*   categorical       ALL      
21 OHCA             39 1.000    categorical       ALL      
22 CPR.duration     39 0.055*   continuous normal ALL      
23 CA.ECMO.Time     39 0.076*   continuous normal ALL      
24 EKG              75 0.026**  categorical       ALL      
25 Weaning          75 <0.001** categorical       ALL      
26 Vent             75 0.026**  categorical       ALL      
27 HD               75 <0.001** continuous normal ALL      
28 ICU.stay         75 0.001**  continuous normal ALL      
29 CRRT             75 0.218    categorical       ALL      
30 CRRT.duration    29 0.175    continuous normal ALL      
31 CAG              75 0.281    categorical       ALL      
32 PCI              75 0.600    categorical       ALL      
33 GCS              75 0.001**  continuous normal ALL      
34 Bleeding         75 0.535    categorical       ALL      
35 Tamponade        75 1.000    categorical       ALL      
36 Pulmo            75 0.713    categorical       ALL      
37 AKI              75 0.308    categorical       ALL      
38 Liver.failure    75 0.306    categorical       ALL      
39 MOF              75 0.026**  categorical       ALL      
40 Sepsis           75 1.000    categorical       ALL      
41 Stroke           75 0.091*   categorical       ALL      
42 Leg.ischemia     75 0.582    categorical       ALL      
43 WBC              75 0.558    continuous normal ALL      
44 Hct              75 0.583    continuous normal ALL      
45 PLT              75 0.432    continuous normal ALL      
46 Sodium           75 0.404    continuous normal ALL      
47 K                75 0.015**  continuous normal ALL      
48 HCO3             73 0.042**  continuous normal ALL      
49 BUN              75 0.014**  continuous normal ALL      
50 Cr               75 0.003**  continuous normal ALL      
51 Glucose          73 0.437    continuous normal ALL      
52 Albumin          75 0.668    continuous normal ALL      
53 Total..bilirubin 74 0.499    continuous normal ALL      
54 CRP              74 0.187    continuous normal ALL      
55 D.dimer          49 0.796    continuous normal ALL      
56 CK.MB            66 0.784    continuous normal ALL      
57 Troponin         64 0.522    continuous normal ALL      
58 PRC              74 0.187    continuous normal ALL      
59 FFP              74 0.338    continuous normal ALL      
60 PC               74 0.074*   continuous normal ALL      
61 Cryo             74 0.753    continuous normal ALL      
-----
Signif. codes:  0 '**' 0.05 '*' 0.1 ' ' 1 
createTable(res)

--------Summary descriptives table by 'Survive'---------

_________________________________________________________ 
                           0            1       p.overall 
                          N=52         N=23               
<U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF> 
Sex:                                              0.878   
    Female             20 (38.5%)   10 (43.5%)            
    Male               32 (61.5%)   13 (56.5%)            
Age                   63.8 (13.6)  52.5 (19.6)    0.018   
Year                  2012 (1.40)  2013 (1.43)    0.015   
Holiday:                                          1.000   
    Holiday            8 (15.4%)    3 (13.0%)             
    Not Holiday        44 (84.6%)   20 (87.0%)            
Insertion.time:                                   0.673   
    Elective           29 (56.9%)   15 (65.2%)            
    Emergency          22 (43.1%)   8 (34.8%)             
ECMO.duration         6.75 (6.51)  8.36 (4.84)    0.241   
Death:                                           <0.001   
    Non-survivor       52 (100%)    0 (0.00%)             
    Survivor           0 (0.00%)    23 (100%)             
Height                1.64 (0.08)  1.63 (0.19)    0.708   
Weight                60.6 (8.69)  60.5 (14.3)    0.980   
BMI                   26.7 (3.50)  26.3 (3.96)    0.694   
Diagnosis:                                        0.237   
    ACS                27 (51.9%)   6 (26.1%)             
    ARDS               8 (15.4%)    5 (21.7%)             
    Drug Intoxication  1 (1.92%)    0 (0.00%)             
    LA myxoma          1 (1.92%)    0 (0.00%)             
    Myocarditis        2 (3.85%)    4 (17.4%)             
    PTE                4 (7.69%)    4 (17.4%)             
    Pump failure       5 (9.62%)    2 (8.70%)             
    Trauma             1 (1.92%)    1 (4.35%)             
    VF, VT             3 (5.77%)    1 (4.35%)             
Smoking:                                          0.728   
    No                 7 (13.5%)    4 (17.4%)             
    Yes                45 (86.5%)   19 (82.6%)            
CAD:                                              0.369   
    No                 40 (76.9%)   20 (87.0%)            
    Yes                12 (23.1%)   3 (13.0%)             
CVA:                                              1.000   
    No                 49 (94.2%)   22 (95.7%)            
    Yes                3 (5.77%)    1 (4.35%)             
PAD: No                52 (100%)    23 (100%)       .     
DM:                                               0.180   
    No                 36 (69.2%)   20 (87.0%)            
    Yes                16 (30.8%)   3 (13.0%)             
HBP:                                              0.126   
    No                 27 (51.9%)   17 (73.9%)            
    Yes                25 (48.1%)   6 (26.1%)             
Indication:                                       0.750   
    Cardiology         43 (82.7%)   18 (78.3%)            
    Respiratory        9 (17.3%)    5 (21.7%)             
Mode:                                             0.750   
    VA                 43 (82.7%)   18 (78.3%)            
    VV                 9 (17.3%)    5 (21.7%)             
CA:                                               0.083   
    No                 21 (40.4%)   15 (65.2%)            
    Yes                31 (59.6%)   8 (34.8%)             
OHCA:                                             1.000   
    IHCA               27 (87.1%)   7 (87.5%)             
    OHCA               4 (12.9%)    1 (12.5%)             
CPR.duration          33.8 (36.7)  16.9 (15.0)    0.055   
CA.ECMO.Time           109 (221)   35.2 (16.3)    0.076   
EKG:                                              0.026   
    Asystole           15 (28.8%)   1 (4.35%)             
    ETC                33 (63.5%)   21 (91.3%)            
    PEA                4 (7.69%)    1 (4.35%)             
Weaning:                                         <0.001   
    Fail               44 (84.6%)   0 (0.00%)             
    Success            8 (15.4%)    23 (100%)             
Vent:                                             0.026   
    No                 0 (0.00%)    3 (13.0%)             
    Yes                52 (100%)    20 (87.0%)            
HD                    14.7 (15.6)  40.4 (18.1)   <0.001   
ICU.stay              9.75 (9.43)  18.9 (10.1)    0.001   
CRRT:                                             0.218   
    No                 29 (55.8%)   17 (73.9%)            
    Yes                23 (44.2%)   6 (26.1%)             
CRRT.duration         5.78 (4.65)  9.17 (5.00)    0.175   
CAG:                                              0.281   
    No                 23 (44.2%)   14 (60.9%)            
    Yes                29 (55.8%)   9 (39.1%)             
PCI:                                              0.600   
    No                 36 (69.2%)   18 (78.3%)            
    Yes                16 (30.8%)   5 (21.7%)             
GCS                   6.81 (4.24)  10.9 (4.50)    0.001   
Bleeding:                                         0.535   
    No                 33 (63.5%)   17 (73.9%)            
    Yes                19 (36.5%)   6 (26.1%)             
Tamponade:                                        1.000   
    No                 49 (94.2%)   22 (95.7%)            
    Yes                3 (5.77%)    1 (4.35%)             
Pulmo:                                            0.713   
    No                 45 (86.5%)   21 (91.3%)            
    Yes                7 (13.5%)    2 (8.70%)             
AKI:                                              0.308   
    No                 28 (53.8%)   16 (69.6%)            
    Yes                24 (46.2%)   7 (30.4%)             
Liver.failure:                                    0.306   
    No                 48 (92.3%)   23 (100%)             
    Yes                4 (7.69%)    0 (0.00%)             
MOF:                                              0.026   
    No                 36 (69.2%)   22 (95.7%)            
    Yes                16 (30.8%)   1 (4.35%)             
Sepsis:                                           1.000   
    No                 46 (88.5%)   20 (87.0%)            
    Yes                6 (11.5%)    3 (13.0%)             
Stroke:                                           0.091   
    No                 52 (100%)    21 (91.3%)            
    Yes                0 (0.00%)    2 (8.70%)             
Leg.ischemia:                                     0.582   
    No                 50 (96.2%)   21 (91.3%)            
    Yes                2 (3.85%)    2 (8.70%)             
WBC                   13751 (8254) 12786 (5610)   0.558   
Hct                   34.4 (7.31)  35.5 (7.86)    0.583   
PLT                    167 (100)    151 (77.1)    0.432   
Sodium                 142 (5.23)   141 (5.78)    0.404   
K                     4.36 (1.07)  3.93 (0.44)    0.015   
HCO3                  13.7 (5.02)  17.0 (6.54)    0.042   
BUN                   28.7 (17.7)  20.4 (10.3)    0.014   
Cr                    1.62 (1.23)  1.04 (0.41)    0.003   
Glucose                207 (108)    192 (59.6)    0.437   
Albumin               2.99 (0.88)  3.07 (0.78)    0.668   
Total..bilirubin      1.19 (1.36)  1.42 (1.26)    0.499   
CRP                   6.12 (7.73)  9.41 (10.5)    0.187   
D.dimer               9.86 (12.2)  8.94 (11.4)    0.796   
CK.MB                 36.2 (47.2)  31.4 (71.1)    0.784   
Troponin              1.53 (3.48)  1.07 (2.12)    0.522   
PRC                   26.5 (27.6)  19.2 (18.3)    0.187   
FFP                   11.9 (14.2)  8.91 (11.3)    0.338   
PC                    53.6 (62.1)  33.4 (33.3)    0.074   
Cryo                  1.87 (4.83)  1.55 (3.54)    0.753   
<U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF><U+00AF> 
out=lm(Survive~Age+EKG+GCS+K+HCO3+BUN+Cr,data=data1)
anova(out)
Analysis of Variance Table

Response: Survive
          Df Sum Sq Mean Sq F value Pr(>F)   
Age        1   1.30   1.303    8.25 0.0055 **
EKG        2   0.95   0.475    3.01 0.0564 . 
GCS        1   1.63   1.626   10.29 0.0021 **
K          1   0.86   0.861    5.45 0.0227 * 
HCO3       1   0.21   0.213    1.35 0.2500   
BUN        1   0.29   0.286    1.81 0.1831   
Cr         1   0.02   0.017    0.11 0.7410   
Residuals 64  10.11   0.158                  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
confint(out)
               2.5 %    97.5 %
(Intercept) -0.09000  1.324465
Age         -0.01372 -0.001852
EKGETC      -0.06618  0.463505
EKGPEA      -0.28030  0.620177
GCS          0.01068  0.052544
K           -0.19917  0.067662
HCO3        -0.01018  0.029433
BUN         -0.01406  0.002781
Cr          -0.11749  0.164317
out2=lm(Survive~Age+GCS+K,data=data1)
anova(out2)
Analysis of Variance Table

Response: Survive
          Df Sum Sq Mean Sq F value  Pr(>F)    
Age        1   1.62   1.619    10.2 0.00205 ** 
GCS        1   2.31   2.306    14.6 0.00028 ***
K          1   0.81   0.806     5.1 0.02697 *  
Residuals 71  11.22   0.158                    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
confint(out2)
               2.5 %    97.5 %
(Intercept)  0.41150  1.619705
Age         -0.01448 -0.003158
GCS          0.01696  0.056448
K           -0.21032 -0.013095
summary(out)

Call:
lm(formula = Survive ~ Age + EKG + GCS + K + HCO3 + BUN + Cr, 
    data = data1)

Residuals:
    Min      1Q  Median      3Q     Max 
-0.6108 -0.3021 -0.0819  0.2735  0.7925 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)   
(Intercept)  0.61723    0.35402    1.74   0.0860 . 
Age         -0.00778    0.00297   -2.62   0.0109 * 
EKGETC       0.19866    0.13257    1.50   0.1389   
EKGPEA       0.16994    0.22538    0.75   0.4536   
GCS          0.03161    0.01048    3.02   0.0037 **
K           -0.06576    0.06678   -0.98   0.3285   
HCO3         0.00963    0.00991    0.97   0.3353   
BUN         -0.00564    0.00421   -1.34   0.1857   
Cr           0.02341    0.07053    0.33   0.7410   
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.397 on 64 degrees of freedom
  (2 observations deleted due to missingness)
Multiple R-squared:  0.342, Adjusted R-squared:  0.26 
F-statistic: 4.16 on 8 and 64 DF,  p-value: 0.000475
summary(out2)

Call:
lm(formula = Survive ~ Age + GCS + K, data = data1)

Residuals:
   Min     1Q Median     3Q    Max 
-0.642 -0.248 -0.122  0.317  0.856 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  1.01560    0.30297    3.35  0.00129 ** 
Age         -0.00882    0.00284   -3.11  0.00272 ** 
GCS          0.03670    0.00990    3.71  0.00041 ***
K           -0.11171    0.04946   -2.26  0.02697 *  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.397 on 71 degrees of freedom
Multiple R-squared:  0.297, Adjusted R-squared:  0.267 
F-statistic: 9.98 on 3 and 71 DF,  p-value: 1.43e-05
out=glm(Survive~Age+EKG+GCS+K+HCO3+BUN+Cr,family=binomial,data=data1)
summary(out)

Call:
glm(formula = Survive ~ Age + EKG + GCS + K + HCO3 + BUN + Cr, 
    family = binomial, data = data1)

Deviance Residuals: 
   Min      1Q  Median      3Q     Max  
-1.499  -0.665  -0.168   0.277   2.365  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)   
(Intercept)   3.0202     3.0746    0.98   0.3259   
Age          -0.0692     0.0277   -2.49   0.0126 * 
EKGETC        2.4616     1.3212    1.86   0.0624 . 
EKGPEA        1.7998     1.8365    0.98   0.3271   
GCS           0.2876     0.0962    2.99   0.0028 **
K            -0.5689     0.5481   -1.04   0.2992   
HCO3          0.0412     0.0724    0.57   0.5689   
BUN          -0.0535     0.0369   -1.45   0.1470   
Cr           -1.1460     0.8622   -1.33   0.1838   
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 89.355  on 72  degrees of freedom
Residual deviance: 52.918  on 64  degrees of freedom
  (2 observations deleted due to missingness)
AIC: 70.92

Number of Fisher Scoring iterations: 6
exp(coef(out))
(Intercept)         Age      EKGETC      EKGPEA         GCS           K 
    20.4957      0.9331     11.7234      6.0484      1.3332      0.5661 
       HCO3         BUN          Cr 
     1.0421      0.9479      0.3179 
confint(out)
              2.5 %   97.5 %
(Intercept) -3.0326  9.35893
Age         -0.1320 -0.02174
EKGETC       0.2193  5.74244
EKGPEA      -2.0306  5.69066
GCS          0.1182  0.50296
K           -1.7263  0.46948
HCO3        -0.1025  0.18679
BUN         -0.1324  0.01500
Cr          -3.0337  0.43597
exp(confint(out))
              2.5 %    97.5 %
(Intercept) 0.04819 1.160e+04
Age         0.87638 9.785e-01
EKGETC      1.24516 3.118e+02
EKGPEA      0.13126 2.961e+02
GCS         1.12548 1.654e+00
K           0.17795 1.599e+00
HCO3        0.90259 1.205e+00
BUN         0.87596 1.015e+00
Cr          0.04814 1.546e+00
out2=glm(Survive~Age+GCS,family=binomial,data=data1)
summary(out2)

Call:
glm(formula = Survive ~ Age + GCS, family = binomial, data = data1)

Deviance Residuals: 
   Min      1Q  Median      3Q     Max  
-1.695  -0.685  -0.435   0.772   1.962  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)   
(Intercept)   0.0160     1.1395    0.01   0.9888   
Age          -0.0454     0.0179   -2.54   0.0111 * 
GCS           0.2071     0.0642    3.23   0.0012 **
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 92.461  on 74  degrees of freedom
Residual deviance: 72.907  on 72  degrees of freedom
AIC: 78.91

Number of Fisher Scoring iterations: 4
exp(coef(out2))
(Intercept)         Age         GCS 
     1.0161      0.9556      1.2301 
confint(out2)
               2.5 %  97.5 %
(Intercept) -2.21905  2.3195
Age         -0.08318 -0.0121
GCS          0.08737  0.3418
exp(confint(out2))
             2.5 % 97.5 %
(Intercept) 0.1087 10.171
Age         0.9202  0.988
GCS         1.0913  1.408
m<-lm(Survive~K,data=data1)
m

Call:
lm(formula = Survive ~ K, data = data1)

Coefficients:
(Intercept)            K  
      0.754       -0.106  
plot(m)

plot of chunk unnamed-chunk-1plot of chunk unnamed-chunk-1plot of chunk unnamed-chunk-1plot of chunk unnamed-chunk-1