데이타를 읽어온다.

require(ztable)
## Loading required package: ztable
## Welcome to package ztable ver 0.1.5
require(moonBook)
## Loading required package: moonBook
options(ztable.type="html")
#options(ztable.type="latex")
MI=read.csv("inf.csv",stringsAsFactors=FALSE)
#MI$AVB
MI$CAVB=0
MI$CAVB[MI$AVB=="Advanced"]=1
MI$CAVB[MI$AVB=="Complete"]=1
#str(MI$CAVB)

표를 만든다.

res=mytable(Temporary1~.,data=MI)
ztable(res)
0 1 p
(N=274) (N=21)
Age 61.0 ± 12.6 66.5 ± 11.1 0.053
Sex 1.000
    F 74 (27.0%) 6 (28.6%)
    M 200 (73.0%) 15 (71.4%)
death 0.526
    0 263 (96.0%) 19 (90.5%)
    1 11 ( 4.0%) 2 ( 9.5%)
MI 0.907
    0 269 (98.2%) 20 (95.2%)
    1 5 ( 1.8%) 1 ( 4.8%)
TVR 1.000
    0 269 (98.2%) 21 (100.0%)
    1 5 ( 1.8%) 0 ( 0.0%)
MACE 0.421
    0 255 (93.1%) 18 (85.7%)
    1 19 ( 6.9%) 3 (14.3%)
Stroke 0.923
    0 266 (97.1%) 21 (100.0%)
    1 8 ( 2.9%) 0 ( 0.0%)
Fuday 397.2 ± 315.8 398.1 ± 322.4 0.990
Attempted 0.167
    1 201 (76.1%) 11 (55.0%)
    2 53 (20.1%) 7 (35.0%)
    3 7 ( 2.7%) 2 (10.0%)
    4 2 ( 0.8%) 0 ( 0.0%)
    5 1 ( 0.4%) 0 ( 0.0%)
SuccessfullyDilated 1.3 ± 0.6 1.5 ± 0.7 0.080
CASSNo 0.832
    1st OM 2 ( 0.7%) 0 ( 0.0%)
    1st RPL 4 ( 1.5%) 0 ( 0.0%)
    dLCX 10 ( 3.6%) 0 ( 0.0%)
    dRCA 66 (24.1%) 2 ( 9.5%)
    mLCX 16 ( 5.8%) 1 ( 4.8%)
    mRCA 86 (31.4%) 7 (33.3%)
    OM1 1 ( 0.4%) 0 ( 0.0%)
    OM2 1 ( 0.4%) 0 ( 0.0%)
    pLAD 1 ( 0.4%) 0 ( 0.0%)
    pLCX 4 ( 1.5%) 1 ( 4.8%)
    pRCA 74 (27.0%) 10 (47.6%)
    rPDA 5 ( 1.8%) 0 ( 0.0%)
    rPL 3 ( 1.1%) 0 ( 0.0%)
    RPL 1 ( 0.4%) 0 ( 0.0%)
Vessel 0.927
    L 35 (12.8%) 2 ( 9.5%)
    R 239 (87.2%) 19 (90.5%)
N2BMin 13.2 ± 10.2 16.2 ± 13.4 0.201
D2BMin 113.9 ± 164.6 109.8 ± 95.8 0.859
FT 15.4 ± 11.9 16.6 ± 9.9 0.690
Contrast 183.4 ± 78.3 168.3 ± 59.6 0.434
Entry 0.000
    Femoral 60 (21.9%) 17 (81.0%)
    Radial 214 (78.1%) 4 (19.0%)
EF 56.2 ± 7.9 55.9 ± 10.6 0.898
Ht 165.2 ± 8.5 165.0 ± 8.4 0.905
WT 65.7 ± 12.0 65.4 ± 10.8 0.913
FBS 181.4 ± 76.9 219.9 ± 117.4 0.164
HbA1c 6.8 ± 1.6 6.5 ± 1.5 0.505
Cr 1.0 ± 0.3 1.5 ± 0.9 0.041
CrMax 1.1 ± 0.6 1.6 ± 1.0 0.052
CPKMax 2108.4 ± 2782.6 3337.9 ± 2830.1 0.058
CKMB 151.3 ± 111.9 247.4 ± 170.6 0.022
Trop 1.6 ± 3.1 4.5 ± 5.0 0.019
TC 178.1 ± 45.7 166.4 ± 43.0 0.270
LDL 115.7 ± 41.2 113.9 ± 37.6 0.849
HDL 37.5 ± 11.2 40.6 ± 19.1 0.482
TG 117.2 ± 85.2 77.0 ± 38.5 0.000
hs.CRP 1.8 ± 9.9 1.6 ± 1.7 0.766
DM 0.512
    0 193 (70.4%) 16 (80.0%)
    1 81 (29.6%) 4 (20.0%)
HBP 0.716
    0 128 (46.9%) 8 (40.0%)
    1 145 (53.1%) 12 (60.0%)
IHD 1.000
    0 238 (87.2%) 17 (85.0%)
    1 35 (12.8%) 3 (15.0%)
CVA 0.794
    0 257 (94.1%) 18 (90.0%)
    1 16 ( 5.9%) 2 (10.0%)
