Search

Background: It is difficult to detect the outbreak of emergency infectious disease based on the exiting surveillance system. Here we investigate the utility of the Baidu Search Index, an indicator of how large of a keyword is in Baidu's search volume, in the early warning and predicting the epidemic trend of COVID-19., Methods: The daily number of cases and the Baidu Search Index of 8 keywords (weighted by population) from December 1, 2019 to March 15, 2020 were collected and analyzed with times series and Spearman correlation with different time lag. To predict the daily number of COVID-19 cases using the Baidu Search Index, Zero-inflated negative binomial regression was used in phase 1 and negative binomial regression model was used in phase 2 and phase 3 based on the characteristic of independent variable., Results: The Baidu Search Index of all keywords in Wuhan was significantly higher than Hubei (excluded Wuhan) and China (excluded Hubei). Before the causative pathogen was identified, the search volume of "Influenza" and "Pneumonia" in Wuhan increased with the number of new onset cases, their correlation coefficient was 0.69 and 0.59, respectively. After the pathogen was public but before COVID-19 was classified as a notifiable disease, the search volume of "SARS", "Pneumonia", "Coronavirus" in all study areas increased with the number of new onset cases with the correlation coefficient was 0.69 ~ 0.89, while "Influenza" changed to negative correlated (r s : -0.56 ~ -0.64). After COVID-19 was closely monitored, the Baidu Search Index of "COVID-19", "Pneumonia", "Coronavirus", "SARS" and "Mask" could predict the epidemic trend with 15 days, 5 days and 6 days lead time, respectively in Wuhan, Hubei (excluded Wuhan) and China (excluded Hubei). The predicted number of cases would increase 1.84 and 4.81 folds, respectively than the actual number of cases in Wuhan and Hubei (excluded Wuhan) from 21 January to 9 February., Conclusion: The Baidu Search Index could be used in the early warning and predicting the epidemic trend of COVID-19, but the search keywords changed in different period. Considering the time lag from onset to diagnosis, especially in the areas with medical resources shortage, internet search data can be a highly effective supplement of the existing surveillance system., (© 2024. The Author(s).)

Background: By analysing the deaths of inpatients in a tertiary hospital in Hangzhou, this study aimed to understand the epidemiological distribution characteristics and the composition of the causes of death. Additionally, this study aimed to predict the changing trend in the number of deaths, providing valuable insights for hospitals to formulate relevant strategies and measures aimed at reducing mortality rates., Methods: In this study, data on inpatient mortality at a tertiary hospital in Hangzhou from 2015 to 2022 were obtained via the population information registration system of the Chinese Center for Disease Control and Prevention. The death data of inpatients were described and analysed through a retrospective study. Excel 2016 was utilized for data sorting, and SPSS 22.0 software was employed for data analysis. The statistical inference of single factor differences was conducted via χ2 tests. The SARIMA model was established via the forecast, aTSA, and tseries software packages (version 4.3.0) to forecast future changes in the number of deaths., Results: A total of 1938 inpatients died at the tertiary hospital in Hangzhou, with the greatest number of deaths occurring in 2022 (262, 13.52%). The sex ratio was 2.22:1, and there were significant differences between sexes in terms of age, marital status, educational level, and place of residence (P < 0.05). The percentage of males in the groups aged of 20 to 29 and 30 to 39 years was significantly greater than that of females (χ 2  = 46.905, P < 0.001). More females than males died in the widowed group, and divorced and married males experienced a greater number of deaths than divorced and married females did (χ 2  = 61.130, P < 0.001). The proportions of male students with a junior college and senior high school education were significantly greater than that of female students (χ 2  = 12.310, P < 0.05). The primary causes of mortality within the hospital setting included circulatory system diseases, injury, poisoning, tumours, and respiratory system diseases. These leading factors accounted for 86.12% of all recorded deaths. Finally, the SARIMA (2, 1, 1) (1, 1, 1) 12 model was determined to be the optimal model, with an AIC of 380.23, a BIC of 392.79, and an AICc of 381.81. The MAPE was 14.99%, indicating a satisfactory overall fit of this model. The relative error between the predicted and actual number of deaths in 2022 was 8.02%. Therefore, the SARIMA (2, 1, 1) (1, 1, 1) 12 model demonstrates good predictive performance., Conclusions: Hospitals should enhance the management of sudden cardiac death, acute myocardial infarction, severe craniocerebral injury, lung cancer, and lung infection to reduce the mortality rate. The SARIMA model can be employed for predicting the number of deaths., (© 2024. The Author(s).)

Background: Survey studies in medical and health sciences predominantly apply a conventional direct questioning (DQ) format to gather private and highly personal information. If the topic under investigation is sensitive or even stigmatizing, such as COVID-19-related health behaviors and adherence to non-pharmaceutical interventions in general, DQ surveys can lead to nonresponse and untruthful answers due to the influence of social desirability bias (SDB). These effects seriously threaten the validity of the results obtained, potentially leading to distorted prevalence estimates for behaviors for which the prevalence in the population is unknown. While this issue cannot be completely avoided, indirect questioning techniques (IQTs) offer a means to mitigate the harmful influence of SDB by guaranteeing the confidentiality of individual responses. The present study aims at assessing the validity of a recently proposed IQT, the Cheating Detection Triangular Model (CDTRM), in estimating the prevalence of COVID-19-related health behaviors while accounting for cheaters who disregard the instructions., Methods: In an online survey of 1,714 participants in Taiwan, we obtained CDTRM prevalence estimates via an Expectation-Maximization algorithm for three COVID-19-related health behaviors with different levels of sensitivity. The CDTRM estimates were compared to DQ estimates and to available official statistics provided by the Taiwan Centers for Disease Control. Additionally, the CDTRM allowed us to estimate the share of cheaters who disregarded the instructions and adjust the prevalence estimates for the COVID-19-related health behaviors accordingly., Results: For a behavior with low sensitivity, CDTRM and DQ estimates were expectedly comparable and in line with official statistics. However, for behaviors with medium and high sensitivity, CDTRM estimates were higher and thus presumably more valid than DQ estimates. Analogously, the estimated cheating rate increased with higher sensitivity of the behavior under study., Conclusions: Our findings strongly support the assumption that the CDTRM successfully controlled for the validity-threatening influence of SDB in a survey on three COVID-19-related health behaviors. Consequently, the CDTRM appears to be a promising technique to increase estimation validity compared to conventional DQ for health-related behaviors, and sensitive attributes in general, for which a strong influence of SDB is to be expected., (© 2024. The Author(s).)

Background: Experimental studies of wound healing often use survival analysis and time to event outcomes or differences in wound area at a specific time point. However, these methods do not use a potentially large number of observations made over the course of a trial and may be inefficient. A model-based approach can leverage all trial data, but there is little guidance on appropriate models and functional forms to describe wound healing., Methods: We derive a general statistical model and review a wide range of plausible mathematical models to describe wound healing. We identify a range of possible derived estimands and their derivation from the models. Using data from a trial of an intervention to promote ulcer healing in patients affected by leprosy that included three measurement methods repeated across the course of the study, we compare the goodness-of-fit of the models using a range of methods and estimate treatment effects and healing rate functions with the best-fitting models., Results: Overall, we included 5,581 ulcer measurements of 1,578 unique images from 130 patients. We examined the performance of a range of models. The square root, log square root, and log quadratic models were the best fitting models across all outcome measurement methods. The estimated treatment effects magnitude and sign varied by time post-randomisation, model type, and outcome type, but across all models there was little evidence of effectiveness. The estimated effects were significantly more precise than non-parametric alternatives. For example, estimated differences from the three outcome measurements at 42-days post-randomisation were - 0.01 cm 2 (-0.77, 0.74), -0.44 cm 2 (-1.64, 0.76), and 0.11 cm 2 (-0.87, 1.08) using a non-parametric method versus - 0.03 cm 2 (-0.14, 0.06), 0.06 cm 2 (-0.05, 0.17), and 0.03 cm 2 (-0.07, 0.17) using a square-root model., Conclusions: Model-based analyses can dramatically improve the precision of estimates but care must be taken to carefully compare and select the best fitting models. The (log) square-root model is strongly recommended reflecting advice from a century ago., (© 2024. The Author(s).)

Background: The aim of this study is to analyze the trend of acute onset of chronic cor pulmonale at Chenggong Hospital of Kunming Yan'an Hospital between January 2018 and December 2022.Additionally, the study will compare the application of the ARIMA model and Holt-Winters model in predicting the number of chronic cor pulmonale cases., Methods: The data on chronic cor pulmonale cases from 2018 to 2022 were collected from the electronic medical records system of Chenggong Hospital of Kunming Yan'an Hospital. The ARIMA and Holt-Winters models were constructed using monthly case numbers from January 2018 to December 2022 as training data. The performance of the model was tested using the monthly number of cases from January 2023 to December 2023 as the test set., Results: The number of acute onset of chronic cor pulmonale in Chenggong Hospital of Kunming Yan'an Hospital exhibited a downward trend overall from 2018 to 2022. There were more cases in winter and spring, with peaks observed in November to December and January of the following year. The optimal ARIMA model was determined to be ARIMA (0,1,1) (0,1,1) 12 , while for the Holt-Winters model, the optimal choice was the Holt-Winters multiplicative model. It was found that the Holt-Winters multiplicative model yielded the lowest error., Conclusion: The Holt-Winters multiplicative model predicts better accuracy. The diagnosis of acute onset of chronic cor pulmonale is related to many risk factors, therefore, when using temporal models to fit and predict the data, we must consider such factors' influence and try to incorporate them into the models., (© 2024. The Author(s).)

