Open access peer-reviewed chapter

Variance Balanced Design

Written By

D.K. Ghosh

Submitted: 13 November 2021 Reviewed: 01 December 2021 Published: 30 January 2022

DOI: 10.5772/intechopen.101847

From the Edited Volume

Computational Statistics and Applications

Edited by Ricardo López-Ruiz

Chapter metrics overview

193 Chapter Downloads

View Full Metrics

Abstract

In this chapter binary, ternary and n-ary variance balanced design is constructed using balanced incomplete block, resolvable balanced incomplete block, semi regular group divisible, factorial, fractional factorial designs. Constructed variance balanced designs are with v, (v + 1), (v + 2) and (v + r) treatments. Method of construction of variance balanced designs are supported by suitable examples. It is found that all most all variance balanced designs are with high efficiency factors.

Keywords

  • incidence matrix
  • C – Matrix
  • resolvable balanced incomplete block designs
  • eigen values
  • balanced and group divisible designs

1. Introduction

In literature balanced incomplete block designs are either variance balanced (VB), efficiency balanced (EB) or pairwise balanced. Raghvarao ([1], Theorem 4.5.2) discussed that among the class of connected designs the balanced designs are the most efficient designs. A design is said to be variance balanced, if the variance of the estimate of each of the possible elementary treatment contrast is the same, i.e., if ti denotes the estimate of ith treatment effects, then Var (ti - tj) is constant for all i ≠ j.

Chakrabarti [2] gave useful concept of C – matrix of design. It is known that balanced incomplete block designs are the most efficient but do not exist for all parametric specifications, and they are equi replicated and have equal block sizes. In some situations, balanced block designs with equal replicates or unequal block size or both are needed. The variance balanced designs can have both equal and unequal number of replications and block sizes. The importance of variance balanced designs in the context of experimental material is well known, as it yields optimal designs apart from ensuring simplicity in the analysis. Many practical situations demand designs with varying block sizes (Pearce, [3], or resolvable VB designs with unequal replications Mukerjee and Kageyama [4]). Rao [5], Headyat and Federer [6], Raghavarao [7] and Puri and Nigam [8] defined that a design is said to be variance balanced, if every normalized estimable linear function of treatment effect can be estimated with the same precision. They also discussed the necessary and sufficient conditions for the existence of such designs. John [9], Jones et al. [10], Kageyama [11, 12], Kageyama et al. [13], Pal and Pal [14], Roy [15], Sinha [16, 17], Sinha and Jones [18] and Tyagi [19] gave some more methods for constructing block designs with unequal treatment replications and unequal block sizes. Khatri [20], along with a method of construction of VB designs, gave a formula to measure over-all A-efficiency of variance balanced designs. Das and Ghosh [21] gave the methods of construction of variance balanced designs with augmented blocks and treatments. Mukerjee and Kageyama [22] introduced resolvable variance balanced designs. A technique for constructing variance balanced designs, which is based on the unionizing block principle of Headayat and Federer [6], was described in Calvin [23]. Calvin and Sinha [24] extended his technique to produce designs with more than two distinct block sizes that permit fewer replications. Agarwal and Kumar [25] gave a method of construction of variance balanced designs which is associated with group divisible (GD) designs. Rao [5] observed that, if the information matrix C of a block design satisfied

C=θIv1vEvv

where, θ is non zero eigen value of C matrix, Iv is an identity matrix of order v, Evv is the matrix with v rows and v columns where, all the elements are unity, then such design is called Variance balanced designs. Since balanced incomplete block design (BIBD) satisfies this property and hence, balanced incomplete block design is a particular case of Variance balanced designs.

Das and Ghosh [21] defined generalized efficiency balanced (GEB) design, which include both VB as well as EB designs. Ghosh [26], Ghosh and Karmaker [27], Ghosh and Devecha [28], Ghosh, Divecha, and Kageyama [29], Ghosh et al. [30], Ghosh, et al. [31, 32] obtained several methods for construction of VB designs. Ghosh and Joshi [30] constructed VB design through GD design. Again, Ghosh and Joshi [33] Constructed VB Design through Triangular design. Kageyama [10] recommended the use of non-binary VB design, when binary VB designs are not available for given values of parameters. Ghosh and Ahuja [34] carried out VB design using fractional factorial designs. Agarwal and Kumar [35, 36] developed some methods of constructing ternary VB designs with v+s treatments s1, having blocks of unequal sizes, through block designs with v treatment. Ghosh, Kageyama and Joshi [37] developed Ternary VB designs using BIB and GD design. Ghosh et al. [37] further obtained more VB designs using Latin square type PBIB design. Ghosh [38] studied the robustness of variance balanced design against the loss of k treatments and one block. Ghosh et al. [39] discuss construction of VB design using factorial designs. Hedayat and Stufken [40] established a relation between pair wise balanced and variance balanced designs. Jones [41] discussed the property of incomplete block designs. Gupta and Jones [42] constructed equal replicated VB designs.

Advertisement

2. Method of construction

Method of construction of Variance balanced design with equal/unequal replication sizes and equal/unequal block sizes is carried pot in this chapter. Section 3 discusses the construction of variance balanced design using Hadamard matrix. While construction of variance balanced design using semi regular group divisible design is discussed in Section 4. Variance balanced design is constructed by augmenting n more blocks which is discussed in Section 5. Construction of variance balanced design with (v + 1) treatments using unreduced balanced incomplete block design is shown in Section 6. Section 7 discusses the construction of variance balanced design using 2n symmetrical factorial experiments. Variance balanced design is constructed using incidence matrix also, and is shown in section – 8.

Advertisement

3. Variance balanced design using Hadamard matrix

