STATISTICAL ANALYSIS OF SPLIT-PLOT DESIGN USING SPECIAL TYPE OF GRAPHS
V. Saranya1, S. Kavitha2*, M. Pachamuthu3 and S. Vijayan4
Research Scholar, Department of Statistics, Periyar University, Salem, Tamil Nadu, India
saranya88stat@gmail.com
2'
'Assistant Professor, Department of Statistics, Periyar University, Salem, Tamil Nadu, India
pustatkavitha@gmail.com 3Assistant Professor, Department of Statistics, Periyar University, Salem, Tamil Nadu, India
pachamuthu@periyaruniversity.ac.in "Research Scholar, Department of Statistics, Periyar University, Salem, Tamil Nadu, India
vijaystatistician1210@gmail.com
Abstract
When all experimental runs cannot be done under homogeneous conditions, blocking can be utilized to increase the power for testing treatment effects. In many real-life environments, there is at least one factor that is hard to change, leading to a split-plot structure. This paper demonstrates how to generate certain graphs using main-plot and sub-plot analyses, as well as providing a catalog. As a result, during situations where the candidate set is too huge to be tractable, the design of split-plot experiments becomes computationally feasible. The designs were considered ideal because they were capable and efficient in estimating the fixed effects of the suitable statistical model given the split-plot design structure. The Split-Plot Design (SPD) is the complete block design which plays an important role in the fields of agriculture, medicine, and industries. This SPD is specifically suited for a two-factor experiment that has more treatments than can be accommodated by a complete block design. In an SPD, one factor is assigned to the main- plot. The assigned first factor is called the main - plot factor. The main-plot is then divided into subplots and the second factor is called the sub - plot factor. SPD is most used for (i) few experimental materials may be rare while the other experimental materials may be available in large quantity, (ii) the levels of one or more treatment factor or easy to change and the alteration of levels of other treatment factors are costly or time-consuming. Given the extensive study done in graph theory, it has developed to be a very broad subject in mathematics. Graphs are important because they are a visual way of expressing information. A graph shows data that is equivalent to many words. A graph can convey information that is difficult to express in words. A bipartite graph is a type of graph in which the entire graph may be divided into two bipartite sets, with edges connecting vertices in one set to vertices in the other. Vertex coloring is the procedure of assigning labels or colors to each vertex in a graph. The data set was also manually analyzed to validate the software-analyzed outcomes. R gave the same results as the manual analysis, showing that they were both correct. R is mainly command-based. The proposed approach is demonstrated using agricultural and industrial examples.
Key Word: Split-plot design, complete bipartite graph, colored graph.
1. Introduction
Tamil Nadu is one of the leading rice-growing states in India and has been successfully cultivating rice since ancient times as the state has all the favorable climatic conditions suitable for rice cultivation. Rice research was initiated in Madras State to increase rice production and productivity (mixed).1902 at Samalkota in East Godavari district was extended to 12 more places. This study uses primary data to determine the rice production of the Salem district of Tamil Nadu. The present study aims to analyze rice cultivation using different levels of nitrogen. The constraint analysis is applied to find out the problems of paddy cultivation.
A split plot design and some graphs were adopted for this present study. The split-plot design originated in the field of agriculture. Experimenters applied one treatment to a large area of land, called a whole plot, and other treatments to smaller areas of land within the whole plot called a subplot. Split plots have two types of factors Hard-to-change (HTC) applied to the whole plots and Easy-to-change (ETC) applied to the subplots. In such case different sizes of plots are required and the resulting design is known as Split Plot Design (SPD). In 1925, Fisher developed these designs for the purpose of agricultural experiments.
One of the fastest-expanding sciences in modern technology is graph theory. Graphs are commonly used in applications of many fields to represent different objects and their relationships. The declaration of an object vertex serves as the graph's visual representation, while an edge represents the relationships between objects. Graph theory has recently become established as a significant mathematical tool in a wide range of fields, including functional research, chemistry, genetics, and linguistics, as well as electrical engineering, geography, sociology, and architecture of themselves.
Wooding W M [1] has discussed split-plot designs characteristics and applications. To design the first section, models and least squares are reviewed. The main part shows how a fundamental split-plot design is created through a process of "evolution," starting with a completely random model and progressing through a randomized blocks design to a split-plot while using the same set of runs. George Box and Stephen Jones [2] have evaluated the applicability of split-plot designs for the experimental setting and have concentrated on the use of statistical experimental designs in designing goods that are robust to environmental factors. They conclude that the split-plot and strip-block designs are valuable for creating strong products. Peter Goos and Martina Vandebroek [3] have developed an exchange algorithm for constructing D-optimal split-plot designs and the resulting designs are analyzed. Natalino Calegario et al. [4] have analyzed the split-split-plot design and established the impact of fertilizer concentration on the establishment of Begonia and Petunia. Then they draw the conclusion that the pH values declined with fertilizer concentration over time and the EC values increased over time, resulting in values that limited nutrient availability and plant growth.
Bradley Jones and Peter Goos [5] have suggested a fresh technique for producing ideal split-plot designs. These split-plot designs are best when they are effective at estimating the fixed effects of the proper statistical model, given the structure of the design. Pwasong A D and Choji D N [6] have analyzed the rabbit feeds data obtained from the Department of Agricultural Science, Federal College of Education Pankshin and determined that there is any significant variation in the categories of feeds given. The result illustrates that there was no significant different between the various types of feed utilized to feed the rabbits.
Bradley Jones [7] has suggested for the use of split-plot designs in industrial applications are provided after an examination of current developments. Johannes Ledolter [8] has reviewed the factorial split-plot design and fractional factorial split-plot designs experiments and uses several illustrative examples to illustrate why they frequently occur in industrial investigations. Abhishek K. Shrivastava [9] has presented an effective method for constructing split-plot design catalogues by transforming the design isomorphism problem to a graph isomorphism problem utilizing a new graph
representation. Derya Dogan and Pinar Dundar [10] have introduced the new concept of average covering number of a graph and establish the brief relationship between the average covering number and some other graph parameters.
David I J and Adehi M U [11] have utilized a 21x52 split-plot experiment with three replicates for comparison. Here they review the sorghum thresher's improved threshing efficiency.Vahide Hajihassani and Yadollah Rajaei [12] have used five agricultural machinery companies that have existed accepted in the Tehran Stock Exchange since (1388-1390) and a sample that is representative of society, conduct a split-plot design model study on the factors impacting liquidity accepted in stock exchange Agricultural Machinery companies. David I J et al. [13] have presented the steps for the estimated generalized least square (EGLS) technique, which estimates the parameters of a nonlinear split-plot design (SPD) model utilizing theoretical iterative Gauss Newton via Taylor Series expansion. Yoshimi Egawaa et al. [14] have discussed the 4-connected graph in triangles and let G be a 4-connected graph, and let E~(G) denote the set of those edges of G which are not contained in a triangle, and let Ec(G) denote the set of 4-contractible edges of G. We show that if 3 < I E~ (G)l < 4 or I E~ (G) I > 7, then I Ec(G) I > (I E~(G) I +8)/4 unless G has one of the three specified configurations.
Table 1: Background of this research
The related articles of Split-plot design, application of split-plot design, graph theory, vertex coloring and split-plot design with colored graph are given.
Yield of paddy in different level of nitrogen and the given data are
collected from the agriculture filed of salem district._
Example 2 Application method of paint in different mixing level and the given data
are collected from different hardware's in salem district.
Review
Example 1
2. Preliminaries
2.1 Split-plot Design
A randomized complete block design with two factors is no longer a randomized complete block design because the order of experiments is controlled to obtain observations in each treatment under each block. Splitting the randomization of an experiment to obtain observations under the treatment of one factor is called a split-plot design.
2.2 Complete Bipartite Graph
A complete bipartite graph is a graph whose vertex set V can be divided into two subsets V1 and V2 such that no edge has both endpoints in the same subset and every edge is connected to every vertex of the first subset and every vertex of the second subset.
2.3 Colored Graph
In a graph, the procedure for assigning the labels (colors) to the nodes or edges or areas is known as graph coloring. In this assignment no two adjacent vertices or adjacent edges or adjacent areas are getting the same color.
3. Statistical Analysis of Split Plot Design
The liner model for Split Plot Design is.
Yijk = ^ + n + tj + sk + tsjk + du + eijk, V i=1,2,...r; j=1,2,...v; k=1,2...n.
Where, Ytjk is the observation corresponding to kth level of sub plot factor (B), jth level
of main plot factor (A) and ith replication.
I is general mean effect.
rt is ith replication effect.
tj is jth main - plot treatment effect.
sk is kth sub - plot treatment effect.
tSjk is interaction effect.
The error components dtj and £ijk are independently and normally distributed with mean zero and respective variance a2 d andc2e.
3.1 Main-Plot Analysis
This analysis part is based on the comparisons of main plot totals:
The levels of A are assigned to the main plots within blocks based on RBD and the sum of squares are given below,
G2
• Correction factor (CF) = —
rvn
• Total sum of square (SST) = EX2 — CF
%r2
• Replication sum of square (SSR) =--CF
EA2
Main-plot sum of square (SS (MP)) =--CF
1 1 rn
• Main-plot error sum of square (SSEi) :
3.2 Sub-Plot Analysis
E(AR)2
— CF—SSR—SS(MP)
This analysis part is based on the comparisons of sub plot totals:
Sub-plot sum of square (SS (SP))
IB2
— CF
Interaction effect sum of square (AxB) =
E(AB)2
■ - CF — SSA — SSB
• Sub-plot error sum of square (SSE2) = SST- (SSR + SS (MP) + SSEi + SS (SP) + (AxB)) The analysis of the variance table is outlined as follows
Table 2: ANOVA for split-plot designs
Sv Df Ss Mss F-Ratio
Replication (r-1) S2 SR SSR FR - SR
SSR (r — 1) SE1 F(r-1),(r-1)(v—1)
MP(A) (v-1) S2 -SA - SSA fa = SA
SSA (v — 1) SE1 F(v-1),(r-1)(v—1)
MPE (Ei) (r-1)(v-1) SSEi 2 SE - SSE±
(r — 1)(v — 1) -
SP (B) (n-1) SSB 2 SB SSB FB = SB
(n — 1) SE2 F(v-1),(r-1)(n—1)
IE(AB) (v-1)(n-1) SSAB 2 SAB SSAB fab c2 _ SAB
(v — 1)(n — 1) SE2~F(v-l),(r-1)(n-1)
n
rv
r
SPE(E2) v(r-1)(n-1) SSE2 „2 = SSE2 -
_A v(r - 1)(n - 1)_
Total rvn-1 SST - -
3.3 Flow - Chart
Figure 1: Flow chart
4. Construction of Split-Plot Design using Complete Bipartite Graph
4.1 Method for Construction of Complete Bipartite Graphs
Let us consider the main-plot and sub-plot as vertex set S . This vertex set can be divided into subsets of S1and S2.
In main-plot, the replication is considered as first subset S1and variety as second subset S2.
Now consider the first vertex of first subset and then R1 is connected to all the vertices of second subset through edges.
Next consider the second vertex and it is connected to all the vertices of the second subset through the edges.
Similarly, all the remaining vertices of the first subset are connected to all the vertices
of second subset through the corresponding edges.
Finally, we get the complete bipartite graph for main-plot and sub-plot.
4.1.1 Application
This example is to determine the yield response in N fertilization between different paddy varieties, three varieties of Paddy (Vi = ADT 36, V2=ASD 16, V3 = IR50) are the treatments of main plot, nitrogen rates such as 0, 30 and 60 Kg/ha are the sub- plot treatments. The study was replicated four times and the primary data gathered for this experiment from the agricultural field of salem district of tamil nadu in India and shown in table 3.
Table 3: Replication wise data for yield of paddy (Kg/ha)
Replication Ri R2 R3 R4
Variety Nitrogen (N1)
Vi 15.8 19.2 13.2 13.2
V2 20.8 15.3 20.5 13.8
V3 15.9 16.3 16.2 12.8
Nitrogen (N2)
Vi 17.8 20.5 14.8 13.8
V2 24.8 20.8 18.8 17.8
V3 18.5 16.1 20.8 12.2
Nitrogen (N3)
Vi 21.1 24.8 13.8 18.8
V2 30.5 19.2 25.7 15.2
V3 18.3 18.2 22.8 10.8
Table 4: Replication x variety (R x V) for main - - plot
Vi V2 V3 Replication Total
Ri 54.7 76.1 52.7 183.5
R2 64.5 55.3 51 170.8
R3 41.8 65 59.8 166.6
R4 45.8 46.8 35.8 128.4
Variety Total 206.8 243.2 199.3 649.3
The procedure for constructing the complete bipartite graph mentioned in section 4.1 is followed for the main-plot and sub-plot for the above experiments and then the finalized complete bipartite graph.
From the above table 4 as vertex is fixed as S, which is divided into two subsets, figure 2 shows that S1 (replication) and S2 (variety). Figure 3 shows that the first replication vertex (Ri) and it is connected to all the vertices of variety (V1, V2 and V3) through the edge values 54.7(Yi), 76.1(Y2) and 54.7(Y3).
Figure 2: Graph of subsets Figure 3: Graph for first replication (Ri)
Figure 4 shows that the second replication vertex (R2) and it is connected to all the vertices of variety (V1,V2 and V3) through the edge values 64.5(V1) 55.3(y2) and 51(Y3). Similarly, figure 5 shows that the third and fourth replication vertices R3 and R4 are connected to all the vertices of variety (V1, V2and V3) through the corresponding edge values (Yv Y2and Y3) 41.8,65, and 59.8 (Yv Y2and Y3) 45.8, 46.8 and 35.8.
Figure 4: Graph for second replication (R2) Figure 5: Graph for third and fourth replication (R3 and R4)
Finally, figure 6 shows that the complete bipartite graph of variety and replication for main - plot.
Replication Variety
Ms^iJV/
1^55.3 \/ j
sK /
A ^•Vl
Rj< li fiS V \
/
Ri^T--5W- —^V'3
Fig. 6: Graph for complete bipartite graph of main - plot Table 5: Variety x nitrogen (V x N) for sub-plot
Ni N2 N3 Variety Total
Vi 61.4 66.9 78.5 206.8
V2 70.4 82.2 90.6 243.2
V3 61.2 67.6 70.5 199.2
Nitrogen Total 193 216.7 239.6 649.3
The construction of complete bipartite graph for the sub-plot are given below.
From the above table 5 as vertex is fixed as G, which is divided into two subsets, figure 7 shows that G1(variety) and G2 (nitrogen). Figure 8 shows that the first variety (VI) is connected to all nitrogen (N1t N2and N3) through the values 61.401), 66.901) and 78.501).
Figure 7: Graph for vertex subset Figure 8: Graph for first variety (Vi)
Similarly, figure shows that the second and third variety V2 and V3 is connected to all the nitrogen (W1, N2 and N3) through the corresponding values 70.4, 82.2 and 90.6 (Y1,Y2 and Y3), 61.2, 67.6 and 70.7 (Y1,Y2 and Y3).Finally, figure 10 shows the complete bipartite graph for variety and nitrogen.
The sum of squares for main- plot:
• Correction factor (CF) = 17910.8469
• Total sum of square (SST) = 637.46
• Replication sum of square (SSR) = 187.709
• Variety sum of square (SSV) = 91.9006
• Main-plot error sum of square (SSE1) = 175.4527 The sum of squares for sub-plot:
• Nitrogen sum of square (SSN) = 90.4906
• Interaction effect sum of square (V*N) = 10.4194
• Sub-plot error sum of square (SSE2) = 81.4869 The ANOVA table for split-plot design is shown in below table:
Table 6: ANOVA table for split-plot design
Sv Df Ss Mss F-Ratio P-Value
Replication 3 187.7098 62.570 2.1397 0.196473
Variety (A) 2 91.9006 45.950 1.5714 0.282634
Main - plot error(Ei) 6 175.4527 29.242 - -
Nitrogen(B) 2 90.4906 45.245 9.9945 0.001204**
Interaction (AB) 4 10.4194 2.605 0.5754 0.684090
Sub - plot error(E2) 18 81.49 4.527 - -
Total 35 - - - -
The table value of replication and variety are greater than the calculated values. So, the null hypothesis is accepted. There is no significant difference between the four replications and three varieties. The table value of nitrogen level is greater than the calculated value. So, the null hypothesis is accepted. There is no significant difference between the three nitrogen levels. The table value of the interaction effect is
V. Saranya, S. Kavitha, M. Pachamuthu, S. Vijayan RT&A, No 1 (77)
SPLIT-PLOT MODEL USING SPECIAL TYPE OF GRAPHS_Volume 19, March 2024
also greater than the calculated value. So, the null hypothesis is accepted.
There is no significant difference between the interaction effects. The P-value of the above experiment is greater than the 5% level of significant. Therefore, the null hypothesis is accepted. There is no significant difference that occurred in the above experiment.
4.2 Method for Construction of Colored Graph
• Let us consider the main-plot and sub-plot factors as the vertex set 51 and S2 . Here the common factor assigned as 51 and the other factor assigned S2.
• The vertices of set 51 and S2 are colored using the vertex coloring and the vertices are differentiate with different colors. Now consider the first vertex of 51 and it is connected to the corresponding vertices of S2 through edges.
• Next the second vertex of 51 which is connected to the corresponding vertices of S2 through edges.
• Similarly, all the remaining vertices of the first set are connected to the corresponding vertices of second set through the corresponding edges.
• Finally, we get the colored graph (vertex coloring graph) for main-plot and sub-plot.
4.2.1 Application
The test is designed to examine pigment dispersion in paint. Three different mixing levels of a particular pigment are studied. The procedure consists of three application methods (brushing, sparing, and rolling) and measured the percentage reflectance of a pigment. Four days required running the experiment from hardware shops in salem district and the data obtained below.
Table 7: Replication wise data form mixes level and application method ofpaint
Replication Ri R2 Rs R4
Application Method Mixing level (Mi)
Ai 65.8 70.2 65.2 69.2
A2 69.8 65.3 70.5 63.8
As 70.8 67.3 68.2 69.8
Mixing level (M2)
Ai 68.7 73.5 69.9 66.8
Ä2 74.8 70.8 68.8 67.8
As 50.8 69.1 71.8 63.2
Mixing level (M3)
Ai 72.2 77.8 71.6 70.8
A2 81.5 69.2 75.7 65.2
As 69.3 71.6 77.8 60.8
Table 8: Replication x application method (RxA) for main - plot
Ai A2 As Replication Total
Ri 206.8 226.1 190.9 623.8
R2 221.5 205.3 208 634.8
R3 206.7 215 217.8 639.5
R4 206.8 196.8 193.8 597.4
Application Total 841.8 843.2 810.5 2495.5
Table 9: Application method x mixing (A xM) for sub-plot
M, M2 M3 Application Total
A, 270.4 219 292.4 841.8
A2 269.4 282.2 291.6 843.2
A3 276.1 254.9 279.5 810.5
Mixing Total 815.9 816.1 863.5 2495.5
The procedure for constructing the colored graph mentioned in section 4.2 for main-plot and sub-plot for the above experiments and then the finalize.
Here we take main-plot factors (replication and application) and sub-plot factors (application method and mixing level) as the set 5X and S2. Figure 11 shows that the Here 5X consists of the common factor which is application and S2 consists of factors such as replication and mixing level. Next, figure 12 shows that the first vertex A1 of first set and then A1 is connected to the corresponding (replication and mixing) vertices of the second set through edges.
Figure 11: Colored graph of vertex set Figure 12: Colored graph for first application method
Next figure 13 shows that the second vertex of first set and it is connected to the corresponding vertices of second set. Similarly, figure 14 shows that the third vertex of first set are connected to the corresponding vertices of second set and finally, we get the colored graph (vertex coloring graph) for main-plot and sub-plot.
V. Saranya, S. Kavitha, M. Pachamuthu, S. Vijayan RT&A, No 1 (77)
SPLIT-PLOT MODEL USING SPECIAL TYPE OF GRAPHS Volume 19, March 2024
Figure 13: Colored graph for second application method Figure 14: Colored graph for main - plot and sub - plot
The sum of squares for main-plot:
• Correction factor (CF) = 172986.6738
• Total sum of square (SST) = 975.3762
• Replication sum of square (SSR) = 118.269533
• Application method sum of square (SSA) = 56.970367
• Main-plot error sum of square (SSE1) = 173161.9137 The sum of squares for main-plot:
• Sub-plot sum of square (SSM) = 125.348
• Interaction effect sum of square (A*M) = 87.0449
• Sub-plot error sum of square (SSE2) = 334.2397
Table 10: ANOVA table for split- plot design
Sv Df Ss Mss F-Ratio P-Value
Replication 3 118.2695 39.4217 0.9330 0.48062
Application method(A) 2 56.97037 28.4852 0.6742 0.54435
Main Plot Error(Ei) 6 253.503 42.2505 - -
Mixing(M) 2 125.3487 62.674 3.3752 0.05691
Interaction (AM) 4 87.0449 21.7612 1.1719 0.35616
Sub-plot Error(E2) 18 334.2397 18.5688 - -
Total 35 - - - -
The table values of replication and application method are greater than the calculated values. So, the null hypothesis is accepted. There is no significant difference between the four replications and three application methods. The table value of mixing level is greater than the calculated value. So, the null hypothesis is accepted. There is no significant difference between the three application methods. The table value of the interaction effect is also greater than the calculated value. So, the null hypothesis is accepted. There is no significant difference between the interaction effects.
The P-value of the above experiment is greater than the 5% level of significant. Therefore, the null hypothesis is accepted. There is no significant difference that occurred in the above experiment.
V. Saranya, S. Kavitha, M. Pachamuthu, S. Vijayan
SPLIT-PLOT MODEL USING SPECIAL TYPE OF GRAPHS_
Table 11: Table for comparative study
ANOVA Example 1 Example 2
Null hypothesis is Null hypothesis is
Traditional method
accepted accepted
R-Software method Null hypothesis is accepted Null hypothesis is accepted
5. Conclusion
Many of these real-world agricultural and industrial experiments involve factors called HTC. In these situations, experimenters have realized that the most efficient way to conduct an experiment is to fix the level of the hard-to-change factor and then run all or some combination of the easily changeable factors. This is repeated a few times. As we have seen, this leads to a split-plot design. Accounting for the split-plot nature of the design is equally important in the analysis of the data because the split-plot test contains two error terms. The present paper is classified into three parts namely rice production, nitrogen level and variety of rice in salem district. Constraint analysis is used to increase rice production. To construct and analyze the SPD using some special type of graphs through numerical esxamples from different field and the hypothesis testing is compared by the split-plot ANOVA method with software using method. When comparing the results of these methods, they produce the same results. Here some special type of graphs is used to construct the SPD. In future, there is an idea to expanding this procedure to other experimental designs such as strip-plot design and incomplete block designs etc.
List of Abbreviation
Sv - Sources of variance
Df - Degrees of freedom
Ss - Sum of squares
Mss - Mean sum of squares
MP(A) - Main -Plot(A)
MPE(E1) - Main - Plot Error (E1)
SP(B) - Sup -Plot(B)
SPE(E2) - Sup - Plot Error (E2)
IE(AB) - Interaction Effect (AB)
References
[1] George Box and Stephen Jones. (1992). Split-Plot Designs for Robust Product Experimentation, Journal of Applied Statistics, (1):3-26.
[2] Robin J. Wilson. Introduction to Graph Theory, Fourth Edition, Addison Wesley Longman Limited, England, 1996.
[3] Peter Goos, and Martina Vandebroek. (2001). Optimal Split-Plot Designs, Journal of Quality Technology, 33(4): 436-450.
[4] Calegario, C. L. L., Calegario, N., Lemos, P. C., and de Mello, A. A. (2006). Effect of different fertilization practices on the growth of two plant species: Begonia sp. AND Petunia sp. Cerne, i2(2):107-112.
[5] Bradley Jones, and Peter Goos. (2007). A Candidate-Set-Free Algorithm for Generating D-Optimal Split-Plot Designs, Applied Statistics, 56(3):347-364.
[6] Pwasong, A D and Choji D N. (2008). The Application of the Split-Plot Design in the Analysis of Experimental Rabbit Feeds Data, Science World Journal, 3(2):73-77.
[7] Bkggradley Jones. (2009). Split-Plot Designs: What, Why, and How, Journal of Quality Technology, 41(4): 340-361.
[8] Johannes Ledolter. (2010). Split-Plot Designs: Discussion and Examples", International Journal Quality Engineering and Technology, 1(4):441-457.
[9] Abhishek, K. Shrivastava. (2013). Efficient Construction of Split-Plot Design Catalogs using Graphs, IIE Transactions,45, 1137-1152.
[10] Derya Dogan, and Pinar Dundar. (2013). The Average Covering Number of a Graph Hindawi Publishing Corporation, Journal of Applied Mathematics, 1-4.
[11] A. Bondy, U.S.R. Murty. Graph theory, Springer International Edition, 20i3.
[12] David, I J and Adehi, M U. (2014). Effectiveness of Split-Plot Design over Randomized Complete Block Design in Some Experiments", Journal of Biology, Agriculture and Healthcare, 4 (19): 75-80.
[13] Vahide Hajihassani, and Yadollah Rajaei. (2015). Using Model of Split-Plot Design in the Study of Factors Affecting on Accepted in Stock Exchange Agricultural Machinery Companies Liquidity, Indian Journal of Science and Technology, 8:526-529.
[14] David I J, Asiribo, O E and Dikko, H G. (2018). On Parameters Estimation of Nonlinear Split-Plot Design Model with EGLS-MLE, Professional Statisticians Society of Nigeria, 2:381-386.
[15] Douglas, C. Montgomery. Design and analysis of experiments, Fifth Edition, International Student Version, Arizona State University, 2016.
[16] Egawa, Y., Kotani, K., and Nakamura, S. (2018). Structure of edges in a 4-connected graph not contained in triangles and the number of contractible edges. AKCE International Journal of Graphs and Combinatorics, i5(2):202-210.
[17] David, I. J., Asiribo, O. E., and Dikko, H. G. (2019). Parameter estimation of nonlinear split-plot design models: a theoretical framework. Journal of Reliability and Statistical Studies, 117-129.
[18] Santoso, K. A., Agustin, I. H., Prihandini, R. M. and Alfarisi, R. (2019, April). Vertex colouring using the adjacency matrix. In Journal of Physics: Conference Series 1211(1): 012019.