Open access peer-reviewed chapter

Constant Block-Size with Constant Sum-Block Partially Balanced Incomplete Block Designs

Written By

Babatunde L. Adeleke, Gabriel O. Adebayo and Kazeem A. Osuolale

Submitted: 23 January 2024 Reviewed: 24 January 2024 Published: 10 July 2024

DOI: 10.5772/intechopen.1004363

From the Edited Volume

Response Surface Methods - Theory, Applications and Optimization Techniques

Valter Silva and João Sousa Cardoso

Chapter metrics overview

18 Chapter Downloads

View Full Metrics

Abstract

This chapter delves into the introduction, methodology, and application of Constant Block-Size with Constant Sum-Block Partially Balanced Incomplete Block Designs (CBS-CSB PBIBDs) as an innovative class of experimental designs. The study defines key concepts, including constant block-size, constant block-sum, and constant sum-block, and outlines conditions for implementing CBS-CSB PBIBDs. The research presents experimental designs for various scenarios, such as t = 21, k = 3; t = 15, k = 3; t = 27, k = 3, and illustrates their respective cases through line charts. Notably, the study demonstrates that as the number of replicates (r) increases, the efficiency factors also increase, and when r = 1, efficiency (E) equals 0. In conclusion, CBS-CSB PBIBDs emerge as a valuable tool for experimental prioritization, offering a new perspective on partially balanced incomplete block designs. The research highlights the significance of CBS-CSB PBIBDs across diverse fields of experimentation. The findings equip readers with the ability to define, design, and analyze experiments using CBS-CSB PBIBDs, while acknowledging challenges and limitations associated with this approach. This study contributes to the broader understanding of experimental design and provides a foundation for future research directions in this innovative field.

Keywords

  • experimental design
  • partially balanced incomplete block designs
  • constant block-size
  • constant sum-block
  • prioritization of treatments

1. Introduction

Constant Block-Size with Constant Sum-Block Partially Balanced Incomplete Block Designs, CBS-CSB PBIBDs is a new class of partially balanced incomplete block designs that is use when it requires prioritizing the treatments combination. This new class of the design is significant in every aspect of experiment in physical, life, clinical, engineering, education, social science, psychology, pharmaceutical etc. once there is need to prioritize the treatment combination then CBS-CSB PBIBDs will be used.

The goal of this chapter is to define constant block-size partially balanced incomplete block designs and constant sum-block partially balanced incomplete block designs, design experiments using constant block-size with constant sum-block partially balanced incomplete block designs, CBS-CSB PBIBDs, state the Conditions for Constant Block-Size with Constant Sum-Block Partially Balanced Incomplete Block Designs, CBS-CSB PBIBDs, perform simple analysis on CBS-CSB PBIBDs and state some challenges and limitation on CBS-CSB PBIBDs.

In the remaining part of this chapter, we will discuss the background and theory, design principles and methodology, data analysis, discussion of results and conclusion of study’s findings.

1.1 Background and theory

The optimum design is balanced incomplete block designs (BIBDs) which gives 100% efficiency factor, e = 1. But the shortcoming of BIBDs leads to the introduction of partially balanced incomplete block designs (PBIBDs) introduced by the authors in [1]. This chapter introduces partially balanced incomplete block designs and its Efficiency Factors but its limitation was that it took into consideration the intra-block information only. The important properties of partially balanced incomplete block designs with two associate classes with both the intra-block and inter-block analysis in a comprehensive was introduced by the author [2]. Analysis of partially balanced incomplete block designs using simple square and rectangular lattices [3], construction of partially balanced incomplete block designs with two treatments per block and less or equal to ten (10) replications of each treatment [4], the class of partially balanced incomplete block designs with more than two associate classes comprehensively [5], m-classes partially balanced incomplete block designs and association scheme [6] and the table of 2-associate classes of partially balanced incomplete block designs were constructed [7].

There are several patterns and strategies in the construction of PBIBDs for the purpose of estimating higher efficiency factors. Some of the strategies are complicated and complex in designs and the efficiency factors is not satisfactory and comprehensive. The PBIBDs with higher associate classes were constructed by authors in [8, 9, 10]. Triangular and four associate classes PBIBDs with two replicates using dualization method were constructed by other scholars that contributed to the construction of PBIBDs with their different and respective pattern and strategies [11, 12, 13, 14]. In the physical and life experiment, there may be need to prioritize some treatment combinations. It was observed that development of PBIBDs through the incorporation of a constant block-sum strategy [15, 16] and later, in a subsequent work [17] delved deeper into and expanded upon this strategy in the construction of PBIBDs. This design changed the new pattern and strategy towards the construction of PBIBDs by authors cited in [18, 19, 20, 21, 22]. The limitations of the design was that r = 1, λ1 = 1, λ2 = 0. If r = 1, then E = 0 where ‘E’ is the efficiency factor, ‘r’ is the number of replicate and ‘λ’ is the number of treatment pair; this does not fit for analysis. To design an experiment, r 2, for r is positive integer. Instead, the design may be categorized into two or more to fit for analysis [23]. The key concepts for this work are Partially Balanced Incomplete Block Design; Constant Block-Size; Constant Sum-Block; Constant Block-Sum; Efficiency Factors.

Advertisement

2. Design principles and methodology

In this section, the paper defines, give conditions and design CBS-CSB PBIBDs.

2.1 Definition of terms

This sub-section defines Constant Block-Size, Constant Block-Sum and Constant Sum-Block of partially balanced incomplete block designs;

Definition 1: A design is said to have a constant block-size if the number of experimental unit, k, are the same irrespective of the number of blocks, b. For example;

  • t = 21, k = 3, r = 1, b = 7

  • t = 21, k = 3, r = 2, b = 14

  • t = 21, k = 3, r = 3, b = 21

Definition 2: A design is said to be constant block-sum if the number of blocks, b, are the same, b=tk, and the sum of each corresponding number of experimental units are the same, irrespective of the number of treatments and replicates [24]. It is determined by the sum of all the treatments divided by number of blocks, CBS=tb. Furthermore, Khattree model the design as:

t = pq, k = min (p,q), b = max (p,q), r = 1, λ1 = 1, λ2 = 0 where p and q must be odd numbers.

For example;

t=21,k=3,r=1,b=7,λ1=1,λ2=0,n1=2,n2=18E1
CSB=1+2++217=2317=33.E2

Table 1 shows a Constant Block-Sum table, where each row (b1, b2, b3, etc.) represents a block, and the numbers within each block represent the treatments. The key feature of this table is that the sum of values in each block is constant and equal to 33.

Sum
b11112133
b22151633
b33102033
b44121733
b5591933
b66131433
b7781833

Table 1.

Constant block-sum table.

Definition 3: A design is said to be a constant sum-block, if the following conditions exist;

  • There is constant block size.

  • Each corresponding blocks has constant sum.

  • The number of blocks are not the same, b=trk

  • The number of treatments, t, are constant.

  • The number of replicates, r, are at least one, r1

For example;

t=21,k=3,r=2,b=14,λ1=1,λ2=0,n1=4,n2=16.

In Table 2, each row represents a block labeled as b1, b2, b3, and so on. The numbers within each block represent different treatments. For example, in block b1, the treatments are 1, 11, 21, and the sum of these treatments is 33. The same pattern follows for the other blocks. The sum of treatments in each block is the same and equal to 33. This is evident by the “Sum” column at the end of each row, indicating that the values in each row (treatments within a block) add up to 33. This design, with a constant sum within each block is used in experimental design to ensure balance and control for variability. It allows for a systematic distribution of treatments across different blocks while maintaining a consistent sum within each block.

Sum
b11112133
b22151633
b33102033
b44121733
b5591933
b66131433
b7781833
b81122033
b92141733
b10392133
b114111833
b125131533
b13681933
b147101633

Table 2.

Constant sum-block table.

2.2 Conditions for constant block-size with constant sum-block partially balanced incomplete block designs, CBS-CSB PBIBDs

There are conditions for CBS-CSB PBIBDs;

  • The treatments combinations must be prioritized, that is, there must be number of zero associates, λ2 = 0.

  • The maximum number of replicates must be necessarily determine from treatment one.

  • The number of treatments and block sizes remain fixed.

For example;

The designs;

t=21,k=3,λ2=0,

The three numerical values in each row in Table 3 represent different treatments or conditions. For example, in the first row (SN = 1), the values are 1, 21, and 11 and their sum is 33. Similarly, the sum of values in each subsequent row is consistently 33. The key feature here, as in the previous tables, is that the sum of the numerical values in each row is constant and equal to 33. This ensures a balanced distribution of the treatments or conditions. In this context, the table structure represents a specific design where the sum of values in each row is held constant. The associate for treatment 1 is given in Table 4.

SNSum
11211133
21201233
31191333
41181433
51171533

Table 3.

Maximum number of replicates table.

Tλ1 = 1λ2 = 0
111, 21,12,20,13,19,14,18,15,172,3,4,5,6,7,8,9,10,16

Table 4.

Associate table for treatment 1.

Table 4 associates different sets of values with two conditions, denoted as λ1 and λ2, for a specific treatment labeled as T. The row labeled as T corresponds to a specific treatment condition while λ1 and λ2 are conditions or parameters associated with the treatment. For treatment 1, the associated values are presented under λ1 = 1 and λ2 = 0. Under each condition (λ1 and λ2), there are sets of values associated with the treatment. For example, under λ1 = 1, the associated values are 11, 21, 12, 20, 13, 19, 14, 18, 15, 17. Under λ2 = 0, the associated values are 2, 3, 4, 5, 6, 7, 8, 9, 10, 16. Table 4 shows that 1 < r 5 and λ2 = 0. This met the conditions for CBS-CSB PBIBDs.

2.3 Case studies on the design constant block-size with constant sum-block partially balanced incomplete block designs, CBS-CSB PBIBDs

This section shows the construction of CBS-CSB PBIBDs. The researcher and the user will determine the number of replicates based on the priority of the treatment combination and the cost of the experiment (Table 5). This means 1<rmaxr.

Sum
b11112133
b22151633
b33102033
b44121733
b5591933
b66131433
b7781833

Table 5.

Case 1 design.

2.3.1 Construction of design 1

The design t = 21, k = 3, λ2 = 0 will be constructed in five cases since r = 5.

Case 1: t = 21, k = 3, r = 1, b = 7, λ1 = 1, λ2 = 0, n1 = 2, n2 = 18.

The table represents the design with 7 blocks (b1 to b7), and each block contains 3 treatments (k = 3). The values in each block are arranged in a way that the sum of treatments within each block is constant and equal to 33. This design follows a structure for the Constant Block-Sum table provided earlier, with the same sum of treatments in each block.

P1=10018,P2=02215.
E1=0,E2=0,E=0

The efficiency factor (E) is calculated using R software based on a method proposed by authors in [25]. Table 6 associates each treatment (T) with sets of values for the 1st and 2nd conditions. For example, treatment 1 (T = 1) is associated with values 11 and 21 under the 1st condition, and with values 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, 18, 19, 20 under the 2nd condition. The structure of the table indicates how each treatment is linked to specific sets of values under different conditions. This indicates that the efficiency is not directly given but needs to be calculated using specific software and methodology. Case 2 designs are shown in Tables 7 and 8, respectively.

