In clinical research, study outcomes usually consist of various patients’ information corresponding to the treatment. To have a better understanding of the effects of different treatments, one often needs to analyze multiple clinical outcomes simultaneously, while the data are usually mixed with both continuous and discrete variables. We propose the multivariate mixed response model to implement statistical inference based on the conditional grouped continuous model through a pairwise composite-likelihood approach. It can simplify the multivariate model by dealing with three types of bivariate models and incorporating the asymptotical properties of the composite likelihood via the Godambe information. We demonstrate the validity and the statistic power of the multivariate mixed response model through simulation studies and clinical applications. This composite-likelihood method is advantageous for statistical inference on correlated multivariate mixed outcomes.
The Asian Emerald Ash Borer beetle (EAB, Agrilus planipennis Fairmaire) can cause damage to all species of Ash trees (Fraxinus), and rampant, unchecked infestations of this insect can cause significant damage to forests. It is thus critical to assess and model the spread of the EAB in a manner that allows authorities to anticipate likely areas of future tree infestation. In this study, a generalized linear mixed model (GLMM), combining the features of the commonly used generalized linear model (GLM) and a random effects model, was developed to predict future EAB spread patterns in Southern Ontario, Canada. The GLMM was designed to deal with autocorrelation in the data. Two random effects were established based on the geographic information provided with the EAB data, and a method based on statistical inference was proposed to identify the most significant factors associated with the distribution of the EAB. The results of the model showed that 95% of the testing data were correctly classified. The predictive performance of the GLMM was substantially enhanced in comparison with that obtained by the GLM. The influence of climatic factors, such as wind speed and anthropogenic activities, had the most significant influence on the spread of the EAB.
10 September 2025 in Canadian Journal of Statistics
When data from multiple tasks have outlier contamination, existing multitask learning methods perform less efficiently. To address this issue, we propose a robust multitask feature learning method by combining the adaptive Huber regression tasks with mixed regularization. The robustification parameters can be chosen to adapt to the sample size, model dimension, and moments of the error distribution while striking a balance between unbiasedness and robustness. We consider heavy‐tailed distributions for multiple datasets that have bounded th moment for any . Our method can achieve estimation and sign recovery consistency. Additionally, we propose a robust information criterion to conduct joint inference on related tasks, which can be used for consistent model selection. Through different simulation studies and real data applications, we illustrate that the performance of the proposed model can provide smaller estimation errors and higher feature selection accuracy than the non robust multitask learning and robust single‐task methods.
Background: In clinical trials, evaluating de-intensified oncologic treatment strategies can help reduce treatment-related toxicities while preserving patients’ quality of life. However, de-intensification is typically evaluated in cancers with a low relapse rate, and if the cancer type is uncommon, a randomized trial may require an impractically extended period to accumulate sufficient events for reliable inferential conclusions. Method: This paper introduces a Bayesian adaptive method for the single-arm trial design that provides efficient analysis of survival data under these constraints. By incorporating data from previous studies to establish prior knowledge and a historical control arm, this approach enables robust and accurate estimations and predictions for trial design, sample size determination, and inferential decision-making. To support the implementation of this method, we developed an R package called “BayesAT,” which offers significant flexibility in modelling and supports multi-stage interim analyses, particularly for evaluating de-intensified oncologic treatments. Result: Our approach is validated through comprehensive simulation studies and sensitivity analyses. Additionally, this algorithm has been applied to a pediatric Hodgkin lymphoma trial, showcasing its capability to effectively leverage information from previous studies and conduct interim analyses that expedite conclusions regarding treatment efficacy.