Space-time smoothing of complex survey data: Small area estimation for child mortality

Many people living in low- and middle-income countries are not covered by civil registration and vital statistics systems. Consequently, a wide variety of other types of data, including many household sample surveys, are used to estimate health and population indicators. In this paper we combine data from sample surveys and demographic surveillance systems to produce small area estimates of child mortality through time. Small area estimates are necessary to understand geographical heterogeneity in health indicators when full-coverage vital statistics are not available. For this endeavor spatio-temporal smoothing is beneficial to alleviate problems of data sparsity. The use of conventional hierarchical models requires careful thought since the survey weights may need to be considered to alleviate bias due to nonrandom sampling and nonresponse. The application that motivated this work is an estimation of child mortality rates in five-year time intervals in regions of Tanzania. Data come from Demographic and Health Surveys conducted over the period 1991-2010 and two demographic surveillance system sites. We derive a variance estimator of under five years child mortality that accounts for the complex survey weighting. For our application, the hierarchical models we consider include random effects for area, time and survey and we compare models using a variety of measures including the conditional predictive ordinate (CPO). The method we propose is implemented via the fast and accurate integrated nested Laplace approximation (INLA).Comment: Published at http://dx.doi.org/10.1214/15-AOAS872 in the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Maximum Likelihood Estimation of Individual Inbreeding Coefficients and Null Allele Frequencies

In this paper, we developed and compared several expectation-maximization (EM) algorithms to find maximum likelihood estimates of individual inbreeding coefficients using molecular marker information. The first method estimates the inbreeding coefficient for a single individual and assumes that allele frequencies are known without error. The second method jointly estimates inbreeding coefficients and allele frequencies for a set of individuals that have been genotyped at several loci. The third method generalizes the second method to include the case in which null alleles may be present. In particular, it is able to jointly estimate individual inbreeding coefficients and allele frequencies, including the frequencies of null alleles, and accounts for missing data. We compared our methods with several other estimation procedures using simulated data and found that our methods perform well. The maximum likelihood estimators consistently gave among the lowest root-mean-square-error (RMSE) of all the estimators that were compared. Our estimator that accounts for null alleles performed particularly well and was able to tease apart the effects of null alleles, randomly missing genotypes and differing degrees of inbreeding among members of the datasets we analysed. To illustrate the performance of our estimators, we analysed previously published datasets on mice (Mus musculus) and white-tailed deer (Odocoileus virginianus)

Is the Positive Association Between Middle-Income and Rich Household Wealth and Adult Sub-Saharan African Women\u27s Overweight Status Modified by the Level of Education Attainment? A Cross-Sectional Study of 22 Countries

BACKGROUND: Previous studies show a positive association between household wealth and overweight in sub-Saharan African (SSA) countries; however, the manner in which this relationship differs in the presence of educational attainment has not been well-established. This study examined the multiplicative effect modification of educational attainment on the association between middle-income and rich household wealth and overweight status among adult females in 22 SSA countries. We hypothesized that household wealth was associated with a greater likelihood of being overweight among middle income and rich women with lower levels of educational attainment compared to those with higher levels of educational attainment. METHODS: Demographic and Health Survey (DHS) data from 2006 to 2016 for women aged 18-49 years in SSA countries were used for the study. Overweight was defined as a body mass index (BMI) ≥ 25 kg/m2. Household wealth index tertile was the exposure and educational attainment, the effect modifier. Potential confounders included age, ethnicity, place of residence, and parity. Descriptive analysis was conducted, and separate logistic regression models were fitted for each of the 22 SSA countries to compute measures of effect modification and 95% confidence intervals. Analysis of credibility (AnCred) methods were applied to assess the intrinsic credibility of the study findings and guide statistical inference. RESULTS: The prevalence of overweight ranged from 12.6% in Chad to 56.6% in Swaziland. Eighteen of the 22 SSA countries had measures of effect modification below one in at least one wealth tertile. This included eight of the 12 low-income countries and all 10 middle income countries. This implied that the odds of overweight were greater among middle-income and rich women with lower levels of educational attainment than those with higher educational attainment. On the basis of the AnCred analysis, it was found that the majority of the study findings across the region provided some support for the study hypothesis. CONCLUSIONS: Women in higher wealth strata and with lower levels of educational attainment appear to be more vulnerable to overweight compared to those in the same wealth strata but with higher levels of educational attainment in most low- and middle- income SSA countries

Evaluating the impact of programmatic mass drug administration for malaria in Zambia using routine incidence data.

BACKGROUND NlmCategory: BACKGROUND content: In 2016, the Zambian National Malaria Elimination Centre started programmatic mass drug administration (pMDA) campaigns with dihydroartemisinin-piperaquine as a malaria elimination tool in Southern Province. Two rounds were administered, two months apart (coverage 70% and 57% respectively). We evaluated the impact of one year of pMDA on malaria incidence using routine data. - Label: METHODS NlmCategory: METHODS content: We conducted an interrupted time series with comparison group analysis on monthly incidence data collected at the health facility catchment area (HFCA) level, with a negative binomial model using generalized estimating equations. pMDA was conducted in HFCAs with greater than 50 cases/1,000 people/year. Ten HFCAs with incidence rates marginally above this threshold (pMDA group) were compared to 20 HFCAs marginally below (comparison group). - Label: RESULTS NlmCategory: RESULTS content: "The pMDA HFCAs saw a 46% greater decrease in incidence at the time of intervention than the comparison areas (incidence rate ratio: 0.536 [0.337-0.852]); however, incidence increased toward the end of the season. No HFCAs saw a transmission interruption." - Label: CONCLUSION NlmCategory: CONCLUSIONS content: pMDA, implemented during one year with imperfect coverage in low transmission areas with sub-optimal vector control coverage, significantly reduced incidence. However, elimination will require additional tools. Routine data are important resources for programmatic impact evaluations and should be considered for future analyses

