Graphical Tests for the Frailty Distribution in the Shared Frailty Model

We extend a diagnostic plot for the frailty distribution in proportional hazards models to the case of shared frailty. The plot is based on a closure property of exponential family failure distributions with canonical statistics z and g(z), namely that the frailty distribution among survivors at time t has the same form, with the same values of the parameters associated with g(z). We extend this property to shared frailty, considering various definitions of a “surviving” cluster at time t. We illustrate the effectiveness of the method in the case where the “death” of the cluster is defined by the first death among its members.     

The use of length‐biased distributions in statistical monitoring

Statistical process monitoring (SPM) has been used extensively recently in order to assure the quality of the output of industrial processes. Techniques of SPM have been efficiently applied during the last two decades in non‐industrial processes. A field of application with great interest is public health monitoring, where a pitfall with which we have to deal is the fact that available samples are not random in all cases. In the majority of cases, we monitor measurements derived from patient admissions to a hospital against control limits that were calculated using a sample of data taken from an epidemiological survey. In this work, we bridge the gap of a change in the sampling scheme from Phase I to Phase II, studying the case where the sampling during Phase II is biased. We present the appropriate methodology and then apply extensive numerical simulation in order to explore the performance of the proposed methodology, for measurements following various asymmetrical distributions. As the simulations show, the proposed methodology has a significantly better performance than the standard procedure.     

Join the Union and Be Safe: The Effects of Unionization on Occupational Safety and Health in the European Union

This paper investigates the effect of unionization on occupational safety and health, as measured by the fatal and non‐fatal work accidents, after controlling for the country's gross domestic product. It uses a panel sample of 10 European Union countries, for the period 1982–2006. The study takes into account the time persistence in work injuries and the endogenous nature of the work injuries–union density relationship. In addition, the effect of union density is decomposed into a temporary and permanent effect. It is shown that increasing union density is associated with a decrease in the number of both fatal and non‐fatal work injuries.     

Modelling survival data using mixtures of frailties

Frailty models are often used to describe the extra heterogeneity in survival data by introducing an individual random, unobserved effect. The frailty term is usually assumed to act multiplicatively on a baseline hazard function common to all individuals. In order to apply the frailty model, a specific frailty distribution has to be assumed. If at least one of the latent variables is continuous, the frailty must follow a continuous distribution. In this paper, a finite mixture of continuous frailty distributions is used in order to describe situations in which one (or more) of the latent variables separates the population in study into two (or more) subpopulations. Closure properties of the unobserved quantity are given along with the maximum-likelihood estimates under the most common choices of frailty distributions. The model is illustrated on a set of lifetime data.     

Sample Tests for Detection of Size-Biased Sampling Mechanism

The aim of this article is to present a new test for the detection of size-bias in a sample with or without censored observations. The test is simple in the form and demands only the knowledge of consistent estimators of any nuisance parameters appeared in the model. With the use of simulated samples from the Weibull distribution, we show the advantages of the new test compared to the Likelihood Ratio and the Wald test.     

On the generalized cumulative residual entropy weighted distributions

Recently, Feizjavadian and Hashemi (2015 Feizjavadian, S.H., Hashemi, R. (2015). Mean residual weighted versus the length-biased Rayleigh distribution. J. Stat. Comput. Simul. 85:2823–2838.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) introduced and studied the mean residual weighted (MRW) distribution as an alternative to the length-biased distribution, by using the concepts of the mean residual lifetime and the cumulative residual entropy (CRE). In this article, a new sequence of weighted distributions is introduced based on the generalized CRE. This sequence includes the MRW distribution. Properties of this sequence are obtained generalizing and extending previous results on the MRW distribution. Moreover, expressions for some known distributions are given, and finite mixtures between the new sequence of weighted distributions and the length-biased distribution are studied. Numerical examples are given to illustrate the new results.     

Fitting parametric frailty and mixture models under biased sampling

Biased sampling from an underlying distribution with p.d.f. f(t), t>0, implies that observations follow the weighted distribution with p.d.f. f ^w (t)=w(t)f(t)/E[w(T)] for a known weight function w. In particular, the function w(t)=t ^α has important applications, including length-biased sampling (α=1) and area-biased sampling (α=2). We first consider here the maximum likelihood estimation of the parameters of a distribution f(t) under biased sampling from a censored population in a proportional hazards frailty model where a baseline distribution (e.g. Weibull) is mixed with a continuous frailty distribution (e.g. Gamma). A right-censored observation contributes a term proportional to w(t)S(t) to the likelihood; this is not the same as S ^w (t), so the problem of fitting the model does not simply reduce to fitting the weighted distribution. We present results on the distribution of frailty in the weighted distribution and develop an EM algorithm for estimating the parameters of the model in the important Weibull–Gamma case. We also give results for the case where f(t) is a finite mixture distribution. Results are presented for uncensored data and for Type I right censoring. Simulation results are presented, and the methods are illustrated on a set of lifetime data.     

On Small Samples Testing for Frailty Through Homogeneity Test

We derive a test in order to examine the need of modeling survival data using frailty models based on the likelihood ratio (LR) test for homogeneity. Test is developed for both complete and censored samples from a family of baseline distributions that satisfy a closure property. Approach motivated by I-divergence distance is used in order to determine “credible” regions for all parameters of baseline distribution for which homogeneity hypothesis is not rejected. Proposed test outperforms the usual asymptotic LR test both in very small samples with known frailty and for all small sample sizes under misspecified frailty.     

On the Bayesian analysis of the mixture of power function distribution using the complete and the censored sample

The power function distribution is often used to study the electrical component reliability. In this paper, we model a heterogeneous population using the two-component mixture of the power function distribution. A comprehensive simulation scheme including a large number of parameter points is followed to highlight the properties and behavior of the estimates in terms of sample size, censoring rate, parameters size and the proportion of the components of the mixture. The parameters of the power function mixture are estimated and compared using the Bayes estimates. A simulated mixture data with censored observations is generated by probabilistic mixing for the computational purposes. Elegant closed form expressions for the Bayes estimators and their variances are derived for the censored sample as well as for the complete sample. Some interesting comparison and properties of the estimates are observed and presented. The system of three non-linear equations, required to be solved iteratively for the computations of maximum likelihood (ML) estimates, is derived. The complete sample expressions for the ML estimates and for their variances are also given. The components of the information matrix are constructed as well. Uninformative as well as informative priors are assumed for the derivation of the Bayes estimators. A real-life mixture data example has also been discussed. The posterior predictive distribution with the informative Gamma prior is derived, and the equations required to find the lower and upper limits of the predictive intervals are constructed. The Bayes estimates are evaluated under the squared error loss function.