Properties of doubly-truncated gamma variables

The truncated gamma distribution has been widely studied, primarily in life-testing and reliability settings. Most work has assumed an upper bound on the support of the random variable, i.e. the space of the distribution is (0,u). We consider a doubly-truncated gamma random variable restricted by both a lower (l) and upper (u) truncation point, both of which are considered known. We provide simple forms for the density, cumulative distribution function (CDF), moment generating function, cumulant generating function, characteristic function, and moments. We extend the results to describe the density, CDF, and moments of a doubly-truncated noncentral chi-square variable.     

Tests for generalized exponential laws based on the empirical Mellin transform

Recently, many standard families of distributions have been generalized by exponentiating their cumulative distribution function (CDF). In this paper, test statistics are constructed based on CDF–transformed observations and the corresponding moments of arbitrary positive order. Simulation results for generalized exponential distributions show that the proposed test compares well with standard methods based on the empirical distribution function.     

The Correlated Bivariate Noncentral F Distribution and Its Application

This article proposes the singly and doubly correlated bivariate noncentral F (BNCF) distributions. The probability density function (pdf) and the cumulative distribution function (cdf) of the distributions are derived for arbitrary values of the parameters. The pdf and cdf of the distributions for different arbitrary values of the parameters are computed, and their graphs are plotted by writing and implementing new R codes. An application of the correlated BNCF distribution is illustrated in the computations of the power function of the pre-test test for the multivariate simple regression model (MSRM).     

Enhancement for two commonly-used approximations for the inverse cumulative function of the normal distribution

Two commonly used approximations for the inverse distribution function of the normal distribution are Schmeiser's and Shore's. Both approximations are based on a power transformation of either the cumulative density function (CDF) or a simple function of it. In this note we demonstrate, that if these approximations are presented in the form of the classical one-parameter Box-Cox transformation, and the exponent of the transformation is expressed as a simple function of the CDF, then the accuracy of both approximations may be considerably enhanced, without losing much in algebraic simplicity. Since both approximations are special cases of more general four-parameter systems of distributions, the results presented here indicate that the accuracy of the latter, when used to represent non-normal density functions, may also be considerably enhanced.     

On singular elliptical models

The probability density function (pdf) ofsingular elliptical distributions is represented as an integralseries of singular normal distributions. Explicit formulas for the pdf and the cdf of the generalized Chi-square distribution are derived under singular elliptical assumptions extending the result of Díaz-García [(2002). Singular elliptical distribution: density and applications. Commun. Stat.—Theory Methods 31:665–681]. Applications are given of the proposed result for singular mixedmodels.     

On estimation of the PDF and CDF of the Lindley distribution

This article addresses two methods of estimation of the probability density function (PDF) and cumulative distribution function (CDF) for the Lindley distribution. Following estimation methods are considered: uniformly minimum variance unbiased estimator (UMVUE) and maximum likelihood estimator (MLE). Since the Lindley distribution is more flexible than the exponential distribution, the same estimators have been found out for the exponential distribution and compared. Monte Carlo simulations and a real data analysis are performed to compare the performances of the proposed methods of estimation.     

Estimation of the distribution function under hybrid ranked set sampling

Recently, a hybrid ranked set sampling (HRSS) scheme has been proposed in the literature. The HRSS scheme encompasses several existing ranked set sampling (RSS) schemes, and it is a cost-effective alternative to the classical RSS and double RSS schemes. In this paper, we propose an improved estimator for estimating the cumulative distribution function (CDF) using HRSS. It is shown, both theoretically and numerically, that the CDF estimator under HRSS scheme is unbiased and its variance is always less than the variance of the CDF estimator with simple random sampling (SRS). An unbiased estimator of the variance of CDF estimator using HRSS is also derived. Using Monte Carlo simulations, we also study the performances of the proposed and existing CDF estimators under both perfect and imperfect rankings. It turns out that the proposed CDF estimator is by far a superior alternative to the existing CDF estimators with SRS, RSS and L-RSS schemes. For a practical application, a real data set is considered on the bilirubin level of babies in neonatal intensive care.     