PAD 1.000
    0 269 (98.5%) 20 (100.0%)
    1 4 ( 1.5%) 0 ( 0.0%)
CKD 0.878
    0 264 (96.7%) 20 (100.0%)
    1 9 ( 3.3%) 0 ( 0.0%)
smoking 0.063
    0 92 (33.7%) 11 (55.0%)
    1 39 (14.3%) 4 (20.0%)
    2 142 (52.0%) 5 (25.0%)
Statin 1.000
    0 232 (85.0%) 17 (85.0%)
    1 41 (15.0%) 3 (15.0%)
HD 6.3 ± 3.8 8.8 ± 4.7 0.005
DX 0.811
    Aborted MI 5 ( 1.8%) 0 ( 0.0%)
    NSTEMI 6 ( 2.2%) 0 ( 0.0%)
    old MI 1 ( 0.4%) 0 ( 0.0%)
    STEMI 262 (95.6%) 21 (100.0%)
EKGDx 0.142
2 ( 0.7%) 0 ( 0.0%)
    anterior, inferior 5 ( 1.8%) 0 ( 0.0%)
    anterior, inferior, posterior 1 ( 0.4%) 0 ( 0.0%)
    inferior 194 (70.8%) 16 (76.2%)
    inferior, lateral 41 (15.0%) 1 ( 4.8%)
    inferior, lateral, posterior 4 ( 1.5%) 0 ( 0.0%)
    inferior, lateral, RV 1 ( 0.4%) 0 ( 0.0%)
    inferior, posterior 14 ( 5.1%) 0 ( 0.0%)
    inferior, posterior, RV 2 ( 0.7%) 0 ( 0.0%)
    inferior, RV 10 ( 3.6%) 4 (19.0%)
AI 1.000
    0 270 (98.5%) 21 (100.0%)
    1 4 ( 1.5%) 0 ( 0.0%)
LWI 0.453
    0 223 (81.4%) 19 (90.5%)
    1 51 (18.6%) 2 ( 9.5%)
precordialSTdep 0.981
    0 188 (68.6%) 15 (71.4%)
    1 86 (31.4%) 6 (28.6%)
RVI 0.174
    0 261 (95.3%) 18 (85.7%)
    1 13 ( 4.7%) 3 (14.3%)
RHYTHMBPCI 0.180
    AF 16 ( 5.8%) 4 (19.0%)
    Junctional escape 2 ( 0.7%) 1 ( 4.8%)
    junctional rhythm 3 ( 1.1%) 1 ( 4.8%)
    RSR 239 (87.2%) 14 (66.7%)
    Sinus 1 ( 0.4%) 0 ( 0.0%)
    sinus bradycardia 10 ( 3.6%) 1 ( 4.8%)
    ventricular escape 1 ( 0.4%) 0 ( 0.0%)
    Ventricular escape 1 ( 0.4%) 0 ( 0.0%)
    ventricular rhythm 1 ( 0.4%) 0 ( 0.0%)
Rate 0.000
    <40bpm 19 ( 6.9%) 9 (42.9%)
    > 60 1 ( 0.4%) 0 ( 0.0%)
    >60bpm 170 (62.0%) 7 (33.3%)
    40-60 2 ( 0.7%) 0 ( 0.0%)
    40-60bpm 82 (29.9%) 5 (23.8%)
AVB 0.000
    0 AV block 4 ( 1.5%) 0 ( 0.0%)
    0 AV Block 179 (65.3%) 5 (23.8%)
    1st degree 32 (11.7%) 0 ( 0.0%)
    Advanced 13 ( 4.7%) 1 ( 4.8%)
    Complete 46 (16.8%) 15 (71.4%)
BBB 0.030
    0ne 255 (93.1%) 18 (85.7%)
    LBBB 6 ( 2.2%) 0 ( 0.0%)
    RBBB 7 ( 2.6%) 3 (14.3%)
    Ventricular 6 ( 2.2%) 0 ( 0.0%)
VTVF 0.008
    0 258 (94.2%) 16 (76.2%)
    1 16 ( 5.8%) 5 (23.8%)
CPR 0.005
    0 265 (96.7%) 17 (81.0%)
    1 9 ( 3.3%) 4 (19.0%)
DC 0.018
    0 262 (95.6%) 17 (81.0%)
    1 12 ( 4.4%) 4 (19.0%)
IABP 0.000
    0 269 (98.2%) 15 (71.4%)
    1 5 ( 1.8%) 6 (28.6%)
ECMO 0.095
    0 274 (100.0%) 20 (95.2%)
    1 0 ( 0.0%) 1 ( 4.8%)
Shock 0.000
    0 242 (88.3%) 4 (19.0%)
    1 32 (11.7%) 17 (81.0%)
Noreflow 0.007
    0 269 (98.2%) 18 (85.7%)
    1 5 ( 1.8%) 3 (14.3%)
Temporary 0.000
    0t Done 274 (100.0%) 0 ( 0.0%)
    before PCI, removed at Cath Room 0 ( 0.0%) 7 (33.3%)
    performed after PCI 0 ( 0.0%) 1 ( 4.8%)
    performed after PCI, removed at ward 0 ( 0.0%) 5 (23.8%)
    performed before PCI,removed at ward 0 ( 0.0%) 8 (38.1%)
