Widely Applicable Information Criterion Research Articles

Cross-Cohort Mixture Analysis: A Data Integration Approach With Applications on Gestational Age and DNA-Methylation-Derived Gestational Age Acceleration Metrics.

Data integration of multiple studies can provide enhanced exposure contrast and statistical power to examine associations between environmental exposure mixtures and health outcomes. Extant research has combined populations and identified an overall mixture-outcome association, without accounting for differences across studies. We extended the Bayesian Weighted Quantile Sum (BWQS) regression to a hierarchical framework to analyze mixtures across cohorts. The hierarchical BWQS (HBWQS) approach aggregates sample size of multiple cohorts to calculate an overall mixture index, thereby identifying the most harmful exposure(s) across cohorts; and provides cohort-specific associations between the overall mixture index and the outcome. We showed results from 10 simulated scenarios including four mixture components in three, eight, and ten populations, and two real-case examples on the association between prenatal metal mixture exposure-comprising arsenic, cadmium, and lead-and both gestational age and epigenetic-derived gestational age acceleration metrics. Simulated scenarios showed good empirical coverage and little bias for all HBWQS-estimated parameters. The Watanabe-Akaike information criterion showed a better average performance for the HBWQS regression than the BWQS across scenarios. HBWQS results incorporating cohorts within the national Environmental influences on Child Health Outcomes (ECHO) program from three different sites showed that the environmental mixture was negatively associated with gestational age in a single site. The HBWQS approach facilitates the combination of multiple cohorts and accounts for individual cohort differences in mixture analyses. HBWQS findings can be used to develop regulations, policies, and interventions regarding multiple co-occurring environmental exposures and it will maximize the use of extant publicly available data.

Hematocrits of red blood cell units increase during storage due to changes in mean corpuscular volume, impacting outcomes of red cell exchange procedures

Abstract Introduction Hematocrits of red blood cell (RBC) units are generally not measured by blood centers or transfusion services. Red blood cell exchange (RCE) is commonly used to treat sickle cell disease, and a key parameter needed in RCE is the average hematocrit of RBC units. However, there is no consensus on the average hematocrit to use or how it is affected by storage. The goal of our study is to determine the average hematocrit of RBC units by measuring hematologic parameters RBC units during storage and performing statistical modeling of hematocrit changes over time. Methods We sampled aliquots of leukoreduced RBC units in AS-1 (n=10) using the Terumo TSCD-II Sterile Tubing Welder to sterilely collect ~15mL of RBCs in transfer-pack containers. Aliquots were collected weekly from Day 7 to Day 42 of storage. Complete blood counts were performed using the Sysmex XN-3100 Hematology Analyzer. We statistically modeled changes in RBC units during storage by utilizing Bayesian multilevel linear regression models, allowing for capturing effects at the individual RBC-unit level. Posterior distributions of model parameters were estimated and processed using Stan and the rethinking and rstan packages in R. Posterior distributions are summarized with posterior means and 95%-credible intervals in parentheses. Model selection was done by calculating Akaike weights using the widely applicable information criterion (WAIC) to signify statistical support for each model. Results Multilevel linear regression models of hematologic parameters as a function of storage time were developed. Hematocrits increased with a posterior mean of 1.4% (0.98–1.75%) per week of storage, with posterior mean hematocrits of: 61.9% (60.2–63.6%) on Day 7, 63.9% (61.5–66.4%) on Day 21, and 66.0% (62.8–69.7%) on Day 42. These results indicate that a single average cannot be used for RBC units since hematocrits dramatically increase during storage. Mean corpuscular volume (MCV) increased with a posterior mean of 2.31 fL (1.4–3.71) per week whereas hemoglobin remained stable with a posterior mean increase of 0.008 g/dL (-0.016–0.031) per week. Using these results, we subsequently developed a multilevel regression model of hematocrit as a function of MCV. The model using MCV (WAIC=100.2, weight=100%) was favored over the model using time (WAIC=120.0, weight=0%), suggesting that changes in MCV are driving the increasing hematocrits. Conclusion Our study is the first to leverage Bayesian statistics to carefully quantify the average hematocrits of RBC units during storage time and to demonstrate that increasing hematocrits are caused by changes in MCV. These results indicate that transfusion protocols must consider the age of RBC units for RCE calculations, and provisions (e.g., requiring use of fresh units for RCE) may help minimize the variability in post-exchange testing. Moreover, our findings suggest that transfusion services should regularly measure hematocrits of RBC units in their inventory.