Approximating Symmetric Distributions via Sampling and Coefficient of Variation

Abstract

The Coefficient of Variation is one of the most commonly used statistical tool across various scientific fields. This paper proposes a use of the Coefficient of Variation, obtained by Sampling, to define the polynomial probability density function (pdf) of a continuous and symmetric random variable on the interval [a, b]. The basic idea behind the first proposed algorithm is the transformation of the interval from [a, b] to [0, b-a]. The chi-square goodness-of-fit test is used to compare the proposed (observed) sample distribution with the expected probability distribution. The experimental results show that the collected data are approximated by the proposed pdf. The second algorithm proposes a new method to get a fast estimate for the degree of the polynomial pdf when the random variable is normally distributed. Using the known percentages of values that lie within one, two and three standard deviations of the mean, respectively, the so-called three-sigma rule of thumb, we conclude that the degree of the polynomial pdf takes values between 1.8127 and 1.8642. In the case of a Laplace (μ, b) distribution, we conclude that the degree of the polynomial pdf takes values greater than 1. All calculations and graphs needed are done using statistical software R.     

Skewed Bessel function distributions with application to rainfall data

Skewed distributions have attracted significant attention in the last few years. In this article, a skewed Bessel function distribution with the probability density function (pdf) f(x)=2 g (x) G (λ x) is introduced, where g (·) and G (·) are taken, respectively, to be the (pdf) and the cumulative distribution function of the Bessel function distribution [McKay, A.T., 1932, A Bessel function distribution, Biometrica, 24, 39–44]. Several particular cases of this distribution are identified and various representations for its moments derived. Estimation procedures by the method of maximum likelihood are also derived. Finally, an application is provided to rainfall data from Orlando, Florida.     

Using multiple imputation to estimate cumulative distribution functions in longitudinal data analysis with data missing at random

In longitudinal clinical studies, after randomization at baseline, subjects are followed for a period of time for development of symptoms. The interested inference could be the mean change from baseline to a particular visit in some lab values, the proportion of responders to some threshold category at a particular visit post baseline, or the time to some important event. However, in some applications, the interest may be in estimating the cumulative distribution function (CDF) at a fixed time point post baseline. When the data are fully observed, the CDF can be estimated by the empirical CDF. When patients discontinue prematurely during the course of the study, the empirical CDF cannot be directly used. In this paper, we use multiple imputation as a way to estimate the CDF in longitudinal studies when data are missing at random. The validity of the method is assessed on the basis of the bias and the Kolmogorov–Smirnov distance. The results suggest that multiple imputation yields less bias and less variability than the often used last observation carried forward method. Copyright © 2013 John Wiley & Sons, Ltd.     

A Dimensional CLT for Non-Central Wilks' Lambda in Multivariate Analysis

Abstract.  We consider the non-central distribution of the classical Wilks' lambda statistic for testing the general linear hypothesis in MANOVA. We prove that as the dimension of the observation vector goes to infinity, Wilks' lambda obeys a central limit theorem under simple growth conditions on the non-centrality matrix. In one case we also prove a stronger result: the saddlepoint cumulative distribution function (CDF) approximation for the standardized version of Wilks' lambda converges uniformly on compact sets to the standard normal CDF. These theoretical results go some way towards explaining why saddlepoint approximations to the distribution of Wilks' lambda retain excellent accuracy in high-dimensional cases.     

Further Results on Order Statistics Generated by Two Simulation Methods

Uniform order statistics generated by two simulation methods are compared by means of Pitman’s measure of closeness. This measure, as a probability, is shown to be asymptotically 1/2. Some results are also established for fixed points of the cumulative distribution function (CDF) for a uniform order statistic. These fixed points are important for calculations involving the joint distribution of these order statistics.     

On the Computation of the Distribution for the Analysis of Means

Abstract

We discuss the accuracy of the computation and present a fortran program to compute the cumulative distribution function (CDF) for the analysis of means (ANOM).     

The distribution of the product of two independent generalized trapezoidal random variables

ABSTRACT