CAVB 0.000
    0 215 (78.5%) 5 (23.8%)
    1 59 (21.5%) 16 (76.2%) 
res1=mytable(Temporary1~.,data=MI[MI$CAVB==1,])
ztable(res1)
0 1 p
(N=59) (N=16)
Age 62.8 ± 10.8 67.4 ± 11.1 0.130
Sex 1.000
    F 17 (28.8%) 5 (31.2%)
    M 42 (71.2%) 11 (68.8%)
death 1.000
    0 56 (94.9%) 15 (93.8%)
    1 3 ( 5.1%) 1 ( 6.2%)
MI 0.481
    0 59 (100.0%) 15 (93.8%)
    1 0 ( 0.0%) 1 ( 6.2%)
TVR 0.000
    0 59 (100.0%) 16 (100.0%)
MACE 0.624
    0 56 (94.9%) 14 (87.5%)
    1 3 ( 5.1%) 2 (12.5%)
Stroke 0.000
    0 59 (100.0%) 16 (100.0%)
Fuday 386.8 ± 270.9 465.2 ± 338.5 0.334
Attempted 0.084
    1 40 (72.7%) 8 (53.3%)
    2 15 (27.3%) 6 (40.0%)
    3 0 ( 0.0%) 1 ( 6.7%)
SuccessfullyDilated 0.117
    1 41 (74.5%) 9 (60.0%)
    2 14 (25.5%) 5 (33.3%)
    3 0 ( 0.0%) 1 ( 6.7%)
CASSNo 0.969
    dRCA 6 (10.2%) 2 (12.5%)
    mLCX 1 ( 1.7%) 0 ( 0.0%)
    mRCA 24 (40.7%) 6 (37.5%)
    OM1 1 ( 1.7%) 0 ( 0.0%)
    pRCA 25 (42.4%) 8 (50.0%)
    rPL 1 ( 1.7%) 0 ( 0.0%)
    RPL 1 ( 1.7%) 0 ( 0.0%)
Vessel 1.000
    L 2 ( 3.4%) 0 ( 0.0%)
    R 57 (96.6%) 16 (100.0%)
N2BMin 11.8 ± 11.9 16.5 ± 15.1 0.194
D2BMin 96.2 ± 105.2 113.6 ± 107.8 0.561
FT 13.5 ± 11.2 18.1 ± 10.9 0.182
Contrast 157.3 ± 68.3 177.4 ± 61.2 0.332
Entry 0.000
    Femoral 14 (23.7%) 12 (75.0%)
    Radial 45 (76.3%) 4 (25.0%)
EF 57.0 ± 6.6 56.9 ± 10.6 0.965
Ht 164.9 ± 8.0 165.6 ± 8.6 0.779
WT 63.2 ± 10.2 64.5 ± 11.7 0.666
FBS 189.2 ± 83.6 212.9 ± 128.6 0.492
HbA1c 6.5 ± 1.3 6.5 ± 1.5 0.940
Cr 1.1 ± 0.5 1.5 ± 1.0 0.156
CrMax 1.3 ± 1.0 1.6 ± 1.1 0.256
CPKMax 2859.4 ± 4826.4 2533.1 ± 1996.7 0.686
CKMB 172.8 ± 113.7 199.5 ± 151.1 0.441
Trop 1.1 ± 2.3 4.7 ± 5.3 0.017
TC 173.3 ± 45.3 171.5 ± 43.4 0.887
LDL 115.3 ± 44.2 118.1 ± 39.3 0.822
HDL 36.7 ± 9.0 42.4 ± 20.5 0.292
TG 120.6 ± 105.8 81.5 ± 39.8 0.025
hs.CRP 1.4 ± 2.3 1.5 ± 1.9 0.924
DM 1.000
    0 46 (78.0%) 12 (75.0%)
    1 13 (22.0%) 4 (25.0%)
HBP 1.000
    0 28 (47.5%) 7 (43.8%)
    1 31 (52.5%) 9 (56.2%)
IHD 1.000
    0 55 (93.2%) 15 (93.8%)
    1 4 ( 6.8%) 1 ( 6.2%)
CVA 0.819
    0 55 (93.2%) 14 (87.5%)
    1 4 ( 6.8%) 2 (12.5%)
PAD 1.000
    0 58 (98.3%) 16 (100.0%)
    1 1 ( 1.7%) 0 ( 0.0%)
CKD 0.658
    0 55 (93.2%) 16 (100.0%)
    1 4 ( 6.8%) 0 ( 0.0%)
smoking 0.036
    0 18 (30.5%) 9 (56.2%)
    1 5 ( 8.5%) 3 (18.8%)
    2 36 (61.0%) 4 (25.0%)
Statin 1.000
    0 55 (93.2%) 15 (93.8%)
    1 4 ( 6.8%) 1 ( 6.2%)
HD 6.5 ± 3.7 8.6 ± 4.8 0.059
DX 0.000
    STEMI 59 (100.0%) 16 (100.0%)
EKGDx 0.101
    inferior 37 (62.7%) 12 (75.0%)
    inferior, lateral 11 (18.6%) 0 ( 0.0%)
    inferior, posterior 4 ( 6.8%) 0 ( 0.0%)
    inferior, posterior, RV 2 ( 3.4%) 0 ( 0.0%)
    inferior, RV 5 ( 8.5%) 4 (25.0%)