Modelling of the time to death of breast cancer patients at Hiwot Fana Specialized University Hospital

Breast cancer is the most common cause of cancer death and is a frequently diagnosed cancer among women worldwide. It is becoming a challenging health condition in Ethiopia with a high rate of morbidity and mortality. The main aim of this study was to model the time to death in breast cancer patients at Hiwot Fana Specialized University Hospital. A retrospective cohort study was carried out from April 1st, 2020, to April 1st, 2023, and 296 women were included in the study. We used nonparametric methods and Bayesian accelerated failure time models (with Laplace approximation) to identify risk factors and choose a model fitting breast cancer patient data. Model comparison was performed using the marginal likelihood, deviance information criterion and Watanabe Akaike information criterion. From the total of 296 patients in the study, 56 (18.9%) died. The estimated median survival time was 33 months. The log-rank test showed that age group, stage, alcohol consumption, smoking habit, and comorbidity were potential risk factors associated with the time to death in breast cancer patients at the 5% level of significance. The Bayesian Weibull accelerated failure time model was found to be the best fitted model for predicting the survival time of patients with minimum DIC (520.39) and WAIC (521.59) values. The final Bayesian Weibull AFT model with the integrated nested Laplace approximation estimation technique revealed that age group, stage, alcohol consumption, smoking habit, and comorbidity were significantly associated with the time to death in breast cancer patients. Individuals older than 65 years, with stage IV disease, drinking alcohol, smoking cigarettes and having comorbidities had shortened survival times in patients with breast cancer. Hence, Hiwot Fana Specialized University Hospital and related bodies should work on awareness creation to reduce smoking habits and alcohol use as well as give due attention to elderly and stage IV breast cancer patients during intervention.

Are information criteria good enough to choose the right the number of regimes in hidden Markov models?

Selecting the number of regimes in Hidden Markov models is an important problem. Many criteria are used to do so, such as Akaike information criterion (AIC), Bayesian information criterion (BIC), integrated completed likelihood (ICL), deviance information criterion (DIC), and Watanabe-Akaike information criterion (WAIC), to name a few. In this article, we introduced goodness-of-fit tests for general Hidden Markov models with covariates, where the distribution of the observations is arbitrary, i.e. continuous, discrete, or a mixture of both. A selection procedure is proposed based on this goodness-of-fit test. The main aim of this article is to compare the classical information criterion with the new criterion when the outcome is either continuous, discrete or zero-inflated. Numerical experiments assess the finite sample performance of the goodness-of-fit tests, and comparisons between the different criteria are made.

Patient- and Areal-Level Risk Factors Associated With Lung Cancer Mortality in Victoria, Australia: A Bayesian Spatial Survival Analysis.

In Australia, lung cancer is the leading cause of cancer-related deaths. In Victoria, the mortality risk is assumed to vary across Local Government Areas (LGAs) due to variations in socioeconomic advantage, remoteness, and healthcare accessibility. Thus, we applied Bayesian spatial survival models to examine the geographic variation in lung cancer survival in Victoria. Data on lung cancer cases were extracted from the Victorian Lung Cancer Registry (VLCR). To account for spatial dependence and risk factors of survival in lung cancer patients, we employed a Bayesian spatial survival model. Conditional Autoregressive (CAR) prior was assigned to model the spatial dependence. Deviance Information Criterion (DIC), Watanabe Akaike Information Criterion (WAIC), and Log Pseudo Marginal Likelihood (LPML) were used for model comparison. In the final best-fitted model, the Adjusted Hazard Ratio (AHR) with the 95% Credible Interval (CrI) was reported. The outcome variable was the survival status of lung cancer patients, defined as whether they survived or died during the follow-up period (death was our interest). Our study revealed substantial variations in lung cancer mortality in Victoria. Poor Eastern Cooperative Oncology Group (ECOG) performance status, diagnosed at a regional hospital, Small Cell Lung Cancer (SCLC), advanced age, and advanced clinical stage were associated with a higher risk of mortality, whereas being female, presented at Multidisciplinary Team (MDT) meeting, and diagnosed at a metropolitan private hospital were significantly associated with a lower risk of mortality. Identifying geographical disparities in lung cancer survival may help shape healthcare policy to implement more targeted and effective lung cancer care services.

Learning under singularity: an information criterion improving WBIC and sBIC

