Search designs for 2m factorials derived from balanced arrays of strength 2(l+1) and ad-optimal search designs

We consider a search design for the 2^m type such that at most knonnegative effects can be searched among (l+1)-factor interactions and estimated along with the effects up to l- factor interactions, provided (l+1)-factor and higher order interactions are negligible except for the k effects. We investigate some properties of a search design which is yielded by a balanced 2^m design of resolution 2l+1 derived from a balanced array of strength 2(l+1). A necessary and sufficient condition for the balanced design of resolution 2l+1 to be a search design for k=1 is given. Optimal search designs for k=1 in the class of the balanced 2^m designs of resolution V (l=2), with respect to the AD-optimality criterion given by Srivastava (1977), with N assemblies are also presented, where the range of (m,N) is (m=6; 28≤N≤41), (m=7; 35≤N≤63) and (m=8; 44≤N≤74).     

Balanced fractional factorial designs of resolution 2l+1 for interesting effects orthogonal to some nuisance parameters: 2m1+m2 series

For 2^m₁+m₂ factorial designs, this paper investigates balanced fractional 2^m₁ factorial designs of resolution 2l+1 with some nuisance parameters concerning the second factors. They are derivable from partially balanced arrays and further permit estimation of the effects up to the l-factor interactions concerning the first factors orthogonally to the nuisance parameters.     

Balanced arrays of strength 4 and balanced fractional 3m factorial designs

A connection between a balanced fractional 2^m factorial design of resolution 2l + 1 and a balanced array of strength 2l with index set {μ₀, μ₁,…, μ_2l} was established by Yamamoto, Shirakura and Kuwada (1975). The main purpose of this paper is to give a connection between a balanced fractional 3^m factorial design of resolution V and a balanced array of strength 4, size N, m constraints, 3 levels and index set {λ_l0l1l2}.     

Existence Conditions for Balanced Fractional 2m Factorial Designs of Resolution 2l + 1 Derived from Simple Arrays

We consider a fractional 2^m factorial design derived from a simple array (SA) such that the (? + 1)-factor and higher-order interactions are negligible, where 2? ? m. The purpose of this article is to give a necessary and sufficient condition for an SA to be a balanced fractional 2^m factorial design of resolution 2? + 1. Such a design is concretely characterized by the suffixes of the indices of an SA.     

Some results on constructing two-level block designs with general minimum lower order confounding

Block designs are widely used in experimental situations where the experimental units are heterogeneous. The blocked general minimum lower order confounding (B-GMC) criterion is suitable for selecting optimal block designs when the experimenters have some prior information on the importance of ordering of the treatment factors. This paper constructs B-GMC 2^{n ? m}: 2^r designs with 5 × 2^l/16 + 1 ? n ? (N ? 2^l) < 2^{l ? 1} for l(r + 1 ? l ? n ? m ? 1), where 2^{n ? m}: 2^r denotes a two-level regular block design with N = 2^{n ? m} runs, n treatment factors, and 2^r blocks. With suitable choice of the blocking factors, each B-GMC block design has a common specific structure. Some examples illustrate the simple and effective construction method.     

Nonparametric comparison of curves with dependent errors

Suppose that data {(x _l,i,n , y _l,i,n ): l?=?1, …, k; i?=?1, …, n} are observed from the regression models: Y _l,i,n ?=?m _l (x _l,i,n )?+?? _l,i,n , l?=?1, …, k, where the regression functions {m _l } _l=1 ^k are unknown and the random errors {? _l,i,n } are dependent, following an MA(∞) structure. A new test is proposed for testing the hypothesis H ₀: m ₁?=?·?·?·?=?m _k , without assuming that {m _l } _l=1 ^k are in a parametric family. The criterion of the test derives from a Crámer-von-Mises-type functional based on different distances between {[mcirc]} _l and {[mcirc]} _s , l?≠?s, l, s?=?1, …, k, where {[mcirc] _l } _l=1 ^k are nonparametric Gasser–Müller estimators of {m _l } _l=1 ^k . A generalization of the test to the case of unequal design points, with different sample sizes {n _l } _l=1 ^k and different design densities {f _l } _l=1 ^k , is also considered. The asymptotic normality of the test statistic is obtained under general conditions. Finally, a simulation study and an analysis with real data show a good behavior of the proposed test.     