AI 0.000
    0 59 (100.0%) 16 (100.0%)
LWI 1.000
    0 51 (86.4%) 14 (87.5%)
    1 8 (13.6%) 2 (12.5%)
precordialSTdep 1.000
    0 39 (66.1%) 11 (68.8%)
    1 20 (33.9%) 5 (31.2%)
RVI 0.761
    0 52 (88.1%) 13 (81.2%)
    1 7 (11.9%) 3 (18.8%)
RHYTHMBPCI 0.981
    AF 8 (13.6%) 2 (12.5%)
    Junctional escape 2 ( 3.4%) 1 ( 6.2%)
    junctional rhythm 3 ( 5.1%) 1 ( 6.2%)
    RSR 42 (71.2%) 11 (68.8%)
    sinus bradycardia 2 ( 3.4%) 1 ( 6.2%)
    ventricular escape 1 ( 1.7%) 0 ( 0.0%)
    Ventricular escape 1 ( 1.7%) 0 ( 0.0%)
Rate 0.177
    <40bpm 16 (27.1%) 9 (56.2%)
    >60bpm 12 (20.3%) 2 (12.5%)
    40-60 1 ( 1.7%) 0 ( 0.0%)
    40-60bpm 30 (50.8%) 5 (31.2%)
AVB 0.282
    Advanced 13 (22.0%) 1 ( 6.2%)
    Complete 46 (78.0%) 15 (93.8%)
BBB 0.538
    0ne 49 (83.1%) 14 (87.5%)
    LBBB 3 ( 5.1%) 0 ( 0.0%)
    RBBB 4 ( 6.8%) 2 (12.5%)
    Ventricular 3 ( 5.1%) 0 ( 0.0%)
VTVF 0.469
    0 54 (91.5%) 13 (81.2%)
    1 5 ( 8.5%) 3 (18.8%)
CPR 1.000
    0 55 (93.2%) 15 (93.8%)
    1 4 ( 6.8%) 1 ( 6.2%)
DC 0.819
    0 55 (93.2%) 14 (87.5%)
    1 4 ( 6.8%) 2 (12.5%)
IABP 0.039
    0 58 (98.3%) 13 (81.2%)
    1 1 ( 1.7%) 3 (18.8%)
ECMO 0.000
    0 59 (100.0%) 16 (100.0%)
Shock 0.000
    0 45 (76.3%) 4 (25.0%)
    1 14 (23.7%) 12 (75.0%)
Noreflow 0.105
    0 57 (96.6%) 13 (81.2%)
    1 2 ( 3.4%) 3 (18.8%)
Temporary 0.000
    0t Done 59 (100.0%) 0 ( 0.0%)
    before PCI, removed at Cath Room 0 ( 0.0%) 7 (43.8%)
    performed after PCI 0 ( 0.0%) 1 ( 6.2%)
    performed after PCI, removed at ward 0 ( 0.0%) 5 (31.2%)
    performed before PCI,removed at ward 0 ( 0.0%) 3 (18.8%)
CAVB 0.000
    1 59 (100.0%) 16 (100.0%) 
ztable(mytable(CAVB~.,data=MI))
0 1 p
(N=220) (N=75)
Age 60.5 ± 13.0 63.8 ± 10.9 0.055
Sex 0.727
    F 58 (26.4%) 22 (29.3%)
    M 162 (73.6%) 53 (70.7%)
death 0.899
    0 211 (95.9%) 71 (94.7%)
    1 9 ( 4.1%) 4 ( 5.3%)
MI 0.981
    0 215 (97.7%) 74 (98.7%)
    1 5 ( 2.3%) 1 ( 1.3%)
TVR 0.424
    0 215 (97.7%) 75 (100.0%)
    1 5 ( 2.3%) 0 ( 0.0%)
MACE 0.962
    0 203 (92.3%) 70 (93.3%)
    1 17 ( 7.7%) 5 ( 6.7%)
Stroke 0.207
    0 212 (96.4%) 75 (100.0%)
    1 8 ( 3.6%) 0 ( 0.0%)
Fuday 395.2 ± 325.8 403.5 ± 286.0 0.843
Attempted 0.215
    1 164 (76.6%) 48 (68.6%)
    2 39 (18.2%) 21 (30.0%)
    3 8 ( 3.7%) 1 ( 1.4%)
    4 2 ( 0.9%) 0 ( 0.0%)
    5 1 ( 0.5%) 0 ( 0.0%)
SuccessfullyDilated 1.3 ± 0.6 1.3 ± 0.5 0.602
CASSNo 0.001
    1st OM 2 ( 0.9%) 0 ( 0.0%)
    1st RPL 4 ( 1.8%) 0 ( 0.0%)
    dLCX 10 ( 4.5%) 0 ( 0.0%)
    dRCA 60 (27.3%) 8 (10.7%)
    mLCX 16 ( 7.3%) 1 ( 1.3%)
    mRCA 63 (28.6%) 30 (40.0%)
    OM1 0 ( 0.0%) 1 ( 1.3%)
    OM2 1 ( 0.5%) 0 ( 0.0%)
    pLAD 1 ( 0.5%) 0 ( 0.0%)
    pLCX 5 ( 2.3%) 0 ( 0.0%)
    pRCA 51 (23.2%) 33 (44.0%)
    rPDA 5 ( 2.3%) 0 ( 0.0%)
    rPL 2 ( 0.9%) 1 ( 1.3%)
    RPL 0 ( 0.0%) 1 ( 1.3%)
