data-raw/R/make_charts_preterm.R

In growthcharts/chartbox: Collection of Growth Charts

# make_charts_preterm.R
#
# Writes selected preterm growth chart series

library(chartdesigner)
library(chartbox)

preterms <- c(
  paste0("PJAAN", 25:36), paste0("PJABN", 25:36),
  paste0("PMAAN", 25:36), paste0("PMABN", 25:36),
  paste0("PJEAN", 25:36),
  paste0("PMEAN", 25:36))

solo <-
  c(paste0("PJAON", 25:36), paste0("PMAON", 25:36),
    paste0("PJAWN", 25:36), paste0("PMAWN", 25:36),
    paste0("PJAHN", 25:36), paste0("PMAHN", 25:36),
    paste0("PJEHN", 25:36), paste0("PMEHN", 25:36),
    paste0("PJEWN", 25:36), paste0("PMEWN", 25:36),
    paste0("PJADN", 25:36), paste0("PMADN", 25:36),
    paste0("PJEDN", 25:36), paste0("PMEDN", 25:36))

sides <- c(preterms, solo)

# # Update only selected charts for design E
# sides <- c(paste0("PJEAN", 25:36), paste0("PMEAN", 25:36),
#            paste0("PJEHN", 25:36), paste0("PMEHN", 25:36),
#            paste0("PJEWN", 25:36), paste0("PMEWN", 25:36))
#            # paste0("PJEDN", 25:36), paste0("PMEDN", 25:36))

# write to chartbox
project <- path.expand("~/package/chartbox/chartbox")
chartbox <- file.path(project, "inst/library/preterm")
for (chartcode in sides) {
  outfile <- file.path(chartbox, paste(chartcode, "rds", sep = "."))
  g <- growthchart(chartcode, paper = "A4")
  saveRDS(g, file = outfile, compress = "xz")
}

# check - reload every file
files <- dir(chartbox)
chartcodes <- unlist(strsplit(files, split = "[.]"))[seq(1, 2 * length(files), 2)]

# check - which sides do not have files?
sides[!sides %in% chartcodes]

# export chart to pdf
pdfdir <- file.path(project, "pdf")
pdf.options(useDingbats = FALSE, colormodel = "srgb", family = "Helvetica")

setpalet <- function(chartcode) {
  palettes <- chartbox::palettes
  parsed <- chartcatalog::parse_chartcode(chartcode)
  palette(palettes[parsed$population,])
}

for (chartcode in chartcodes) {
  infile <- file.path(chartbox, paste(chartcode, "rds", sep = "."))
  g <- readRDS(file = infile)
  filename <- file.path(pdfdir, "preterm", paste(chartcode, "pdf", sep = "."))
  pdf(filename, height = 29.7/2.54, width = 21/2.54)
  setpalet(chartcode)
  grid::grid.draw(g)
  dev.off()
}