Testing for a Unit Root in a Time Series With a Changing Mean: Corrections and Extensions

This article explores the relation between nonexponential waiting times between events and the distribution of the number of events in a fixed time interval. It is shown that within this framework the frequently observed phenomenon of overdispersion—that is, a variance that exceeds the mean—is caused by a decreasing hazard function of the waiting times, whereas an increasing hazard function leads to underdispersion. Using the assumption of iid gamma-distributed waiting times, a new count-data model is derived. Its use is illustrated in two applications, the number of births and the number of doctor consultations.     

The Kumaraswamy generalized gamma distribution with application in survival analysis

We introduce and study the so-called Kumaraswamy generalized gamma distribution that is capable of modeling bathtub-shaped hazard rate functions. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a large number of well-known lifetime special sub-models such as the exponentiated generalized gamma, exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma, generalized Rayleigh, among others. Some structural properties of the new distribution are studied. We obtain two infinite sum representations for the moments and an expansion for the generating function. We calculate the density function of the order statistics and an expansion for their moments. The method of maximum likelihood and a Bayesian procedure are adopted for estimating the model parameters. The usefulness of the new distribution is illustrated in two real data sets.     

The log-exponentiated generalized gamma regression model for censored data

For the first time, we introduce a generalized form of the exponentiated generalized gamma distribution [Cordeiro et al. The exponentiated generalized gamma distribution with application to lifetime data, J. Statist. Comput. Simul. 81 (2011), pp. 827–842.] that is the baseline for the log-exponentiated generalized gamma regression model. The new distribution can accommodate increasing, decreasing, bathtub- and unimodal-shaped hazard functions. A second advantage is that it includes classical distributions reported in the lifetime literature as special cases. We obtain explicit expressions for the moments of the baseline distribution of the new regression model. The proposed model can be applied to censored data since it includes as sub-models several widely known regression models. It therefore can be used more effectively in the analysis of survival data. We obtain maximum likelihood estimates for the model parameters by considering censored data. We show that our extended regression model is very useful by means of two applications to real data.     

Detection of frailty in weibull lifetime data using outlier tests

Heterogeneity in lifetime data may be modelled by multiplying an individual's hazard by an unobserved frailty. We test for the presence of frailty of this kind in univariate and bivariate data with Weibull distributed lifetimes, using statistics based on the ordered Cox–Snell residuals from the null model of no frailty. The form of the statistics is suggested by outlier testing in the gamma distribution. We find through simulation that the sum of the k largest or k smallest order statistics, for suitably chosen k, provides a powerful test when the frailty distribution is assumed to be gamma or positive stable, respectively. We provide recommended values of k for sample sizes up to 100 and simple formulae for estimated critical values for tests at the 5% level.     

A compound beta distribution with applications in finance

If X and Y are gamma distributed independent random variables then it is well known that the ratio X / (X + Y) has the beta distribution. In this note, the distribution of W = X / (X + Y) is considered when X and Y have the compound gamma distribution. We refer to the distribution of W as compound beta and describe an application to consumer price indices to show that compound beta is a better model than one based on the standard beta distribution. We derive various properties of W, including its probability density function, cumulative distribution function, hazard rate function and moments.     

Some gamma distributions

Three new generalizations of the standard gamma distribution introduced by the author are reviewed. Various properties are derived for each distribution, including its hazard rate function and moments. An application is illustrated to drought data.     

The exponentiated generalized gamma distribution with application to lifetime data

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.     

Performance analysis of the retrial queues with finite number of sources and service interruptions

This paper aims at presenting an analytic approach for investigating a single-server retrial queue with finite population of customers where the server is subject to interruptions. A free source may generate a primary call to request service. If the server is free upon arrival, the call starts to be served and the service times are independent, generally distributed random variables. During the service time the source cannot generate a new primary call. After service the source moves into the free state and can generate a new primary call. There is no waiting space in front of the server, and a call who finds the server unavailable upon arrival joins an orbit of unsatisfied customers. The server is subject to interruptions during the service processes. When the server is interrupted, the call being served just before server interruption goes to the retrial orbit and will retry its luck after a random amount of time until it finds the server available. The recovery times of the interrupted server are assumed to be generally distributed. Our analysis extends previous work on this topic and includes the analysis of the arriving customer’s distribution, the busy period, and the waiting time process.     

On a recent generalization of gamma distribution

In this paper, we investigate a generalized gamma distribution recentIy developed by Agarwal and Kalla (1996). Also, we show that such generalized distribution, like the ordinary gamma distribution, has a unique mode and, unlike the ordinary gamma distribution, may have a hazard rate (mean residual life) function which is upside-down bathtub (bathtub) shaped.     

Multi-class M/PH/1 queues with deterministic impatience times