Background: Multi-Arm, Multi-Stage (MAMS) clinical trial designs allow for multiple therapies to be compared across a spectrum of clinical trial phases. MAMS designs fall under several overarching design groups, including adaptive designs (AD) and multi-arm (MA) designs. Factorial clinical trials designs represent a combination of factorial and MAMS trial designs and can provide increased efficiency relative to fixed, traditional designs. We explore design choices associated with Factorial Adaptive Multi-Arm Multi-Stage (FAST) designs, which represent the combination of factorial and MAMS designs., Methods: Simulation studies were conducted to assess the impact of the type of analyses, the timing of analyses, and the effect size observed across multiple outcomes on trial operating characteristics for a FAST design. Given multiple outcomes types assessed within the hypothetical trial, the primary analysis approach for each assessment varied depending on data type., Results: The simulation studies demonstrate that the proposed class of FAST trial designs can offer a framework to potentially provide improvements relative to other trial designs, such as a MAMS or factorial trial. Further, we note that the design implementation decisions, such as the timing and type of analyses conducted throughout trial, can have a great impact on trial operating characteristics., Conclusions: Motivated by a trial currently under design, our work shows that the FAST category of trial can potentially offer benefits similar to both MAMS and factorial designs; however, the chosen design aspects which can be included in a FAST trial need to be thoroughly explored during the planning phase., (© 2024. The Author(s).)

Background: Chemical bioproduction has attracted attention as a key technology in a decarbonized society. In computational design for chemical bioproduction, it is necessary to predict changes in metabolic fluxes when up-/down-regulating enzymatic reactions, that is, responses of the system to enzyme perturbations. Structural sensitivity analysis (SSA) was previously developed as a method to predict qualitative responses to enzyme perturbations on the basis of the structural information of the reaction network. However, the network structural information can sometimes be insufficient to predict qualitative responses unambiguously, which is a practical issue in bioproduction applications. To address this, in this study, we propose BayesianSSA, a Bayesian statistical model based on SSA. BayesianSSA extracts environmental information from perturbation datasets collected in environments of interest and integrates it into SSA predictions., Results: We applied BayesianSSA to synthetic and real datasets of the central metabolic pathway of Escherichia coli. Our result demonstrates that BayesianSSA can successfully integrate environmental information extracted from perturbation data into SSA predictions. In addition, the posterior distribution estimated by BayesianSSA can be associated with the known pathway reported to enhance succinate export flux in previous studies., Conclusions: We believe that BayesianSSA will accelerate the chemical bioproduction process and contribute to advancements in the field., (© 2024. The Author(s).)

Background: In sub-Saharan African countries, preventable and manageable diseases such as diarrhea and acute respiratory infections still claim the lives of children. Hence, this study aims to estimate the rate of change in the log expected number of days a child suffers from Diarrhea (NOD) and flu/common cold (NOF) among children aged 6 to 11 months at the baseline of the study., Methodology: This study used secondary data which exhibit a longitudinal and multilevel structure. Based on the results of exploratory analysis, a multilevel zero-inflated Poisson regression model with a rate of change in the log expected NOD and NOF described by a quadratic trend was proposed to efficiently analyze both outcomes accounting for correlation between observations and individuals through random effects. Furthermore, residual plots were used to assess the goodness of fit of the model., Results: Considering subject and cluster-specific random effects, the results revealed a quadratic trend in the rate of change of the log expected NOD. Initially, low dose iron Micronutrient Powder (MNP) users exhibited a higher rate of change compared to non-users, but this trend reversed over time. Similarly, the log expected NOF decreased for children who used MNP and exclusively breastfed for six months, in comparison to their counterparts. In addition, the odds of not having flu decreased with each two-week increment for MNP users, as compared to non-MNP users. Furthermore, an increase in NOD resulted in an increase in the log expected NOF. Region and exclusive breastfeeding also have a significant relationships with both NOD and NOF., Conclusion: The findings of this study underscore the importance of commencing analysis of data generated from a study with exploratory analysis. The study highlights the critical role of promoting EBF for the first six months and supporting children with additional food after six months to reduce the burden of infectious diseases., (© 2024. The Author(s).)

The inability to correctly account for unmeasured confounding can lead to bias in parameter estimates, invalid uncertainty assessments, and erroneous conclusions. Sensitivity analysis is an approach to investigate the impact of unmeasured confounding in observational studies. However, the adoption of this approach has been slow given the lack of accessible software. An extensive review of available R packages to account for unmeasured confounding list deterministic sensitivity analysis methods, but no R packages were listed for probabilistic sensitivity analysis. The R package unmconf implements the first available package for probabilistic sensitivity analysis through a Bayesian unmeasured confounding model. The package allows for normal, binary, Poisson, or gamma responses, accounting for one or two unmeasured confounders from the normal or binomial distribution. The goal of unmconf is to implement a user friendly package that performs Bayesian modeling in the presence of unmeasured confounders, with simple commands on the front end while performing more intensive computation on the back end. We investigate the applicability of this package through novel simulation studies. The results indicate that credible intervals will have near nominal coverage probability and smaller bias when modeling the unmeasured confounder(s) for varying levels of internal/external validation data across various combinations of response-unmeasured confounder distributional families., (© 2024. The Author(s).)

Background: Early identification of children at high risk of developing myopia is essential to prevent myopia progression by introducing timely interventions. However, missing data and measurement error (ME) are common challenges in risk prediction modelling that can introduce bias in myopia prediction., Methods: We explore four imputation methods to address missing data and ME: single imputation (SI), multiple imputation under missing at random (MI-MAR), multiple imputation with calibration procedure (MI-ME), and multiple imputation under missing not at random (MI-MNAR). We compare four machine-learning models (Decision Tree, Naive Bayes, Random Forest, and Xgboost) and three statistical models (logistic regression, stepwise logistic regression, and least absolute shrinkage and selection operator logistic regression) in myopia risk prediction. We apply these models to the Shanghai Jinshan Myopia Cohort Study and also conduct a simulation study to investigate the impact of missing mechanisms, the degree of ME, and the importance of predictors on model performance. Model performance is evaluated using the receiver operating characteristic curve (AUROC) and the area under the precision-recall curve (AUPRC)., Results: Our findings indicate that in scenarios with missing data and ME, using MI-ME in combination with logistic regression yields the best prediction results. In scenarios without ME, employing MI-MAR to handle missing data outperforms SI regardless of the missing mechanisms. When ME has a greater impact on prediction than missing data, the relative advantage of MI-MAR diminishes, and MI-ME becomes more superior. Furthermore, our results demonstrate that statistical models exhibit better prediction performance than machine-learning models., Conclusion: MI-ME emerges as a reliable method for handling missing data and ME in important predictors for early-onset myopia risk prediction., (© 2024. The Author(s).)

Background: Chronic kidney disease (CKD) is an important public health problem worldwide; therefore, forecasting CKD mortality rates and death numbers globally is vital for planning CKD prevention programs. This study aimed to characterize the temporal trends in CKD mortality at the international level from 1990 to 2019 and predict CKD mortality rates and numbers until 2030., Methods: Data were obtained from the Global Burden of Disease 2019 Study. A joinpoint regression model was used to estimate the average annual percentage change in CKD mortality rates and numbers. Finally, we used a generalized additive model to predict CKD mortality through 2030., Results: The number of CKD-related deaths worldwide increased from 591.80 thousand in 1990 to 1425.67 thousand in 2019. The CKD age-adjusted mortality rate increased from 15.95 per 100,000 people to 18.35 per 100,000 people during the same period. Between 2020 and 2030, the number of CKD deaths is forecasted to increase further to 1812.85 thousand by 2030. The CKD age-adjusted mortality rate is expected to decrease slightly to 17.76 per 100,000 people (95% credible interval (CrI): 13.84 to 21.68). Globally, it is predicted that in the next decade, the CKD mortality rate will decrease in men, women, all subgroups of disease etiology except glomerulonephritis, people younger than 40 years old, and all groupings of countries based on the sociodemographic index (SDI) except high-middle-SDI countries., Conclusions: The CKD mortality rate is predicted to decrease in the next decade. However, more attention should be given to people with glomerulonephritis, people over 40 years old, and people in high- to middle-income countries because the mortality rate due to CKD in these subgroups is expected to increase until 2030., (© 2024. The Author(s).)

Background: Distributed statistical analyses provide a promising approach for privacy protection when analyzing data distributed over several databases. Instead of directly operating on data, the analyst receives anonymous summary statistics, which are combined into an aggregated result. Further, in discrimination model (prognosis, diagnosis, etc.) development, it is key to evaluate a trained model w.r.t. to its prognostic or predictive performance on new independent data. For binary classification, quantifying discrimination uses the receiver operating characteristics (ROC) and its area under the curve (AUC) as aggregation measure. We are interested to calculate both as well as basic indicators of calibration-in-the-large for a binary classification task using a distributed and privacy-preserving approach., Methods: We employ DataSHIELD as the technology to carry out distributed analyses, and we use a newly developed algorithm to validate the prediction score by conducting distributed and privacy-preserving ROC analysis. Calibration curves are constructed from mean values over sites. The determination of ROC and its AUC is based on a generalized linear model (GLM) approximation of the true ROC curve, the ROC-GLM, as well as on ideas of differential privacy (DP). DP adds noise (quantified by the ℓ 2 sensitivity Δ 2 ( f ^ ) ) to the data and enables a global handling of placement numbers. The impact of DP parameters was studied by simulations., Results: In our simulation scenario, the true and distributed AUC measures differ by Δ AUC < 0.01 depending heavily on the choice of the differential privacy parameters. It is recommended to check the accuracy of the distributed AUC estimator in specific simulation scenarios along with a reasonable choice of DP parameters. Here, the accuracy of the distributed AUC estimator may be impaired by too much artificial noise added from DP., Conclusions: The applicability of our algorithms depends on the ℓ 2 sensitivity Δ 2 ( f ^ ) of the underlying statistical/predictive model. The simulations carried out have shown that the approximation error is acceptable for the majority of simulated cases. For models with high Δ 2 ( f ^ ) , the privacy parameters must be set accordingly higher to ensure sufficient privacy protection, which affects the approximation error. This work shows that complex measures, as the AUC, are applicable for validation in distributed setups while preserving an individual's privacy., (© 2024. The Author(s).)

