-
Notifications
You must be signed in to change notification settings - Fork 61
/
Copy pathREADME.Rmd
167 lines (133 loc) · 9.95 KB
/
README.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
---
title: icd README
output:
md_document:
variant: gfm
date: "`r format(Sys.time(), '%B %d, %Y')`"
author: Jack O. Wasey
---
<!--
Copyright (C) 2014 - 2020 Jack O. Wasey
This file is part of icd.
icd is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
icd is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with icd. If not, see <http:#www.gnu.org/licenses/>.
-->
<!-- README.md is generated from README.Rmd. Please edit that file and render with rmarkdown::render("README.Rmd")
-->
```{r, echo = FALSE}
suppressPackageStartupMessages({
library(icd, warn.conflicts = FALSE)
library(magrittr, warn.conflicts = FALSE)
})
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = "man/figures/README-"
)
patient_data <- data.frame(
visit_id = c(1000, 1000, 1000, 1000, 1001, 1001, 1002, 1000),
icd9 = as.icd9cm(c("40201", "2258", "7208", "25001",
"34400", "4011", "4011", NA)),
stringsAsFactors = FALSE
)
```
# icd
<!-- badges: start -->
[](https://cran.r-project.org/package=icd)
[](https://www.tidyverse.org/lifecycle/#stable)
[](https://www.repostatus.org/#active)
[](https://travis-ci.org/jackwasey/icd)
[](https://ci.appveyor.com/project/jackwasey/icd/branch/main)
[](https://codecov.io/github/jackwasey/icd?branch=main)
[](https://bestpractices.coreinfrastructure.org/projects/2092)
[](https://cran.r-project.org/package=icd)
[](https://cran.r-project.org/package=icd)
<!-- badges: end -->
<!-- [](https://cran.r-project.org/package=icd) -->
# Fast comorbidities from ICD-9 and ICD-10 codes, decoding, manipulation and validation
## Introduction
```{r pkgdesc, results = "asis", echo = FALSE}
cat(packageDescription("icd")$Description)
```
`icd` is used by many researchers around the world who work in public health, epidemiology, clinical research, nutrition, journalism, health administration, insurance, and more. I'm grateful for contact from people in these fields for their feedback and code contributions, and I'm pleased to say that `icd` has been used in works like the [Pulitzer finalist](http://www.pulitzer.org/finalists/staff-propublica) work on [maternal death](http://www.pulitzer.org/finalists/staff-propublica) by [ProPublica](https://www.propublica.org).
## Features
* find comorbidities of patients based on ICD-9 or ICD-10 codes, e.g. Cancer, Heart Disease
+ several standard mappings of ICD codes to comorbidities are included (Quan, Deyo, Elixhauser, AHRQ, PCCC)
+ *very fast* assignment of ICD codes to comorbidities (using novel matrix multiplication algorithm and C++ internally -- see 'efficiency' vignette for details)
* use your existing wide or long data format, icd can guess which columns are ICD-9 or ICD-10 codes.
* explain and summarize groups of ICD codes in natural language, using ICD editions from the WHO, USA, France and Belgium. Many different annual editions of these data are available, and these may be downloaded automatically when used, or in bulk with `download_all_icd_data()`.
* Charlson and Van Walraven score calculations
* Hierarchical Condition Codes (HCC) from CMS
* Clinical Classifications Software (CCS) comorbidities from AHRQ
* Pediatric Complex Chronic Condition comorbidities
* AHRQ ICD-10 procedure code classification
* correct conversion between different representations of ICD codes, with and without a decimal points, leading and trailing characters (this is not trivial for ICD-9-CM). ICD-9 to ICD-10 cross-walk is not yet implemented
* comprehensive test suite to increase confidence in accurate processing of ICD codes
## Examples
See also the vignettes and examples embedded in the help for each function for more. Here's a taste:
```{r example}
# install.packages("icd")
library(icd)
# Typical diagnostic code data, with many-to-many relationship
patient_data
# get comorbidities using Quan's application of Deyo's Charlson comorbidity groups
comorbid_charlson(patient_data)
# or go straight to the Charlson scores:
charlson(patient_data)
# plot summary of Uranium Cancer Registry sample data using AHRQ comorbidities
plot_comorbid(uranium_pathology)
```
## Comorbidities example: make "Table 1" summary data
A common requirement for medical research involving patients is determining new or existing comorbidities. This is often reported in *Table 1* of research papers to demonstrate the similarity or differences of groups of patients. This package is focussed on fast and accurate generation of this comorbidity information from raw lists of ICD-9 and ICD-10 codes.
Here we are using the US National Hospital Discharge Survey 2010 data from the [nhds](https://github.com/jackwasey/nhds) package. For the sake of example, let us compare emergency to other admissions. A real table would have more patient features; this primarily demonstrates how to get ICD codes into your Table 1.
NHDS 2010 comorbidities to demonstrate Table One creation. Presented as counts (percentage prevalence in group).
```{r tableone}
nhds <- nhds::nhds2010
# get the comorbidities using the Quan-Deyo version of the Charlson categories
cmb <- icd::comorbid_quan_deyo(nhds, abbrev_names = FALSE)
nhds <- cbind(nhds, cmb, stringsAsFactors = FALSE)
Y <- nhds$adm_type == "emergency"
tab_dat <- vapply(
unname(unlist(icd_names_charlson)),
function(x) {
c(sprintf("%i (%.2f%%)",
sum(nhds[Y, x]),
100 * mean(nhds[Y, x])),
sprintf("%i (%.2f%%)",
sum(nhds[!Y, x]),
100 * mean(nhds[!Y, x])))
},
character(2)
)
knitr::kable(t(tab_dat), col.names = c("Emergency", "Not emergency"))
```
## How to get help
Look at the help files for details and examples of almost every function in this package. There are several vignettes showing the main features (See list with `vignette(package = "icd")`):
- Introduction `vignette("introduction", package = "icd")`
- Charlson scores `vignette("charlson-scores", package = "icd")`
- Examples using ICD-10 codes `vignette("ICD-10", package = "icd")`
- CMS Hierarchical Condition Codes (HCC) `vignette("CMS-HCC", package = "icd")`
- Pediatric Complex Chronic Conditions (PCCC) `vignette("PCCC", package = "icd")`
- Working with ICD code ranges `vignette("ranges", package = "icd")`
- Comparing comorbidity maps `vignette("compare-maps", package = "icd")`
- Paper detailing efficient matrix method of comorbidities `vignette("efficiency", package = "icd")`
Many users have emailed me directly for help, and I'll do what I can, but it is often better to examine or add to the list of [issues](https://github.com/jackwasey/icd) so we can help each other. Advanced users may look at the [source code](https://github.com/jackwasey/icd), particularly the extensive test suite which exercises all the key functions.
``` r
?comorbid
?comorbid_hcc
?explain_code
?is_valid
```
## ICD-9 codes
ICD-9 codes are still in heavy use around the world, particularly in the USA where the ICD-9-CM (Clinical Modification) was in widespread use until the end of 2015. ICD-10 has been used worldwide for reporting cause of death for more than a decade, and ICD-11 is due to be released in 2019. ICD-10-CM is now the primary coding scheme for US hospital admission and discharge diagnoses used for regulatory purposes and billing. A vast amount of electronic patient data is recorded with ICD-9 codes of some kind: this package enables their use in R alongside ICD-10.
ICD-9 codes are not numbers, and great care is needed when matching individual codes and ranges of codes. It is easy to make mistakes, hence the need for this package. ICD-9 codes can be presented in *short* 5 character format, or *decimal* format, with a decimal place separating the code into two groups. There are also codes beginning with V and E which have different validation rules. Zeroes after a decimal place are meaningful, so numeric ICD-9 codes cannot be used in most cases. In addition, most clinical databases contain invalid codes, and even decimal and non-decimal format codes in different places. This package primarily deals with ICD-9-CM (Clinical Modification) codes, but should be applicable or easily extendable to the original WHO ICD-9 system.
## ICD-10 codes
ICD-10 has a somewhat simpler format, with consistent use of a letter, then two alphanumeric characters. However, especially for ICD-10-CM, there are a multitude of qualifiers, e.g. specifying recurrence, laterality, which vastly increase the number of possible codes. This package recognizes validity of codes by syntax alone, or whether the codes appear in a canonical list. There is not yet the capability of converting between ICD-9 and ICD-10, but comorbidities can be generated from older ICD-9 codes and newer ICD-10 codes in parallel, and the comorbidities can then be compared.