AbstractWe introduce a novel information criterion, termed learning under singularity (LS), designed to enhance the functionality of the widely applicable Bayes information criterion (WBIC) and the singular Bayesian information criterion (sBIC). LS is effective without regularity constraints and demonstrates stability. Watanabe defined a statistical model or a learning machine as regular if the mapping from a parameter to a probability distribution is one-to-one and its Fisher information matrix is positive definite. In contrast, models not meeting these conditions are termed singular. Over the past decade, several information criteria for singular cases have been proposed, including WBIC and sBIC. WBIC is applicable in non-regular scenarios but faces challenges with large sample sizes and redundant estimation of known learning coefficients. Conversely, sBIC is limited in its broader application due to its dependence on maximum likelihood estimates. LS addresses these limitations by enhancing the utility of both WBIC and sBIC. It incorporates the empirical loss from the widely applicable information criterion (WAIC) to represent the goodness of fit to the statistical model, along with a penalty term similar to that of sBIC. This approach offers a flexible and robust method for model selection, free from regularity constraints.

Bayesian Spatio-Temporal Multilevel Modelling of Patient-Reported Quality of Life following Prostate Cancer Surgery.

Globally, prostate cancer is the second leading cause of cancer deaths among males. It is the most commonly diagnosed cancer in Australia. The quality of life of prostate cancer patients is poorer when compared to the general population due to the disease itself and its related complications. However, there is limited research on the geographic pattern of quality of life and its risk factors in Victoria. Therefore, an examination of the spatio-temporal pattern and risk factors of poor quality of life, along with the impact of spatial weight matrices on estimates and model performance, was conducted. A retrospective study was undertaken based on the Prostate Cancer Outcome Registry-Victoria data. Patient data (n = 5238) were extracted from the Prostate Cancer Outcome Registry, a population-based clinical quality outcome assessment from 2015 to 2021. A Bayesian spatio-temporal multilevel model was fitted to identify risk factors for poor quality of life. This study also evaluated the impact of distance- and adjacency-based spatial weight matrices. Model convergence was assessed using Gelman-Rubin statistical plots, and model comparison was based on the Watanabe-Akaike Information Criterion. A total of 1906 (36.38%) prostate cancer patients who had undergone surgery experienced poor quality of life in our study. Belonging to the age group between 76 and 85 years (adjusted odds ratio (AOR) = 2.90, 95% credible interval (CrI): 1.39, 2.08), having a prostate-specific antigen level between 10.1 and 20.0 (AOR = 1.33, 95% CrI: 1.12, 1.58), and being treated in a public hospital (AOR = 1.35, 95% CrI: 1.17, 1.53) were significantly associated with higher odds of poor quality of life. Conversely, residing in highly accessible areas (AOR = 0.60, 95% CrI: 0.38, 0.94) was significantly associated with lower odds of poor prostate-specific antigen levels. Variations in estimates and model performance were observed depending on the choice of spatial weight matrices. Belonging to an older age group, having a high prostate-specific antigen level, receiving treatment in public hospitals, and remoteness were statistically significant factors linked to poor quality of life. Substantial spatio-temporal variations in poor quality of life were observed in Victoria across local government areas. The distance-based weight matrix performed better than the adjacency-based matrix. This research finding highlights the need to reduce geographical disparities in quality of life. The statistical methods developed in this study may also be useful to apply to other population-based clinical registry settings.

Optimizing spatio-temporal correlation structures for modeling food security in Africa: a simulation-based investigation

This study investigates the impact of spatio- temporal correlation using four spatio-temporal models: Spatio-Temporal Poisson Linear Trend Model (SPLTM), Poisson Temporal Model (TMS), Spatio-Temporal Poisson Anova Model (SPAM), and Spatio-Temporal Poisson Separable Model (STSM) concerning food security and nutrition in Africa. Evaluating model goodness of fit using the Watanabe Akaike Information Criterion (WAIC) and assessing bias through root mean square error and mean absolute error values revealed a consistent monotonic pattern. SPLTM consistently demonstrates a propensity for overestimating food security, while TMS exhibits a diverse bias profile, shifting between overestimation and underestimation based on varying correlation settings. SPAM emerges as a beacon of reliability, showcasing minimal bias and WAIC across diverse scenarios, while STSM consistently underestimates food security, particularly in regions marked by low to moderate spatio-temporal correlation. SPAM consistently outperforms other models, making it a top choice for modeling food security and nutrition dynamics in Africa. This research highlights the impact of spatial and temporal correlations on food security and nutrition patterns and provides guidance for model selection and refinement. Researchers are encouraged to meticulously evaluate the biases and goodness of fit characteristics of models, ensuring their alignment with the specific attributes of their data and research goals. This knowledge empowers researchers to select models that offer reliability and consistency, enhancing the applicability of their findings.

