Order statistics from discrete distributions

Several recurrence relations and identities available for single and product moments of order1 statistics in a sample size n from an arbitrary continuous distribution are extended for the discrete case,, Making use of these recurrence relations it is shown that it is sufficient to evaluate just two single moments and (n-l)/2 product moments when n is odd and two single moments and {n-2)/2 product moments when n is even, in order to evaluate the first, second and product moments of order statistics in a sample of size n drawn from an arbitrary discrete distribution, given these moments in samples of sizes n-1 and less.. A series representation for the product moments of order statistics is derived.. Besides enabling us to obtain an exact and explicit expression for the product moments of order statistics from the geometric distribution, it. makes the computation of the product moments of order statistics from other discrete distributions easy too.     

Moments of Order Statistics—From Even to Odd Sample Sizes

Abstract

A method is demonstrated to compute the complete set of first moments of order statistics for an arbitrary distribution, given only the first moments of the maximal order statistics either for all even sample sizes, or for all odd samples sizes.     

Minimax optimality of the Shiryayev–Roberts change-point detection rule

The result of Pollak [1985. Optimal detection of a change in distribution. Ann. Statist. 13, 206–227] proving the asymptotic optimality in sequential change-point detection of a suitable Shirayayev–Roberts stopping rule up to terms that vanish in the limit is generalized from the case of two completely specified distributions to that of a composite alternative hypothesis in a multidimensional exponential family. An explicit asymptotic lower bound on the expected Kullback–Leibler information required to detect a change-point is derived and is shown to be attained by a Shirayayev–Roberts stopping rule.     

Further results on the trace of a noncentral wishart matrix