Introduction: While there is an interest in defining longitudinal change in people with chronic illness like Parkinson's disease (PD), statistical analysis of longitudinal data is not straightforward for clinical researchers. Here, we aim to demonstrate how the choice of statistical method may influence research outcomes, (e.g., progression in apathy), specifically the size of longitudinal effect estimates, in a cohort., Methods: In this retrospective longitudinal analysis of 802 people with typical Parkinson's disease in the Luxembourg Parkinson's study, we compared the mean apathy scores at visit 1 and visit 8 by means of the paired two-sided t-test. Additionally, we analysed the relationship between the visit numbers and the apathy score using linear regression and longitudinal two-level mixed effects models., Results: Mixed effects models were the only method able to detect progression of apathy over time. While the effects estimated for the group comparison and the linear regression were smaller with high p-values (+ 1.016/ 7 years, p = 0.107, -0.056/ 7 years, p = 0.897, respectively), effect estimates for the mixed effects models were positive with a very small p-value, indicating a significant increase in apathy symptoms by + 2.345/ 7 years (p < 0.001)., Conclusion: The inappropriate use of paired t-tests and linear regression to analyse longitudinal data can lead to underpowered analyses and an underestimation of longitudinal change. While mixed effects models are not without limitations and need to be altered to model the time sequence between the exposure and the outcome, they are worth considering for longitudinal data analyses. In case this is not possible, limitations of the analytical approach need to be discussed and taken into account in the interpretation., (© 2024. The Author(s).)

Background: Rifampicin resistant tuberculosis (RR-TB) poses a growing threat to individuals and communities. This study utilized a seasonal autoregressive integrated moving average (SARIMA) model to quantitatively predict the monthly incidence of RR-TB in Yunnan Province which could guide government health administration departments and the centers for disease control and prevention (CDC) in preventing and controlling the RR-TB epidemic., Methods: The study utilized routine surveillance reporting data from the infectious Disease Network Surveillance and Reporting System. Monthly incidence rates of RR-TB were collected from January 2019 to December 2022. A time series SARIMA model was used to predict the number of monthly RR-TB cases in Yunnan Province in 2023, and the model was validated using time series plots, seasonal and non-seasonal differencing, autocorrelation and partial autocorrelation analysis, and white noise tests., Results: From 2019 to 2022, the incidence of RR-TB decreases as the incidence of all TB decreases (P < 0.05). There was no significant change in the proportion of RR-TB among all TB cases, which remained within 2.5% (P>0.05). The time series decomposition shows that it presented obvious seasonality, periodicity and randomness after being decomposed. Time series analysis was performed on the original series after 1 non-seasonal difference and 1 seasonal difference, the ADF test showed P < 0.05. According to ACF and PACF, the SARIMA (1, 1, 1) (1, 1, 0) 12 model was chosen and statistically significant model parameter estimates (P < 0.05). The predicted seasonal trend of RR-TB incidence in 2019 to 2023 was similar to the actual data. The percentage accuracy in the prediction excesses 80% in 2019 to 2022 and is all within 95% CI. However there was a certain gap between the actual incidence and the predicted value in 2023, and the acutual incidence had increased by 12.4% compared to 2022. The percentage of accuracy in the prediction was only 70% in 2023., Conclusions: We found the incidence of RR-TB was based on that of all TB in Yunnan. The SARIMA model successfully predicted the seasonal incidence trend of RR-TB in Yunnan Province in 2019 to 2023, but the prediction precision could be influenced by factors such as new infectious disease outbreaks or pandemics, social issues, environmental challenges or other unknown risks. Hence CDCs should pay special attention to the post epidemic effects of new infectious disease outbreaks or pandemics, carry out monitoring and early warning, and better optimize disease prediction models., (© 2024. The Author(s).)

Background: Describing the transmission dynamics of infectious diseases across different regions is crucial for effective disease surveillance. The multivariate time series (MTS) model has been widely adopted for constructing cross-regional infectious disease transmission networks due to its strengths in interpretability and predictive performance. Nevertheless, the assumption of constant parameters frequently disregards the dynamic shifts in disease transmission rates, thereby compromising the accuracy of early warnings. This study investigated the applicability of time-varying MTS models in multi-regional infectious disease monitoring and explored strategies for model selection., Methods: This study focused on two prominent time-varying MTS models: the time-varying parameter-stochastic volatility-vector autoregression (TVP-SV-VAR) model and the time-varying VAR model using the generalized additive framework (tvvarGAM), and intended to explore and verify their applicable conditions for the surveillance of infectious diseases. For the first time, this study proposed the time delay coefficient and spatial sparsity indicators for model selection. These indicators quantify the temporal lags and spatial distribution of infectious disease data, respectively. Simulation study adopted from real-world infectious disease surveillance was carried out to compare model performances under various scenarios of spatio-temporal variation as well as random volatility. Meanwhile, we illustrated how the modelling process could help the surveillance of infectious diseases with an application to the influenza-like case in Sichuan Province, China., Results: When the spatio-temporal variation was small (time delay coefficient: 0.1-0.2, spatial sparsity:0.1-0.3), the TVP-SV-VAR model was superior with smaller fitting residuals and standard errors of parameter estimation than those of the tvvarGAM model. In contrast, the tvvarGAM model was preferable when the spatio-temporal variation increased (time delay coefficient: 0.2-0.3, spatial sparsity: 0.6-0.9)., Conclusion: This study emphasized the importance of considering spatio-temporal variations when selecting appropriate models for infectious disease surveillance. By incorporating our novel indicators-the time delay coefficient and spatial sparsity-into the model selection process, the study could enhance the accuracy and effectiveness of infectious disease monitoring efforts. This approach was not only valuable in the context of this study, but also has broader implications for improving time-varying MTS analyses in various applications., (© 2024. The Author(s).)

Background: Simulation is an important tool for assessing the performance of statistical methods for the analysis of data and for the planning of studies. While methods are available for the simulation of correlated binary random variables, all have significant practical limitations for simulating outcomes from longitudinal cluster randomised trial designs, such as the cluster randomised crossover and the stepped wedge trial designs. For these trial designs as the number of observations in each cluster increases these methods either become computationally infeasible or their range of allowable correlations rapidly shrinks to zero., Methods: In this paper we present a simple method for simulating binary random variables with a specified vector of prevalences and correlation matrix. This method allows for the outcome prevalence to change due to treatment or over time, and for a 'nested exchangeable' correlation structure, in which observations in the same cluster are more highly correlated if they are measured in the same time period than in different time periods, and where different individuals are measured in each time period. This means that our method is also applicable to more general hierarchical clustered data contexts, such as students within classrooms within schools. The method is demonstrated by simulating 1000 datasets with parameters matching those derived from data from a cluster randomised crossover trial assessing two variants of stress ulcer prophylaxis., Results: Our method is orders of magnitude faster than the most well known general simulation method while also allowing a much wider range of correlations than alternative methods. An implementation of our method is available in an R package NestBin., Conclusions: This simulation method is the first to allow for practical and efficient simulation of large datasets of binary outcomes with the commonly used nested exchangeable correlation structure. This will allow for much more effective testing of designs and inference methods for longitudinal cluster randomised trials with binary outcomes., (© 2024. The Author(s).)

Statistical regression models are used for predicting outcomes based on the values of some predictor variables or for describing the association of an outcome with predictors. With a data set at hand, a regression model can be easily fit with standard software packages. This bears the risk that data analysts may rush to perform sophisticated analyses without sufficient knowledge of basic properties, associations in and errors of their data, leading to wrong interpretation and presentation of the modeling results that lacks clarity. Ignorance about special features of the data such as redundancies or particular distributions may even invalidate the chosen analysis strategy. Initial data analysis (IDA) is prerequisite to regression analyses as it provides knowledge about the data needed to confirm the appropriateness of or to refine a chosen model building strategy, to interpret the modeling results correctly, and to guide the presentation of modeling results. In order to facilitate reproducibility, IDA needs to be preplanned, an IDA plan should be included in the general statistical analysis plan of a research project, and results should be well documented. Biased statistical inference of the final regression model can be minimized if IDA abstains from evaluating associations of outcome and predictors, a key principle of IDA. We give advice on which aspects to consider in an IDA plan for data screening in the context of regression modeling to supplement the statistical analysis plan. We illustrate this IDA plan for data screening in an example of a typical diagnostic modeling project and give recommendations for data visualizations., (© 2024. The Author(s).)

Background: Prediction models are often externally validated with data from a single study or cohort. However, the interpretation of performance estimates obtained with single-study external validation is not as straightforward as assumed. We aimed to illustrate this by conducting a large number of external validations of a prediction model for functional outcome in subarachnoid hemorrhage (SAH) patients., Methods: We used data from the Subarachnoid Hemorrhage International Trialists (SAHIT) data repository (n = 11,931, 14 studies) to refit the SAHIT model for predicting a dichotomous functional outcome (favorable versus unfavorable), with the (extended) Glasgow Outcome Scale or modified Rankin Scale score, at a minimum of three months after discharge. We performed leave-one-cluster-out cross-validation to mimic the process of multiple single-study external validations. Each study represented one cluster. In each of these validations, we assessed discrimination with Harrell's c-statistic and calibration with calibration plots, the intercepts, and the slopes. We used random effects meta-analysis to obtain the (reference) mean performance estimates and between-study heterogeneity (I 2 -statistic). The influence of case-mix variation on discriminative performance was assessed with the model-based c-statistic and we fitted a "membership model" to obtain a gross estimate of transportability., Results: Across 14 single-study external validations, model performance was highly variable. The mean c-statistic was 0.74 (95%CI 0.70-0.78, range 0.52-0.84, I 2  = 0.92), the mean intercept was -0.06 (95%CI -0.37-0.24, range -1.40-0.75, I 2  = 0.97), and the mean slope was 0.96 (95%CI 0.78-1.13, range 0.53-1.31, I 2  = 0.90). The decrease in discriminative performance was attributable to case-mix variation, between-study heterogeneity, or a combination of both. Incidentally, we observed poor generalizability or transportability of the model., Conclusions: We demonstrate two potential pitfalls in the interpretation of model performance with single-study external validation. With single-study external validation. (1) model performance is highly variable and depends on the choice of validation data and (2) no insight is provided into generalizability or transportability of the model that is needed to guide local implementation. As such, a single single-study external validation can easily be misinterpreted and lead to a false appreciation of the clinical prediction model. Cross-validation is better equipped to address these pitfalls., (© 2024. The Author(s).)