Theorem – 3.1: Equi-replicated Variance balanced design with parameters v = n – 1, b = n, r = n/2, k = {n – 1, n/2–1} and Cim = 3n4n1n2 can always be constructed from a Hadamard matrix of size n by deleting its first row and then considering rows as treatments and columns as blocks.

Proof: Consider a Hadamard matrix of size n. Delete its first row. The size of this matrix become (n – 1) x n. We replace −1 by 0, and call this matrix by N. This matrix contains (n-1) rows and n columns where, element “1” occurs n/2 times in each row, (n – 1) times in first column, and (n/2–1) times in the remaining columns. Consider matrix N as an incidence matrix of a variance balanced design, where rows are treatments and columns are blocks, so, v = n-1, b = n, r = n/2 and k = {n - 1, n/2–1} .

For variance balanced design, Cim = jbnijnmjn.j, where, i ≠ m = 1 to v.

Cim is computed as C1m = 1n1 + 1n21 = 3n4n1n2 .

We can verify that, Cim gives same constant value for each pair of treatments. Now a block design Is said to variance balanced design, if C matrix satisfies, C = θ (Iv – Evv/v), where, θ is non zero eigen value of C matrix with multiplicity (v – 1), where,

C=diag(r1,r2,,rv)NK−1N.
C=n/20..00n2..0:0:0:..:n/2nn23n4..3n43n4nn2..3n4:3n4:3n4....:nn2/n1n2

After simplification C reduces to

C=nn2n3+6n82n1n2IvEvvvE1

Where, θ = nn2n3+6n82n1n2 denotes the non-zero eigen value of C matrix with multiplicity (n – 2).

Eq. (1) satisfy the condition of variance balanced design. Hence, this is an equi-replicated and two unequal block sizes variance balanced design.

3.1 Efficiency factor of a variance balanced design

The efficiency factor of a variance balanced design is defined as

E=Vartîtm̂RBDVartîtm̂VB

Where,Vartîtm̂RBD = (2/r) σ2 = 2n/2 σ2 and

Vartîtm̂VB=2/θσ2=2nn2n3+6n82n1n2σ2
E=nn2n3+6n8nn1(n2

Example–3.1 Construct a variance balanced design from a Hadamard matrix of size 8.

Using Theorem – 3.1, we construct a variance balanced design from a Hadamard matrix of size 8 as following:

Hadamard Matrix of size 8 Incidence matrix of a Variance balanced design

1 1 1 1 1 1 1 1.

1 -1 1 -1 1 -1 1 -11 0 1 0 1 0 1 0

1 1 -1 -1 1 1 -1 -11 1 0 0 1 1 0 0

1 -1 -1 1 1 -1 -1 1 N = 10 0 1 1 0 0 1

1 1 1 1 -1 -1 -1 -11 1 1 1 0 0 0 0

1 -1 1 -1 -1 1 -1 11 0 1 0 0 1 0 1

1 1 -1 -1 -1 -1 1 11 1 0 0 0 0 1 1

1 -1 -1 1 -1 1 1 -11 0 0 1 0 1 1 0

N gives the incidence matrix of an equi-replicated and un equal block sizes variance balanced design with parameters v = 7, b = 8, r = 4, k = {7, 3}, Cim = 10/21 and information matrix,

C=40..004..0:0:0:..:44820..202048..20:20:20....:48/42

After simplification, C reduces to

C=28084I7E777=θI7E777E2

Where, θ = 103,is the non zero eigen value of C matrix with multiplicity 6. Hence, it is a variance balanced design, with tî = (1/ θ)Qi = (3/10)Qi,

Vartîtm̂VB=2/θσ2=6/10σ2.Vartîtm̂RBD=2/rσ2=24σ2

and Efficiency factor, E = 5/6. This shows that efficiency factor is very high.

Advertisement

4. Variance balanced design through semi regular group divisible designs

In this section, we discuss the construction of variance balanced design by adding the blocks of semi-regular group divisible design with its groups, provided the following conditions (i) block sizes, k = λ2, (ii) λ1 = 0 and (iii) number of groups are considered as number of blocks, are satisfied.

Theorem – 4.1 Let the parameters of a semi regular group divisible design are v, b, r, k, λ1 = 0, λ2, m and n, where k = λ2. By adding the b blocks of this semi regular group divisible design with number of groups as blocks, an equi-replicated and un-equal block sizes variance balanced design is constructed with parameters v1 = v, b1 = b + mn, r1 = r + n, k1 = {k, n} and Cim = λ2 / k or Cim = λ1 / k + n/n.

Proof: Consider a semi regular group divisible design with parameters v, b, r, k, λ1 = 0, λ2 = k, m and n, where, m denotes number of groups and n number of treatments per group. Denote N as the incidence matrix of the resulting design. Consider one group as one block. Here, there are m groups and hence, we have m more blocks. Add b blocks of the semi-regular group divisible design with its m more blocks, provided m blocks are repeated n times. Hence, v1 = v, b1 = b + mn, r1 = r + n and k1 = {k, n}. We can check, Cim = jbnijnmjn.j, where i ≠ m = 1 to v, for each pair of treatment as following. C1m = λ2k, for those pair of treatments, which occur λ2 times. Since, λ2 = k and hence, Cim = 1. Again, for those pair of treatments for which λ1 = 0, C1m = λ1k + nn = 1. For variance balanced design, Cim should be the same for each pair of treatments. Hence, λ1k + nn = λ2k . This implies that, using this method, we can construct a variance balanced design from those semi - regular group divisible designs in which (λ2 – λ1) = k holds true.

Again, a block design Is said to variance balanced design, if C matrix satisfies,

C = θ (Iv – Evv/v), where, θ is non zero eigen value of C matrix with multiplicity (v – 1), and, C = diag (r1, r2, …,rv) – N K−1N.

C=r+n0..00r+n..0:0:0:..:r+nr+kk1..11r+kk..1:1:1....:r+kk

Diagonal elements = [k (r + n) – (r + k)]/k, and off diagonal elements = − k/k = −1. After simplification, C reduces to

C=kr+nrkIvEvvvE3

Where, θ = kr+nrk denotes the non-zero eigen value of C matrix with multiplicity (v – 1).

Eq. (3) satisfy the condition of variance balanced design. Hence, this is equi replicated and two unequal block sizes variance balanced design with parameters v1 = v, b1 = b + mn, r1 = r + n, k1 = {k, n}.

4.1 Efficiency factor

The efficiency factor of a variance balanced design is defined as

E=Vartîtm̂RBDVartîtm̂VB,where,Vartîtm̂RBD=2/rσ2=2r+nσ2,
Vartîtm̂VB=2/θσ2=2kr+nrkσ2andE=kr+n)rkr+n

