1
|
Keller EP, Lawson AB, Wagner CL, Reed SG. Bayesian modeling of spatially differentiated multivariate enamel defects of the children's primary maxillary central incisor teeth. BMC Med Res Methodol 2024; 24:88. [PMID: 38622506 PMCID: PMC11017560 DOI: 10.1186/s12874-024-02211-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2023] [Accepted: 04/02/2024] [Indexed: 04/17/2024] Open
Abstract
BACKGROUND The analysis of dental caries has been a major focus of recent work on modeling dental defect data. While a dental caries focus is of major importance in dental research, the examination of developmental defects which could also contribute at an early stage of dental caries formation, is also of potential interest. This paper proposes a set of methods which address the appearance of different combinations of defects across different tooth regions. In our modeling we assess the linkages between tooth region development and both the type of defect and associations with etiological predictors of the defects which could be influential at different times during the tooth crown development. METHODS We develop different hierarchical model formulations under the Bayesian paradigm to assess exposures during primary central incisor (PMCI) tooth development and PMCI defects. We evaluate the Bayesian hierarchical models under various simulation scenarios to compare their performance with both simulated dental defect data and real data from a motivating application. RESULTS The proposed model provides inference on identifying a subset of etiological predictors of an individual defect accounting for the correlation between tooth regions and on identifying a subset of etiological predictors for the joint effect of defects. Furthermore, the model provides inference on the correlation between the regions of the teeth as well as between the joint effect of the developmental enamel defects and dental caries. Simulation results show that the proposed model consistently yields steady inferences in identifying etiological biomarkers associated with the outcome of localized developmental enamel defects and dental caries under varying simulation scenarios as deemed by small mean square error (MSE) when comparing the simulation results to real application results. CONCLUSION We evaluate the proposed model under varying simulation scenarios to develop a model for multivariate dental defects and dental caries assuming a flexible covariance structure that can handle regional and joint effects. The proposed model shed new light on methods for capturing inclusive predictors in different multivariate joint models under the same covariance structure and provides a natural extension to a nested hierarchical model.
Collapse
|
2
|
Rotejanaprasert C, Chinpong K, Lawson AB, Chienwichai P, Maude RJ. Evaluation and comparison of spatial cluster detection methods for improved decision making of disease surveillance: a case study of national dengue surveillance in Thailand. BMC Med Res Methodol 2024; 24:14. [PMID: 38243198 PMCID: PMC10797994 DOI: 10.1186/s12874-023-02135-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2023] [Accepted: 12/21/2023] [Indexed: 01/21/2024] Open
Abstract
BACKGROUND Dengue is a mosquito-borne disease that causes over 300 million infections worldwide each year with no specific treatment available. Effective surveillance systems are needed for outbreak detection and resource allocation. Spatial cluster detection methods are commonly used, but no general guidance exists on the most appropriate method for dengue surveillance. Therefore, a comprehensive study is needed to assess different methods and provide guidance for dengue surveillance programs. METHODS To evaluate the effectiveness of different cluster detection methods for dengue surveillance, we selected and assessed commonly used methods: Getis Ord [Formula: see text], Local Moran, SaTScan, and Bayesian modeling. We conducted a simulation study to compare their performance in detecting clusters, and applied all methods to a case study of dengue surveillance in Thailand in 2019 to further evaluate their practical utility. RESULTS In the simulation study, Getis Ord [Formula: see text] and Local Moran had similar performance, with most misdetections occurring at cluster boundaries and isolated hotspots. SaTScan showed better precision but was less effective at detecting inner outliers, although it performed well on large outbreaks. Bayesian convolution modeling had the highest overall precision in the simulation study. In the dengue case study in Thailand, Getis Ord [Formula: see text] and Local Moran missed most disease clusters, while SaTScan was mostly able to detect a large cluster. Bayesian disease mapping seemed to be the most effective, with adaptive detection of irregularly shaped disease anomalies. CONCLUSIONS Bayesian modeling showed to be the most effective method, demonstrating the best accuracy in adaptively identifying irregularly shaped disease anomalies. In contrast, SaTScan excelled in detecting large outbreaks and regular forms. This study provides empirical evidence for the selection of appropriate tools for dengue surveillance in Thailand, with potential applicability to other disease control programs in similar settings.
Collapse
|
3
|
Reed SG, Fan S, Wagner CL, Lawson AB. Predictors of Developmental Defects of Enamel in Primary Maxillary Central Incisors Using Bayesian Model Selection. Caries Res 2023; 58:30-38. [PMID: 37918363 PMCID: PMC10922907 DOI: 10.1159/000534793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2022] [Accepted: 10/22/2023] [Indexed: 11/04/2023] Open
Abstract
INTRODUCTION Localized non-inheritable developmental defects of tooth enamel (DDE) are classified as enamel hypoplasia (EH), opacity (OP), and post-eruptive breakdown (PEB) using the enamel defects index. To better understand the etiology of DDE, we assessed the linkages amongst exposome variables for these defects during the specific time duration for enamel mineralization of the human primary maxillary central incisor enamel crowns. In general, these two teeth develop between 13 and 14 weeks in utero and 3-4 weeks' postpartum of a full-term delivery, followed by tooth eruption at about 1 year of age. METHODS We utilized existing datasets for mother-child dyads that encompassed 12 weeks' gestation through birth and early infancy, and child DDE outcomes from digital images of the erupted primary maxillary central incisor teeth. We applied a Bayesian modeling paradigm to assess the important predictors of EH, OP, and PEB. RESULTS The results of Gibbs variable selection showed a key set of predictors: mother's prepregnancy body mass index (BMI); maternal serum concentrations of calcium and phosphorus at gestational week 28; child's gestational age; and both mother's and child's functional vitamin D deficiency (FVDD). In this sample of healthy mothers and children, significant predictors for OP included the child having a gestational period >36 weeks and FVDD at birth, and for PEB included a mother's prepregnancy BMI <21.5 and higher serum phosphorus concentration at week 28. CONCLUSION In conclusion, our methodology and results provide a roadmap for assessing timely biomarker measures of exposures during specific tooth development to better understand the etiology of DDE for future prevention.
Collapse
|
4
|
Lawson AB, Kim J, Johnson C, Ratnapradipa KL, Alberg AJ, Akonde M, Hastert T, Bandera EV, Terry P, Mandle H, Cote ML, Bondy M, Marks J, Peres LC, Schildkraut J, Peters ES. The Association between Mediated Deprivation and Ovarian Cancer Survival among African American Women. Cancers (Basel) 2023; 15:4848. [PMID: 37835542 PMCID: PMC10571563 DOI: 10.3390/cancers15194848] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 09/25/2023] [Accepted: 10/02/2023] [Indexed: 10/15/2023] Open
Abstract
BACKGROUND Deprivation indices are often used to adjust for socio-economic disparities in health studies. Their role has been partially evaluated for certain population-level cancer outcomes, but examination of their role in ovarian cancer is limited. In this study, we evaluated a range of well-recognized deprivation indices in relation to cancer survival in a cohort of self-identified Black women diagnosed with ovarian cancer. This study aimed to determine if clinical or diagnostic characteristics lie on a mediating pathway between socioeconomic status (SES) and deprivation and ovarian cancer survival in a minority population that experiences worse survival from ovarian cancer. METHODS We used mediation analysis to look at the direct and indirect causal effects of deprivation indices with main mediators of the SEER stage at diagnosis and residual disease. The analysis employed Bayesian structural equation models with variable selection. We applied a joint Bayesian structural model for the mediator, including a Weibull mixed model for the vital outcome with deprivation as exposure. We selected modifiers via a Monte Carlo model selection procedure. RESULTS The results suggest that high SES-related indices, such as Yost, Kolak urbanicity (URB), mobility (MOB) and SES dimensions, and concentrated disadvantage index (CDI), all have a significant impact on improved survival. In contrast, area deprivation index (ADI)/Singh, and area level poverty (POV) did not have a major impact. In some cases, the indirect effects have very wide credible intervals, so the total effect is not well estimated despite the estimation of the direct effect. CONCLUSIONS First, it is clear that commonly used indices such as Yost, or CDI both significantly impact the survival experience of Black women diagnosed with epithelial ovarian cancer. In addition, the Kolak dimension indices (URB, MOB, mixed immigrant: MICA and SES) also demonstrate a significant association, depending on the mediator. Mediation effects differ according to the mediator chosen.
Collapse
|
5
|
Lawson AB, Kim J, Johnson C, Hastert T, Bandera EV, Alberg AJ, Terry P, Akonde M, Mandle H, Cote ML, Bondy M, Marks J, Peres L, Ratnapradipa KL, Xin Y, Schildkraut J, Peters ES. Deprivation and segregation in ovarian cancer survival among African American women: a mediation analysis. Ann Epidemiol 2023; 86:57-64. [PMID: 37423270 PMCID: PMC10538403 DOI: 10.1016/j.annepidem.2023.07.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Revised: 07/03/2023] [Accepted: 07/05/2023] [Indexed: 07/11/2023]
Abstract
PURPOSE Deprivation and segregation indices are often examined as possible explanations for observed health disparities in population-based studies. In this study, we assessed the role of recognized deprivation and segregation indices specifically as they affect survival in a cohort of self-identified Black women diagnosed with ovarian cancer who enrolled in the African American Cancer Epidemiology Study. METHODS Mediation analysis was used to examine the direct and indirect effects between deprivation or segregation and overall survival via a Bayesian structural equation model with Gibbs variable selection. RESULTS The results suggest that high socioeconomic status-related indices have an association with increased survival, ranging from 25% to 56%. In contrast, index of concentration at the extremes-race does not have a significant impact on overall survival. In many cases, the indirect effects have very wide credible intervals; consequently, the total effect is not well estimated despite the estimation of the direct effect. CONCLUSIONS Our results show that Black women living in higher socioeconomic status neighborhoods are associated with increased survival with ovarian cancer using area-level economic indices such as Yost or index of concentration at the extremes-income. In addition, the Kolak urbanization index has a similar impact and highlights the importance of area-level deprivation and segregation as potentially modifiable social factors in ovarian cancer survival.
Collapse
|
6
|
Johnson CE, Alberg AJ, Bandera EV, Peres LC, Akonde M, Collin LJ, Cote ML, Hastert TA, Hébert JR, Peters ES, Qin B, Terry P, Schwartz AG, Bondy M, Epstein MP, Mandle HB, Marks JR, Lawson AB, Schildkraut JM. Association of inflammation-related exposures and ovarian cancer survival in a multi-site cohort study of Black women. Br J Cancer 2023; 129:1119-1125. [PMID: 37537254 PMCID: PMC10539498 DOI: 10.1038/s41416-023-02385-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2023] [Revised: 07/25/2023] [Accepted: 07/26/2023] [Indexed: 08/05/2023] Open
Abstract
BACKGROUND An association was observed between an inflammation-related risk score (IRRS) and worse overall survival (OS) among a cohort of mostly White women with invasive epithelial ovarian cancer (EOC). Herein, we evaluated the association between the IRRS and OS among Black women with EOC, a population with higher frequencies of pro-inflammatory exposures and worse survival. METHODS The analysis included 592 Black women diagnosed with EOC from the African American Cancer Epidemiology Study (AACES). Cox proportional hazards models were used to compute hazard ratios (HRs) and 95% confidence intervals (CIs) for the association of the IRRS and OS, adjusting for relevant covariates. Additional inflammation-related exposures, including the energy-adjusted Dietary Inflammatory Index (E-DIITM), were evaluated. RESULTS A dose-response trend was observed showing higher IRRS was associated with worse OS (per quartile HR: 1.11, 95% CI: 1.01-1.22). Adding the E-DII to the model attenuated the association of IRRS with OS, and increasing E-DII, indicating a more pro-inflammatory diet, was associated with shorter OS (per quartile HR: 1.12, 95% CI: 1.02-1.24). Scoring high on both indices was associated with shorter OS (HR: 1.54, 95% CI: 1.16-2.06). CONCLUSION Higher levels of inflammation-related exposures were associated with decreased EOC OS among Black women.
Collapse
|
7
|
Lawson AB. Evaluation of predictive capability of Bayesian spatio-temporal models for Covid-19 spread. BMC Med Res Methodol 2023; 23:182. [PMID: 37568119 PMCID: PMC10422743 DOI: 10.1186/s12874-023-01997-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Accepted: 07/20/2023] [Indexed: 08/13/2023] Open
Abstract
BACKGROUND Bayesian models have been applied throughout the Covid-19 pandemic especially to model time series of case counts or deaths. Fewer examples exist of spatio-temporal modeling, even though the spatial spread of disease is a crucial factor in public health monitoring. The predictive capabilities of infectious disease models is also important. METHODS In this study, the ability of Bayesian hierarchical models to recover different parts of the variation in disease counts is the focus. It is clear that different measures provide different views of behavior when models are fitted prospectively. Over a series of time horizons one step predictions have been generated and compared for different models (for case counts and death counts). These Bayesian SIR models were fitted using MCMC at 28 time horizons to mimic prospective prediction. A range of goodness of prediction measures were analyzed across the different time horizons. RESULTS A particularly important result is that the peak intensity of case load is often under-estimated, while random spikes in case load can be mimicked using time dependent random effects. It is also clear that during the early wave of the pandemic simpler model forms are favored, but subsequently lagged spatial dependence models for cases are favored, even if the sophisticated models perform better overall. DISCUSSION The models fitted mimic the situation where at a given time the history of the process is known but the future must be predicted based on the current evolution which has been observed. Using an overall 'best' model for prediction based on retrospective fitting of the complete pandemic waves is an assumption. However it is also clear that this case count model is well favored over other forms. During the first wave a simpler time series model predicts case counts better for counties than a spatially dependent one. The picture is more varied for morality. CONCLUSIONS From a predictive point of view it is clear that spatio-temporal models applied to county level Covid-19 data within the US vary in how well they fit over time and also how well they predict future events. At different times, SIR case count models and also mortality models with cumulative counts perform better in terms of prediction. A fundamental result is that predictive capability of models varies over time and using the same model could lead to poor predictive performance. In addition it is clear that models addressing the spatial context for case counts (i.e. with lagged neighborhood terms) and cumulative case counts for mortality data are clearly better at modeling spatio-temporal data which is commonly available for the Covid-19 pandemic in different areas of the globe.
Collapse
|
8
|
Kim J, Lawson AB, Neelon B, Korte JE, Eberth JM, Chowell G. Evaluation of Bayesian spatiotemporal infectious disease models for prospective surveillance analysis. BMC Med Res Methodol 2023; 23:171. [PMID: 37481553 PMCID: PMC10363300 DOI: 10.1186/s12874-023-01987-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Accepted: 07/11/2023] [Indexed: 07/24/2023] Open
Abstract
BACKGROUND COVID-19 brought enormous challenges to public health surveillance and underscored the importance of developing and maintaining robust systems for accurate surveillance. As public health data collection efforts expand, there is a critical need for infectious disease modeling researchers to continue to develop prospective surveillance metrics and statistical models to accommodate the modeling of large disease counts and variability. This paper evaluated different likelihoods for the disease count model and various spatiotemporal mean models for prospective surveillance. METHODS We evaluated Bayesian spatiotemporal models, which are the foundation for model-based infectious disease surveillance metrics. Bayesian spatiotemporal mean models based on the Poisson and the negative binomial likelihoods were evaluated with the different lengths of past data usage. We compared their goodness of fit and short-term prediction performance with both simulated epidemic data and real data from the COVID-19 pandemic. RESULTS The simulation results show that the negative binomial likelihood-based models show better goodness of fit results than Poisson likelihood-based models as deemed by smaller deviance information criteria (DIC) values. However, Poisson models yield smaller mean square error (MSE) and mean absolute one-step prediction error (MAOSPE) results when we use a shorter length of the past data such as 7 and 3 time periods. Real COVID-19 data analysis of New Jersey and South Carolina shows similar results for the goodness of fit and short-term prediction results. Negative binomial-based mean models showed better performance when we used the past data of 52 time periods. Poisson-based mean models showed comparable goodness of fit performance and smaller MSE and MAOSPE results when we used the past data of 7 and 3 time periods. CONCLUSION We evaluate these models and provide future infectious disease outbreak modeling guidelines for Bayesian spatiotemporal analysis. Our choice of the likelihood and spatiotemporal mean models was influenced by both historical data length and variability. With a longer length of past data usage and more over-dispersed data, the negative binomial likelihood shows a better model fit than the Poisson likelihood. However, as we use a shorter length of the past data for our surveillance analysis, the difference between the Poisson and the negative binomial models becomes smaller. In this case, the Poisson likelihood shows robust posterior mean estimate and short-term prediction results.
Collapse
|
9
|
Rotejanaprasert C, Lawson AB, Maude RJ. Spatiotemporal reproduction number with Bayesian model selection for evaluation of emerging infectious disease transmissibility: an application to COVID-19 national surveillance data. BMC Med Res Methodol 2023; 23:62. [PMID: 36915077 PMCID: PMC10010957 DOI: 10.1186/s12874-023-01870-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Accepted: 02/20/2023] [Indexed: 03/16/2023] Open
Abstract
BACKGROUND To control emerging diseases, governments often have to make decisions based on limited evidence. The effective or temporal reproductive number is used to estimate the expected number of new cases caused by an infectious person in a partially susceptible population. While the temporal dynamic is captured in the temporal reproduction number, the dominant approach is currently based on modeling that implicitly treats people within a population as geographically well mixed. METHODS In this study we aimed to develop a generic and robust methodology for estimating spatiotemporal dynamic measures that can be instantaneously computed for each location and time within a Bayesian model selection and averaging framework. A simulation study was conducted to demonstrate robustness of the method. A case study was provided of a real-world application to COVID-19 national surveillance data in Thailand. RESULTS Overall, the proposed method allowed for estimation of different scenarios of reproduction numbers in the simulation study. The model selection chose the true serial interval when included in our study whereas model averaging yielded the weighted outcome which could be less accurate than model selection. In the case study of COVID-19 in Thailand, the best model based on model selection and averaging criteria had a similar trend to real data and was consistent with previously published findings in the country. CONCLUSIONS The method yielded robust estimation in several simulated scenarios of force of transmission with computing flexibility and practical benefits. Thus, this development can be suitable and practically useful for surveillance applications especially for newly emerging diseases. As new outbreak waves continue to develop and the risk changes on both local and global scales, our work can facilitate policymaking for timely disease control.
Collapse
|
10
|
Schildkraut JM, Johnson C, Dempsey LF, Qin B, Terry P, Akonde M, Peters ES, Mandle H, Cote ML, Peres L, Moorman P, Schwartz AG, Epstein M, Marks J, Bondy M, Lawson AB, Alberg AJ, Bandera EV. Survival of epithelial ovarian cancer in Black women: a society to cell approach in the African American cancer epidemiology study (AACES). Cancer Causes Control 2023; 34:251-265. [PMID: 36520244 PMCID: PMC9753020 DOI: 10.1007/s10552-022-01660-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Accepted: 12/05/2022] [Indexed: 12/23/2022]
Abstract
PURPOSE The causes for the survival disparity among Black women with epithelial ovarian cancer (EOC) are likely multi-factorial. Here we describe the African American Cancer Epidemiology Study (AACES), the largest cohort of Black women with EOC. METHODS AACES phase 2 (enrolled 2020 onward) is a multi-site, population-based study focused on overall survival (OS) of EOC. Rapid case ascertainment is used in ongoing patient recruitment in eight U.S. states, both northern and southern. Data collection is composed of a survey, biospecimens, and medical record abstraction. Results characterizing the survival experience of the phase 1 study population (enrolled 2010-2015) are presented. RESULTS Thus far, ~ 650 patients with EOC have been enrolled in the AACES. The five-year OS of AACES participants approximates those of Black women in the Surveillance Epidemiology and End Results (SEER) registry who survive at least 10-month past diagnosis and is worse compared to white women in SEER, 49 vs. 60%, respectively. A high proportion of women in AACES have low levels of household income (45% < $25,000 annually), education (51% ≤ high school education), and insurance coverage (32% uninsured or Medicaid). Those followed annually differ from those without follow-up with higher levels of localized disease (28 vs 24%) and higher levels of optimal debulking status (73 vs 67%). CONCLUSION AACES is well positioned to evaluate the contribution of social determinants of health to the poor survival of Black women with EOC and advance understanding of the multi-factorial causes of the ovarian cancer survival disparity in Black women.
Collapse
|
11
|
Lawson AB, Kim J. Bayesian space-time SIR modeling of Covid-19 in two US states during the 2020-2021 pandemic. PLoS One 2022; 17:e0278515. [PMID: 36548256 PMCID: PMC9778953 DOI: 10.1371/journal.pone.0278515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2022] [Accepted: 11/17/2022] [Indexed: 12/24/2022] Open
Abstract
This paper describes the Bayesian SIR modeling of the 3 waves of Covid-19 in two contrasting US states during 2020-2021. A variety of models are evaluated at the county level for goodness-of-fit and an assessment of confounding predictors is also made. It is found that models with three deprivation predictors and neighborhood effects are important. In addition, the work index from Google mobility was also found to provide an increased explanation of the transmission dynamics.
Collapse
|
12
|
Lawson AB, Kim J. Issues in Bayesian prospective surveillance of spatial health data. Spat Spatiotemporal Epidemiol 2022; 41:100431. [PMID: 35691635 DOI: 10.1016/j.sste.2021.100431] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/01/2020] [Revised: 04/30/2021] [Accepted: 04/30/2021] [Indexed: 10/21/2022]
Abstract
In this paper I review some of the major issues that arise when geo-referenced health data are to be the subject of prospective surveillance. The review focusses on modelbased approaches to this activity, and proposes the Bayesian paradigm as a convenient vehicle for modeling. Various posterior functional measures are discussed including the SCPO and SKL and a number of extensions to these are considered. Overall the value of Bayesian Hierarchical Modeling (BHM) with surveillance functionals is stressed in its relevance to early warning of adverse risk scenarios.
Collapse
|
13
|
Bozigar M, Lawson AB, Pearce JL, Svendsen ER, Vena JE. Using Bayesian time-stratified case-crossover models to examine associations between air pollution and "asthma seasons" in a low air pollution environment. PLoS One 2021; 16:e0260264. [PMID: 34879071 PMCID: PMC8654232 DOI: 10.1371/journal.pone.0260264] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2021] [Accepted: 11/05/2021] [Indexed: 11/18/2022] Open
Abstract
Many areas of the United States have air pollution levels typically below Environmental Protection Agency (EPA) regulatory limits. Most health effects studies of air pollution use meteorological (e.g., warm/cool) or astronomical (e.g., solstice/equinox) definitions of seasons despite evidence suggesting temporally-misaligned intra-annual periods of relative asthma burden (i.e., “asthma seasons”). We introduce asthma seasons to elucidate whether air pollutants are associated with seasonal differences in asthma emergency department (ED) visits in a low air pollution environment. Within a Bayesian time-stratified case-crossover framework, we quantify seasonal associations between highly resolved estimates of six criteria air pollutants, two weather variables, and asthma ED visits among 66,092 children ages 5–19 living in South Carolina (SC) census tracts from 2005 to 2014. Results show that coarse particulates (particulate matter <10 μm and >2.5 μm: PM10-2.5) and nitrogen oxides (NOx) may contribute to asthma ED visits across years, but are particularly implicated in the highest-burden fall asthma season. Fine particulate matter (<2.5 μm: PM2.5) is only associated in the lowest-burden summer asthma season. Relatively cool and dry conditions in the summer asthma season and increased temperatures in the spring and fall asthma seasons are associated with increased ED visit odds. Few significant associations in the medium-burden winter and medium-high-burden spring asthma seasons suggest other ED visit drivers (e.g., viral infections) for each, respectively. Across rural and urban areas characterized by generally low air pollution levels, there are acute health effects associated with particulate matter, but only in the summer and fall asthma seasons and differing by PM size.
Collapse
|
14
|
Babatunde OA, Zahnd WE, Eberth JM, Lawson AB, Adams SA, Boakye EA, Jefferson MS, Allen CG, Pearce JL, Li H, Halbert CH. Association between Neighborhood Social Deprivation and Stage at Diagnosis among Breast Cancer Patients in South Carolina. INTERNATIONAL JOURNAL OF ENVIRONMENTAL RESEARCH AND PUBLIC HEALTH 2021; 18:ijerph182211824. [PMID: 34831579 PMCID: PMC8625868 DOI: 10.3390/ijerph182211824] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/17/2021] [Revised: 11/01/2021] [Accepted: 11/04/2021] [Indexed: 11/25/2022]
Abstract
The purpose of this study was to examine the association between neighborhood social deprivation and individual-level characteristics on breast cancer staging in African American and white breast cancer patients. We established a retrospective cohort of patients with breast cancer diagnosed from 1996 to 2015 using the South Carolina Central Cancer Registry. We abstracted sociodemographic and clinical variables from the registry and linked these data to a county-level composite that captured neighborhood social conditions-the social deprivation index (SDI). Data were analyzed using chi-square tests, Student's t-test, and multivariable ordinal regression analysis to evaluate associations. The study sample included 52,803 female patients with breast cancer. Results from the multivariable ordinal regression model demonstrate that higher SDI (OR = 1.06, 95% CI: 1.02-1.10), African American race (OR = 1.35, 95% CI: 1.29-1.41), and being unmarried (OR = 1.17, 95% CI: 1.13-1.22) were associated with a distant stage at diagnosis. Higher tumor grade, younger age, and more recent year of diagnosis were also associated with distant-stage diagnosis. As a proxy for neighborhood context, the SDI can be used by cancer registries and related population-based studies to identify geographic areas that could be prioritized for cancer prevention and control efforts.
Collapse
|
15
|
Sartorius B, Lawson AB, Pullan RL. Author Correction: Modelling and predicting the spatio‑temporal spread of COVID‑19, associated deaths and impact of key risk factors in England. Sci Rep 2021; 11:17699. [PMID: 34465862 PMCID: PMC8406384 DOI: 10.1038/s41598-021-97282-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
16
|
Baer DR, Lawson AB, Joseph JE. Joint space-time Bayesian disease mapping via quantification of disease risk association. Stat Methods Med Res 2021; 30:35-61. [PMID: 33595403 DOI: 10.1177/0962280220938975] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Alzheimer's disease is an increasingly prevalent neurological disorder with no effective therapies. Thus, there is a need to characterize the progression of Alzheimer's disease risk in order to preclude its inception in patients. Characterizing Alzheimer's disease risk can be accomplished at the population-level by the space-time modeling of Alzheimer's disease incidence data. In this paper, we develop flexible Bayesian hierarchical models which can borrow risk information from conditions antecedent to Alzheimer's disease, such as mild cognitive impairment, in an effort to better characterize Alzheimer's disease risk over space and time. From an application of these models to real-world Alzheimer's disease and mild cognitive impairment spatiotemporal incidence data, we found that our novel models provided improved model goodness of fit, and via a simulation study, we demonstrated the importance of diagnosing the label-switching problem for our models as well as the importance of model specification in order to best capture the contribution of time in modeling Alzheimer's disease risk.
Collapse
|
17
|
Lawson AB, Prates M, Anderson C. GEOMED 2019 Editorial. Stat Methods Med Res 2021; 30:5. [PMID: 33595404 DOI: 10.1177/0962280220930177] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
18
|
Sartorius B, Lawson AB, Pullan RL. Modelling and predicting the spatio-temporal spread of COVID-19, associated deaths and impact of key risk factors in England. Sci Rep 2021; 11:5378. [PMID: 33686125 PMCID: PMC7940626 DOI: 10.1038/s41598-021-83780-2] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2020] [Accepted: 02/08/2021] [Indexed: 12/13/2022] Open
Abstract
COVID-19 caseloads in England have passed through a first peak, and at the time of this analysis appeared to be gradually increasing, potentially signalling the emergence of a second wave. To ensure continued response to the epidemic is most effective, it is imperative to better understand both retrospectively and prospectively the geographical evolution of COVID-19 caseloads and deaths at small-area resolution, identify localised areas in space-time at significantly higher risk, quantify the impact of changes in localised population mobility (or movement) on caseloads, identify localised risk factors for increased mortality and project the likely course of the epidemic at high spatial resolution in coming weeks. We applied a Bayesian hierarchical space-time SEIR model to assess the spatiotemporal variability of COVID-19 caseloads (transmission) and deaths at small-area scale in England [Middle Layer Super Output Area (MSOA), 6791 units] and by week (using observed data from week 5 to 34 of 2020), including key determinants, the modelled transmission dynamics and spatial-temporal random effects. We also estimate the number of cases and deaths at small-area resolution with uncertainty projected forward in time by MSOA (up to week 51 of 2020), the impact mobility reductions (and subsequent easing) have had on COVID-19 caseloads and quantify the impact of key socio-demographic risk factors on COVID-19 related mortality risk by MSOA. Reductions in population mobility during the course of the first lockdown had a significant impact on the reduction of COVID-19 caseloads across England, however local authorities have had a varied rate of reduction in population movement which our model suggest has substantially impacted the geographic heterogeneity in caseloads at small-area scale. The steady gain in population mobility, observed from late April, appears to have contributed to a slowdown in caseload reductions towards late June and subsequent start of the second wave. MSOA with higher proportions of elderly (70+ years of age) and elderly living in deprivation, both with very distinct geographic distributions, have a significantly elevated COVID-19 mortality rates. While non-pharmaceutical interventions (that is, reductions in population mobility and social distancing) had a profound impact on the trajectory of the first wave of the COVID-19 outbreak in England, increased population mobility appears to have significantly contributed to the second wave. A number of contiguous small-areas appear to be at a significant elevated risk of high COVID-19 transmission, many of which are also at increased risk for higher mortality rates. A geographically staggered re-introduction of intensified social distancing measures is advised and limited cross MSOA movement if the magnitude and geographic extent of the second wave is to be reduced.
Collapse
|
19
|
Davis ML, Neelon B, Nietert PJ, Burgette LF, Hunt KJ, Lawson AB, Egede LE. Propensity score matching for multilevel spatial data: accounting for geographic confounding in health disparity studies. Int J Health Geogr 2021; 20:10. [PMID: 33639940 PMCID: PMC7913404 DOI: 10.1186/s12942-021-00265-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Accepted: 02/10/2021] [Indexed: 11/20/2022] Open
Abstract
BACKGROUND Diabetes is a public health burden that disproportionately affects military veterans and racial minorities. Studies of racial disparities are inherently observational, and thus may require the use of methods such as Propensity Score Analysis (PSA). While traditional PSA accounts for patient-level factors, this may not be sufficient when patients are clustered at the geographic level and thus important confounders, whether observed or unobserved, vary by geographic location. METHODS We employ a spatial propensity score matching method to account for "geographic confounding", which occurs when the confounding factors, whether observed or unobserved, vary by geographic region. We augment the propensity score and outcome models with spatial random effects, which are assigned scaled Besag-York-Mollié priors to address spatial clustering and improve inferences by borrowing information across neighboring geographic regions. We apply this approach to a study exploring racial disparities in diabetes specialty care between non-Hispanic black and non-Hispanic white veterans. We construct multiple global estimates of the risk difference in diabetes care: a crude unadjusted estimate, an estimate based solely on patient-level matching, and an estimate that incorporates both patient and spatial information. RESULTS In simulation we show that in the presence of an unmeasured geographic confounder, ignoring spatial heterogeneity results in increased relative bias and mean squared error, whereas incorporating spatial random effects improves inferences. In our study of racial disparities in diabetes specialty care, the crude unadjusted estimate suggests that specialty care is more prevalent among non-Hispanic blacks, while patient-level matching indicates that it is less prevalent. Hierarchical spatial matching supports the latter conclusion, with a further increase in the magnitude of the disparity. CONCLUSIONS These results highlight the importance of accounting for spatial heterogeneity in propensity score analysis, and suggest the need for clinical care and management strategies that are culturally sensitive and racially inclusive.
Collapse
|
20
|
Lawson AB, Hughes-Halbert C, Babatunde OA, Zahnd WE, Eberth JM. Abstract PO-167: Area-level social deprivation and stage at diagnosis among breast cancer patients in South Carolina. Cancer Epidemiol Biomarkers Prev 2020. [DOI: 10.1158/1538-7755.disp20-po-167] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Abstract
Background: Cancer incidence, staging and mortality rates differ across geographic areas; however, there is a need for a better understanding of how neighborhood level socioeconomic and access to care factors impact cancer burden to tailor cancer control interventions in ways that appropriately target geographic determinants of cancer health disparities. The purpose of this study was to characterize the distribution of neighborhood deprivation in a cohort of breast cancer patients and examine the effect of social deprivation, healthcare professional shortage Area (HPSA) designation, and individual-level characteristics on breast cancer staging.
Methods: We established a retrospective cohort of breast cancer patients diagnosed from 1996 to 2015 using the South Carolina Central Cancer Registry.
Sociodemographic (e.g., race, age) and clinical variables were abstracted from the registry. We linked registry data to county-level variables to determine levels of social deprivation and residence in a health care professional shortage area using the Robert Graham Center’s Social Deprivation Index (SDI) and the Health Resources and Services Administration primary care HPSA designation. Bivariate analyses and multivariate regression analyses were used to examine associations. Results: The sample included 54,501 female breast cancer patients. Overall, the mean for SDI was 54.2 (+18.1) and the range was 76 (19-95). Approximately 44.4% of women lived in areas with high levels of social deprivation (e.g., SDI score of 52 to 95). In the logistic regression model, living in a geographic area with high social deprivation was significantly associated with African American race (OR=2.3, 95% C.I. 2.2-2.4), being unmarried (OR=1.2, 95% C.I. 1.1-1.3), and HPSA designation (OR=14.0, 95% C.I. 13.5- 14.6). Higher tumor grade (OR=1.2; 95% CI.1.2-1.3) and later stage (OR=1.1, 95% C.I. 1.1-1.2) were also significantly associated with neighborhood deprivation. Conclusion: This study shows that SDI differs by race and clinical characteristics among breast cancer patients. The SDI could be integrated into tumor registries and cancer research to understand the effects of neighborhood level variables on cancer health disparities to improve the precision of cancer control interventions that are developed to address geographic determinants.
Citation Format: Andrew B. Lawson, Chanita Hughes-Halbert, Oluwole A. Babatunde, Whitney E. Zahnd, Jan M. Eberth. Area-level social deprivation and stage at diagnosis among breast cancer patients in South Carolina [abstract]. In: Proceedings of the AACR Virtual Conference: Thirteenth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2020 Oct 2-4. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(12 Suppl):Abstract nr PO-167.
Collapse
|
21
|
Bozigar M, Lawson AB, Pearce JL, King K, Svendsen ER. A Bayesian spatio-temporal analysis of neighborhood pediatric asthma emergency department visit disparities. Health Place 2020; 66:102426. [PMID: 33011491 DOI: 10.1016/j.healthplace.2020.102426] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 07/17/2020] [Accepted: 08/17/2020] [Indexed: 11/25/2022]
Abstract
Asthma disparities have complex, neighborhood-level drivers that are not well understood. Consequently, identifying particular contextual factors that contribute to disparities is a public health goal. We study pediatric asthma emergency department (ED) visit disparities and neighborhood factors associated with them in South Carolina (SC) census tracts from 1999 to 2015. Leveraging a Bayesian framework, we identify risk clusters, spatially-varying relationships, and risk percentile-specific associations. Clusters of high risk occur in both rural and urban census tracts with high probability, with neighborhood-specific associations suggesting unique risk factors for each locale. Bayesian methods can help clarify the neighborhood drivers of health disparities.
Collapse
|
22
|
Boaz RM, Lawson AB, Pearce JL. Multivariate Air Pollution Prediction Modeling with partial Missingness. ENVIRONMETRICS 2019; 30:e2592. [PMID: 31983873 PMCID: PMC6980235 DOI: 10.1002/env.2592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2018] [Accepted: 05/23/2019] [Indexed: 06/10/2023]
Abstract
Missing observations from air pollution monitoring networks have posed a longstanding problem for health investigators of air pollution. Growing interest in mixtures of air pollutants has further complicated this problem, as many new challenges have arisen that require development of novel methods. The objective of this study is to develop a methodology for multivariate prediction of air pollution. We focus specifically on tackling different forms of missing data, such as: spatial (sparse sites), outcome (pollutants not measured at some sites), and temporal (varieties of interrupted time series). To address these challenges, we develop a novel multivariate fusion framework, which leverages the observed inter-pollutant correlation structure to reduce error in the simultaneous prediction of multiple air pollutants. Our joint fusion model employs predictions from the Environmental Protection Agency's Community Multiscale Air Quality (CMAQ) model along with spatio-temporal error terms. We have implemented our models on both simulated data and a case study in South Carolina for 8 pollutants over a 28-day period in June 2006. We found that our model, which uses a multivariate correlated error in a Bayesian framework, showed promising predictive accuracy particularly for gaseous pollutants.
Collapse
|
23
|
Reed SG, Miller CS, Wagner CL, Hollis BW, Lawson AB. Toward Preventing Enamel Hypoplasia: Modeling Maternal and Neonatal Biomarkers of Human Calcium Homeostasis. Caries Res 2019; 54:55-67. [PMID: 31665727 PMCID: PMC7299520 DOI: 10.1159/000502793] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 08/15/2019] [Indexed: 01/01/2023] Open
Abstract
AIM The aim of this study was to assess biomarkers of calcium homeostasis and tooth development, in mothers during pregnancy and their children at birth, for enamel hypoplasia (EH) in the primary maxillary central incisor teeth. METHODS Bayesian methodology was used for secondary data analyses from a randomized, controlled trial of prenatal vitamin D3 supplementation in healthy mothers (N = 350) and a follow-up study of a subset of the children. The biomarkers were serum calcium (Ca), phosphorus (P), intact parathyroid hormone (iPTH), total circulating 25-dihydroxyvitamin D (25(OH)D), and 1,25-dihydroxyvitamin D (1,25(OH)2D). The maternal biomarkers were assayed monthly during pregnancy, and the child's biomarkers were derived from cord blood. Digital images of the child's 2 teeth were scored for EH using Enamel Defects Index criteria for each of the incisal, middle, and cervical regions for an EH extent score. RESULTS The child EH prevalence was 41% (60/145), with most defects present in the incisal and middle tooth regions. Cord blood iPTH and 1,25(OH)2D levels were significantly associated with EH extent after controlling for maternal factors. For every 1 pg/mL increase in cord blood iPTH, the EH extent decreased by approximately 6%. For every 10 pg/mL increase in cord blood 1,25(OH)2D, the EH extent increased by almost 30% (holding all other terms constant and adjusting for subject-level heterogeneity). The relationship between maternal 25(OH)D and maternal mean iPTH varied significantly by EH extent. CONCLUSION The results suggest possible modifiable relationships of maternal and neonatal factors of calcium homeostasis during pregnancy and at birth for EH, contributing to the frontier of knowledge regarding sound tooth development for dental caries prevention.
Collapse
|
24
|
Rotejanaprasert C, Lawson AB, Iamsirithaworn S. Spatiotemporal multi-disease transmission dynamic measure for emerging diseases: an application to dengue and zika integrated surveillance in Thailand. BMC Med Res Methodol 2019; 19:200. [PMID: 31655546 PMCID: PMC6815359 DOI: 10.1186/s12874-019-0833-6] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2018] [Accepted: 09/12/2019] [Indexed: 11/26/2022] Open
Abstract
Background New emerging diseases are public health concerns in which policy makers have to make decisions in the presence of enormous uncertainty. This is an important challenge in terms of emergency preparation requiring the operation of effective surveillance systems. A key concept to investigate the dynamic of infectious diseases is the basic reproduction number. However it is difficult to be applicable in real situations due to the underlying theoretical assumptions. Methods In this paper we propose a robust and flexible methodology for estimating disease strength varying in space and time using an alternative measure of disease transmission within the hierarchical modeling framework. The proposed measure is also extended to allow for incorporating knowledge from related diseases to enhance performance of surveillance system. Results A simulation was conducted to examine robustness of the proposed methodology and the simulation results demonstrate that the proposed method allows robust estimation of the disease strength across simulation scenarios. A real data example is provided of an integrative application of Dengue and Zika surveillance in Thailand. The real data example also shows that combining both diseases in an integrated analysis essentially decreases variability of model fitting. Conclusions The proposed methodology is robust in several simulated scenarios of spatiotemporal transmission force with computing flexibility and practical benefits. This development has potential for broad applicability as an alternative tool for integrated surveillance of emerging diseases such as Zika.
Collapse
|
25
|
Carroll R, Lawson AB, Zhao S. Temporally dependent accelerated failure time model for capturing the impact of events that alter survival in disease mapping. Biostatistics 2019; 20:666-680. [PMID: 29939209 DOI: 10.1093/biostatistics/kxy023] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2017] [Revised: 03/08/2018] [Accepted: 04/24/2018] [Indexed: 11/15/2022] Open
Abstract
The introduction of spatial and temporal frailty parameters in survival models furnishes a way to represent unmeasured confounding in the outcome of interest. Using a Bayesian accelerated failure time model, we are able to flexibly explore a wide range of spatial and temporal options for structuring frailties as well as examine the benefits of using these different structures in certain settings. A setting of particular interest for this work involved using temporal frailties to capture the impact of events of interest on breast cancer survival. Our results suggest that it is important to include these temporal frailties when there is a true temporal structure to the outcome and including them when a true temporal structure is absent does not sacrifice model fit. Additionally, the frailties are able to correctly recover the truth imposed on simulated data without affecting the fixed effect estimates. In the case study involving Louisiana breast cancer-specific mortality, the temporal frailty played an important role in representing the unmeasured confounding related to improvements in knowledge, education, and disease screenings as well as the impacts of Hurricane Katrina and the passing of the Affordable Care Act. In conclusion, the incorporation of temporal, in addition to spatial, frailties in survival analysis can lead to better fitting models and improved inference by representing both spatially and temporally varying unmeasured risk factors and confounding that could impact survival. Specifically, we successfully estimated changes in survival around the time of events of interest.
Collapse
|