A flexible Bayesian tool for CoDa mixed models: logistic-normal distribution with Dirichlet covariance

Compositional Data Analysis (CoDa) has gained popularity in recent years. This type of data consists of values from disjoint categories that sum up to a constant. Both Dirichlet regression and logistic-normal regression have become popular as CoDa analysis methods. However, fitting this kind of multivariate models presents challenges, especially when structured random effects are included in the model, such as temporal or spatial effects. To overcome these challenges, we propose the logistic-normal Dirichlet Model (LNDM). We seamlessly incorporate this approach into the R-INLA package, facilitating model fitting and model prediction within the framework of Latent Gaussian Models. Moreover, we explore metrics like Deviance Information Criteria, Watanabe Akaike information criterion, and cross-validation measure conditional predictive ordinate for model selection in R-INLA for CoDa. Illustrating LNDM through two simulated examples and with an ecological case study on Arabidopsis thaliana in the Iberian Peninsula, we underscore its potential as an effective tool for managing CoDa and large CoDa databases.

Estimating the Distribution of True Rates of Visual Field Progression in Glaucoma.

The purpose of this study was to estimate the distribution of the true rates of progression (RoP) of visual field (VF) loss. We analyzed the progression of mean deviation over time in series of ≥ 10 tests from 3352 eyes (one per patient) from 5 glaucoma clinics, using a novel Bayesian hierarchical Linear Mixed Model (LMM); this modeled the random-effect distribution of RoPs as the sum of 2 independent processes following, respectively, a negative exponential distribution (the "true" distribution of RoPs) and a Gaussian distribution (the "noise"), resulting in a skewed exGaussian distribution. The exGaussian-LMM was compared to a standard Gaussian-LMM using the Watanabe-Akaike Information Criterion (WAIC). The random-effect distributions were compared to the empirical cumulative distribution function (eCDF) of linear regression RoPs using a Kolmogorov-Smirnov test. The WAIC indicated a better fit with the exGaussian-LMM (estimate [standard error]: 192174.4 [721.2]) than with the Gaussian-LMM (192595 [697.4], with a difference of 157.2 [22.6]). There was a significant difference between the eCDF and the Gaussian-LMM distribution (P < 0.0001), but not with the exGaussian-LMM distribution (P = 0.108). The estimated mean (95% credible intervals, CIs) "true" RoP (-0.377, 95% CI = -0.396 to -0.359 dB/year) was more negative than the observed mean RoP (-0.283, 95% CI = -0.299 to -0.268 dB/year), indicating a bias likely due to learning in standard LMMs. The distribution of "true" RoPs can be estimated with an exGaussian-LMM, improving model accuracy. We used these results to develop a fast and accurate analytical approximation for sample-size calculations in clinical trials using standard LMMs, which was integrated in a freely available web application.

Modeling Factors Influencing Covid-19 Cases in South Sulawesi Using Bayesian Conditional Autoregressive Localised

South Sulawesi Province is listed as the province with the highest number of Covid-19 cases in the Sulawes island. Research on Covid-19 modeling has been carried out by many researchers, but until now, there has been no research using the Bayesian spatial Conditional Autoregressive Localized model which involves a combination of factors such as distance to the provincial capital, population density, and the number of elderly people in each district in South Sulawesi Province. The aim of this research is to get the best Bayesian Conditional Autoregressive Localized model. The best model is based on four criteria, namely: Deviance Information Criteria, Watanabe Akaike Information Criteria, residuals from Modified Moran's I, and the number of areas included in a group. It was found that model with G=3 by including population density covariates was the best model. A significant factor influencing the increase in Covid-19 cases is the population density factor which has a positive effect. This shows that the more densely populated an area is, the greater the chance of being infected with Covid-19. Makassar has the highest relative risk value for Covid-19 followed by Toraja district and Pare-Pare City. Meanwhile, Bone district has the lowest relative risk value for Covid-19, followed by Wajo district and Enrekang district.

A study of model selection techniques for modeling synthetic aperture sonar backscatter statistics