This article considers computational procedures for the waiting time and queue length distributions in stationary multi-class first-come, first-served single-server queues with deterministic impatience times. There are several classes of customers, which are distinguished by deterministic impatience times (i.e., maximum allowable waiting times). We assume that customers in each class arrive according to an independent Poisson process and a single server serves customers on a first-come, first-served basis. Service times of customers in each class are independent and identically distributed according to a phase-type distribution that may differ for different classes. We first consider the stationary distribution of the virtual waiting time and then derive numerically feasible formulas for the actual waiting time distribution and loss probability. We also analyze the joint queue length distribution and provide an algorithmic procedure for computing the probability mass function of the stationary joint queue length.     

A relationship between a family of waiting time distributions and the asymptotic residual waiting time distribution

For a class of renewal process waiting time distributions defined herein, one may describe the distribution of asymptotic residual waiting times. The relationship between the two distributions characterizes the class, which includes the gamma distribution. Possible consequences for hypothesis testing are discussed.     

Modeling Multilevel Survival Data Using Frailty Models

Often the dependence in multivariate survival data is modeled through an individual level effect called the frailty. Due to its mathematical simplicity, the gamma distribution is often used as the frailty distribution for hazard modeling. However, it is well known that the gamma frailty distribution has many drawbacks. For example, it weakens the effect of covariates. In addition, in the presence of a multilevel model, overall frailty comes from several levels. To overcome such drawbacks, more heavy-tailed distributions are needed to model the frailty distribution in order to incorporate extra variability. In this article, we develop a class of log-skew-t distributions for the frailty. This class includes the log-normal distribution along with many other heavy tailed distributions, e.g., log-Cauchy, log normal, and log-t as special cases.

Conditional on the frailty, the survival times are assumed to be independent with proportional hazard structure. The modeling process is then completed by assuming multilevel frailty-effects. Instead of tuning a strict parameterization of the baseline hazard function, we consider the partial likelihood approach and thus leave the baseline function unspecified. By eliminating the hazard, the pre-specification and computation are simplified considerably.     

Estimation in Semiparametric Marginal Shared Gamma Frailty Models

The semiparametric marginal shared frailty models in survival analysis have the non–parametric hazard functions multiplied by a random frailty in each cluster, and the survival times conditional on frailties are assumed to be independent. In addition, the marginal hazard functions have the same form as in the usual Cox proportional hazard models. In this paper, an approach based on maximum likelihood and expectation–maximization is applied to semiparametric marginal shared gamma frailty models, where the frailties are assumed to be gamma distributed with mean 1 and variance θ. The estimates of the fixed–effect parameters and their standard errors obtained using this approach are compared in terms of both bias and efficiency with those obtained using the extended marginal approach. Similarly, the standard errors of our frailty variance estimates are found to compare favourably with those obtained using other methods. The asymptotic distribution of the frailty variance estimates is shown to be a 50–50 mixture of a point mass at zero and a truncated normal random variable on the positive axis for θ₀ = 0. Simulations demonstrate that, for θ₀ < 0, it is approximately an x −(100 − x )%, 0 ≤ x ≤ 50, mixture between a point mass at zero and a truncated normal random variable on the positive axis for small samples and small values of θ₀; otherwise, it is approximately normal.     

An extension of the exponential distribution

A generalization of the exponential distribution is presented. The generalization always has its mode at zero and yet allows for increasing, decreasing and constant hazard rates. It can be used as an alternative to the gamma, Weibull and exponentiated exponential distributions. A comprehensive account of the mathematical properties of the generalization is presented. A real data example is discussed to illustrate its applicability.     

ASYMPTOTIC AND SMALL SAMPLE STATISTICAL PROPERTIES OF RANDOM FRAILTY VARIANCE ESTIMATES FOR SHARED GAMMA FRAILTY MODELS

This paper concerns maximum likelihood estimation for the semiparametric shared gamma frailty model; that is the Cox proportional hazards model with the hazard function multiplied by a gamma random variable with mean 1 and variance θ. A hybrid ML-EM algorithm is applied to 26 400 simulated samples of 400 to 8000 observations with Weibull hazards. The hybrid algorithm is much faster than the standard EM algorithm, faster than standard direct maximum likelihood (ML, Newton Raphson) for large samples, and gives almost identical results to the penalised likelihood method in S-PLUS 2000. When the true value θ₀ of θ is zero, the estimates of θ are asymptotically distributed as a 50–50 mixture between a point mass at zero and a normal random variable on the positive axis. When θ₀ > 0, the asymptotic distribution is normal. However, for small samples, simulations suggest that the estimates of θ are approximately distributed as an x ? (100 ? x)% mixture, 0 ≤ x ≤ 50, between a point mass at zero and a normal random variable on the positive axis even for θ₀ > 0. In light of this, p-values and confidence intervals need to be adjusted accordingly. We indicate an approximate method for carrying out the adjustment.     

