Functional Mixed Effects Model for Small Area Estimation.
Ontology highlight
ABSTRACT:
SUBMITTER: Maiti T
PROVIDER: S-EPMC5079448 | biostudies-literature | 2016 Sep
REPOSITORIES: biostudies-literature
Publications
Functional Mixed Effects Model for Small Area Estimation.
Maiti Tapabrata T Sinha Samiran S Zhong Ping-Shou PS
Scandinavian journal of statistics, theory and applications 20160315 3
Functional data analysis has become an important area of research due to its ability of handling high dimensional and complex data structures. However, the development is limited in the context of linear mixed effect models, and in particular, for small area estimation. The linear mixed effect models are the backbone of small area estimation. In this article, we consider area level data, and fit a varying coefficient linear mixed effect model where the varying coefficients are semi-parametricall ...[more]
Similar Datasets
Project description:Obtaining reliable estimates about health outcomes for areas or domains where only few to no samples are available is the goal of small area estimation (SAE). Often, we rely on health surveys to obtain information about health outcomes. Such surveys are often characterised by a complex design, stratification, and unequal sampling weights as common features. Hierarchical Bayesian models are well recognised in SAE as a spatial smoothing method, but often ignore the sampling weights that reflect the complex sampling design. In this paper, we focus on data obtained from a health survey where the sampling weights of the sampled individuals are the only information available about the design. We develop a predictive model-based approach to estimate the prevalence of a binary outcome for both the sampled and non-sampled individuals, using hierarchical Bayesian models that take into account the sampling weights. A simulation study is carried out to compare the performance of our proposed method with other established methods. The results indicate that our proposed method achieves great reductions in mean squared error when compared with standard approaches. It performs equally well or better when compared with more elaborate methods when there is a relationship between the responses and the sampling weights. The proposed method is applied to estimate asthma prevalence across districts.
| S-EPMC5627524 | biostudies-literature
Project description:Bayesian hierarchical regression (BHR) is often used in small area estimation (SAE). BHR conditions on the samples. Therefore, when data are from a complex sample survey, neither survey sampling design nor survey weights are used. This can introduce bias and/or cause large variance. Further, if non-informative priors are used, BHR often requires the combination of multiple years of data to produce sample sizes that yield adequate precision; this can result in poor timeliness and can obscure trends. To address bias and variance, we propose a design assisted model-based approach for SAE by integrating adjusted sample weights. To address timeliness, we use historical data to define informative priors (power prior); this allows estimates to be derived from a single year of data. Using American Community Survey data for validation, we applied the proposed method to Behavioral Risk Factor Surveillance System data. We estimated the prevalence of disability for all U.S. counties. We show that our method can produce estimates that are both more timely than those arising from widely-used alternatives and are closer to ACS' direct estimates, particularly for low-data counties. Our method can be generalized to estimate the county-level prevalence of other health related measurements.
| S-EPMC7182002 | biostudies-literature
Project description:Small area estimation (SAE) entails estimating characteristics of interest for domains, often geographical areas, in which there may be few or no samples available. SAE has a long history and a wide variety of methods have been suggested, from a bewildering range of philosophical standpoints. We describe design-based and model-based approaches and models that are specified at the area-level and at the unit-level, focusing on health applications and fully Bayesian spatial models. The use of auxiliary information is a key ingredient for successful inference when response data are sparse and we discuss a number of approaches that allow the inclusion of covariate data. SAE for HIV prevalence, using data collected from a Demographic Health Survey in Malawi in 2015-2016, is used to illustrate a number of techniques. The potential use of SAE techniques for outcomes related to COVID-19 is discussed.
| S-EPMC9451141 | biostudies-literature
Project description:BackgroundModel-based small area estimation methods can help generate parameter estimates at the district level, where planned population survey sample sizes are not large enough to support direct estimates of HIV prevalence with adequate precision. We computed district-level HIV prevalence estimates and their 95% confidence intervals for districts in Uganda.MethodsOur analysis used direct survey and model-based estimation methods, including Fay-Herriot (area-level) and Battese-Harter-Fuller (unit-level) small area models. We used regression analysis to assess for consistency in estimating HIV prevalence. We use a ratio analysis of the mean square error and the coefficient of variation of the estimates to evaluate precision. The models were applied to Uganda Population-Based HIV Impact Assessment 2016/2017 data with auxiliary information from the 2016 Lot Quality Assurance Sampling survey and antenatal care data from district health information system datasets for unit-level and area-level models, respectively.ResultsEstimates from the model-based and the direct survey methods were similar. However, direct survey estimates were unstable compared with the model-based estimates. Area-level model estimates were more stable than unit-level model estimates. The correlation between unit-level and direct survey estimates was (β1 = 0.66, r2 = 0.862), and correlation between area-level model and direct survey estimates was (β1 = 0.44, r2 = 0.698). The error associated with the estimates decreased by 37.5% and 33.1% for the unit-level and area-level models, respectively, compared to the direct survey estimates.ConclusionsAlthough the unit-level model estimates were less precise than the area-level model estimates, they were highly correlated with the direct survey estimates and had less standard error associated with estimates than the area-level model. Unit-level models provide more accurate and reliable data to support local decision-making when unit-level auxiliary information is available.
| S-EPMC8345831 | biostudies-literature
Project description:Reducing child undernutrition is a key social policy objective of the Ethiopian government. Despite substantial reduction over the last decade and a half, child undernutrition is still high; with 48 percent of children either stunted, underweight or wasted, undernutrition remains an important child health challenge. The existing literature highlights that targeting of efforts to reduce undernutrition in Ethiopia is inefficient, in part due to lack of data and updated information. This paper remedies some of this shortfall by estimating levels of stunting and underweight in each woreda for 2014. The estimates are small area estimations based on the 2014 Demographic and Health Survey and the latest population census. It is shown that small area estimations are powerful predictors of undernutrition, even compared to household characteristics, such as wealth and education, and hence a valuable targeting metric. The results show large variations in share of children undernourished within each region, more than between regions. The results also show that the locations with larger challenges depend on the chosen undernutrition statistic, as the share, number and concentration of undernourished children point to vastly different locations. There is also limited correlation between share of children underweight and stunted across woredas, indicating that different locations face different challenges.
| S-EPMC5391934 | biostudies-literature
Project description:OBJECTIVES:We developed a method to evaluate geographic and temporal variations in community-level obesity prevalence and used that method to identify communities in Massachusetts that should be considered high priority communities for obesity control. METHODS:We developed small-area estimation models to estimate community-level obesity prevalence among community-living adults 18 years or older. Individual-level data from the Behavioral Risk Factors Surveillance System from 1999 to 2005 were integrated with community-level data from the 2000 US Census. Small-area estimation models assessed the associations of obesity (body mass index >or= 30 kg/m(2)) with individual- and community-level characteristics. A classification system based on level and precision of obesity prevalence estimates was then used to identify high-priority communities. RESULTS:Estimates of the prevalence of community-level obesity ranged from 9% to 38% in 2005 and increased in all communities from 1999 to 2005. Fewer than 7% of communities met the Healthy People 2010 objective of prevalence rates below 15%. The highest prevalence rates occurred in communities characterized by lower income, less education, and more blue-collar workers. CONCLUSIONS:Similar to the rest of the nation, Massachusetts faces a great challenge in reaching the national obesity control objective. Targeting high-priority communities identified by small-area estimation may maximize use of limited resources.
| S-EPMC2642528 | biostudies-literature
Project description:There is an urgent need for childhood surveillance systems to design, implement, and evaluate interventions at the local level. We estimated obesity prevalence for individuals aged 5-17 years using a southcentral Wisconsin EHR data repository, Public Health Information Exchange (PHINEX, 2007-2012). The prevalence estimates were calculated by aggregating the estimated probability of each individual being obese, which was obtained via a generalized linear mixed model. We incorporated the random effects at the area level into our model. A weighted procedure was employed to account for missingness in EHR data. A non-parametric kernel smoothing method was used to obtain the prevalence estimates for locations with no or little data (<20 individuals) from the EHR. These estimates were compared to results from newly available obesity atlas (2015-2016) developed from various EHRs with greater statewide representation. The mean of the zip code level obesity prevalence estimates for males and females aged 5-17 years is 16.2% (SD 2.72%); 17.9% (SD 2.14%) for males and 14.4% (SD 2.00%) for females. The results were comparable to the Wisconsin Health Atlas (WHA) estimates, a much larger dataset of local community EHRs in Wisconsin. On average, prevalence estimates were 2.12% lower in this process than the WHA estimates, with lower estimation occurring more frequently for zip codes without data in PHINEX. Using this approach, we can obtain estimates for local areas that lack EHRs data. Generally, lower prevalence estimates were produced for those locations not represented in the PHINEX database when compared to WHA estimates. This underscores the need to ensure that the reference EHRs database can be made sufficiently similar to the geographic areas where synthetic estimates are being created.
| S-EPMC7895416 | biostudies-literature
Project description:IntroductionEstimation of accurate time-integrated activity coefficients (TIACs) and radiation absorbed doses (ADs) is desirable for treatment planning in peptide-receptor radionuclide therapy (PRRT). This study aimed to investigate the accuracy of a simplified dosimetry using a physiologically-based pharmacokinetic (PBPK) model, a nonlinear mixed effect (NLME) model, and single-time-point imaging to calculate the TIACs and ADs of 90Y-DOTATATE in various organs of dosimetric interest and tumors.Materials & methodsBiokinetic data of 111In-DOTATATE in tumors, kidneys, liver, spleen, and whole body were obtained from eight patients using planar scintigraphic imaging at T1 = (2.9 ± 0.6), T2 = (4.6 ± 0.4), T3 = (22.8 ± 1.6), T4 = (46.7 ± 1.7) and T5 = (70.9 ± 1.0) h post injection. Serum activity concentration was measured at 5 and 15 min; 0.5, 1, 2, and 4 h; and 1, 2, and 3 d p.i.. A published PBPK model for PRRT, NLME, and a single-time-point imaging datum at different time points were used to calculate TIACs in tumors, kidneys, liver, spleen, whole body, and serum. Relative deviations (RDs) (median [min, max]) between the calculated TIACs from single-time-point imaging were compared to the TIACs calculated from the all-time-points fit. The root mean square error (RMSE) of the difference between the computed ADs from the single-time-point imaging and reference ADs from the all-time point fittings were analyzed. A joint root mean square error RMSEjoint of the ADs was calculated with the RSME from both the tumor and kidneys to sort the time points concerning accurate results for the kidneys and tumor dosimetry. The calculations of TIACs and ADs from the single-time-point dosimetry were repeated using the sum of exponentials (SOE) approach introduced in the literature. The RDs and the RSME of the PBPK approach in our study were compared to the SOE approach.ResultsUsing the PBPK and NLME models and the biokinetic measurements resulted in a good fit based on visual inspection of the fitted curves and the coefficient of variation CV of the fitted parameters (<50%). T4 was identified being the time point with a relatively low median and range of TIACs RDs, i.e., 5 [1, 21]% and 2 [-15, 21]% for kidneys and tumors, respectively. T4 was found to be the time point with the lowest joint root mean square error RMSEjoint of the ADs. Based on the RD and RMSE, our results show a similar performance as the SOE and NLME model approach.SummaryIn this study, we introduced a simplified calculation of TIACs/ADs using a PBPK model, an NLME model, and a single-time-point measurement. Our results suggest a single measurement might be used to calculate TIACs/ADs in the kidneys and tumors during PRRT.
| S-EPMC10082376 | biostudies-literature
Project description:Genome-wide association studies (GWASs) have successfully identified thousands of genetic variants for many complex diseases; however, these variants explain only a small fraction of the heritability. Recently, genetic association studies that leverage external transcriptome data have received much attention and shown promise for discovering novel variants. One such approach, PrediXcan, is to use predicted gene expression through genetic regulation. However, there are limitations in this approach. The predicted gene expression may be biased, resulting from regularized regression applied to moderately sample-sized reference studies. Further, some variants can individually influence disease risk through alternative functional mechanisms besides expression. Thus, testing only the association of predicted gene expression as proposed in PrediXcan will potentially lose power. To tackle these challenges, we consider a unified mixed effects model that formulates the association of intermediate phenotypes such as imputed gene expression through fixed effects, while allowing residual effects of individual variants to be random. We consider a set-based score testing framework, MiST (mixed effects score test), and propose two data-driven combination approaches to jointly test for the fixed and random effects. We establish the asymptotic distributions, which enable rapid calculation of p values for genome-wide analyses, and provide p values for fixed and random effects separately to enhance interpretability over GWASs. Extensive simulations demonstrate that our approaches are more powerful than existing ones. We apply our approach to a large-scale GWAS of colorectal cancer and identify two genes, POU5F1B and ATF1, which would have otherwise been missed by PrediXcan, after adjusting for all known loci.
| S-EPMC5986723 | biostudies-literature
Project description:Accurately quantifying species' area requirements is a prerequisite for effective area-based conservation. This typically involves collecting tracking data on species of interest and then conducting home-range analyses. Problematically, autocorrelation in tracking data can result in space needs being severely underestimated. Based on the previous work, we hypothesized the magnitude of underestimation varies with body mass, a relationship that could have serious conservation implications. To evaluate this hypothesis for terrestrial mammals, we estimated home-range areas with global positioning system (GPS) locations from 757 individuals across 61 globally distributed mammalian species with body masses ranging from 0.4 to 4000 kg. We then applied block cross-validation to quantify bias in empirical home-range estimates. Area requirements of mammals <10 kg were underestimated by a mean approximately15%, and species weighing approximately100 kg were underestimated by approximately50% on average. Thus, we found area estimation was subject to autocorrelation-induced bias that was worse for large species. Combined with the fact that extinction risk increases as body mass increases, the allometric scaling of bias we observed suggests the most threatened species are also likely to be those with the least accurate home-range estimates. As a correction, we tested whether data thinning or autocorrelation-informed home-range estimation minimized the scaling effect of autocorrelation on area estimates. Data thinning required an approximately93% data loss to achieve statistical independence with 95% confidence and was, therefore, not a viable solution. In contrast, autocorrelation-informed home-range estimation resulted in consistently accurate estimates irrespective of mass. When relating body mass to home range size, we detected that correcting for autocorrelation resulted in a scaling exponent significantly >1, meaning the scaling of the relationship changed substantially at the upper end of the mass spectrum.
| S-EPMC7496598 | biostudies-literature