High resolution acoustic imaging of the seafloor, such as with synthetic aperture sonar, can reveal complex environments due to the ability to resolve scatterers of different types. The statistical distribution of the backscattered field is appropriately modeled using a mixture distribution, which consists of a sum of a K pdf for each scatterer type, weighted by the relative frequency with which they occur. In this type of modeling, the number of distributions must be selected before parameters can be fit, and as the number of components increases, so does the danger of overfitting the data. Several methods of model selection are explored here. Two methods, the Bayesian information criterion and the Akaike information criterion are based on point parameter estimates, and a penalty due to the number of model parameters. The other two, the deviance information criterion and the Watanabe-Akaike information criterion, are based on Monte-Carlo sampling of the model parameter space. These four model selection criteria are compared to each other, and to the authors manual selection of the number of components.

Modelling spatiotemporal patterns of visceral leishmaniasis incidence in two endemic states in India using environment, bioclimatic and demographic data, 2013-2022.

As of 2021, the National Kala-azar Elimination Programme (NKAEP) in India has achieved visceral leishmaniasis (VL) elimination (<1 case / 10,000 population/year per block) in 625 of the 633 endemic blocks (subdistricts) in four states. The programme needs to sustain this achievement and target interventions in the remaining blocks to achieve the WHO 2030 target of VL elimination as a public health problem. An effective tool to analyse programme data and predict/ forecast the spatial and temporal trends of VL incidence, elimination threshold, and risk of resurgence will be of use to the programme management at this juncture. We employed spatiotemporal models incorporating environment, climatic and demographic factors as covariates to describe monthly VL cases for 8-years (2013-2020) in 491 and 27 endemic and non-endemic blocks of Bihar and Jharkhand states. We fitted 37 models of spatial, temporal, and spatiotemporal interaction random effects with covariates to monthly VL cases for 6-years (2013-2018, training data) using Bayesian inference via Integrated Nested Laplace Approximation (INLA) approach. The best-fitting model was selected based on deviance information criterion (DIC) and Watanabe-Akaike Information Criterion (WAIC) and was validated with monthly cases for 2019-2020 (test data). The model could describe observed spatial and temporal patterns of VL incidence in the two states having widely differing incidence trajectories, with >93% and 99% coverage probability (proportion of observations falling inside 95% Bayesian credible interval for the predicted number of VL cases per month) during the training and testing periods. PIT (probability integral transform) histograms confirmed consistency between prediction and observation for the test period. Forecasting for 2021-2023 showed that the annual VL incidence is likely to exceed elimination threshold in 16-18 blocks in 4 districts of Jharkhand and 33-38 blocks in 10 districts of Bihar. The risk of VL in non-endemic neighbouring blocks of both Bihar and Jharkhand are less than 0.5 during the training and test periods, and for 2021-2023, the probability that the risk greater than 1 is negligible (P<0.1). Fitted model showed that VL occurrence was positively associated with mean temperature, minimum temperature, enhanced vegetation index, precipitation, and isothermality, and negatively with maximum temperature, land surface temperature, soil moisture and population density. The spatiotemporal model incorporating environmental, bioclimatic, and demographic factors demonstrated that the KAMIS database of the national programmme can be used for block level predictions of long-term spatial and temporal trends in VL incidence and risk of outbreak / resurgence in endemic and non-endemic settings. The database integrated with the modelling framework and a dashboard facility can facilitate such analysis and predictions. This could aid the programme to monitor progress of VL elimination at least one-year ahead, assess risk of resurgence or outbreak in post-elimination settings, and implement timely and targeted interventions or preventive measures so that the NKAEP meet the target of achieving elimination by 2030.

The Best Fit Bayesian Hierarchical Generalized Linear Model Selection Using Information Complexity Criteria in the MCMC Approach