T1st2nd
111, 212,3,4,5,6,7,8,9,10,12,13,14,15,16,17,18,19,20
215,161,3,4,5,6,7,8,9,10,11,12,13,14,17,18,19,20,21
310, 201,2,4,5,6,7,8,9,11,12,13,14,15,16,17,18,19,21
412, 171,2,3,5,6,7,8,9,10,11,13,14,15,16,18,19,20,21
59, 191,2,3,4,6,7,8,10,11,12,13,14,15,16,17,18,20,21
613, 141,2,3,4,5,7,8,9,10,11,12,15,16,17,18,19,20,21
78, 181,2,3,4,5,6,9,10,11,12,13,14,15,16,17,19,20,21
87, 181,2,3,4,5,6,9,10,11,12,13,14,15,16,17,19,20,21
95, 191,2,3,4,6,7,8,10,11,12,13,14,15,16,17,18,20,21
103, 201,2,4,5,6,7,8,9,11,12,13,14,15,16,17,18,19,21
111, 212,3,4,5,6,7,8,9,10,12,13,14,15,16,17,18,19,20
124,171,2,3,5,6,7,8,9,10,11,13,14,15,16,18,19,20,21
136,141,2,3,4,5,7,8,9,10,11,12,15,16,17,18,19,20,21
146,131,2,3,4,5,7,8,9,10,11,12,15,16,17,18,19,20,21
152,161,3,4,5,6,7,8,9,10,11,12,13,14,17,18,19,20,21
162,151,3,4,5,6,7,8,9,10,11,12,13,14,17,18,19,20,21
174,121,2,3,5,6,7,8,9,10,11,13,14,15,16,18,19,20,21
187,81,2,3,4,5,6,9,10,11,12,13,14,15,16,17,19,20,21
195,91,2,3,4,6,7,8,10,11,12,13,14,15,16,17,18,20,21
203,101,2,4,5,6,7,8,9,11,12,13,14,15,16,17,18,19,21
211,112,3,4,5,6,7,8,9,10,12,13,14,15,16,17,18,19,20

Table 6.

Case 1 associates.

Sum
b11112133
b22151633
b33102033
b44121733
b5591933
b66131433
b7781833
b81122033
b92141733
b10392133
b114111833
b125131533
b13681933
b147101633

Table 7.

Case 2 designs.

T1st2nd
111, 21,12,202,3,4,5,6,7,8,9,10,13,14,15,16,17,18,19
215,16,14,171,3,4,5,6,7,8,9,10,11,12,13,18,19,20,21
310, 20,9,211,2,4,5,6,7,8,11,12,13,14,15,16,17,18,19
412, 17,11,181,2,3,5,6,7,8,9,10,13,14,15,16,19,20,21
59, 19,13,151,2,3,4,6,7,8,10,11,12,14,16,17,18,20,21
613, 14,8,191,2,3,4,5,7,9,10,11,12,15,16,17,18,20,21
78, 18,10,161,2,3,4,5,6,9,11,12,13,14,15,17,19,20,21
87, 18,6,191,2,3,4,5,9,10,11,12,13,14,15,16,17,20,21
95, 19,3,211,2,4,6,7,8,10,11,12,13,14,15,16,17,18,20
103, 20,7,161,2,4,5,6,8,9,11,12,13,14,15,17,18,19,21
111, 21,4,182,3,5,6,7,8,9,10,12,13,14,15,16,17,19,20
124,17,1,202,3,5,6,7,8,9,10,11,13,14,15,16,18,19,21
136,14,5,151,2,3,4,7,8,9,10,11,12,16,17,18,19,20,21
146,13,2,171,3,4,5,7,8,9,10,11,12,15,16,18,19,20,21
152,16,5,131,3,4,6,7,8,9,10,11,12,14,17,18,19,20,21
162,15,7,101,3,4,5,6,8,9,11,12,13,14,17,18,19,20,21
174,12,2,141,3,5,6,7,8,9,10,11,13,15,16,18,19,20,21
187,8,4,111,2,3,5,6,9,10,12,13,14,15,16,17,19,20,21
195,9,6,81,2,3,4,7,10,11,12,13,14,15,16,17,18,20,21
203,10,1,122,4,5,6,7,8,9,11,13,14,15,16,17,18,19,21
211,11,3,92,4,5,6,7,8,10,12,13,14,15,16,17,18,19,20

Table 8.

Case 2 associates.

Case 2: t = 21, k = 3, r = 2, b = 14, λ1 = 1, λ2 = 0, n1 = 4, n2 = 16

P1=12214,P2=04411
E1=0.6667,E2=0.4444,E=0.4762

Each row in Table 7 represents a block, and there are 14 blocks in total. The numbers within each block represent different treatments. The sum of treatments in each block is constant and equal to 33, maintaining a similar structure to the design in Case 1. Table 8 associates each treatment (T) with sets of values for the 1st and 2nd conditions. For instance, treatment 1 (T = 1) is associated with values 11, 21, 12, and 20 under the 1st condition, and with values 2, 3, 4, 5, 6, 7, 8, 9, 10, 13, 14, 15, 16, 17, 18, 19 under the 2nd condition.

Case 3: t = 21, k = 3, r = 3, b = 21, λ1 = 1, λ2 = 0, n1 = 6, n2 = 14

P1=23311,P2=2449
E1=0.7111,E2=0.5926,E=0.6238

The table represents the design with 21 blocks (b1 to b21), and each block contains 3 treatments (k = 3). The values in each block are arranged in a way that the sum of treatments within each block is constant and equal to 33. This design follows a structure similar to the Constant Block-Sum tables provided earlier, with the same sum of treatments in each block (Table 9). Table 10 is similar to the previous associates tables. It associates each treatment (T) with sets of values for the 1st and 2nd conditions. For example, treatment 1 (T = 1) is associated with values 11, 21, 12, 20, 13, and 19 under the 1st condition, and with values 2, 3, 4, 5, 6, 7, 8, 9, 10, 14, 15, 16, 17, 18 under the 2nd condition.

Sum
b11112133
b22151633
b33102033
b44121733
b5591933
b66131433
b7781833
b81122033
b92141733
b10392133
b114111833
b125131533
b13681933
b147101633
b151131933
b162112033
b173141633
b18482133
b195101833
b206121533
b21791733

Table 9.

Case 3 design.

