Moderate deviations for the moment estimators in Rayleigh distribution with two parameters

ABSTRACT

We study the moderate deviations of the moment estimators in Rayleigh distribution with two parameters. The moderate deviations are obtained by the delta method in large deviation principle.     

A General Class of Distributions: Properties and Applications

ABSTRACT

In this article, we derive a general class of distributions and establish its relationship to χ² distribution. The proposed class includes normal, inverse Gaussian, lognormal, gamma, Rayleigh, and Maxwell distributions. Various statistical properties of the class are discussed. Some applications of the class are given.     

Bayesian prediction of the total time on test using doubly censored rayleigh data

In a Bayesian setting, and on the basis of a doubly censored random sample of failure times drawn from a Rayleigh distribution, Fernandez (2000, Statist. Probab. Lett. , 48 , 393-399) considered the problem of predicting an independent future sample from the same distribution. In this article, we extend his work to include the estimation of the predictive distribution of the total time on test up to a certain failure in a future sample, as well as that of the remaining testing time time until all the items in the original sample have failed. Two examples are used to illustrate the prediction procedure.     

Mean residual weighted versus the length-biased Rayleigh distribution

In some statistical applications, data may not be considered as a random sample of the whole population and some subjects have less probability of belonging to the sample. Consequently, statistical inferences for such data sets, usually yields biased estimation. In such situations, the length-biased version of the original random variable as a special weighted distribution often produces better inferences. An alternative weighted distribution based on the mean residual life is suggested to treat the biasedness. The Rayleigh distribution is applied in many real applications, hence the proposed method of weighting is performed to produce a new lifetime distribution based on the Rayleigh model. In addition, statistical properties of the proposed distribution is investigated. A simulation study and a real data set are prepared to illustrate that the mean residual weighted Rayleigh distribution gives a better fit than the original and also the length-biased Rayleigh distribution.     

Nonparametric Estimation of the Conditional Distribution at Regression Boundary Points

Abstract

Nonparametric regression is a standard statistical tool with increased importance in the Big Data era. Boundary points pose additional difficulties but local polynomial regression can be used to alleviate them. Local linear regression, for example, is easy to implement and performs quite well both at interior and boundary points. Estimating the conditional distribution function and/or the quantile function at a given regressor point is immediate via standard kernel methods but problems ensue if local linear methods are to be used. In particular, the distribution function estimator is not guaranteed to be monotone increasing, and the quantile curves can “cross.” In the article at hand, a simple method of correcting the local linear distribution estimator for monotonicity is proposed, and its good performance is demonstrated via simulations and real data examples. Supplementary materials for this article are available online.     

A note on the most powerful invariant test of Rayleigh against exponential distribution

Abstract

The problem of testing Rayleigh distribution against exponentiality, based on a random sample of observations is considered. This problem arises in survival analysis, when testing a linearly increasing hazard function against a constant hazard function. It is shown that for this problem the most powerful invariant test is equivalent to the “ratio of maximized likelihoods” (RML) test. However, since the two families are separate, the RML test statistic does not have the usual asymptotic chi-square distribution. Normal and saddlepoint approximations to the distribution of the RML test statistic are derived. Simulations show that saddlepoint approximation is more accurate than the normal approximation, especially for tail probabilities that are the main values of interest in hypothesis testing.     

Bayesian inference for the generalized exponential distribution

The two-parameter generalized exponential (GE) distribution was introduced by Gupta and Kundu [Gupta, R.D. and Kundu, D., 1999, Generalized exponential distribution. Australian and New Zealand Journal of Statistics, 41(2), 173–188.]. It was observed that the GE can be used in situations where a skewed distribution for a nonnegative random variable is needed. In this article, the Bayesian estimation and prediction for the GE distribution, using informative priors, have been considered. Importance sampling is used to estimate the parameters, as well as the reliability function, and the Gibbs and Metropolis samplers data sets are used to predict the behavior of further observations from the distribution. Two data sets are used to illustrate the Bayesian procedure.     

ON BAYES PREDICTION OF FUTURE MEDIAN

An expression for the Bayesian predictive survival function of the median of a set of future observations is obtained whether its size is assumed to be odd or even. Both of the informative and future samples are drawn from a population whose distribution is a general class that includes several distributions used in life testing (and other areas as well) such as the Weibull (including the exponential and Rayleigh), compound Weibull (including the compound exponential and compound Rayleigh), Pareto, beta, Gompertz and compound Gompertz, among other distributions. A general proper (conjugate) prior density function is used to cover most prior distributions that have been used in literature. Applications to the Weibull, exponential and Rayleigh models are illustrated.     

