The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

Using Rdiagnosislist functions with custom hierarchies

Anoop Shah

2023-04-18

Using Rdiagnosislist functions with custom hierarchies

The Rdiagnosislist functions can be used with a 'fake' SNOMED environment to show hierarchical relationships with other coding systems, such as ICD-10.

In order to do this, you can modify the sample SNOMED environment by replacing the following three tables:

This example models a small section of the ICD-10 hierarchy for heart failure: - IX Diseases of the circulatory system (icdchapter) - I50 Heart failure (icd3) - I50.0 Congestive heart failure (leaf) - I50.1 Left ventricular failure (leaf) - I50.9 Heart failure, unspecified (leaf) - I11.0 Hypertensive heart disease with (congestive) heart failure (leaf) - Heart failure (phenotype)

The concepts are of the following semantic types: - leaf - a 4-character ICD-10 code which can be used for coding - icd3 - a 3-character ICD-10 category - icdchapter - lettered ICD-10 chapter - phenotype - a set of leaf codes that constitute a research phenotype

As there is only one description per concept in this example, the data can be specified by two tables: a 'description' table and a 'relations' table.

# Load package
require(Rdiagnosislist)

CUSTOM <- sampleSNOMED()

NEW_DESCRIPTION <- data.table::data.table(
  conceptId = bit64::as.integer64(1:7),
  term = c(
  'IX Diseases of the circulatory system (icdchapter)',
  'I50 Heart failure (icd3)',
  'I50.0 Congestive heart failure (leaf)',
  'I50.1 Left ventricular failure (leaf)',
  'I50.9 Heart failure, unspecified (leaf)',
  'I11.0 Hypertensive heart disease with (congestive) heart failure (leaf)',
  'Heart failure (phenotype)'
  ))
NEW_RELATIONSHIP <- data.table::data.table(
  parent = bit64::as.integer64(c(1,2,2,2,2,7,7)),
  child =  bit64::as.integer64(c(2,3,4,5,6,2,6)))

CUSTOM$CONCEPT <- NEW_DESCRIPTION[, list(id = conceptId,
  effectiveTime = as.Date(Sys.Date()), active = TRUE,
  moduleId = 0, definitionStatusId = 0)]
CUSTOM$DESCRIPTION <- NEW_DESCRIPTION[, list(id = conceptId,
  effectiveTime = as.Date(Sys.Date()), active = TRUE,
  moduleId = 0, conceptId = conceptId, languageCode = 'en',
  typeId = bit64::as.integer64('900000000000003001'), # fully specified name
  term, caseSignificanceId = 0)]
CUSTOM$STATEDRELATIONSHIP <- NEW_RELATIONSHIP[, list(id = 1:.N,
  effectiveTime = as.Date(Sys.Date()), active = TRUE,
  moduleId = 0, sourceId = child, destinationId = parent,
  relationshipGroup = 0, typeId = bit64::as.integer64('116680003'), # is a
  characteristicTypeId = 0, modifierId = 0)]
CUSTOM$RELATIONSHIP <- CUSTOM$STATEDRELATIONSHIP[0]

# Using the new dictionaries
myconcept <- SNOMEDconcept('Diseases of the circulatory system',
  exact = FALSE, SNOMED = CUSTOM)
myphenotype <- SNOMEDconcept('Heart failure (phenotype)',
  SNOMED = CUSTOM)
  
# Show the concept using the new DESCRIPTION table
description(myconcept, SNOMED = CUSTOM)
##    id conceptId                                               term
## 1:  1         1 IX Diseases of the circulatory system (icdchapter)
# Create and view a codelist based on ICD10 and the phenotype
mycodelist <- SNOMEDcodelist(c(myconcept, myphenotype),
  include_desc = TRUE, SNOMED = CUSTOM)
## Converting 2 concept(s) to a codelist
# Create HTML codelist (not run)
# temp = paste0(tempdir(), 'test.html')
# htmlCodelistHierarchy(mycodelist, file = temp, SNOMED = CUSTOM)
# system(paste0('firefox ', temp, ' &')) # open in firefox (on Linux)

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.