Random r-content of an r-simplex from beta-type-2 random points

We study the r-content Δ of the r -simplex generated by r+ 1 independent random points in R”. Each random point Z_j is isotropic and distributed according to λ||Z_j||² ～ beta-type-2(n/2, v), λ, v > 0. We provide an asymptotic normality result which is analogous to the conjecture made by Miles (1971). A method is introduced to work out the exact density of W = (rλ)^r(r!Δ)²/(r + |)^r+l and hence that of Δ. The distribution of W is also related to some hypothesis-testing problems in multivariate analysis. Furthermore, by using this method, the distribution of W or Δ can easily be simulated.     

Maximum Likelihood Estimation in the Bivariate Binomial (0,1) Distribution:Application TO 2×2 Tables

We consider a 2×2 contingency table, with dichotomized qualitative characters (A,A) and (B,B), as a sample of size n drawn from a bivariate binomial (0,1) distribution. Maximum likelihood estimates p?₁p?₂ and p? are derived for the parameters of the two marginals p₁p₂ and the coefficient of correlation. It is found that p? is identical to Pearson's (1904)?=(χ²/n)½, where ?² is Pearson's usual chi-square for the 2×2 table. The asymptotic variance-covariance matrix of p?_lp?₂and p is also derived.     

Invariance In the linear model – an argument for χ2 and F in nonnormal situations

For a suitable norm, conservation of the distance between expectation and hypothesis may furnish a basis for data reduction by invariance in the linear, not neces-sarily normal, model. If the norm is Euclidean (i.e. based on some inner product), the maximal invariant is a pair of sums of squares. This provides support for traditional χ² (or F) - methods also in nonnormal cases. If the norm is l^p p≠2, or the supnorm, the maximal invariant is, at the best a air of order statistics.     

A class of non-abelian 2-groups containing Menon difference sets

We show that every group in a certain class of 2-groups contains a Menon difference set. This provides further positive evidence for a conjecture of Dillon concerning 2-groups of order 2^2l which contain a normal subgroup isomorphic to Z^l₂. The conjecture, however, remains open.     

On the norm of alias matrices in balanced fractional 2m factorial designs of resolution 2l + 1

The norm $\text{6A6 = {tr(A′A)}}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ of the alias matrix A of a design can be used as a measure for selecting a design. In this paper, an explicit expression for 6A6 will be given for a balanced fractional 2^m factorial design of resolution 2l + 1 which obtained from a simple array with parameters (m; λ₀, λ₁,…, λ_m). This array is identical with a balanced array of strength m, m constraints and index set {λ₀, λ₁,…, λ_m}. In the class of the designs of resolution V (l = 2) obtained from S-arrays, ones which minimize 6A6 will be presented for any fixed N assemblies satisfying (i) m = 4, 11 ? N ? 16, (ii) m = 5, 16 ? N ? 32, and (iii) m = 6, 22 ? N ? 40.     

Monotonic minimax estimators of a 2×2 covariance matrix

Let S : 2 × 2 have a nonsingular Wishart distribution with unknown matrix σ and n degrees of freedom. For estimating σ two families of mimmax estimators, with respect to the entropy loss, are presented. These estimators are of the form σ(S) = Rø(L)R^t where R is orthogonal, L and Φ are diagonal, and RLR^T = S. Conditions under which the components of Φ and L follow the same order relation [i.e., writing Φ = diag(Φ₁,Φ₂) and L = diag(l₁,/₂) with l₁ ≥ l₂, we have Φ₁ ≥ Φ₂] are established. Comparisons with Stein's estimators and other orthogonally invariant estimators are discussed.     

Characterization of Balanced Fractional 3 m Factorial Designs of Resolution III

We consider a fractional 3 ^m factorial design derived from a simple array (SA), which is a balanced array of full strength, where the non negligible factorial effects are the general mean and the linear and quadratic components of the main effect, and m ≥ 2. In this article, we give a necessary and sufficient condition for an SA to be a balanced fractional 3 ^m factorial design of resolution III. Such a design is characterized by the suffixes of indices of an SA.     

Likelihood ratio inference for mean residual life

