Optimality aspects of 3-concurrence most balanced designs

Takcuchi (1961,1963) established E-optimality of Group Divisible Designs (GDDs) with λ₂=λ₁+1. Much later, Cheng (1980) and Jacroux (1980,1983) demonstrated E-optimality property of the GDDs with n=2,λ₁=λ₂+1 or with m=2,λ₂=λ₁+2. The purpose of this paper is to provide a unified approach for identifying certain classes of designs as E-optimal. In the process, we come up with a complete characterization of all E-optimal designs attaining a specific bound for the smallest non-zero eigenvalue of the underlying C-matrices. This establishes E-optimality of a class of 3-concurrence most balanced designs with suitable intra- and inter-group balancing. We also discuss the MV-optimality aspect of such designs.     

Super-simple group divisible designs with block size 4 and index 9

A design is said to be super-simple if the intersection of any two blocks has at most two elements. In statistical planning of experiments, super-simple designs are the ones providing samples with maximum intersection as small as possible. Super-simple GDDs are useful in constructing super-simple BIBDs. The existence of super-simple (4,λ)‐GDDs has been determined for λ=2-6. In this paper, we investigate the existence of a super-simple (4,9)-GDD of group type g^u and show that such a design exists if and only if u≥4, g(u−2)≥18 and u(u−1)g²≡0 (mod 4).     

Aliasing in 2n–k fractions in the absence of within group interactions

It is often known in advance that certain subsets of factors act independently upon a response. Such information can be used to estimational advantage by aliasing low-order effects with such zero interactions. We find the best 2^n–k fractions for the case when the factors can be partitioned into two classes such that non-zero interactions may exist only between classes but not within a class.     

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.     

Addition of runs to a two-level supersaturated design

The purpose of this article is to introduce a new class of extended E(s²)-optimal two level supersaturated designs obtained by adding runs to an existing E(s²)-optimal two level supersaturated design. The extended design is a union of two optimal SSDs belonging to different classes. New lower bound to E(s²) has been obtained for the extended supersaturated designs. Some examples and a small catalogue of E(s²)-optimal SSDs are also included.     

Partition of a query set into minimal number of subsets having consecutive retrieval property

This paper is concerned with information retrieval. The basic problem is how to store large masses of data in such a way that whenever information regarding some particular aspect of the data is needed, such information is easily and efficiently retrieved. Work in this field is thus very important for organizations dealing with large classes of data.The consecutive retrieval (C-R) property defined by S.P. Ghosh is an important relation between a set of queries and a set of records. Its existence enables the design of information retrieval system with a minimal search time and no redundant storage in that the records can be organized in such a way that those pertinent to any query are stored in consecutive storage locations. The C-R property, however, can not exist between every arbitrary query set and every record set.A subset of the query set Q having the C-R property is called a C-R subset and a C-R subset having the maximum cardinality is called the maximal C-R subset. A partition of Q is called the C-R partition if every subset has the C-R property. A C-R partition with minimum number of subsets is called the minimal C-R partition. With respect to the set of all binary queries and the set of all binary records, it is shown that the maximal cardinality of a C-R subset is 2l-1 where l is the number of attributes concerned. A combinatorial characterization of a maximal C-R subset is also given. A lower bound on the number of subsets in a C-R partition and several examples which attain the lower bound are given. A general procedure for obtaining a minimal C-R partition which attains the lower bound is given provided the number of attributes is even.     

A new ranked set sampling protocol for the signed rank test

In this paper, we construct a new ranked set sampling protocol that maximizes the Pitman asymptotic efficiency of the signed rank test. The new sampling design is a function of the set size and independent order statistics. If the set size is odd and the underlying distribution is symmetric and unimodal, then the new sampling protocol quantifies only the middle observation. On the other hand, if the set size is even, the new sampling design quantifies the two middle observations. This data collection procedure for use in the signed rank test outperforms the data collection procedure in the standard ranked set sample. We show that the exact null distribution of the signed rank statistic W_RSS⁺ based on a data set generated by the new ranked set sample design for odd set sizes is the same as the null distribution of the simple random sample signed rank statistic W_SRS⁺ based on the same number of measured observations. For even set sizes, the exact null distribution of W_RSS⁺ is simulated.     

Fast nonparametric classification based on data depth

A new procedure, called D D α-procedure, is developed to solve the problem of classifying d-dimensional objects into q ≥ 2 classes. The procedure is nonparametric; it uses q-dimensional depth plots and a very efficient algorithm for discrimination analysis in the depth space [0,1] ^q . Specifically, the depth is the zonoid depth, and the algorithm is the α-procedure. In case of more than two classes several binary classifications are performed and a majority rule is applied. Special treatments are discussed for ‘outsiders’, that is, data having zero depth vector. The D Dα-classifier is applied to simulated as well as real data, and the results are compared with those of similar procedures that have been recently proposed. In most cases the new procedure has comparable error rates, but is much faster than other classification approaches, including the support vector machine.     

