Characterizations of exponentiality within the HNBUE class and related tests

A harmonic new better than used in expectation (HNBUE) variable is a random variable which is dominated by an exponential distribution in the convex stochastic order. We use a recently obtained condition on stochastic equality under convex domination to derive characterizations of the exponential distribution and bounds for HNBUE variables based on the mean values of the order statistics of the variable. We apply the results to generate discrepancy measures to test if a random variable is exponential against the alternative that is HNBUE, but not exponential.     

Consistency and power in tests with shape-restricted alternatives

Likelihood ratio tests of constant vs. monotone regression function, as well as linear vs. convex regression function and other tests with shape-restricted alternatives, are known to have null distributions equivalent to mixtures of beta random variates. The monotone and convex regression estimators are known to be inconsistent at the endpoints, where there is “spiking.” This spiking affects the critical values of the test statistic. Modified versions of the monotone and convex regression estimators are proposed that are consistent everywhere; when the modified versions are used in hypothesis testing, the null distribution is again an exact mixture of beta distributions, with different mixing parameters. Simulations show that the power of the test using the modified version is larger for the examples chosen.     

A new nonparametric estimator for the mean of the selected population

A modified bootstrap estimator of the mean of the population selected from two populations is proposed which is a convex combination of the two sample means, where the weights are random quantities. The estimator is shown to be strongly consistent. The small sample behavior of the estimator is investigated and compared with some competitors by means of Monte Carlo studies. It is found that the newly proposed estimator has smaller mean squared error for a wide range of parameter values.     

Bayesian revision of a prior given prior-data conflict,expert opinion,or a similar insight: a large-deviation approach

Learning from model diagnostics that a prior distribution must be replaced by one that conflicts less with the data raises the question of which prior should instead be used for inference and decision. The same problem arises when a decision maker learns that one or more reliable experts express unexpected beliefs. In both cases, coherence of the solution would be guaranteed by applying Bayes's theorem to a distribution of prior distributions that effectively assigns the initial prior distribution a probability arbitrarily close to 1. The new distribution for inference would then be the distribution of priors conditional on the insight that the prior distribution lies in a closed convex set that does not contain the initial prior. A readily available distribution of priors needed for such conditioning is the law of the empirical distribution of sufficiently large number of independent parameter values drawn from the initial prior. According to the Gibbs conditioning principle from the theory of large deviations, the resulting new prior distribution minimizes the entropy relative to the initial prior. While minimizing relative entropy accommodates the necessity of going beyond the initial prior without departing from it any more than the insight demands, the large-deviation derivation also ensures the advantages of Bayesian coherence. This approach is generalized to uncertain insights by allowing the closed convex set of priors to be random.     

Projection Mean–Variance Bounds on Expectations of kth Record Values from Restricted Families

We present sharp mean–variance bounds for expectations of kth record values based on distributions coming from restricted families of distributions. These families are defined in terms of convex or star ordering with respect to generalized Pareto distribution. The bounds for expectations of kth record values from DD, DFR, DDA, and DFRA families are special cases of our results. The bounds are derived by application of the projection method.     

On functionals of random convex hulls

The completeness of the induced distribution of the convex hull of a random set of points drawn uniformly from a convex region seems not to have been noticed before. The result generalizes a well-known result for dimension one. As a consequence, there exists a theory of best unbiased estimation for certain functionals of the convex region. For example, a best estimator of the centroid of the convex region can be constructed. This estimator is distinct from the centroid of the convex hull. Therefore another generalization of a property holding in dimension one, namely the unbiased-ness of the centroid of the convex hull, is seen to fail.     

Testing non-oblique hypotheses

This paper considers a likelihood ratio test for testing hypotheses defined by non-oblique closed convex cones, satisfying the so called iteration projection property, in a set of k normal means. We obtain the critical values of the test using the Chi-Bar-Squared distribution. The obtuse cones are introduced as a particular class of cones which are non-oblique with every one of their faces. Examples with the simple tree order cone and the total order cone are given to illustrate the results.     

Testing uniformity for the case of a planar unknown support

