Star-shaped distributions and their generalizations

Elliptically contoured distributions can be considered to be the distributions for which the contours of the density functions are proportional ellipsoids. We generalize elliptically contoured densities to “star-shaped distributions” with concentric star-shaped contours and show that many results in the former case continue to hold in the more general case. We develop a general theory in the framework of abstract group invariance so that the results can be applied to other cases as well, especially those involving random matrices.     

On smoothness of Tukey depth contours

The smoothness of Tukey depth contours is a regularity condition often encountered in asymptotic theory, among others. This condition ensures that the Tukey depth fully characterizes the underlying multivariate probability distribution. In this paper we demonstrate that this regularity condition is rarely satisfied. It is shown that even well-behaved probability distributions with symmetrical, smooth and (strictly) quasi-concave densities may have non-smooth Tukey depth contours, and that the smoothness behaviour of depth contours is fairly unpredictable.     

Distributional results for special cases of the jolly-seber model

The multinomial-binomial approach to the Jolly-Seber capture- recapture model is used as a basis to derive explicit probability distributions for special cases of the Jolly-Seber model:no recruitment, or no mortality. Also given are the residual distributions that allow tests of these restricted models compared to the general Jolly-Seber model. Losses on capture are allowed. The special case distribution is also derived for no recruitment and no mortality, but allowing losses on capture; this is a generalized version of Darroch's closed capture-recapture model. Here, however, it was not possible to obtain a closed form residual distribution.     

Obtaining a Trapezoidal Distribution

Given a most believed value for a quantity together with upper and lower possible deviations from that value, a rectangular distribution might be used to represent state-of-knowledge about the quantity. If the deviations are themselves known by probability distributions, and the value conditioned on the deviations is rectangular, then the marginal distribution of the value is determined by the distributions of the deviations. Here we show under quite general conditions that conversely, given the marginal distribution, the distributions of the deviations are uniquely determined. The case in which the marginal distribution is trapezoidal is studied in some detail.     

A generalized binomial distribution

A new generalization of the binomial distribution is introduced that allows dependence between trials, nonconstant probabilities of success from trial to trial, and which contains the usual binomial distribution as a special case. Along with the number of trials and an initial probability of ‘success’, an additional parameter that controls the degree of correlation between trials is introduced. The resulting class of distributions includes the binomial, unirnodal distributions, and bimodal distributions. Formulas for the moments, mean, and variance of this distribution are given along with a method for fitting the distribution to sample data.     

Weighted bivariate logarithmic series distributions

In this paper, we are interested in the weighted distributions of a bivariate three parameter logarithmic series distribution studied by Kocherlakota and Kocherlakota (1990). The weighted versions of the model are derived with weight W(x,y) = x^[r] y^[s]. Explicit expressions for the probability mass function and probability generating functions are derived in the case r = s = l. The marginal and conditional distributions are derived in the general case. The maximum likelihood estimation of the parameters, in both two parameter and three parameter cases, is studied. A procedure for computer generation of bivariate data from a discrete distribution is described. This enables us to present two examples, in order to illustrate the methods developed, for finding the maximum likelihood estimates.     

Joint distribution of simultaneous exposures to several carcinogens in a case–control study: sample size determination

Consider a population of individuals who are free of a disease under study, and who are exposed simultaneously at random exposure levels, say X,Y,Z,… to several risk factors which are suspected to cause the disease in the populationm. At any specified levels X=x, Y=y, Z=z, …, the incidence rate of the disease in the population ot risk is given by the exposure–response relationship r(x,y,z,…) = P(disease|x,y,z,…). The present paper examines the relationship between the joint distribution of the exposure variables X,Y,Z, … in the population at risk and the joint distribution of the exposure variables U,V,W,… among cases under the linear and the exponential risk models. It is proven that under the exponential risk model, these two joint distributions belong to the same family of multivariate probability distributions, possibly with different parameters values. For example, if the exposure variables in the population at risk have jointly a multivariate normal distribution, so do the exposure variables among cases; if the former variables have jointly a multinomial distribution, so do the latter. More generally, it is demonstrated that if the joint distribution of the exposure variables in the population at risk belongs to the exponential family of multivariate probability distributions, so does the joint distribution of exposure variables among cases. If the epidemiologist can specify the differnce among the mean exposure levels in the case and control groups which are considered to be clinically or etiologically important in the study, the results of the present paper may be used to make sample size determinations for the case–control study, corresponding to specified protection levels, i.e., size α and 1–β of a statistical test. The multivariate normal, the multinomial, the negative multinomial and Fisher's multivariate logarithmic series exposure distributions are used to illustrate our results.     

