L'utilisation des méthodes géostatistiques reste rare dans le domaine de la santé. Afin d’évaluer l'utilité de l'approche géostatistique en épidémiologie, nous avons choisi d'appliquer ces méthodes à la couverture vaccinale contre les papillomavirus humains (HPV) en France. En effet, la couverture vaccinale contre les HPV reste faible en France et les disparités géographiques sont importantes. L'objectif de cette étude était d'identifier les facteurs socioécologiques qui peuvent expliquer ces variations géographiques.

Des données sociologiques, économiques et comportementales pour l'année 2016 ont été rassemblées (bases de données démographiques et de santé publique, web et réseaux sociaux) et ont été corrélées au taux de couverture vaccinale sur l'ensemble du territoire français. Des zones géographiques homogènes définies par des corrélations élevées pour des groupes de variables ont été sélectionnées. Dans chaque zone homogène, une analyse en composantes principales a été réalisée et une approche de modélisation géostatistique a permis d'estimer la couverture vaccinale contre le papillomavirus à l’échelle de la commune.

Un modèle unique ne permet pas d'expliquer les variations spatiales de la couverture vaccinale en France. Dans les zones urbaines, un faible taux de couverture vaccinale est préférentiellement associé à des facteurs socioéconomiques défavorables (pauvreté, chômage, immigration). Dans les zones rurales, un faible taux de couverture vaccinale est préférentiellement associé à des facteurs socioculturels (catégories socioprofessionnelles, niveau scolaire, intérêt pour les médecines alternatives et le mouvement antivaccin). Deux zones géographiques secondaires ont été définies : Île-de-France et 12 départements du quart nord-est. En Île-de-France, l'association avec les facteurs économiques est retrouvée comme dans la zone urbaine dans son ensemble. La zone nord-est est un cas particulier car le taux de couverture vaccinale y est relativement élevé en dépit d'indicateurs économiques de pauvreté.

Ce modèle géostatistique permet d'identifier de nouvelles explications potentielles des disparités territoriales de la couverture vaccinale du vaccin contre le papillomavirus en France. Ces résultats pourraient permettre d'adapter ou de développer de futurs programmes de vaccination sur des territoires spécifiques en tenant compte des caractéristiques sociologiques, économiques et comportementales de la population.

The use of geostatistical methods remains rare in health studies. In order to assess the usefulness of the geostatistical approach in epidemiology, we chose to apply these methods to the vaccination coverage rate (VCR) against human papillomavirus (HPV) in France. Indeed, HPV vaccine coverage remains low in France and geographical disparities are sizable. The objective of this study was to identify the socioecological factors that may explain these geographical variations.

Sociological, economic and behavioral data for 2016 have been gathered (demographics and public health database, web and social networks) and were correlated with the HPV VCR vaccine coverage over the French territory. Homogeneous geographical areas defined by strong correlations for groups of variables were selected. In each homogeneous area, principal component analysis was performed and a geostatistical approach provided an estimate predicting vaccine coverage at a given scale.

HPV VCR spatial variations in France cannot be fully explained by a single model. In urban areas, a low rate of HPV VCR is preferentially associated with unfavorable socioeconomic factors (poverty, unemployment, immigration). In rural areas, HPV VCR is preferentially associated with sociocultural factors (socio-professional categories, education level, interest in alternative medicines the anti-vaccine movement). Two secondary geographical areas were defined: the Île-de-France region and 12 departments in northeastern France. In the Île-de-France region, the association with the economic factors one again appears as in urban areas in general. The northeasteran departments represent a particular case insofar as HPV VCR is relatively high, notwithstanding economic poverty indicators.

Geostatistical modeling successfully identifies new potential explanations for HPV VCR geographical disparities in France. These results could help to adapt or develop future vaccination programs in specific areas by taking into account the sociological, economic and behavioral characteristics of their populations.