T1st2nd
111, 21,12,20,13,192,3,4,5,6,7,8,9,10,14,15,16,17,18
215,16,14,17,11,201,3,4,5,6,7,8,9,10,12,13,18,19,21
310, 20,9,21,14,161,2,4,5,6,7,8,11,12,13,15,17,18,19
412, 17,11,18,8,211,2,3,5,6,7,9,10,13,14,15,16,19,20
59, 19,13,15,10,181,2,3,4,6,7,8,11,12,14,16,17,20,21
613, 14,8,19,12,151,2,3,4,5,7,9,10,11,16,17,18,20,21
78, 18,10,16,9,171,2,3,4,5,6,11,12,13,14,15,19,20,21
87, 18,6,19,4,211,2,3,5,9,10,11,12,13,14,15,16,17,20
95, 19,3,21,7,171,2,4,6,8,10,11,12,13,14,15,16,18,20
103, 20,7,16,5,181,2,4,6,8,9,11,12,13,14,15,17,19,21
111, 21,4,18,2,203,5,6,7,8,9,10,12,13,14,15,16,17,19
124,17,1,20,6,152,3,5,7,8,9,10,11,13,14,16,18,19,21
136,14,5,15,1,192,3,4,7,8,9,10,11,12,16,17,18,20,21
146,13,2,17,3,161,4,5,7,8,9,10,11,12,15,18,19,20,21
152,16,5,13,6,121,3,4,7,8,9,10,11,14,17,18,19,20,21
162,15,7,10,3,141,4,5,6,8,9,11,12,13,17,18,19,20,21
174,12,2,14,7,91,3,5,6,8,10,11,13,15,16,18,19,20,21
187,8,4,11,5,101,2,3,6,9,12,13,14,15,16,17,19,20,21
195,9,6,8,1,132,3,4,7,10,11,12,14,15,16,17,18,20,21
203,10,1,12,2,114,5,6,7,8,9,13,14,15,16,17,18,19,21
211,11,3,9,4,82,5,6,7,10,12,13,14,15,16,17,18,19,20

Table 10.

Case 3 associates.

Case 4: t = 21, k = 3, r = 4, b = 21, λ1 = 1, λ2 = 0, n1 = 8, n2 = 12

P1=3448,P2=5338
E1=0.713,E2=0.6417,E=0.6684

Table 11 represents the design for Case 4, with specific values assigned to each block. Each row corresponds to a block (b1 to b28), and within each block, there are different treatments arranged in a way that the sum of treatments is constant and equal to 33. Table 12 is similar to the previous associate tables. It associates each treatment (T) with sets of values for the 1st and 2nd conditions. For example, treatment 1 (T = 1) is associated with values 11, 21, 12, 20, 13, 19, 14, and 18 under the 1st condition, and with values 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 16, 17 under the 2nd condition. The notations P1 and P2 are partitions with specific characteristics. The E1 = 0.713, E2 = 0.6417, E = 0.6684 values represent efficiency factors. In experimental design, efficiency is a measure of how well an experimental design utilizes the available resources. Here, E1 and E2 are specific efficiency values for P1 and P2, and E is an overall efficiency factor.

Sum
b11112133
b22151633
b33102033
b44121733
b5591933
b66131433
b7781833
b81122033
b92141733
b10392133
b114111833
b125131533
b13681933
b147101633
b151131933
b162112033
b173141633
b18482133
b195101833
b206121533
b21791733
b221141833
b232121933
b243131733
b25492033
b26572133
b276111633
b288101533

Table 11.

Case 4 design.

T1st2nd
111, 21,12,20,13,19,14,182,3,4,5,6,7,8,9,10,15,16,17
215,16,14,17,11,20,12,191,3,4,5,6,7,8,9,10,13,18,21
310, 20,9,21,14,16,13,171,2,4,5,6,7,8,11,13,15,18,19
412, 17,11,18,8,21,9,201,2,3,5,6,7,10,13,14,15,16,19
59, 19,13,15,10,18,7,211,2,3,4,6,8,11,12,14,16,17,20
613, 14,8,19,12,15,11,161,2,3,4,5,7,9,10,17,18,20,21
78, 18,10,16,9,17,5,211,2,3,4,6,11,12,13,14,15,19,20
87, 18,6,19,4,21,10,151,2,3,5,9,11,12,13,14,16,17,20
95, 19,3,21,7,17,4,201,2,6,8,10,11,12,13,14,15,16,18
103, 20,7,16,5,18,8,151,2,4,6,9,11,12,13,14,17,19,21
111, 21,4,18,2,20,6,163,5,7,8,9,10,12,13,14,15,17,19
124,17,1,20,6,15,2,193,5,7,8,9,10,11,13,14,16,18,21
136,14,5,15,1,19,3,172,4,7,8,9,10,11,12,16,18,20,21
146,13,2,17,3,16,1,184,5,7,8,9,10,11,12,15,19,20,21
152,16,5,13,6,12,8,101,3,4,7,9,11,14,17,18,19,20,21
162,15,7,10,3,14,6,111,4,5,8,9,12,13,17,18,19,20,21
174,12,2,14,7,9,3,131,5,6,8,10,11,15,16,18,19,20,21
187,8,4,11,5,10,1,142,3,6,9,12,13,15,16,17,19,20,21
195,9,6,8,1,13,2,123,4,7,10,11,14,15,16,17,18,20,21
203,10,1,12,2,11,4,95,6,7,8,13,14,15,16,17,18,19,21
211,11,3,9,4,8,5,72,6,10,12,13,14,15,16,17,18,19,20

Table 12.

Case 4 associates.

Case 5: t = 21, k = 3, r = 5, b = 21, λ1 = 1, λ2 = 0, n1 = 10, n2 = 10

P1=4555,P2=8227
E1=0.7077,E2=0.6571,E=0.6815

Table 13 presents a Case 5 design. It represents the design for Case 5, with specific values assigned to each block. Each row corresponds to a block (b1 to b35), and within each block, there are different treatments arranged in a way that the sum of treatments is constant and equal to 33. Table 14 contained Case 5 associates. This table shows the associations between the treatments in Case 5. Each row represents a treatment (T1 to T21), and the “1st” and “2nd” columns indicate which other treatments are associated with the corresponding treatment. For example, the first row indicates that for Treatment 1 (T1), the associated treatments in the 1st group are 11, 21, 12, 20, 13, 19, 14, 18, 15, 17, and in the 2nd group are 2, 3, 4, 5, 6, 7, 8, 9, 10, 16. Similarly, each subsequent row provides information about the associations for the corresponding treatment. The efficiency values (E1, E2, and E) are also provided as 0.7077, 0.6571, and 0.6815, respectively (Figure 1).

