The simultaneous clustering of observations and features of datasets (known as co-clustering) has recently emerged as a central topic in machine learning applications. However, most models focus on continuous data in stationary scenarios, where cluster assignments do not evolve over time. We propose in this paper the dynamic latent block model (dLBM), which extends the classical binary latent block model, making amenable such analysis to dynamic cases where data are counts. Our approach operates on temporal count matrices allowing to detect abrupt changes in the way existing clusters interact with each other. The time breaks detection is performed through clustering of time instants, that allows for better model parsimony. The time dependent counting data are modeled via non-homogeneous Poisson processes (HHPPs), conditionally to the latent variables. In order to handle the model inference, we rely on a SEM-Gibbs algorithm and the ICL criterion is used for model selection. Numerical experiments on simulated data highlight the main features of the proposed approach and show the interest of dLBM with respect to related works. An application to adverse drug reaction in pharmacovigilance is also proposed, where dLBM was able to recognize clusters in a meaningful way that identified safety events that were consistent with retrospective knowledge. Hence, our aim is to propose this dynamic co-clustering method as a tool for automatic safety signal detection, to support medical authorities.

(1) The aim of our study is to evaluate the capacity of the Visually AcceSAble Rembrandt Images (VASARI) scoring system in discerning between the different degrees of glioma and Isocitrate Dehydrogenase (IDH) status predictions, with a possible application in machine learning. (2) A retrospective study was conducted on 126 patients with gliomas (M/F = 75/51; mean age: 55.30), from which we obtained their histological grade and molecular status. Each patient was analyzed with all 25 features of VASARI, blinded by two residents and three neuroradiologists. The interobserver agreement was assessed. A statistical analysis was conducted to evaluate the distribution of the observations using a box plot and a bar plot. We then performed univariate and multivariate logistic regressions and a Wald test. We also calculated the odds ratios and confidence intervals for each variable and the evaluation matrices with receiver operating characteristic (ROC) curves in order to identify cut-off values that are predictive of a diagnosis. Finally, we did the Pearson correlation test to see if the variables grade and IDH were correlated. (3) An excellent ICC estimate was obtained. For the grade and IDH status prediction, there were statistically significant results by evaluation of the degree of post-contrast impregnation (F4) and the percentage of impregnated area (F5), not impregnated area (F6), and necrotic (F7) tissue. These models showed good performances according to the area under the curve (AUC) values (>70%). (4) Specific MRI features can be used to predict the grade and IDH status of gliomas, with important prognostic implications. The standardization and improvement of these data (aim: AUC > 80%) can be used for programming machine learning software.