A new test is proposed for the hypothesis of uniformity on bi‐dimensional supports. The procedure is an adaptation of the “distance to boundary test” (DB test) proposed in Berrendero, Cuevas, & Vázquez‐Grande (2006). This new version of the DB test, called DBU test, allows us (as a novel, interesting feature) to deal with the case where the support S of the underlying distribution is unknown. This means that S is not specified in the null hypothesis so that, in fact, we test the null hypothesis that the underlying distribution is uniform on some support S belonging to a given class ${\cal C}$ . We pay special attention to the case that ${\cal C}$ is either the class of compact convex supports or the (broader) class of compact λ‐convex supports (also called r‐convex or α‐convex in the literature). The basic idea is to apply the DB test in a sort of plug‐in version, where the support S is approximated by using methods of set estimation. The DBU method is analysed from both the theoretical and practical point of view, via some asymptotic results and a simulation study, respectively. The Canadian Journal of Statistics 40: 378–395; 2012 © 2012 Statistical Society of Canada     

Discrete regularized discriminant analysis

A method of regularized discriminant analysis for discrete data, denoted DRDA, is proposed. This method is related to the regularized discriminant analysis conceived by Friedman (1989) in a Gaussian framework for continuous data. Here, we are concerned with discrete data and consider the classification problem using the multionomial distribution. DRDA has been conceived in the small-sample, high-dimensional setting. This method has a median position between multinomial discrimination, the first-order independence model and kernel discrimination. DRDA is characterized by two parameters, the values of which are calculated by minimizing a sample-based estimate of future misclassification risk by cross-validation. The first parameter is acomplexity parameter which provides class-conditional probabilities as a convex combination of those derived from the full multinomial model and the first-order independence model. The second parameter is asmoothing parameter associated with the discrete kernel of Aitchison and Aitken (1976). The optimal complexity parameter is calculated first, then, holding this parameter fixed, the optimal smoothing parameter is determined. A modified approach, in which the smoothing parameter is chosen first, is discussed. The efficiency of the method is examined with other classical methods through application to data.     

Non-parametric maximum likelihood estimation of interval-censored failure time data subject to misclassification

The paper considers non-parametric maximum likelihood estimation of the failure time distribution for interval-censored data subject to misclassification. Such data can arise from two types of observation scheme; either where observations continue until the first positive test result or where tests continue regardless of the test results. In the former case, the misclassification probabilities must be known, whereas in the latter case, joint estimation of the event-time distribution and misclassification probabilities is possible. The regions for which the maximum likelihood estimate can only have support are derived. Algorithms for computing the maximum likelihood estimate are investigated and it is shown that algorithms appropriate for computing non-parametric mixing distributions perform better than an iterative convex minorant algorithm in terms of time to absolute convergence. A profile likelihood approach is proposed for joint estimation. The methods are illustrated on a data set relating to the onset of cardiac allograft vasculopathy in post-heart-transplantation patients.     

Detection of a pair of outliers in a sample from a Gumbel distribution with known scale parameter

Test procedures on outlier detection problems for Gumbel distribution are rarely available. Hence, a test statistic is proposed here for detection of a pair of upper and lower outliers from a Gumbel distribution with known scale parameter. The critical values of the statistic are obtained and some examples are also given to highlight the use of the statistic. The advantage of the proposed statistic is that the scale parameter, though assumed to be known is not explicitly involved in the determination of the critical values.     

Signatures in Markov-Modulated Processes

□ In recent years, signatures are widely used for analysis of coherent systems consisting of unreliable components. If component lifetimes are independent and identically distributed, then system lifetime distribution function is a convex combination of distribution functions of order statistics for component lifetimes. Coefficients of this convex combination are called signatures. This article considers the case when a system operates in a so-called random environment, i.e., component failure rates are jointly modulated by a finite-state continuous-time Markov chain. In this model, component lifetimes remain exchangeable. An expression for distribution function of time to system failure is derived. Here, a crucial role is played by an elaborated procedure of deriving a distribution function of order statistics for system component lifetimes. A numerical example illustrates the suggested approach and analyzes the influence of random environment on the distribution function of system lifetime.     

Linearity of convex mean residual life

It is shown that linearity of convex mean residual life identifies the distribution up to a scale in the class of distributions with zero mean. On this basis new characterizations of the uniform and, unexpectedly, the Student distribution with two degrees of freedom are obtained. It is observed that, strictly speaking, the conjecture of Nagaraja and Nevzerov about uniqueness property of the convex mean residual life time is incorrect. The condition under study is in obvious relations with linearity of regression of observations with respect to order statistics.     

