Title: PROFILING, RANKING
1- PROFILING, RANKING
- League Tables
2(No Transcript)
3RANKING IN THE NEWS
4LETTERMANS TOP 10 LIST
5NEW YORKS MOST DEADLY CARDIAC SURGEONS!!!!
6(No Transcript)
7HOPKINS IS THELEADING SPH!!!
8PROFILING(League Tables)
- The process of comparing units on an outcome
measure with relative or normative standards - Quality of care, use of services, cost
- Educational quality
- Disease rates in small areas
- Gene expression
- Best of breed livestock
- Developing and implementing performance indices
to compare physicians, hospitals, schools,
teachers, genes, ........
9PROFILING OBJECTIVES(in health services)
- Estimate and compare provider-specific
performance measures - Utilization/cost
- Process measures
- Clinical outcomes
- Patient satisfaction/QoL
- Compare using a normative (external) or
- a relative (internal) standard
10(No Transcript)
11RANKING IS EASY
- Just compute estimates order them
12MLE ESTIMATED SMRs
13RANKING IS DIFFICULT
- Need to trade-off the
- estimates and uncertainties
14MLE ESTIMATED SMRs 95 CIs
15Statistical Challenges
- Need a valid method of adjusting for
- case mix and other features
- Patient, physician and hospital characteristics
- But, beware of over adjustment
- Need a valid model for stochastic properties
- Account for variation at all levels
- Account for within-hospital, within-patient
correlations - Need to
- Adjust for systematic variation
- Estimate and account for statistical variation
16PROPER USE OF STATISTICAL SUMMARIES
- The challenge
- Differences in standard errors of
hospital-specific estimates invalidate direct
comparisons - In any case, large SEs make comparisons imprecise
- Consequence
- Even after valid case mix adjustment, differences
in directly estimated performance are due, in
part, to sampling variability - (Partial) Solution, use
- Shrinkage estimates to balance and reduce
variability - Goal-specific estimates to hit the right target
17Comparing performance measures
- Ranks/percentiles, of
- Direct estimates (MLEs)
- Shrunken estimates (BLUPs, Posterior Means)
- Z-scores testing H0 that a unit is just like
others - Optimal (best) ranks or percentiles
- Other measures
- Probability of a large difference between
unit-specific true and H0-generated event rates - Probability of excess mortality
- For the typical patient, on average or for a
specific patient type - Z-score/P-value declarations
- ....
18(No Transcript)
19USRDS
20USRDS
21(No Transcript)
22(No Transcript)
23(No Transcript)
24MLE ESTIMATED SMRs CIs
25Poisson-Normal Model(N, Yk , emortk) are
inputs
- model
-
- precdgamma(0.00001,0.00001)
- for (k in 1N)
- logsmrkdnorm(0,prec)
- smrklt-exp(logsmrk)
- rateklt-emortksmrk
- Yk dpois(ratek)
-
-
- Monitor the SMRk
26MLE, SE POSTERIOR MEAN SMRs (using a
log-normal/Poisson model)
SE
MLE
PM
27Posterior Mean estimated SMRs CIs using a
log-normal/Poisson model (original scale)
28Posterior Mean estimated SMRs CIs using a
Gamma/Poisson model (expanded scale)
29Caterpillar Plot (Hofer et al. JAMA 1999)
- Estimated relative, physician-specific visit
rate and 95 CI - Adjusted for patient demographic and case-mix
- (1.0 is the typical rate)
-
30- Amount that physician-specific, laboratory costs
for diabetic - patients deviates from the mean for all
physicians /(pt. yr.) - Lines show the path from the direct estimate
(the MLE) to the - shrunken estimate (Hofer et al JAMA 1999)
DIRECT
ADJUSTED
31Example using BUGS forhospital performance
ranking
32(No Transcript)
33BUGS Model specification
- model
-
- for k in 1K
- bkdnorm(0, prec)
- rkdbin(pk, nk)
- logit(pk) lt- mu bk
-
- pop.meanlt-exp(mu bb)/(1exp(mu bb))
- mudnorm(0, 1E-6)
- precdgamma(.0001,.0001)
- tausqlt-1/prec
- adddnorm(0, prec)
- bblt- mu add
-
- Monitor the pk and ask for ranks
34Summary Statistics
35Posterior distributions of the ranks
36(No Transcript)
37(No Transcript)
38(No Transcript)
39LOS
X (Posterior Mean-Based Ranks) (Optimal Ranks)
?
40LOS
41BACK TO THE USRDS, SMR EXAMPLE
42Relations among percentiling methods 1998 USRDS
Percentiles
43False detection and non-detection
44Minimize OC
45Advantages of Pk
- Relates percentiles to a substantive scale
- Far easier to compute than
- Shows that the Normand et al. (JASA 1997)
approach is loss function based
46(No Transcript)
47(No Transcript)
48(No Transcript)
49K is large and we can use a completely
non-parametric prior
50? (1-B) ?2/(?2 ?2) ICC
51Probability of being in the upper 10 as a
function of true percentile ? intra-class
correlation
52(No Transcript)