We have introduced the R package jmBIG to facilitate the analysis of large healthcare datasets and the development of predictive models. This package provides a comprehensive set of tools and functions specifically designed for the joint modelling of longitudinal and survival data in the context of big data analytics. The jmBIG package offers efficient and scalable implementations of joint modelling algorithms, allowing for integrating large-scale healthcare datasets.By utilizing the capabilities of jmBIG, researchers and analysts can effectively handle the challenges associated with big healthcare data, such as high dimensionality and complex relationships between multiple outcomes.With the support of jmBIG, analysts can seamlessly fit Bayesian joint models, generate predictions, and evaluate the performance of the models. The package incorporates cutting-edge methodologies and harnesses the computational capabilities of parallel computing to accelerate the analysis of large-scale healthcare datasets significantly. In summary, jmBIG empowers researchers to gain deeper insights into disease progression and treatment response, fostering evidence-based decision-making and paving the way for personalized healthcare interventions that can positively impact patient outcomes on a larger scale., (© 2024. The Author(s).)

Background: Mediation analysis, often completed as secondary analysis to estimating the main treatment effect, investigates situations where an exposure may affect an outcome both directly and indirectly through intervening mediator variables. Although there has been much research on power in mediation analyses, most of this has focused on the power to detect indirect effects. Little consideration has been given to the extent to which the strength of the mediation pathways, i.e., the intervention-mediator path and the mediator-outcome path respectively, may affect the power to detect the total effect, which would correspond to the intention-to-treat effect in a randomized trial., Methods: We conduct a simulation study to evaluate the relation between the mediation pathways and the power of testing the total treatment effect, i.e., the intention-to-treat effect. Consider a sample size that is computed based on the usual formula for testing the total effect in a two-arm trial. We generate data for a continuous mediator and a normal outcome using the conventional mediation models. We estimate the total effect using simple linear regression and evaluate the power of a two-sided test. We explore multiple data generating scenarios by varying the magnitude of the mediation paths whilst keeping the total effect constant., Results: Simulations show the estimated total effect is unbiased across the considered scenarios as expected, but the mean of its standard error increases with the magnitude of the mediator-outcome path and the variability in the residual error of the mediator, respectively. Consequently, this affects the power of testing the total effect, which is always lower than planned when the mediator-outcome path is non-trivial and a naive sample size was employed. Analytical explanation confirms that the intervention-mediator path does not affect the power of testing the total effect but the mediator-outcome path. The usual effect size consideration can be adjusted to account for the magnitude of the mediator-outcome path and its residual error., Conclusions: The sample size calculation for studies with efficacy and mechanism evaluation should account for the mediator-outcome association or risk the power to detect the total effect/intention-to-treat effect being lower than planned., (© 2024. The Author(s).)

