How To Use Propensity Score Analysis: http://www.mc.vanderbilt.edu/crc/workshop_files/2008-04-11.pdf
Vanderbilt Biostatistics Wiki>Main Web>DataSets: http://biostat.mc.vanderbilt.edu/twiki/bin/view/Main/DataSets
Illustration of PS uses: http://rpubs.com/kaz_yos/use-ps
Propensity score demos: http://rpubs.com/kaz_yos/epi271-lab3
AHRQ Summary Variables in Observational Research: Propensity Scores and Disease Risk Scores: http://effectivehealthcare.ahrq.gov/index.cfm/search-for-guides-reviews-and-reports/?productid=1084&pageaction=displayproduct (quotes are from this article).
Colors in ggplot2: http://www.cookbook-r.com/Graphs/Colors_(ggplot2)/
Holmes: Using Propensity Scores in Quasi-Experimental Designs.
## Load Matching for PS matching
library(Matching)
## Right heart cath dataset
rhc <- read.csv("http://biostat.mc.vanderbilt.edu/wiki/pub/Main/DataSets/rhc.csv")
## first 6 rows
head(rhc)
X cat1 cat2 ca sadmdte dschdte dthdte lstctdte death cardiohx chfhx dementhx psychhx
1 1 COPD <NA> Yes 11142 11151 NA 11382 No 0 0 0 0
2 2 MOSF w/Sepsis <NA> No 11799 11844 11844 11844 Yes 1 1 0 0
3 3 MOSF w/Malignancy MOSF w/Sepsis Yes 12083 12143 NA 12400 No 0 0 0 0
4 4 ARF <NA> No 11146 11183 11183 11182 Yes 0 0 0 0
5 5 MOSF w/Sepsis <NA> No 12035 12037 12037 12036 Yes 0 0 0 0
6 6 COPD <NA> No 12389 12396 NA 12590 No 0 1 0 0
chrpulhx renalhx liverhx gibledhx malighx immunhx transhx amihx age sex edu surv2md1 das2d3pc t3d30 dth30
1 1 0 0 0 1 0 0 0 70.25 Male 12.000 0.641 23.50 30 No
2 0 0 0 0 0 1 1 0 78.18 Female 12.000 0.755 14.75 30 No
3 0 0 0 0 1 1 0 0 46.09 Female 14.070 0.317 18.14 30 No
4 0 0 0 0 0 1 0 0 75.33 Female 9.000 0.441 22.93 30 No
5 0 0 0 0 0 0 0 0 67.91 Male 9.945 0.437 21.05 2 Yes
6 1 0 0 0 0 0 0 0 86.08 Female 8.000 0.665 17.50 30 No
aps1 scoma1 meanbp1 wblc1 hrt1 resp1 temp1 pafi1 alb1 hema1 bili1 crea1 sod1 pot1 paco21 ph1 swang1 wtkilo1
1 46 0 41 22.10 124 10 38.7 68.0 3.5 58.0 1.01 1.2 145 4.000 40 7.359 No RHC 64.7
2 50 0 63 28.90 137 38 38.9 218.3 2.6 32.5 0.70 0.6 137 3.300 34 7.329 RHC 45.7
3 82 0 57 0.05 130 40 36.4 275.5 3.5 21.1 1.01 2.6 146 2.900 16 7.359 RHC 0.0
4 48 0 55 23.30 58 26 35.8 156.7 3.5 26.3 0.40 1.7 117 5.800 30 7.460 No RHC 54.6
5 72 41 65 29.70 125 27 34.8 478.0 3.5 24.0 1.01 3.6 126 5.800 17 7.229 RHC 78.4
6 38 0 115 18.00 134 36 39.2 184.2 3.1 30.5 1.01 1.4 138 5.399 68 7.300 No RHC 54.9
dnr1 ninsclas resp card neuro gastr renal meta hema seps trauma ortho adld3p urin1 race income
1 No Medicare Yes Yes No No No No No No No No 0 NA white Under $11k
2 No Private & Medicare No No No No No No No Yes No No NA 1437 white Under $11k
3 No Private No Yes No No No No No No No No NA 599 white $25-$50k
4 No Private & Medicare Yes No No No No No No No No No NA NA white $11-$25k
5 Yes Medicare No Yes No No No No No No No No NA 64 white Under $11k
6 No Medicare Yes No No No No No No No No No 0 242 white Under $11k
ptid
1 5
2 7
3 9
4 10
5 11
6 12
## Show outcome (death) and exposure (swang1)
addmargins(table(rhc[,c("swang1", "death")]))
death
swang1 No Yes Sum
No RHC 1315 2236 3551
RHC 698 1486 2184
Sum 2013 3722 5735
## Create weights for plotting
rhc$one <- 1
## PS formula
psFormula <- (rhc$swang1 == "RHC") ~ age + sex + race + edu + income + ninsclas + cat1 + das2d3pc + dnr1 + ca + surv2md1 + aps1 + scoma1 + wtkilo1 + temp1 + meanbp1 + resp1 + hrt1 + pafi1 + paco21 + ph1 + wblc1 + hema1 + sod1 + pot1 + crea1 + bili1 + alb1 + resp + card + neuro + gastr + renal + meta + hema + seps + trauma + ortho + cardiohx + chfhx + dementhx + psychhx + chrpulhx + renalhx + liverhx + gibledhx + malighx + immunhx + transhx + amihx
## PS model
psModel <- glm(formula = psFormula,
family = binomial(link = "logit"),
data = rhc)
## PS (predicted probability of treatment)
rhc$ps <- predict(psModel, type = "response")
## PS model diagnostic by ROC
rocPsModel <- roc(swang1 ~ ps, data = rhc)
plot(rocPsModel, legacy.axes = TRUE)
Call:
roc.formula(formula = swang1 ~ ps, data = rhc)
Data: ps in 3551 controls (swang1 No RHC) < 2184 cases (swang1 RHC).
Area under the curve: 0.798
## Matching
listMatch <- Match(Tr = (rhc$swang1 == "RHC"), # Need to be in 0,1
X = log(rhc$ps / (1 - rhc$ps)), # logit of PS,i.e., log(PS/(1-PS)) as matching variable
M = 1, # 1:1 match
caliper = 0.05, # caliper = x * SD(each matching variable)
replace = FALSE,
ties = TRUE,
version = "fast")
## Summary
summary(listMatch)
Estimate... 0
SE......... 0
T-stat..... NaN
p.val...... NA
Original number of observations.............. 5735
Original number of treated obs............... 2184
Matched number of observations............... 1540
Matched number of observations (unweighted). 1540
Caliper (SDs)........................................ 0.05
Number of obs dropped by 'exact' or 'caliper' 644
MatchBalance(formul = psFormula,
data = rhc,
match.out = listMatch)
***** (V1) age *****
Before Matching After Matching
mean treatment........ 60.75 60.663
mean control.......... 61.761 60.855
std mean diff......... -6.4686 -1.2154
mean raw eQQ diff..... 1.684 1.5129
med raw eQQ diff..... 1.582 1.256
max raw eQQ diff..... 4.797 3.696
mean eCDF diff........ 0.028456 0.023341
med eCDF diff........ 0.022275 0.021429
max eCDF diff........ 0.070272 0.062987
var ratio (Tr/Co)..... 0.81749 0.82544
T-test p-value........ 0.022435 0.74462
KS Bootstrap p-value.. < 2.22e-16 0.004
KS Naive p-value...... 0.0000031678 0.0044423
KS Statistic.......... 0.070272 0.062987
***** (V2) sexMale *****
Before Matching After Matching
mean treatment........ 0.58516 0.56948
mean control.......... 0.539 0.55844
std mean diff......... 9.3671 2.2287
mean raw eQQ diff..... 0.046245 0.011039
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.023081 0.0055195
med eCDF diff........ 0.023081 0.0055195
max eCDF diff........ 0.046162 0.011039
var ratio (Tr/Co)..... 0.9771 0.99427
T-test p-value........ 0.00061028 0.54405
***** (V3) raceother *****
Before Matching After Matching
mean treatment........ 0.065018 0.063636
mean control.......... 0.059983 0.06039
std mean diff......... 2.0417 1.3296
mean raw eQQ diff..... 0.0050366 0.0032468
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0025176 0.0016234
med eCDF diff........ 0.0025176 0.0016234
max eCDF diff........ 0.0050352 0.0032468
var ratio (Tr/Co)..... 1.0783 1.0501
T-test p-value........ 0.44645 0.70388
***** (V4) racewhite *****
Before Matching After Matching
mean treatment........ 0.78159 0.77922
mean control.......... 0.77527 0.78442
std mean diff......... 1.529 -1.252
mean raw eQQ diff..... 0.0064103 0.0051948
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0031594 0.0025974
med eCDF diff........ 0.0031594 0.0025974
max eCDF diff........ 0.0063188 0.0051948
var ratio (Tr/Co)..... 0.97998 1.0173
T-test p-value........ 0.57544 0.72825
***** (V5) edu *****
Before Matching After Matching
mean treatment........ 11.856 11.798
mean control.......... 11.569 11.788
std mean diff......... 9.1037 0.30618
mean raw eQQ diff..... 0.29396 0.15741
med raw eQQ diff..... 0 0
max raw eQQ diff..... 3 3
mean eCDF diff........ 0.018112 0.0071587
med eCDF diff........ 0.0075068 0.0058442
max eCDF diff........ 0.051076 0.018831
var ratio (Tr/Co)..... 1.0147 0.95258
T-test p-value........ 0.00079428 0.93422
KS Bootstrap p-value.. < 2.22e-16 0.794
KS Naive p-value...... 0.0017251 0.94767
KS Statistic.......... 0.051076 0.018831
***** (V6) income$11-$25k *****
Before Matching After Matching
mean treatment........ 0.20696 0.21364
mean control.......... 0.20079 0.20974
std mean diff......... 1.5229 0.95026
mean raw eQQ diff..... 0.0064103 0.0038961
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0030856 0.0019481
med eCDF diff........ 0.0030856 0.0019481
max eCDF diff........ 0.0061712 0.0038961
var ratio (Tr/Co)..... 1.023 1.0136
T-test p-value........ 0.57384 0.79131
***** (V7) income$25-$50k *****
Before Matching After Matching
mean treatment........ 0.17995 0.16688
mean control.......... 0.14081 0.16364
std mean diff......... 10.187 0.87046
mean raw eQQ diff..... 0.039377 0.0032468
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.01957 0.0016234
med eCDF diff........ 0.01957 0.0016234
max eCDF diff........ 0.03914 0.0032468
var ratio (Tr/Co)..... 1.22 1.0159
T-test p-value........ 0.00010535 0.80971
***** (V8) incomeUnder $11k *****
Before Matching After Matching
mean treatment........ 0.52427 0.53961
mean control.......... 0.58603 0.54545
std mean diff......... -12.365 -1.1721
mean raw eQQ diff..... 0.061813 0.0058442
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.030882 0.0029221
med eCDF diff........ 0.030882 0.0029221
max eCDF diff........ 0.061765 0.0058442
var ratio (Tr/Co)..... 1.0283 1.002
T-test p-value........ 0.0000049842 0.74524
***** (V9) ninsclasMedicare *****
Before Matching After Matching
mean treatment........ 0.23397 0.23896
mean control.......... 0.26669 0.23766
std mean diff......... -7.7248 0.30444
mean raw eQQ diff..... 0.032509 0.0012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.016356 0.00064935
med eCDF diff........ 0.016356 0.00064935
max eCDF diff........ 0.032711 0.0012987
var ratio (Tr/Co)..... 0.91664 1.0038
T-test p-value........ 0.0052471 0.93254
***** (V10) ninsclasMedicare & Medicaid *****
Before Matching After Matching
mean treatment........ 0.056319 0.061688
mean control.......... 0.070684 0.06039
std mean diff......... -6.23 0.53963
mean raw eQQ diff..... 0.014194 0.0012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0071828 0.00064935
med eCDF diff........ 0.0071828 0.00064935
max eCDF diff........ 0.014366 0.0012987
var ratio (Tr/Co)..... 0.80922 1.0201
T-test p-value........ 0.028239 0.88019
***** (V11) ninsclasNo insurance *****
Before Matching After Matching
mean treatment........ 0.062271 0.057143
mean control.......... 0.05238 0.054545
std mean diff......... 4.0924 1.1187
mean raw eQQ diff..... 0.010073 0.0025974
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0049457 0.0012987
med eCDF diff........ 0.0049457 0.0012987
max eCDF diff........ 0.0098915 0.0025974
var ratio (Tr/Co)..... 1.1766 1.0447
T-test p-value........ 0.12124 0.75335
***** (V12) ninsclasPrivate *****
Before Matching After Matching
mean treatment........ 0.33471 0.31818
mean control.......... 0.27232 0.31688
std mean diff......... 13.218 0.27874
mean raw eQQ diff..... 0.062271 0.0012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.031195 0.00064935
med eCDF diff........ 0.031195 0.00064935
max eCDF diff........ 0.062389 0.0012987
var ratio (Tr/Co)..... 1.1239 1.0022
T-test p-value........ 0.00000070868 0.93786
***** (V13) ninsclasPrivate & Medicare *****
Before Matching After Matching
mean treatment........ 0.22436 0.22987
mean control.......... 0.21008 0.23247
std mean diff......... 3.4217 -0.61713
mean raw eQQ diff..... 0.014194 0.0025974
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0071387 0.0012987
med eCDF diff........ 0.0071387 0.0012987
max eCDF diff........ 0.014277 0.0025974
var ratio (Tr/Co)..... 1.0488 0.99217
T-test p-value........ 0.20429 0.86157
***** (V14) cat1CHF *****
Before Matching After Matching
mean treatment........ 0.095696 0.11623
mean control.......... 0.069558 0.10325
std mean diff......... 8.8832 4.0507
mean raw eQQ diff..... 0.026099 0.012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.013069 0.0064935
med eCDF diff........ 0.013069 0.0064935
max eCDF diff........ 0.026138 0.012987
var ratio (Tr/Co)..... 1.3374 1.1095
T-test p-value........ 0.00059649 0.24819
***** (V15) cat1Cirrhosis *****
Before Matching After Matching
mean treatment........ 0.022436 0.030519
mean control.......... 0.049282 0.02987
std mean diff......... -18.123 0.37738
mean raw eQQ diff..... 0.026557 0.00064935
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.013423 0.00032468
med eCDF diff........ 0.013423 0.00032468
max eCDF diff........ 0.026846 0.00064935
var ratio (Tr/Co)..... 0.46819 1.0211
T-test p-value........ 0.000000026922 0.9156
***** (V16) cat1Colon Cancer *****
Before Matching After Matching
mean treatment........ 0.00045788 0.00064935
mean control.......... 0.0016897 0.00064935
std mean diff......... -5.7566 0
mean raw eQQ diff..... 0.0013736 0
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 0
mean eCDF diff........ 0.00061589 0
med eCDF diff........ 0.00061589 0
max eCDF diff........ 0.0012318 0
var ratio (Tr/Co)..... 0.27137 1
T-test p-value........ 0.13667 1
***** (V17) cat1Coma *****
Before Matching After Matching
mean treatment........ 0.043498 0.05
mean control.......... 0.096029 0.051948
std mean diff......... -25.748 -0.89354
mean raw eQQ diff..... 0.052656 0.0019481
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.026266 0.00097403
med eCDF diff........ 0.026266 0.00097403
max eCDF diff........ 0.052531 0.0019481
var ratio (Tr/Co)..... 0.47937 0.96448
T-test p-value........ 1.9984e-15 0.80715
***** (V18) cat1COPD *****
Before Matching After Matching
mean treatment........ 0.026557 0.037013
mean control.......... 0.11236 0.034416
std mean diff......... -53.355 1.3753
mean raw eQQ diff..... 0.085623 0.0025974
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.042903 0.0012987
med eCDF diff........ 0.042903 0.0012987
max eCDF diff........ 0.085806 0.0025974
var ratio (Tr/Co)..... 0.25924 1.0726
T-test p-value........ < 2.22e-16 0.6921
***** (V19) cat1Lung Cancer *****
Before Matching After Matching
mean treatment........ 0.0022894 0.0032468
mean control.......... 0.0095748 0.0025974
std mean diff......... -15.24 1.1411
mean raw eQQ diff..... 0.007326 0.00064935
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0036427 0.00032468
med eCDF diff........ 0.0036427 0.00032468
max eCDF diff........ 0.0072854 0.00064935
var ratio (Tr/Co)..... 0.24091 1.2492
T-test p-value........ 0.00015944 0.73892
***** (V20) cat1MOSF w/Malignancy *****
Before Matching After Matching
mean treatment........ 0.072344 0.082468
mean control.......... 0.067868 0.07987
std mean diff......... 1.7275 0.94394
mean raw eQQ diff..... 0.0045788 0.0025974
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0022381 0.0012987
med eCDF diff........ 0.0022381 0.0012987
max eCDF diff........ 0.0044761 0.0025974
var ratio (Tr/Co)..... 1.061 1.0296
T-test p-value........ 0.5207 0.79288
***** (V21) cat1MOSF w/Sepsis *****
Before Matching After Matching
mean treatment........ 0.32051 0.24675
mean control.......... 0.14841 0.24481
std mean diff......... 36.87 0.45171
mean raw eQQ diff..... 0.17216 0.0019481
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.086052 0.00097403
med eCDF diff........ 0.086052 0.00097403
max eCDF diff........ 0.1721 0.0019481
var ratio (Tr/Co)..... 1.7235 1.0054
T-test p-value........ < 2.22e-16 0.8972
***** (V22) das2d3pc *****
Before Matching After Matching
mean treatment........ 20.701 20.587
mean control.......... 20.371 20.652
std mean diff......... 6.5434 -1.2794
mean raw eQQ diff..... 0.49681 0.43556
med raw eQQ diff..... 0.28516 0.31641
max raw eQQ diff..... 3.8516 3.8516
mean eCDF diff........ 0.024966 0.017295
med eCDF diff........ 0.013749 0.015584
max eCDF diff........ 0.066851 0.041558
var ratio (Tr/Co)..... 0.8429 0.83397
T-test p-value........ 0.020104 0.73467
KS Bootstrap p-value.. < 2.22e-16 0.112
KS Naive p-value...... 0.000011267 0.13988
KS Statistic.......... 0.066851 0.041558
***** (V23) dnr1Yes *****
Before Matching After Matching
mean treatment........ 0.070971 0.082468
mean control.......... 0.14052 0.083117
std mean diff......... -27.081 -0.23599
mean raw eQQ diff..... 0.069597 0.00064935
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.034777 0.00032468
med eCDF diff........ 0.034777 0.00032468
max eCDF diff........ 0.069553 0.00064935
var ratio (Tr/Co)..... 0.54601 0.99289
T-test p-value........ < 2.22e-16 0.94777
***** (V24) caNo *****
Before Matching After Matching
mean treatment........ 0.79075 0.76169
mean control.......... 0.74683 0.75649
std mean diff......... 10.794 1.2189
mean raw eQQ diff..... 0.043956 0.0051948
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.02196 0.0025974
med eCDF diff........ 0.02196 0.0025974
max eCDF diff........ 0.043919 0.0051948
var ratio (Tr/Co)..... 0.87528 0.98539
T-test p-value........ 0.00011208 0.73164
***** (V25) caYes *****
Before Matching After Matching
mean treatment........ 0.15293 0.17403
mean control.......... 0.17967 0.18182
std mean diff......... -7.427 -2.0546
mean raw eQQ diff..... 0.026557 0.0077922
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.013369 0.0038961
med eCDF diff........ 0.013369 0.0038961
max eCDF diff........ 0.026737 0.0077922
var ratio (Tr/Co)..... 0.87908 0.96626
T-test p-value........ 0.0077861 0.56643
***** (V26) surv2md1 *****
Before Matching After Matching
mean treatment........ 0.56846 0.58686
mean control.......... 0.60721 0.58037
std mean diff......... -19.539 3.292
mean raw eQQ diff..... 0.038596 0.008629
med raw eQQ diff..... 0.039978 0.0041146
max raw eQQ diff..... 0.064 0.03302
mean eCDF diff........ 0.047529 0.01003
med eCDF diff........ 0.041643 0.0084416
max eCDF diff........ 0.095698 0.031169
var ratio (Tr/Co)..... 1.0663 0.941
T-test p-value........ 4.1656e-13 0.36799
KS Bootstrap p-value.. < 2.22e-16 0.438
KS Naive p-value...... 3.5003e-11 0.44297
KS Statistic.......... 0.095698 0.031169
***** (V27) aps1 *****
Before Matching After Matching
mean treatment........ 60.739 56.976
mean control.......... 50.934 57.346
std mean diff......... 48.371 -1.8942
mean raw eQQ diff..... 9.8324 1.0013
med raw eQQ diff..... 10 1
max raw eQQ diff..... 14 28
mean eCDF diff........ 0.079748 0.0085099
med eCDF diff........ 0.057146 0.0071429
max eCDF diff........ 0.21274 0.022727
var ratio (Tr/Co)..... 1.1609 1.0256
T-test p-value........ < 2.22e-16 0.57798
KS Bootstrap p-value.. < 2.22e-16 0.736
KS Naive p-value...... < 2.22e-16 0.82128
KS Statistic.......... 0.21274 0.022727
***** (V28) scoma1 *****
Before Matching After Matching
mean treatment........ 18.973 18.746
mean control.......... 22.253 19.071
std mean diff......... -11.604 -1.1508
mean raw eQQ diff..... 3.266 0.83377
med raw eQQ diff..... 0 0
max raw eQQ diff..... 34 28
mean eCDF diff........ 0.03459 0.0085006
med eCDF diff........ 0.027772 0.0071429
max eCDF diff........ 0.068259 0.018182
var ratio (Tr/Co)..... 0.81163 0.94251
T-test p-value........ 0.000043792 0.75777
KS Bootstrap p-value.. < 2.22e-16 0.592
KS Naive p-value...... 0.0000067334 0.96098
KS Statistic.......... 0.068259 0.018182
***** (V29) wtkilo1 *****
Before Matching After Matching
mean treatment........ 72.36 70.662
mean control.......... 65.04 69.398
std mean diff......... 26.398 4.6031
mean raw eQQ diff..... 7.428 2.5227
med raw eQQ diff..... 6.3 2
max raw eQQ diff..... 47.8 72.1
mean eCDF diff........ 0.066108 0.023683
med eCDF diff........ 0.060428 0.016234
max eCDF diff........ 0.13874 0.071429
var ratio (Tr/Co)..... 0.88347 1.0822
T-test p-value........ < 2.22e-16 0.17925
KS Bootstrap p-value.. < 2.22e-16 < 2.22e-16
KS Naive p-value...... < 2.22e-16 0.00077396
KS Statistic.......... 0.13874 0.071429
***** (V30) temp1 *****
Before Matching After Matching
mean treatment........ 37.595 37.622
mean control.......... 37.633 37.64
std mean diff......... -2.0869 -1.0039
mean raw eQQ diff..... 0.07974 0.13597
med raw eQQ diff..... 0 0.10156
max raw eQQ diff..... 1.5977 3
mean eCDF diff........ 0.0069645 0.012262
med eCDF diff........ 0.0037106 0.0087662
max eCDF diff........ 0.028031 0.038312
var ratio (Tr/Co)..... 1.1024 0.83604
T-test p-value........ 0.43463 0.79236
KS Bootstrap p-value.. 0.156 0.13
KS Naive p-value...... 0.23842 0.20838
KS Statistic.......... 0.028031 0.038312
***** (V31) meanbp1 *****
Before Matching After Matching
mean treatment........ 68.198 73.275
mean control.......... 84.869 73.662
std mean diff......... -48.685 -1.0884
mean raw eQQ diff..... 16.673 1.9256
med raw eQQ diff..... 11 1
max raw eQQ diff..... 49 49
mean eCDF diff........ 0.092555 0.010774
med eCDF diff........ 0.06949 0.0071429
max eCDF diff........ 0.21172 0.044805
var ratio (Tr/Co)..... 0.77589 0.97611
T-test p-value........ < 2.22e-16 0.75394
KS Bootstrap p-value.. < 2.22e-16 0.064
KS Naive p-value...... < 2.22e-16 0.090854
KS Statistic.......... 0.21172 0.044805
***** (V32) resp1 *****
Before Matching After Matching
mean treatment........ 26.652 28.112
mean control.......... 28.978 28.286
std mean diff......... -16.413 -1.2436
mean raw eQQ diff..... 2.3471 0.93571
med raw eQQ diff..... 2 0
max raw eQQ diff..... 12 20
mean eCDF diff........ 0.032485 0.012967
med eCDF diff........ 0.021114 0.0097403
max eCDF diff........ 0.091046 0.046104
var ratio (Tr/Co)..... 1.033 1.0439
T-test p-value........ 0.0000000013589 0.7266
KS Bootstrap p-value.. < 2.22e-16 0.046
KS Naive p-value...... 3.6682e-10 0.075753
KS Statistic.......... 0.091046 0.046104
***** (V33) hrt1 *****
Before Matching After Matching
mean treatment........ 118.93 117.77
mean control.......... 112.87 117.21
std mean diff......... 14.601 1.3955
mean raw eQQ diff..... 7.0481 2.0065
med raw eQQ diff..... 5 1
max raw eQQ diff..... 38 33
mean eCDF diff........ 0.034418 0.0098074
med eCDF diff........ 0.032739 0.0077922
max eCDF diff........ 0.07816 0.029221
var ratio (Tr/Co)..... 1.026 0.90805
T-test p-value........ 0.000000071816 0.70741
KS Bootstrap p-value.. < 2.22e-16 0.398
KS Naive p-value...... 0.0000001335 0.5266
KS Statistic.......... 0.07816 0.029221
***** (V34) pafi1 *****
Before Matching After Matching
mean treatment........ 192.43 212.52
mean control.......... 240.63 208.04
std mean diff......... -45.664 4.1538
mean raw eQQ diff..... 48.13 5.9221
med raw eQQ diff..... 46.312 6.3125
max raw eQQ diff..... 241.88 94.25
mean eCDF diff........ 0.11148 0.014006
med eCDF diff........ 0.12874 0.014286
max eCDF diff........ 0.18158 0.033766
var ratio (Tr/Co)..... 0.8184 1.013
T-test p-value........ < 2.22e-16 0.23222
KS Bootstrap p-value.. < 2.22e-16 0.342
KS Naive p-value...... < 2.22e-16 0.34375
KS Statistic.......... 0.18158 0.033766
***** (V35) paco21 *****
Before Matching After Matching
mean treatment........ 36.792 37.462
mean control.......... 39.953 37.477
std mean diff......... -28.804 -0.13271
mean raw eQQ diff..... 3.1849 0.56345
med raw eQQ diff..... 2 0
max raw eQQ diff..... 21 24
mean eCDF diff........ 0.043943 0.008979
med eCDF diff........ 0.041576 0.0045455
max eCDF diff........ 0.10809 0.029221
var ratio (Tr/Co)..... 0.59347 1.1714
T-test p-value........ < 2.22e-16 0.96987
KS Bootstrap p-value.. < 2.22e-16 0.38
KS Naive p-value...... 3.7859e-14 0.5266
KS Statistic.......... 0.10809 0.029221
***** (V36) ph1 *****
Before Matching After Matching
mean treatment........ 7.3802 7.3905
mean control.......... 7.3935 7.3889
std mean diff......... -11.625 1.5437
mean raw eQQ diff..... 0.013163 0.0062075
med raw eQQ diff..... 0.009766 0.000584
max raw eQQ diff..... 0.13086 0.13086
mean eCDF diff........ 0.014537 0.0071353
med eCDF diff........ 0.0043743 0.0038961
max eCDF diff........ 0.087984 0.02987
var ratio (Tr/Co)..... 1.1322 0.99486
T-test p-value........ 0.000012523 0.6643
KS Bootstrap p-value.. < 2.22e-16 0.344
KS Naive p-value...... 0.0000000016156 0.49798
KS Statistic.......... 0.087984 0.02987
***** (V37) wblc1 *****
Before Matching After Matching
mean treatment........ 16.266 15.898
mean control.......... 15.263 15.863
std mean diff......... 7.9869 0.2768
mean raw eQQ diff..... 1.0778 0.49585
med raw eQQ diff..... 0.89844 0.19922
max raw eQQ diff..... 71 92
mean eCDF diff........ 0.019048 0.0062007
med eCDF diff........ 0.012898 0.0051948
max eCDF diff........ 0.056097 0.023377
var ratio (Tr/Co)..... 1.2088 1.1395
T-test p-value........ 0.0023905 0.93513
KS Bootstrap p-value.. < 2.22e-16 0.712
KS Naive p-value...... 0.00040252 0.79406
KS Statistic.......... 0.056097 0.023377
***** (V38) hema1 *****
Before Matching After Matching
mean treatment........ 30.509 30.952
mean control.......... 32.7 30.855
std mean diff......... -29.541 1.2878
mean raw eQQ diff..... 2.3215 0.66275
med raw eQQ diff..... 2 0.5
max raw eQQ diff..... 7 9.5996
mean eCDF diff........ 0.050813 0.015861
med eCDF diff........ 0.024711 0.016234
max eCDF diff........ 0.14789 0.04026
var ratio (Tr/Co)..... 0.7109 0.87021
T-test p-value........ < 2.22e-16 0.72979
KS Bootstrap p-value.. < 2.22e-16 0.12
KS Naive p-value...... < 2.22e-16 0.16472
KS Statistic.......... 0.14789 0.04026
***** (V39) sod1 *****
Before Matching After Matching
mean treatment........ 136.33 136.64
mean control.......... 137.04 136.57
std mean diff......... -9.2669 1.0017
mean raw eQQ diff..... 0.72756 0.42922
med raw eQQ diff..... 1 0
max raw eQQ diff..... 12 7
mean eCDF diff........ 0.0099576 0.0069512
med eCDF diff........ 0.0029821 0.0045455
max eCDF diff........ 0.056846 0.027922
var ratio (Tr/Co)..... 0.97993 0.86576
T-test p-value........ 0.00069316 0.79157
KS Bootstrap p-value.. < 2.22e-16 0.4
KS Naive p-value...... 0.0003202 0.58562
KS Statistic.......... 0.056846 0.027922
***** (V40) pot1 *****
Before Matching After Matching
mean treatment........ 4.0495 4.0466
mean control.......... 4.0773 4.0218
std mean diff......... -2.741 2.4806
mean raw eQQ diff..... 0.048535 0.050469
med raw eQQ diff..... 0 0
max raw eQQ diff..... 2.5996 0.80078
mean eCDF diff........ 0.0058546 0.0069174
med eCDF diff........ 0.0049272 0.0045455
max eCDF diff........ 0.029526 0.023377
var ratio (Tr/Co)..... 0.9501 0.94903
T-test p-value........ 0.31831 0.49442
KS Bootstrap p-value.. 0.122 0.592
KS Naive p-value...... 0.18911 0.79406
KS Statistic.......... 0.029526 0.023377
***** (V41) crea1 *****
Before Matching After Matching
mean treatment........ 2.4734 2.2616
mean control.......... 1.9236 2.3038
std mean diff......... 26.781 -2.1901
mean raw eQQ diff..... 0.58035 0.2237
med raw eQQ diff..... 0.40015 0.1001
max raw eQQ diff..... 8.0977 4.8984
mean eCDF diff........ 0.039684 0.015114
med eCDF diff........ 0.0037401 0.012987
max eCDF diff........ 0.20108 0.053247
var ratio (Tr/Co)..... 1.0276 0.66442
T-test p-value........ < 2.22e-16 0.58678
KS Bootstrap p-value.. < 2.22e-16 0.018
KS Naive p-value...... < 2.22e-16 0.025398
KS Statistic.......... 0.20108 0.053247
***** (V42) bili1 *****
Before Matching After Matching
mean treatment........ 2.7057 2.5418
mean control.......... 1.9973 2.4608
std mean diff......... 13.291 1.5963
mean raw eQQ diff..... 0.73006 0.17997
med raw eQQ diff..... 0.10986 0.0097656
max raw eQQ diff..... 9.4004 6.6016
mean eCDF diff........ 0.024365 0.0058442
med eCDF diff........ 0.019348 0.0038961
max eCDF diff........ 0.11186 0.085714
var ratio (Tr/Co)..... 1.4504 0.92538
T-test p-value........ 0.00000020365 0.66286
KS Bootstrap p-value.. < 2.22e-16 < 2.22e-16
KS Naive p-value...... 3.9968e-15 0.000024395
KS Statistic.......... 0.11186 0.085714
***** (V43) alb1 *****
Before Matching After Matching
mean treatment........ 2.9776 3.043
mean control.......... 3.1635 3.0429
std mean diff......... -20.095 0.0016884
mean raw eQQ diff..... 0.20252 0.058434
med raw eQQ diff..... 0.2998 0
max raw eQQ diff..... 18 23.2
mean eCDF diff........ 0.035867 0.0084798
med eCDF diff........ 0.021534 0.0058442
max eCDF diff........ 0.09896 0.038961
var ratio (Tr/Co)..... 1.8947 1.9305
T-test p-value........ 4.4409e-16 0.99956
KS Bootstrap p-value.. < 2.22e-16 0.062
KS Naive p-value...... 6.2834e-12 0.19293
KS Statistic.......... 0.09896 0.038961
***** (V44) respYes *****
Before Matching After Matching
mean treatment........ 0.28938 0.33247
mean control.......... 0.41707 0.34935
std mean diff......... -28.151 -3.5826
mean raw eQQ diff..... 0.12775 0.016883
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.063844 0.0084416
med eCDF diff........ 0.063844 0.0084416
max eCDF diff........ 0.12769 0.016883
var ratio (Tr/Co)..... 0.84597 0.97637
T-test p-value........ < 2.22e-16 0.31659
***** (V45) cardYes *****
Before Matching After Matching
mean treatment........ 0.42308 0.38377
mean control.......... 0.28358 0.38247
std mean diff......... 28.229 0.26697
mean raw eQQ diff..... 0.13965 0.0012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.069747 0.00064935
med eCDF diff........ 0.069747 0.00064935
max eCDF diff........ 0.13949 0.0012987
var ratio (Tr/Co)..... 1.2016 1.0013
T-test p-value........ < 2.22e-16 0.94001
***** (V46) neuroYes *****
Before Matching After Matching
mean treatment........ 0.054029 0.07013
mean control.......... 0.16193 0.071429
std mean diff......... -47.715 -0.5084
mean raw eQQ diff..... 0.1076 0.0012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.053948 0.00064935
med eCDF diff........ 0.053948 0.00064935
max eCDF diff........ 0.1079 0.0012987
var ratio (Tr/Co)..... 0.37669 0.98319
T-test p-value........ < 2.22e-16 0.88405
***** (V47) gastrYes *****
Before Matching After Matching
mean treatment........ 0.19231 0.18506
mean control.......... 0.147 0.17208
std mean diff......... 11.493 3.3431
mean raw eQQ diff..... 0.04533 0.012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.022653 0.0064935
med eCDF diff........ 0.022653 0.0064935
max eCDF diff........ 0.045307 0.012987
var ratio (Tr/Co)..... 1.2389 1.0586
T-test p-value........ 0.000011569 0.34036
***** (V48) renalYes *****
Before Matching After Matching
mean treatment........ 0.067766 0.061039
mean control.......... 0.041397 0.057792
std mean diff......... 10.489 1.3558
mean raw eQQ diff..... 0.026557 0.0032468
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.013184 0.0016234
med eCDF diff........ 0.013184 0.0016234
max eCDF diff........ 0.026369 0.0032468
var ratio (Tr/Co)..... 1.5922 1.0525
T-test p-value........ 0.000032073 0.70388
***** (V49) metaYes *****
Before Matching After Matching
mean treatment........ 0.042582 0.048701
mean control.......... 0.048437 0.049351
std mean diff......... -2.8989 -0.30158
mean raw eQQ diff..... 0.0059524 0.00064935
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0029273 0.00032468
med eCDF diff........ 0.0029273 0.00032468
max eCDF diff........ 0.0058546 0.00064935
var ratio (Tr/Co)..... 0.88469 0.98752
T-test p-value........ 0.29815 0.93428
***** (V50) hemaYes *****
Before Matching After Matching
mean treatment........ 0.052656 0.061688
mean control.......... 0.067305 0.06039
std mean diff......... -6.5575 0.53963
mean raw eQQ diff..... 0.014652 0.0012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0073247 0.00064935
med eCDF diff........ 0.0073247 0.00064935
max eCDF diff........ 0.014649 0.0012987
var ratio (Tr/Co)..... 0.79477 1.0201
T-test p-value........ 0.021435 0.88152
***** (V51) sepsYes *****
Before Matching After Matching
mean treatment........ 0.23626 0.21104
mean control.......... 0.14503 0.21883
std mean diff......... 21.473 -1.909
mean raw eQQ diff..... 0.091117 0.0077922
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.045617 0.0038961
med eCDF diff........ 0.045617 0.0038961
max eCDF diff........ 0.091234 0.0077922
var ratio (Tr/Co)..... 1.4555 0.97401
T-test p-value........ < 2.22e-16 0.59155
***** (V52) traumaYes *****
Before Matching After Matching
mean treatment........ 0.015568 0.0084416
mean control.......... 0.005069 0.0097403
std mean diff......... 8.4788 -1.4191
mean raw eQQ diff..... 0.010531 0.0012987
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0052494 0.00064935
med eCDF diff........ 0.0052494 0.00064935
max eCDF diff........ 0.010499 0.0012987
var ratio (Tr/Co)..... 3.0393 0.8678
T-test p-value........ 0.00030679 0.7055
***** (V53) orthoYes *****
Before Matching After Matching
mean treatment........ 0.0018315 0.00064935
mean control.......... 0.00084483 0.00064935
std mean diff......... 2.3071 0
mean raw eQQ diff..... 0.00091575 0
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 0
mean eCDF diff........ 0.00049333 0
med eCDF diff........ 0.00049333 0
max eCDF diff........ 0.00098667 0
var ratio (Tr/Co)..... 2.1661 1
T-test p-value........ 0.3414 1
***** (V54) cardiohx *****
Before Matching After Matching
mean treatment........ 0.20421 0.2026
mean control.......... 0.15967 0.20325
std mean diff......... 11.046 -0.1615
mean raw eQQ diff..... 0.044414 0.00064935
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.02227 0.00032468
med eCDF diff........ 0.02227 0.00032468
max eCDF diff........ 0.044539 0.00064935
var ratio (Tr/Co)..... 1.2114 0.99762
T-test p-value........ 0.000026753 0.96437
***** (V55) chfhx *****
Before Matching After Matching
mean treatment........ 0.1946 0.20325
mean control.......... 0.16784 0.1974
std mean diff......... 6.7572 1.4518
mean raw eQQ diff..... 0.027015 0.0058442
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.013379 0.0029221
med eCDF diff........ 0.013379 0.0029221
max eCDF diff........ 0.026757 0.0058442
var ratio (Tr/Co)..... 1.1223 1.0221
T-test p-value........ 0.01118 0.68282
***** (V56) dementhx *****
Before Matching After Matching
mean treatment........ 0.069139 0.072727
mean control.......... 0.11631 0.077273
std mean diff......... -18.588 -1.7498
mean raw eQQ diff..... 0.047161 0.0045455
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.023583 0.0022727
med eCDF diff........ 0.023583 0.0022727
max eCDF diff........ 0.047166 0.0045455
var ratio (Tr/Co)..... 0.6263 0.94581
T-test p-value........ 7.3188e-10 0.62495
***** (V57) psychhx *****
Before Matching After Matching
mean treatment........ 0.045788 0.054545
mean control.......... 0.080541 0.055195
std mean diff......... -16.623 -0.28585
mean raw eQQ diff..... 0.034799 0.00064935
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.017377 0.00032468
med eCDF diff........ 0.017377 0.00032468
max eCDF diff........ 0.034753 0.00064935
var ratio (Tr/Co)..... 0.59009 0.98891
T-test p-value........ 0.000000056902 0.93599
***** (V58) chrpulhx *****
Before Matching After Matching
mean treatment........ 0.14423 0.1539
mean control.......... 0.21797 0.1539
std mean diff......... -20.983 0
mean raw eQQ diff..... 0.073718 0
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 0
mean eCDF diff........ 0.036868 0
med eCDF diff........ 0.036868 0
max eCDF diff........ 0.073736 0
var ratio (Tr/Co)..... 0.72423 1
T-test p-value........ 6.3505e-13 1
***** (V59) renalhx *****
Before Matching After Matching
mean treatment........ 0.048535 0.05
mean control.......... 0.04196 0.057143
std mean diff......... 3.0589 -3.2763
mean raw eQQ diff..... 0.0064103 0.0071429
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0032874 0.0035714
med eCDF diff........ 0.0032874 0.0035714
max eCDF diff........ 0.0065748 0.0071429
var ratio (Tr/Co)..... 1.149 0.88163
T-test p-value........ 0.24869 0.37696
***** (V60) liverhx *****
Before Matching After Matching
mean treatment........ 0.062271 0.066883
mean control.......... 0.074627 0.067532
std mean diff......... -5.112 -0.25984
mean raw eQQ diff..... 0.012363 0.00064935
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0061779 0.00032468
med eCDF diff........ 0.0061779 0.00032468
max eCDF diff........ 0.012356 0.00064935
var ratio (Tr/Co)..... 0.84572 0.99107
T-test p-value........ 0.069159 0.94043
***** (V61) gibledhx *****
Before Matching After Matching
mean treatment........ 0.024725 0.029221
mean control.......... 0.036891 0.024026
std mean diff......... -7.8326 3.0833
mean raw eQQ diff..... 0.011905 0.0051948
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0060829 0.0025974
med eCDF diff........ 0.0060829 0.0025974
max eCDF diff........ 0.012166 0.0051948
var ratio (Tr/Co)..... 0.67881 1.2097
T-test p-value........ 0.0080415 0.36504
***** (V62) malighx *****
Before Matching After Matching
mean treatment........ 0.2033 0.23117
mean control.......... 0.24556 0.23636
std mean diff......... -10.5 -1.2318
mean raw eQQ diff..... 0.042125 0.0051948
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.021134 0.0025974
med eCDF diff........ 0.021134 0.0025974
max eCDF diff........ 0.042268 0.0051948
var ratio (Tr/Co)..... 0.87441 0.98468
T-test p-value........ 0.00017202 0.7302
***** (V63) immunhx *****
Before Matching After Matching
mean treatment........ 0.29121 0.27792
mean control.......... 0.25542 0.28961
std mean diff......... 7.8754 -2.6083
mean raw eQQ diff..... 0.035714 0.011688
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.017894 0.0058442
med eCDF diff........ 0.017894 0.0058442
max eCDF diff........ 0.035788 0.011688
var ratio (Tr/Co)..... 1.0855 0.97543
T-test p-value........ 0.0032935 0.46761
***** (V64) transhx *****
Before Matching After Matching
mean treatment........ 0.14973 0.12013
mean control.......... 0.09434 0.11948
std mean diff......... 15.519 0.19967
mean raw eQQ diff..... 0.055403 0.00064935
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.027693 0.00032468
med eCDF diff........ 0.027693 0.00032468
max eCDF diff........ 0.055386 0.00064935
var ratio (Tr/Co)..... 1.4903 1.0047
T-test p-value........ 0.0000000011484 0.95508
***** (V65) amihx *****
Before Matching After Matching
mean treatment........ 0.043498 0.031169
mean control.......... 0.029569 0.033117
std mean diff......... 6.8272 -1.1207
mean raw eQQ diff..... 0.014194 0.0019481
med raw eQQ diff..... 0 0
max raw eQQ diff..... 1 1
mean eCDF diff........ 0.0069645 0.00097403
med eCDF diff........ 0.0069645 0.00097403
max eCDF diff........ 0.013929 0.0019481
var ratio (Tr/Co)..... 1.4502 0.94307
T-test p-value........ 0.0075345 0.7607
Before Matching Minimum p.value: < 2.22e-16
Variable Name(s): age edu cat1COPD cat1MOSF w/Sepsis das2d3pc dnr1Yes surv2md1 aps1 scoma1 wtkilo1 meanbp1 resp1 hrt1 pafi1 paco21 ph1 wblc1 hema1 sod1 crea1 bili1 alb1 respYes cardYes neuroYes sepsYes Number(s): 1 5 18 21 22 23 26 27 28 29 31 32 33 34 35 36 37 38 39 41 42 43 44 45 46 51
After Matching Minimum p.value: < 2.22e-16
Variable Name(s): wtkilo1 bili1 Number(s): 29 42
## Extract matched data
psMatchData <- rhc[unlist(listMatch[c("index.treated","index.control")]), ]