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.     

Environmental risk assessment model for low dose extrapolation ∗

Use of experimental data from animal studies to estimate human risk due to long-term exposure to very low doses of chemicals in the environment poses a number of biological and statistical problems. One of the statistical problems is to extrapolate the animal dose-response relation from the high dose levels where data are available to low dose, which humans might encounter. Here, a quantal dose-response model is developed based on a multi-hit theory of toxic response. The development of the model utilizes a weighted Lagrange-Poisson distribution for the number of hits. When spontaneous background toxic response is included, the model involves three unknown parameters. The maximum likelihood estimators for these parameters are given as the solution of a nonlinear iterative algorithm. The use of this model for low-dose extrapolation is indicated. The results are applied to nine sets of toxic response data.     

Nonlinear voxel-based modelling of the haemodynamic response in fMRI

A common assumption for data analysis in functional magnetic resonance imaging is that the response signal can be modelled as the convolution of a haemodynamic response (HDR) kernel with a stimulus reference function. Early approaches modelled spatially constant HDR kernels, but more recently spatially varying models have been proposed. However, convolution limits the flexibility of these models and their ability to capture spatial variation. Here, a range of (nonlinear) parametric curves are fitted by least squares minimisation directly to individual voxel HDRs (i.e., without using convolution). A ‘constrained gamma curve’ is proposed as an efficient form for fitting the HDR at each individual voxel. This curve allows for spatial variation in the delay of the HDR, but places a global constraint on the temporal spread. The approach of directly fitting individual parameters of HDR shape is demonstrated to lead to an improved fit of response estimates.     

On univariate and bivariate generalized gamma convolutions

This paper has two parts. In the first part some results for generalized gamma convolutions (GGCs) are reviewed. A GGC is a limit distribution for sums of independent gamma variables. In the second part, bivariate gamma distributions and bivariate GGCs are considered. New bivariate gamma distributions are derived from shot-noise models. The remarkable property hyperbolic complete monotonicity (HCM) for a function is considered both in the univariate case and in the bivariate case.     

A Family of Near-Exact Distributions Based on Truncations of the Exact Distribution for the Generalized Wilks Lambda Statistic

For the case where at least two sets have an odd number of variables we do not have the exact distribution of the generalized Wilks Lambda statistic in a manageable form, adequate for manipulation. In this article, we develop a family of very accurate near-exact distributions for this statistic for the case where two or three sets have an odd number of variables. We first express the exact characteristic function of the logarithm of the statistic in the form of the characteristic function of an infinite mixture of Generalized Integer Gamma distributions. Then, based on truncations of this exact characteristic function, we obtain a family of near-exact distributions, which, by construction, match the first two exact moments. These near-exact distributions display an asymptotic behaviour for increasing number of variables involved. The corresponding cumulative distribution functions are obtained in a concise and manageable form, relatively easy to implement computationally, allowing for the computation of virtually exact quantiles. We undertake a comparative study for small sample sizes, using two proximity measures based on the Berry-Esseen bounds, to assess the performance of the near-exact distributions for different numbers of sets of variables and different numbers of variables in each set.     

The gamma-Lomax distribution

For any continuous baseline G distribution, Zografos and Balakrishnan [On families of beta- and generalized gamma-generated distributions and associated inference. Statist Methodol. 2009;6:344–362] proposed a generalized gamma-generated distribution with an extra positive parameter. A new three-parameter continuous distribution called the gamma-Lomax distribution, which extends the Lomax distribution is proposed and studied. Various structural properties of the new distribution are derived including explicit expressions for the moments, generating and quantile functions, mean deviations and Rényi entropy. The estimation of the model parameters is performed by maximum likelihood. We also determine the observed information matrix. An application illustrates the usefulness of the proposed model.     

A testing approach to estimation

Results from the theory of uniformly most powerful invariant tests are used to develop a new parameter estimation procedure. The procedure is used to derive parameter estimators for several important distributions. Results of simulation studies comparing the performances of the new estimators and maximum likelihood estimators are presented.     

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.     

Exact Sampling from a Continuous State Space

Propp & Wilson (1996) described a protocol, called coupling from the past, for exact sampling from a target distribution using a coupled Markov chain Monte Carlo algorithm. In this paper we extend coupling from the past to various MCMC samplers on a continuous state space; rather than following the monotone sampling device of Propp & Wilson, our approach uses methods related to gamma-coupling and rejection sampling to simulate the chain, and direct accounting of sample paths.     

Inference procedures for the variance gamma model and applications

Goodness-of-fit tests for the family of the four-parameter normal–variance gamma distribution are constructed. The tests are based on a weighted integral incorporating the empirical characteristic function of suitably standardized data. Non-standard algorithms are employed for the computation of the maximum-likelihood estimators of the parameters involved in the test statistic, while Monte Carlo results are used in order to compare the new test with some classical goodness-of-fit methods. A real-data application is also included.