A prior for the variance in hierarchical models

The choice of prior distributions for the variances can be important and quite difficult in Bayesian hierarchical and variance component models. For situations where little prior information is available, a ‘nonin-formative’ type prior is usually chosen. ‘Noninformative’ priors have been discussed by many authors and used in many contexts. However, care must be taken using these prior distributions as many are improper and thus, can lead to improper posterior distributions. Additionally, in small samples, these priors can be ‘informative’. In this paper, we investigate a proper ‘vague’ prior, the uniform shrinkage prior (Strawder-man 1971; Christiansen & Morris 1997). We discuss its properties and show how posterior distributions for common hierarchical models using this prior lead to proper posterior distributions. We also illustrate the attractive frequentist properties of this prior for a normal hierarchical model including testing and estimation. To conclude, we generalize this prior to the multivariate situation of a covariance matrix.     

Relation of modified power series distributions to lagrangian probability distributions

The class of Lagrangian probability distributions ‘LPD’, given by the expansion of a probability generating function f‘t’ under the transformation u = t/g‘t’ where g ‘t’ is also a p.g.f., has been substantially widened by removing the restriction that the defining functions g ‘t’ and f‘t’ be probability generating functions. The class of modified power series distributions defined by Gupta ‘1974’ has been shown to be a sub-class of the wider class of LPDs     

When there are outliers in the carriers:the univariate case

Robust regression estimators studied to date are robust against non-normal distributions of the errors only If the carriers ‘Independent variables’ do not also contain outliers. Several alternative estimators that are robust even 1f there are outliers in the carriers are studied. Two estimators seem to be preferable, but even these can be very Inefficient ‘relative to least squares’ If the errors are normally distributed.     

Distributional properties of a model for the spread or drug abuse

We examine the properties of the distribution of the number of drug abusers in a previously free community assuming that initiators enter the community during a ‘latent’ period in which they randomly infect other members of the community, and that during the subsequent ‘control’ period the spread of abuse follows a linear birth and death process. The form of distribution is shown to be unaltered by a series of steady changes in the birth and death parameters. The distribution can be regarded as the convolution of three distributions:pseudo-binomial or binomial, negative binomial, and Polya-Aeppli. Special cases include the Laguerre series distribution.     

Household budget-share distributions and welfare implications: an application of multivariate distributional statistics

In this paper the consequences of considering the household ‘food share’ distribution as a welfare measure, in isolation from the joint distribution of itemized budget shares, is examined through the unconditional and conditional distribution of ‘food share’ both parametrically and nonparametrically. The parametric framework uses Dirichlet and Beta distributions, while the nonparametric framework uses kernel smoothing methods. The analysis, in a three commodity setup (‘food’, ‘durables’, ‘others’), based on household level rural data for West Bengal, India, for the year 2009–2010 shows significant underrepresentation of households by the conventional unconditional ‘food share’ distribution in the higher range of food budget shares that correspond to the lower end of the income profile. This may have serious consequences for welfare measurement.     

A comparison of analysis procedures for correlated binary data in dedicated multi‐rater imaging trials

In this paper, three analysis procedures for repeated correlated binary data with no a priori ordering of the measurements are described and subsequently investigated. Examples for correlated binary data could be the binary assessments of subjects obtained by several raters in the framework of a clinical trial. This topic is especially of relevance when success criteria have to be defined for dedicated imaging trials involving several raters conducted for regulatory purposes. First, an analytical result on the expectation of the ‘Majority rater’ is presented when only the marginal distributions of the single raters are given. The paper provides a simulation study where all three analysis procedures are compared for a particular setting. It turns out that in many cases, ‘Average rater’ is associated with a gain in power. Settings were identified where ‘Majority significant’ has favorable properties. ‘Majority rater’ is in many cases difficult to interpret. Copyright © 2014 John Wiley & Sons, Ltd.     

A new family of generalized distributions

