 Research
 Open Access
 Published:
US countylevel estimation for maternal and infant healthrelated behavior indicators using pregnancy risk assessment monitoring system data, 2016–2018
Population Health Metrics volume 20, Article number: 14 (2022)
Abstract
Background
There is a critical need for maternal and child health data at the local level (for example, county), yet most counties lack sustainable resources or capabilities to collect locallevel data. In such case, modelbased small area estimation (SAE) could be a feasible approach. SAE for maternal or infant healthrelated behaviors at small areas has never been conducted or evaluated.
Methods
We applied multilevel regression with poststratification approach to produce countylevel estimates using Pregnancy Risk Assessment Monitoring System (PRAMS) data, 2016–2018 (n = 65,803 from 23 states) for 2 key outcomes, breastfeeding at 8 weeks and infant nonsupine sleeping position.
Results
Among the 1,471 counties, the median model estimate of breastfeeding at 8 weeks was 59.8% (ranged from 34.9 to 87.4%), and the median of infant nonsupine sleeping position was 16.6% (ranged from 10.3 to 39.0%). Strong correlations were found between model estimates and direct estimates for both indicators at the state level. Model estimates for both indicators were close to direct estimates in magnitude for Philadelphia County, Pennsylvania.
Conclusion
Our findings support this approach being potentially applied to other maternal and infant health and behavioral indicators in PRAMS to facilitate public health decisionmaking at the local level.
Maternal and infant health indicators include a wide range of health outcomes and healthrelated behaviors before, during, and shortly after pregnancy. The demand for locallevel (e.g., county) data related to maternal and infant health indicators is growing among state and local health departments. Data at the local level can help improve the understanding of geographic disparities in these indicators and lead to more effective planning of program and policy strategies to promote maternal and infant health. Locallevel data can also help guide local health departments on which maternal and infant health issues most need to be addressed; state estimates may not reflect local needs. However, most maternal and infant data, particularly healthrelated behavioral indicators, are usually only available at the national or state level. Some local health departments have developed specific surveys to obtain locallevel data, such as the Philadelphia Pregnancy Risk Assessment Monitoring System (PhillyPRAMS) and the Los Angeles Mommy and Baby Project [1]; however, most counties lack sustainable resources or capabilities to collect locallevel data. In such cases, modelbased estimation for small areas could be a feasible approach.
A variety of modelbased estimation methods have been developed for small area estimation (SAE) over the past decades [2, 3], and consequently there has been a rapid increase in the use of SAE in estimating maternal and infant health indicators in recent years [4]. Most of these studies have focused on infant mortality and some infectious diseases in developing countries; most estimation was based on Bayesian hierarchical models; and estimation techniques overall are complex [4]. Estimation for healthrelated behaviors, such as breastfeeding duration and infant sleeping positions, is challenging, partly because the sources of survey data for SAE are limited, and independent sources of local data used for validation are rarely available.
In this study, we focused on 2 maternal and infant healthrelated behaviors and used a multilevel regression with poststratification (MRP) approach to estimate countylevel prevalence from the Pregnancy Risk Assessment Monitoring System (PRAMS) [5] survey data. The multilevel regression model has been increasingly used for small area estimation because it accounts for individual factors and the hierarchical nesting of data simultaneously [6,7,8,9]. MRP was originally developed by Gelman and Little to estimate statelevel voter preference using national polls [10,11,12]. The Centers for Disease Control and Prevention (CDC) expanded this approach to produce prevalence estimates for adult chronic conditions at 4 geographic levels (county, place [incorporated and Census designated], ZIP Code Tabulation Area, and census tract levels in the United States [www.cdc.gov/PLACES]) by using health data from the Behavioral Risk Factor Surveillance System [13,14,15]. The objective of this study is to estimate countylevel prevalence of maternal breastfeeding and infant nonsupine sleep position (a risk factor for sudden unexpected infant death) in the US and evaluate how well the MRP method performs.
Methods
Data source
PRAMS is an ongoing statebased surveillance system of selfreported maternal behaviors, attitudes, and experiences before, during, and shortly after pregnancy [5]. It is administered by the Division of Reproductive Health (DRH) at CDC in collaboration with state health departments. The PRAMS protocol is approved by the Institutional Review Board of the CDC and by each participating PRAMS site. From each participating site, a stratified, random sample of women with a recent live birth is selected monthly from birth certificate files. Women are surveyed 2–6 months postpartum (average = 4 months) using a standardized protocol and questionnaire. In the current study, we examined 2016–2018 PRAMS data (n = 65,803) from 23 states (a total of 1,471 counties) that agreed to provide their data for this research. The dataset included a geographic identifier of county (the number of respondents in 1,405 counties with nonzero respondents ranged from 1 to 2,335 with a median of 10, and 66 counties did not contain any respondents) for each respondent and individual level demographic and health variables. Annually, PRAMS data for each site are weighted for sampling design, nonresponse, and noncoverage to produce data representative of the site’s birth population for the year.
Two healthrelated behaviors, infant sleeping position and maternal breastfeeding duration, were assessed using 2 questions from PRAMS. Respondents were asked, “In which one position do you most often lay your baby down to sleep now?” and responses were classified into 2 groups, supine sleep position (on back) and nonsupine (on side or on stomach). Respondents were also asked several questions related to maternal breastfeeding. If they were currently (at the time of survey completion) breastfeeding (or feeding pumped milk) or, for those not breastfeeding at the time of survey completion, if they ever breastfed (or pumped breast milk) for ≥ 8 weeks or ≥ 2 months, then they were coded as any breastfeeding at 8 weeks; otherwise, they were coded as not. Missing values or responses with “do not know” were excluded for each indicator (n = 2,680 for infant sleeping position and n = 2,637 for maternal breastfeeding duration). We calculated the statelevel weighted direct estimates for both outcomes using SUDAAN (Research Triangle Institute, North Carolina).
We obtained countylevel direct estimates from PhillyPRAMS 2018, which is not a part of PRAMS but an independent survey sponsored by the Philadelphia Department of Public Health (https://www.phila.gov/departments/departmentofpublichealth). It provides estimates for Philadelphia County (which includes the city of Philadelphia) for indicators similar to those in PRAMS with a sample size of 1,489. The data were acquired by DRH through a data use agreement. Given identical questions on breastfeeding and infant sleep position on the 2 surveys, we defined infant nonsupine sleeping position and breastfeeding for 8 weeks for PhillyPRAMS in the same way as for PRAMS data and calculated their weighted direct estimates for Philadelphia County using SUDAAN (Research Triangle Institute, North Carolina).
Model construction
We constructed a multilevel logistic regression model for each of the binary outcomes, Y, respectively.
where
Y_{ijk}: Survey response of infant nonsupine sleeping position (yes or no) or breastfeeding at 8 weeks (yes or no) from respondent k in county j and state i
\({\mathrm{X}}_{\mathrm{ijk}}\): A matrix of respondent k’s demographic and socioeconomic variables.
\(\beta\): A fixed but unknown parameter vector for \({X}_{\mathrm{ijk}}\)
\({\mathrm{re}}_{j (i)}\): Random effect of county j nested in state i.
\({\mathrm{re}}_{i}\): Random effect of state i.
\({X}_{ijk}\) was a set of respondent k’s covariates. It initially included variables that were associated with Y_{ijk} based on literature review, such as maternal age, maternal race, maternal marital status, maternal education level, paternal education level, and infant sex and that were evaluated based on Pickering et al.’s empirical approach [16, 17]. Only those variables that showed significant associations (the square of the parameter estimate divided by the square of its standard error greater than 4) entered the model. The final variables included in the model were maternal education (less than high school, high school, some college, college and above), maternal Hispanic ethnicity (yes or no), and maternal race (white, black, Asian, and other; other included American Indian, Hawaiian, nonwhite other, Alaska Natives, or mixed race) for both indicators. As the countylevel random effect, \({\mathrm{re}}_{j(i)}\), was not significant (p > 0.05) in the model for infant nonsupine sleeping position, it was excluded from the final model for this outcome. If the betweenstate variance (variance in null model minus variance in final model [model includes random effects only] and then divided by variance in the null model, %) was greater than 40%, then we moved on to the next step of prediction and produced the estimates. We performed modeling using the SAS 9.4 GLIMMIX procedure (SAS Institute, Cary, North Carolina). The residual subjectspecific pseudolikelihood method was used to produce the parameter estimates \(\widehat{\beta }\) (with variance \({\widehat{\upsigma }}_{\widehat{\beta }}\)) and empirical best linear unbiased predictors, \({\widehat{\mathrm{re}}}_{j(i)}\)(with variance \({\widehat{\upsigma }}_{{\widehat{\mathrm{re}}}_{j(i)}}\)) and \({\widehat{\mathrm{re}}}_{i} (\mathrm{with variance }{\widehat{\upsigma }}_{{\widehat{\mathrm{re}}}_{i}})\). The covariance structure was specified as variance components. For the 66 counties in the 23 states included in this study without any respondents in PRAMS, we also generated their \({\widehat{\mathrm{re}}}_{j(i)}\) by averaging their neighboring \({\widehat{\mathrm{re}}}_{j(i)}\) so that we could provide estimated prevalence for all the counties.
Model parameters were then applied to live birth counts from National Center for Health Statistics birth certificate files, 2016–2018. To poststratify, we categorized the live birth counts by maternal education level, maternal Hispanic status, and maternal race for each county; thus, each county had a total of 32 (4 × 2 × 4) categories. The predicted probability (\({\widehat{p}}_{\mathrm{ijm}}\)) of breastfeeding at 8 weeks for the m^{th} category in county j, state i was calculated based on the following formula:
where X is a matrix of demographic variables, and \({X}_{\mathrm{ijm}}\) is the row of category m, and \({X}_{ijm}\widehat{\beta }={\widehat{\beta }}_{\mathrm{maternal\,education}}+{\widehat{\beta }}_{\mathrm{maternal}\,\mathrm{hispanic}}+{\widehat{\beta }}_{\mathrm{maternal}\,\mathrm{race}}\). The formula for the probability of infant nonsupine sleeping position was similar to (2) but omitting \({\widehat{\mathrm{re}}}_{j(i)}\). With \({\widehat{p}}_{ijm}\), we could calculate the estimated prevalence of the indicator through poststratification for state i and for county j as below:
where \(\sum {N}_{ijm}\) is the live birth counts of state i, and \(\sum {N}_{jm}\) is the live birth counts of county j. To take into account the uncertainty arising from the models, we adopted Monte Carlo simulation, a tool to generate sample statistics by using point estimates of model parameters and their asymptotic covariance matrix of these estimates [18] to simulate the distribution for the estimate. Thus formula (2) became the following:
where \({\widehat{\beta }}^{*}\) is a normal variate with mean \(\widehat{\beta }\) and variance \({\widehat{\upsigma }}_{\widehat{\beta }}\), \({\widehat{\mathrm{re}}}_{j(i)}^{*}\) is a normal variate with mean \({\widehat{\mathrm{re}}}_{j(i)}\) and variance \({\widehat{\upsigma }}_{{\widehat{\mathrm{re}}}_{j(i)}}\), and \({\widehat{\mathrm{re}}}_{i}^{*}\) is a normal variate with mean \({\widehat{\mathrm{re}}}_{i}\) and variance \({\widehat{\upsigma }}_{{\widehat{\mathrm{re}}}_{i}}\). Model (4) was repeated 1000 times, and then we applied each \({\widehat{p}}_{ijm}\) to formula (3), which generated 1,000 \({\widehat{P}}_{j}\). The mean and 95% confidence interval (CI, a range between the 2.5^{th} and 97.5^{th} values) were determined for breastfeeding at 8 weeks and similarly for infant nonsupine sleeping position but without \({\widehat{\mathrm{re}}}_{j(i)}^{*}\) in formula (4). In a sensitivity analysis, we fit the same model but incorporated the sampling weights from the PRAMS data. We checked the performance of the models by comparing the model estimates with weighted estimates observed from the PRAMS survey file (statelevel only) and PhillyPRAMS survey (countylevel). All the above analysis was conducted in SAS 9.4 (SAS Institute, Cary, North Carolina).
Results
Table 1 shows the state and countylevel random effects for breastfeeding at 8 weeks and statelevel random effect for infant nonsupine sleeping position in null and final models. The final model for breastfeeding at 8 weeks explained 47% of betweenstate variation and 51% of betweencounty variation. In the model for infant nonsupine sleeping position, the covariates explained 65% of the statelevel variance. The results were tested further by the comparison with direct estimates at both state and county levels.
In Table 2, we listed 2 sets of estimates for the indicators by state. The median model estimate of breastfeeding at 8 weeks was 67.8% and ranged from 46.8 to 82.0% among 23 states. Model estimates were very close to the direct survey estimates for each state. The Pearson correlation coefficient (Pearson’s r) between the 2 sets of estimates was 93%. The 95% CIs of model estimates were wider than those of the corresponding direct estimates. The median model estimate of infant nonsupine sleeping position was 18.4%, ranged from 12.6 to 29.9%, and was very close to direct estimates for each state as well, but the 95% CIs were narrower than the direct estimates. The Pearson’s r between the 2 sets of estimates was 93%. We also observed that the 95% confidence intervals of maternal breastfeeding duration state estimates were wider than those of infant sleeping position.
Among the 1,471 counties from 23 states, the median modelbased estimate of breastfeeding at 8 weeks was 59.8% with a range from 34.9 to 87.4%, and the median of infant nonsupine sleeping position was 16.6% with a range from 10.3 to 39.0%. Figure 1 shows ordered model estimates and their 95% CIs of each outcome. The 95% CIs for maternal breastfeeding duration county estimates were found to be wider than those for infant sleeping position. In Philadelphia County (Table 2), the estimate for breastfeeding at 8 weeks was 58.8% (model, 95% CI: 53.1%, 64.1%) versus 65.9% (direct, 95% CI: 61.6%, 70.0%), and the estimate for infant nonsupine sleeping position was 24.0% (model, 95% CI: 22.1%, 25.9%) versus 23.9% (direct, 95% CI: 20.3%, 28.0%). The 2 sets of estimates’ 95% CIs overlapped for both indicators.
Figure 2 illustrates the countylevel geographic distribution of the 2 indicators among 23 states. The top map shows model estimates of breastfeeding at 8 weeks. For most of the states, the variations in prevalence among counties were small. For example, the prevalence of maternal breastfeeding at 8 weeks in Oregon was relatively high across the counties and was low across the counties in Louisiana, and the same patterns were observed in almost all the counties in the respective states. In some other states, such as Alabama, Georgia, and South Carolina, variations in the prevalence of breastfeeding at 8 weeks were seen across the counties. For example, countylevel estimates varied from 38.7 to 64.2% in Alabama. Similar patterns were observed for infant nonsupine sleeping position. In the sensitivity analysis, we used the sampling weights when modelling the data. Differences between the weighted and unweighted estimates were within 1% for 98% of the counties and thus were not reported.
Discussion
In this study, we generated estimates of 2 maternal and infant healthrelated behaviors for counties from 23 states via an MRP approach and used direct estimates at both state and county levels to assess the validity of the model estimates.
Good performance of MRP relies on multiple factors, such as the prevalence of the indicator in the population, data sources, model specification, availability of matching populations, and quality of the data used for validation. One of the major challenges of MRP is how to specify the multilevel regression model adequately and accurately. There are no “gold standard” criteria on variable selection. Methodologists suggest including all variables that have an important impact on sampling and nonresponse and are also potentially predictive of the outcome of interest [19, 20]. However, increasing the number of variables in the model does not necessarily increase the accuracy of the model estimates, because adding variables may cause nonconvergence or overfitting problems [21]. The variables considered in the model should be not only correlated with the outcome, but also able to explain much of the arealevel variance [16]. The approach we adopted, which was based on Pickering et al.’s empirical experience, was straightforward and performed well according to the validation analysis.
Another challenge of MRP is assessing the uncertainty around the model estimates. Because bias and variance cannot be derived in closed form [22], the measure of uncertainty relies on resampling techniques. In Bayesian estimation a posterior distribution of small area estimates is produced, however, large datasets with complex sampling can present convergence and computational challenges [3]. Therefore, frequentist approaches, such as bootstrapping and Monte Carlo simulation [18, 23], are also used to estimate mean square error or construct 95% confidence intervals. Maternal breastfeeding estimates had somewhat wider confidence intervals than infant sleeping position estimates. Besides different outcomes and different specification of models, the extent of the betweenarea variation that remained unexplained in the model also played a role [24]. The confidence intervals can be used by local health departments that may wish to compare rates and estimated precision across counties, or to evaluate subgroup differences on maternal and infant health indicators [21].
External validation is typically done by comparing similarities in the distribution and correlations between model and direct estimates rather than the absolute values because multiple factors may contribute to systematic differences in absolute values. However, independent, reliable, and concurrent local surveys are rare, which makes external validation a big challenge in small area estimation. In this study, PhillyPRAMS data was available for external data validation, and absolute values and overlapping CIs were used for comparison purposes. The model estimates were quite consistent with the direct estimates for the county of Philadelphia, which is primarily due to 2 reasons. First, the variables included in the model captured important covariates for maternal and infant outcomes, which yielded good estimation results. Second, PhillyPRAMS is similar to PRAMS, including identical questions and weighting. If more countylevel samples are drawn by local districts, future studies might be able to offer additional validation evidence of this methodology across a range of outcomes.
Other limitations should also be noted. First, the small area estimation process does not serve to reduce any nonsampling errors in the survey data, such as nonresponse or recall biases; to the extent they exist these errors remain within the modeling process. Second, although we pooled 3 survey years of data, the sample size in some counties remained small (< 10), which might pose a challenge in modelling. We are not clear how these small sample sizes influenced the predictions of the countylevel random effects and subsequently the model estimates. However, a study found that although increases in sample size help improve precision of the model estimates, including relevant covariates helps even more as sample sizes become very small [25]. Further study is warranted to better understand the influence of small sample sizes related to the categorization of the random effect variable. Finally, we did not account for geographic covariates, such as health care access or hospital density, as some other studies have done when estimating maternal and infant health indictors [5]. Although there was little variation among counties after demographic and socioeconomic variables were introduced, if a strong statelevel covariate was included, the model performance might be improved.
The multilevel model estimation process implemented in this study has important public health implications. While our estimation process did not explicitly include policy or program intervention effects, the analysis of county and statelevel random effects could provide insights into potential policy and practice implications. For example, infant sleeping position tended to have more variability at the state level than the county level, suggesting that policies, programs, or practices applied statewide may be sufficient to meet the needs of populations at higher levels of risk. Maternal breastfeeding may potentially be affected by a wide range of local factors, such as lactation service or support, given there was relatively more heterogeneity among counties than there was for infant sleeping position; this suggests that these data may be more helpful for determining where local interventions are most needed.
Conclusion
Our findings support this approach being potentially applied to other maternal and infant health and behavioral indicators in PRAMS for estimates at substate levels, such as county, ZIP Code, and census tract; however, the model specification should be tailored specifically to the outcome of interest.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from CDC on reasonable request.
Abbreviations
 CDC:

Centers for Disease Control and Prevention
 CI:

Confidence interval
 DRH:

Division of Reproductive Health
 MRP:

Multilevel regression with poststratification
 PhillyPRAMS:

Philadelphia Pregnancy Risk Assessment Monitoring System
 PRAMS:

Pregnancy Risk Assessment Monitoring System
 SAE:

Small area estimation
References
Archuleta J, Chao SM. Maternal characteristics that impact postpartum weight retention: results from the 2016 Los Angeles Mommy and Baby (LAMB) followup study. Matern Child Health J. 2021;25:151–61.
Ghosh M, Rao JNK. Small area estimation: an appraisal. Stat Sci. 1994;9:65–93.
Pfeffermann D. New important developments in small area estimation. Stat Sci. 2013;28:40–68.
Ferreira LZ, Blumenberg C, Utazi CE, et al. Geospatial estimation of reproductive, maternal, newborn and child health indicators: a systematic review of methodological aspects of studies based on household surveys. Int J Health Geogr. 2020;19:41.
Shulman HB, D’Angelo DV, Harrison L, Smith RA, Warner L. The Pregnancy risk assessment monitoring system (PRAMS): overview of design and methodology. Am J Public Health. 2018;108:1305–13.
Congdon P. A multilevel model for cardiovascular disease prevalence in the US and its application to micro area prevalence estimates. Int J Health Geogr. 2009;8:6.
Berkowitz Z, Zhang X, Richards TB, et al. Multilevel smallarea estimation of multiple cigarette smoking status categories using the 2012 behavioral risk factor surveillance system. Cancer Epidemiol Biomark Prev. 2016;25:1402–10.
Taylor J, Moon G, Twigg L. Using geocoded survey data to improve the accuracy of multilevel small area synthetic estimates. Soc Sci Res. 2016;56:108–16.
Knutson K, Zhang W, Tabnak F. Applying the smallarea estimation method to estimate a population eligible for breast cancer detection services. Prev Chronic Dis. 2008;5(1):A10.
Gelman A, Little TC. Poststratification into many categories using hierarchical logistic regression. Surv Methodol. 1997;23:127–35.
Park DK, Gelman A, Bafumi J. Bayesian multilevel estimation with poststratification: statelevel estimates from national polls. Polit Anal. 2004;12:375–85.
Park DK, Gelman A, Bafumi J. Statelevel opinions from national surveys: Poststratification using multilevel logistic regression. In: Cohen JE, editor. Public opinion in state politics. Palo Alto: Stanford University Press; 2006.
Zhang X, Holt JB, Lu H, et al. Multilevel regression and postratification for smallarea estimation of population health outcomes: a case study of chronic obstructive pulmonary disease prevalence using the behavioral risk factor surveillance system. Am J Epidemiol. 2014;179:1025–33.
Zhang X, Holt JB, Yun S, et al. Validation of multilevel regression and postratification methodology for small area estimation of health indicators from the behavioral risk factor surveillance system. Am J Epidemiol. 2015;182:127–37.
Wang Y, Holt JB, Zhang X, et al. Comparison of methods for estimating prevalence of chronic diseases and health behaviors for small geographic areas: Boston validation study, 2013. Prev Chronic Dis. 2017;14:170281.
Pickering K, Scholes S, Bajekal M. Synthetic estimation of healthy lifestyles indicators: Stage 2 report. https://www.publichealth.ie/files/file/Synthetic_Estimation_Stage_2_Report.pdf. (Accessed June 29, 2021)
Pickering K, Scholes S, Bajekal M. Synthetic estimation of healthy lifestyles indicators: Stage 3 report. https://www.publichealth.ie/files/file/Synthetic_Estimation_Stage_3_Report.pdf. (Accessed June 29, 2021)
Preacher K. Advantages of Monte Carlo confidence intervals for indirect effects. Commun Methods Meas. 2012;6:77–98.
Gelman A. Struggles with survey weighting and regression modeling. Stat Sci. 2007;22:153–64.
Downes M, Carlin JB. Multilevel regression and poststratification as a modeling approach for estimating population quantities in large population health studies: a simulation study. Biom J. 2020;62:479–91.
Haslett S, Jones G. Small area estimation using surveys and censuses: some practical and statistical issues. Stat Transit. 2005;7:541–55.
Skrondal A, RabeHesketh S. Prediction in multilevel generalized linear models. J R Stat Soc A. 2009;172:659–87.
Chatterjee S, Lahiri P, Li H. Parameteric bootstrap approximation to the distribution of EBLUP and related prediction intervals in linear mixed models. Ann Stat. 2008;36:1221–45.
Bajekal M, Scholes S, Pickering K, Purdon S. Synthetic estimation of healthy lifestyle indicators: Stage 1 report. https://www.publichealth.ie/files/file/Synthetic_Estimation_Stage_1_Report.pdf. (Accessed June 29, 2021).
Srebotnjak T, Mokdad AH, Murray CJL. A novel framework for validating and applying standardized small area measurement strategies. Popul Health Metr. 2010;8:26.
Acknowledgements
We acknowledge Vanessa Mensah, formerly of the Philadelphia Department of Public Health, and Sara Thuma, Pennsylvania Department of Health.
Disclaimer
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
YW designed this study, analyzed data, and drafted the manuscript. HT, HL, SC, SC, RL and WB participated in the study design, directed the implementation, and reviewed drafts. HL made the maps. HS, BM, and PH helped prepare the data and reviewed the draft. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare they have no conflict of interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wang, Y., Tevendale, H., Lu, H. et al. US countylevel estimation for maternal and infant healthrelated behavior indicators using pregnancy risk assessment monitoring system data, 2016–2018. Popul Health Metrics 20, 14 (2022). https://doi.org/10.1186/s12963022002916
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12963022002916
Keywords
 Infant nonsupine sleeping position
 Maternal breastfeeding
 Multilevel regression
 Pregnancy risk assessment monitoring system
 Small area estimation