Sum
b11112133
b22151633
b33102033
b44121733
b5591933
b66131433
b7781833
b81122033
b92141733
b10392133
b114111833
b125131533
b13681933
b147101633
b151131933
b162112033
b173141633
b18482133
b195101833
b206121533
b21791733
b221141833
b232121933
b243131733
b25492033
b26572133
b276111633
b288101533
b291151733
b302102133
b313111933
b324131633
b33582033
b34691833
b357121433

Table 13.

Case 5 design.

T1st2nd
111, 21,12,20,13,19,14,18,15,172,3,4,5,6,7,8,9,10,16
215,16,14,17,11,20,12,19,10,211,3,4,5,6,7,8,9,13,18
310, 20,9,21,14,16,13,17,11,191,2,4,5,6,7,8,13,15,18
412, 17,11,18,8,21,9,20,13,161,2,3,5,6,7,10,14,15,19
59, 19,13,15,10,18,7,21,8,201,2,3,4,6,11,12,14,16,17
613, 14,8,19,12,15,11,16,9,181,2,3,4,5,7,10,17,20,21
78, 18,10,16,9,17,5,21,12,141,2,3,4,6,11,13,15,19,20
87, 18,6,19,4,21,10,15,5,201,2,3,9,11,12,13,14,16,17
95, 19,3,21,7,17,4,20,6,181,2,8,10,11,12,13,14,15,16
103, 20,7,16,5,18,8,15,2,211,4,6,9,11,12,13,14,17,19
111, 21,4,18,2,20,6,16,3,195,7,8,9,10,12,13,14,15,17
124,17,1,20,6,15,2,19,7,143,5,8,9,10,11,13,16,18,21
136,14,5,15,1,19,3,17,4,162,7,8,9,10,11,12,18,20,21
146,13,2,17,3,16,1,18,7,124,5,8,9,10,11,15,19,20,21
152,16,5,13,6,12,8,10,1,173,4,7,9,11,14,18,19,20,21
162,15,7,10,3,14,6,11,4,131,5,8,9,12,17,18,19,20,21
174,12,2,14,7,9,3,13,1,155,6,8,10,11,16,18,19,20,21
187,8,4,11,5,10,1,14,6,92,3,12,13,15,16,17,19,20,21
195,9,6,8,1,13,2,12,3,114,7,10,14,15,16,17,18,20,21
203,10,1,12,2,11,4,9,5,86,7,13,14,15,16,17,18,19,21
211,11,3,9,4,8,5,7,2,106,12,13,14,15,16,17,18,19,20

Table 14.

Case 5 associates.

Figure 1.

Design parameters (t = 21, k = 3).

2.3.2 Construction of design 2

The design t = 15, k = 3, λ2 = 0 will be constructed in five cases since r = 4.

Case 1: t = 15, k = 3, r = 1, b = 10, λ1 = 1, λ2 = 0, n1 = 2, n2 = 12

P1=10012,P2=0229
E1=0,E2=0,E=0

Case 2

t=15,k=3,r=2,b=10,λ1=1,λ2=0,n1=4,n2=10
P1=1228,P2=2227
E1=0.75,E2=0.6,E=0.6364

Case 3

t=15,k=3,r=3,b=15,λ1=1,λ2=0,n1=6,n2=8
P1=2335,P2=4225
E1=0.7302,E2=0.6389,E=0.6751

Case 4

t=15,k=3,r=4,b=20,λ1=1,λ2=0,n1=8,n2=6
P1=4333,P2=7114
E1=0.725,E2=0.6591,E=0.6952

For Case 1 an efficiency of E = 0 implies a 100% efficient design. By achieving E1 = E2 = E = 0, the allocations defined in the P1 and P2 matrices result in a maximally efficient balanced incomplete block design for estimating the main effects of the 3 factors over the 15 runs. Case 2 defines a constant block-size with constant sum-block partially balanced incomplete block (PBIB) design. This achieves variances and overall design efficiency: E1 = 0.75 E2 = 0.6 E = 0.6364. An efficiency of 0.6364 implies the partial blocking and subsampling still allows reasonably efficient estimation of the main effects for the 3 factors at strength 2. The allocations balance runs across subsets to optimize estimation of 2-factor interactions at the cost of some design efficiency compared to a fully balanced design (E = 0).

Case 3 also defines a constant block-size with constant sum-block partially balanced incomplete block (PBIB) design. This allocation achieves variances and efficiency: E1 = 0.7302, E2 = 0.6389 and E = 0.6751. The partial blocking balances runs across subsets to allow estimation of 3-factor interactions. An efficiency of 0.6751 indicates the subsample allocations still allow moderately efficient estimation of main effects, although less precisely than Case 1 or 2 constructions. The matrices P1 and P2 systematically distribute runs to balance representation across subsets at the cost of some variance inflation. In Case 4, the allocation of runs to the two subsamples allows estimation of 3-factor interactions between the variables. The partial blocking distributes runs systematically though unequally across subsets. This achieves reasonable variance for assessing main effects, traded off against balance. Appropriate analysis can still estimate factor impacts. The matrices P1 and P2 define the subdivision rules that split the 15 runs to optimize estimation of 3-way interactions between the 3 factors. The efficiency factor, E, increases as the number of case increases (Figures 2 and 3).

Figure 2.

Design parameters (t = 15, k = 3).

Figure 3.

Design parameters (t = 27, k = 3).

2.3.3 Construction of design 3

The design t = 27, k = 3, λ2 = 0 will be constructed in five cases since r = 7.

Case 1

t=27,k=3,r=1,b=9,λ1=1,λ2=0,n1=2,n2=24
P1=10024,P2=02221
E1=0,E2=0,E3=0

