Multivariate normal mean–variance mixture distribution based on Birnbaum–Saunders distribution

A multivariate normal mean–variance mixture based on a Birnbaum–Saunders (NMVMBS) distribution is introduced and several properties of this new distribution are discussed. A new robust non-Gaussian ARCH-type model is proposed in which there exists a relation between the variance of the observations, and the marginal distributions are NMVMBS. A simple EM-based maximum likelihood estimation procedure to estimate the parameters of this normal mean–variance mixture distribution is given. A simulation study and some real data are used to demonstrate the modelling strength of this new model.     

Estimation and simulation of conditional hazard function in the quasi-associated framework when the observations are linked via a functional single-index structure

We consider the estimation of the conditional hazard function of a scalar response variable Y given a Hilbertian random variable X when the observations are linked via a single-index structure in the quasi-associated framework. We establish the pointwise almost complete convergence and the uniform almost complete convergence (with the rate) of the estimate of this model. A simulation is given to illustrate the good behavior in the practice of our methodology.     

Conditional Akaike Information Criteria for a Class of Poisson Mixture Models with Random Effects

Focusing on the model selection problems in the family of Poisson mixture models (including the Poisson mixture regression model with random effects and zero‐inflated Poisson regression model with random effects), the current paper derives two conditional Akaike information criteria. The criteria are the unbiased estimators of the conditional Akaike information based on the conditional log‐likelihood and the conditional Akaike information based on the joint log‐likelihood, respectively. The derivation is free from the specific parametric assumptions about the conditional mean of the true data‐generating model and applies to different types of estimation methods. Additionally, the derivation is not based on the asymptotic argument. Simulations show that the proposed criteria have promising estimation accuracy. In addition, it is found that the criterion based on the conditional log‐likelihood demonstrates good model selection performance under different scenarios. Two sets of real data are used to illustrate the proposed method.     

The beta generalized linear exponential distribution

In this paper, a new five-parameter lifetime distribution called beta generalized linear exponential distribution (BGLED) is introduced. It includes at least 17 popular sub-models as special cases such as the beta linear exponential, the beta generalized exponential, and the exponentiated generalized linear distributions. Mathematical and statistical properties of the proposed distribution are discussed in details. In particular, explicit expression for the density function, moments, asymptotics distributions for sample extreme statistics, and other statistical measures are obtained. The estimation of the parameters by the method of maximum-likelihood is discussed and the finite sample properties of the maximum-likelihood estimators (MLEs) are investigated numerically. A real data set is used to demonstrate its superior performance fit over several existing popular lifetime models.     

Sum,product and ratio of Pareto and gamma variables

The exact distributions of X+Y, X Y and X/(X+Y) are studied when X and Y are independent Pareto and gamma random variables. Applications are discussed, to real problems in clinical trials, computer networks and economics.     

Generalised K-S confidence regions: some exact results

One-sided confidence regions for continuous cumulative distribution function are constructed using empirical cumulative distribution functions and the generalized Kolmogorov-Smirnov distance. The band width of such regions becomes narrower in the right or left tail of the distribution. Significance levels necessary for implementation are given. Some other K-S type distances useful in constructing a confidence region with nonconstant width are also included.     

Simulating random variables using moment-generating functions and the saddlepoint approximation

When we are given only a transform such as the moment-generating function of a distribution, it is rare that we can efficiently simulate random variables. Possible approaches such as the inverse transform using numerical inversion of the transform are computationally very expensive. However, the saddlepoint approximation is known to be exact for the Normal, Gamma, and inverse Gaussian distribution and remarkably accurate for a large number of others. We explore the efficient use of the saddlepoint approximation for simulating distributions and provide three examples of the accuracy of these simulations.     

Double shrinkage estimators for large sparse covariance matrices

Covariance matrices play an important role in many multivariate techniques and hence a good covariance estimation is crucial in this kind of analysis. In many applications a sparse covariance matrix is expected due to the nature of the data or for simple interpretation. Hard thresholding, soft thresholding, and generalized thresholding were therefore developed to this end. However, these estimators do not always yield well-conditioned covariance estimates. To have sparse and well-conditioned estimates, we propose doubly shrinkage estimators: shrinking small covariances towards zero and then shrinking covariance matrix towards a diagonal matrix. Additionally, a richness index is defined to evaluate how rich a covariance matrix is. According to our simulations, the richness index serves as a good indicator to choose relevant covariance estimator.     

Rank repeated measures analysis of covariance

When analyzing a response variable at the presence of both factors and covariates, with potentially correlated responses and violated assumptions of the normal residual or the linear relationship between the response and the covariates, rank-based tests can be an option for inferential procedures instead of the parametric repeated measures analysis of covariance (ANCOVA) models. This article derives a rank-based method for multi-way ANCOVA models with correlated responses. The generalized estimating equations (GEE) technique is employed to construct the proposed rank tests. Asymptotic properties of the proposed tests are derived. Simulation studies confirmed the performance of the proposed tests.     

Decomposition of Variance for Spatial Cox Processes

Abstract. Spatial Cox point processes is a natural framework for quantifying the various sources of variation governing the spatial distribution of rain forest trees. We introduce a general criterion for variance decomposition for spatial Cox processes and apply it to specific Cox process models with additive or log linear random intensity functions. We moreover consider a new and flexible class of pair correlation function models given in terms of normal variance mixture covariance functions. The proposed methodology is applied to point pattern data sets of locations of tropical rain forest trees.