Kumaraswamy [Generalized probability density-function for double-bounded random-processes, J. Hydrol. 462 (1980), pp. 79–88] introduced a distribution for double-bounded random processes with hydrological applications. For the first time, based on this distribution, we describe a new family of generalized distributions (denoted with the prefix ‘Kw’) to extend the normal, Weibull, gamma, Gumbel, inverse Gaussian distributions, among several well-known distributions. Some special distributions in the new family such as the Kw-normal, Kw-Weibull, Kw-gamma, Kw-Gumbel and Kw-inverse Gaussian distribution are discussed. We express the ordinary moments of any Kw generalized distribution as linear functions of probability weighted moments (PWMs) of the parent distribution. We also obtain the ordinary moments of order statistics as functions of PWMs of the baseline distribution. We use the method of maximum likelihood to fit the distributions in the new class and illustrate the potentiality of the new model with an application to real data.     

A note on the unimddality of discrete distributions

It is noted that the unimadility property is very inportant and necessary in many probabilistic-statistical models. In this paper, we consider the definition of discrete uniirodality such that the mode may be unique integer or a sequence of consecutive integers. It will be shown that the necessary and sufficient condition for the discrete distribution to be uninnodal about a can be given through some canonical representation of its characteristic function ‘ch.f.’ Further, characterization results for some well-known distributions are established.     

Dissimilarity measures for time trajectories

In this paper, by considering one of the possible geometric representations of time arrays in the ‘object space’ (space of the units), we analyze different dissimilarity measures between multivariate time trajectories of the units, which are classified, systematically, in various approaches, by taking into account their features. In particular, we define three classes of dissimilarities: the ‘geometric’ class, in which dissimilarities are built according to the geometrical features of the trajectories (instantaneous position, slope of the inter-temporal segments (velocity), concavity and convexity of each pair of inter-temporal segments (acceleration), polygonal line (shape), portion of area between each pair of trajectories); the ‘correlative’ class, in which dissimilarity measures that analyze the autocorrelation and cross-correlation functions of the univariate components of each multivariate trajectory are classified; the ‘structural’ class, containing dissimilarities which consider the structural aspects of the trajectories, such as the linear or polynomial trends and the seasonality of each univariate component. An empirical comparison is also included.     

Tests for bivariate mean residual life

In this paper we have developed tests for bivariate exponentiaIity against the ‘bivariate decreasing mean residual life (BDMRL)’ and ‘bivariate new better than used in expectation (BNBUE)’ classes of non-exponentia1 probability distributions. We have also obtained a large-sample approximation to make the test readily applicable.     

Stochastic ordering and a class of multivariate new better than used distributions

A new characterization for the univariate class of new better than used ‘NBU’ distributions in terms of stochastic ordering is introduced. A multivariate version of this characterization is then used to define a multivariate class of NBU distributions. Basic properties of this class are derived. Comparisons and relationships of this new class with earlier classes are developed. Two multivariate new worse than used (NWU) classes of life distributions are defined and compared and their basic properties are studied.     

Residual life estimation under time-varying conditions based on a Wiener process

Residual life (RL) estimation plays an important role in prognostics and health management. In operating conditions, components usually experience stresses continuously varying over time, which have an impact on the degradation processes. This paper investigates a Wiener process model to track and predict the RL under time-varying conditions. The item-to-item variation is captured by the drift parameter and the degradation characteristic of the whole population is described by the diffusion parameter. The bootstrap method and Bayesian theorem are employed to estimate and update the distribution parameters of ‘a’ and ‘b’, which are the coefficients of the linear drifting process in the degradation model. Once new degradation information becomes available, the RL distributions considering the future operating condition are derived. The proposed method is tested on Lithium-ion battery devices under three levels of charging/discharging rates. The results are further validated by a simulation method.     

On the Breakdown Properties of Some Multivariate M‐Functionals*

Abstract. For probability distributions on ? ^q, a detailed study of the breakdown properties of some multivariate M‐functionals related to Tyler's [Ann. Statist. 15 (1987) 234] ‘distribution‐free’ M‐functional of scatter is given. These include a symmetrized version of Tyler's M‐functional of scatter, and the multivariate t M‐functionals of location and scatter. It is shown that for ‘smooth’ distributions, the (contamination) breakdown point of Tyler's M‐functional of scatter and of its symmetrized version are 1/q and , respectively. For the multivariate t M‐functional which arises from the maximum likelihood estimate for the parameters of an elliptical t distribution on ν ≥ 1 degrees of freedom the breakdown point at smooth distributions is 1/( q + ν). Breakdown points are also obtained for general distributions, including empirical distributions. Finally, the sources of breakdown are investigated. It turns out that breakdown can only be caused by contaminating distributions that are concentrated near low‐dimensional subspaces.     

