On the Reversible Jump Markov Chain Monte Carlo (RJMCMC) Algorithm for Extreme Value Mixture Distribution as a Location-Scale Transformation of the Weibull Distribution

Data with a multimodal pattern can be analyzed using a mixture model. In a mixture model, the most important step is the determination of the number of mixture components, because finding the correct number of mixture components will reduce the error of the resulting model. In a Bayesian analysis, one method that can be used to determine the number of mixture components is the reversible jump Markov chain Monte Carlo (RJMCMC). The RJMCMC is used for distributions that have location and scale parameters or location-scale distribution, such as the Gaussian distribution family. In this research, we added an important step before beginning to use the RJMCMC method, namely the modification of the analyzed distribution into location-scale distribution. We called this the non-Gaussian RJMCMC (NG-RJMCMC) algorithm. The following steps are the same as for the RJMCMC. In this study, we applied it to the Weibull distribution. This will help many researchers in the field of survival analysis since most of the survival time distribution is Weibull. We transformed the Weibull distribution into a location-scale distribution, which is the extreme value (EV) type 1 (Gumbel-type for minima) distribution. Thus, for the mixture analysis, we call this EV-I mixture distribution. Based on the simulation results, we can conclude that the accuracy level is at minimum 95%. We also applied the EV-I mixture distribution and compared it with the Gaussian mixture distribution for enzyme, acidity, and galaxy datasets. Based on the Kullback–Leibler divergence (KLD) and visual observation, the EV-I mixture distribution has higher coverage than the Gaussian mixture distribution. We also applied it to our dengue hemorrhagic fever (DHF) data from eastern Surabaya, East Java, Indonesia. The estimation results show that the number of mixture components in the data is four; we also obtained the estimation results of the other parameters and labels for each observation. Based on the Kullback–Leibler divergence (KLD) and visual observation, for our data, the EV-I mixture distribution offers better coverage than the Gaussian mixture distribution.

Modelling the patient accrual with truncated Gaussian mixture distribution for the accurate estimation of sample size in survival trials

In survival trials with fixed trial length, the patient accrual rate has a significant impact on the sample size estimation or equivalently, on the power of trials. A larger sample size is required for the staggered patient entry. During enrollment, the patient accrual rate changes with the recruitment publicity effect, disease incidence and many other factors and fluctuations of the accrual rate occur frequently. However, the existing accrual models are either over-simplified for the constant rate assumption or complicated in calculation for the subdivision of the accrual period. A more flexible accrual model is required to represent the fluctuant patient accrual rate for accurate sample size estimation. In this paper, inspired by the flexibility of the Gaussian mixture distribution in approximating continuous densities, we propose the truncated Gaussian mixture distribution accrual model to represent different variations of accrual rate by different parameter configurations. The sample size calculation formula and the parameter setting of the proposed accrual model are discussed further.