Nested Hadamard difference sets

A Hadamard difference set (HDS) has the parameters (4N², 2N² − N, N² − N). In the abelian case it is equivalent to a perfect binary array, which is a multidimensional matrix with elements ±1 such that all out-of-phase periodic autocorrelation coefficients are zero. We show that if a group of the form H × Z²_p^r contains a (hp^2r, √hp^r(2√hp^r − 1), √hp^r(√hp^r − 1)) HDS (HDS), p a prime not dividing |H| = h and p^j ≡ −1 (mod exp(H)) for some j, then H × Z²_p^t has a (hp^2t, √hp^t(2√hp^t − 1), √hp^t(√hp^t − 1)) HDS for every 0⩽t⩽r. Thus, if these families do not exist, we simply need to show that H × Z²_p does not support a HDS. We give two examples of families that are ruled out by this procedure.     

Multistage ranked set sampling

The superiority of using ranked set sampling, for estimating the mean of a population, over simple random sampling, is well established. This technique is useful when visual ordering of a small set of size (m) can be done easily and fairly accurately, but exact measurement of an observation is difficult and expensive. It is noted that for many distributions, an increase in the efficiency of ranked set sampling can be achieved by increasing the set size m. However, in practice, m should be kept very small so that visual ranking errors will not destroy the gain in efficiency. In this paper, multistage ranked set sampling is considered as a generalization of ranked set sampling, that results in an increase of the efficiency for fixed value of m. Steady state efficiency, the limiting efficiency as the number of stages approaches infinity, varies from one distribution to another. It is shown that this efficiency is always larger than 1, close to m² for symmetric distributions and equal to m² for the uniform distribution. Some real applications of the technique are discussed. Data on olive yield of olive trees is collected to illustrate the technique.     

Existence of certain class of duadic group algebra codes

Duadic codes are defined in terms of idempotents of a group algebra GF(q)G, where G is a finite group and gcd(q,|G|)=1. Under the conditions of (1) q=2^m, and (2) the idempotents are taken to be central and (3) the splitting is μ₋₁, we show that such duadic codes exist if and only if q has odd-order modulo |G|.     

One-factor-at-a-time designs of resolution V

Two classes of designs of resolution V are constructed using one-factor-at-a-time techniques. They facilitate sequential learning and are more economical in run size than regular 2_V^k−p designs. Comparisons of D efficiency with other designs are given to assess the suitability of the proposed designs. Their D_s efficiencies can be dramatically improved with the addition of a few runs.     

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}.     

Two-dimensional density estimation using smooth invertible transformations

We investigate the problem of estimating a smooth invertible transformation f when observing independent samples X₁,…,X_n∼P°f where P is a known measure. We focus on the two-dimensional case where P and f are defined on R². We present a flexible class of smooth invertible transformations in two dimensions with variational equations for optimizing over the classes, then study the problem of estimating the transformation f by penalized maximum likelihood estimation. We apply our methodology to the case when P°f has a density with respect to Lebesgue measure on R² and demonstrate improvements over kernel density estimation on three examples.     

On computer generation of balanced arrays

This paper provides an algebraic (and hence computing) procedure for generation of balanced arrays having two symbols, m rows, specified minimum and maximum column weights, arbitrary strength t≦m, and index set parameters μ^t₁, μ^t₂,…, μ^t_t. μ^t₀ is unspecified, and calculated as part of the algorithm, although the procedure for specifying it is straightforward and can be used if desired. Array generation is herein reduced to finding integral solutions to a linear programming problem. It is shown that the integral solutions of the system of equations comprise all balanced arrays with the given set of parameters.A computing algorithm is provided which constructs the system of equations to be solved; it has been interfaced with a standard linear programming package to provide some preliminary results.Additional algorithms whose development should result in substantial decreases in computing costs are discussed.     

Theory of factorial designs of the parallel flats type. I: the coefficient matrix