Case 2

t=27,k=3,r=2,b=18,λ1=1,λ2=0,n1=4,n2=22
P1=12220,P2=31120
E1=0.7667,E2=0.6389,E3=0.6557

Case 3

t=27,k=3,r=3,b=27,λ1=1,λ2=0,n1=6,n2=20
P1=14416,P2=51118
E1=0.7284,E2=0.6556,E3=0.671

Case 4

t=27,k=3,r=4,b=36,λ1=1,λ2=0,n1=8,n2=18
P1=16612,P2=71116
E1=0.7115,E2=0.6607,E3=0.6756

Case 5

t=27,k=3,r=5,b=45,λ1=1,λ2=0,n1=10,n2=16
P1=1888,P2=91114
E1=0.702,E2=0.663,E3=0.6774

Case 6

t=27,k=3,r=6,b=54,λ1=1,λ2=0,n1=12,n2=14
P1=110104,P2=111112
E1=0.6958,E2=0.6641,E3=0.6784

Case 7

t=27,k=3,r=7,b=63,λ1=1,λ2=0,n1=14,n2=12
P1=111112,P2=131110
E1=0.6914,E2=0.6648,E3=0.6789
Advertisement

3. Data analysis

3.1 Dataset description

This dataset was collected on the impact of Body Mass Index (BMI) on impact oscillometry (IOS) measures on children with sickle cell diseases. The source of the dataset [26].

For the purpose of this chapter, 45 observations were extracted which has the following factors; asthma, hydroxyurea, ICS and LABA. Furthermore, each treatment was replicated 3 times alongside with the gender. The BMI results were recorded.

Table 15 dataset will be addressed in two cases (2 and 3) based on the priorities. The priority is based on the fixed sum of 24 (Tables 1621).

SNSubject IDHydroxyureaAsthmaICSLABAREPGenderBMI
11YesYesYesYes1Male23.6
21YesYesYesYes1Male23.74
38YesYesYesYes1Male18.39
41YesYesYesNo2Male22.84
51YesYesYesNo2Male23.53
63YesYesYesNo2Male20.95
73YesYesNoNo3Male20.79
821YesYesNoNo3Male15.87
930YesYesNoNo3Male13.98
1013YesNoNoNo4Male16.36
1113YesNoNoNo4Male17.52
1213YesNoNoNo4Male19.07
1325NoYesYesYes5Male17.68
1459NoYesYesYes5male16.71
1564NoYesYesYes5male16.01
1624NoYesYesNo6Male15.36
1725NoYesYesNo6Male16.03
1825NoYesYesNo6Male16.26
1924NoYesNoNo7Male15.25
2037NoYesNoNo7Male15.09
2137NoYesNoNo7Male15.83
2215NoNoNoNo8Male22.03
2336NoNoNoNo8Male14.81
2436NoNoNoNo8Male14.83
2510YesYesYesYes9Female20.19
2610YesYesYesYes9Female20.49
2710YesYesYesYes9Female21.66
2814YesYesNoNo10Female16.02
2914YesYesNoNo10Female16.79
3014YesYesNoNo10Female17.8
3157NoYesYesYes11Female17.68
3257NoYesYesYes11Female16.82
3357NoYesYesYes11Female16.75
342NoNoNoNo12Female27.82
354NoNoNoNo12Female22.29
364NoNoNoNo12Female22.35
3712YesYesYesNo13Female16.41
3812YesYesYesNo13Female17.39
3912YesYesYesNo13Female18.16
4011YesNoNoNo14Female16.13
4111YesNoNoNo14Female17.26
4211YesNoNoNo14Female21
437NoYesYesNo15Female20.98
4416NoYesYesNo15Female20.8
4544NoYesYesNo15Female14.93

Table 15.

Table showing the impact of BMI on IOS.

Sum
b11(23.6)8(22.03)15(20.98)24
b21(23.74)11(17.68)12(27.82)24
b32(22.84)10(16.02)12(22.29)24
b42(23.53)7(15.25)15(20.8)24
b53(20.79)7(15.09)14(16.13)24
b6315.87)8(14.81)13(16.41)24
b74(16.36)9(20.19)11(16.82)24
b84(17.52)6(15.36)14(17.26)24
b95(17.68)6(16.03)13(17.39)24
b105(16.71)9(20.49)10(16.79)24

Table 16.

The dataset arranged in blocks when r = 2.

BlockBBkB2kB2
b11815
23.622.0320.9866.6122.214794436.89
b211112
23.7417.6827.8269.2423.081598.14794.18
b321012
22.8416.0222.2961.1520.381246.43739.32
b42715
23.5315.2520.859.5819.861183.33549.78
b53714
20.7915.0916.1352.0117.34901.682705.04
b63813
15.8714.8116.4147.0915.7739.162217.47
b74911
16.3620.1916.8253.3717.79949.452848.36
b84614
17.5215.3617.2650.1416.71838.012514.02
b95613
17.6816.0317.3951.117.03870.42611.21
b105910
16.7120.4916.7953.9918971.642914.92
Total564.310,77732331.2

Table 17.

Block estimation.

Tb1b2b3b4b5b6b7b8b9b10TQq
123.623.7447.34
22.223.0845.2832.067.66
222.8423.5346.37
20.3819.8640.2436.132.62
320.7915.8736.66
17.3415.733.0333.63−4.6
416.3617.5233.88
17.7916.7134.503−0.62−3.1
517.6816.7134.39
17.031835.03−0.64−2.6
615.3616.0331.39
16.7117.0333.747−2.36−3.9
715.2515.0930.34
19.8617.3437.197−6.86−0.4
822.0314.8136.84
22.215.737.9−1.060.28
920.1920.4940.68
17.791835.7874.89−1.8
1016.0216.7932.81
20.381838.38−5.570.76
1117.6816.8234.5
23.0817.7940.87−6.373.25
1227.8222.2950.11
23.0820.3843.4636.655.84
1316.4117.3933.8
15.717.0332.731.07−4.9
1416.1317.2633.39
17.3416.7134.05−0.66−3.6
1520.9820.841.78
22.219.8642.063−0.284.44