Example – 4.1 Construct a variance balanced design with parameters v1 = 6, b1 = 18, r1 = 8, k1 = {3, 2} from a semi regular group divisible design SR – 20, having parameters v = 6, b = 12, r = 6, k = 3, λ1 = 0, λ2 = 3, m = 3 and n = 2. Where, group is (3,2).

Three groups each with 2 treatments are (1 4), (2 5), (3 6).

Blocks of the semi-regular group divisible design, SR – 20 are.

(1 2 3), (2 4 6), (3 4 5), (1 5 6), (1 2 3), (2 4 6), (3 4 5), (1 5 6), (1 2 6), (1 3 5), (2 3 4), (4 5 6),

Using Theorem – 4.1, incidence matrix of the variance balanced design is given as.

  1 0 0 1 1 0 0 1 1 1 0 0 1 0 0 1 0 0

  1 1 0 0 1 1 0 0 1 0 1 0 0 1 0 0 1 0

N1 = 1 0 1 0 1 0 1 0 0 1 1 0 0 0 1 0 0 1

  0 1 1 0 0 1 1 0 0 0 1 1 1 0 0 1 0 0

  0 0 1 1 0 0 1 1 0 1 0 1 0 1 0 0 1 0

  0 1 0 1 0 1 0 1 1 0 0 1 0 0 1 0 0 1

N1 gives the incidence matrix of an equi-replicated and un-equal block sizes variance balanced design with parameters v1 = 6, b1 = 18, r1 = 8, k1 = {3, 2} Cim = 1 and information matrix,

C=80..008..0:0:0:..:8186..6618..6:6:6....:18/6

After simplification, C reduces to

C=6I7E777=θI6E66

Where, θ = 6is the non zero eigen value of C matrix with multiplicity 5. Hence, it is a variance balance design with tî = (1/ θ) Qi = (1/6) Qi, Vartîtm̂VB=(2/θ)σ2 = (2/6)σ2,Vartîtm̂RBD =(2/r)σ2 = (2/8) σ2, and Efficiency factor, E = 3/4. This shows that efficiency factor is very high.

Advertisement

5. Variance balanced design through augmenting n (≥ 1) blocks

In this section, variance balanced designs are obtained through balanced incomplete block design by augmenting one and more than one blocks, such that each augmented block contains each of the v treatments. The resulting design is an un - equal replicated and un - equal blocks sizes variance balanced design.

Theorem – 5.1 Let N be the incidence matrix of a balanced incomplete block design with parameters v, b, r, k and λ. Let n blocks are added with the blocks of the given balanced incomplete block design. The incidence matrix N1 defined as

N1=Nvxb1vx1

gives variance balanced design with parameters v1 = v, b1 = b + n, r1 = {(r + n), b}, k1 = {k, v}, where, N1 is the incidence matrix of Variance balanced design.

Proof: Consider a balanced incomplete block design with parameters v, b, r, k and λ, whose incidence matrix is denoted by N. Next n more blocks are augmented, hence, for resulting design, v1 = v, b1 = b + n, r1 = (r + n), k1 = {k, v}. Cim = λk + nv = λv+nkvk.

Again, a block design Is said to variance balanced design, if C matrix satisfies,

C = θ (Iv – Evv/v), where, θ is non zero eigen value of C matrix with multiplicity (v – 1) and C = diag(r1, r2, …,rv) – N K−1N.

C=r+n0..00r+n..0:0:0:..:r+nkb+nλv+nk..λv+nkλv+nkkb+n..λv+nk:λv+nk:λv+nk....:kb+n/vk

Diagonal elements are [vk(r + n) – k(b + n)]/vk, and off diagonal elements are -λv+nkvk . After simplification, C reduces to C =

vkr+nkb+nλv+nk..λv+nkλv+nkvkr+nkb+n..λv+nk:λv+nk:λv+nk....:vkr+nkb+n/vk
Finally,C=λv+nkkIvEvvvE4

where, θ = λv+nkk is the non zero eigen value of C matrix with multiplicity (v – 1). Eq. (4) satisfy the condition of variance balanced design. Hence, this is equi replicated and two unequal block sizes variance balanced design with parameters v1 = v, b1 = b + n, r1 = (r + n), k1 = {k, v}.

5.1 Efficiency factor

The efficiency factor of a variance balanced design is defined as

E=Vartîtm̂RBDVartîtm̂VB,where,Vartîtm̂RBD=2/rσ2=2r+nσ2,
Vartîtm̂VB=2/θσ2=2λv+nkkσ2=2kλv+nkandE=λv+nkkr+n.

Example – 5.1 Construct a variance balanced design with parameters v1 = 9, b1 = 15, r1 = 7, k1 = {3, 9} from a balanced incomplete block design having parameters v = 9, b = 12, r = 4, k = 3, λ = 1.

Blocks of the balanced incomplete block design are (1 2 3), (4 5 6), (7 8 9), (1 4 7), (2 5 8), (3 6 9), (1 6 8), (2 4 9), (3 5 7), (1 5 9), (2 6 7), (3 4 8). Let n = 3.

