Effectiveness of Interventions for Improving HPV Vaccination Cover

The inclusion criteria are: public and private educational institutions recognized by the Ministry of National Education, located in municipalities with a reported HPV vaccination coverage rate of less than 25% (using the national 2021 vaccination coverage statistics as a baseline). Parents or guardians, the educational community, and children attending the selected institution will participate. Children must be between the ages of 9 and 17 years old and boys between 9 and 14 years old, in accordance with the guidelines of the Ministry of Health and Social Protection. To estimate the sample size, the method suggested by Campbell and van Brukelen was implemented, which requires the participation of 12 groups of 100 participants (n=1,200). For each municipality, schools will be selected using the official list of educational institutions provided by the Ministry of Education as a sampling framework. Simple random sampling will be used. The randomization and analysis units will be educational institutions randomly assigned in a 1:1:1 ratio to the control or intervention group. All study arms will receive psychosocial training with training from health professionals and teachers from the selected municipalities, defined within the community setting. Within the school-based intervention, the behavioral component will use two strategies implemented with an approach to biopsychosocial models: the health belief model, the theory of planned behavior, social cognitive theory, or the theory of reasoned action. The outcomes to be measured will be: the proportion of participants (boys and girls aged 9-17 years and boys aged 9-14 years, in accordance with the guidelines of the Ministry of Health and Social Protection) exposed to the intervention who receive a dose of vaccine during the first 30 days after the intervention, the degree of acceptability of vaccination, and the effectiveness of the training strategies in increasing HPV vaccination coverage. Statistical analysis will be performed following the "Intent-to-Treat (ITT)" principle. The primary outcome will be assessed between the intervention and control groups using a two-dimensional multilevel model, with participants nested within each educational institution. The model will be adjusted for covariates at the individual level (Level 1: student sex and age) and the educational setting (Level 2: nature of the educational institution, religious education, and geographic area). Data will be analyzed using R software version 4.4.1 (RStudio 2021.09.0).

Expected Impact The development of this project will allow for the design and implementation of a program to improve vaccination coverage with the quadrivalent human papillomavirus (HPV) vaccine in Colombia.

Analysis plan: the primary outcome will be assessed during the first 30 days after administering the intervention, and statistical analysis will be performed following the "Intent-to-Treat (ITT)" principle. The intention-to-treat population is defined as all participants registered for the training intervention, regardless of their compliance. The primary outcome will be measured during the first 30 days after the intervention, emphasizing that receiving at least one dose of HPV vaccination will be considered the outcome of interest. The population exposed to the intervention who has not received the vaccine at the end of the observation period will be followed up by a member of the research group, who will determine, through home visits, whether the non-vaccination was due to lack of consent or assent, or because they were deemed no longer needed (e.g., initiation of sexual activity, teenage pregnancy). If home contact is not feasible, telephone and written (email) communication will be established to determine whether the non-vaccination was a consequence of migration or displacement from the territory. If contact is not feasible, loss to follow-up will be considered. Sensitivity analyses will be performed to assess the robustness of the conclusions based on the percentage of losses to follow-up.

The primary outcome will be assessed in each intervention group and then compared to the control group using a two-dimensional multilevel model, with participants nested within each cluster (educational institution). The model will be adjusted for covariates at the individual level (level 1: student sex and age) and the educational setting (level 2: nature of the educational institution, religious education, and geographic area). Data will be analyzed using R software version 4.4.1 (RStudio 2021.09.0). Descriptive statistics will be applied to the sociodemographic variables. Measures of central tendency and dispersion will be estimated for continuous data, and absolute and relative frequencies will be used for qualitative data. Frequency and distribution will be examined for categorical variables, and normal distribution and homogeneity of variance will be analyzed for continuous variables. For categorical variables, the differences between each intervention group and the control group will be analyzed using Fisher's exact test or Pearson's chi-square test. Continuous variables will be compared between each intervention group and the control group using the Mann-Whitney test or Student's t test, depending on the presence or absence of normal distribution of the data. Multivariate logistic regression will be implemented to adjust the association of interest for the presence of the previously mentioned level one and two covariates.

A Bonferroni correction will be performed for the level of significance. Confidence intervals will be estimated, and the adjusted and unadjusted odds ratios (OR and OR) will be implemented as measures of association. The study protocol must be previously approved by the Ethics Committee of the National Cancer Institute as the executing entity.