Background: Dimension reduction methods do not always reduce their underlying indicators to a single composite score. Furthermore, such methods are usually based on optimality criteria that require discarding some information. We suggest, under some conditions, to use the joint probability density function (joint pdf or JPD) of p-dimensional random variable (the p indicators), as an index or a composite score. It is proved that this index is more informative than any alternative composite score. In two examples, we compare the JPD index with some alternatives constructed from traditional methods., Methods: We develop a probabilistic unsupervised dimension reduction method based on the probability density of multivariate data. We show that the conditional distribution of the variables given JPD is uniform, implying that the JPD is the most informative scalar summary under the most common notions of information. B. We show under some widely plausible conditions, JPD can be used as an index. To use JPD as an index, in addition to having a plausible interpretation, all the random variables should have approximately the same direction(unidirectionality) as the density values (codirectionality). We applied these ideas to two data sets: first, on the 7 Brief Pain Inventory Interference scale (BPI-I) items obtained from 8,889 US Veterans with chronic pain and, second, on a novel measure based on administrative data for 912 US Veterans. To estimate the JPD in both examples, among the available JPD estimation methods, we used its conditional specifications, identified a well-fitted parametric model for each factored conditional (regression) specification, and, by maximizing the corresponding likelihoods, estimated their parameters. Due to the non-uniqueness of conditional specification, the average of all estimated conditional specifications was used as the final estimate. Since a prevalent common use of indices is ranking, we used measures of monotone dependence [e.g., Spearman's rank correlation (rho)] to assess the strength of unidirectionality and co-directionality. Finally, we cross-validate the JPD score against variance-covariance-based scores (factor scores in unidimensional models), and the "person's parameter" estimates of (Generalized) Partial Credit and Graded Response IRT models. We used Pearson Divergence as a measure of information and Shannon's entropy to compare uncertainties (informativeness) in these alternative scores., Results: An unsupervised dimension reduction was developed based on the joint probability density (JPD) of the multi-dimensional data. The JPD, under regularity conditions, may be used as an index. For the well-established Brief Pain Interference Inventory (BPI-I (the short form with 7 Items) and for a new mental health severity index (MoPSI) with 6 indicators, we estimated the JPD scoring. We compared, assuming unidimensionality, factor scores, Person's scores of the Partial Credit model, the Generalized Partial Credit model, and the Graded Response model with JPD scoring. As expected, all scores' rankings in both examples were monotonically dependent with various strengths. Shannon entropy was the smallest for JPD scores. Pearson Divergence of the estimated densities of different indices against uniform distribution was maximum for JPD scoring., Conclusions: An unsupervised probabilistic dimension reduction is possible. When appropriate, the joint probability density function can be used as the most informative index. Model specification and estimation and steps to implement the scoring were demonstrated. As expected, when the required assumption in factor analysis and IRT models are satisfied, JPD scoring agrees with these established scores. However, when these assumptions are violated, JPD scores preserve all the information in the indicators with minimal assumption., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Background: Pocock-Simon's minimisation method has been widely used to balance treatment assignments across prognostic factors in randomised controlled trials (RCTs). Previous studies focusing on the survival outcomes have demonstrated that the conservativeness of asymptotic tests without adjusting for stratification factors, as well as the inflated type I error rate of adjusted asymptotic tests conducted in a small sample of patients, can be relaxed using re-randomisation tests. Although several RCTs using minimisation have suggested the presence of non-proportional hazards (non-PH) effects, the application of re-randomisation tests has been limited to the log-rank test and Cox PH models, which may result in diminished statistical power when confronted with non-PH scenarios. To address this issue, we proposed two re-randomisation tests based on a maximum combination of weighted log-rank tests (MaxCombo test) and the difference in restricted mean survival time (dRMST) up to a fixed time point τ , both of which can be extended to adjust for randomisation stratification factors., Methods: We compared the performance of asymptotic and re-randomisation tests using the MaxCombo test, dRMST, log-rank test, and Cox PH models, assuming various non-PH situations for RCTs using minimisation, with total sample sizes of 50, 100, and 500 at a 1:1 allocation ratio. We mainly considered null, and alternative scenarios featuring delayed, crossing, and diminishing treatment effects., Results: Across all examined null scenarios, re-randomisation tests maintained the type I error rates at the nominal level. Conversely, unadjusted asymptotic tests indicated excessive conservatism, while adjusted asymptotic tests in both the Cox PH models and dRMST indicated inflated type I error rates for total sample sizes of 50. The stratified MaxCombo-based re-randomisation test consistently exhibited robust power across all examined scenarios., Conclusions: The re-randomisation test is a useful alternative in non-PH situations for RCTs with minimisation using the stratified MaxCombo test, suggesting its robust power in various scenarios., (© 2024. The Author(s).)

Background: A platform trial approach allows adding arms to on-going trials to speed up intervention discovery programs. A control arm remains open for recruitment in a platform trial while intervention arms may be added after the onset of the study and could be terminated early for efficacy and/or futility when early stopping is allowed. The topic of utilising non-concurrent control data in the analysis of platform trials has been explored and discussed extensively. A less familiar issue is the presence of heterogeneity, which may exist for example due to modification of enrolment criteria and recruitment strategy., Method: We conduct a simulation study to explore the impact of heterogeneity on the analysis of a two-stage platform trial design. We consider heterogeneity in treatment effects and heteroscedasticity in outcome data across stages for a normally distributed endpoint. We examine the performance of some hypothesis testing procedures and modelling strategies. The use of non-concurrent control data is also considered accordingly. Alongside standard regression analysis, we examine the performance of a novel method that was known as the pairwise trials analysis. It is similar to a network meta-analysis approach but adjusts for treatment comparisons instead of individual studies using fixed effects., Results: Several testing strategies with concurrent control data seem to control the type I error rate at the required level when there is heteroscedasticity in outcome data across stages and/or a random cohort effect. The main parameter of treatment effects in some analysis models correspond to overall treatment effects weighted by stage wise sample sizes; while others correspond to the effect observed within a single stage. The characteristics of the estimates are not affected significantly by the presence of a random cohort effect and/ or heteroscedasticity., Conclusion: In view of heterogeneity in treatment effect across stages, the specification of null hypotheses in platform trials may need to be more subtle. We suggest employing testing procedure of adaptive design as opposed to testing the statistics from regression models; comparing the estimates from the pairwise trials analysis method and the regression model with interaction terms may indicate if heterogeneity is negligible., (© 2024. The Author(s).)

The Platform trial In COVID-19 priming and BOOsting (PICOBOO) is a multi-site, adaptive platform trial designed to generate evidence of the immunogenicity, reactogenicity, and cross-protection of different booster vaccination strategies against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants, specific for the Australian context. The PICOBOO trial randomises participants to receive one of three COVID-19 booster vaccine brands (Pfizer, Moderna, Novavax) available for use in Australia, where the vaccine brand subtypes vary over time according to the national vaccine roll out strategy, and employs a Bayesian hierarchical modelling approach to efficiently borrow information across consecutive booster doses, age groups and vaccine brand subtypes. Here, we briefly describe the PICOBOO trial structure and report the statistical considerations for the estimands, statistical models and decision making for trial adaptations. This paper should be read in conjunction with the PICOBOO Core Protocol and PICOBOO Sub-Study Protocol 1: Booster Vaccination. PICOBOO was registered on 10 February 2022 with the Australian and New Zealand Clinical Trials Registry ACTRN12622000238774., (© 2024. The Author(s).)

Systematic reviews and meta-analyses are essential tools in contemporary evidence-based medicine, synthesizing evidence from various sources to better inform clinical decision-making. However, the conclusions from different meta-analyses on the same topic can be discrepant, which has raised concerns about their reliability. One reason is that the result of a meta-analysis is sensitive to factors such as study inclusion/exclusion criteria and model assumptions. The arm-based meta-analysis model is growing in importance due to its advantage of including single-arm studies and historical controls with estimation efficiency and its flexibility in drawing conclusions with both marginal and conditional effect measures. Despite its benefits, the inference may heavily depend on the heterogeneity parameters that reflect design and model assumptions. This article aims to evaluate the robustness of meta-analyses using the arm-based model within a Bayesian framework. Specifically, we develop a tipping point analysis of the between-arm correlation parameter to assess the robustness of meta-analysis results. Additionally, we introduce some visualization tools to intuitively display its impact on meta-analysis results. We demonstrate the application of these tools in three real-world meta-analyses, one of which includes single-arm studies., (© 2024. The Author(s).)

Background: Electronic Health Records (EHR) are widely used to develop clinical prediction models (CPMs). However, one of the challenges is that there is often a degree of informative missing data. For example, laboratory measures are typically taken when a clinician is concerned that there is a need. When data are the so-called Not Missing at Random (NMAR), analytic strategies based on other missingness mechanisms are inappropriate. In this work, we seek to compare the impact of different strategies for handling missing data on CPMs performance., Methods: We considered a predictive model for rapid inpatient deterioration as an exemplar implementation. This model incorporated twelve laboratory measures with varying levels of missingness. Five labs had missingness rate levels around 50%, and the other seven had missingness levels around 90%. We included them based on the belief that their missingness status can be highly informational for the prediction. In our study, we explicitly compared the various missing data strategies: mean imputation, normal-value imputation, conditional imputation, categorical encoding, and missingness embeddings. Some of these were also combined with the last observation carried forward (LOCF). We implemented logistic LASSO regression, multilayer perceptron (MLP), and long short-term memory (LSTM) models as the downstream classifiers. We compared the AUROC of testing data and used bootstrapping to construct 95% confidence intervals., Results: We had 105,198 inpatient encounters, with 4.7% having experienced the deterioration outcome of interest. LSTM models generally outperformed other cross-sectional models, where embedding approaches and categorical encoding yielded the best results. For the cross-sectional models, normal-value imputation with LOCF generated the best results., Conclusion: Strategies that accounted for the possibility of NMAR missing data yielded better model performance than those did not. The embedding method had an advantage as it did not require prior clinical knowledge. Using LOCF could enhance the performance of cross-sectional models but have countereffects in LSTM models., (© 2024. The Author(s).)

Background: Society is burdened with stroke-associated pneumonia (SAP) after intracerebral haemorrhage (ICH). Cerebral small vessel disease (CSVD) complicates clinical manifestations of stroke. In this study, we redefined the CSVD burden score and incorporated it into a novel radiological-clinical prediction model for SAP., Materials and Methods: A total of 1278 patients admitted to a tertiary hospital between 1 January 2010 and 31 December 2019 were included. The participants were divided into training and testing groups using fivefold cross-validation method. Four models, two traditional statistical models (logistic regression and ISAN) and two machine learning models (random forest and support vector machine), were established and evaluated. The outcomes and baseline characteristics were compared between the SAP and non-SAP groups., Results: Among the of 1278 patients, 281(22.0%) developed SAP after their first ICH. Multivariate analysis revealed that the logistic regression (LR) model was superior in predicting SAP in both the training and testing groups. Independent predictors of SAP after ICH included total CSVD burden score (OR, 1.29; 95% CI, 1.03-1.54), haematoma extension into ventricle (OR, 2.28; 95% CI, 1.87-3.31), haematoma with multilobar involvement (OR, 2.14; 95% CI, 1.44-3.18), transpharyngeal intubation operation (OR, 3.89; 95% CI, 2.7-5.62), admission NIHSS score ≥ 10 (OR, 2.06; 95% CI, 1.42-3.01), male sex (OR, 1.69; 95% CI, 1.16-2.52), and age ≥ 67 (OR, 2.24; 95% CI, 1.56-3.22). The patients in the SAP group had worse outcomes than those in the non-SAP group., Conclusion: This study established a clinically combined imaging model for predicting stroke-associated pneumonia and demonstrated superior performance compared with the existing ISAN model. Given the poor outcomes observed in patients with SAP, the use of individualised predictive nomograms is vital in clinical practice., (© 2024. The Author(s).)

Background: In randomized clinical trials, treatment effects may vary, and this possibility is referred to as heterogeneity of treatment effect (HTE). One way to quantify HTE is to partition participants into subgroups based on individual's risk of experiencing an outcome, then measuring treatment effect by subgroup. Given the limited availability of externally validated outcome risk prediction models, internal models (created using the same dataset in which heterogeneity of treatment analyses also will be performed) are commonly developed for subgroup identification. We aim to compare different methods for generating internally developed outcome risk prediction models for subject partitioning in HTE analysis., Methods: Three approaches were selected for generating subgroups for the 2,441 participants from the United States enrolled in the ASPirin in Reducing Events in the Elderly (ASPREE) randomized controlled trial. An extant proportional hazards-based outcomes predictive risk model developed on the overall ASPREE cohort of 19,114 participants was identified and was used to partition United States' participants by risk of experiencing a composite outcome of death, dementia, or persistent physical disability. Next, two supervised non-parametric machine learning outcome classifiers, decision trees and random forests, were used to develop multivariable risk prediction models and partition participants into subgroups with varied risks of experiencing the composite outcome. Then, we assessed how the partitioning from the proportional hazard model compared to those generated by the machine learning models in an HTE analysis of the 5-year absolute risk reduction (ARR) and hazard ratio for aspirin vs. placebo in each subgroup. Cochran's Q test was used to detect if ARR varied significantly by subgroup., Results: The proportional hazard model was used to generate 5 subgroups using the quintiles of the estimated risk scores; the decision tree model was used to generate 6 subgroups (6 automatically determined tree leaves); and the random forest model was used to generate 5 subgroups using the quintiles of the prediction probability as risk scores. Using the semi-parametric proportional hazards model, the ARR at 5 years was 15.1% (95% CI 4.0-26.3%) for participants with the highest 20% of predicted risk. Using the random forest model, the ARR at 5 years was 13.7% (95% CI 3.1-24.4%) for participants with the highest 20% of predicted risk. The highest outcome risk group in the decision tree model also exhibited a risk reduction, but the confidence interval was wider (5-year ARR = 17.0%, 95% CI= -5.4-39.4%). Cochran's Q test indicated ARR varied significantly only by subgroups created using the proportional hazards model. The hazard ratio for aspirin vs. placebo therapy did not significantly vary by subgroup in any of the models. The highest risk groups for the proportional hazards model and random forest model contained 230 participants each, while the highest risk group in the decision tree model contained 41 participants., Conclusions: The choice of technique for internally developed models for outcome risk subgroups influences HTE analyses. The rationale for the use of a particular subgroup determination model in HTE analyses needs to be explicitly defined based on desired levels of explainability (with features importance), uncertainty of prediction, chances of overfitting, and assumptions regarding the underlying data structure. Replication of these analyses using data from other mid-size clinical trials may help to establish guidance for selecting an outcomes risk prediction modelling technique for HTE analyses., (© 2024. The Author(s).)

Background: Mental health is an important fundamental element of health that influences different dimensions of an older individual's life. There are various factors that affect the mental health of older adults, such as resilience and social support. However, the relationship between various factors and mental health is unclear. Purposes of the study were to investigates whether resilience and social support affect the mental health of older adults and to analyze the structural equation model of mental health among the older adults with resilience as a mediating variable., Methods: A cross-sectional study was conducted among 964 older adults by multi-stage random sampling. The data were collected using the following questionnaires: Thai Elderly Resilience Scale, Social Support Scale, and Thai Geriatric Mental Health Assessment. Data analysis used descriptive statistics, Pearson's correlation, and Path analysis., Results: The older adults had a high level of resilience, social support, and mental health (MH). The best model we obtained demonstrated that social support and resilience had directly statistically significant influences the MH. Social support also had statistically significant indirect effect via resilience to the MH., Conclusion: Understanding the factors that influence mental health of older adults can be beneficial for health professionals to provide appropriate social support and resilience, including helping the older adult can join with others, be confident in life, have necessities, live with spiritual security, and be able to de-stress and management problems were recommended., (© 2024. The Author(s).)

When different researchers study the same research question using the same dataset they may obtain different and potentially even conflicting results. This is because there is often substantial flexibility in researchers' analytical choices, an issue also referred to as "researcher degrees of freedom". Combined with selective reporting of the smallest p-value or largest effect, researcher degrees of freedom may lead to an increased rate of false positive and overoptimistic results. In this paper, we address this issue by formalizing the multiplicity of analysis strategies as a multiple testing problem. As the test statistics of different analysis strategies are usually highly dependent, a naive approach such as the Bonferroni correction is inappropriate because it leads to an unacceptable loss of power. Instead, we propose using the "minP" adjustment method, which takes potential test dependencies into account and approximates the underlying null distribution of the minimal p-value through a permutation-based procedure. This procedure is known to achieve more power than simpler approaches while ensuring a weak control of the family-wise error rate. We illustrate our approach for addressing researcher degrees of freedom by applying it to a study on the impact of perioperative p a O 2 on post-operative complications after neurosurgery. A total of 48 analysis strategies are considered and adjusted using the minP procedure. This approach allows to selectively report the result of the analysis strategy yielding the most convincing evidence, while controlling the type 1 error-and thus the risk of publishing false positive results that may not be replicable., (© 2024. The Author(s).)

Background: Identification of human leukocyte antigen (HLA) types from DNA-sequenced human samples is important in organ transplantation and cancer immunotherapy and remains a challenging task considering sequence homology and extreme polymorphism of HLA genes., Results: We present Orthanq, a novel statistical model and corresponding application for transparent and uncertainty-aware quantification of haplotypes. We utilize our approach to perform HLA typing while, for the first time, reporting uncertainty of predictions and transparently observing mutations beyond reported HLA types. Using 99 gold standard samples from 1000 Genomes, Illumina Platinum Genomes and Genome In a Bottle projects, we show that Orthanq can provide overall superior accuracy and shorter runtimes than state-of-the-art HLA typers., Conclusions: Orthanq is the first approach that allows to directly utilize existing pangenome alignments and type all HLA loci. Moreover, it can be generalized for usages beyond HLA typing, e.g. for virus lineage quantification. Orthanq is available under https://orthanq.github.io ., (© 2024. The Author(s).)

Background: Smoking is the major risk factor for tracheal, bronchus, and lung (TBL) cancers. We investigated the feasibility of projecting TBL cancer incidence using smoking incidence rates by incorporating a range of latent periods from the main risk factor exposure to TBL cancer diagnosis., Methods: In this ecological study, we extracted data on TBL cancer incidence rates in Iran from 1990 to 2018 from the Global Burden of Disease (GBD) database. We also collected data on Iranian cigarette smoking patterns over the past 40 years through a literature review. The weighted average smoking incidence was calculated using a fixed-effects model with Comprehensive Meta-Analysis (CMA) software. Using these data, the five-year TBL cancer incidence in Iran was projected through time series modeling with IT Service Management (ITSM) 2000 software. A second model was developed based on cigarette smoking incidence using linear regression with SPSS (version 22), incorporating different latent periods. The results of these two models were compared to determine the best latent periods., Results: An increasing trend in TBL cancer incidence was observed from 2019 to 2023 (first model: 10.30 [95% CI: 9.62, 10.99] to 11.42 [95% CI: 10.85, 11.99] per 100,000 people). In the second model, the most accurate prediction was obtained with latent periods of 17 to 20 years, with the best prediction using a 17-year latent period (10.13 to 11.40 per 100,000 people) and the smallest mean difference of 0.08 (0.84%) per 100,000 people using the standard forecasting model (the ARIMA model)., Conclusion: Projecting an increase in TBL cancer incidence rates in the future, an optimal latent period of 17 to 20 years between exposure to cigarette smoke and TBL cancer incidence has implications for macrolevel preventive health policymaking to help reduce the burden of TBL cancer in upcoming years., (© 2024. The Author(s).)

Background: We aimed to estimate the age-specific and age-standardized incidence rate of diabetes for men and women in Mexico between 2003 and 2015, and to assess the relative change in incidence of diabetes between 2003 and 2015., Methods: We use a partial differential equation describing the illness-death model to estimate the incidence rate (IR) of diabetes for the years 2003, 2009 and 2015 based on prevalence data from National Health Surveys conducted in Mexico, the mortality rate of the Mexican general population and plausible input values for age-specific mortality rate ratios associated with diabetes., Results: The age-standardized IR of diabetes per 1000 person years (pryr) was similar among men (IRm) and women (IRw) in the year 2003 (IRm 6.1 vs. IRw 6.5 1000/pryr), 2009 (IRm: 7.0 vs. IRw: 8.4 1000/pryr), and in 2015 (IRm 8.0 vs. IRw 10.6 1000/pryr). The highest incident rates were observed among men and women in the 60-69 age group., Conclusions: Overall, the incidence rate of diabetes in Mexico between the years 2003 and 2015 remained stable. However, rates were markedly higher among women in the age group 40-49 and 50-59 in the year 2015 compared with rates in 2003., (© 2024. The Author(s).)

We propose a compartmental model for investigating smoking dynamics in an Italian region (Tuscany). Calibrating the model on local data from 1993 to 2019, we estimate the probabilities of starting and quitting smoking and the probability of smoking relapse. Then, we forecast the evolution of smoking prevalence until 2043 and assess the impact on mortality in terms of attributable deaths. We introduce elements of novelty with respect to previous studies in this field, including a formal definition of the equations governing the model dynamics and a flexible modelling of smoking probabilities based on cubic regression splines. We estimate model parameters by defining a two-step procedure and quantify the sampling variability via a parametric bootstrap. We propose the implementation of cross-validation on a rolling basis and variance-based Global Sensitivity Analysis to check the robustness of the results and support our findings. Our results suggest a decrease in smoking prevalence among males and stability among females, over the next two decades. We estimate that, in 2023, 18% of deaths among males and 8% among females are due to smoking. We test the use of the model in assessing the impact on smoking prevalence and mortality of different tobacco control policies, including the tobacco-free generation ban recently introduced in New Zealand., (© 2024. The Author(s).)

Background: Risk prediction models are routinely used to assist in clinical decision making. A small sample size for model development can compromise model performance when the model is applied to new patients. For binary outcomes, the calibration slope (CS) and the mean absolute prediction error (MAPE) are two key measures on which sample size calculations for the development of risk models have been based. CS quantifies the degree of model overfitting while MAPE assesses the accuracy of individual predictions., Methods: Recently, two formulae were proposed to calculate the sample size required, given anticipated features of the development data such as the outcome prevalence and c-statistic, to ensure that the expectation of the CS and MAPE (over repeated samples) in models fitted using MLE will meet prespecified target values. In this article, we use a simulation study to evaluate the performance of these formulae., Results: We found that both formulae work reasonably well when the anticipated model strength is not too high (c-statistic < 0.8), regardless of the outcome prevalence. However, for higher model strengths the CS formula underestimates the sample size substantially. For example, for c-statistic = 0.85 and 0.9, the sample size needed to be increased by at least 50% and 100%, respectively, to meet the target expected CS. On the other hand, the MAPE formula tends to overestimate the sample size for high model strengths. These conclusions were more pronounced for higher prevalence than for lower prevalence. Similar results were drawn when the outcome was time to event with censoring. Given these findings, we propose a simulation-based approach, implemented in the new R package 'samplesizedev', to correctly estimate the sample size even for high model strengths. The software can also calculate the variability in CS and MAPE, thus allowing for assessment of model stability., Conclusions: The calibration and MAPE formulae suggest sample sizes that are generally appropriate for use when the model strength is not too high. However, they tend to be biased for higher model strengths, which are not uncommon in clinical risk prediction studies. On those occasions, our proposed adjustments to the sample size calculations will be relevant., (© 2024. The Author(s).)

Background: Viral hepatitis imposes a heavy disease burden worldwide and is also one of the most serious public health problems in China. We aimed to describe the epidemiological characteristics of hepatitis in China and to investigate the influencing factors., Methods: We first used the JoinPoint model to analyze the percentage change (APC) and average annual percentage change (AAPC) of hepatitis in Chinese provinces from 2002 to 2021. We then explored the influencing factors by using the time-series global principal component analysis (GPCA) and the panel fixed-effects model., Results: The disease burden varied across different provinces from 2002 to 2021. The AAPC of the total HAV incidence decreased by 10.39% (95% CI: [-12.70%, -8.02%]) from 2002 to 2021. Yet the AAPC of HBV, HCV, and HEV increased by 1.50% (95% CI: [0.23%, 2.79%]), 13.99% (95% CI: [11.28%, 16.77%]), and 7.10% (95% CI: [0.90%, 13.69%]), respectively. The hotspots of HAV, HBV, HCV, and HEV moved from the west to the center, from the northwest to the southeast, from the northeast to the whole country, and from the northeast to the southeast, respectively. Different types of viral hepatitis infections were associated with hygiene, pollutant, and meteorological factors. Their roles in spatial-temporal incidence were expressed by panel regression functions., Conclusions: Viral hepatitis infection in China showed spatiotemporal heterogeneity. Interventions should be tailored to its epidemiological characteristics and determinants of viral hepatitis., (© 2024. The Author(s).)

This paper introduces a novel approach to modeling malaria incidence in Nigeria by integrating clustering strategies with regression modeling and leveraging meteorological data. By decomposing the datasets into multiple subsets using clustering techniques, we increase the number of explanatory variables and elucidate the role of weather in predicting different ranges of incidence data. Our clustering-integrated regression models, accompanied by optimal barriers, provide insights into the complex relationship between malaria incidence and well-established influencing weather factors such as rainfall and temperature.We explore two models. The first model incorporates lagged incidence and individual-specific effects. The second model focuses solely on weather components. Selection of a model depends on decision-makers priorities. The model one is recommended for higher predictive accuracy. Moreover, our findings reveal significant variability in malaria incidence, specific to certain geographic clusters and beyond what can be explained by observed weather variables alone.Notably, rainfall and temperature exhibit varying marginal effects across incidence clusters, indicating their differential impact on malaria transmission. High rainfall correlates with lower incidence, possibly due to its role in flushing mosquito breeding sites. On the other hand, temperature could not predict high-incidence cases, suggesting that other factors other than temperature contribute to high cases.Our study addresses the demand for comprehensive modeling of malaria incidence, particularly in regions like Nigeria where the disease remains prevalent. By integrating clustering techniques with regression analysis, we offer a nuanced understanding of how predetermined weather factors influence malaria transmission. This approach aids public health authorities in implementing targeted interventions. Our research underscores the importance of considering local contextual factors in malaria control efforts and highlights the potential of weather-based forecasting for proactive disease management., (© 2024. The Author(s).)

Background: Hospital outbreaks of multidrug resistant Pseudomonas aeruginosa are often caused by Pseudomonas aeruginosa clones which produce metallo-β-lactamases, such as Verona Integron-encoded Metallo-β-lactamase (VIM). Although different sources have been identified, the exact transmission routes often remain unknown. However, quantifying the role of different transmission routes of VIM-PA is important for tailoring infection prevention and control measures. The aim of this study is to quantify the relative importance of different transmission routes by applying a mathematical transmission model using admission and discharge dates as well as surveillance culture data of patients. Methods: We analyzed VIM-PA surveillance data collected between 2010 and 2018 of two intensive-care unit (ICU) wards for adult patients of the Erasmus University Medical Center Rotterdam using a mathematical transmission model. We distinguished two transmission routes: direct cross-transmission and a persistent environmental route. Based on admission, discharge dates, and surveillance cultures, we estimated the proportion of transmissions assigned to each of the routes. Results: Our study shows that only 13.7% (95% CI 1.4%, 29%) of the transmissions that occurred in these two ICU wards were likely caused by cross-transmission, leaving the vast majority of transmissions (86.3%, 95% CI 71%, 98.6%) due to persistent environmental contamination. Conclusions: Our results emphasize that persistent contamination of the environment may be an important driver of nosocomial transmissions of VIM-PA in ICUs. To minimize the transmission risk from the environment, potential reservoirs should be regularly and thoroughly cleaned and disinfected, or redesigned. [ABSTRACT FROM AUTHOR]

Background: Longitudinal ordinal data are commonly analyzed using a marginal proportional odds model for relating ordinal outcomes to covariates in the biomedical and health sciences. The generalized estimating equation (GEE) consistently estimates the regression parameters of marginal models even if the working covariance structure is misspecified. For small-sample longitudinal binary data, recent studies have shown that the bias of regression parameters may result from the GEE and have addressed the issue by applying Firth's adjustment for the likelihood score equation to the GEE as if generalized estimating functions were likelihood score functions. In this manuscript, for the proportional odds model for longitudinal ordinal data, the small-sample properties of the GEE were investigated, and a bias-reduced GEE (BR-GEE) was derived., Methods: By applying the adjusted function originally derived for the likelihood score function of the proportional odds model to the GEE, we produced the BR-GEE. We investigated the small-sample properties of both GEE and BR-GEE through simulation and applied them to a clinical study dataset., Results: In simulation studies, the BR-GEE had a bias closer to zero, smaller root mean square error than the GEE with coverage probability of confidence interval near or above the nominal level. The simulation also showed that BR-GEE maintained a type I error rate near or below the nominal level., Conclusions: For the analysis of longitudinal ordinal data involving a small number of subjects, the BR-GEE is advantageous for obtaining estimates of the regression parameters of marginal proportional odds models., (© 2024. The Author(s).)

Background and Aim: Pancreatic cancer possesses a high prevalence and mortality rate among other cancers. Despite the low survival rate of this cancer type, the early prediction of this disease has a crucial role in decreasing the mortality rate and improving the prognosis. So, this study., Materials and Methods: In this retrospective study, we used 654 alive and dead PC cases to establish the prediction model for PC. The six chosen machine learning algorithms and prognostic factors were utilized to build the prediction models. The importance of the predictive factors was assessed using the relative importance of a high-performing algorithm., Results: The XG-Boost with AU-ROC of 0.933 (95% CI= [0.906-0.958]) and AU-ROC of 0.836 (95% CI= [0.789-0.865] in internal and external validation modes were considered as the best-performing model for predicting the mortality risk of PC. The factors, including tumor size, smoking, and chemotherapy, were considered the most influential for prediction., Conclusion: The XG-Boost gained more performance efficiency in predicting the mortality risk of PC patients, so this model can promote the clinical solutions that doctors can achieve in healthcare environments to decrease the mortality risk of these patients., (© 2024. The Author(s).)

Background: HALE is now a regular strategic planning indicator for all levels of the Chinese government. However, HALE measurements necessitate comprehensive data collection and intricate technology. Therefore, effectively converting numerous diseases into the years lived with disability (YLD) rate is a significant challenge for HALE measurements. Our study aimed to construct a simple YLD rate measurement model with high applicability based on the current situation of actual data resources within China to address challenges in measuring HALE target values during planning., Methods: First, based on the Chinese YLD rate in the Global Burden of Disease (GBD) 2019, Pearson correlation analysis, the global optimum method, etc., was utilized to screen the best predictor variables from the current Chinese data resources. Missing data for predictor variables were filled in via spline interpolation. Then, multiple linear regression models were fitted to construct the YLD rate measurement model. The Sullivan method was used to measure HALE. The Monte Carlo method was employed to generate 95% uncertainty intervals. Finally, model performances were assessed using the mean absolute error (MAE) and mean absolute percentage error (MAPE)., Results: A three-input-parameter model was constructed to measure the age-specific YLD rates by sex in China, directly using the incidence of infectious diseases, the incidence of chronic diseases among persons aged 15 and older, and the addition of an under-five mortality rate covariate. The total MAE and MAPE for the combined YLD rate were 0.0007 and 0.5949%, respectively. The MAE and MAPE of the combined HALE in the 0-year-old group were 0.0341 and 0.0526%, respectively. There were slightly fewer males (0.0197, 0.0311%) than females (0.0501, 0.0755%)., Conclusion: We constructed a high-accuracy model to measure the YLD rate in China by using three monitoring indicators from the Chinese national routine as predictor variables. The model provides a realistic and feasible solution for measuring HALE at the national and especially regional levels, considering limited data., (© 2024. The Author(s).)

Investigators often conduct randomized controlled trials (RCTs) at multiple centers/sites when determining the effect of a treatment or an intervention. Diversifying recruitment across multiple institutions allows investigators to make recruitment go faster within a shorter timeframe and allows generalizing the study results across diverse populations. Despite having a common study protocol across multiple centers, the eligible participants may be heterogeneous, site policies and practices may vary, and the investigators' experience, training, and expertise may also vary across sites. These factors may contribute to the heterogeneity in effect estimates across centers. As a result, we usually observe some degree of heterogeneity in effect estimates across centers, despite all centers following the same study protocol. During the analysis of such a trial, investigators typically ignore center effects, but some have suggested considering centers as fixed or random effects in the model. It is not clear how considering the effects of centers, either as fixed or random effects, impacts the test of the primary hypothesis. In this article, we first review the practice of accounting for center effects in the analyses of published RCTs and illustrate the extent of heterogeneity observed in a few preexisting multicenter RCTs. To determine the impact of heterogeneity on the test of a primary hypothesis of an RCT, we considered continuous and binary outcomes and the corresponding appropriate model, namely, a simple linear regression model for a continuous outcome and a logistic regression model for the binary outcome. For each model type, we considered three methods: (a) ignore the center effect, (b) account for centers as fixed effects, or (c) account for centers as random effects. Based on simulation studies of these models, we then examine whether considering the center as a fixed or random effect in the model helps to preserve or reduce the type I and type II error rates during the analysis phase of an RCT. Finally, we outline the threshold at which center-level effects are negligible and thus negligible and provide recommendations on when it may be necessary to account for center effects during the analyses of multicenter randomized controlled trials., (© 2024. The Author(s).)

Objective: To compare how different imputation methods affect the estimates and performance of a prediction model for premature mortality., Study Design and Setting: Sex-specific Weibull accelerated failure time survival models were run on four separate datasets using complete case, mode, single and multiple imputation to impute missing values. Six performance measures were compared to access predictive accuracy (Nagelkerke R 2 , integrated brier score), discrimination (Harrell's c-index, discrimination slope) and calibration (calibration in the large, calibration slope)., Results: The highest proportion of missingness for a single variable was 10.86% for the female model and 8.24% for the male model. Comparing the performance measures for complete case, mode, single and multiple imputation: the Nagelkerke R 2 values for the female model was 0.1084, 0.1116, 0.1120 and 0.111-0.1120 with the male model exhibited similar variation of 0.1050, 0.1078, 0.1078 and 0.1078-0.1081. Harrell's c-index also demonstrated small variation with values of 0.8666, 0.8719, 0.8719 and 0.8711-0.8719 for the female model and 0.8549, 0.8548, 0.8550 and 0.8550-0.8553 for the male model., Conclusion: In the scenarios examined in this study, mode imputation performed well when using a population health survey compared to single and multiple imputation when predictive performance measures is the main model goal. To generate unbiased hazard ratios, multiple imputation methods were superior. This study shows the need to consider the best imputation approach for a predictive model development given the conditions of missing data and the goals of the analysis., (© 2024. The Author(s).)

Background: Causal mediation analysis plays a crucial role in examining causal effects and causal mechanisms. Yet, limited work has taken into consideration the use of sampling weights in causal mediation analysis. In this study, we compared different strategies of incorporating sampling weights into causal mediation analysis., Methods: We conducted a simulation study to assess 4 different sampling weighting strategies-1) not using sampling weights, 2) incorporating sampling weights into mediation "cross-world" weights, 3) using sampling weights when estimating the outcome model, and 4) using sampling weights in both stages. We generated 8 simulated population scenarios comprising an exposure (A), an outcome (Y), a mediator (M), and six covariates (C), all of which were binary. The data were generated so that the true model of A given C and the true model of A given M and C were both logit models. We crossed these 8 population scenarios with 4 different sampling methods to obtain 32 total simulation conditions. For each simulation condition, we assessed the performance of 4 sampling weighting strategies when calculating sample-based estimates of the total, direct, and indirect effects. We also applied the four sampling weighting strategies to a case study using data from the National Survey on Drug Use and Health (NSDUH)., Results: Using sampling weights in both stages (mediation weight estimation and outcome models) had the lowest bias under most simulation conditions examined. Using sampling weights in only one stage led to greater bias for multiple simulation conditions., Discussion: Using sampling weights in both stages is an effective approach to reduce bias in causal mediation analyses under a variety of conditions regarding the structure of the population data and sampling methods., (© 2024. The Author(s).)

