Confidence intervals and significance tests for a single trial

Confidence intervals are developed for the location parameter of a continuous, symmetric, unimodal distribution in the casev where only a single observation from the distribution is available. These intervals are similar to those given by Abbott and Rosenblatt (1963), but shorter. The result is extended to include distributions which can be standardized to have unit scale. The procedure is exemplified for the normal distribution and the power of one- and two-sided significance tests are computed under normality.     

The inverted Nadarajah–Haghighi distribution: estimation methods and applications

The inverted (or inverse) distributions are sometimes very useful to explore additional properties of the phenomenons which non-inverted distributions cannot. We introduce a new inverted model called the inverted Nadarajah–Haghighi distribution which exhibits decreasing and unimodal (right-skewed) density while the hazard rate shapes are decreasing and upside-down bathtub. Our main focus is the estimation (from both frequentist and Bayesian points of view) of the unknown parameters along with some mathematical properties of the new model. The Bayes estimators and the associated credible intervals are obtained using Markov Chain Monte Carlo techniques under squared error loss function. The gamma priors are adopted for both scale and shape parameters. The potentiality of the distribution is analysed by means of two real data sets. In fact, it is found to be superior in its ability to sufficiently model the data as compared to the inverted Weibull, inverted Rayleigh, inverted exponential, inverted gamma, inverted Lindley and inverted power Lindley models.     

Gamma-minimax results for the class of unimodal priors

Olman and Shmundak proved 1985 that in estimating a bounded normal mean under squared error loss the Bayes estimator with respect to the uniform distribution on the parameter interval is gamma-minimax when the parameter interval is sufficiently small and the class of priors consists of all symmetric and unimodal distributions. Recently, one of the authors showed that this result remains valid for quite general families of distributions which satisfy some regularity conditions. In the present paper a generalization to the class of unimodal priors with fixed mode is derived. It is proved that the Bayes estimator with respect to a suitable mixture of two uniform distributions is gamma-minimax for sufficiently small parameter intervals. To that end appropriate characterizations of a saddle point in the corresponding statistical games are established. Some results of a numerical study are presented.     

Modeling with a large class of unimodal multivariate distributions

In this paper we introduce a new class of multivariate unimodal distributions, motivated by Khintchine's representation for unimodal densities on the real line. We start by introducing a new class of unimodal distributions which can then be naturally extended to higher dimensions, using the multivariate Gaussian copula. Under both univariate and multivariate settings, we provide MCMC algorithms to perform inference about the model parameters and predictive densities. The methodology is illustrated with univariate and bivariate examples, and with variables taken from a real data set.     

Confidence curves and improved exact confidence intervals for discrete distributions

The author describes a method for improving standard “exact” confidence intervals in discrete distributions with respect to size while retaining correct level. The binomial, negative binomial, hypergeometric, and Poisson distributions are considered explicitly. Contrary to other existing methods, the author's solution possesses a natural nesting condition: if α < α', the 1 ‐ α' confidence interval is included in the 1 ‐ α interval. Nonparametric confidence intervals for a quantile are also considered.     

A likelihood ratio test to discriminate exponential–Poisson and gamma distributions

The exponential–Poisson (EP) distribution with scale and shape parameters β>0 and λ∈?, respectively, is a lifetime distribution obtained by mixing exponential and zero-truncated Poisson models. The EP distribution has been a good alternative to the gamma distribution for modelling lifetime, reliability and time intervals of successive natural disasters. Both EP and gamma distributions have some similarities and properties in common, for example, their densities may be strictly decreasing or unimodal, and their hazard rate functions may be decreasing, increasing or constant depending on their shape parameters. On the other hand, the EP distribution has several interesting applications based on stochastic representations involving maximum and minimum of iid exponential variables (with random sample size) which make it of distinguishable scientific importance from the gamma distribution. Given the similarities and different scientific relevance between these models, one question of interest is how to discriminate them. With this in mind, we propose a likelihood ratio test based on Cox's statistic to discriminate the EP and gamma distributions. The asymptotic distribution of the normalized logarithm of the ratio of the maximized likelihoods under two null hypotheses – data come from EP or gamma distributions – is provided. With this, we obtain the probabilities of correct selection. Hence, we propose to choose the model that maximizes the probability of correct selection (PCS). We also determinate the minimum sample size required to discriminate the EP and gamma distributions when the PCS and a given tolerance level based on some distance are before stated. A simulation study to evaluate the accuracy of the asymptotic probabilities of correct selection is also presented. The paper is motivated by two applications to real data sets.     

