Negative dependence and the Simes inequality

It is shown that the Simes inequality is reversed for a broad class of negatively dependent distributions.     

Two sets of isotones for comparing tests of exponentiality

Isotones   are a deterministic graphical device introduced by Mudholkar et al. [1991. A graphical procedure for comparing goodness-of-fit tests. J. Roy. Statist. Soc. B 53, 221–232], in the context of comparing some tests of normality. An isotone of a test is a contour of p$p$ values of the test applied to “ideal samples”, called profiles, from a two-shape-parameter family representing the null and the alternative distributions of the parameter space. The isotone is an adaptation of Tukey's sensitivity curves, a generalization of Prescott's stylized sensitivity contours, and an alternative to the isodynes   of Stephens. The purpose of this paper is two fold. One is to show that the isotones can provide useful qualitative information regarding the behavior of the tests of distributional assumptions other than normality. The other is to show that the qualitative conclusions remain the same from one two-parameter family of alternatives to another. Towards this end we construct and interpret the isotones of some tests of the composite hypothesis of exponentiality, using the profiles of two Weibull extensions, the generalized Weibull and the exponentiated Weibull families, which allow IFR, DFR, as well as unimodal and bathtub failure rate alternatives. Thus, as a by-product of the study, it is seen that a test due to Csörg? et al. [1975. Application of characterizations in the area of goodness-of-fit. In: Patil, G.P., Kotz, S., Ord, J.K. (Eds.), Statistical Distributions in Scientific Work, vol. 2. Reidel, Boston, pp. 79–90], and Gnedenko's Q(r)$Q(r)$ test [1969. Mathematical Methods of Reliability Theory. Academic Press, New York], are appropriate for detecting monotone failure rate alternatives, whereas a bivariate F$F$ test due to Lin and Mudholkar [1980. A test of exponentiality based on the bivariate F$F$ distribution. Technometrics 22, 79–82] and their entropy test [1984. On two applications of characterization theorems to goodness-of-fit. Colloq. Math. Soc. Janos Bolyai 45, 395–414] can detect all alternatives, but are especially suitable for nonmonotone failure rate alternatives.     

An optimality criterion for supersaturated designs with quantitative factors

A supersaturated design (SSD) is a factorial design in which the degrees of freedom for all its main effects exceed the total number of distinct factorial level-combinations (runs) of the design. Designs with quantitative factors, in which level permutation within one or more factors could result in different geometrical structures, are very different from designs with nominal ones which have been treated as traditional designs. In this paper, a new criterion is proposed for SSDs with quantitative factors. Comparison and analysis for this new criterion are made. It is shown that the proposed criterion has a high efficiency in discriminating geometrically nonisomorphic designs and an advantage in computation.     

D-optimal designs for logistic models with three and four parameters

Logistic functions are used in different applications, including biological growth studies and assay data analysis. Locally D-optimal designs for logistic models with three and four parameters are investigated. It is shown that these designs are minimally supported. Efficiencies are computed for equally spaced and uniform designs.     

Prediction intervals for general balanced linear random models

The main interest of prediction intervals lies in the results of a future sample from a previously sampled population. In this article, we develop procedures for the prediction intervals which contain all of a fixed number of future observations for general balanced linear random models. Two methods based on the concept of a generalized pivotal quantity (GPQ) and one based on ANOVA estimators are presented. A simulation study using the balanced one-way random model is conducted to evaluate the proposed methods. It is shown that one of the two GPQ-based and the ANOVA-based methods are computationally more efficient and they also successfully maintain the simulated coverage probabilities close to the nominal confidence level. Hence, they are recommended for practical use. In addition, one example is given to illustrate the applicability of the recommended methods.     

Cost-cautious designs for confirmatory bioassay

Confirmatory bioassay experiments take place in late stages of the drug discovery process when a small number of compounds have to be compared with respect to their properties. As the cost of the observations may differ considerably, the design problem is well specified by the cost of compound used rather than by the number of observations. We show that cost-efficient designs can be constructed using useful properties of the minimum support designs. These designs are particularly suited for studies where the parameters of the model to be estimated are known with high accuracy prior to the experiment, although they prove to be robust against typical inaccuracies of these values. When the parameters of the model can only be specified with ranges of values or by a probability distribution, we use a Bayesian criterion of optimality to construct the required designs. Typically, the number of their support points depends on the prior knowledge for the model parameters. In all cases we recommend identifying a set of designs with good statistical properties but different potential costs to choose from.     

Predictability of designs which adjust for imbalances in prognostic factors

Minimisation is a method often used in clinical trials to balance the treatment groups with respect to some prognostic factors. In the case of two treatments, the predictability of this method is calculated for different numbers of factors, different numbers of levels of each factor and for different proportions of the population at each level. It is shown that if we know nothing about the previous patients except the last treatment allocation, the next treatment can be correctly guessed more than 60% of the time if no biased coin is used. If the two previous assignments are known to have been the same, the next treatment can be guessed correctly around 80% of the time. Therefore, it is suggested that a biased coin should always be used with minimisation. Different choices of biased coin are investigated in terms of the reduction in predictability and the increase in imbalance that they produce. An alternative design to minimisation which makes use of optimum design theory is also investigated, by means of simulation, and does not appear to have any clear advantages over minimisation with a biased coin.     

Defining equations for two-level factorial designs

Defining equations are introduced in the context of two-level factorial designs and they are shown to provide a concise specification of both regular and nonregular designs. The equations are used to find orthogonal arrays of high strength and some optimal designs. The latter optimal designs are formed in a new way by augmenting notional orthogonal arrays which are allowed to have some runs with a negative number of replicates before augmentation. Defining equations are also shown to be useful when the factorial design is blocked.     