One of the basic parameters in survival analysis is the mean residual life M ₀. For right censored observation, the usual empirical likelihood based log-likelihood ratio leads to a scaled c₁²{\chi_1^2} limit distribution and estimating the scaled parameter leads to lower coverage of the corresponding confidence interval. To solve the problem, we present a log-likelihood ratio l(M ₀) by methods of Murphy and van der Vaart (Ann Stat 1471–1509, 1997). The limit distribution of l(M ₀) is the standard c₁²{\chi_1^2} distribution. Based on the limit distribution of l(M ₀), the corresponding confidence interval of M ₀ is constructed. Since the proof of the limit distribution does not offer a computational method for the maximization of the log-likelihood ratio, an EM algorithm is proposed. Simulation studies support the theoretical result.     

Bayesian Optimal Adaptive Estimation Using a Sieve Prior

We derive rates of contraction of posterior distributions on non‐parametric models resulting from sieve priors. The aim of the study was to provide general conditions to get posterior rates when the parameter space has a general structure, and rate adaptation when the parameter is, for example, a Sobolev class. The conditions employed, although standard in the literature, are combined in a different way. The results are applied to density, regression, nonlinear autoregression and Gaussian white noise models. In the latter we have also considered a loss function which is different from the usual l ² norm, namely the pointwise loss. In this case it is possible to prove that the adaptive Bayesian approach for the l ² loss is strongly suboptimal and we provide a lower bound on the rate.     

MULTI-FACTOR BALANCED BLOCK DESIGNS WITH COMPLETE ADJUSTED ORTHOGONALITY FOR ALL PAIRS OF TREATMENT FACTORS

A new series of multi-factor balanced block designs is introduced. Each of these designs has the following properties: (i) each of its k– 1 treatment factors is disposed in a cyclic or multi-cyclic balanced incomplete block design with parameters (v,b,r,k,Λ) = (a(k-l) + 1,a²(k-1) +a, ak, k, k) (a > 1); (ii) the incidence of any one of the treatment factors on any other is balanced; and (iii) after adjustment for blocks only, the relationship between any two of the treatment factors is that of adjusted orthogonality. The treatment factors are thus orthogonal to one another in the within-blocks stratum of the analysis of variance. The designs provide a benchmark with which other designs may be compared.     

Probability generating function of GPED<Subscript>2</Subscript>

The Probability generating function of a random variable which has Generalized Polya Eggenberger Distribution of the second kind (GPED ₂) is obtained. The probability density function of the range R, in random sampling from a uniform distribution on (k, l) and exponential distribution with parameter λ is obtained, when the sample size is a random variable from GPED ₂. The results of Bazargan-Lari (2004) follow as special cases.     

The estimation of R 2 and adjusted R 2 in incomplete data sets using multiple imputation

The coefficient of determination, known also as the R ², is a common measure in regression analysis. Many scientists use the R ² and the adjusted R ² on a regular basis. In most cases, the researchers treat the coefficient of determination as an index of ‘usefulness’ or ‘goodness of fit,’ and in some cases, they even treat it as a model selection tool. In cases in which the data is incomplete, most researchers and common statistical software will use complete case analysis in order to estimate the R ², a procedure that might lead to biased results. In this paper, I introduce the use of multiple imputation for the estimation of R ² and adjusted R ² in incomplete data sets. I illustrate my methodology using a biomedical example.     

New Section on Social Statistics

The analysis of a general k-factor factorial experiment having unequal numbers of observations per cell is complex. For the special case of a 2 ^k experiment with unequal numbers of observations per cell, the method of unweighted means provides a simple vehicle for analysis that requires no matrix inversion and can be used with existing software programs for the analysis of balanced data. All numerator sums of squares for testing main effects and interactions are χ² with one degree of freedom. In addition, for tests having one degree of freedom in any factorial experiment, the method of unweighted means may be modified to yield exact tests.     

Generating generalized cyclic designs with factorial balance

An algorithm for generating paired comparison factorially balanced generalized cyclic designs is described. The algorithm is based upon the 2 ⁿ ? 1 class association scheme defined by Shah (1960) for n-factor experiments. The algorithm is highly successful in achieving its objective. Firstorder designs with block size greater than two can also be obtained using the algorithm.