A Bayesian piecewise exponential phase II design for monitoring a time-to-event endpoint

A robust Bayesian design is presented for a single-arm phase II trial with an early stopping rule to monitor a time to event endpoint. The assumed model is a piecewise exponential distribution with non-informative gamma priors on the hazard parameters in subintervals of a fixed follow up interval. As an additional comparator, we also define and evaluate a version of the design based on an assumed Weibull distribution. Except for the assumed models, the piecewise exponential and Weibull model based designs are identical to an established design that assumes an exponential event time distribution with an inverse gamma prior on the mean event time. The three designs are compared by simulation under several log-logistic and Weibull distributions having different shape parameters, and for different monitoring schedules. The simulations show that, compared to the exponential inverse gamma model based design, the piecewise exponential design has substantially better performance, with much higher probabilities of correctly stopping the trial early, and shorter and less variable trial duration, when the assumed median event time is unacceptably low. Compared to the Weibull model based design, the piecewise exponential design does a much better job of maintaining small incorrect stopping probabilities in cases where the true median survival time is desirably large.     

Compound patterns and generating functions: From basic waiting times to counts of occurrence

We derive an explicit, closed form expression for the double generating function of the corresponding counts of occurrence, within a finite time horizon, of the single patterns contained in a compound pattern. The expression is in terms of a basic single, and a basic joint, generating functions for which exact solutions exist in the literature. The single generating function is associated with the basic waiting time for the first occurrence of the compound pattern. The joint generating function is that for the waiting time to reach a given single pattern and the associated counts of occurrence, within that waiting time, of the single patterns contained in the compound pattern. The literature on patterns is huge. Also, there are results that establish links between generating functions for counts of occurrence of the single patterns contained in a compound pattern with generating functions of some more complex waiting times associated with that compound pattern. The latter waiting times are known in the literature with names such as sooner, or later waiting times, or generalisations of such. On the other hand, our result fills a gap in the literature by providing a neat link connecting the generating functions of the basic quantities associated with occurrence of compound patterns.     

Generalized exponential distributions

The three-parameter gamma and three-parameter Weibull distributions are commonly used for analysing any lifetime data or skewed data. Both distributions have several desirable properties, and nice physical interpretations. Because of the scale and shape parameters, both have quite a bit of flexibility for analysing different types of lifetime data. They have increasing as well as decreasing hazard rate depending on the shape parameter. Unfortunately both distributions also have certain drawbacks. This paper considers a three-parameter distribution which is a particular case of the exponentiated Weibull distribution originally proposed by Mudholkar, Srivastava & Freimer (1995) when the location parameter is not present. The study examines different properties of this model and observes that this family has some interesting features which are quite similar to those of the gamma family and the Weibull family, and certain distinct properties also. It appears this model can be used as an alternative to the gamma model or the Weibull model in many situations. One dataset is provided where the three-parameter generalized exponential distribution fits better than the three-parameter Weibull distribution or the three-parameter gamma distribution.     

The generalized symmetric Laplace distribution on the sphere

A distribution on the unit sphere is generated by conditioning a normal mixture distribution with an inverse gamma distributed weighting function. It can be regarded as the generalized symmetric Laplace distribution on the unit sphere. The density involves a modified Bessel function of the third kind which can be approximated by other simpler functions in certain limiting cases. As a consequence, the von Mises–Fisher, cardioid and Jones–Pewsey distributions are limiting cases of the new distribution. No closed form expressions exist for the roots of the likelihood equations. However, given the normal mixture structure of the distribution, we propose an E–M-algorithm-based approach for finding the maximum-likelihood estimates of the parameters which assumes the weights in the mixture to be missing data. The modeling capabilities of the proposed distribution are illustrated by fitting it and some of its competitors to two circular data sets.     

BAYESIAN HIDDEN MARKOV MODELS FOR LONGITUDINAL COUNTS

This paper describes a Bayesian approach to modelling carcinogenity in animal studies where the data consist of counts of the number of tumours present over time. It compares two autoregressive hidden Markov models. One of them models the transitions between three latent states: an inactive transient state, a multiplying state for increasing counts and a reducing state for decreasing counts. The second model introduces a fourth tied state to describe non‐zero observations that are neither increasing nor decreasing. Both these models can model the length of stay upon entry of a state. A discrete constant hazards waiting time distribution is used to model the time to onset of tumour growth. Our models describe between‐animal‐variability by a single hierarchy of random effects and the within‐animal variation by first‐order serial dependence. They can be extended to higher‐order serial dependence and multi‐level hierarchies. Analysis of data from animal experiments comparing the influence of two genes leads to conclusions that differ from those of Dunson (2000). The observed data likelihood defines an information criterion to assess the predictive properties of the three‐ and four‐state models. The deviance information criterion is appropriately defined for discrete parameters.