Evolution of mutation rates in hypermutable populations of Escherichia coli propagated at very small effective population size
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Travel burdens are a major barrier to healthcare access among US Veteran patient populations, particularly those residing in rural areas. Spatial accessibility to points of care for US Veteran populations is commonly assessed in two ways. The first approach uses open data from the US Census to represent collective travel burdens, for example the distance between population-weighted census tract centroids and VHA points of care. The second approach uses restricted-access VHA patient data to measure travel costs (e.g., distance, time) for accessing points of care with respect to geolocated patient addresses and real or approximated transportation networks. While the advantage of the open data approach lies in its reproducibility, it has notable limitations in its tendency to infer individual travel behavior from aggregate population characteristics, a problem known as ecological fallacy. Conversely, while the patient data approach is able to account for individual travel behavior, its ability to account for localized access disparities (e.g., a neighborhood with exceptionally high transportation costs) and patient demographics is limited as protecting individual patient data requires their storage in closed systems with limited capacity for adequately modeling real-world travel patterns or for supplementing patient attributes. Additionally, the patient data approach cannot account for veterans who are not enrolled in the VHA system but who may be eligible for care. These challenges limit the ability to perform “what if” analyses on the effects of place-specific interventions on veteran populations with high access barriers to healthcare. To address these challenges, we explore the application of realistic synthetic populations to examine travel burdens and spatial accessibility issues among veteran patient populations. Synthetic populations provide a virtual, individually-resolved and cross-sectional representation of the veteran patient population that enables investigation of spatial access to points of care in ways in which aggregate data and patient data do not. First, synthetic populations allow one to directly assess how individuals access points of care, from synthesized residential locations to outpatient facilities on real-world transportation networks. Modeling access to points of care at the individual scale addresses the ecological fallacy problem associated with using aggregated census data to represent veteran populations and patterns of movement. Second, synthetic populations provide a means of completely representing an area’s veteran population using only publicly available, anonymized census microdata from the American Community Survey (ACS) to ensure the privacy of real-world individuals. Generating synthetic populations from the ACS also expands descriptive characteristics beyond what patient data typically offers to include socio-demographic, economic, housing, and mobility attributes. More detailed profiles of both VHA patient populations and veterans not enrolled in the VA system will provide a comprehensive picture of groups that may benefit from interventions or outreach. As an initial exercise for using synthetic populations to measure veteran travel burdens to VA care, we apply Oak Ridge National Laboratory’s (ORNL) UrbanPop capability to generate a series of synthetic VHA patient populations for 9 Veterans Integrated Services Networks (VISN) market areas in 9 Census Divisions across the continental United States, which are listed in Table 1. We use UrbanPop to produce synthetic populations for the VISN markets selected for each US Census Division, then assign VA outpatient clinic destinations to synthetic VHA patients based on travel about each VISN market’s road network. To demonstrate using the synthetic populations to evaluate healthcare travel burdens, we compare the time-based impedance between simulated home locations and VA outpatient clinics in each VISN market. We then perform validation exercises on the synthetic populations with respect to neighborhood (block group) demographic composition as well as patient mobility, comparing aggregate origin-destination statistics for the synthetic population to outpatient visits available in restricted patient data from the VA’s Corporate Data Warehouse (CDW) database.
Abstract Grouping populations for management may overlook the fine-scale diversity underpinning the stability and resilience of meta-populations and fisheries. A bimodal timing distribution of summer-run steelhead Oncorhynchus mykiss (anadromous Rainbow Trout) historically was observed at Bonneville Dam (BON), the first barrier to upstream migration in the Columbia River basin. Early mode fish (A-run) tended to be younger and smaller (<78 cm) than later fish (B-run). While A-run fish spawn throughout the Columbia River basin, B-run fish spawn primarily in the Snake River basin. Managers used indices of these modes to make fishery decisions, and later these criteria were adopted for conservation. It is still unclear how life history and body size differences among wild Snake River populations are related to the categories at BON. We examined population parameters characterizing the two categories (date of passage at BON, length) and parameters directly affecting population dynamics (age composition, sex ratio). The life history portfolio of Snake River steelhead is quite diverse. There was broad overlap among populations in several respects, forming a gradient in life history characteristics rather than a dichotomous break. All populations produced adults <78 cm and adults returning after August 25. Median lengths of putative B-run populations were close to the criterion that was supposed to be a defining characteristic. In contrast, few A-run populations produced many adults ≥78 cm. Mean percentage of two-ocean fish was 52.1% for A-run populations and 82.0% for B-run populations. Mean age at spawn was greater in populations producing older smolts. Sex ratio was female biased, and older populations had greater percentages of females. Although the run-type dichotomy was useful for management of fisheries in the past, it is not useful for conservation. A combination of genetic stock identification at main-stem dams and population-specific monitoring in natal streams provides a unified framework for the assessment of fisheries management and conservation objectives. Received March 2, 2016; accepted November 9, 2016 Published online March 10, 2017
Genomic prediction is a useful tool to accelerate genetic gain in selection using DNA marker information. However, this technology typically relies on standard prediction procedures, such as genomic BLUP, that are not designed to accommodate population heterogeneity resulting from differences in marker effects across populations. In this study, we assayed different prediction procedures to capture marker-by-population interactions in genomic prediction models. Prediction procedures included genomic BLUP and two kernel-based extensions of genomic BLUP which explicitly accounted for population heterogeneity. To model population heterogeneity, dissemblance between populations was either depicted by a unique coefficient (as previously reported), or a more flexible function of genetic distance between populations (proposed herein). Models under investigation were applied in a diverse switchgrass sample under two validation schemes: whole-sample calibration, where all individuals except selection candidates are included in the calibration set, and cross-population calibration, where the target population is entirely excluded from the calibration set. First, we showed that using fixed effects, from principal components or putative population groups, appeared detrimental to prediction accuracy, especially in cross-population calibration. Then we showed that modelling population heterogeneity by our proposed procedure resulted in highly significant improvements in model fit. In such cases, gains in accuracy were often positive. These results suggest that population heterogeneity may be parsimoniously captured by kernel methods. However, in cases where improvement in model fit by our proposed procedure is null-to-moderate, ignoring heterogeneity should probably be preferred due to the robustness and simplicity of the standard genomic BLUP model.
Abstract Trapping is an effective strategy for control of invasive wild pig (Sus scrofa), but eradication of populations via trapping alone typically is impractical throughout much of the species’ introduced range. Consequently, recovery of wild pig populations to pre-trapping levels following control is a major impediment to success of long-term control. In this study, we conducted camera-trap surveys to monitor recovery of two wild pig populations immediately following density reductions achieved via trapping on two sites on the Savannah River Site, SC. We used spatial capture-recapture models to estimate wild pig density during seven 3 week periods spanning 12 months to: (1) determine the time required for populations to recover to pre-trapping levels; and (2) evaluate trends in population growth, density, and distribution during recovery. Recovery of both wild pig populations started immediately after trapping ceased, with populations reaching pre-trapping levels in 3 and 5 months. Threeweek population growth rates for each wild pig population during recovery averaged 1.28 (range = 1.11–1.39 on site 1) and 1.24 (range = 0.64–1.92 on site 2), which corresponded to a 345% and 132% increase in abundance (an increase of 141 and 126 individual wild pigs on site 1 and 2, respectively) of each population by the last survey period. Rapid recovery and growth beyond pre-trapping levels of wild pig populations likely were driven by residual wild pigs not removed by trapping and compensatory reproduction and immigration. Our study is one of the first to provide information on how quickly wild pig populations can recover following overall density reductions of 54–68% and indicates practitioners should be prepared to reinitiate trapping efforts in as little as 3 months to prevent recovery of wild pig populations following control.
The building of the NASA Orbital Debris Engineering Model, whether ORDEM2000 or its recently updated version ORDEM2010, uses as its foundation a number of model debris populations, each truncated at a minimum object-size ranging from 10 micron to 1 m. This paper discusses the development of the ORDEM2010 model debris populations in LEO (low Earth orbit), focusing on centimeter (smaller than 10 cm) and millimeter size regimes. Primary data sets used in the statistical derivation of the cm- and mm-size model populations are from the Haystack radar operated in a staring mode. Unlike cataloged objects of sizes greater than approximately 10 cm, ground-based radars monitor smaller-size debris only in a statistical manner instead of tracking every piece. The mono-static Haystack radar can detect debris as small as approximately 5 mm at moderate LEO altitudes. Estimation of millimeter debris populations (for objects smaller than approximately 6 mm) rests largely on Goldstone radar measurements. The bi-static Goldstone radar can detect 2- to 3-mm objects. The modeling of the cm- and mm-debris populations follows the general approach to developing other ORDEM2010-required model populations for various components and types of debris. It relies on appropriate reference populations to provide necessary prior information on the orbital structures and other important characteristics of the debris objects. NASA's LEO-to-GEO Environment Debris (LEGEND) model is capable of furnishing such reference populations in the desired size range. A Bayesian statistical inference process, commonly adopted in ORDEM2010 model-population derivations, changes a priori distribution into a posteriori distribution and thus refines the reference populations in terms of data. This paper describes key elements and major steps in the statistical derivations of the cm- and mm-size debris populations and presents results. Due to lack of data for near 1-mm sizes, the model populations of 1- to 3.16-mm objects are an empirical extension from larger debris. The extension takes into account the results of micro-debris (from 10 micron to 1 mm) population modeling that is based on shuttle impact data, in the hope of making a smooth transition between micron and millimeter size regimes. This paper also includes a brief discussion on issues and potential future work concerning the analysis and interpretation of Goldstone radar data.
This chapter discusses the basis for these population estimates, scope, and limitations based on experiences in development of massive global population datasets and usage of these datasets as a basis for sampling design. It presents tools and approaches for using these georeferenced population estimates for complex household survey sampling. The chapter provides a unified resource for understanding current gridded population datasets, their use in survey research, and promising areas of future work toward improved population estimates. It includes an overview of popular gridded population datasets, common methodological frameworks for developing gridded population estimates, and discusses pros and cons for selecting a gridded population dataset from a survey perspective. The chapter focuses on survey methods that use gridded population estimates specifically. It highlights differences between census population estimates and gridded population estimates that may be relevant when planning data collection. The chapter presents a case study of gridded population data and sampling methods in Nigeria.
ABSTRACT We look at the contribution of different galaxy populations to the atomic hydrogen (H i) mass function (HIMF) and the H i density parameter, ΩH i, in the local Universe. Our analysis is based on a sample of 7857 H i selected galaxies selected from a volume common to the Sloan Digital Sky Survey (SDSS) and Arecibo Fast Legacy ALFA (40 per cent catalogue – α.40) surveys. We define different populations of galaxies in the colour (u − r)–magnitude (Mr) plane and compute the HIMF for each of them. Additionally we compute the HIMF for dark galaxies; these are undetected in SDSS and represent $\sim 2{{\ \rm per\ cent}}$ of the total sample. We find that the luminous red population dominates the total HIMF for $\log _{10}(M_{\rm{H\,{\small I}}}h^2_{70}/M_{\odot }) \ge 10.4$. The full red population – luminous and faint – represents about $\sim 17{{\ \rm per\ cent}}$ of the ΩH i budget, while that of the dark population is $\sim 3{{\ \rm per\ cent}}$. The HIMF about the knee, $\log _{10}(M_{\rm{H\,{\small I}}}h^2_{70}/M_{\odot }) \in [8,10.4]$, is dominated by the faint and luminous blue populations, the latter dominating at larger masses in this interval. Their total contribution to ΩH i is $\sim 55\!-\!70{{\ \rm per\ cent}}$, the variation depending on the definition of population. The dominant populations at the low-mass end, $\log _{10}(M_{\rm{H\,{\small I}}}h^2_{70}/M_{\odot }) \le 8.0$ are the faint blue and faint bluer populations, the latter’s dominance being sensitive to its definition. The full blue (blue–bluer luminous and faint) population represents $\sim 80{{\ \rm per\ cent}}$ of ΩH i. A bimodal HIMF suggested by our results is however not seen since the amplitude of the HIMF of the luminous red population is small compared to that of the luminous blue population.
Population data are normally collected at various census administrative levels, and areal interpolation of population is often required to transform population to the desired spatial resolution. Building footprint datasets, such as Microsoft building footprints, have proven to be useful in estimating population distribution and can therefore be used for areal interpolation of population. In addition to Microsoft building footprints, the recently released USA Structures dataset provides additional information such as building type and building height for some regions, which may provide valuable information for a better depiction of population distribution and improved population areal interpolation accuracy. In this study, we have conducted areal interpolation of population projections consistent with three different Shared Socioeconomic Pathways (SSP2, SSP3, and SSP5) from 1-km grid cells to block level in Washington state for every ten years from 2020 to 2040 based on USA Structures. We assessed USA Structures-based population downscaling accuracy using U.S. decennial survey data in 2020 under three different downscaling schemes, including population downscaling from census tracts to block groups, from census tracts to blocks, and from block groups to blocks. The resulting accuracies were compared with those based on Microsoft building footprints. The comparison showed that USA Structures achieved higher accuracies across different population density regions and areas with different urbanization extent within our study area.
Datasets supporting findings and visualizations behind McManamay et al. (2024) Divergent urban land trajectories under alternative population projections within the shared socioeconomic pathways. Environmental Research Letters, DOI: 10.1088/1748-9326/ad2eec Using a spatial modeling experiment at high resolution (1-km), this study compared how two alternative US population projections, varying in the spatially explicit nature of demographic patterns and migration, affect urban land dynamics simulated by the Spatially Explicit, Long-term, Empirical City development (SELECT) model for SSP2, SSP3, and SSP5. The numerical experiment, inputs and outputs are fully described in McManamay et al. (2024). The datasets summarize SELECT model simulations within urban areas and rural areas of the conterminous United States. For code to reproduce the results, please refer to https://github.com/IMMM-SFA/mcmanamay_etal_2024_erl Please refer to the README file provided in Files for more details. Descriptions of the datasets are provided below. Dataset(s) Descriptions: Urban_Land_delta_UA_County.csv: Urban land area for original (default) population, updated population, and urban land delta (difference between the two) according to Shared Socioeconomic Pathways (SSPs) and year. Urban land is summarized for U.S. counties (FIPS code) and urban areas (GEOID) as unique areas within various counties. JO_Gao_pop_ufdelta.csv: Comparison of population differences and urban land differences between the origial (default) and updated population projections. City_Case_studies.csv: Differences in urban land areas arising from different population projections summarized for selected cities (Atlanta, Los Angeles, Houston, and New York) and their surrounding rural areas UA_county_FID.zip: Zipped folder of .tiff file depicting unique combination of urban areas and US counties used for summarizing urban land differences in urban and rural areas based on different population projections. ULD_clusters.zip: Zipped folder of .csv files depicting normalized changes in urban land delta (difference in urban land arising from population projections) and the respective Ward's hierarchical clusters, which group urban and rural areas based on similarities in temporal trends. Urban_Percent_change_2100.zip: Zipped folder of .csv files depicting percent changes in urban land area for year 2100 based on differences in the original (default) population and updated population. Each file is for SSPs and urban or rural areas.
There is overwhelming evidence that forest trees are locally adapted to climate. Thus, genecological models based on population phenotypes have been used to measure local adaptation, infer genetic maladaptation to climate, and guide assisted migration. However, instead of phenotypes, there is increasing interest in using genomic data for gene resource management. We used whole-genome resequencing and common-garden experiments to understand the genetic architecture of adaptive traits in black cottonwood. We studied the potential of using genome-wide association studies (GWAS) and genomic prediction to detect causal loci, identify climate-adapted phenotypes, and inform gene resource management. We analyzed population structure by partitioning phenotypic and genomic (single-nucleotide polymorphism) variation among 840 genotypes collected from 91 stands along 16 rivers. Most phenotypic variation (60 to 81%) occurred among populations and was strongly associated with climate. Population phenotypes were predicted well using genomic data (e.g., predictive abilityr> 0.9) but almost as well using climate or geography (r> 0.8). In contrast, genomic prediction within populations was poor (r< 0.2). We identified many GWAS associations among populations, but most appeared to be spurious based on pooled within-population analyses. Hierarchical partitioning of linkage disequilibrium and haplotype sharing suggested that within-population genomic prediction and GWAS were poor because allele frequencies of causal loci and linked markers differed among populations. Given the urgent need to conserve natural populations and ecosystems, our results suggest that climate variables alone can be used to predict population phenotypes, delineate seed zones and deployment zones, and guide assisted migration.
Population change is a main driver behind global environmental change, including urban land expansion. In future scenario modeling, assumptions regarding how populations will change locally, despite identical global constraints of Shared Socioeconomic Pathways (SSPs), can have dramatic effects on subsequent regional urbanization. Using a spatial modeling experiment at high resolution (1 km), this study compared how two alternative US population projections, varying in the spatially explicit nature of demographic patterns and migration, affect urban land dynamics simulated by the Spatially Explicit, Long-term, Empirical City development (SELECT) model for SSP2, SSP3, and SSP5. The population projections included: (1) newer downscaled state-specific population (SP) projections inclusive of updated international and domestic migration estimates, and (2) prevailing downscaled national-level projections (NP) agnostic to localized demographic processes. Our work shows that alternative population inputs, even those under the same SSP, can lead to dramatic and complex differences in urban land outcomes. Under the SP projection, urbanization displays more of an extensification pattern compared to the NP projection. This suggests that recent demographic information supports more extreme urban extensification and land pressures on existing rural areas in the US than previously anticipated. Urban land outcomes to population inputs were spatially variable where areas in close spatial proximity showed divergent patterns, reflective of the spatially complex urbanization processes that can be accommodated in SELECT. Although different population projections and assumptions led to divergent outcomes, urban land development is not a linear product of population change but the result of complex relationships between population, dynamic urbanization processes, stages of urban development maturity, and feedback mechanisms. These findings highlight the importance of accounting for spatial variations in the population projections, but also urbanization process to accurately project long-term urban land patterns.
With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.
Ecological risk assessment (ERA) is charged with assessing the likelihood a chemical will have adverse environmental or ecological effects. When assessing the risk of a potential contaminant to biological organisms, ecologists are most concerned with the sustainability of populations of organisms, rather than protecting every individual. However, ERA most commonly relies on data on the effect of a potential contaminant on individuals because these experiments are more feasible than costly population-level exposures. In this work, we address the challenge of extrapolating these individual-level results to predict population-level effects. Previous per-capita population growth rate estimates calculated from individual-level exposures of Daphnia pulicaria to silver nanoparticles (AgNPs) at different food rations predict a critical daily food requirement for daphnid populations exposed to 200 μg/L AgNPs to avoid extinction. To test this, we exposed daphnid populations to the same AgNP concentration at three different food inputs, with the lowest ration close to the extinction threshold predicted from data on individuals. The two populations with the higher food inputs persisted, and the population with the lowest food input went extinct after 50 days but did persist through two generations. Further, we demonstrate that we can extrapolate between these levels of biological organization by parameterizing an individual-level biomass model with data on individuals’ response to AgNPs and using these parameters to predict the outcome for control and AgNP-exposed populations. Key to successful extrapolation is careful modeling of temporal changes in resource density, driven by both the experimental protocols and feedback from the consumer. The implication for ecotoxicology is that estimates of extinction thresholds based on studies of individuals may be reliable predictors of population outcomes, but only with careful treatment of resource dynamics.
To understand the forces driving differentiation and diversification in wild bacterial populations, we must be able to delineate and track ecologically relevant units through space and time. Mapping metagenomic sequences to reference genomes derived from the same environment can reveal genetic heterogeneity within populations, and in some cases, be used to identify boundaries between genetically similar, but ecologically distinct, populations. Here we examine population-level heterogeneity within abundant and ubiquitous freshwater bacterial groups such as the acI Actinobacteria and LD12 Alphaproteobacteria (the freshwater sister clade to the marine SAR11) using 33 single-cell genomes and a 5-year metagenomic time series. The single-cell genomes grouped into 15 monophyletic clusters (termed "tribes") that share at least 97.9% 16S rRNA identity. Distinct populations were identified within most tribes based on the patterns of metagenomic read recruitments to single-cell genomes representing these tribes. Genetically distinct populations within tribes of the acI Actinobacterial lineage living in the same lake had different seasonal abundance patterns, suggesting these populations were also ecologically distinct. In contrast, sympatric LD12 populations were less genetically differentiated. This suggests that within one lake, some freshwater lineages harbor genetically discrete (but still closely related) and ecologically distinct populations, while other lineages are composed of less differentiated populations with overlapping niches. Our results point at an interplay of evolutionary and ecological forces acting on these communities that can be observed in real time.
Gridded population datasets represent high-resolution distributions of human occupancy, enabling informed decision-making across a broad range of fields. These data products are valuable for assessing environmental risk, urban development, disaster preparedness and resource allocation—areas where accurate population estimates directly enhance policy effectiveness and optimize resource distribution. Despite the importance of gridded population datasets, traditional population modeling approaches often overlook inherent uncertainties in the estimation process. This limitation can create a false sense of certainty in population estimates, potentially leading to flawed decisions by those who rely on the data. To address this methodological gap, we introduce a probabilistic machine learning modeling framework, LandScan Mosaic, that explicitly incorporates uncertainty into the population modeling process. Our approach systematically quantifies uncertainty in three key modeling parameters of the LandScan HD gridded population dataset: building use types, floor counts, and occupancy rates. By employing Monte Carlo simulations, we propagate these uncertainties through the modeling process, yielding probability distributions of population counts in place of deterministic point estimates. We demonstrate the practical application of this framework in Iloilo City, Philippines, using structured decision-making techniques and our probabilistic estimates to identify and prioritize areas most affected by projected flooding, supporting targeted interventions that address both economic and social risks. In doing so, we propose a population-specific approach for incorporating confidence into structured decision making processes. Through a comparative analysis with conventional deterministic approaches and point estimate approaches, including LandScan HD and WorldPop, we evaluate how the incorporation of machine learning and uncertainty influences decision rankings. This research advances population distribution modeling by offering a robust, quantitative approach that explicitly accounts for uncertainty in the underlying data, along with guidance for how users can apply uncertainty in their decision-making.
Population estimates are critical for government services, development projects, and public health campaigns. Such data are typically obtained through a national population and housing census. However, population estimates can quickly become inaccurate in localized areas, particularly where migration or displacement has occurred. Some conflict-affected and resource-poor countries have not conducted a census in over 10 y. We developed a hierarchical Bayesian model to estimate population numbers in small areas based on enumeration data from sample areas and nationwide information about administrative boundaries, building locations, settlement types, and other factors related to population density. We demonstrated this model by estimating population sizes in every 10- m grid cell in Nigeria with national coverage. These gridded population estimates and areal population totals derived from them are accompanied by estimates of uncertainty based on Bayesian posterior probabilities. The model had an overall error rate of 67 people per hectare (mean of absolute residuals) or 43% (using scaled residuals) for predictions in out-of-sample survey areas (approximately 3 ha each), with increased precision expected for aggregated population totals in larger areas. This statistical approach represents a significant step toward estimating populations at high resolution with national coverage in the absence of a complete and recent census, while also providing reliable estimates of uncertainty to support informed decision making.
The purpose of this study was to assess the validity of a practical diabetes risk score amongst two heterogenous populations, a working population and a non-working population. Study population 1 (n = 2,089) participated in a large-scale screening program offered to retired workers to discover previously undetected/incipient chronic illness. Study population 2 (n = 3,293) was part of a Colorado worksite wellness program health risk assessment. We assessed the relationship between a continuous diabetes risk score at baseline and development of diabetes in the future using logistic regression. Receiver operating curves and sensitivity/specificity of the models were calculated. Across both study populations, we observed that participants with diabetes at follow-up had higher diabetes risk scores at baseline than participants who did not have diabetes at follow-up. On average, the odds ratio of developing diabetes in the future was 1.38 (95% CI: 1.26–1.50, p < 0.0001) for study population 1 and 1.68 (95% CI: 1.45–1.95, p-value < 0.0001) for study population 2. These findings indicate that the diabetes risk score may be generalizable to diverse individuals, and thus potentially a population level diabetes screening tool. Minimally-invasive diabetes risk scores can aid in the identification of sub-populations of individuals at risk for diabetes.