In this article, we derive the probability density function (pdf) of the product of two independent generalized trapezoidal random variables having different supports, in closed form, by considering all possible cases. We also show that the results for the product of two triangular and uniform random variables follow as special cases of our main result. As an illustration, we obtain pdf of product for a suitably constrained set of parameters and plot some graphs using MATLAB, which express variation in pdf with change in different parameters of the generalized trapezoidal distribution.     

Noninformative priors for the generalized half-normal distribution

In this paper, we develop noninformative priors for the generalized half-normal distribution when scale and shape parameters are of interest, respectively. Especially, we develop the first and second order matching priors for both parameters. For the shape parameter, we reveal that the second order matching prior is a highest posterior density (HPD) matching prior and a cumulative distribution function (CDF) matching prior. In addition, it matches the alternative coverage probabilities up to the second order. For the scale parameter, we reveal that the second order matching prior is neither a HPD matching prior nor a CDF matching prior. Also, it does not match the alternative coverage probabilities up to the second order. For both parameters, we present that the one-at-a-time reference prior is a second order matching prior. However, Jeffreys’ prior is neither a first nor a second order matching prior. Methods are illustrated with both a simulation study and a real data set.     

Efficient Estimation of the PDF and the CDF of the Exponentiated Gumbel Distribution

The exponentiated Gumbel model has been shown to be useful in climate modeling including global warming problem, flood frequency analysis, offshore modeling, rainfall modeling, and wind speed modeling. Here, we consider estimation of the probability density function (PDF) and the cumulative distribution function (CDF) of the exponentiated Gumbel distribution. The following estimators are considered: uniformly minimum variance unbiased (UMVU) estimator, maximum likelihood (ML) estimator, percentile (PC) estimator, least-square (LS) estimator, and weighted least-square (WLS) estimator. Analytical expressions are derived for the bias and the mean squared error. Simulation studies and real data applications show that the ML estimator performs better than others.     

Bayesian estimation of the mean of an autoregressive process

The marginal posterior probability density function (pdf) for the mean of a stationary pth order Gaussian autoregressive process is derived using the conditional likelihood function. While the posterior pdf provides a small sample analysis, the pdf is not well known and must be analyzed numerically. This is relatively easy since it is a function of only one variable. Two sets of examples are presented. The first set involves synthetic data generated by computer, and the second set deals with energy expenditure data on a bum patient.     

A bivariate distribution with gamma and beta marginals with application to drought data

The first known bivariate distribution with gamma and beta marginals is introduced. Various representations are derived for its joint probability density function (pdf), joint cumulative distribution function (cdf), product moments, conditional pdfs, conditional cdfs, conditional moments, joint moment generating function, joint characteristic function and entropies. The method of maximum likelihood and the method of moments are used to derive the associated estimation procedures as well as the Fisher information matrix, variance–covariance matrix and the profile likelihood confidence intervals. An application to drought data from Nebraska is provided. Some other applications are also discussed. Finally, an extension of the bivariate distribution to the multivariate case is proposed.     

A weighted bootstrap approximation for comparing the error distributions in nonparametric regression

Several procedures have been proposed for testing the equality of error distributions in two or more nonparametric regression models. Here we deal with methods based on comparing estimators of the cumulative distribution function (CDF) of the errors in each population to an estimator of the common CDF under the null hypothesis. The null distribution of the associated test statistics has been approximated by means of a smooth bootstrap (SB) estimator. This paper proposes to approximate their null distribution through a weighted bootstrap. It is shown that it produces a consistent estimator. The finite sample performance of this approximation is assessed by means of a simulation study, where it is also compared to the SB. This study reveals that, from a computational point of view, the proposed approximation is more efficient than the one provided by the SB.     

Exact distributions of two non-central generalized Wilks's statistics

We derive the exact expressions of the probability density function (pdf) and the cumulative distribution function (cdf) of Wilks's likelihood ratio criterion Λ and Wilks-Lawley's statistic U in the non-central linear and the non-central planar cases. Those expressions are given in rapidly converging infinite series and can be used for numerical computation. For applications, we compute the exact power of these statistics in a multivariate analysis of variance exercise, and show by simulation the precision of our analytic formulae.