Effect Modification by Baseline Mortality in the MORDOR Azithromycin Trial.

We examined whether baseline mortality risk, as a function of child age and site, modified the azithromycin mortality-reduction effect in the Macrolide Oraux pour Réduire les Décès avec un Oeil sur la Résistance (MORDOR) clinical trial. We used the Cox proportional hazards model with an interaction term. Three models were examined representing three sources for the baseline-risk covariate: two using sources external to MORDOR and the third leveraging data within MORDOR. All three models provided moderate evidence for the effect becoming stronger with increasing baseline mortality (P = 0.02, 0.02, and 0.07, respectively) at the rate of approximately 6-12% additional mortality reduction per doubling of baseline mortality. Etiological and programmatic implications of these findings are discussed

Intraperitoneal drain placement and outcomes after elective colorectal surgery: international matched, prospective, cohort study

Despite current guidelines, intraperitoneal drain placement after elective colorectal surgery remains widespread. Drains were not associated with earlier detection of intraperitoneal collections, but were associated with prolonged hospital stay and increased risk of surgical-site infections.Background Many surgeons routinely place intraperitoneal drains after elective colorectal surgery. However, enhanced recovery after surgery guidelines recommend against their routine use owing to a lack of clear clinical benefit. This study aimed to describe international variation in intraperitoneal drain placement and the safety of this practice. Methods COMPASS (COMPlicAted intra-abdominal collectionS after colorectal Surgery) was a prospective, international, cohort study which enrolled consecutive adults undergoing elective colorectal surgery (February to March 2020). The primary outcome was the rate of intraperitoneal drain placement. Secondary outcomes included: rate and time to diagnosis of postoperative intraperitoneal collections; rate of surgical site infections (SSIs); time to discharge; and 30-day major postoperative complications (Clavien-Dindo grade at least III). After propensity score matching, multivariable logistic regression and Cox proportional hazards regression were used to estimate the independent association of the secondary outcomes with drain placement. Results Overall, 1805 patients from 22 countries were included (798 women, 44.2 per cent; median age 67.0 years). The drain insertion rate was 51.9 per cent (937 patients). After matching, drains were not associated with reduced rates (odds ratio (OR) 1.33, 95 per cent c.i. 0.79 to 2.23; P = 0.287) or earlier detection (hazard ratio (HR) 0.87, 0.33 to 2.31; P = 0.780) of collections. Although not associated with worse major postoperative complications (OR 1.09, 0.68 to 1.75; P = 0.709), drains were associated with delayed hospital discharge (HR 0.58, 0.52 to 0.66; P < 0.001) and an increased risk of SSIs (OR 2.47, 1.50 to 4.05; P < 0.001). Conclusion Intraperitoneal drain placement after elective colorectal surgery is not associated with earlier detection of postoperative collections, but prolongs hospital stay and increases SSI risk

Space-Time Smoothing Models for Surveillance and Complex Survey Data

Thesis (Ph.D.)--University of Washington, 2016-06Area and time-specific estimates of disease rates, cause-specific mortality rates and other key health indicators are of great interest for health care and policy purposes. Such estimates provide the information needed to identify areas with increased risk, effectively allocate resources, and target interventions. A wide variety of data, such as vital statistics, complex surveys, demographic surveillance sites, and disease registries, are used for these purposes. Unfortunately, the sample size of data available at a granular space-time scale is often too small to provide reliable estimates and uncertainty intervals. Using data from multiple sources and spatial and temporal smoothing is beneficial to alleviate problems of data scarcity. The purpose of the work described herein is to use Bayesian space-time models, to combine data from multiple sources to provide reliable area-based estimates. This work is motivated by estimating rates of health indicators (e.g. diabetes, smoking) by health reporting areas in King County from the Behavioral Risk Factor Surveillance Survey, child mortality by regions in Tanzania from Demographic and Health Surveys and demographic surveillance sites, and cancer-specific incidence and mortality rates in Europe from government data and local registries

Estimating the levels and trends of family planning indicators in 436 sub-national areas across 26 countries in sub-Saharan Africa

AbstractBackgroundScaling up access to safe, effective, and voluntary family planning (FP) services to achieve universal access for women and families will require increased commitment by countries and international organizations. On the way toward universal access, quantitative family planning goals have also been established by the United Nations through the sustainable development goals (SDGs). Here, we present a model-based framework that can help monitor progress toward these goals at the sub-national and demographic subgroup scale.MethodsWe utilize 90 demographic health surveys for 26 countries that contain associated geographic position system data. We extract survey cluster level data to fit multiple small area estimation models that estimate FP indicators by administrative unit one and two regions as well as different demographic subgroups.FindingsWe find significant variation of modern contraceptive prevalence rates (mCPR), unmet need, and demand satisfied by country, sub-national region, and demographic subgroup. Our model-based estimates show that on average for 436 administrative unit one regions, mCPR has increased at a rate of 0.75% per year and unmet need has decreased by 0.26% per year. There are also striking differences of FP indicators by demographic subgroup; for example, unmet need can be up to 40% different based on age and parity.InterpretationWe have developed a framework to help monitor progress, provide insights about the inequitable progress by region and demographic groups, and account for the sparsity of available data. These results and framework can help policy-makers better allocate and target interventions to help achieve family planning goals.FundingThis publication is based on models and data analysis performed by the Institute for Disease Modeling at the Bill &amp; Melinda Gates Foundation.</jats:sec