Data Summary

outcome <- read.csv("data/outcome-of-care-measures.csv", colClasses = "character")
head(outcome)
str(outcome[,c(2,7,11,17,19,23)])

'data.frame':   4706 obs. of  6 variables:
 $ Hospital.Name                                                                        : chr  "SOUTHEAST ALABAMA MEDICAL CENTER" "MARSHALL MEDICAL CENTER SOUTH" "ELIZA COFFEE MEMORIAL HOSPITAL" "MIZELL MEMORIAL HOSPITAL" ...
 $ State                                                                                : chr  "AL" "AL" "AL" "AL" ...
 $ Hospital.30.Day.Death..Mortality..Rates.from.Heart.Attack                            : num  14.3 18.5 18.1 NA NA NA 17.7 18 15.9 NA ...
 $ Hospital.30.Day.Death..Mortality..Rates.from.Heart.Failure                           : chr  "11.4" "15.2" "11.3" "13.6" ...
 $ Lower.Mortality.Estimate...Hospital.30.Day.Death..Mortality..Rates.from.Heart.Failure: chr  "9.5" "12.2" "9.1" "10.0" ...
 $ Hospital.30.Day.Death..Mortality..Rates.from.Pneumonia                               : chr  "10.9" "13.9" "13.4" "14.9" ...
summary(outcome[,c(2,7, 11,17,19,23)])

 Hospital.Name         State           Hospital.30.Day.Death..Mortality..Rates.from.Heart.Attack
 Length:4706        Length:4706        Min.   :10.10                                            
 Class :character   Class :character   1st Qu.:14.50                                            
 Mode  :character   Mode  :character   Median :15.40                                            
                                       Mean   :15.45                                            
                                       3rd Qu.:16.40                                            
                                       Max.   :21.90                                            
                                       NA's   :1986                                             
 Hospital.30.Day.Death..Mortality..Rates.from.Heart.Failure Lower.Mortality.Estimate...Hospital.30.Day.Death..Mortality..Rates.from.Heart.Failure
 Length:4706                                                Length:4706                                                                          
 Class :character                                           Class :character                                                                     
 Mode  :character                                           Mode  :character                                                                     
                                                                                                                                                 
                                                                                                                                                 
                                                                                                                                                 
                                                                                                                                                 
 Hospital.30.Day.Death..Mortality..Rates.from.Pneumonia
 Length:4706                                           
 Class :character                                      
 Mode  :character

30-day mortality rates plot for heart attack

outcome[, 11] <- as.numeric(outcome[, 11])
## You may get a warning about NAs being introduced; that is okay
hist(outcome[, 11])

Finding the best hospital in a state