Background: Dynamical mathematical models defined by a system of differential equations are typically not easily accessible to non-experts. However, forecasts based on these types of models can help gain insights into the mechanisms driving the process and may outcompete simpler phenomenological growth models. Here we introduce a friendly toolbox, SpatialWavePredict, to characterize and forecast the spatial wave sub-epidemic model, which captures diverse wave dynamics by aggregating multiple asynchronous growth processes and has outperformed simpler phenomenological growth models in short-term forecasts of various infectious diseases outbreaks including SARS, Ebola, and the early waves of the COVID-19 pandemic in the US., Results: This tutorial-based primer introduces and illustrates a user-friendly MATLAB toolbox for fitting and forecasting time-series trajectories using an ensemble spatial wave sub-epidemic model based on ordinary differential equations. Scientists, policymakers, and students can use the toolbox to conduct real-time short-term forecasts. The five-parameter epidemic wave model in the toolbox aggregates linked overlapping sub-epidemics and captures a rich spectrum of epidemic wave dynamics, including oscillatory wave behavior and plateaus. An ensemble strategy aims to improve forecasting performance by combining the resulting top-ranked models. The toolbox provides a tutorial for forecasting time-series trajectories, including the full uncertainty distribution derived through parametric bootstrapping, which is needed to construct prediction intervals and evaluate their accuracy. Functions are available to assess forecasting performance, estimation methods, error structures in the data, and forecasting horizons. The toolbox also includes functions to quantify forecasting performance using metrics that evaluate point and distributional forecasts, including the weighted interval score., Conclusions: We have developed the first comprehensive toolbox to characterize and forecast time-series data using an ensemble spatial wave sub-epidemic wave model. As an epidemic situation or contagion occurs, the tools presented in this tutorial can facilitate policymakers to guide the implementation of containment strategies and assess the impact of control interventions. We demonstrate the functionality of the toolbox with examples, including a tutorial video, and is illustrated using daily data on the COVID-19 pandemic in the USA., (© 2024. The Author(s).)

