Code status proposal

concepts/Code_status/Readme.md

@@ -0,0 +1,9 @@
+The files here create a patient table that tracks the changes to the patient code status. The careplangeneral table was used instead of careplaneol because careplaneol had much less recorded patients, and the type of code status change was recorded in careplangeneral. There are two versions of the files: .Rmd and .sql. The Rmd version is the supported version, while the SQL version is not supported and was added just as a template.
+
+Rmd:
+  - carePlan_getItemValues creates a csv that contains all the cplitemvalues for inspection, in order to choose the values that correspond to patient code status.
+  - carePlan_getPatientCode creates a table that lists the recorded patient code status at a certain time (offset).
+
+SQL:
+  - carePlan_getItemValues creates a materialized view that lists all the cplitemvalues (possible statuses) in descending order of the counts
+  - carePlan_getPatientCode creates a materialized view that lists the recorded patient code status at a certain time (offset).

concepts/Code_status/carePlan_getItemValues.Rmd

@@ -0,0 +1,34 @@
+---
+title: "carePlan_getItemValues"
+output: html_document
+---
+
+```{r setup, include=FALSE}
+knitr::opts_chunk$set(echo = TRUE)
+```
+
+Connect to database.
+
+```{r}
+require("RPostgreSQL")
+drv <- dbDriver("PostgreSQL")
+con <- dbConnect(drv, dbname = "eicu",
+                 host = "localhost", port = 5432,
+                 user = "josephpark", password = rstudioapi::askForPassword("Database password"),
+                 options = "--search_path=eicu")
+dbExistsTable(con, "patient")
+```
+
+Write cplitemvalues to file "possible_status" for inspection on which items indicate patient code status.
+All cplitemvalues that were chosen were in cplgroup = 'Care Limitation' (all cplitemvalues in this cplgroup was used),
+  except cplitemvalue = 'End of life', which was in cplgroup = 'Ordered Protocols'.
+
+```{r, connection = con}
+library(dplyr)
+carePlanGeneral_tbl <- tbl(con, "careplangeneral")
+possible_status <- carePlanGeneral_tbl %>% 
+  group_by(cplitemvalue) %>%
+  summarize(n=n()) %>%
+  arrange(desc(n))
+write.csv(possible_status, "possible_status.csv")
+```

concepts/Code_status/carePlan_getItemValues.sql

@@ -0,0 +1,9 @@
+-- NOT SUPPORTED
+
+DROP MATERIALIZED VIEW IF EXISTS possibleCodeStatus;
+CREATE MATERIALIZED VIEW possibleCodeStatus AS
+
+SELECT "cplitemvalue", count(*) AS "n"
+FROM "careplangeneral"
+GROUP BY "cplitemvalue"
+ORDER BY "n" DESC;

concepts/Code_status/carePlan_getPatientCode.Rmd

@@ -0,0 +1,64 @@
+---
+title: "carePlan_getPatientCode"
+output: html_document
+---
+
+```{r setup, include=FALSE}
+knitr::opts_chunk$set(echo = TRUE)
+```
+
+Connect to database.
+
+```{r}
+require("RPostgreSQL")
+drv <- dbDriver("PostgreSQL")
+con <- dbConnect(drv, dbname = "eicu",
+                 host = "localhost", port = 5432,
+                 user = "josephpark", password = rstudioapi::askForPassword("Database password"),
+                 options = "--search_path=eicu")
+dbExistsTable(con, "patient")
+```
+
+After inspecting table generated from careplan_getitemvalues.Rmd, the values in cpitemvalue_code were selected to relate to patient code status.
+
+```{r, connection = con}
+library(dplyr)
+carePlanGeneral_tbl <- tbl(con, "careplangeneral")
+cplitemvalue_code = c('Full therapy', 'Do not resuscitate', 'No CPR', 'No intubation', 'Comfort measures only',
+                      'No cardioversion', 'No vasopressors/inotropes', 'No augmentation of care', 
+                      'End of life', 'No blood products', 'No blood draws', 'Advance directives')
+```
+
+Generate table that lists the patient along with code status.
+
+```{r, connection = con}
+patientCodeStatus <- carePlanGeneral_tbl %>%
+  filter(cplitemvalue %in% cplitemvalue_code) %>%
+  select(patientunitstayid, cplitemoffset, cplgroup, cplitemvalue) %>%
+  arrange(patientunitstayid, cplitemoffset, cplitemvalue) %>%
+  group_by(patientunitstayid, cplitemoffset) %>%
+  summarize(cplitemvalue = paste(cplitemvalue, collapse=", "))
+patientCodeStatus
+```
+
+Output the file as a csv.
+
+```{r, connection = con}
+write.csv(patientCodeStatus, "patient_status_table.csv", row.names=FALSE)
+```
+
+Optional: sorts patients by which patient had the most number of status updates. 
+
+```{r, connection = con,eval=FALSE}
+patientCodeStatus %>%
+  group_by(patientunitstayid) %>%
+  summarize(n = n()) %>%
+  arrange(desc(n))
+```
+
+Example of a patient. This patient was strange with both full therapy and no CPR at the same time.
+
+```{r, connection = con,eval=FALSE}
+patientCodeStatus %>%
+  filter(patientunitstayid == 266999)
+```

concepts/Code_status/carePlan_getPatientCode.sql

@@ -0,0 +1,13 @@
+-- NOT SUPPORTED
+
+DROP MATERIALIZED VIEW IF EXISTS patientCodeStatus;
+CREATE MATERIALIZED VIEW patientCodeStatus AS
+
+SELECT "patientunitstayid", "cplitemoffset", string_agg("cplitemvalue", ', ') AS "cplitemvalue"
+FROM (SELECT *
+FROM (SELECT "patientunitstayid" AS "patientunitstayid", "cplitemoffset" AS "cplitemoffset", "cplgroup" AS "cplgroup", "cplitemvalue" AS "cplitemvalue"
+FROM (SELECT *
+FROM "careplangeneral"
+WHERE ("cplitemvalue" IN ('Full therapy', 'Do not resuscitate', 'No CPR', 'No intubation', 'Comfort measures only', 'No cardioversion', 'No vasopressors/inotropes', 'No augmentation of care', 'End of life', 'No blood products', 'No blood draws', 'Advance directives'))) "ivtmaaxwzk") "oytldaespk"
+ORDER BY "patientunitstayid", "cplitemoffset", "cplitemvalue") "vrxthajevr"
+GROUP BY "patientunitstayid", "cplitemoffset";