A population-based study of sociodemographic and geographic variation in hpv vaccination

Background: Human papillomavirus (HPV) vaccination rates in the United States remain low and lag behind other recommended adolescent vaccines. Studies evaluating the association of geographic and area-level characteristics with HPV vaccination rates provide a valuable resource for public health planning. Method: We used the Rochester Epidemiology Project data linkage system to ascertain HPV vaccination rates between 2010 and 2015 in a 7-county region of southern Minnesota. Geocoded individual patient data were spatially linked to socioeconomic data from the American Community Survey at the census block group level. Bayesian hierarchical logistic regression was used to model incident vaccination rates, adjusting for individual-and area-level sociodemographic characteristics, and geolocation. Geolocation was modeled as an approximated Gaussian field using a Stochastic Partial Differential Equations approach. All models were estimated using Integrated Nested Laplace Approximations. Results: In adjusted models, increasing age and female sex were associated with increased HPV vaccination. Lower socioeconomic status was associated with decreased rates of initiation [adjusted odds ratio (AOR); 95% confidence interval 0.90 (0.86-0.95)], completion of the second dose [AOR 0.88 (0.83-0.93)], and completion of the third dose [AOR 0.85 (0.80-0.92)]. Geographic spatial analysis demonstrated increased odds of vaccination for the eastern region and in the greater Rochester metropolitan area, showing significant spatial variation not explained by individual level characteristics and ACS block group-level data. Conclusions: HPV vaccination rates varied geographically and by individual and geographically indexed sociodemographic characteristics. Impact: Identifying geographic regions with low HPV vaccination rates can help target clinical and community efforts to improve vaccination rates.