Bayes and robust Bayes predictions in a subfamily of scale parameters under a precautionary loss function

ABSTRACT

This paper deals with Bayes, robust Bayes, and minimax predictions in a subfamily of scale parameters under an asymmetric precautionary loss function. In Bayesian statistical inference, the goal is to obtain optimal rules under a specified loss function and an explicit prior distribution over the parameter space. However, in practice, we are not able to specify the prior totally or when a problem must be solved by two statisticians, they may agree on the choice of the prior but not the values of the hyperparameters. A common approach to the prior uncertainty in Bayesian analysis is to choose a class of prior distributions and compute some functional quantity. This is known as Robust Bayesian analysis which provides a way to consider the prior knowledge in terms of a class of priors Γ for global prevention against bad choices of hyperparameters. Under a scale invariant precautionary loss function, we deal with robust Bayes predictions of Y based on X. We carried out a simulation study and a real data analysis to illustrate the practical utility of the prediction procedure.     

On the matrix-variate generalized hyperbolic distribution and its Bayesian applications

In the first part of the paper, we introduce the matrix-variate generalized hyperbolic distribution by mixing the matrix normal distribution with the matrix generalized inverse Gaussian density. The p-dimensional generalized hyperbolic distribution of [Barndorff-Nielsen, O. (1978). Hyperbolic distributions and distributions on hyperbolae. Scand. J. Stat., 5, 151–157], the matrix-T distribution and many well-known distributions are shown to be special cases of the new distribution. Some properties of the distribution are also studied. The second part of the paper deals with the application of the distribution in the Bayesian analysis of the normal multivariate linear model.     

Goodness of fit test for the generalized Rayleigh distribution with unknown parameters

The use of goodness-of-fit test based on Anderson–Darling (AD) statistic is discussed, with reference to the composite hypothesis that a sample of observations comes from a generalized Rayleigh distribution whose parameters are unspecified. Monte Carlo simulation studies were performed to calculate the critical values for AD test. These critical values are then used for testing whether a set of observations follows a generalized Rayleigh distribution when the scale and shape parameters are unspecified and are estimated from the sample. Functional relationship between the critical values of AD is also examined for each shape parameter (α), sample size (n) and significance level (γ). The power study is performed with the hypothesized generalized Rayleigh against alternate distributions.     

A spatial scan statistic for survival data based on generalized life distribution

ABSTRACT

For many years, detection of clusters has been of great public health interest and widely studied. Several methods have been developed to detect clusters and their performance has been evaluated in various contexts. Spatial scan statistics are widely used for geographical cluster detection and inference. Different types of discrete or continuous data can be analyzed using spatial scan statistics for Bernoulli, Poisson, ordinal, exponential, and normal models. In this paper, we propose a scan statistic for survival data which is based on generalized life distribution model that provides three important life distributions, viz. Weibull, exponential, and Rayleigh. The proposed method is applied to the survival data of tuberculosis patients in Nainital district of Uttarakhand, India, for the year 2004–05. The Monte Carlo simulation studies reveal that the proposed method performs well for different survival distributions.     

Estimation of R=P[Y<X] for three-parameter generalized Rayleigh distribution

Surles and Padgett [Inference for reliability and stress–strength for a scaled Burr type X distribution. Lifetime Data Anal. 2001;7:187–200] introduced a two-parameter Burr-type X distribution, which can be described as a generalized Rayleigh distribution. In this paper, we consider the estimation of the stress–strength parameter R=P[Y<X], when X and Y are both three-parameter generalized Rayleigh distributions with the same scale and locations parameters but different shape parameters. It is assumed that they are independently distributed. It is observed that the maximum-likelihood estimators (MLEs) do not exist, and we propose a modified MLE of R. We obtain the asymptotic distribution of the modified MLE of R, and it can be used to construct the asymptotic confidence interval of R. We also propose the Bayes estimate of R and the construction of the associated credible interval based on importance sampling technique. Analysis of two real data sets, (i) simulated and (ii) real, have been performed for illustrative purposes.     

Spatial Regression Models Using Inter-Region Distances in a Non-Random Context

This article considers spatial data z( s ₁), z( s ₂),…, z( s _n ) collected at n locations, with the objective of predicting z( s ₀) at another location. The usual method of analysis for this problem is kriging, but here we introduce a new signal-plus-noise model whose essential feature is the identification of hot spots. The signal decays in relation to distance from hot spots. We show that hot spots can be located with high accuracy and that the decay parameter can be estimated accurately. This new model compares well to kriging in simulations.     