Table 18.

Treatments estimation.

G123Gi
11611
2.06−2.36−6.37−6.67
22712
6.13−6.866.655.92
33813
3.63−1.061.073.64
44914
−0.624.89−0.663.61
551015
−0.64−5.57−0.28−6.49

Table 19.

Group for treatment adjustment.

TQiGiW = krk1QiZ = 1rk1GiT = W-ZTQi
12.06−6.671.545−1.66753.21256.61775
26.135.924.59751.483.117519.11028
33.633.642.72250.911.81256.579375
4−0.623.61−0.4650.9025−1.36750.84785
5−0.64−6.49−0.48−1.62251.1425−0.7312
6−2.36−6.67−1.77−1.6675−0.10250.2419
7−6.865.92−5.1451.48−6.62545.4475
8−1.063.64−0.7950.91−1.7051.8073
94.893.613.66750.90252.76513.52085
10−5.57−6.49−4.1775−1.6225−2.55514.23135
11−6.37−6.67−4.7775−1.6675−3.1119.8107
126.655.924.98751.483.507523.32488
131.073.640.80250.91−0.1075−0.11503
14−0.663.61−0.4950.9025−1.39750.92235
15−0.28−6.49−0.21−1.62251.4125−0.3955
Total151.2204

Table 20.

Treatment adjustment estimates.

(1)G2N = 564.3x564.33010614.483
(2)x210939.08
(3)B2k10,777

Table 21.

Sum of square table.

Case 1

When r = 2

t=15,k=3,r=2,b=10
Letp=rGN=2x564.2830=37.6187

Then q1 = 45.283–37.62 = 7.66

Q1=47.3445.283=2.06
SSB=31
SSTotal=21

The model assumption

Yij=μ+αi+βj+eij

Where α is the treatment, β is the block, eij is the error term and μ is the mean.

The Anova Table 22 shows that both the blocks and the treatments adjustment are significant (p-value ˂0.05). Therefore, we can identify the treatments responsible for the significance (Figure 4).

SVSSdfMSFP-value
blocks162.517918.05749.9770.00557***
Treatments adj151.22041410.80155.9680.0186***
Error10.859661.8099
Total324.59729

Table 22.

ANOVA for the Design when r = 2.

Note: Estimates with “***” means they are significant since the p < 0.05 is considered significant in the study.

Figure 4.

Chart showing the treatment and block means.

The Table 23 shows that treatments 1,2,12 and 15 contributed to the significant of the results (Table 24)

TTiBiMean
147.3445.28323.67
246.3740.24323.185
336.6633.03318.33
433.8834.50316.94
534.3935.0317.195
631.3933.74715.695
730.3437.19715.17
836.8437.918.42
940.6835.78720.34
1032.8138.3816.405
1134.540.8717.25
1250.1143.46325.055
1333.832.7316.9
1433.3934.0516.695
1541.7842.06320.89

Table 23.

Treatment, block and mean.