In an earlier article Mathai (1980) has given compact representations for the moments and cumulants of the trace of a noncentral Wishart matrix. He has also shown that (trA-ntr;∑)/(2ntri∑²)172. is asymptotically standard normal where A is a noncentral wishart matrix with n degrees of freedom and covariance matrix [0, In the present article explicit expressions for the exact density of the trace are given in terms of confluent hypergeometric functions and in terms of zonal polynomials for the general case and as finite sums when the sample size is odd. As a consequence of some of these representations some summation formulae for zonal polynomials are also given     

On the distribution of the maximum of bivariate normal random variables with general means and variances

A great amount of effort has been devoted to achieving exact expressions for moments of order statistics of independent normal random variables, as well as the dependent case with the same correlation coefficients, means and variances. It does not seem as if there are handy formulae for the order statistics of even the simple bivariate normal random variables when the means and variances are allowed to be different. In this paper we give an explicit formula for the Lanl ace-Stielties Transform of the maximum of bivariate normal random variables by which we obtain formulae for the first two moments in the standard way.     

Inference Problems Involving Moment Determinacy of Distributions

We describe diverse stochastic inference problems whose solution essentially depends on the moment determinacy of some distributions involved. For a variety of stochastic models we ask questions such as “how to identify a distribution if knowing its moments?” “how asymmetric can be a distribution with zero odd order moments?” “is any mixture model identifiable?” For specific models we provide answers, motivating arguments, and illustrations. Some challenging open questions are outlined.     

Sojourn time in an M/M/1 processor sharing queue with permanent customers

We explicitly compute the sojourn time distribution of an arbitrary customer in an M/M/1 processor sharing (PS) queue with permanent customers. We notably exhibit the orthogonal structure associated with this queuing system and we show how sieved Pollaczek polynomials and their associated orthogonality measure can be used to obtain an explicit representation for the complementary cumulative distribution function of the sojourn time of a customer. This explicit formula subsequently allows us to compute the two first moments of this random variable and to study the asymptotic behavior of its distribution. The most salient result is that the decay rate depends on the load of the system and the number K of permanent customers. When the load is above a certain threshold depending on K, the decay rate is identical to that of a regular M/M/1 PS queue.     

Uniform robustness against nonnormality of the t and f tests

The size of the two-sample t test is generally thought to be robust against nonnormal distributions if the sample sizes are large. This belief is based on central limit theory, and asymptotic expansions of the moments of the t statistic suggest that robustness may be improved for moderate sample sizes if the variance, skewness, and kurtosis of the distributions are matched, particularly if the sample sizes are also equal.

It is shown that asymptotic arguments such as these can be misleading and that, in fact, the size of the t test can be as large as unity if the distributions are allowed to be completely arbitrary. Restricting the distributions to be identical or symmetric (but otherwise arbitrary) does not guarantee that the size can be controlled either, but controlling the tail-heaviness of the distributions does. The last result is proved more generally for the k-sample F test.     

Explicit formulas for the cumulants and the vector-valued odd moments of the multivariate linearly skewed elliptical distributions

In this letter explicit expressions are derived for the cumulants and the vector-valued odd moments of the multivariate linearly skewed elliptical family of distributions. The general calculations of such moments are described by multivariate integrals which complicate the calculations. We show how such multivariate computations can be projected into a univariate framework, which extremely simplifies the computations.     

Another Proof of Borel's Strong Law of Large Numbers

Let X be a continuous nonnegative random variable with finite first and second moments and a continuous pdf that is positive on the interior of its support. A nonzero limiting density at the origin and a coefficient of variation (CV) greater than 1 are shown to be sufficient conditions for the distribution truncated below at t > 0 to have a variance greater than the variance of the full distribution. Distributions that satisfy these conditions include those with decreasing hazard rates (e.g., the gamma and Weibull distributions with shape parameters less than 1) and the beta distribution with parameter values p and q for which q > p(p + q + 1). The bound T for which truncation at 0 < t < T increases the variance relative to the full distribution is shown to be greater than the (1 — 1/CV)th percentile of the full distribution.     

A generalized skew two-piece skew-elliptical distribution

We present a new generalized family of skew two-piece skew-elliptical (GSTPSE) models and derive some its statistical properties. It is shown that the new family of distribution may be written as a mixture of generalized skew elliptical distributions. Also, a new representation theorem for a special case of GSTPSE-distribution is given. Next, we will focus on t kernel density and prove that it is a scale mixture of the generalized skew two-piece skew normal distributions. An explicit expression for the central moments as well as a recurrence relations for its cumulative distribution function and density are obtained. Since, this special case is a uni-/bimodal distribution, a sufficient condition for each cases is given. A real data set on heights of Australian females athletes is analysed. Finally, some concluding remarks and open problems are discussed.     

AN INFINITE-PHASE QUASI-BIRTH-AND-DEATH MODEL FOR THE NON-PREEMPTIVE PRIORITY M/PH/1 QUEUE

This paper considers a single server queue that handles arrivals from N classes of customers on a non-preemptive priority basis. Each of the N classes of customers features arrivals from a Poisson process at rate λ _i and class-dependent phase type service. To analyze the queue length and waiting time processes of this queue, we derive a matrix geometric solution for the stationary distribution of the underlying Markov chain. A defining characteristic of the paper is the fact that the number of distinct states represented within the sub-level is countably infinite, rather than finite as is usually assumed. Among the results we obtain in the two-priority case are tractable algorithms for the computation of both the joint distribution for the number of customers present and the marginal distribution of low-priority customers, and an explicit solution for the marginal distribution of the number of high-priority customers. This explicit solution can be expressed completely in terms of the arrival rates and parameters of the two service time distributions. These results are followed by algorithms for the stationary waiting time distributions for high- and low-priority customers. We then address the case of an arbitrary number of priority classes, which we solve by relating it to an equivalent three-priority queue. Numerical examples are also presented.     

A study of moments and likelihood estimators of the odd Weibull distribution

The odd Weibull distribution is a three-parameter generalization of the Weibull and the inverse Weibull distributions having rich density and hazard shapes for modeling lifetime data. This paper explored the odd Weibull parameter regions having finite moments and examined the relation to some well-known distributions based on skewness and kurtosis functions. The existence of maximum likelihood estimators have shown with complete data for any sample size. The proof for the uniqueness of these estimators is given only when the absolute value of the second shape parameter is between zero and one. Furthermore, elements of the Fisher information matrix are obtained based on complete data using a single integral representation which have shown to exist for any parameter values. The performance of the odd Weibull distribution over various density and hazard shapes is compared with generalized gamma distribution using two different test statistics. Finally, analysis of two data sets has been performed for illustrative purposes.     

The Concept of Majorization in Experimental Design

Row and column designs are examined and the concept of majorization is applied to find optimal designs, when the observations are either independent or dependent. The dependence follows a first-order autoregression with parameter α. The case of two treatments is examined and the universally optimal or Φ-optimal designs are given, for different values of a, when the number of experimental units is even or odd. A filtering procedure is followed to reduce the number of competing designs.     

Risk comparison of the Stein-rule estimator in a linear regression model with omitted relevant regressors and multivariate<Emphasis Type="Italic">t</Emphasis> errors under the Pitman nearness criterion

In this paper we consider a linear regression model with omitted relevant regressors and multivariatet error terms. The explicit formula for the Pitman nearness criterion of the Stein-rule (SR) estimator relative to the ordinary least squares (OLS) estimator is derived. It is shown numerically that the dominance of the SR estimator over the OLS estimator under the Pitman nearness criterion can be extended to the case of the multivariatet error distribution when the specification error is not severe. It is also shown that the dominance of the SR estimator over the OLS estimator cannot be extended to the case of the multivariatet error distribution when the specification error is severe. This research is partially supported by the Grants-in-Aid for 21st Century COE program.     

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.     

On some study of skew-t distributions

Abstract

In this note, through ratio of independent random variables, new families of univariate and bivariate skew-t distributions are introduced. Probability density function for each skew-t distribution will be given. We also derive explicit forms of moments of the univariate skew-t distribution and recurrence relations for its cumulative distribution function. Finally we illustrate the flexibility of this class of distributions with applications to a simulated data and the volcanos heights data.     

Useful moment and CDF formulations for the COM–Poisson distribution

The Poisson distribution is as important for discrete events as the normal distribution is to large sample data. In this note, we discuss a generalized Poisson distribution recently introduced in the statistics literature. We derive—for the first time—exact and explicit expressions for its moments and the cumulative distribution function for the case of over-dispersion. Computational issues are discussed to show the real value of these expressions.     

On moments of dual generalized order statistics from Topp-Leone distribution

In this paper, we have derived exact and explicit expressions for the ratio and inverse moments of dual generalized order statistics from Topp-Leone distribution. This result includes the single and product moments of order statistics and lower records . Further, based on n dual generalized order statistics, we have deduced the expression for Maximum likelihood estimator (MLE) and Uniformly minimum variance unbiased estimator (UMVUE) for the shape parameter of Topp-Leone distribution. Finally, based on order statistics and lower records, a simulation study is being carried out to check the efficiency of these estimators.     

Parameter Estimation for the Discrete Stable Family

The discrete stable family constitutes an interesting two-parameter model of distributions on the non-negative integers with a Paretian tail. The practical use of the discrete stable distribution is inhibited by the lack of an explicit expression for its probability function. Moreover, the distribution does not possess moments of any order. Therefore, the usual tools—such as the maximum-likelihood method or even the moment method—are not feasible for parameter estimation. However, the probability generating function of the discrete stable distribution is available in a simple form. Hence, we initially explore the application of some existing estimation procedures based on the empirical probability generating function. Subsequently, we propose a new estimation method by minimizing a suitable weighted L ²-distance between the empirical and the theoretical probability generating functions. In addition, we provide a goodness-of-fit statistic based on the same distance.