Results

This report summarizes results from the current retained analysis dataset.

Table 1. Distribution of doses by delivery modality

knitr::kable(
  table1,
  align = c("l", "r", "r", "r", "r", "r", "r"),
  caption = "Distribution of first doses, second doses, and all doses by delivery modality."
)
Distribution of first doses, second doses, and all doses by delivery modality.
modality First Doses First Doses % Second Doses Second Doses % All Doses All Doses %
School 3081 33% 2333 34% 5414 33%
Health facility 5873 63% 4188 61% 10061 62%
Outreach 409 4% 291 4% 700 4%
Total 9363 6812 16175

Table 2. Within-girl patterns of care seeking for second doses

knitr::kable(
  table2,
  align = c("l", "r", "r"),
  caption = "Within-girl patterns of care seeking for girls who received a second HPV vaccine dose."
)
Within-girl patterns of care seeking for girls who received a second HPV vaccine dose.
Pattern of care seeking Simple average across sites Weighted by total second doses
Both doses in school 23.88% 25.50%
Mixed location 34.29% 31.92%
Both doses out of school 41.84% 42.58%

Figure 1. Distribution of dose 1 delivered in schools

suppressPackageStartupMessages(library(ggplot2))

ggplot(plot_data, aes(x = dose1_school_prop_pct)) +
  geom_histogram(binwidth = 5, boundary = 0, closed = "left", color = "white", fill = "#2c7fb8") +
  labs(
    x = "Dose 1 delivered in schools (%)",
    y = "Number of health areas",
    title = "Distribution of proportion of dose 1 delivered in schools",
    subtitle = "Across retained health-area records (n = 20)"
  ) +
  theme_minimal(base_size = 12)
plot of chunk fig-1
plot of chunk fig-1

Table 3. Proportion of dose 1 delivered in schools by health area

table3 <- plot_data %>%
  transmute(
    `Health Area` = health_area_name,
    `Dose 1 in schools (%)` = round(dose1_school_prop_pct, 2)
  ) %>%
  arrange(desc(`Dose 1 in schools (%)`))

knitr::kable(
  table3,
  align = c("l", "r"),
  caption = "One row per health area, sorted in descending order by proportion of dose 1 delivered in schools."
)
One row per health area, sorted in descending order by proportion of dose 1 delivered in schools.
Health Area Dose 1 in schools (%)
Guápiles - Limón 99.55
Curridabat - San José 67.49
San Rafael de Heredia - Heredia 47.95
Santa Cruz - Guanacaste 47.50
Liberia - Guanacaste 34.38
Heredia Cubujuquí - Heredia 32.99
Aguas Zarcas - Alajuela 30.24
Zapote-Catedral - San José 28.12
La Unión - Cartago 26.92
Ciudad Quesada - Alajuela 26.01
Santa Ana - San José 25.50
Escazú - San José 15.56
Cartago - Cartago 13.54
Oreamuno - Cartago 13.30
Limón - Limón 10.25
Siquirres - Limón 10.12
Montes de Oro - Puntarenas 9.45
San Juan San Diego Concepción - Cartago 7.14
Desamparados 3 - San José 4.79
Turrialba - Cartago 1.41

Notes

  • Table 1 is based on 20 retained health-area records.
  • Table 2 simple averages weight each health area equally.
  • Table 2 weighted estimates use total second doses as the weight.
  • Total second-dose weight for Table 2 is 6,812.