-1

I am reading data from a csv file which has 3 columns (hospital name- character, state-character, mortality rate-numeric) by:

datafile <- read.csv("outcome-of-care-measures.csv", 
    na.strings = "Not Available",
    colClasses = c("character","character","numeric"))

Now I split the data based on state:

## split data based on state name
data_split <- split(datafile,datafile$State)

My problem is to find the “worst” hospital (highest mortality rate) in each state and display the result. For this, first I sorted the data”: (rate is a list) 

for (i in 1:length(data_split)){
  ## remove all rows with NA
  rate[[i]] <- data_split[[i]][complete.cases(data_split[[i]][ ,3]), ]  
  ##sort by mortality and remove
  ## conflict by hospital name
  rate[[i]] <- rate[[i]][order(rate[[i]][, 3],rate[[i]][ ,1]), ]  

}

Program is working but I am getting the wrong hospital name for many states. I am unable to find error in the program.

MrFlick
  • 195,160
  • 17
  • 277
  • 295
bipvan
  • 77
  • 7
  • 3
    Hi. It would be good to have some sort of reproducible example of your data. Check [this topic](http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example). Beside that, a better title for you question may improve your chances to get a better answer, allowing it to be easier to find while searching for topics. – Paulo E. Cardoso May 10 '14 at 13:17

1 Answers1

1

Why split the data.frame?

something like this would help?

df <- data.frame('hospital' = LETTERS[1:6],
                 'state' = rep(c('state1', 'state2', 'state3'),2),
                 'mr' = c(1:6))
df
  hospital  state mr
1        A state1  1
2        B state2  2
3        C state3  3
4        D state1  4
5        E state2  5
6        F state3  6

df2 <- df[with(df, order(-mr, state)), ]

df2[!duplicated(df2$state), ]
  hospital  state mr
6        F state3  6
5        E state2  5
4        D state1  4

You could do it with your approach but maintaining the list with all initial complete entries. But why?

ds <- split(df, df$state) 
rate <- list()
for (i in 1:length(ds)){
  ## remove all rows with NA
  rate[[i]] <- ds[[i]][complete.cases(ds[[i]][ ,3]), ]  
  ##sort by mortality and remove
  ## conflict by hospital name
  rate[[i]] <- rate[[i]][order(- rate[[i]][, 3], rate[[i]][ ,1]), ]  

}
rate

 [[1]]
  hospital  state mr
4        D state1  4
1        A state1  1

[[2]]
  hospital  state mr
5        E state2  5
2        B state2  2

[[3]]
  hospital  state mr
6        F state3  6
3        C state3  3
Paulo E. Cardoso
  • 5,778
  • 32
  • 42
  • Thanks for the suggested changes but I would appreciate if anyone can elaborate on what is wrong in my code. Also, the name of one of the column is, "Hospital 30-Day Death (Mortality) Rates from Heart Failure." How can I access that column with a $ sign? – bipvan May 11 '14 at 06:40
  • @bipvan. You can call it with df$"var name" if necessary or df[ ,1]. I'd say to avoid such names. Use short, lower case, names instead and change them when necessary for outputs from within or outside R. You'll see that while importing objects R will try to eliminate or replace any space on variable names. – Paulo E. Cardoso May 12 '14 at 08:37