Exploring length of stay in the ICU

In [48]:

    

# install.packages("RPostgreSQL")
require("RPostgreSQL")
require("ggplot2")


    

In [4]:

    

# load the PostgreSQL driver
drv <- dbDriver("PostgreSQL")


    

In [19]:

    

# create a connection to the postgres database
# set the search path to the mimiciii schema
con <- dbConnect(drv, dbname = "mimic",
                 host = "localhost", port = 5432,
                 user = "mimic")
dbSendQuery(con, 'set search_path to mimiciii')


    

In [34]:

    

# show a list of tables
dbListTables(con)


    

    
Out[34]:





	
'admissions'
	
'callout'
	
'caregivers'
	
'chartevents'
	
'chartevents_1'
	
'chartevents_2'
	
'chartevents_3'
	
'chartevents_4'
	
'chartevents_5'
	
'chartevents_6'
	
'chartevents_7'
	
'chartevents_8'
	
'chartevents_9'
	
'chartevents_10'
	
'chartevents_11'
	
'chartevents_12'
	
'chartevents_13'
	
'chartevents_14'
	
'cptevents'
	
'datetimeevents'
	
'diagnoses_icd'
	
'drgcodes'
	
'd_cpt'
	
'd_icd_diagnoses'
	
'd_icd_procedures'
	
'd_items'
	
'd_labitems'
	
'icustays'
	
'inputevents_cv'
	
'inputevents_mv'
	
'labevents'
	
'microbiologyevents'
	
'noteevents'
	
'outputevents'
	
'prescriptions'
	
'procedureevents_mv'
	
'procedures_icd'
	
'services'
	
'transfers'
	
'noteevents_deid'
	
'patients'

In [90]:

    

# assign ICU length of stay to variable
iculos = dbGetQuery(con, "select los from icustays where los >= 0 and los <=20")


    

In [93]:

    

# plot histogram of length of stay in the ICU
qplot(iculos$los, geom="histogram", binwidth = 1, main = "Length of stay in the ICU", 
      xlab = "Length of stay, days", fill=I("#9ebcda"), col=I("#FFFFFF"))


    

In [ ]:

    

# close the connection
dbDisconnect(con)
dbUnloadDriver(drv)