Vessel 0.005
    L 35 (15.9%) 2 ( 2.7%)
    R 185 (84.1%) 73 (97.3%)
N2BMin 13.6 ± 9.6 12.8 ± 12.7 0.603
D2BMin 118.3 ± 175.6 99.9 ± 105.2 0.281
FT 15.9 ± 12.0 14.3 ± 11.2 0.331
Contrast 189.8 ± 79.3 160.9 ± 67.1 0.006
Entry 0.071
    Femoral 51 (23.2%) 26 (34.7%)
    Radial 169 (76.8%) 49 (65.3%)
EF 55.9 ± 8.3 57.0 ± 7.5 0.324
Ht 165.2 ± 8.7 165.1 ± 8.1 0.894
WT 66.5 ± 12.3 63.5 ± 10.4 0.063
FBS 180.6 ± 75.1 194.2 ± 94.5 0.258
HbA1c 6.9 ± 1.6 6.5 ± 1.3 0.099
Cr 1.0 ± 0.3 1.2 ± 0.7 0.010
CrMax 1.0 ± 0.4 1.4 ± 1.0 0.009
CPKMax 1986.7 ± 1965.3 2789.8 ± 4368.5 0.127
CKMB 150.8 ± 117.3 178.5 ± 122.0 0.082
Trop 1.7 ± 3.3 1.8 ± 3.5 0.840
TC 178.8 ± 45.9 172.9 ± 44.6 0.334
LDL 115.5 ± 40.3 115.9 ± 43.0 0.941
HDL 37.6 ± 11.7 37.9 ± 12.4 0.858
TG 115.1 ± 78.6 112.2 ± 96.7 0.817
hs.CRP 2.0 ± 11.1 1.4 ± 2.2 0.588
DM 0.217
    0 151 (68.9%) 58 (77.3%)
    1 68 (31.1%) 17 (22.7%)
HBP 1.000
    0 101 (46.3%) 35 (46.7%)
    1 117 (53.7%) 40 (53.3%)
IHD 0.092
    0 185 (84.9%) 70 (93.3%)
    1 33 (15.1%) 5 ( 6.7%)
CVA 0.619
    0 206 (94.5%) 69 (92.0%)
    1 12 ( 5.5%) 6 ( 8.0%)
PAD 1.000
    0 215 (98.6%) 74 (98.7%)
    1 3 ( 1.4%) 1 ( 1.3%)
CKD 0.353
    0 213 (97.7%) 71 (94.7%)
    1 5 ( 2.3%) 4 ( 5.3%)
smoking 0.515
    0 76 (34.9%) 27 (36.0%)
    1 35 (16.1%) 8 (10.7%)
    2 107 (49.1%) 40 (53.3%)
Statin 0.031
    0 179 (82.1%) 70 (93.3%)
    1 39 (17.9%) 5 ( 6.7%)
HD 6.4 ± 3.9 6.9 ± 4.0 0.276
DX 0.234
    Aborted MI 5 ( 2.3%) 0 ( 0.0%)
    NSTEMI 6 ( 2.7%) 0 ( 0.0%)
    old MI 1 ( 0.5%) 0 ( 0.0%)
    STEMI 208 (94.5%) 75 (100.0%)
EKGDx 0.009
2 ( 0.9%) 0 ( 0.0%)
    anterior, inferior 5 ( 2.3%) 0 ( 0.0%)
    anterior, inferior, posterior 1 ( 0.5%) 0 ( 0.0%)
    inferior 161 (73.2%) 49 (65.3%)
    inferior, lateral 31 (14.1%) 11 (14.7%)
    inferior, lateral, posterior 4 ( 1.8%) 0 ( 0.0%)
    inferior, lateral, RV 1 ( 0.5%) 0 ( 0.0%)
    inferior, posterior 10 ( 4.5%) 4 ( 5.3%)
    inferior, posterior, RV 0 ( 0.0%) 2 ( 2.7%)
    inferior, RV 5 ( 2.3%) 9 (12.0%)
AI 0.550
    0 216 (98.2%) 75 (100.0%)
    1 4 ( 1.8%) 0 ( 0.0%)
LWI 0.300
    0 177 (80.5%) 65 (86.7%)
    1 43 (19.5%) 10 (13.3%)
precordialSTdep 0.749
    0 153 (69.5%) 50 (66.7%)
    1 67 (30.5%) 25 (33.3%)
RVI 0.001
    0 214 (97.3%) 65 (86.7%)
    1 6 ( 2.7%) 10 (13.3%)