Exponentiated Sinh Cauchy Distribution with Applications

The exponentiated sinh Cauchy distribution is characterized by four parameters: location, scale, symmetry, and asymmetry. The symmetry parameter preserves the symmetry of the distribution by producing both bimodal and unimodal densities having coefficient of kurtosis values ranging from one to positive infinity. The asymmetry parameter changes the symmetry of the distribution by producing both positively and negatively skewed densities having coefficient of skewness values ranging from negative infinity to positive infinity. Bimodality, skewness, and kurtosis properties of this regular distribution are presented. In addition, relations to some well-known distributions are examined in terms of skewness and kurtosis by constructing aliases of the proposed distribution on the symmetry and asymmetry parameter plane. The maximum likelihood parameter estimation technique is discussed, and examples are provided and analyzed based on data from astronomy and medical sciences to illustrate the flexibility of the distribution for modeling bimodal and unimodal data.     

Dependent models for observations which include angular ones

This paper discusses some stochastic models for dependence of observations which include angular ones. First, we provide a theorem which constructs four-dimensional distributions with specified bivariate marginals on certain manifolds such as two tori, cylinders or discs. Some properties of the submodel of the proposed models are investigated. The theorem is also applicable to the construction of a related Markov process, models for incomplete observations, and distributions with specified marginals on the disc. Second, two maximum entropy distributions on the cylinder are discussed. The circular marginal of each model is distributed as the generalized von Mises distribution which represents a symmetric or asymmetric, unimodal or bimodal shape. The proposed cylindrical model is applied to two data sets.     

Bimodal Birnbaum–Saunders distribution with applications to non negative measurements

In this article, we introduce a new extension of the Birnbaum–Saunders (BS) distribution as a follow-up to the family of skew-flexible-normal distributions. This extension produces a family of BS distributions including densities that can be unimodal as well as bimodal. This flexibility is important in dealing with positive bimodal data, given the difficulties experienced by the use of mixtures of distributions. Some basic properties of the new distribution are studied including moments. Parameter estimation is approached by the method of moments and also by maximum likelihood, including a derivation of the Fisher information matrix. Three real data illustrations indicate satisfactory performance of the proposed model.     

On Variance Upper Bounds for Unimodal Distributions

We obtain upper bounds on the variance of discrete unimodal distributions. The alternative proofs of the corresponding bounds for the continuous unimodal distributions are also given.     

Characterization of generalized markov-polya and generalized polya-eggenberger distributions

A discrete model is considered where the original observation is subjected to partial destruction according to the Generalized Markov-Polya (GMP) damage model. A characterization of the Generalized Polya-Eggenberger distribution (GPED) is given in the context of the Rao-Rubin condition. More specifically, if the probability that an observation n of a non-negative integer valued r.v.X is reduced to an integer k during a damage, process is given by the GMPD, and if the resulting r.v.Y is such thatrit satisfies the RR-conditlon, then X has a GPED. Secondly, if N = A + B, where B is the missing part and A is the recorded part such that the conditional distribution P(A= x|N=n) is the GMPD, then the r.v.'s A and B are independent if, and only if, N has a GPED. Several other characterizations are also given for these two distributions. The results of Rao-Rubin ‘1964’, Patil-Ratnaparkhi (1977) and Consul (1975) follow as special cases.     

Prediction Intervals for an Ordered Observation from Weibull Distribution Based on Censored Samples

In this article, we provide some suitable pivotal quantities for constructing prediction intervals for the jth future ordered observation from the two-parameter Weibull distribution based on censored samples. Our method is more general in the sense that it can be applied to any data scheme. We present a simulation of our method to analyze its performance. Two illustrative examples are also included. For further study, our method is easily applied to other location and scale family distributions.     

Bivariate kurtotic distributions of garment fibre data