Let GF(s) be the finite field with s elements.(Thus, when s=3, the elements of GF(s) are 0, 1 and 2.)Let A(r×n), of rank r, and c_i(i=1,…,f), (r×1), be matrices over GF(s). (Thus, for n=4, r=2, f=2, we could have A=[¹¹¹⁰₀₁₂₁], c₁=[¹₀], c₂=[⁰₂].) Let T_i (i=1,…,f) be the flat in EG(n, s) consisting of the set of all the s^n?r solutions of the equations At=c_i, wheret′=(t₁,…,t_n) is a vector of variables.(Thus, EG(4, 3) consists of the 3⁴=81 points of the form (t₁,t₂,t₃,t₄), where t's take the values 0,1,2 (in GF(3)). The number of solutions of the equations At=c_i is s^n?r, where r=Rank(A), and the set of such solutions is said to form an (n?r)-flat, i.e. a flat of (n?r) dimensions. In our example, both T₁ and T₂ are 2-flats consisting of 3^4?2=9 points each. The flats T₁,T₂,…,T_f are said to be parallel since, clearly, no two of them can have a common point. In the example, the points of T₁ are (1000), (0011), (2022), (0102), (2110), (1121), (2201), (1212) and (0220). Also, T₂ consists of (0002), (2010), (1021), (2101), (1112), (0120), (1200), (0211) and (2222).) Let T be the fractional design for a sⁿ symmetric factorial experiment obtained by taking T₁,T₂,…,T_f together. (Thus, in the example, 3⁴=81 treatments of the 3⁴ factorial experiment correspond one-one with the points of EG(4,3), and T will be the design (i.e. a subset of the 81 treatments) consisting of the 18 points of T₁ and T₂ enumerated above.)In this paper, we lay the foundation of the general theory of such ‘parallel’ types of designs. We define certain functions of A called the alias component matrices, and use these to partition the coefficient matrix X (n×v), occuring in the corresponding linear model, into components X._j(j=0,1,…,g), such that the information matrix X is the direct sum of the X′._jX._j. Here, v is the total number of parameters, which consist of (possibly μ), and a (general) set of (geometric) factorial effects (each carrying (s?1) degrees of freedom as usual). For j≠0, we show that the spectrum of X′._jX._j does not change if we change (in a certain important way) the usual definition of the effects. Assuming that such change has been adopted, we consider the partition of the X._j into the X_ij (i=1,…,f). Furthermore, the X_ij are in turn partitioned into smaller matrices (which we shall here call the) X_ijh. We show that each X_ijh can be factored into a product of 3 matrices J, ζ (not depending on i,j, and h) and Q(j,h,i)where both the Kronecker and ordinary product are used. We introduce a ring R using the additive groups of the rational field and GF(s), and show that the Q(j,h,i) belong to a ring isomorphic to R. When s is a prime number, we show that R is the cyclotomic field. Finally, we show that the study of the X._j and X′._jX._j can be done in a much simpler manner, in terms of certain relatively small sized matrices over R.     

Bayesian inference for Gaussian excursion set generated Cox processes with set-marking

This work considers spatial Cox point processes where the random intensity is defined by a random closed set such that different point intensities appear in the two phases formed by the random set and its complement. The point pattern is observed as a set of point coordinates in a bounded region W?? ^d together with the information on the phase of the location of each point. This phase information, called set-marking, is not a representative sample from the random set, and hence it cannot be directly used for deducing properties of the random set. Excursion sets of continuous-parameter Gaussian random fields are applied as a flexible model for the random set. Fully Bayesian method and Markov chain Monte Carlo (MCMC) simulation is adopted for inferring the parameters of the model and estimating the random set. The performance of the new approach is studied by means of simulation experiments. Further, two forestry data sets on point patterns of saplings are analysed. The saplings grow in a clear-cut forest area where, before planting and natural seeding, the soil has been mounded forming a blotched soil structure. The tree densities tend to be different in the tilled patches and in the area outside the patches. The coordinates of each sapling have been measured and it is known whether this location is in a patch or outside. This example has been a motivation for the study.     

MBMUDs: a combinatorial extension of BIBDs showing good optimality behaviour

The construction of a balanced incomplete block design (BIBD) is formulated in terms of combinatorial optimization by defining a cost function that reaches its lower bound on all and only those configurations corresponding to a BIBD. This cost function is a linear combination of distribution measures for each of the properties of a block design (number of plots, uniformity of rows, uniformity of columns, and balance). The approach generalizes naturally to a super-class BIBDs, which we call maximally balanced maximally uniform designs (MBMUDs), that allow two consecutive values for their design parameters [r,r+1;k,k+1;λ,λ+1]. In terms of combinatorial balance, MBMUDs are the closest possible approximation to BIBDs for all experimental settings where no set of admissible parameters exists. Thus, other design classes previously proposed with the same approximation aim—such as RDGs, SRDGs and NBIBDs of type I—can be viewed as particular cases of MBMUDs. Interestingly, experimental results show that the proposed combinatorial cost function has a monotonic relation with A- and D-statistical optimality in the space of designs with uniform rows and columns, while its computational cost is much lower.     

New classes of D-optimal edge designs

Edge designs are screening experimental designs that allow a model-independent estimate of the set of relevant variables, thus providing more robustness than traditional designs. In this paper we construct new classes of D-optimal edge designs. This construction uses weighing matrices of order n and weight k together with permutation matrices of order n to obtain D-optimal edge designs. One linear and one quadratic simulated screening scenarios are studied and compared using linear regression and edge designs analysis.