RHYTHMBPCI 0.000
    AF 10 ( 4.5%) 10 (13.3%)
    Junctional escape 0 ( 0.0%) 3 ( 4.0%)
    junctional rhythm 0 ( 0.0%) 4 ( 5.3%)
    RSR 200 (90.9%) 53 (70.7%)
    Sinus 1 ( 0.5%) 0 ( 0.0%)
    sinus bradycardia 8 ( 3.6%) 3 ( 4.0%)
    ventricular escape 0 ( 0.0%) 1 ( 1.3%)
    Ventricular escape 0 ( 0.0%) 1 ( 1.3%)
    ventricular rhythm 1 ( 0.5%) 0 ( 0.0%)
Rate 0.000
    <40bpm 3 ( 1.4%) 25 (33.3%)
    > 60 1 ( 0.5%) 0 ( 0.0%)
    >60bpm 163 (74.1%) 14 (18.7%)
    40-60 1 ( 0.5%) 1 ( 1.3%)
    40-60bpm 52 (23.6%) 35 (46.7%)
AVB 0.000
    0 AV block 4 ( 1.8%) 0 ( 0.0%)
    0 AV Block 184 (83.6%) 0 ( 0.0%)
    1st degree 32 (14.5%) 0 ( 0.0%)
    Advanced 0 ( 0.0%) 14 (18.7%)
    Complete 0 ( 0.0%) 61 (81.3%)
BBB 0.012
    0ne 210 (95.5%) 63 (84.0%)
    LBBB 3 ( 1.4%) 3 ( 4.0%)
    RBBB 4 ( 1.8%) 6 ( 8.0%)
    Ventricular 3 ( 1.4%) 3 ( 4.0%)
VTVF 0.261
    0 207 (94.1%) 67 (89.3%)
    1 13 ( 5.9%) 8 (10.7%)
CPR 0.436
    0 212 (96.4%) 70 (93.3%)
    1 8 ( 3.6%) 5 ( 6.7%)
DC 0.398
    0 210 (95.5%) 69 (92.0%)
    1 10 ( 4.5%) 6 ( 8.0%)
IABP 0.620
    0 213 (96.8%) 71 (94.7%)
    1 7 ( 3.2%) 4 ( 5.3%)
ECMO 1.000
    0 219 (99.5%) 75 (100.0%)
    1 1 ( 0.5%) 0 ( 0.0%)
Shock 0.000
    0 197 (89.5%) 49 (65.3%)
    1 23 (10.5%) 26 (34.7%)
Noreflow 0.042
    0 217 (98.6%) 70 (93.3%)
    1 3 ( 1.4%) 5 ( 6.7%)
Temporary 0.000
    0t Done 215 (97.7%) 59 (78.7%)
    before PCI, removed at Cath Room 0 ( 0.0%) 7 ( 9.3%)
    performed after PCI 0 ( 0.0%) 1 ( 1.3%)
    performed after PCI, removed at ward 0 ( 0.0%) 5 ( 6.7%)
    performed before PCI,removed at ward 5 ( 2.3%) 3 ( 4.0%)
Temporary1 0.000
    0 215 (97.7%) 59 (78.7%)
    1 5 ( 2.3%) 16 (21.3%)

생존분석

require(survival)
Loading required package: survival
Loading required package: splines
MI$Sur=Surv(MI$Fuday,MI$MACE==1)
fit=survfit(Sur~1,data=MI)
fit
Call: survfit(formula = Sur ~ 1, data = MI)

records   n.max n.start  events  median 0.95LCL 0.95UCL 
    295     295     295      22      NA      NA      NA 
plot(fit)

plot(fit,conf.int=FALSE)

fit1=survfit(Sur~Temporary1,data=MI)
plot(fit1,col=1:2,lty=1:2,xlim=c(0,600))
legend("bottomleft",legend=c("No","Yes"),col=1:2,lty=1:2)

survdiff(Sur~Temporary1,data=MI)
Call:
survdiff(formula = Sur ~ Temporary1, data = MI)

               N Observed Expected (O-E)^2/E (O-E)^2/V
Temporary1=0 274       19    20.35     0.089      1.18
Temporary1=1  21        3     1.65     1.094      1.18

 Chisq= 1.2  on 1 degrees of freedom, p= 0.277 
out=coxph(Sur~Temporary1,data=MI)
summary(out)
Call:
coxph(formula = Sur ~ Temporary1, data = MI)

  n= 295, number of events= 22 

             coef exp(coef) se(coef)     z Pr(>|z|)
Temporary1 0.6642    1.9429   0.6217 1.068    0.285

           exp(coef) exp(-coef) lower .95 upper .95
Temporary1     1.943     0.5147    0.5745     6.571

Concordance= 0.519  (se = 0.031 )
Rsquare= 0.003   (max possible= 0.529 )
Likelihood ratio test= 0.97  on 1 df,   p=0.3243
Wald test            = 1.14  on 1 df,   p=0.2853
Score (logrank) test = 1.18  on 1 df,   p=0.2766
out1=mycph(Sur~.,data=MI)

 mycph : perform coxph of individual expecting variables

 Call: Sur ~ ., data= MI 

