The Double Burden of Malnutrition and Pediatric Sepsis in Tanzania

1. Defining malnutrition

1.1. Prerequisites

library(tidyverse)
library(zscorer)
library(readr)
library(dplyr)
library (lubridate)
library (ggplot2)
library (rmarkdown)
library(openxlsx)

PSFI_df_malnutrition <- read_csv("Ben_cut_4_27.csv")

  PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(
    date_present = as.POSIXct(date_present, format = "%d%b%Y:%H:%M:%S", tz = "UTC"),
    dob = mdy(dob),
    age_days_exact = as.numeric(difftime(date_present, dob, units = "days"))
  )
  
  
summary(PSFI_df_malnutrition$age_days_exact)

   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  32.03  289.19  604.61 1124.44 1461.64 5452.40 

View(PSFI_df_malnutrition)

1.2. Anthroprometric analysis

PSFI_df_malnutrition %>%
  filter(case_control == 1) %>%
  select(ht, wt, age_years, muac) %>%
  summary()

Error in `select()`:
! Can't select columns that don't exist.
✖ Column `age_years` doesn't exist.
Run `]8;;x-r-run:rlang::last_trace()rlang::last_trace()]8;;` to see where the error occurred.

# Layout to split the screen
layout(mat = matrix(c(1,2),2,1, byrow=TRUE),  height = c(1,8))
 
# Draw the boxplot and the histogram 
par(mar=c(0, 3.1, 1.1, 2.1))
boxplot(PSFI_df_malnutrition$ht[PSFI_df_malnutrition$case_control == 1] , horizontal=TRUE , ylim=c(0,200), xaxt="n" , col=rgb(0.8,0.8,0,0.5) , frame=F)
par(mar=c(4, 3.1, 1.1, 2.1))

hist(PSFI_df_malnutrition$ht[PSFI_df_malnutrition$case_control == 1]
     , breaks=40 , col=rgb(1,0.8,0.8,1) , border=F , main="" , xlab="Height (cm)", xlim=c(0,200))

layout(mat = matrix(c(1,2),2,1, byrow=TRUE),  height = c(1,8))
 par(mar=c(0, 3.1, 1.1, 2.1))
boxplot(PSFI_df_malnutrition$wt[PSFI_df_malnutrition$case_control == 1] , horizontal=TRUE , ylim=c(0,80), xaxt="n" , col=rgb(0.8,0.8,0,0.5) , frame=F)
par(mar=c(4, 3.1, 1.1, 2.1))

hist(PSFI_df_malnutrition$wt[PSFI_df_malnutrition$case_control == 1]
     , breaks=40 , col=rgb(1,0.8,0.8,1) , border=F , main="" , xlab="Weight (kg)", xlim=c(0,80))

layout(mat = matrix(c(1,2),2,1, byrow=TRUE),  height = c(1,8))
par(mar=c(0, 3.1, 1.1, 2.1))
boxplot(PSFI_df_malnutrition$muac[PSFI_df_malnutrition$case_control == 1] , horizontal=TRUE , ylim=c(5,30), xaxt="n" , col=rgb(0.8,0.8,0,0.5) , frame=F)
par(mar=c(4, 3.1, 1.1, 2.1))

hist(PSFI_df_malnutrition$muac[PSFI_df_malnutrition$case_control == 1]
     , breaks=20 , col=rgb(1,0.8,0.8,1) , border=F , main="" , xlab="Mid-upper arm circumference (cm)", xlim=c(5,30))
abline(v = 11.5, col = "red", lwd = 2, lty = 2)   # severe

abline(v = 12.5, col = "blue", lwd = 2, lty = 2)  # moderate

PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(age_years = age_days_exact / 365.25)

summary(PSFI_df_malnutrition$age_years)

    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
 0.08768  0.79176  1.65535  3.07854  4.00176 14.92786 

layout(mat = matrix(c(1,2),2,1, byrow=TRUE),  height = c(1,8))
 par(mar=c(0, 3.1, 1.1, 2.1))
boxplot(PSFI_df_malnutrition$age_years[PSFI_df_malnutrition$case_control == 1] , horizontal=TRUE , ylim=c(0,15), xaxt="n" , col=rgb(0.8,0.8,0,0.5) , frame=F)
par(mar=c(4, 3.1, 1.1, 2.1))

hist(PSFI_df_malnutrition$age_years[PSFI_df_malnutrition$case_control == 1]
     , breaks=40 , col=rgb(1,0.8,0.8,1) , border=F , main="" , xlab="Age (years)", xlim=c(0,15))

