chemist wishes to test the effect of four chemical agents on the strength of a particular type of cloth. Because there might be variability from one bolt to another, the chemist decides to use a randomized block design, with the boltsof cloth considered as blocks. She selects five bolts and applies all four chemicals in random order to each bolt. The resulting tensile strengths follow. Analyze the data from this experiment (use α=0.15) and draw appropriate conclusions. Be sure to state the linear effects model and hypotheses being tested.
library(GAD)## Loading required package: matrixStats## Loading required package: R.methodsS3## R.methodsS3 v1.8.2 (2022-06-13 22:00:14 UTC) successfully loaded. See ?R.methodsS3 for help.Chemical_1<-c(73,68,74,71,67)
Chemical_2<-c(73,67, 75,72,70)
Chemical_3<-c(75,68,78,73,68)
Chemical_4<-c(73,71,75,75,69)
obs<-c(Chemical_1,Chemical_2,Chemical_3,Chemical_4)
Blots<-c(rep(seq(1,5),4))
chemical<-c(rep(1,5),rep(2,5),rep(3,5),rep(4,5))
chemical<-as.fixed(chemical)
Blots<-as.fixed(Blots)
model<-lm(obs~chemical+Blots)
gad(model)## Analysis of Variance Table
##
## Response: obs
## Df Sum Sq Mean Sq F value Pr(>F)
## chemical 3 12.95 4.317 2.3761 0.1211
## Blots 4 157.00 39.250 21.6055 2.059e-05 ***
## Residual 12 21.80 1.817
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Comment:
With Alpha =.15, we came to the conclusion that there was insufficient data to reject the null hypothesis that the chemicals had the same effect on the fabric.
Assume now that she didn’t block on Bolt and rather ran the experiment at a completely randomized design on random pieces of cloth, resulting in the following data. Analyze the data from this experiment (use α=0.15) and draw appropriate conclusions. Be sure to state the linear effects model and hypotheses being tested.
chemical<-c(rep(1,5),rep(2,5),rep(3,5),rep(4,5))
chemical<-as.fixed(chemical)
model<-lm(obs~chemical)
gad(model)## Analysis of Variance Table
##
## Response: obs
## Df Sum Sq Mean Sq F value Pr(>F)
## chemical 3 12.95 4.3167 0.3863 0.7644
## Residual 16 178.80 11.1750Linear effects model observations :
\[ y_{i,j}=\mu+\tau_{i,j}+\epsilon_{i,j} \]
Since the p- value is greater than 0.15, we concluded that we accept the null hypotheses that the chemicals have the same effect.
Comment on any differences in the findings from questions 1 and 2. Do you believe that the Bolt of cloth represents a significant amount of nuisance variability.
Comment: Blocking causes the MSE to rise, therefore removing the barrier eliminates the annoyance and makes a big impact. Blocking has a p-value of 0.1211, whereas not blocking has a p-value of 0.7644.