An economical acceptance–rejection algorithm for uniform random variate generation over constrained simplexes

This paper develops an algorithm for uniform random generation over a constrained simplex, which is the intersection of a standard simplex and a given set. Uniform sampling from constrained simplexes has numerous applications in different fields, such as portfolio optimization, stochastic multi-criteria decision analysis, experimental design with mixtures and decision problems involving discrete joint distributions with imprecise probabilities. The proposed algorithm is developed by combining the acceptance–rejection and conditional methods along with the use of optimization tools. The acceptance rate of the algorithm is analytically compared to that of a crude acceptance–rejection algorithm, which generates points over the simplex and then rejects any points falling outside the intersecting set. Finally, using convex optimization, the setup phase of the algorithm is detailed for the special cases where the intersecting set is a general convex set, a convex set defined by a finite number of convex constraints or a polyhedron.     

Ordered Regressions

There are often situations where two or more regression functions are ordered over a range of covariate values. In this paper, we develop efficient constrained estimation and testing procedures for such models. Specifically, necessary and sufficient conditions for ordering generalized linear regressions are given and shown to unify previous results obtained for simple linear regression, for polynomial regression and in the analysis of covariance models. We show that estimating the parameters of ordered linear regressions requires either quadratic programming or semi‐infinite programming, depending on the shape of the covariate space. A distance‐type test for order is proposed. Simulations demonstrate that the proposed methodology improves the mean square error and power compared with the usual, unconstrained, estimation and testing procedures. Improvements are often substantial. The methodology is extended to order generalized linear models where convex semi‐infinite programming plays a role. The methodology is motivated by, and applied to, a hearing loss study.     

MVE, MCD, and MZE: A simulation study comparing convex body minimizers

Summary: In a recent work (Paris Scholz, 2002), a new robust estimator for convex bodies has been proposed, based on the estimation of a zonoid of a distribution. This so–called minimum volume zonoid estimator (MZE) is similar in type to the well–known robust approaches of the minimum volume ellipsoid (MVE) and the minimum covariance determinant (MCD), all three seeking for a subset of given data for which some criteria are minimized. We investigate the similarity between these three concepts by comparing which subsets are chosen to be the optimal ones.     

Depth and outliers for samples of sets and random sets distributions

We suggest several constructions suitable to define the depth of set-valued observations with respect to a sample of convex sets or with respect to the distribution of a random closed convex set. With the concept of a depth, it is possible to determine if a given convex set should be regarded an outlier with respect to a sample of convex closed sets. Some of our constructions are motivated by the known concepts of half-space depth and band depth for function-valued data. A novel construction derives the depth from a family of non-linear expectations of random sets. Furthermore, we address the role of positions of sets for evaluation of their depth. Two case studies concern interval regression for Greek wine data and detection of outliers in a sample of particles.     

A General Model of Random Variation

A statistical distribution of a random variable is uniquely represented by its normal-based quantile function. For a symmetrical distribution it is S-shaped (for negative kurtosis) and inverted S-shaped (otherwise). As skewness departs from zero, the quantile function gradually transforms into a monotone convex function (positive skewness) or concave function (otherwise). Recently, a new general modeling platform has been introduced, response modeling methodology, which delivers good representation to monotone convex relationships due to its unique “continuous monotone convexity” property. In this article, this property is exploited to model the normal-based quantile function, and explored using a set of 27 distributions.     

A flexible approach for multivariate mixed-effects models with non-ignorable missing values

We propose a flexible model approach for the distribution of random effects when both response variables and covariates have non-ignorable missing values in a longitudinal study. A Bayesian approach is developed with a choice of nonparametric prior for the distribution of random effects. We apply the proposed method to a real data example from a national long-term survey by Statistics Canada. We also design simulation studies to further check the performance of the proposed approach. The result of simulation studies indicates that the proposed approach outperforms the conventional approach with normality assumption when the heterogeneity in random effects distribution is salient.     

Convex transformation on survival functions and related dependence concepts

Some concepts of stochastic dependence for continuous bivariate distribution functions are investigated by defining a convex transformation on their reliability or survival functions. We also study notions of bivariate hazard rate and hazard dependence. Some dependence orderings are characterized by using convex transformation. To clarify the discussions, illustrative examples are given.