1.3. Z-scorer

PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(sex_who = if_else(sex == 1, 1, 2))

summary(PSFI_df_malnutrition$sex_who)

   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  1.000   1.000   1.000   1.427   2.000   2.000 

PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(
    age_group = case_when(
      age_days_exact < 5 * 365.25 ~ 0L,
      age_days_exact >= 5 * 365.25 ~ 1L,
      TRUE ~ NA_integer_
    )
  )

PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(
    wflz = addWGSR(
      data = .,
      sex = "sex_who",
      firstPart = "wt",
      secondPart = "ht",
      index = "wfl"
    )$wflz,
    
    wfhz = addWGSR(
      data = .,
      sex = "sex_who",
      firstPart = "wt",
      secondPart = "ht",
      index = "wfh"
    )$wfhz,
    
    baz = addWGSR(
      data = .,
      sex = "sex_who",
      firstPart = "wt",
      secondPart = "ht",
      thirdPart = "age_days_exact",
      index = "bfa"
    )$bfaz,
  )

==============================================================================================
==============================================================================================
==============================================================================================

PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(
    muacz = case_when(
      muac >= 12.5 ~ 0,
      muac >= 11.5 & muac < 12.5 ~ -2.5,
      muac < 11.5 ~ -4,
      TRUE ~ NA_real_
    )
  )

PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(
    zscore_unified = case_when(
      age_group == 0 & ht >= 45 & ht < 65 ~ wflz,
      age_group == 0 & ht >= 65 & ht < 120 ~ wfhz,
      age_group == 0 & (ht < 45 | ht >= 120 | is.na(ht)) ~ muacz,
      age_group == 1 ~ baz,
      TRUE ~ NA_real_
    )
  )

PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(
    malnutrition = case_when(
  is.na(zscore_unified) ~ NA_integer_,
  zscore_unified < -3 ~ 2L,
  zscore_unified >= -3 & zscore_unified < -2 ~ 1L,
  TRUE ~ 0L
    )
  )

PSFI_df_malnutrition <- PSFI_df_malnutrition %>%
  mutate(
    malnutrition_source = case_when(
      age_group == 0 & ht >= 45 & ht < 65  ~ "WFL",
      age_group == 0 & ht >= 65 & ht < 120 ~ "WFH",
      age_group == 0 & (ht < 45 | ht >= 120 | is.na(ht)) ~ "MUAC",
      age_group == 1 ~ "BFA",
      TRUE ~ NA_character_
    )
  )

1.4 Z score check

PSFI_df_malnutrition %>%
  filter(age_group == 1 & case_control == 1) %>%
  pull(baz) %>%
  summary()

    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
-11.0000  -2.1325  -0.4700  -0.2723   1.1300  32.6700 

PSFI_df_malnutrition %>%
  filter(age_group == 0 & ht >= 45 & ht < 65 & case_control == 1) %>%
  pull(wflz) %>%
  summary()

   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
-4.9000 -1.0075  0.8800  0.9719  2.4400 13.7700 

PSFI_df_malnutrition %>%
  filter(age_group == 0 & ht >= 65 & ht < 120 & case_control == 1) %>%
  pull(wfhz) %>%
  summary()

    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
-10.5900  -2.4175  -0.8900  -0.9841   0.7500  12.8800 

layout(mat = matrix(c(1,2),2,1, byrow=TRUE),  height = c(1,8))
 par(mar=c(0, 3.1, 1.1, 2.1))
boxplot(PSFI_df_malnutrition$baz[PSFI_df_malnutrition$age_group == 1 & PSFI_df_malnutrition$case_control == 1] , horizontal=TRUE , ylim=c(-10,10), xaxt="n" , col=rgb(0.8,0.8,0,0.5) , frame=F)
par(mar=c(4, 3.1, 1.1, 2.1))

hist(PSFI_df_malnutrition$baz[PSFI_df_malnutrition$age_group == 1 & PSFI_df_malnutrition$case_control == 1]
     , breaks=50 , col=rgb(1,0.8,0.8,1) , border=F , main="" , xlab="BMI for age (z-score)", xlim=c(-10,10))

layout(mat = matrix(c(1,2),2,1, byrow=TRUE),  height = c(1,8))
 par(mar=c(0, 3.1, 1.1, 2.1))
