Matta et al., 10.1001/jamainternmed.2021.6454, Table 2

Matta-T2.R

# Matta et al., 10.1001/jamainternmed.2021.6454

# Table 2

n.s0b0 <- 25271
n.s0b1 <- 461
n.s1b0 <- 638
n.s1b1 <- 453

symptoms <- list(
  list(name="Sleep problems", nn=c(2580, 49, 55, 45))
, list(name="Joint pain", nn=c(1802, 30, 26, 36))
, list(name="Back pain", nn=c(1525, 32, 33, 40))
, list(name="Digestive", nn=c(838, 33, 20, 18))
, list(name="Muscular", nn=c(808, 22, 18, 19))
, list(name="Fatigue_", nn=c(625, 57, 22, 62))
, list(name="Attention", nn=c(555, 34, 17, 38))
, list(name="Skin problems", nn=c(598, 17, 6, 11))
, list(name="Other symptoms", nn=c(463, 17, 8, 26))
, list(name="Sensory_", nn=c(463, 16, 8, 5))
, list(name="Hearing_", nn=c(456, 7, 6, 10))
, list(name="Headache", nn=c(323, 13, 8, 16))
, list(name="Breathing", nn=c(192, 29, 9, 26))
, list(name="Palpitations", nn=c(175, 17, 6, 15))
, list(name="Dizziness", nn=c(158, 7, 5, 8))
, list(name="Chest pain", nn=c(138, 14, 2, 20))
, list(name="Cough___", nn=c(144, 10, 2, 11))
, list(name="Anosmia_", nn=c(75, 20, 7, 44))
)

for (s in symptoms) {
  mat <- matrix(c(
    s$nn[1], n.s0b0 - s$nn[1]
  , s$nn[2], n.s0b1 - s$nn[2]
  , s$nn[3], n.s1b0 - s$nn[3]
  , s$nn[4], n.s1b1 - s$nn[4]
  ), ncol=4)
  p <- chisq.test(mat)$p.value

  cat(s$name
    , s$nn[1], paste0("(", sprintf("%.1f", (s$nn[1] * 100) / n.s0b0), ")")
    , s$nn[2], paste0("(", sprintf("%.1f", (s$nn[2] * 100) / n.s0b1), ")")
    , s$nn[3], paste0("(", sprintf("%.1f", (s$nn[3] * 100) / n.s1b0), ")")
    , s$nn[4], paste0("(", sprintf("%.1f", (s$nn[4] * 100) / n.s1b1), ")")
    , sprintf("%.3f", p), "\n", sep="\t")
}