Multilevel (individual, county, public health region) random-intercept logit models were fit to HPV vaccination data (¡Ý1 dose Gardasil) from the 2008 Behavioral Risk Factor Surveillance System. Using the parameter estimates from the final model, we simulated 10,000 data sets for each regression coefficient from the normal distribution and applied them to the logit model to estimate HPV vaccine coverage in each county.
County-level coverage estimates ranged from 7 % to 29 % , compared with the state average of 18 % (95 % confidence interval [CI], 13.59-21.88). Many Southwestern border and metropolitan counties exhibited high coverage estimates. Low coverage estimates were noted in the Panhandle, Southeastern border region, and Northeast. Significant correlations were observed between HPV vaccination and Hispanic ethnicity, county poverty, and public health region poverty.
Harnessing the flexibility of multilevel small area models to estimate HPV vaccine coverage at the county level, we have provided data that may inform the development of health education programs/policies, the provision of health services, and the planning of new research studies. Additionally, we have provided a framework for modeling other health outcomes at the county level using national survey data.