boxplot(PSFI_df_malnutrition$wfhz[PSFI_df_malnutrition$age_group == 0 & PSFI_df_malnutrition$ht >= 65 & PSFI_df_malnutrition$ht < 120 & PSFI_df_malnutrition$case_control == 1] , horizontal=TRUE , ylim=c(-10,10), xaxt="n" , col=rgb(0.8,0.8,0,0.5) , frame=F)
par(mar=c(4, 3.1, 1.1, 2.1))

hist(PSFI_df_malnutrition$wfhz[PSFI_df_malnutrition$age_group == 0 & PSFI_df_malnutrition$ht >= 65 & PSFI_df_malnutrition$ht < 120 & PSFI_df_malnutrition$case_control == 1]
     , breaks=50 , col=rgb(1,0.8,0.8,1) , border=F , main="" , xlab="Weight-for-height (z-score)", xlim=c(-10,10))

layout(mat = matrix(c(1,2),2,1, byrow=TRUE),  height = c(1,8))
 par(mar=c(0, 3.1, 1.1, 2.1))
boxplot(PSFI_df_malnutrition$wflz[PSFI_df_malnutrition$age_group == 0 & PSFI_df_malnutrition$ht >= 45 & PSFI_df_malnutrition$ht < 65 & PSFI_df_malnutrition$case_control == 1] , horizontal=TRUE , ylim=c(-10,10), xaxt="n" , col=rgb(0.8,0.8,0,0.5) , frame=F)
par(mar=c(4, 3.1, 1.1, 2.1))

hist(PSFI_df_malnutrition$wflz[PSFI_df_malnutrition$age_group == 0 & PSFI_df_malnutrition$ht >= 45 & PSFI_df_malnutrition$ht < 65 & PSFI_df_malnutrition$case_control == 1]
     , breaks=50 , col=rgb(1,0.8,0.8,1) , border=F , main="" , xlab="Weight for length (z-score)", xlim=c(-10,10))

## muac is not based on z-score, but rather if it is <115mm, between 115-125 or >125mm
PSFI_df_malnutrition %>%
  filter(case_control == 1) %>%
  drop_na(muacz) %>%
  mutate(
    malnutrition_status = case_when(
      muacz < -3 ~ "Severe malnutrition",
      muacz >= -3 & muacz < -2 ~ "Moderate malnutrition",
      muacz >= -2 ~ "No malnutrition",
      TRUE ~ NA_character_
    )
  ) %>%
  ggplot(aes(x = malnutrition_status)) +
  geom_bar() +
  labs(
    x = "Malnutrition status based on MUAC",
    y = "N"
  )

1.5 Data check

PSFI_df_malnutrition%>%
  group_by(malnutrition) %>%
  summarize(
    malnutrition_missing = sum(is.na(malnutrition)),
    mort_inhosp_missing = sum(is.na(mort_inhosp)), #103 mort_inhosp_miissing  = controls
  )

PSFI_df_malnutrition %>%
  filter (case_control == 1) %>%
  group_by(malnutrition_source) %>%
  summarize(
    count_malnut = n()
  )

# NA due to child with height of 20 cm (under WHO curve for wfh/wfl) and 38 days old (too young for MUAC)

PSFI_df_malnutrition %>%
  filter (case_control == 1) %>%
  group_by(malnutrition) %>%
  summarize(
    count_malnut = n()
  )

PSFI_df_malnutrition %>%
  drop_na(malnutrition, mort_inhosp) %>%
  ggplot(aes(x = factor(malnutrition))) +
  geom_bar() +
  labs(
    x = "Malnutrition",
    y = "N"
  )

PSFI_df_malnutrition%>%
  filter (case_control == 1) %>%
  count(malnutrition, mort_inhosp) %>%
  group_by(malnutrition) %>%
  mutate(prop = n / sum(n))

PSFI_df_malnutrition %>%
  drop_na(mort_inhosp, malnutrition) %>% #only cases & minus the one previously mentioned NA
  ggplot +
  (aes(x = factor(malnutrition), fill = factor(mort_inhosp))) +
  geom_bar(stat = "count", position = "fill") +
  labs(
    x = "Malnutrition",
    y = "Proportion",
    fill = "Mortality"
  )

PSFI_df_malnutrition %>%
  drop_na(malnutrition) %>%
  ggplot +
  (aes(x = factor(malnutrition), fill = factor(case_control))) +
  geom_bar(stat = "count", position = "fill") +
  labs(
    x = "Malnutrition",
    y = "Proportion",
    fill = "Case/Control"
  )

write.xlsx(PSFI_df_malnutrition, file = "PSFI_final_malnutrition.xlsx")