Background: Time-varying exposures like pet ownership pose challenges for identifying critical windows due to multicollinearity when modeled simultaneously. The Distributed Lag Model (DLM) estimates critical windows for time-varying exposures, which are mainly continuous variables. However, applying complex functions such as high-order splines and nonlinear functions within DLMs may not be suitable for situations with limited time points or binary exposure, such as in questionnaire surveys., Objectives: (1) We examined the estimation performance of a simple DLM with fractional polynomial function for time-varying binary exposures through simulation experiments. (2) We evaluated the impact of pet ownership on childhood wheezing onset and estimate critical windows., Methods: (1) We compared logistic regression including time-varying exposure in separate models, in one model simultaneously, and using DLM. For evaluation, we employed bias, empirical standard error (EmpSE), and mean squared error (MSE). (2) The Japan Environment and Children's Study (JECS) is a prospective birth cohort study of approximately 100,000 parent-child pairs, registered across Japan from 2011 to 2014. We applied DLM to the JECS data up to age 3. The estimated odds ratios (OR) were considered to be within critical windows when they were significant at the 5% level., Results: (1) DLM and the separate model exhibited lower bias compared to the simultaneously model. Additionally, both DLM and the simultaneously model demonstrated lower EmpSEs than the separate model. In all scenarios, DLM had lower MSEs than the other methods. Specifically, where critical windows is clearly present and exposure correlation is high, DLM showed MSEs about 1/2 to 1/200 of those of other models. (2) Application of DLM to the JECS data showed that, unlike other models, a significant exposure effect was observed only between the ages of 0 and 6 months. During that periods, the highest ORs were 1.07 (95% confidence interval, 1.01 to 1.14) , observed between the ages of 2 and 5 months., Conclusions: (1) A simple DLM improves the accuracy of exposure effect and critical windows estimation. (2) 0-6 months may be the critical windows for the effect of pet ownership on the wheezing onset at 3 years., (© 2024. The Author(s).)