Prediction intervals for generalized-order statistics with random sample size

The problem of predicting future generalized-order statistics, by assuming the future sample size is a random variable, is discussed. A general expression for the coverage probability of the prediction intervals is derived. Since k-records and progressively type-II censored-order statistics are contained in the model of generalized-order statistics, the corresponding results for them can be deduced as special cases. When the future sample size has degenerate, binomial, Poisson and geometric distributions, numerical computations are given. The procedure for finding an optimal prediction interval is presented for each case. Finally, we apply our results to a real data set in life testing given in Lee and Wang [Statistical methods for survival data analysis. Hoboken, NJ: John Wiley and Sons; 2003, p. 58, Table 3.4] for illustrative the proposed procedure in this paper.     

New class of location-parameter discrete probability distributions and their characterizations

A new class of location-parameter discrete probability distributions (LDPD) has been defined where the population mean is the location parameter. It has been shown that some single parameter discrete distributions do not belong to this class and all discrete probability distributions belonging to this class can be characterized by their variances only. Expressions are given for the first four central moments and a recurrence formula for higher central moments has been obtained. Eight theorems are given to characterize the various distributions in the LDPD class.     

Simultaneous closeness among order statistics to population quantiles

We derive expressions for the probability that an individual order statistic is closest to the target parameter among the order statistics from a complete random sample. Results are given for random variables with bounded and complete support. We then apply these general results to location-scale parameter families of distributions with specific applications to estimation of percentiles. In this case, simultaneous-closeness probabilities depend upon the parameters through the value of p in the percentile and the sample size, n. Results are finally illustrated with the estimation of percentiles for normal and exponential distributions.     

A class of weighted multivariate elliptical models useful for robust analysis of nonnormal and bimodal data

A class of weighted elliptical models useful for analyzing nonnormal and bimodal multivariate data is introduced. It is obtained from the marginal distribution of a centrally truncated multivariate elliptical distribution. As a special case, a finite mixture of weighted multinormal distribution is examined in detail, establishing connections with the multinormal and the finite mixture of multinormal. The special class of distributions is studied from several aspects such as weighting of probability density functions, association with centrally truncated distributions, and a finite scale mixture scheme. The relationships among these aspects are given, and various properties of the class are also discussed. For the inference of the class, an MCMC procedure and its numerical example are provided.     

A note on bounding monte carlo variances michael lavine

A density bounded class P of probability distributions on a space χ is the set of all probability distributions corresponding to probability densities bounded below by a given subprob-ability density and bounded above by a given superprobability density. Density bounded classes arise in robust Bayesian analysis (Lavine 1991) and also in Monte Carlo integration (Fishman Granovsky and Rubin 1989). Finding upper and lower bounds on the variance over all p? P allows one to bound the Monte Carlo variance. Fishman Granovsky and Rubin (1989) find bounds on the variance over all p ? P and also find the densities in P achieving those bounds in the case where χ is discrete; that is, where P is actually a set of probability mass functions. This article generalizes their result by showing how to bound the variance and find the densities achieving the bounds when χ is continuous.     

Some aspects of decision making under uncertainty

Two discrete-time insurance models are studied in the framework of cost approach. The models being non-deterministic one deals with decision making under uncertainty. Three different situations are investigated: (1) underlying processes are stochastic however their probability distributions are given; (2) information concerning the distribution laws is incomplete; (3) nothing is known about the processes under consideration. Mathematical methods useful for establishing the (asymptotically) optimal control are demonstrated in each case. Algorithms for calculation of critical levels are proposed. Numerical results are presented as well.     

Some Geometric Methods for Constructing Decision Criteria Based On Two-Dimensional Parameters

This paper reviews two types of geometric methods proposed in recent years for defining statistical decision rules based on 2-dimensional parameters that characterize treatment effect in a medical setting. A common example is that of making decisions, such as comparing treatments or selecting a best dose, based on both the probability of efficacy and the probability toxicity. In most applications, the 2-dimensional parameter is defined in terms of a model parameter of higher dimension including effects of treatment and possibly covariates. Each method uses a geometric construct in the 2-dimensional parameter space based on a set of elicited parameter pairs as a basis for defining decision rules. The first construct is a family of contours that partitions the parameter space, with the contours constructed so that all parameter pairs on a given contour are equally desirable. The partition is used to define statistical decision rules that discriminate between parameter pairs in term of their desirabilities. The second construct is a convex 2-dimensional set of desirable parameter pairs, with decisions based on posterior probabilities of this set for given combinations of treatments and covariates under a Bayesian formulation. A general framework for all of these methods is provided, and each method is illustrated by one or more applications.     