MACE  was excluded : infinite
CASSNo1st RPL CASSNodLCX CASSNodRCA CASSNomLCX CASSNomRCA CASSNoOM1 CASSNoOM2 CASSNopLAD CASSNopLCX CASSNopRCA CASSNorPDA CASSNorPL CASSNoRPL  was excluded : infinite
DXNSTEMI DXold MI DXSTEMI  was excluded : infinite
EKGDxanterior, inferior EKGDxanterior, inferior, posterior EKGDxinferior EKGDxinferior, lateral EKGDxinferior, lateral, posterior EKGDxinferior, lateral, RV EKGDxinferior, posterior EKGDxinferior, posterior, RV EKGDxinferior, RV  was excluded : infinite
AI  was excluded : infinite
RVI  was excluded : infinite
RHYTHMBPCIJunctional escape RHYTHMBPCIjunctional rhythm RHYTHMBPCIRSR RHYTHMBPCISinus RHYTHMBPCIsinus bradycardia RHYTHMBPCIventricular escape RHYTHMBPCIVentricular escape RHYTHMBPCIventricular rhythm  was excluded : infinite
Rate> 60 Rate>60bpm Rate40-60 Rate40-60bpm  was excluded : infinite
AVB0 AV Block AVB1st degree AVBAdvanced AVBComplete  was excluded : infinite
Temporarybefore PCI, removed at Cath Room Temporaryperformed after PCI Temporaryperformed after PCI, removed at ward Temporaryperformed before PCI,removed at ward  was excluded : infinite
out1
                       HR   lcl    ucl     p
Age                  1.07  1.03   1.11 0.001
SexM                 0.58  0.25   1.35 0.205
death               35.40 14.97  83.71 0.000
MI                  19.60  7.43  51.72 0.000
TVR                 18.63  6.72  51.63 0.000
Stroke              14.91  4.85  45.80 0.000
Fuday                0.93  0.90   0.96 0.000
Attempted            1.41  0.79   2.51 0.250
SuccessfullyDilated  1.22  0.64   2.32 0.556
VesselR              1.68  0.39   7.27 0.486
N2BMin               0.99  0.94   1.04 0.676
D2BMin               1.00  1.00   1.00 0.883
FT                   0.99  0.96   1.03 0.733
Contrast             0.99  0.99   1.00 0.140
EntryRadial          0.69  0.29   1.66 0.412
EF                   0.94  0.91   0.98 0.007
Ht                   0.95  0.91   0.99 0.020
WT                   0.94  0.91   0.98 0.001
FBS                  1.00  1.00   1.01 0.046
HbA1c                1.20  0.93   1.54 0.160
Cr                   2.14  1.34   3.43 0.001
CrMax                1.73  1.29   2.33 0.000
CPKMax               1.00  1.00   1.00 0.248
CKMB                 1.00  1.00   1.00 0.640
Trop                 0.99  0.87   1.12 0.826
TC                   0.98  0.97   0.99 0.001
LDL                  0.98  0.96   0.99 0.000
HDL                  0.99  0.96   1.03 0.695
TG                   1.00  0.99   1.00 0.246
hs.CRP               1.01  0.93   1.11 0.776
DM                   1.44  0.60   3.46 0.413
HBP                  2.85  1.05   7.74 0.040
IHD                  4.57  1.95  10.72 0.000
CVA                  1.65  0.38   7.13 0.504
PAD                  4.29  0.57  32.15 0.157
CKD                  1.16  0.16   8.68 0.885
smoking              0.71  0.44   1.12 0.139
Statin               3.74  1.49   9.40 0.005
HD                   1.04  0.95   1.13 0.390
LWI                  0.81  0.24   2.75 0.734
precordialSTdep      0.76  0.30   1.96 0.575
BBBLBBB              2.89  0.38  21.87 0.304
BBBRBBB              2.55  0.59  11.04 0.211
BBBVentricular       2.33  0.31  17.58 0.410
VTVF                 2.19  0.73   6.53 0.162
CPR                  5.41  1.82  16.09 0.002
DC                   3.32  1.12   9.86 0.031
IABP                 3.28  0.96  11.20 0.057
ECMO                40.29  4.96 327.55 0.001
Shock                4.88  2.11  11.27 0.000
Noreflow             1.02  0.13   7.67 0.988
Temporary1           1.94  0.57   6.57 0.285
CAVB                 0.89  0.33   2.40 0.811
HRplot(out1,sig.level=0.05)

out2=coxph(Sur~Temporary1+Shock+CPR+Cr+Age+LDL,data=MI)
HRplot(out2,show.CI=TRUE,type=2)