Background: Influenza is a highly contagious respiratory disease that presents a significant challenge to public health globally. Therefore, effective influenza prediction and prevention are crucial for the timely allocation of resources, the development of vaccine strategies, and the implementation of targeted public health interventions., Method: In this study, we utilized historical influenza case data from January 2013 to December 2021 in Fuzhou to develop four regression prediction models: SARIMA, Prophet, Holt-Winters, and XGBoost models. Their predicted performance was assessed by using influenza data from the period from January 2022 to December 2022 in Fuzhou. These models were used for fitting and prediction analysis. The evaluation metrics, including Mean Squared Error (MSE), Root Mean Squared Error (RMSE), and Mean Absolute Error (MAE), were employed to compare the performance of these models., Results: The results indicate that the epidemic of influenza in Fuzhou exhibits a distinct seasonal and cyclical pattern. The influenza cases data displayed a noticeable upward trend and significant fluctuations. In our study, we employed SARIMA, Prophet, Holt-Winters, and XGBoost models to predict influenza outbreaks in Fuzhou. Among these models, the XGBoost model demonstrated the best performance on both the training and test sets, yielding the lowest values for MSE, RMSE, and MAE among the four models., Conclusion: The utilization of the XGBoost model significantly enhances the prediction accuracy of influenza in Fuzhou. This study makes a valuable contribution to the field of influenza prediction and provides substantial support for future influenza response efforts., (© 2024. The Author(s).)

Background: Extended illness-death models (a specific class of multistate models) are a useful tool to analyse situations like hospital-acquired infections, ventilation-associated pneumonia, and transfers between hospitals. The main components of these models are hazard rates and transition probabilities. Calculation of different measures and their interpretation can be challenging due to their complexity., Methods: By assuming time-constant hazards, the complexity of these models becomes manageable and closed mathematical forms for transition probabilities can be derived. Using these forms, we created a tool in R to visualize transition probabilities via stacked probability plots., Results: In this article, we present this tool and give some insights into its theoretical background. Using published examples, we give guidelines on how this tool can be used. Our goal is to provide an instrument that helps obtain a deeper understanding of a complex multistate setting., Conclusion: While multistate models (in particular extended illness-death models), can be highly complex, this tool can be used in studies to both understand assumptions, which have been made during planning and as a first step in analysing complex data structures. An online version of this tool can be found at https://eidm.imbi.uni-freiburg.de/ ., (© 2024. The Author(s).)

Background: Nested case-control (NCC) designs are efficient for developing and validating prediction models that use expensive or difficult-to-obtain predictors, especially when the outcome is rare. Previous research has focused on how to develop prediction models in this sampling design, but little attention has been given to model validation in this context. We therefore aimed to systematically characterize the key elements for the correct evaluation of the performance of prediction models in NCC data., Methods: We proposed how to correctly evaluate prediction models in NCC data, by adjusting performance metrics with sampling weights to account for the NCC sampling. We included in this study the C-index, threshold-based metrics, Observed-to-expected events ratio (O/E ratio), calibration slope, and decision curve analysis. We illustrated the proposed metrics with a validation of the Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA version 5) in data from the population-based Rotterdam study. We compared the metrics obtained in the full cohort with those obtained in NCC datasets sampled from the Rotterdam study, with and without a matched design., Results: Performance metrics without weight adjustment were biased: the unweighted C-index in NCC datasets was 0.61 (0.58-0.63) for the unmatched design, while the C-index in the full cohort and the weighted C-index in the NCC datasets were similar: 0.65 (0.62-0.69) and 0.65 (0.61-0.69), respectively. The unweighted O/E ratio was 18.38 (17.67-19.06) in the NCC datasets, while it was 1.69 (1.42-1.93) in the full cohort and its weighted version in the NCC datasets was 1.68 (1.53-1.84). Similarly, weighted adjustments of threshold-based metrics and net benefit for decision curves were unbiased estimates of the corresponding metrics in the full cohort, while the corresponding unweighted metrics were biased. In the matched design, the bias of the unweighted metrics was larger, but it could also be compensated by the weight adjustment., Conclusions: Nested case-control studies are an efficient solution for evaluating the performance of prediction models that use expensive or difficult-to-obtain biomarkers, especially when the outcome is rare, but the performance metrics need to be adjusted to the sampling procedure., (© 2024. The Author(s).)

Background: Identifying the impact of environmental mixtures on human health is an important topic. However, such studies face challenges when exposure measurements lie below limit of detection (LOD). While various approaches for accommodating a single exposure subject to LOD have been used, their impact on mixture analysis has not been thoroughly investigated. Our study aims to understand the impact of five popular LOD accommodation approaches on mixture analysis results with multiple exposures subject to LOD, including omitting subjects with any exposures below LOD (complete case analysis); single imputations by LOD/ 2 , and by estimates from a censored accelerated failure time (AFT) model; and multiple imputation (MI) with or without truncation based on LOD., Methods: In extensive simulation studies with high-dimensional and highly correlated exposures and a continuous health outcome, we examined the performance of each LOD approach on three mixture analysis methods: elastic net regression, weighted quantile sum regression (WQS) and Bayesian kernel machine regression (BKMR). We further analyzed data from the National Health and Nutrition Examination Survey (NHANES) on how persistent organic pollutants (POPs) influenced leukocyte telomere length (LTL)., Results: Complete case analysis was inefficient and could result in severe bias for some mixture methods. Imputation by LOD/ 2 showed unstable performance across mixture methods. Conventional MI was associated with consistent mild biases, which can be reduced by using a truncated distribution for imputation. Estimating censored values by AFT models had a minimal impact on the results. In the NHANES analysis, imputation by LOD/ 2 , truncated MI and censored AFT models performed similarly, with a positive overall effect of POPs on LTL while PCB126, PCB169 and furan 2,3,4,7,8-pncdf being the most important exposures., Conclusions: Our study favored using truncated MI and censored AFT models to accommodate values below LOD for the stability of downstream mixture analysis., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)