1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
|
hist.pscore <- function(x, numdraws=5000, xlab="Propensity Score", main=NULL, freq=F, xlim = NULL,...){
treat <- x$treat
pscore <- x$distance
weights <- x$weights
matched <- weights!=0
q.cut <- x$q.cut
cwt <- sqrt(weights)
ratio <- x$call$ratio
if(is.null(ratio)){ratio <- 1}
## For full or ratio matching, sample numdraws observations using the weights
if(identical(x$call$method,"full") | (ratio!=1)) {
pscore.treated.matched <- sample(names(treat)[treat==1],
numdraws/2, replace=TRUE,
prob=x$weights[treat==1])
pscore.treated.matched <- pscore[pscore.treated.matched]
pscore.control.matched <- sample(names(treat)[treat==0],
numdraws/2, replace=TRUE,
prob=x$weights[treat==0])
pscore.control.matched <- pscore[pscore.control.matched]
} else {
pscore.treated.matched <- pscore[treat==1 & weights!=0]
pscore.control.matched <- pscore[treat==0 & weights!=0]
}
par(mfrow=c(2,2))
if(!is.null(xlim)){warning("xlim may not be user specified. xlim returned to default.")}
xlim <- range(na.omit(pscore))
if(is.null(main)){
hist(pscore[treat==1],xlim=xlim,
xlab=xlab, freq=freq,
main="Raw Treated", ...)
hist(pscore.treated.matched,xlim=xlim,
xlab=xlab, freq=freq,
main="Matched Treated",...)
if(!is.null(q.cut)){abline(v=q.cut,col="grey",lty=1)}
hist(pscore[treat==0],xlim=xlim,
xlab=xlab, freq=freq,
main="Raw Control",...)
hist(pscore.control.matched,xlim=xlim,
xlab=xlab, freq=freq,
main="Matched Control",...)
if(!is.null(q.cut)){abline(v=q.cut,col="grey",lty=1)}
}else{
hist(pscore[treat==1],xlim=xlim,
xlab=xlab, freq=freq,
main=main, ...)
hist(pscore.treated.matched,xlim=xlim,
xlab=xlab, freq=freq,
main=main,...)
if(!is.null(q.cut)){abline(v=q.cut,col="grey",lty=1)}
hist(pscore[treat==0],xlim=xlim,
xlab=xlab, freq=freq,
main=main,...)
hist(pscore.control.matched,xlim=xlim,
xlab=xlab, freq=freq,
main=main,...)
if(!is.null(q.cut)){abline(v=q.cut,col="grey",lty=1)}
}
}
|