Using Theorem – 5.1, incidence matrix of the variance balanced design is given as.

  1 0 0 1 0 0 1 0 0 1 0 0 1 1 1.

  1 0 0 0 1 0 0 1 0 0 1 0 1 1 1

N1 = 1 0 0 0 0 1 0 0 1 0 0 1 1 1 1

  0 1 0 1 0 0 0 1 0 0 0 1 1 1 1

  0 1 0 0 1 0 0 0 1 1 0 0 1 1 1

  0 1 0 0 0 1 1 0 0 0 1 0 1 1 1

  0 0 1 1 0 0 0 0 1 0 1 0 1 1 1

  0 0 1 0 1 0 1 0 0 0 0 1 1 1 1

  0 0 1 0 0 1 0 1 0 1 0 0 1 1 1

N1 gives the incidence matrix of an equi replicated and un - equal block sizes variance balanced design with parameters v1 = 9, b1 = 15, r1 = 7, k1 = {3, 9} Cim = 2/3 and information matrix,

C=70..007..0:0:0:..:74518..181845..18:18:18....:45/27

After simplification, C reduces to C = 14418..1818144..18:18:18....:144/27

Finally,

C=16227I9E999=183I9E999=θI9E999E5

Where, θ = 6is the non zero eigen value of C matrix with multiplicity 8. Hence, it is a variance balance design with tî = (1/ θ) Qi = (1/6) Qi,Vartîtm̂VB =(2/θ)σ2 = (2/6) σ2. Vartîtm̂RBD= (2/r) σ2 = (2/7) σ2 and Efficiency factor, E = 6/7. This shows that efficiency factor is very high.

Advertisement

6. Variance balanced design with (v + 1) treatments

Variance balanced design with (v + 1) treatments is constructed by reinforcing one treatment in each block of a balanced incomplete block design.

6.1 Variance balanced designs with (v + 1) treatments from a series of balanced incomplete block design with parameters v, b=vC2, r = v1C21, k = 2 and λ = 1

In this section, method of the construction of variance balanced design with (v + 1) treatments is discussed. Variance balanced design with (v + 1) treatments can always be constructed through a balanced incomplete block design by reinforcing one treatment and augmenting n blocks. Let the parameters of a balanced incomplete block design are v, b=vC2, r = v1C21, k = 2 and λ = 1, provided v (r – 1) must be divisible by (k + 1) = 3. This is shown in Theorem – 6.1.

Theorem – 6.1: Let the parameters of an unreduced balanced incomplete block design are v, b=vC2, r = v1C21, k = 2 and λ = 1, whose incidence matrix is denoted by N. Let balanced incomplete block design is reinforced by one treatment up to b blocks and augmented with n blocks, such that each block contains each of the v treatments, where, n = 1, 2, …. . The incidence matrix N1 defined by

N1=NvxbEvxn11xb01xn

gives the incidence matrix of a Variance balanced design with parameters v1 = (v + 1), b1 = b + n, r1 = {r + n, b}, k = {3, v}, where, Ev x n is a matrix of v rows and n columns with elements as 1, 1 is a vector of one row and b columns, 0 is a vector of one row and n columns, provided n = vr13, n being integer.

Proof: Let us consider a unreduced balanced incomplete block designs with parameters b=vC2, r = v1C21, k = 2 and λ = 1, provided v is divisible by k. This series of balanced incomplete block design is reinforced by one more treatment and augmented by n blocks such that (v + 1)th treatment appears in each of the b blocks and each of the n more blocks contains each of the v treatment once and only once, hence, v1 = v + 1, b1 = b + n, r1 = {r + n, b}, k = {3, v} become the parameters of the resulting variance balanced design. Let us check the Cim (i ≠ m = 1 to v) value for each pair of treatments, Cim value for any pair of treatments among v treatments is computed as

Cim=13+nvE6

Again,

Cim=r3,i=1,,v,andm=v+1E7

For variance balanced design, Cim for each pair of treatment must be same and hence, from (6) and (7), 13 + nv = r3, or, nv = r13, hence, n = vr13.

Again, a block design Is said to variance balanced design, if C matrix satisfies,

C = θ (Iv – Evv/v), where, θ is non zero eigen value of C matrix with multiplicity (v – 1) and C = diag (r1, r2, …,rv) – N K−1N.

C=r+n0..00r+n..0:0:0:..:bvr+3nv+3n..v+3nv+3nvr+3n..v+3n:v+3n:v+3n....:bv/3v