Sum
b11(23.6)8(22.03)15(20.98)24
b21(23.74)11(17.68)12(27.82)24
b31(18.39)10(16.02)13(16.41)24
b42(22.84)10(16.79)12(22.29)24
b52(23.53)7(15.25)15(20.8)24
b62(20.95)8(14.81)14(16.13)24
b73(20.79)7(15.09)14(17.26)24
b83(15.87)8(14.8313(17.39)24
b93(13.98)9(20.19)12(22.35)24
b104(16.36)9(20.49)11(16.82)24
b114(17.52)6(15.36)14(21)24
b124(19.07)5(16.01)15(14.93)24
b135(17.68)6(16.03)13(18.16)24
b145(16.71)9(21.66)10(17.8)24
b156(16.26)7(15.83)11(16.82)24

Table 24.

Table showing r = 3.

Case 2

When r = 3, then

t=15,k=3,r=3andb=15

Following the similar calculation, the Anova table was estimate in Table 25

SVSSdfMSFP-value
Blocks186.5691413.3262.18960.0674
Treatments adj150.091410.7211.76160.1385
Error97.373166.086
Total434.03244

Table 25.

ANOVA for the design when r = 3.

Table 25 showed that both the blocks and the treatments are not significant since the p-value is greater than 0.05

Advertisement

4. Discussion

This study has introduced Constant Block-Size with Constant Sum-Block Partially Balanced Incomplete Block Designs (CBS-CSB PBIBDs) as a new class of partially balanced incomplete block designs that will be used when there is need to prioritize the treatments combination. This study defined constant block-size; a design is said to have a constant block-size if the number of experimental unit, k, are the same irrespective of the number of blocks, b. Also the study defined constant block-sum; a design is said to be constant block-sum if the number of blocks, b, are the same, b=tk, and the sum of each corresponding number of experimental units are the same, irrespective of the number of treatments and replicates [24]. It is determined by the sum of all the treatments divided by number of blocks, CBS=tb. Furthermore, this study gave conditions on constant sum-block; there is constant block size, each corresponding blocks has constant sum, the number of blocks are not the same, b=trk, the number of treatments, t, are constant and the number of replicates, r, are at least one, r1. This study also gave conditions for Constant Block-Size with Constant Sum-Block Partially Balanced Incomplete Block Designs, CBS-CSB PBIBDs; the treatments combinations must be prioritized, that is, there must be number of zero associates, λ2 = 0, the maximum number of replicates must be necessarily determine from treatment one and the number of treatments and block sizes remain fixed. The study designed t = 21, k = 3; t = 15, k = 3; t = 27, k = 3. And their respective cases were designed and presented with line charts. It was observed when the number of replicate, r, increases the efficiency factors increases. It was also observed that when r = 1, then E = 0.

After the design of the experiments, the study analyzes the experiments using real-life application dataset. The dataset was collected on the impact of Body Mass Index (BMI) on impact oscillometry (IOS) measures on children with sickle cell diseases. [26]. Forty-five (45) observations were extracted which has the following factors; asthma, hydroxyurea, ICS and LABA. Furthermore, each treatment was replicated 3 times alongside with gender. The BMI results were recorded. The dataset was addressed in two cases (2 and 3) based on the priorities. The priority is based on the fixed sum of 24. For the case one, r = 2, the study shows that both the blocks and the treatments adjustment are significant since the p-value is less than 0.05. But the case two, r = 3, the study shows that both the blocks and the treatments are not significant since the p-value is greater than 0.05.

Consequently, the study is not applicable on the even number of treatments, i.e., when t = 2, 4, 6 … Even number of treatments will never give a constant block size [15]. In that case, there is need to research on constant sum-block with even number of treatments.

Advertisement

5. Conclusion

In conclusion, Constant Block-Size with Constant Sum-Block Partially Balanced Incomplete Block Designs (CBS-CSB PBIBDs) has been constructed in this study as a new class of partially balanced incomplete block designs that should be used when it requires prioritizing the treatments combination. The CBS-CSB PBIBDs are designed to provide efficient and precise estimates of treatment effects and they are useful in obtaining maximum information from the available resources, leading to increased experimental efficiency. The CBS-CSB PBIBDs constructed are applicable for testing the impact of various factors on product quality, manufacturing processes, or other parameters of interest. Constant sum-block with even number of treatments will be considered for future research.

References

  1. 1. Bose RC, Nair KR. Partially balanced incomplete block designs. Sankhyā: The Indian Journal of Statistics. 1939;4:337-372
  2. 2. Bose RC, Shimamoto T. Classification and analysis of partially balanced incomplete block designs with two associate classes. Journal of the American Statistical. 1952;47:151-184
  3. 3. Nair KR. Analysis of partially balanced incomplete block designs illustrated on the simple square and rectangular lattices. Biometrics. 1952;8(2):122-155
  4. 4. Clatworthy WH. Partially balanced incomplete block designs with two associate classes and two treatments per block. Journal of Research of the National Bureau of Standards. 1955;54(2):177-190
  5. 5. Hinkelmann K. Extended group divisible partially balanced incomplete block designs. The Annals of Mathematical Statistics. 1964;35(2):681-695
  6. 6. Mesner DM. A new family of partially balanced incomplete block designs with some Latin square design properties. The Annals of Mathematical Statistics. 1967;38(2):571-581
  7. 7. Clatworthy WH, Cameron JM, Speckman JA. Tables of two-associate class partially balanced designs. In: NBS Applied Mathematics Series 63. Washington, DC: National Bureau of Standards; 1973
  8. 8. Garg DK. Pseudo new modified Latin Square (NML m (m)) type PBIB designs. Communications in Statistics—Theory and Methods. 2010;39(19):3485-3491
  9. 9. Kaur P, Garg DK. Construction of some higher associate class PBIB designs using symmetrically repeated differences. Arya Bhatta Journal of Mathematics and Informatics. 2016;8(2):65-78
  10. 10. Sharma K, Garg DK. Construction of three associate PBIB designs using some sets of initial blocks. International Journal of Agricultural and Statistical Sciences. 2017;13(1):55-60
  11. 11. Garg DK, Singh GP. Construction of four associate class PBIB designs using pairing in triplets system. Advances and Applications in Statistics. 2014;41(1):11
  12. 12. Garg DK, Farooq SA. Construction of PBIB designs through chosen lines and triangles of graphs. International Journal of Mathematics Trends and Technology. 2014;8(1):25-32
  13. 13. Saurabh S, Singh MK. Construction of some series of PBIB designs. Journal of Computer and Mathematical Sciences. 2018;9(12):2031-2036
  14. 14. Singh GP, Garg DK. Construction of PBIB designs by using Boolean algebra. Journal of Computer and Mathematical Sciences. 2018;9(2):56-62
  15. 15. Khattree R. A note on the nonexistence of the constant block-sum balanced incomplete block designs. Communications in Statistics-Theory and Methods. 2019;48(20):5165-5168
  16. 16. Khattree R. The Parshvanath yantram yields a constant-sum partially balanced incomplete block design. Communications in Statistics-Theory and Methods. 2019;48(4):841-849
  17. 17. Khattree R. On construction of constant block-sum partially balanced incomplete block designs. Communications in Statistics-Theory and Methods. 2020;49(11):2585-2606
  18. 18. Kaur P, Sharma K. Constant block-sum Youden-m square and PBIB designs using galois field. Communications in Some - Simulation and Computation. 2023;52(11):5396-5400
  19. 19. Bansal N, Garg DK. Construction of some new Youden-m squares, incomplete Sudoku squares and PBIB designs. Communications in Statistics-Simulation and Computation. 2022:1-15. DOI: 10.1080/03610918.2022.2109674
  20. 20. Gupta S. Constant block-sum designs through confounded factorials. Statistics and Applications. 2021;20:93-101
  21. 21. Gupta S. Constant block-sum two-associate class group divisible designs. Statistics and Applications. 2021;19(1):141-148
  22. 22. Khattree R. On construction of equireplicated constant block-sum designs. Communications in Statistics-Theory and Methods. 2022;51(13):4434-4450
  23. 23. Bansal N, Garg DK, Khattree R. Some results on the existence or nonexistence of constant block-sum triangular and other related designs. Journal of Statistical Theory and Practice. 2023;17(1):1
  24. 24. Khattree R. Constant and nearly constant block-sum partially balanced incomplete block designs and magic rectangles. Journal of Statistical Theory and Practice. 2023;17(1):8
  25. 25. Kaur P, Sharma K, Garg DK, Sharma MK. Package ‘PBIBD’. Comprehensive R Archive Network (CRAN): O’Reilly Media, Inc.; 2017
  26. 26. Utkarsh Singh. Impact of BMI on IOS measures on children with sickle cell diseases. 2023. Available from: Kaggle.com

Written By

Babatunde L. Adeleke, Gabriel O. Adebayo and Kazeem A. Osuolale

Submitted: 23 January 2024 Reviewed: 24 January 2024 Published: 10 July 2024