Both frequentist and Bayesian statistics schools have improved statistical tools and model choices for the collected data or measurements. Model selection approaches have advanced due to the difficulty of comparing complicated hierarchical models in which linear predictors vary by grouping variables, and the number of model parameters is not distinct. Many regression model selection criteria are considered, including the maximum likelihood (ML) point estimation of the parameter and the logarithm of the likelihood of the dataset. This paper demonstrates the information complexity (ICOMP), Bayesian deviance information, or the widely applicable information criterion (WAIC) of the BRMS to hierarchical linear models fitted with repeated measures with a simulation and two real data examples. The Fisher information matrix for the Bayesian hierarchical model considering fixed and random parameters under maximizing a posterior estimation is derived. Using Gibbs sampling and Hybrid Hamiltonian Monte Carlo approaches, six different models were fitted for three distinct application datasets. The best-fitted candidate models were identified under each application dataset with the two MCMC approaches. In this case, the Bayesian hierarchical (mixed effect) linear model with random intercepts and random slopes estimated using the Hamiltonian Monte Carlo method best fits the two application datasets. Information complexity (ICOMP) is a better indicator of the best-fitted models than DIC and WAIC. In addition, the information complexity criterion showed that hierarchical models with gradient-based Hamiltonian Monte Carlo estimation are the best fit and have supper convergence relative to the gradient-free Gibbs sampling methods.

Incorporating Road Network Connectivity in Neighboring Structures for Crash Prediction Models at the Area Level

Spatial correlation models have been traditionally used in road safety to account for spatial effects resulting from unmeasured or unknown risk factors that induce spatial correlation between neighboring areas. In transportation, the interaction between neighboring areas is highly influenced by the number of roads that connect those areas and the importance of those roads. This paper proposes an approach in which the weights of the spatial interaction (and therefore the spatial correlation) between areas depends on the number of road connections between those areas and the importance of those connections. The results using districts in Costa Rica show that the inclusion of road network connectivity in the models of spatial correlation significantly improves model fit, even after accounting for model complexity using the deviance information criterion (DIC) and widely applicable information criterion (WAIC). The inclusion of higher weights for national roads compared to municipal or local roads further improved the model fit. The best three models with respect to the posterior deviance, DIC, and WAIC are those that give at least three times more weight to national roads compared to local roads. With respect to site ranking, those three models present similar results, which also highlights the consistency among those models.

Automatic detection of myocardial ischaemia using generalisable spatio-temporal hierarchical Bayesian modelling of DCE-MRI

Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI) can be used as a non-invasive method for the assessment of myocardial perfusion. The acquired images can be utilised to analyse the spatial extent and severity of myocardial ischaemia (regions with impaired microvascular blood flow). In the present paper, we propose a novel generalisable spatio-temporal hierarchical Bayesian model (GST-HBM) to automate the detection of ischaemic lesions and improve the in silico prediction accuracy by systematically integrating spatio-temporal context information. We present a computational inference procedure with an adequate trade-off between accuracy and computational efficiency, whereby model parameters are sampled from the posterior distribution with Gibbs sampling, while lower-level hyperparameters are selected using model selection strategies based on the Watanabe Akaike information criterion (WAIC). We have assessed our method on both synthetic (in silico) data with known gold-standard and 12 sets of clinical first-pass myocardial perfusion DCE-MRI datasets. We have also carried out a comparative performance evaluation with four established alternative methods: Gaussian mixture model (GMM), opening and closing operations based on Gaussian mixture model (GMMC&Omax), Markov random field constrained Gaussian mixture model (GMM-MRF) and model-based hierarchical Bayesian model (M-HBM). Our results show that the proposed GST-HBM method achieves much higher in silico prediction accuracy than the established alternative methods. Furthermore, this method appears to provide a more robust delineation of ischaemic lesions in datasets affected by spatially variant noise.

Prediction for Risk of Hypouricemia in Hospitalized Patients: A Single-center, Retrospective Cohort Study.

Background Hypouricemia, defined as a serum uric acid (SUA) level ≤2 mg/dL, could be a risk factor for death in hospitalized patients. However, how explanatory variables can explain hypouricemia as an objective variable in a logistic regression analysis remains unknown. Purpose To predict the risk factors for hypouricemia in hospitalized patients using a robust Bayesian logistic (RBL) model. Methods This study retrospectively enrolled patients who visited Yonago Medical Center between April 2020 and March 2021. The association between potential risk factors and hypouricemia was analyzed using the RBL model in Python-modulated PyMC3. The final model was selected based on the lowest Watanabe-Akaike information criterion (WAIC). Results Of the 618 patients, 64 (10.4%) had hypouricemia. Based on the model according to the lowest WAIC, independent risk factors for hypouricemia were febuxostat [odds ratio (OR) 5.46, 95% confidence interval (CI) 2.32-13.4], amino acids in parenteral nutrition (OR 5.19, 95% CI 1.62-15.1), TMP-SMX (OR 4.20, 95% CI 1.66-10.9), emaciation (OR 3.48, 95% CI 1.75-7.21), and serum sodium level (OR 0.90, 95% CI 0.84-0.96). Conclusion The RBL model predicted amino acids in parenteral nutrition, TMP-SMX, emaciation, and low serum sodium levels for hypouricemia, in addition to the authentic risk factor febuxostat.