out1=mycph(Sur~.-CAVB-TVR-Stroke-DX-AI-ECMO,data=MI[MI$CAVB==1,])

 mycph : perform coxph of individual expecting variables

 Call: Sur ~ . - CAVB - TVR - Stroke - DX - AI - ECMO, data= [ Sur ~ . - CAVB - TVR - Stroke - DX - AI - ECMO, data= MI Sur ~ . - CAVB - TVR - Stroke - DX - AI - ECMO, data= MI$CAVB == 1 Sur ~ . - CAVB - TVR - Stroke - DX - AI - ECMO, data=  

SexM  was excluded : infinite
MACE  was excluded : infinite
CASSNomLCX CASSNomRCA CASSNoOM1 CASSNopRCA CASSNorPL CASSNoRPL  was excluded : infinite
VesselR  was excluded : infinite
CVA  was excluded : infinite
PAD  was excluded : infinite
Statin  was excluded : infinite
EKGDxinferior, lateral EKGDxinferior, posterior EKGDxinferior, posterior, RV EKGDxinferior, RV  was excluded : infinite
precordialSTdep  was excluded : infinite
RVI  was excluded : infinite
RHYTHMBPCIJunctional escape RHYTHMBPCIjunctional rhythm RHYTHMBPCIRSR RHYTHMBPCIsinus bradycardia RHYTHMBPCIventricular escape RHYTHMBPCIVentricular escape  was excluded : infinite
Rate>60bpm Rate40-60 Rate40-60bpm  was excluded : infinite
AVBComplete  was excluded : infinite
BBBLBBB BBBRBBB BBBVentricular  was excluded : infinite
IABP  was excluded : infinite
Shock  was excluded : infinite
Temporarybefore PCI, removed at Cath Room Temporaryperformed after PCI Temporaryperformed after PCI, removed at ward Temporaryperformed before PCI,removed at ward  was excluded : infinite
out1
                       HR  lcl    ucl     p
Age                  1.01 0.92   1.11 0.842
death               42.01 4.36 404.29 0.001
MI                   2.84 0.15  55.62 0.491
Fuday                0.93 0.87   1.00 0.043
Attempted            1.07 0.20   5.64 0.935
SuccessfullyDilated  1.16 0.22   6.03 0.857
N2BMin               0.99 0.90   1.09 0.860
D2BMin               1.00 1.00   1.01 0.371
FT                   1.00 0.93   1.09 0.907
Contrast             1.00 0.99   1.01 0.961
EntryRadial          0.99 0.16   6.01 0.995
EF                   0.95 0.88   1.02 0.147
Ht                   0.97 0.86   1.09 0.612
WT                   0.92 0.83   1.02 0.132
FBS                  1.00 0.99   1.01 0.758
HbA1c                1.14 0.64   2.04 0.647
Cr                   3.89 1.65   9.19 0.002
CrMax                2.42 1.39   4.23 0.002
CPKMax               1.00 1.00   1.00 0.263
CKMB                 1.00 1.00   1.01 0.156
Trop                 1.09 0.93   1.28 0.297
TC                   0.97 0.94   1.00 0.026
LDL                  0.96 0.92   1.00 0.032
HDL                  0.96 0.86   1.08 0.525
TG                   1.00 0.99   1.01 0.971
hs.CRP               1.09 0.73   1.64 0.669
DM                   0.73 0.07   7.17 0.789
HBP                  1.27 0.21   7.62 0.794
IHD                  1.53 0.12  19.42 0.743
CKD                  5.00 0.52  48.30 0.164
smoking              3.01 0.54  16.69 0.208
HD                   1.14 0.96   1.34 0.127
LWI                  2.14 0.22  20.64 0.511
VTVF                 2.43 0.37  16.09 0.356
CPR                  5.00 0.52  48.30 0.164
DC                   3.38 0.49  23.42 0.218
Noreflow             4.10 0.42  40.07 0.225
Temporary1           1.84 0.29  11.73 0.520
HRplot(out1,sig.level=0.05)

out2=coxph(Sur~Cr+LDL+Temporary1,data=MI[MI$CAVB==1,])
out2
Call:
coxph(formula = Sur ~ Cr + LDL + Temporary1, data = MI[MI$CAVB == 
    1, ])

              coef exp(coef) se(coef)      z     p
Cr          1.0069     2.737   0.5031  2.001 0.045
LDL        -0.0339     0.967   0.0238 -1.420 0.160
Temporary1 -0.1894     0.827   1.0918 -0.173 0.860

Likelihood ratio test=11.5  on 3 df, p=0.00927  n= 75, number of events= 5 
HRplot(out2,show.CI=TRUE,type=2)

require(rms)
Loading required package: rms
Loading required package: Hmisc
Loading required package: grid
Loading required package: lattice

Attaching package: 'lattice'

The following object is masked from 'package:moonBook':

    densityplot

Loading required package: Formula

Attaching package: 'Hmisc'

The following objects are masked from 'package:base':

    format.pval, round.POSIXt, trunc.POSIXt, units

Loading required package: SparseM

Attaching package: 'SparseM'

The following object is masked from 'package:base':

    backsolve
f=npsurv(Sur~Temporary1,data=MI)
survplot(f,xlim=c(0,600),col=1:2,xlab="Follow-up days")

survplot(f,xlim=c(0,600),col=1:2,xlab="Follow-up days",conf="none")

survplot(f,xlim=c(0,600),col=1:2,xlab="Follow-up days",conf="none",label.curve=FALSE,
         label.curves=list(keys="lines"))

survplot(f,xlim=c(0,600),col=1:2,xlab="Follow-up days",conf="none",label.curve=FALSE,
         label.curves=list(keys="lines"),n.risk=TRUE)

survplot(f,xlim=c(0,600),ylim=c(0.6,1),col=1:2,xlab="Follow-up days",conf="none",
         label.curve=FALSE,label.curves=list(keys="lines"),n.risk=TRUE)
text(360,0.92,"log rank p=0.239")

require(GGally)
## Loading required package: GGally
fit1=survfit(Sur~Temporary1,data=MI)
ggsurv(fit1)
## Loading required package: scales