Bayesian variable selection using an adaptive powered correlation prior

The problem of selecting the correct subset of predictors within a linear model has received much attention in recent literature. Within the Bayesian framework, a popular choice of prior has been Zellner's g$g$-prior which is based on the inverse of empirical covariance matrix of the predictors. An extension of the Zellner's prior is proposed in this article which allow for a power parameter on the empirical covariance of the predictors. The power parameter helps control the degree to which correlated predictors are smoothed towards or away from one another. In addition, the empirical covariance of the predictors is used to obtain suitable priors over model space. In this manner, the power parameter also helps to determine whether models containing highly collinear predictors are preferred or avoided. The proposed power parameter can be chosen via an empirical Bayes method which leads to a data adaptive choice of prior. Simulation studies and a real data example are presented to show how the power parameter is well determined from the degree of cross-correlation within predictors. The proposed modification compares favorably to the standard use of Zellner's prior and an intrinsic prior in these examples.     

Kernel density estimators for random fields satisfying an interlaced mixing condition

For a sequence of strictly stationary random fields that are uniformly ρ′-mixing$\rho \prime \mathrm{-}\mathrm{mixing}$ and satisfy a Lindeberg condition, a central limit theorem is obtained for sequences of “rectangular” sums from the given random fields. The “Lindeberg CLT” is then used to prove a CLT for some kernel estimators of probability density for some strictly stationary random fields satisfying ρ′-mixing$\rho \prime \mathrm{-}\mathrm{mixing}$, and whose probability density and joint densities are absolutely continuous.     

Statistical analysis for rounded data

When random variables do not take discrete values, observed data are often the rounded values of continuous random variables. Errors caused by rounding of data are often neglected by classical statistical theories. While some pioneers have identified and made suggestions to rectify the problem, few suitable approaches were proposed. In this paper, we propose an approximate MLE (AMLE) procedure to estimate the parameters and discuss the consistency and asymptotic normality of the estimates. For our illustration, we shall consider the estimates of the parameters in AR(p)$\mathit{AR}(p)$ and MA(q)$\mathit{MA}(q)$ models for rounded data.     

Using pilot study information to increase efficiency in clinical trials

It is often necessary to conduct a pilot study to determine the sample size required for a clinical trial. Due to differences in sampling environments, the pilot data are usually discarded after sample size calculation. This paper tries to use the pilot information to modify the subsequent testing procedure when a two-sided t$t$-test or a regression model is used to compare two treatments. The new test maintains the required significance level regardless of the dissimilarity between the pilot and the target populations, but increases the power when the two are similar. The test is constructed based on the posterior distribution of the parameters given the pilot study information, but its properties are investigated from a frequentist's viewpoint. Due to the small likelihood of an irrelevant pilot population, the new approach is a viable alternative to the current practice.     

Simultaneous inference on diversity of biological communities

In environmental statistics, surveys on the structure of biological communities are generally carried out by focusing on diversity indexes. A more complete analysis may be performed by means of an appropriate function giving a spectrum of different measures of diversity: diversity profiles. They can be expressed as a function of the unknown abundance vector of the ecological population. In this paper we develop a non parametric approach based on bootstrap in order to make inference on diversity profiles. The proposed procedure is applied on biological data of four parks in Milan, Italy.     

Asymptotic normality of kernel type regression estimators for random fields

The asymptotic normality of the Nadaraya–Watson regression estimator is studied for α-mixing$\alpha \text{-}\mathrm{mixing}$ random fields. The infill-increasing setting is considered, that is when the locations of observations become dense in an increasing sequence of domains. This setting fills the gap between continuous and discrete models. In the infill-increasing case the asymptotic normality of the Nadaraya–Watson estimator holds, but with an unusual asymptotic covariance structure. It turns out that this covariance structure is a combination of the covariance structures that we observe in the discrete and in the continuous case.     

Optimal block designs for three treatments when observations are correlated

E$E$-optimal designs for comparing three treatments in blocks of size three are identified, where intrablock observations are correlated according to a first order autoregressive error process with parameter ρ∈(0,1)$\rho \in (0,1)$. For number of blocks b   of the form b=3n+1$b=3n+1$, there are two distinct optimal designs depending on the value of ρ$\rho$, with the best design being unequally replicated for large ρ$\rho$. For other values of b$b$, binary, equireplicate designs with specified within-block assignment patterns are best. In many cases, the stronger majorization optimality is established.     

Least squares estimation of the coefficients of a partially explosive model,with polynomial regressions of same degree,generating a pair of related time series

This paper is concerned with a partially explosive linear model with polynomial regression components generating a pair of related time series. The least squares estimates of the coefficients are shown to be √N-consistent and asymptotically singular normal, when the degrees of polynomial regression components are same, thus generalising a result due to Venkataraman (1974).     

A surprising property of uniformly best linear affine estimation in linear regression when prior information is fuzzy

It is already shown in Arnold and Stahlecker (2009) that, in linear regression, a uniformly best estimator exists in the class of all Γ-compatible$\Gamma \text{-compatible}$ linear affine estimators. Here, prior information is given by a fuzzy set Γ$\Gamma$ defined by its ellipsoidal α-cuts$\alpha \text{-cuts}$. Surprisingly, such a uniformly best linear affine estimator is uniformly best not only in the class of all Γ-compatible$\Gamma \text{-compatible}$ linear affine estimators but also in the class of all estimators satisfying a very weak and sensible condition. This property of a uniformly best linear affine estimator is shown in the present paper. Furthermore, two illustrative special cases are mentioned, where a generalized least squares estimator on the one hand and a general ridge or Kuks–Olman estimator on the other hand turn out to be uniformly best.