Bayesian robustness in bidimensional models: Prior independence

When θ is a multidimensional parameter, the issue of prior dependence or independence of coordinates is a serious concern. This is especially true in robust Bayesian analysis; Lavine et al. (J. Amer. Statist. Assoc.86, 964–971 (1991)) show that allowing a wide range of prior dependencies among coordinates can result in near vacuous conclusions. It is sometimes possible, however, to make confidently the judgement that the coordinates of θ are independent a priori and, when this can be done, robust Bayesian conclusions improve dramatically. In this paper, it is shown how to incorporate the independence assumption into robust Bayesian analysis involving -contamination and density band classes of priors. Attention is restricted to the case θ = (θ₁, θ₂) for clarity, although the ideas generalize.     

Analytical proofs of classical inequalities between Spearman's and Kendall's

Short analytical proofs are given for classical inequalities due to Daniels [1950. Rank correlation and population models. J. Roy. Statist. Soc. Ser. B 12, 171–181; 1951. Note on Durbin and Stuart's formula for E(r_s). J. Roy. Statist. Soc. Ser. B 13, 310] and Durbin and Stuart [1951. Inversions and rank correlation coefficients. J. Roy. Statist. Soc. Ser. B 13, 303–309] relating Spearman's ρ and Kendall's τ.     

Order statistics from trivariate normal and -distributions in terms of generalized skew-normal and skew- distributions

We consider here a generalization of the skew-normal distribution, GSN(λ₁,λ₂,ρ), defined through a standard bivariate normal distribution with correlation ρ, which is a special case of the unified multivariate skew-normal distribution studied recently by Arellano-Valle and Azzalini [2006. On the unification of families of skew-normal distributions. Scand. J. Statist. 33, 561–574]. We then present some simple and useful properties of this distribution and also derive its moment generating function in an explicit form. Next, we show that distributions of order statistics from the trivariate normal distribution are mixtures of these generalized skew-normal distributions; thence, using the established properties of the generalized skew-normal distribution, we derive the moment generating functions of order statistics, and also present expressions for means and variances of these order statistics.Next, we introduce a generalized skew-t_ν distribution, which is a special case of the unified multivariate skew-elliptical distribution presented by Arellano-Valle and Azzalini [2006. On the unification of families of skew-normal distributions. Scand. J. Statist. 33, 561–574] and is in fact a three-parameter generalization of Azzalini and Capitanio's [2003. Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t distribution. J. Roy. Statist. Soc. Ser. B 65, 367–389] univariate skew-t_ν form. We then use the relationship between the generalized skew-normal and skew-t_ν distributions to discuss some properties of generalized skew-t_ν as well as distributions of order statistics from bivariate and trivariate t_ν distributions. We show that these distributions of order statistics are indeed mixtures of generalized skew-t_ν distributions, and then use this property to derive explicit expressions for means and variances of these order statistics.     

Design of robust parameter experiments in a continuous space using an evolutionary optimization algorithm

Two common experimental designs used in robust parameter design (RPD) are crossed array and mixed resolution designs. However, the prohibited number of runs, constraints in the design space or special model requirements render some of these designs inadequate. This paper presents the application of an evolutionary strategy to produce nearly optimal design matrices for RPD. The designs are derived by solving a nonlinear optimization problem involving both 𝒟- and 𝒢-efficiency simultaneously. The methodology presented allows the user to obtain new exact designs for a specific number of runs, and a particular experimental region. The combination of 𝒟- and 𝒢-efficiency results in experimental designs that outperform the corresponding benchmarks.     

General upper bound for the number of blocks having a given number of treatments common with a given block

The purpose of this paper is systematically to derive the general upper bound for the number of blocks having a given number of treatments common with a given block of certain incomplete block designs. The approach adopted here is based on the spectral decomposition of N′N for the incidence matrix N of a design, where N' is the transpose of the matrix N. This approach will lead us to upper bounds for incomplete block designs, in particular for a large number of partially balanced incomplete block (PBIB) designs, which are not covered with the standard approach (Shah 1964, 1966), Kapadia (1966)) of using well known relations between blocks of the designs and their association schemes. Several results concerning block structure of block designs are also derived from the main theorem. Finally, further generalizations of the main theorem are discussed with some illustrations.     

Klingenberg structures and partial designs I: Congruence relations and solutions

We consider generalizations of projective Klingenberg and projective Hjelmslev planes, mainly (b, c)-K-structures. These are triples (φ, Π, Π′) where Π and Π′ are incidence structures and φ : Π → Π′ is an epimorphism which satisfies certain lifting axioms for double flags. The congruence relations of such triples are investigated, leading to nice factorizations of φ. Two integer invariants are associated with each congruence relation, generalizing a theorem of Kleinfeld on projective Hjelmslev planes. These invariants are completely characterized for a special class of triples: the balanced, minimally uniform neighbor cohesive (1,1)-K-structures. We show that a balanced neighbor cohesive (1,1)-K-structure Π “over” a PBIBD Π′ is again a PBIBD and compute its invariants. Several methods are given for constructing symmetric “divisible” PBIBD's on arbitrarily many classes.     

Some new classes of multivariate survival distribution functions

The concept of the univariate mean remaining life (m.r.l.) function is generalized to the multivariate case. The multivariate mean remaining life (m.m.r.l.) function is utilized to introduce four new classes of multivariate survival distribution functions (s.d.f.'s). Each of these classes is a new generalization of the univariate decreasing mean remaining life (DMRL) class of s.d.f.'s. The duals of these classes are introduced. Some properties, physical interpretation, and relationships among these classes are investigated. Also for each case, the class of s.d.f.'s common in a class and its dual is characterized.     

DDMA-charts: Nonparametric multivariate moving average control charts based on data depth

Summary: This paper studies the DDMA–chart, a data depth based moving–average control chart for monitoring multivariate data. This chart is nonparametric and it can detect simultaneously location and scale changes in the process. It improves upon the existing r– and Q–chart in the efficiency of detecting location changes. Both theoretical justifications and simulation studies are provided. Comparisons with some existing multivariate control charts via simulation results are also provided. Some applications of the DDMA–chart to the analysis of airline performance data (collected by the FAA) are demonstrated. The results indicate that the DDMA–chart is an effective nonparametric multivariate control chart.*Research supported in part by grants from the National Science Foundation, the National Security Agency, and the Federal Aviation Administration. The discussion on aviation safety in this paper reects the views of the authors, who are solely responsible for the accuracy of the analysis results presented herein, and does not necessarily reect the official view or policy of the FAA. The dataset used in this paper has been partially masked in order to protect confidentiality.     

Bayesian quadrature with non-normal approximating functions

We consider an efficient Bayesian approach to estimating integration-based posterior summaries from a separate Bayesian application. In Bayesian quadrature we model an intractable posterior density function f(·) as a Gaussian process, using an approximating function g(·), and find a posterior distribution for the integral of f(·), conditional on a few evaluations of f (·) at selected design points. Bayesian quadrature using normal g (·) is called Bayes-Hermite quadrature. We extend this theory by allowing g(·) to be chosen from two wider classes of functions. One is a family of skew densities and the other is the family of finite mixtures of normal densities. For the family of skew densities we describe an iterative updating procedure to select the most suitable approximation and apply the method to two simulated posterior density functions.     

Estimating parameters of a selected Pareto population

Let Π₁,…,Π_k be k populations with Π_i being Pareto with unknown scale parameter α_i and known shape parameter β_i;i=1,…,k. Suppose independent random samples (X_i1,…,X_in), i=1,…,k of equal size are drawn from each of k populations and let X_i denote the smallest observation of the ith sample. The population corresponding to the largest X_i is selected. We consider the problem of estimating the scale parameter of the selected population and obtain the uniformly minimum variance unbiased estimator (UMVUE) when the shape parameters are assumed to be equal. An admissible class of linear estimators is derived. Further, a general inadmissibility result for the scale equivariant estimators is proved.     

Nonparametric tests for the mean of a non-negative population

We construct level-α tests for testing the null hypothesis that the mean of a non-negative population falls below a prespecified nominal value. These tests make no assumption about the distribution function other than that it be supported on [0,∞). Simple tests are derived based on either the sample mean or the sample product. The nonparametric likelihood ratio test is also discussed in this context. We also derive the uniformly most powerful monotone (UMP) tests for special cases.     

Robust multivariate Bayesian forecasting under functional distortions in the -metric

This paper investigates robustness of multivariate forecasting in the Bayesian framework. The minimax approach is used to construct robust statistical procedures under deviations from hypothetical assumptions. The deviations are defined as functional distortions using the χ²-pseudo-metric. Two cases of deviations are considered: distortions of parameter distribution and distortions of joint distribution of observations and parameters. Explicit forms for the guaranteed upper risk functional are obtained and integral equations for robust prediction statistics are given for both cases.     

Variance estimation based on blocked 3×2 cross-validation in high-dimensional linear regression

In high-dimensional linear regression, the dimension of variables is always greater than the sample size. In this situation, the traditional variance estimation technique based on ordinary least squares constantly exhibits a high bias even under sparsity assumption. One of the major reasons is the high spurious correlation between unobserved realized noise and several predictors. To alleviate this problem, a refitted cross-validation (RCV) method has been proposed in the literature. However, for a complicated model, the RCV exhibits a lower probability that the selected model includes the true model in case of finite samples. This phenomenon may easily result in a large bias of variance estimation. Thus, a model selection method based on the ranks of the frequency of occurrences in six votes from a blocked 3×2 cross-validation is proposed in this study. The proposed method has a considerably larger probability of including the true model in practice than the RCV method. The variance estimation obtained using the model selected by the proposed method also shows a lower bias and a smaller variance. Furthermore, theoretical analysis proves the asymptotic normality property of the proposed variance estimation.     

Model weights for model choice and averaging

A method is suggested to estimate posterior model probabilities and model averaged parameters via MCMC sampling under a Bayesian approach. The estimates use pooled output for J models (J>1) whereby all models are updated at each iteration. Posterior probabilities are based on averages of continuous weights obtained for each model at each iteration, while samples of averaged parameters are obtained from iteration specific averages that are based on these weights. Parallel sampling of models assists in deriving posterior densities for parameter contrasts between models and in assessing hypotheses regarding model averaged parameters. Four worked examples illustrate application of the approach, two involving fixed effect regression, and two involving random effects.     

Generalized Supremum Tests for the Equality of Cause Specific Hazard Rates

In this paper we propose two new classes of asymptotically distribution-free Renyi-type tests for testing the equality of two risks in a competing risk model with possible censoring. This work extends the work of Aly, Kochar and McKeague [1994, Journal of American Statistical Association, 89, 994–999] and many of the existing tests for this problem belong to these newly proposed classes. The asymptotic properties of the proposed tests are investigated. Simulation studies are done to compare the performance with existing tests. A competing risks data set is analyzed to demonstrate the usefulness of the procedure.     

On the ‘optimal’ density power divergence tuning parameter

The density power divergence, indexed by a single tuning parameter α, has proved to be a very useful tool in minimum distance inference. The family of density power divergences provides a generalized estimation scheme which includes likelihood-based procedures (represented by choice α=0 for the tuning parameter) as a special case. However, under data contamination, this scheme provides several more stable choices for model fitting and analysis (provided by positive values for the tuning parameter α). As larger values of α necessarily lead to a drop in model efficiency, determining the optimal value of α to provide the best compromise between model-efficiency and stability against data contamination in any real situation is a major challenge. In this paper, we provide a refinement of an existing technique with the aim of eliminating the dependence of the procedure on an initial pilot estimator. Numerical evidence is provided to demonstrate the very good performance of the method. Our technique has a general flavour, and we expect that similar tuning parameter selection algorithms will work well for other M-estimators, or any robust procedure that depends on the choice of a tuning parameter.     

A generalized predictive criterion for model selection

Given a random sample from some unknown model belonging to a finite class of parametric models, assume that the estimate of the density of a future observation is of interest San Martini & Spezzaferri (1984) proposed for this problem a predictive criterion based on the logarithmic utility function. The present authors investigate a generalization of this criterion that uses as a loss function an element of the class of α‐divergences discussed by Ali & Silvey (1966) and Csiszár (1967). They also discuss briefly the case in which the class of models considered is not exhaustive.     

Saddlepoint approximations for some models of circular data

In this article we provide saddlepoint approximations for some important models of circular data. The particularity of these saddlepoint approximations is that they do not require solving the saddlepoint equation iteratively, so their evaluation is immediate. We first give very accurate approximations to P-values, critical values and power functions for some optimal tests regarding the concentration parameter under wrapped symmetric α-stable and circular normal models. Then, we consider an approximation to the distribution of a projection of the two-dimensional Pearson random walk with exponential step sizes.     

Erste Ergebnisse faktischer Anonymisierung wirtschaftsstatistischer Einzeldaten

Summary: This paper presents first results of the project Factual anonymization of business microdata. The project aims at creating and providing scientific–use files of business data. For this, appropriate anonymization strategies and methods are created. The anonymization procedures are judged by whether they can guarantee the factual anonymity of company and operating data without limiting their potential value for statistical analysis. As an example the survey of German cost structure in the processing industry is considered. Here, both the impact of different anonymization methods on the potential use of the data as well as their effect on the re–identification risk is examined.

---

Zusammenfassung: Mit diesem Beitrag werden erste Ergebnisse des Projektes Faktische Anonymisierung wirtschaftsstatistischer Einzeldaten vorgestellt. Ziel des Projektes ist es, die Grundlagen für die Erstellung von Scientific–use–files im Bereich der Unternehmens– und Betriebsdaten zu legen. Hierzu werden entsprechende Anonymisierungsverfahren und –strategien entwickelt. Die Anonymisierungsverfahren müssen danach beurteilt werden, ob sie die faktische Anonymität von Unternehmens– und Betriebsdaten sicherstellen können, ohne das Analysepotenzial über Gebühr einzuschränken. Am Beispiel der Kostenstrukturerhebung im Verarbeitenden Gewerbe werden erste Ergebnisse vorgestellt. Dabei wird sowohl untersucht, wie sich verschiedene Anonymisierungsmaßnahmen auf das Analysepotenzial der Daten auswirken, als auch welche Effekte sie auf das Re–Identifikationsrisiko haben.

     

Testing concerning the homogeneity of some predator–prey populations

In this paper, we study an inference problem for a stochastic model where k deterministic Lotka–Volterra systems of ordinary differential equations (ODEs) are perturbed with k pairs of random errors. The k deterministic systems describe the ecological interaction between k predator–prey populations. These k deterministic systems depend on unknown parameters. We consider the testing problem concerning the homogeneity between k pairs of the interaction parameters of the ODEs. We assume that the k pairs of random errors are independent and that, each pair follows correlated Ornstein–Uhlenbeck processes. Thus, we extend the stochastic model suggested in Froda and Colavita [2005. Estimating predator–prey systems via ordinary differential equations with closed orbits. Aust. N.Z. J. Stat. 2, 235–254] as well as in Froda and Nkurunziza [2007. Prediction of predator–prey populations modeled by perturbed ODE. J. Math. Biol. 54, 407–451] where k=1. Under this statistical model, we propose a likelihood ratio test and study the asymptotic properties of this test. Finally, we highlight the performance of our method through some simulations studies.