In the previous discussion we discussed the effect of a single treatment on certain responses . The single treatment is called the factor , and the level or level of the factor is called the level . The factor is symbolized by a capital letter while the level of the factor is symbolized by a lowercase letter . If we simultaneously observe the effect of several factors in the same study, the experiment is called a factorial experiment .
A factorial experiment is an experiment in which the treatment consists of all possible combinations of levels of several factors. Experiments using f factors with t levels for each factor are symbolized by the factorial experiment f t . For example, a 2 2 factorial experiment means that we use 2 factors and the level of each factor consists of 2 levels. Factorial experiment 2 2It is also often written in the form of a 2x2 factorial experiment. The last symbol is often used for factorial experiments where the level of each factor is different, for example 2 levels for factor A and 3 levels for factor B, the experiment is called a 2x3 factorial experiment. A 2x2x3 factorial experiment means a factorial experiment consisting of 3 factors with levels for each factor of 2, 2, and 3. described previously, and the next stage is the selection of the environmental design, which involves the form of experimental designs such as:
- Completely Randomized Design (CRD)
- Randomized Complete Block Design (RCBD)
- Latin Square Design (LSD)
- Split Plot Design
- Split-Split Plot Design
- Strip Plot Design
Sub-discussion:
- Introduction
- Factorial Experiment in Completely Randomized Design
- Randomization and Factorial Trial Plans in a Completely Randomized Design
- Factorial Design Linear Model in RAL
- Example of CRD Factorial
The full discussion can be read in the document below.
Introduction
In the previous discussion we have discussed the effect of single treatment on certain responses. The single treatment is called a factor, and the level or level of the factor is called a level. Factors are symbolized by capital letters while the degree of factors is symbolized by lowercase. If we simultaneously observe the effect of several factors in the same study, then the experiment is called a factorial experiment.
A factorial experiment is an experiment whose treatment consists of all possible combinations of several factors. Experiments using f factors with level for each factor are symbolized by factorial experiments ft. For example, the factorial experiment 22 means that we use 2 factors, and the level of each factor consists of 2 levels. Factorial experiments 22 are also often written in the form of 2x2 factorial experiments. The latter is often used for factorial experiments where the levels of each factor are different, for example 2 levels for factor A and 3 levels for factor B, then the experiment is called a factorial experiment 2x3. A 2x2x3 factorial experiment means a factorial experiment consisting of 3 factors with a level for each of its factors 2, 2, and 3 respectively. Thus, in the factorial experiment, there are two stages that need to be carried out, firstly, namely the treatment design, as previously described, and then the stage of selecting the environmental design, which concerns experimental design forms such as CRD, RCBD, LATIN, Split-plot, Split-Block.
The purpose of the factorial experiment is to look at the interactions between the factors we are trying. Sometimes the two factors synergize with each other to the (positive) response, but sometimes also the existence of one factor hinders the performance of other (negative) factors. The existence of these two mechanisms tends to increase the effect of interactions between the two factors. Interaction measures the failure of the effect of one factor to remain the same at every level of another factor or simply, the interaction between factors is whether the effect of a particular factor depends on the degree of another factor? For example, if the simple effect of N is the same at each level of fertilizer application P, the two factors are mutually independent and it is said that there is no interaction, while if the application of N has a different effect on each level of P, it is said that there is an interaction between Factor N and Factor P.
For example, if we want to observe the effect of nitrogen (N) administration consisting of 2 levels (n0, and n1) and phosphorus administration (P) consisting of 2 levels (p0, p1) on certain responses, with the following results:
Table 19. Simple effect, main effect, and Interaction effect
Factor | Nitrogen (N) | Average P | Simple effect n | |
Phosphorus (P) | n0 | n1 |
| n1-n0 |
p0 | 40 | 48 | 44 | 8 (se N, p0) |
p1 | 42 | 51 | 46.5 | 9 (se N, p1) |
Average N | 41 | 49.5 | 45.25 | 8.5 (me N) |
Simple effect of P (p1-p0) | 2 (se P, n0) | 3 (se P, n1) | 2.5 (me P) |
|
The difference between n1 – n0 and p1 – p0 is called simple effects symbolized by (se N) and (se P). The average of simple effect is called the main effect, symbolized (me N) and (me P).
The approximate effect of the interaction and the main effect of the average treatment can be calculated by the following formula:
Simple effect (se):
$\begin{matrix}se\ N\ \text{ on }\ p_0=n_1p_0-n_0p_0\\=48-40\\=8\\se\ N\ \text{ on } \ p_1=n_1p_1-n_0p_1\\=51-42\\=9\\\end{matrix}$ $\begin{matrix}se\ P\ \text{ on } \ n_0=p_1n_0-p_0n_0\\=42-40\\=2\\se\ P\ \text{ on }\ n_1=p_1n_1-p_0n_1\\=51-48\\=3\\\end{matrix}$
Main Effect (main effect, me):
$\begin{matrix}me\ N=\frac{1}{2}(se\ \ N\ \ on\ \ p0+se\ \ N\ \ on\ \ p1)\\=\frac{1}{2}\left[(n_1p_0-n_0p_0)+(n_1p_1-n_0p_1)\right]\\=\frac{1}{2}\left[(48-40)+(51-42)\right]\\=8.5\\\end{matrix}$
$\begin{matrix}me\ P=\frac{1}{2}(se\ \ P\ \ on\ \ n0+se\ \ P\ \ on\ \ n1)\\=\frac{1}{2}\left[(p_1n_0-p_0n_0)+(p_1n_1-p_0n_1)\right]\\=\frac{1}{2}\left[(42-40)+(51-48)\right]\\=2.5\\\end{matrix}$
Interaction Effects:
$\begin{matrix}interaction\ N\times P=\frac{1}{2}\left[(n_1p_0-n_0p_0)-(n_1p_1-n_0p_1)\right]\\=\frac{1}{2}\left[(48-40)-(51-42)\right]\\=-0.5\\\end{matrix}$
or
$\begin{matrix}=\frac{1}{2}\left[\left(p_1n_0-p_0n_0\right)-\left(p_1n_1-p_0n_1\right)\right]\\=\frac{1}{2}\left[\left(42-40\right)-\left(51-48\right)\right]\\=-0.5\\\end{matrix}$
This simple effect is required by the user (farmer, for example), if he only limits to the use of a certain degree of one of the factors. For example, if farmers want to see the difference in the effect of N on each fertilization level of P, the simple effect of N at the level of p0 = 8 and at the level of p1 = 9.
The difference between simple effect and interactions graphically can be visualized as follows:
Figure 8. The difference between simple effect and interaction
The possibilities that can occur between the main effect and the interaction are presented in the following Figure:
Sources of Variance | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
A | ns | * | ns | * | ns | * | ns | * |
B | ns | ns | * | * | ns | ns | * | * |
AxB | ns | ns | ns | ns | * | * | * | * |
Figure 9. Various possible forms of interaction between factors
Advantage:
- More efficient in using existing sources
- The information obtained is more comprehensive because we can study the main effect and interactions
- The results of the experiment can be applied in a broader range of conditions because we study a combination of various factors
Disadvantage:
- Statistical Analysis becomes more complex
- There are difficulties in providing relatively homogeneous experimental units
- the effect of a certain combination of treatments may not mean anything so that there is a waste of existing resources
Factorial Experiments in a Completely Randomized Design
Factorial experiments in a Completely Randomized Design are factorial experiments using a Completely Randomized Design as their environmental design. In principle it is the same as a Completely Randomized Design, but in this case the factors tried are more than one.
Randomization and Factorial Trial Plan in a Completely Randomized Design
The randomization method is the same as a complete random design. The placement of treatments that are a combination of the level of factors to be tried is carried out in the same way as the Completely Randomized Design. Consider the following case example. An experiment wanted to study the effect of Nitrogen and Variety fertilization on the yield of production carried out in the Greenhouse. Environmental conditions are assumed to be homogeneous. The fertilization factor consists of 2 levels, namely 0 kg N / ha (n0) and 60 kg N / ha (n1). The Variety Factor consists of two levels, namely the IR-64 Variety (v1) and the S-969 Variety (v2). The experiment was designed using the basic CRD design which was repeated 3 times. The experiment was a 2x2 factorial experiment so there were 4 treatment combinations: n0v1; n0v2; n1v1; and n1v2. Because it is repeated 3 times, the unit of experiment consists of 4x3 = 12 units of experiments.
Create 12 plots (experimental units) and the units of experiments are numbered from 1 to 12. The randomization step is the same as randomization on a single CRD. For example, the result of randomization is as follows:
Based on the results of the randomization, the experimental layout is as follows:
1 = n1v1 | 2 = n0v2 | 3 = n0v1 | 4 = n1v2 |
5 = n1v1 | 6 = n1v2 | 7 = n1v2 | 8 = n1v1 |
9 = n0v1 | 10 = n0v2 | 11 = n0v2 | 12 = n0v1 |
Figure 10. Trial Layout of 2 x 2 Factorial Experiment with CRD Environmental Design
Linear Model of Factorial Design In CRD
The linear model of additives for a two-factor factorial design with its environmental design a complete random design is as follows :
Yijk = μ + αi + βj + (αβ)ij + εijk
with i =1.2...,a; j = 1, 2, ..., b; c = 1, 2, ..., r
Yijk = observations on the kth experimental unit that obtained a combination of i-level treatment of factor A and j-th level of factor B
μ = population mean
αi = i-th degree effect of factor A
βj = j-th degree effect of factor B
(αβ) ij = effect of the i-th level of factor A and the j-th level of factor B
εijk = random effect of the kth experimental unit that obtained a combination of ij treatments. εij ~ N(0.σ2).
Assumption:
If all factors (factors A and B) are fixed | If all factors (factors A and B) are random |
$\begin{matrix}\sum{{\alpha}_{i}\ \ =\ {0}\ ;\ \ \ \ \ \sum{\beta}_{j}}\ =\ \mathbf{0}\ ;\ \ \ \ \ \\\sum_{{i}}{({\alpha\beta})_{{ij}}=\sum_{{j}}{({\alpha\beta})_{{ij}}=}}{0}\ ;\ \ \ \ \ {\varepsilon}_{{ijk}}\buildrel~\over~{bsi}{N}(\mathbf{0},{\sigma}^{2})\\\end{matrix}$ | $\begin{matrix}\ \ \alpha_i\buildrel~\over~N(0,{\sigma_\alpha}^2)\ \ ;\ \ \ \ \ \beta_j\buildrel~\over~N(0,{\sigma_\beta}^2)\ ;\ \ \ \ \\\ (\alpha\beta)_{ij}\buildrel~\over~N(0,{\sigma_{\alpha\beta}}^2)\ \ ;\ \ \ \ \ \ \varepsilon_{ijk}\buildrel~\over~bsiN(0,\sigma^2)\\\end{matrix}$ |
Hypothesis:
The hypotheses tested in a factorial design consisting of two factors with a Completely Randomized Design environment design are:
Hypotheses to Be Tested: | Fixed Model (Model I) | Random Model (Model II) |
Effect of AxB Interactions | ||
H0 | (αβ) ij =0 (no effect of interaction on the observed response) | σ2αβ=0 (no variance in the population of combination treatments) |
H1 | there is at least a pair (i,j) so that (αβ)ij ≠0 (there is an effect of the interaction on the observed response) | σ2αβ>0 (there is variance in the combined treatment population) |
Main effect of Factor A | ||
H0 | α1 =α2 =...=αa=0 (no response difference between the levels of factor A attempted) | σ2α=0 (no variance in the population of factor A level) |
H1 | there is at least one i so that αi ≠0 (there is a difference in response among the level of factor A tried) | σ2α>0 (there is variance in the population of factor A level) |
Main effect of Factor B | ||
H0 | β1 =β2 =...=βb=0 (no response difference between the B factor levels attempted) | σ2β=0 (no variance in the population of factor B level) |
H1 | there is at least one j so that βj ≠0 (there is a difference in response between the level of factor B that is tried) | σ2β>0 (there is variance in the population of factor B level)
|
|
|
|
Analysis of variance:
The linear model of factorial experiments with the basic design of CRD is as follows:
$$\begin{matrix}Y_{ijk}=Model+Error\\Y_{ijk}=\mu+\alpha_i+\beta_j+(\alpha\beta)_{ij}+\varepsilon_{ijk}\\Y_{ijk}={\overline{Y}}_{...}+\left({\overline{Y}}_{i..}-{\overline{Y}}_{...}\right)+\left({\overline{Y}}_{.j.}-{\overline{Y}}_{...}\right)+\left({\overline{Y}}_{ij.}-{\overline{Y}}_{i..}-{\overline{Y}}_{.j.}+{\overline{Y}}_{...}\right)+\left(Y_{ijk}-{\overline{Y}}_{ij.}\right)\\\left(Y_{ijk}-{\overline{Y}}_{...}\right)=\left({\overline{Y}}_{i..}-{\overline{Y}}_{...}\right)+\left({\overline{Y}}_{.j.}-{\overline{Y}}_{...}\right)+\left({\overline{Y}}_{ij.}-{\overline{Y}}_{i..}-{\overline{Y}}_{.j.}+{\overline{Y}}_{...}\right)+\left(Y_{ijk}-{\overline{Y}}_{ij.}\right)\\\end{matrix}$$
If the two fields are squared, then we will get:
| Definition | Calculation by Hand |
CF |
| $$\frac{Y...^2}{abr}$$ |
SSTOT | $$\sum_{i=1}\sum_{j=1}\sum_{k=1}{(Y_{ijk}-\bar{Y}...)^2}$$ | $$\sum_{i,j,k}{Y_{ijk}}^2-CF$$ |
SS(A) | $$\sum_{i=1}\sum_{j=1}\sum_{k=1}{({\bar{Y}}_{i..}-\bar{Y}...)^2}$$ | $$\sum_{i}\frac{{Y_{i..}}^2}{br}-CF$$ |
SS(B) | $$\sum_{i=1}\sum_{j=1}\sum_{k=1}{({\bar{Y}}_{.j...}-\bar{Y}...)^2}$$ | $$\sum_{j}\frac{{Y_{.j.}}^2}{ar}-CF$$ |
SS(AB) | $$\sum_{i=1}\sum_{j=1}\sum_{k=1}{({\bar{Y}}_{ij.}-{\bar{Y}}_{i...}-{\bar{Y}}_{..j.}+\bar{Y}...)^2}$$ | $$\sum_{i,j}\frac{{Y_{ij.}}^2}{r}-CF-SS(A)-SS(B)$$ |
SSE | $$\sum_{i=1}\sum_{j=1}\sum_{k=1}{({\bar{Y}}_{ijk}-{\bar{Y}}_{ij.})^2}$$ | SSTOT – SS(A) – SS(B) –JKAB |
Table 20 . Mean Squared Expectation Value of Two-Factorial Design In Complete Random Design
Source Variance | Mean Square | E(MS) | |
|
| Fixed factors A and B | Random factors A and B |
A | MS(A) | $$\sigma^2+rb\sum_{i}{\alpha_i}^2/(a-1)$$ | $$\sigma^2+r{\sigma_{\alpha\beta}}^2+rb{\sigma_\alpha}^2$$ |
B | MS(B) | $$\sigma^2+ra\sum_{j}{\beta_j}^2/(b-1)$$ | $$\sigma^2+r{\sigma_{\alpha\beta}}^2+ra{\sigma_\beta}^2$$ |
AB | MS(AB) | $$\sigma^2+r\sum_{ij}{(\alpha\beta{)_{ij}}^2}/(a-1)(b-1)$$ | $$\sigma^2+r{\sigma_{\alpha\beta}}^2$$ |
Error | MSE | $$\sigma^2$$ | $$\sigma^2$$ |
|
|
|
|
|
| Fixed A factor and random B | Fixed B factor and random A |
A | MS(A) | $$\sigma^2+r{\sigma_{\alpha\beta}}^2+rb\sum_{i}{\alpha_i}^2/(a-1)$$ | $$\sigma^2+rb{\sigma_\alpha}^2$$ |
B | MS(B) | $$\sigma^2+ra{\sigma_\beta}^2$$ | $$\sigma^2+r{\sigma_{\alpha\beta}}^2+ra\sum_{j}{\beta_j}^2/(b-1)$$ |
AB | MS(AB) | $$\sigma^2+r{\sigma_{\alpha\beta}}^2$$ | $$\sigma^2+r{\sigma_{\alpha\beta}}^2$$ |
Error | MSE | $$\sigma^2$$ | $$\sigma^2$$ |
Using the mean-squared expectation value above, we can compile its Analysis of variance Table. The anova table of factorial experiments with two factors in the Completely Randomized Design is as follows:
Table 21. Anova Table of Two-Factorial Designs in a Completely Randomized Design
Sources of variance | Degree of freedom | Sum of Squares | Mean Square | F-stat | F-table |
Treatment | ab-1 | SST | MST | MST/MSE | F(α, df-P, df-G) |
A | a-1 | SS(A) | MS(A) | MS(A)/MSE | F(α, df-A, df-G) |
B | b-1 | SS(B) | MS(B) | MS(B)/MSE | F(α, df-B, df-G) |
AB | (a-1) (b-1) | SS(AB) | MS(AB) | MS(AB)/MSE | F(α, df-AB, df-G) |
Error | ab(r-1) | SS(G) | MSE |
|
|
Total | abr-1 | SSTOT |
|
|
|
If there is an effect of the interaction, then hypothesis testing of the main effect is not necessary. Testing of the main effect will be beneficial if the effect of the interaction is not significant. Reject Ho if the value of F > Fα(df1, df2), and vice versa accept Ho.
Standard Error
The standard error for the difference among treatment averages is calculated by the following formula:
Comparison of the two averages of Factor A:
$$ SED=S_{\bar{Y}}=\sqrt{\frac{2MSE}{rb}}$$
Comparison of the two averages of Factor B:
$$ SED=S_{\bar{Y}}=\sqrt{\frac{2MSE}{ra}}$$
Comparison of the interaction of the two average factors of AxB:
$$ SED=S_{\bar{Y}}=\sqrt{\frac{2MSE}{r}}$$