Diagonal elements are (i) [3v(r + n) – (vr+3n]/3v and (ii) 2bv, and off diagonal elements are - v+3n3v. After simplification, C reduces to

C=3vr+nvr+3nv+3n..v+3nv+3n2vr+nvr+3n..v+3n:v+3n:v+3n....:2bv/3v

For variance balanced design, all the diagonal elements must be same and hence, [3v(r + n) – (vr+3n] = 2bv. This shows that one can use either of diagonal element. In this section, we use [3v(r + n) – (vr+3n] as a diagonal element.

Finally,

C=3n+2r+13IvEvvvE8

Where, θ = 3n+2r+13 is the non zero eigen value of C matrix with multiplicity v. Eq. (8) satisfy the condition of variance balanced design. Hence, this is an unequal replicated and unequal block sizes variance balanced design with parameters v1 = v, b1 = b + n, r1 = {(r + n), b}, k1 = {3, v}.

6.2 Efficiency factor

Since the resulting variance balanced design has two unequal replications and hence, there are two efficiency factors. The efficiency factor of a variance balanced design is defined as.

E1 = Vartîtm̂RBD1Vartîtm̂VB,Vartîtm̂RBD1= (2/r) σ2 = 2r+n σ2, where, ti and tm are any two treatments among v treatments, that is, i ≠ m = 1 to v.Vartîtm̂RBD2 = (1r+n + 1b) σ2, where, i ≠ m = 1 to v and m = (v + 1).

Vartîtm̂VB=2/θσ2=23n+2r+13σ2=63n+2r+1.
E1=3n+2r+13r+n.Again,E2=Vartîtm̂RBD2Vartîtm̂VB=b+r+n3n+2r+16br+n

Example – 6.1: Construct a variance balanced design with parameters v1 = 6, b1 = 15, r1 = {9, 10}, k1 = {3, 5} from a balanced incomplete block design having parameters v = 5, b = 10, r = 4, k = 2, λ = 1.

Blocks of the balanced incomplete block design are.

(1 2), (1 3), (1 4), (1 5), (2 3), (2 4), (2 5), (3 4), (3 5), (4 5). Let n = vr13 = 5. Hence, five blocks are augmented.

Using Theorem −6.1, incidence matrix of the variance balanced design is given as.

  1 1 1 1 0 0 0 0 0 0 1 1 1 1 1

  1 0 0 0 1 1 1 0 0 0 1 1 1 1 1

N1 = 0 1 0 0 1 0 0 1 1 0 1 1 1 1 1

  0 0 1 0 0 1 0 1 0 1 1 1 1 1 1

  0 0 0 1 0 0 1 0 1 1 1 1 1 1 1

  1 1 1 1 1 1 1 1 1 1 0 0 0 0 0

N1 gives the incidence matrix of an unequal replicated and unequal block sizes variance balanced design with parameters v1 = 6, b1 = 15, r1 = {9, 10}, k1 = {3, 5}, Cim = 20/15 and information matrix,

C=9000000900000090000009000000900000010352020202020203520202020202035202020202020352020202020203520202020202050/15

After simplification, C reduces to

C=100-20-20-20-20-20-20100-20-20-20-20-20-20100-20-20-20-20-20-20100-20-20-20-20-20-20100-20-20-20-20-20-20100/15
Finally,C=12015I6666=8I9E999=θI6E666E9

Where, θ = 8is the non zero eigen value of C matrix with multiplicity 5. Hence, it is a variance balanced design with tî = (1/ θ) Qi = (1/8) Qi,

Vartîtm̂VB = (2/θ)σ2 = (2/8) σ2,Vartîtm̂RBD1 = (2/r) σ2 = (2/9) σ2 and Vartîtm̂RBD2= (1r+n + 1b) σ2 = (19 + 110) σ2 = (19/90) σ2 with.

Efficiency factor, E1 = 8/9 and E2 = 38/45. This shows that efficiency factor is very high.

6.3 Variance balanced designs with (v + 1) treatments from a series of balanced incomplete block design with parameters v, b=vCk, r = v1Ck1, k and λ =v2Ck2

In Section 6.1, method of the construction of variance balanced design with (v + 1) treatments is discussed with block sizes k = 2. In this section, we have extended the method of construction of variance balanced designs with (v + 1) treatments through a balanced incomplete block design for any value of k by reinforcing one treatment and augmenting n blocks. Let the parameters of a balanced incomplete block design are v, b=vCk, r = v1Ck1, k and λ = v2Ck2, provided, v(r – λ) must be divisible by (k + 1). This is shown in Theorem – 6.2.

Theorem – 6.2 Let the parameters of a balanced incomplete block design are v, b=vCk, r = v1Ck1, k and λ = v2Ck2. whose incidence matrix is denoted by N. Let balanced incomplete block design is reinforced by one treatment up to b blocks and augmented with n blocks, such that each block contains each of the v treatments, where, n = 1, 2, …. . The incidence matrix N1 defined by

N1=NvxbEvxn11xb01xn

gives the incidence matrix of a Variance balanced design with parameters v1 = (v + 1), b1 = b + n, r1 = {r + n, b}, k1 = {(k + 1), v}, where, Ev x n is a matrix of v rows and n columns with elements as 1, 1 is a vector of one row and b columns, 0 is a vector of one row and n columns, provided n =vrλk+1, n being integers.

Proof: Let us consider a series of balanced incomplete block designs with parameters v,b=vCk, r = v1Ck1, k and λ = v2Ck2, provided n is divisible byvrλk+1 . This series of balanced incomplete block design is reinforced by one more treatment and augmented by n blocks such that (v + 1)th treatment appears in each of the b blocks and each of the n more blocks contains each of the v treatment once and only once, hence, v1 = v + 1, b1 = b + n, r1 = {r + n, b}, k1 = {(k + 1), v} are the parameters of the resulting variance balanced design. Let us check the Cim (i ≠ m = 1 to v) value for each pair of treatments. Cim value for any pair of treatments among v treatments is computed as

Cim=1k+1+nvE10
Again,Cim=rk+1,i=1,,vandm=v+1E11

For variance balanced design, Cim for each pair of treatment, must be same and hence, from (10) and (11), 1k+1 + nv = rk+1, or, nv = rλk+1, hence, n = vrλk+1.

Again, a block design Is said to variance balanced design, if C matrix satisfies,

C = θ (Iv – Evv/v), where, θ is non zero eigen value of C matrix with multiplicity (v – 1) and C = diag(r1, r2, …,rv) – N K−1N.

C=r+n0..00r+n..0:0:0:..:bvr+nk+1λv+nk+1..λv+nk+1λv+nk+1vr+nk+1..λv+nk+1:λv+nk+1:λv+nk+1....:bv/vk+1

Diagonal elements are (i) [v(k + 1)(r + n) – (vr+nk+1]/v(k + 1) and (ii) kbv and off diagonal elements are - λv+nk+1vk+1,. After simplification, C reduces to

C=vkn+r+nvk1λv+nk+1..λv+nk+1λv+nk+1vkn+r+nvk1..v+nk+1:λv+nk+1:λv+nk+1....:kbv/vk+1

For variance balance design, all the diagonal elements must be same and hence,vkn+r+nvk1 = kbv. This shows that we can use either of diagonal element. In this section, we usedvkn+r+nvk1 as a diagonal element.

Finally,C=nk+1+kr+λk+1IvEvvvE12

Where, θ = nk+1+kr+λk+1 is the non zero eigen value of C matrix with multiplicity v. Eq. (12) satisfy the condition of variance balanced design. Hence, this is an unequal replicated and unequal block sizes variance balanced design with parameters v1 = v, b1 = b + n, r1 = {(r + n), b}, k1 = {3(k + 1) v}.

6.3.1 Efficiency factor of this variance balanced design

Since the resulting variance balanced design is a two unequal replicated design and hence, there are two efficiency factors. The efficiency factor of a variance balanced design is defined as.

E1 = Vartîtm̂RBD1Vartîtm̂VB, Vartîtm̂RBD1 = (2/r) σ2 = 2r+n σ2, where ti and tm are any two treatments among v treatments, i ≠ m = 1 to v.Vartîtm̂RBD2 = (1r+n + 1b) σ2, where, i ≠ m = 1 to v and m = (v + 1).

Vartîtm̂VB=2/θσ2=2nk+1+kr+λk+1σ2=2k+1nk+1+kr+λ.
E1=nk+1+kr+λr+nk+1andE2=nk+1+kr+λb+r+n2br+nk+1

Example – 6.2: Construct a variance balanced design with parameters v1 = 6, b1 = 15, r1 = {9, 10}, k1 = {3, 5} from a balanced incomplete block design having parameters v = 6, b = 20, r = 10, k = 3, λ = 4.

Blocks of the balanced incomplete block design are.

(1 2 3), (1 2 4), (1 2 5), (1 2 6), (1 3 4), (1 3 5), (1 3 6), (1 4 5), (1 4 6), (1 5 6), (2 3 4), (2 3 5), (2 3 6), (2 4 5), (2 4 6), (2 5 6), (3 4 5), (3 4 6), (3 5 6), (4 5 6). Let n = vrλk+1 = 9. Hence, nine blocks are augmented. Using Theorem – 6.2, incidence matrix of the variance balanced design is given as.

  1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1.

  1 1 1 1 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 1 1 1 1 1 1 1 1 1.

N1 = 1 0 0 0 1 1 1 0 0 0 1 1 1 0 0 0 1 1 1 0 1 1 1 1 1 1 1 1 1.

  0 1 0 0 1 0 0 1 1 0 1 0 0 1 1 0 1 1 0 1 1 1 1 1 1 1 1 1 1.

  0 0 1 0 0 1 0 1 0 1 0 1 0 1 0 1 1 0 1 1 1 1 1 1 1 1 1 1 1.

  0 0 0 1 0 0 1 0 1 1 0 0 1 0 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1.

  1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0.

N1 gives the incidence matrix of an unequal replicated and un - equal block sizes variance balanced design with parameters v1 = 7, b1 = 29, r1 = {19, 20}, k1 = {4, 6}; Cim = 10/4 and information matrix,

C=190..0019..0:0:0:..:204830..303048..30:30:30....:60/12

After simplification, C reduces to C = 180303030180..30:30:30:..:180/12

Finally,C=21012I7E777=352I7E777=θI7E777E13

Where, θ = 35/2 is the non zero eigen value of C matrix with multiplicity 6. Hence, it is a variance balanced design withtî = (1/ θ) Qi = (2/35) Qi, Vartîtm̂VB = (2/θ)σ2 = (4/35) σ2,Vartîtm̂RBD1 = (2/r) σ2 = (2/19) σ2 and Vartîtm̂RBD2= (1r+n + 1b) σ2 = (19 + 110) σ2 = (39/380) σ2 with efficiency factor, E1 = 35/38 and E2 = 273/304. This shows that efficiency factor is very high.

Nonexistence of variance balanced design by reinforcing (v + t) treatments.

Variance balanced design cannot be constructed from a balanced incomplete block design by reinforcing (v + t) treatments, t = 1, 2, …

Because for variance balanced design Cim must be same for each pair of treatments. In this case Cim = λk+1, i ≠ m = 1 to v and Ci (v + 1) = rk+1, where, t = 1. As per the condition of variance balanced design, rk+1= λk+1, which is not possible as (r – λ) > 0. If two treatments are reinforced, then rk+1= λk+1 = bk+1 must holds true, but b ≠ r ≠ λ.

Advertisement

7. Variance balanced design using 2n symmetrical factorial experiments

This section discusses the construction of variance balanced design using 2n symmetrical factorial experiment.

Theorem: 7.1 Let us consider a 2n factorial experiment. By deleting the control treatment and all the main effects, equi-replicated and unequal block sizes variance balanced design is obtained with parameters v = n, b = 2n – n – 1, r = 2n − 1, k = {2, 3, 4…, n}.

Proof: consider a 2n treatment combination of a 2n factorial experiment. Delete its control treatment and all main effects. Consider each treatment combination as one block. So, we have (2n – 1 – n) blocks with v (= n) treatments, as factors are considered as treatments. Consider this matrix as an incidence matrix of a block design, whose all elements are either zero or 1 only. So, the design is binary.

Since each treatment is repeated (2n − 1– 1) times, design is equi replicated and unequal block sizes with k = {2, 3,…, n}. Let the incidence matrix of the block design is given by

N=011111001110101011111.

The incidence matrix of the block design is a variance balanced, if the C matrix of the block design satisfy C = θIv1vEVV, where, θ is non – zero eigen value of C matrix.

C=2n11002n1100002n11YXXYXXXXY

where, Y =n112+n123+n134+..+1n, and

X = 12+n23+n34++nn1n and r = 2n – 1 -1

After simplification, C reduces to

C=2n11YXX2n11YXXXX2n11Y

Finally,

C=2n11Y+XIvEvvvE14

Where, θ = (2n – 1 –1 - Y + X) is the non-zero eigen value of C matrix with multiplicity (v – 1). Eq. (14) satisfy the condition of variance balanced design. Hence, this is an equal replicated and unequal block sizes variance balanced design with parameters v = n, b = 2n – n – 1, r = 2n − 1 - 1, k = {2, 3, 4…, n.}. Efficiency factor, E = 2n11Y+X2n11 .

Example: 7.1. Construct a variance balanced design with parameters v = 4, b = 11, r = 7, k = {2, 3, …, n.}

Using Theorem 7.1, incidence matrix of the variance balanced design is given by

N=00001111111011100011111011011001111011010101
C=700007000000000733171717173317171717171733171733/12

After simplification, C reduces to

C=51171717175117171717171751171751/12

Finally,

C=6812I4E444=173I4E444=θI4444E15

Where, θ = 17/3 is the non zero eigen value of C matrix with multiplicity 3. Hence, it is a variance balanced design with tî = (1/ θ) Qi = (3/17) Qi, Vartîtm̂VB = (2/θ)σ2 = (6/17) σ2, Vartîtm̂RBD= (2/r) σ2 = (2/7) σ2 and efficiency factor, E = 17/21.

This result is due to Ghosh, Sinojia and Ghosh (2018).

Advertisement

8. Construction of variance balanced designs using some incidence matrix

Theorem 8.1: Let In denotes the identity matrix of order n,jn is a column Vector of one, 0n is the row vectors having all elements zero. An incidence matrix N defined as N=jn01×n2InEnxn/2

gives the incidence matrix of a Variance balanced designs, with parameters

v=n+1,b=n+n2,r=n1+n2 and k=2n, where, n is even.

Proof: Proof is obvious.

Example 8.1. Let = 6. So, the incidence matrix N using Theorem 8.1 is given by

N=111111000100000111010000111001000111000100111000010111000001111
C=60000000400000004000000040000000400000004000000046111111121111111211111112111111121111111211111112/2

After simplification, C reduces to

C=61111111611111116111111161111111611111116111111166/2

Finally,

C=72I7E777=72I7E777=θI7E777E16

where, θ = 7/2 is the non zero eigen value of C matrix with multiplicity 6. Hence, it is a variance balanced design with parameters v = 7, b = 9, r = {6,4}, k = {2, 6}. tî = (1/θ) Qi = (2/7) Qi; Vartîtm̂VB = (2/θ)σ2 = (4/7) σ2, Vartîtm̂RBD1= (2/r) σ2 = (2/4) σ2, Vartîtm̂RBD2= (16 + 14) σ2 = (5/12) σ2 with efficiency factor, E1 = 7/8 and efficiency factor, E2 = 35/48.

This result is due to Ghosh, Sinojia and Ghosh (2018).

Advertisement

9. Conclusions

In this chapter, we have constructed Variance balanced designs using balanced incomplete block, group divisible, resolvable semi - regular group divisible, symmetrical factorial and fractional factorial designs. It is observed that efficiency factor of all most all variance balanced design is high. Variance balanced designs constructed in sections 3 to 6 are new and extended methods, while Section 7 and 8, discuss the review work of Ghosh, Sinojia and Ghosh (2018).

Advertisement

Acknowledgments

Author is very much grateful to University Grants Commission, New Delhi, for providing me an opportunity to work as UGC BSR Faculty Fellow. Author is also thankful to referees for suggesting the important ideas in improving the present chapter.

References

  1. 1. Raghavarao D. Constructions and Combinatorial Problems in Design of Experiments. New York: Wiley; 1971
  2. 2. Chakrabarti MC. On the C-matrix in design of experiments. Journal of the Indian Statistical Association. 1963;1:8-23. 8C23
  3. 3. Pearce SC. Experimenting with blocks of natural size. Biometrics. 1964;20(4):699-706. 699C706. DOI: 10.2307/2528123
  4. 4. Kageyama S. Reduction of associate classes for block designs and related combinatorial arrangements. Hiroshima Mathematical Journal. 1974;4:527C618
  5. 5. Rao VR. A note on balanced designs. The Annals of Mathematical Statistics. 1958;29(1):290-294. 290C294. DOI: 10.1214/aoms/1177706729
  6. 6. Hedayat A, Federer WT. Pairwise and variance balanced incomplete block design. Annals of the Institute of Statistical Mathematics. 1974;26(1):331-338. DOI: 10.1007/BF02479828
  7. 7. Raghavarao D. On balanced unequal block designs. Biometrika. 1962;49(3/4):561-562. 561C562. DOI: 10.1093/biomet/49.3-4.561
  8. 8. Puri PD, Nigam AK. Balanced block designs. Communications in Statistics – Theory and Methods. 1977;6(12):1171-1179. 1171C1179. DOI: 10.1080/03610927708827560
  9. 9. John PWM. Balanced designs with unequal number of replications. Annals of Mathematical Statistics. 1964;35:897C899. DOI: 10.1214/aoms/1177703597
  10. 10. Jones B, sinha, K., & Kageyama, S. Further equireplicated variance balanced designs with unequal block sizes. Utilitas Mathematica. 1987;32:5C10
  11. 11. Kageyama S. Construction of balanced block designs. Utilitas Mathematica. 1976;9:209C229
  12. 12. Kageyama S. Existence of variance-balanced binary designs with fewer experimental units. Statistics & Probability Letters. 1988;7(1):27-28. 27C28. DOI: 10.1016/0167-7152(88)90083-1
  13. 13. Kulshreshtha AC, Dey A, Saha GM. Balanced designs with unequal replications and unequal block sizes. Annals of Mathematical Statistics. 1972;43:1342-1345
  14. 14. Pal S, Pal S. Nonproper variance balanced designs and optimality. Communications in Statistics - Theory and Methods. 1988;17:1685C1695. DOI: 10.1080/03610928808829706
  15. 15. Roy J. On the efficiency factor of block designs. Sankhya. 1958;19:181C188
  16. 16. Sinha K. A note on equireplicated balanced block designs from BIB designs. Journal of the Indian Society of Agricultural Statistics. 1988;42:150C153
  17. 17. Sinha K. Some new equireplicated balanced block designs. Statistics & Probability Letters. 1989;8:89. DOI: 10.1016/0167-7152(89)90089-8
  18. 18. Sinha K, Jones B. Further equireplicated balanced block designs with unequal block sizes. Statistics & Probability Letters. 1988;6:229C330
  19. 19. Tyagi BN. On a class of variance balanced block designs. Journal of Statistical Planning and Inference. 1979;3:333-336. DOI: 10.1016/0378-3758(79)90029-6
  20. 20. Khatri CG. A note on variance balanced designs. Journal of Statistical Planning and Inference. 1982;6(2):173-177. 173C177. DOI: 10.1016/0378-3758(82)90086-6
  21. 21. Das MN, Ghosh DK. Balancing incomplete block designs. Sankhyā: The Indian Journal of Statistics, Series B. 1985;47(1):67-77. 67C77
  22. 22. Mukerjee R, Kageyama S. On resolvable and affine resolvable variance-balanced designs. Biometrika. 1985;72(1):165-172. 165C172. DOI: 10.1093/biomet/72.1.165
  23. 23. Calvin JA. A new class of variance balanced designs. Journal of Statistical Planning and Inference. 1986;14(2–3):251-254. DOI: 10.1016/0378-3758(86)90162-X
  24. 24. Calvin JA, Sinha K. A method of constructing variance balanced designs. Journal of Statistical Planning and Inference. 1989;23(1):127-131. DOI: 10.1016/0378-3758(89)90045-1
  25. 25. Agarwal GG, Kumar S. On a class of variance balanced designs associated with GD designs. Calcutta Statistics Association Bulletin. 1984;33(3–4):187-190. 187C190
  26. 26. Ghosh DK. A series of generalized efficiency balanced designs. Gujarat Statistical Review. 1988;15(1):33-38. 33C38
  27. 27. Ghosh DK, Karmokar PK. Some series of efficiency balanced designs. Australian & New Zealand Journal of Statistics. 1988;30(1):47-51. 47C51
  28. 28. Ghosh DK, Divecha J. Construction of variance balanced designs. TEC.R.No.2, Dept. of Maths. and Stat., Saurashtra Uni. 1989:1-13. 1C13
  29. 29. Ghosh DK, Divecha J, Kageyama S. Equi - replicate variance balanced designs from group divisible designs. Journal of the Japan Statistical Society. 1991;21(1):205-209. 205C209. DOI: 10.1111/j.1467-842X.1988.tb00610.x
  30. 30. Ghosh DK, Joshi K. Note on construction of variance balanced designs through group divisible designs. Utilitas Mathematica. 1991;39:249-253. 249C253
  31. 31. Ghosh DK, Karmokar PK, Divecha J. On Comparison of ternary efficiency balanced and ternary variance balanced designs. Utilitas Mathematica. 1991;40:25-27
  32. 32. Ghosh DK, Anita S, Kageyama S. Construction of variance-balanced and efficiency balanced ternary block designs. Journal of Japan Statistical Association Soc. 1994;24(2):201-208
  33. 33. Ghosh DK, Joshi K. Construction of variance balanced designs through triangular PBIB designs. Calcutta Statistical Association Bulletin. 1995;45(1–2):111-118
  34. 34. Ghosh DK, Sangeeta A. On variance balanced designs. Journal of Modern Applied Statistical Methods. 2017;16(2):124-137
  35. 35. Agarwal GG, Kumar S. A note on construction of variance balanced designs. Journal of the Indian Society of Agricultural Statistics. 1985;37(2):181-183. 181C183
  36. 36. Agarwal GG, Kumar S. Construction of balanced ternary designs. Australian and Newziland Journal of Statistics. 1986;28(2):251-255. 251C255. DOI: 10.1111/j.1467-842X.1986.tb00605.x
  37. 37. Ghosh DK, Joshi K, Kageyama S. Ternary variance balanced designs through BIB and GD designs. Journal of the Japan Statistical Society. 1993;23(1):75-81. 75C81
  38. 38. Ghosh DK. Robustness of Variance balanced designs against the loss of one block. Indian Journal of Statistics and Applications. 2012;1:17-30. 17C30
  39. 39. Ghosh DK, Sinojia CN, Ghosh S. On a class of variance balanced designs. International Journal of Statistics and Probability. 2018;7(3):112-120
  40. 40. Hedayat A, Stufken JA. Relation between pair wise balanced and variance balanced block designs. Journal of the American Statistical Association. 1989;84(407):753-755. DOI: 10.1080/01621459.1989.10478830
  41. 41. Jones RM. On a property of incomplete blocks. Journal of the Royal Statistical Society: Series B. 1959;21:172-179. Available from: http://www.jstor.org/stable/2983939
  42. 42. Gupta SC, Jones B. Equi-replicate balanced block designs with unequal block sizes. Biometrika. 1983;70(2):443-440. 443C440. DOI: 10.1093/biomet/70.2.433

Written By

D.K. Ghosh

Submitted: 13 November 2021 Reviewed: 01 December 2021 Published: 30 January 2022