On a class of nonparametric tests for the treatment vs control problem

The classical problem of testing treatment versus control is revisited by considering a class of test statistics based on a kernel that depends on a constant ‘a’. The proposed class includes the celebrated Wilcoxon-Mann-Whitnet statistics as a special case when ‘a’=1. It is shown that, with optimal choice of ‘a’ depending on the underlying distribution, the optimal member performs better (in terms of Pitman efficiency) than the Wilcoxon-Mann-Whitney and the Median tests for a wide range of underlying distributions. An extended Hodges-Lehmann type point estimator of the shift prameter corresponding to the proposed ‘optimal’ test statistic is also derived.     

Deterministic Evolution of Strength in Multiple Comparisons Models: Who is the Greatest Golfer?

We present a statistical methodology for fitting time‐varying rankings, by estimating the strength parameters of the Plackett–Luce multiple comparisons model at regularly spaced times for each ranked item. We use the little‐known method of barycentric rational interpolation to interpolate between the strength parameters so that a competitor's strength can be evaluated at any time. We chose the time‐varying strengths to evolve deterministically rather than stochastically, a preference that we reason often has merit. There are many statistical and computational problems to overcome on fitting anything beyond ‘toy’ data sets. The methodological innovations here include a method for maximizing a likelihood function for many parameters, approximations for modelling tied data and an approach to the elimination of secular drift of the estimated ‘strengths’. The methodology has obvious applications to fields such as marketing, although we demonstrate our approach by analysing a large data set of golf tournament results, in search of an answer to the question ‘who is the greatest golfer of all time?’     

Joint models for mixed categorical outcomes: a study of HIV risk perception and disease status in Mozambique

Two types of bivariate models for categorical response variables are introduced to deal with special categories such as ‘unsure’ or ‘unknown’ in combination with other ordinal categories, while taking additional hierarchical data structures into account. The latter is achieved by the use of different covariance structures for a trivariate random effect. The models are applied to data from the INSIDA survey, where interest goes to the effect of covariates on the association between HIV risk perception (quadrinomial with an ‘unknown risk’ category) and HIV infection status (binary). The final model combines continuation-ratio with cumulative link logits for the risk perception, together with partly correlated and partly shared trivariate random effects for the household level. The results indicate that only age has a significant effect on the association between HIV risk perception and infection status. The proposed models may be useful in various fields of application such as social and biomedical sciences, epidemiology and public health.     

A re-appraisal of tests for normality

Tests for normality can be divided into two groups - those based upon a function of the empirical distribution function and those based upon a function of the original observations. The latter group of statistics test spherical symmetry and not necessarily normality. If the distribution is completely specified then the first group can be used to test for ‘spherical’ normality. However, if the distribution is incompletely specified and F‘‘x_i - x’/s’ is used these test statistics also test sphericity rather than normality. A Monte Carlo study was conducted for the completely specified case, to investigate the sensitivity of the distance tests to departures from normality when the alternative distributions are non-normal spherically symmetric laws. A “new” test statistic is proposed for testing a completely specified normal distribution     

The behavior of estimated measures of association in small and moderate sample sizes for 2×3 tables

We consider the distributions of Goodman and Kruskal's G, Kendall's tau-b, and correlation coefficients rho and rho-s for sample sizes 10‘10’40 from 2×3 tables. The results are compared with asymptotic theory. It is found that the convergence of G to its asymptotic normal distribution is much slower than the convergence of the other measures to theirs, and that G is more likely to be significantly biased. However, the variances and biases of all four measures come close to their asymptotic values for quite moderate sample sizes.     

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.     

Robust test for means when population variances are unequal

We consider the problem of testing the equality of two population means when the population variances are not necessarily equal. We propose a Welch-type statistic, say T^* _c, based on Tiku!s ‘1967, 1980’ modified maximum likelihood estimators, and show that this statistic is robust to symmetric and moderately skew distributions. We investigate the power properties of the statistic T^* _c; T^* _c clearly seems to be more powerful than Yuen's ‘1974’ Welch-type robust statistic based on the trimmed sample means and the matching sample variances. We show that the analogous statistics based on the ‘adaptive’ robust estimators give misleading Type I errors. We generalize the results to testing linear contrasts among k population means