Comparing DIC and WAIC for multilevel models with missing data.

In Bayesian statistics, the most widely used criteria of Bayesian model assessment and comparison are Deviance Information Criterion (DIC) and Watanabe-Akaike Information Criterion (WAIC). We use a multilevel mediation model as an illustrative example to compare different types of DIC and WAIC. More specifically, we aim to compare the performance of conditional and marginal DICs and WAICs, and investigate their performance with missing data. We focus on two versions of DIC ( and ) and one version of WAIC. In addition, we explore whether it is necessary to include the nuisance models of incomplete exogenous variables in likelihood. Based on the simulation results, whether is better than and WAIC and whether we should include the nuisance models of exogenous variables in likelihood functions depend on whether we use marginal or conditional likelihoods. Overall, we find that the marginal likelihood based- that excludes the likelihood of covariate models generally had the highest true model selection rates.

Importance sampling and Bayesian model comparison in ecology and evolution

Abstract Bayesian approaches to the modelling of ecological systems are increasingly popular, but there are competing methods for formal model comparisons. Here, we focus on the task of performing multimodel inference through estimating posterior model weights, which encompasses uncertainties in the choice of competing model structure into the inference outputs. Model‐based approaches such as reversible‐jump Markov chain Monte Carlo (RJ‐MCMC) are flexible and allow multimodel inference, but can be complex to implement and optimise, and so we translate a model‐based approach for ecological applications using Importance Sampling to estimate the marginal likelihood of the data given a particular model. This approach allows for model comparison through the estimation of Bayes' Factors or interpretable posterior model probabilities, yielding model weights that facilitate multimodel inference through Bayesian model averaging. We demonstrate Importance Sampling with two case study investigations in animal demography: censused analysis of banded mongoose (Mungos mungo) survival where missing data are uncommon, and capture–mark–recapture analysis of European badger (Meles meles) survival where data are commonly missing. We compare outcomes of the model comparison using the Importance Sampling approach to those obtained through single‐model inference approaches using Deviance information criteria and the Watanabe–Akaike information criteria. The results of the Importance Sampling method aligns with RJ‐MCMC model comparisons while often being more straightforward to fit and optimise, particularly if the competing models are non‐nested.

Are ride-hailing services safer than taxis? A multivariate spatial approach with accommodation of exposure uncertainty

Despite many research efforts on ride-hailing services and taxis, limited studies have compared the safety performance of the two modes. A major challenge is the need for reliable mode-specific exposure data to model their safety outcomes. Moreover, crash frequencies of the two modes by injury severities tend to be spatially and inherently correlated. To fully address these issues, this study proposes a novel multivariate conditional autoregressive model considering measurement errors in mode-specific exposures (MVCARME). More specially, a classical measurement error structure accommodates the uncertainty of estimated mode-specific exposures, and a multivariate spatial specification is adopted to capture potential spatial and inherent correlations. The model estimation is accelerated by an integrated nest Laplace approximation method. The census tracts in the city of Chicago are set as the spatial analysis unit. The mode-specific exposures (vehicle-mile-traveled) in each census tract are estimated by trip assignments using ride-hailing and taxi trip data in 2019. The modeling results indicate that both ride-hailing crashes and taxi crashes are positively associated with transportation factors (e.g., vehicle-mile-traveled, mode-specific vehicle-mile-traveled, and traffic signal numbers), land use factors (i.e., number of educational and alcohol-related sites), and demographic factors (e.g., median household income, transit ratio, and walk ratio). By comparison, the proposed model outperforms the others (i.e., negative binomial models and multivariate conditional autoregressive model) by yielding the lowest deviance information criterion (DIC), Watanabe-Akaike information criterion (WAIC), mean absolute error (MAE), and root-mean-square error (RMSE). According to the results of t-tests, ride-hailing services are found to be prone to a higher risk of minor injury crashes compared with taxis, despite no significant difference between the risks of severe injury crashes. Methodologically, this study adds a robust safety evaluation approach for comparing crash risks of different modes to the literature. At the same time, practically, it provides researchers, practitioners, and policy-makers insights into the safety management of various mobility alternatives.