Study of the quality of several asymptotic and near-exact approximations based on moments for the distribution of the Wilks Lambda statistic

In this paper a measure of proximity of distributions, when moments are known, is proposed. Based on cases where the exact distribution is known, evidence is given that the proposed measure is accurate to evaluate the proximity of quantiles (exact vs. approximated). The measure may be applied to compare asymptotic and near-exact approximations to distributions, in situations where although being known the exact moments, the exact distribution is not known or the expression for its probability density function is not known or too complicated to handle. In this paper the measure is applied to compare newly proposed asymptotic and near-exact approximations to the distribution of the Wilks Lambda statistic when both groups of variables have an odd number of variables. This measure is also applied to the study of several cases of telescopic near-exact approximations to the exact distribution of the Wilks Lambda statistic based on mixtures of generalized near-integer gamma distributions.     

A New Markov Binomial Distribution

In this article, we introduce a two-state homogeneous Markov chain and define a geometric distribution related to this Markov chain. We define also the negative binomial distribution similar to the classical case and call it NB related to interrupted Markov chain. The new binomial distribution is related to the interrupted Markov chain. Some characterization properties of the geometric distributions are given. Recursion formulas and probability mass functions for the NB distribution and the new binomial distribution are derived.     

Parameter Estimation and Censored Lives

The usual likelihood used to make inferences about parameters of lifetime distributions when lives may be censored consists of the product of probability densities for completed lives and reliabilities for the censored lives. This works well in many situations. We consider two special cases of a model which assumes a Poisson birth process for the items: in one case the usual method works well, but in the other case it does not, and better estimates are found.     

Tail exponent estimation via broadband log density-quantile regression

Heavy tail probability distributions are important in many scientific disciplines such as hydrology, geology, and physics and therefore feature heavily in statistical practice. Rather than specifying a family of heavy-tailed distributions for a given application, it is more common to use a nonparametric approach, where the distributions are classified according to the tail behavior. Through the use of the logarithm of Parzen's density-quantile function, this work proposes a consistent, flexible estimator of the tail exponent. The approach we develop is based on a Fourier series estimator and allows for separate estimates of the left and right tail exponents. The theoretical properties for the tail exponent estimator are determined, and we also provide some results of independent interest that may be used to establish weak convergence of stochastic processes. We assess the practical performance of the method by exploring its finite sample properties in simulation studies. The overall performance is competitive with classical tail index estimators, and, in contrast, with these our method obtains somewhat better results in the case of lighter heavy-tailed distributions.     

A comparison of two approaches for power and sample size calculations in logistic regression models

Whittemore (1981) proposed an approach for calculating the sample size needed to test hypotheses with specified significance and power against a given alternative for logistic regression with small response probability. Based on the distribution of covariate, which could be either discrete or continuous, this approach first provides a simple closed-form approximation to the asymptotic covariance matrix of the maximum likelihood estimates, and then uses it to calculate the sample size needed to test a hypothesis about the parameter. Self et al. (1992) described a general approach for power and sample size calculations within the framework of generalized linear models, which include logistic regression as a special case. Their approach is based on an approximation to the distribution of the likelihood ratio statistic. Unlike the Whittemore approach, their approach is not limited to situations of small response probability. However, it is restricted to models with a finite number of covariate configurations. This study compares these two approaches to see how accurate they would be for the calculations of power and sample size in logistic regression models with various response probabilities and covariate distributions. The results indicate that the Whittemore approach has a slight advantage in achieving the nominal power only for one case with small response probability. It is outperformed for all other cases with larger response probabilities. In general, the approach proposed in Self et al. (1992) is recommended for all values of the response probability. However, its extension for logistic regression models with an infinite number of covariate configurations involves an arbitrary decision for categorization and leads to a discrete approximation. As shown in this paper, the examined discrete approximations appear to be sufficiently accurate for practical purpose.     

Exact and Approximate Distributions of the Maximum Likelihood Estimate in a Simple Factor Analysis Model

We study a factor analysis model with two normally distributed observations and one factor. In the case when the errors have equal variance, the maximum likelihood estimate of the factor loading is given in closed form. Exact and approximate distributions of the maximum likelihood estimate are considered. The exact distribution function is given in a complex form that involves the incomplete Beta function. Approximations to the distribution function are given for the cases of large sample sizes and small error variances. The accuracy of the approximations is discussed