Health literacy (HLIT) is increasingly recognized as a critical public health capability that mediates individuals’ capacity to access, comprehend, and act on health information. Not only is HLIT important on the individual level, but it also affect...
Health literacy (HLIT) is increasingly recognized as a critical public health capability that mediates individuals’ capacity to access, comprehend, and act on health information. Not only is HLIT important on the individual level, but it also affects population health outcomes and is related to health system efficiencies. However, prevailing research lacks a systematic method to analyze HLIT using population-level survey data and has often treated HLIT as a fixed individual attribute, overlooking its embeddedness within structural and contextual conditions. This study addresses this gap by introducing the Population Health Literacy Assessment through Multilevel Estimation (PHLAME) framework—an integrative analytic strategy designed to model HLIT as a socially stratified and spatially patterned outcome. Drawing on a nationally representative HLIT survey (N = 11,027), the study employs a sequential analytic design combining latent profile analysis, small-area prediction, and multilevel modeling. Individual HLIT scores were downscaled from metropolitan to district and county level using covariate-informed estimation to generate high-resolution spatial predictions. Area-level socioeconomic contexts were classified via latent profile analysis (LPA) of structural indicators (e.g., age dependency, basic livelihood support, educational attainment), yielding four district typologies: Affluent, Moderately Deprived, Severely deprived, and Average SES. Multilevel linear models nested individuals within 255 administrative districts to estimate the independent and interactive effects of place-based deprivation. Intraclass correlation analysis confirmed that 6.6% of HLIT variance was attributable to contextual differences. Incorporating LPA-defined SES profiles significantly improved model fit (ΔAIC = 4477, p < .001) and adding individual-level predictors and cross-level interactions further enhanced explanatory power (ΔAIC = 20, p < .001). Education and employment emerged as strong positive predictors of HLIT, while older age, female sex, and disability were negatively associated. Importantly, the effect of individual education was more pronounced in deprived districts, suggesting a structurally contingent return on personal resources. Findings provide robust evidence that HLIT is co-produced by individual characteristics and the broader socioeconomic environments in which people live. The PHLAME framework offers a scalable and policy-relevant template for mapping HLIT equity, identifying contextual leverage points, and guiding targeted, place-sensitive interventions to reduce HLIT disparities at the population level.