Summary. A bivariate and unimodal distribution is introduced to model an unconventionally distributed data set collected by the Forensic Science Service. This family of distributions allows for a different kurtosis in each orthogonal direction and has a constructive rather than probability density function definition, making conventional inference impossible. However, the construction and inference work well with a Bayesian Markov chain Monte Carlo analysis.     

Bootstrap confidence intervals of CNpk for inverse Rayleigh and log-logistic distributions

In this article bootstrap confidence intervals of process capability index as suggested by Chen and Pearn [An application of non-normal process capability indices. Qual Reliab Eng Int. 1997;13:355–360] are studied through simulation when the underlying distributions are inverse Rayleigh and log-logistic distributions. The well-known maximum likelihood estimator is used to estimate the parameter. The bootstrap confidence intervals considered in this paper consists of various confidence intervals. A Monte Carlo simulation has been used to investigate the estimated coverage probabilities and average widths of the bootstrap confidence intervals. Application examples on two distributions for process capability indices are provided for practical use.     

A 3-parameter Gompertz distribution for survival data with competing risks,with an application to breast cancer data

The cumulative incidence function is of great importance in the analysis of survival data when competing risks are present. Parametric modeling of such functions, which are by nature improper, suggests the use of improper distributions. One frequently used improper distribution is that of Gompertz, which captures only monotone hazard shapes. In some applications, however, subdistribution hazard estimates have been observed with unimodal shapes. An extension to the Gompertz distribution is presented which can capture unimodal as well as monotone hazard shapes. Important properties of the proposed distribution are discussed, and the proposed distribution is used to analyze survival data from a breast cancer clinical trial.     

Multiplicative Strong Unimodality

Multiplicative strong unimodality is defined as the preservation of unimodality in products of independent random variables. An Ibragimov type theorem is proved. As an application, preserving unimodality for scale mixtures of gamma distributions is examined. It is also shown that multiplicative strong unimodal probability measures on R appear as images, by the exponential map, of classical strong unimodal ones. The connection to the star order is also established.     

Shortest prediction intervals for the birnbaum-saunders distribution

Shortest prediction intervals for a future observation from the Birnbaum-Saunders distribution are obtained from both frequentist and Bayesian perspectives. Comparisons are made with alternative intervals obtained via inversion. Monte Carlo simulations are performed to assess the approximate intervals.     

THE WRAPPED t FAMILY OF CIRCULAR DISTRIBUTIONS

This paper considers the three‐parameter family of symmetric unimodal distributions obtained by wrapping the location‐scale extension of Student's t distribution onto the unit circle. The family contains the wrapped normal and wrapped Cauchy distributions as special cases, and can be used to closely approximate the von Mises distribution. In general, the density of the family can only be represented in terms of an infinite summation, but its trigonometric moments are relatively simple expressions involving modified Bessel functions. Point estimation of the parameters is considered, and likelihood‐based methods are used to fit the family of distributions in an illustrative analysis of cross‐bed measurements. The use of the family as a means of approximating the von Mises distribution is investigated in detail, and new efficient algorithms are proposed for the generation of approximate pseudo‐random von Mises variates.     

Bimatrix variate gamma-beta distributions

Two bimatrix distributions with beta and gamma marginals are introduced. Various properties (including product moments of determinants and traces, entropies, marginal distributions) are derived. Parameter estimation by the method of maximum likelihood is discussed. The performance and efficiencies of the maximum likelihood estimators and the associated confidence intervals are assessed by simulation. The efficiencies are compared versus those for the maximum likelihood estimators and the associated confidence intervals based on matrix variate gamma distributions. A discussion of possible applications of the bimatrix distributions is given.     

ESTIMATING THE ANGLE BETWEEN THE MEAN DIRECTIONS OF TWO SPHERICAL DISTRIBUTIONS

This paper considers the problem of calculating a confidence interval for the angular difference between the mean directions of two spherical random variables with rotationally symmetric unimodal distributions. For large sample sizes, it is shown that the asymptotic distribution of 1 – cos α, where α is the sample angular difference, is approximately exponential if the true difference is zero, and approximately normal for a ‘large’ true difference; a scaled beta approximation is determined for the general case. For small sample sizes, a bootstrap approach is recommended. The results are applied to two sets of palaeomagnetic data.