source("best.R")
best("SC", "heart attack")
[1] "MUSC MEDICAL CENTER"
best("NY", "pneumonia")
[1] "MAIMONIDES MEDICAL CENTER"
best("AK", "pneumonia")
[1] "YUKON KUSKOKWIM DELTA REG HOSPITAL"
rankhospital("NC", "heart attack", "worst")
[1] "WAYNE MEMORIAL HOSPITAL"
rankhospital("WA", "heart attack", 7)
[1] "YAKIMA VALLEY MEMORIAL HOSPITAL"
rankhospital("TX", "pneumonia", 10)
[1] "SETON SMITHVILLE REGIONAL HOSPITAL"
rankhospital("NY", "heart attack", 7)
[1] "BELLEVUE HOSPITAL CENTER"
r <- rankall("heart attack", 4)
as.character(subset(r, state == "HI")$hospital)
[1] "CASTLE MEDICAL CENTER"
r <- rankall("pneumonia", "worst")
as.character(subset(r, state == "NJ")$hospital)
[1] "BERGEN REGIONAL MEDICAL CENTER"
r <- rankall("heart failure", 10)
as.character(subset(r, state == "NV")$hospital)
[1] "RENOWN SOUTH MEADOWS MEDICAL CENTER"
LS0tCnRpdGxlOiAiUHJvZ3JhbW1pbmcgQXNzaWdubWVudCAzIgpzdWJ0aXRsZTogIlIgUHJvZ3JhbW1pbmciCmF1dGhvcjogIldhZ25lciBQaW5oZWlybyIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKIyMgRGF0YSBTdW1tYXJ5CgpgYGB7ciwgd2FybmluZz1GQUxTRX0Kb3V0Y29tZSA8LSByZWFkLmNzdigiZGF0YS9vdXRjb21lLW9mLWNhcmUtbWVhc3VyZXMuY3N2IiwgY29sQ2xhc3NlcyA9ICJjaGFyYWN0ZXIiKQpoZWFkKG91dGNvbWUpCmBgYAoKCmBgYHtyLCB3YXJuaW5nPUZBTFNFfQpzdHIob3V0Y29tZVssYygyLDcsMTEsMTcsMTksMjMpXSkKYGBgCgpgYGB7ciwgd2FybmluZz1GQUxTRX0Kc3VtbWFyeShvdXRjb21lWyxjKDIsNywxMSwxNywxOSwyMyldKQpgYGAKCiMjIDMwLWRheSBtb3J0YWxpdHkgcmF0ZXMgcGxvdCBmb3IgaGVhcnQgYXR0YWNrCgpgYGB7ciwgd2FybmluZz1GQUxTRX0Kb3V0Y29tZVssIDExXSA8LSBhcy5udW1lcmljKG91dGNvbWVbLCAxMV0pCiMjIFlvdSBtYXkgZ2V0IGEgd2FybmluZyBhYm91dCBOQXMgYmVpbmcgaW50cm9kdWNlZDsgdGhhdCBpcyBva2F5Cmhpc3Qob3V0Y29tZVssIDExXSkKYGBgCgoKIyMgRmluZGluZyB0aGUgYmVzdCBob3NwaXRhbCBpbiBhIHN0YXRlCgpgYGB7ciwgd2FybmluZz1GQUxTRX0Kc291cmNlKCJiZXN0LlIiKQpgYGAKCmBgYHtyLCB3YXJuaW5nPUZBTFNFfQpiZXN0KCJTQyIsICJoZWFydCBhdHRhY2siKQpgYGAKCmBgYHtyLCB3YXJuaW5nPUZBTFNFfQpiZXN0KCJOWSIsICJwbmV1bW9uaWEiKQpgYGAKCmBgYHtyLCB3YXJuaW5nPUZBTFNFfQpiZXN0KCJBSyIsICJwbmV1bW9uaWEiKQpgYGAKCmBgYHtyLCB3YXJuaW5nPUZBTFNFfQpyYW5raG9zcGl0YWwoIk5DIiwgImhlYXJ0IGF0dGFjayIsICJ3b3JzdCIpCmBgYAoKYGBge3IsIHdhcm5pbmc9RkFMU0V9CnJhbmtob3NwaXRhbCgiV0EiLCAiaGVhcnQgYXR0YWNrIiwgNykKYGBgCgpgYGB7ciwgd2FybmluZz1GQUxTRX0KcmFua2hvc3BpdGFsKCJUWCIsICJwbmV1bW9uaWEiLCAxMCkKYGBgCgpgYGB7ciwgd2FybmluZz1GQUxTRX0KcmFua2hvc3BpdGFsKCJOWSIsICJoZWFydCBhdHRhY2siLCA3KQpgYGAKCmBgYHtyLCB3YXJuaW5nPUZBTFNFfQpyIDwtIHJhbmthbGwoImhlYXJ0IGF0dGFjayIsIDQpCmFzLmNoYXJhY3RlcihzdWJzZXQociwgc3RhdGUgPT0gIkhJIikkaG9zcGl0YWwpCmBgYAoKYGBge3IsIHdhcm5pbmc9RkFMU0V9CnIgPC0gcmFua2FsbCgicG5ldW1vbmlhIiwgIndvcnN0IikKYXMuY2hhcmFjdGVyKHN1YnNldChyLCBzdGF0ZSA9PSAiTkoiKSRob3NwaXRhbCkKYGBgCgpgYGB7ciwgd2FybmluZz1GQUxTRX0KciA8LSByYW5rYWxsKCJoZWFydCBmYWlsdXJlIiwgMTApCmFzLmNoYXJhY3RlcihzdWJzZXQociwgc3RhdGUgPT0gIk5WIikkaG9zcGl0YWwpCmBgYA==