A new lifetime model: the Kumaraswamy generalized Rayleigh distribution

The generalized Rayleigh (GR) distribution [V.G. Vodǎ, Inferential procedures on a generalized Rayleigh variate, I, Appl. Math. 21 (1976), pp. 395–412; V.G. Vodǎ, Inferential procedures on a generalized Rayleigh variate, II, Appl. Math. 21 (1976), pp. 413–419] has been applied in several areas such as health, agriculture, biology and other sciences. For the first time, we propose the Kumaraswamy GR (KwGR) distribution for analysing lifetime data. The new density function can be expressed as a mixture of GR density functions. Explicit formulae are derived for some of its statistical quantities. The density function of the order statistics can be expressed as a mixture of GR density functions. We also propose a linear log-KwGR regression model for analysing data with real support to extend some known regression models. The estimation of parameters is approached by maximum likelihood. The importance of the new models is illustrated in two real data sets.     

Analysis of progressively type-I interval censored competing risks data for a class of an exponential distribution

ABSTRACT

In this paper, we consider some problems of point estimation and point prediction when the competing risks data from a class of exponential distribution are progressive type-I interval censored. The maximum likelihood estimation and mid-point approximation method are proposed for the estimations of parameters. Also several point predictors of censored units such as the maximum likelihood predictor, the best unbiased predictor and the conditional median predictor are obtained. The methods discussed here are applied when the lifetime distributions of the latent failure times are independent and Weibull-distributed. Finally a simulation study is given by using Monte-Carlo simulations to compare the performances of the different methods and one data analysis has been presented for illustrative purposes.     

Bayesian prediction of minimal repair times of a series system based on hybrid censored sample of components’ lifetimes under Rayleigh distribution

In this paper, we develop Bayesian predictive inferential procedures for prediction of repair times of a series system, applying a minimal repair strategy, using the information contained in an independent observed hybrid censored sample of the lifetimes of the components of the system, assuming the underlying distribution of the lifetimes to be Rayleigh distribution. An illustrative real data example and a simulation study are presented for the purpose of illustration and comparison of the proposed predictors.     

Prediction based on linear combinations of order statistics and bivariate concomitants in the case of multivariate elliptical distributions

In this paper, by considering a (3n+1) -dimensional random vector (X₀, X^T, Y^T, Z^T)^T having a multivariate elliptical distribution, we derive the exact joint distribution of (X₀, a^TX_(n), b^TY_[n], c^TZ_[n])^T, where a, b, c∈?ⁿ, X_(n)=(X₍₁₎, …, X_(n))^T, X₍₁₎<···<X_(n), is the vector of order statistics arising from X, and Y_[n]=(Y_[1], …, Y_[n])^T and Z_[n]=(Z_[1], …, Z_[n])^T denote the vectors of concomitants corresponding to X_(n) ((Y_[r], Z_[r])^T, for r=1, …, n, is the vector of bivariate concomitants corresponding to X_(r)). We then present an alternate approach for the derivation of the exact joint distribution of (X₀, X_(r), Y_[r], Z_[r])^T, for r=1, …, n. We show that these joint distributions can be expressed as mixtures of four-variate unified skew-elliptical distributions and these mixture forms facilitate the prediction of X_(r), say, based on the concomitants Y_[r] and Z_[r]. Finally, we illustrate the usefulness of our results by a real data.     

Bayesian predictive distribution for a Poisson model with a parametric restriction

Abstract

Predictive probability estimation for a Poisson distribution is addressed when the parameter space is restricted. The Bayesian predictive probability against the prior on the restricted space is compared with the non-restricted Bayes predictive probability. It is shown that the former predictive probability dominates the latter under some conditions when the predictive probabilities are evaluated by the risk function relative to the Kullback-Leibler divergence. This result is proved by first showing the corresponding dominance result for estimating the restricted parameter and then translating it into the framework of predictive probability estimation.     

Bayesian estimation and prediction based on generalized Type-II hybrid censored sample

ABSTRACT

In this paper, we consider a general form for the underlying distribution and a general conjugate prior, and develop a general procedure for deriving the maximum likelihood and Bayesian estimators based on an observed generalized Type-II hybrid censored sample. The problems of predicting the future order statistics from the same sample and that from a future sample are also discussed from a Bayesian viewpoint. For the illustration of the developed results, the exponential and Pareto distributions are used as examples. Finally, two numerical examples are presented for illustrating all the inferential procedures developed here.