Statadmin: Created page with '==Problem Statement== ''A 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 bo…'

2010-07-31T18:58:16Z

Created page with '==Problem Statement== ''A 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 bo…'

New page

==Problem Statement==
''A 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 bolts of 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 data from this experiment (use <math>\alpha = 0.05</math>) and draw appropriate conclusions.

<table border=1 cellspacing=0 cellpadding=4>
<tr><th bgcolor="#eeeeee"></th><th colspan=5 bgcolor="#eeeeee">Bolt</th></tr>
<tr><th bgcolor="#eeeeee">Chemical</th><th bgcolor="#eeeeee">1</th><th bgcolor="#eeeeee">2</th><th bgcolor="#eeeeee">3</th><th bgcolor="#eeeeee">4</th><th bgcolor="#eeeeee">5</th></tr>
<tr><th bgcolor="#eeeeee">1</th><td>73</td><td>68</td><td>74</td><td>71</td><td>67</td></tr>
<tr><th bgcolor="#eeeeee">2</th><td>73</td><td>67</td><td>75</td><td>72</td><td>70</td></tr>
<tr><th bgcolor="#eeeeee">3</th><td>75</td><td>68</td><td>78</td><td>73</td><td>68</td></tr>
<tr><th bgcolor="#eeeeee">4</th><td>73</td><td>71</td><td>75</td><td>75</td><td>69</td></tr>
</table>

==Solution==
[[Image:4-1dataplot.png|thumb|360px|left|'''Figure 1:''' Illustration of data.]]

In this problem we need to compare <math>a = 4</math> treatments (chemical agents), but the data is being influenced by <math>b = 5</math> blocks (different bolts of fabric). Each datapoint can be thought of as the sum of several factors:<br>
<math> y_{ij} = \mu + \tau_i + \beta_j + \epsilon_{ij} </math><br>
Where:<br>
<math> i = 1,2,...,a </math><br>
<math> j = 1,2,...,b </math><br>
and:<br>
<math> y_{ij} = </math> datapoint from the <math>i</math>th treatment and the <math>j</math>th block<br>
<math> \mu = </math> the grand mean<br>
<math> \tau_i = </math> the effect of the <math>i</math>th treatment<br>
<math> \beta_j = </math> the effect of the <math>j</math>th block<br>
<math> \epsilon_{ij} = </math> a random error affecting the datapoint from the <math>i</math>th treatment and the <math>j</math>th block<br>

In the left plot of Figure 1, this equation is illustrated. The blue line marks the grand sample mean <math>\bar{y}.. \approx \mu</math>, and the red X marks the treatment mean <math>\bar{y}_i.</math> (mean of all data treated by the same chemical). The green number marks where the datapoint would be if the treatment effect and block effect were present, but there was no random error. Finally, the black number marks the location of the datapoint – the sum of the grand mean and each effect. The value of the numbers indicates which block that point was affected by.

In order to examine only the effect of the chemical agent, each bolt was divided into four pieces, each of which was treated by a different chemical agent. This serves to average out the effect of the different bolts. This is displayed in the left plot in Figure 1. The green lines deviate from each treatment mean identically, illustrating that the block effect is constant across chemical treatments.

The plot on the right in Figure 1 illustrates the data similarly, but the squares indicate the block means, and the treatment effects are shown deviating from them. This is merely an alternate visualization of the sum of the grand mean and each effect.

===Hypothesis Testing===

As usual, our null hypothesis will be that all the treatment means are equal (and come from the same distribution). Our alternative hypothesis is that they are not all equal:

<math>H_0\!:~\tau_1 = \tau_2 = ... = \tau_a = 0 </math><br>
<math>H_1\!:~\tau_i \ne 0 </math>

In order to test these hypotheses, we are interested in the sums of squares of all the error effects.

<math> SS_{Total} = \sum_{i=1}^a \sum_{j=1}^b y^2_{ij} - \frac{y^2\!\!..}{N} = 191.75</math><br>
<math> SS_{Treatments} = \frac{1}{b}\sum_{i=1}^a y^2_i\!. - \frac{y^2\!\!..}{N} = 12.95</math><br>
<math> SS_{Blocks} = \frac{1}{a} \sum_{j=1}^b y^2\!\!._j - \frac{y^2\!\!..}{N} = 157.00</math><br>
<math> SS_{Error} = SS_{Total} - SS_{Treatments} - SS_{Blocks} = 21.80</math>

Each of these is a measure of its respective type of error. Graphically, in the left plot of Figure 1 <math>SS_{Total}</math> is the sum of the squares of the distance along the x-axis from the grand mean to each of the black numbers. <math>SS_{Treatments}</math> is similar, but for the x-components of the red lines. For <math>SS_{Blocks}</math> we are summing the squares of the x-component of the green lines. And for <math>SS_{Error}</math> it is the x-component of the black lines that are squared and summed (it is easier to calculate this by arithmetic as in the equation above).

By dividing the statistics above by their respective degrees of freedom, we get a normalized measure of the error effects, or a ''mean sum of squares''.

<math>MS_{Treatments} = \frac{SS_{Treatments}}{a-1} = 4.32</math><br>
<math>MS_{Blocks} = \frac{SS_{Blocks}}{b-1} = 39.25</math><br>
<math>MS_{Error} = \frac{SS_{Error}}{(a-1)(b-1)} = 1.81</math>

In order to test our hypothesis, we must compare the treatment error to random error. We will do this with the F statistic, <math>F_0 = \frac{MS_{Treatments}}{MS_{Error}} = 2.38</math>. We compare this quantity to a <math>F_{\alpha,~a-1,~(a-1)(b-1)} = 3.49</math>, which comes from the F-distribution (this can be calculated in Excel with <tt>FINV(α,a−1,N−a)</tt>; in R use <tt>qf(1-alpha, a-1, (a-1)*(b-1))</tt>. If <math>F_{0} > F_{\alpha,~a-1,~(a-1)(b-1)}</math>, we reject the null hypothesis. As this is not the case, we cannot reject our null hypothesis.

<table cellspacing=0 cellpadding=5 style="border-top: 1px solid black; border-bottom: 1px solid black">
<tr align="center">
<th bgcolor="#eeeeee" style="border-bottom: 1px solid black;"> Source of Variation </th>
<th bgcolor="#eeeeee" style="border-bottom: 1px solid black;"> Sum of Squares </th>
<th bgcolor="#eeeeee" style="border-bottom: 1px solid black;"> Degrees of Freedom </th>
<th bgcolor="#eeeeee" style="border-bottom: 1px solid black;"> Mean Square </th>
<th bgcolor="#eeeeee" style="border-bottom: 1px solid black;"> <math>F_0</math> </th>
<th bgcolor="#eeeeee" style="border-bottom: 1px solid black;"> <math>F_{\alpha}</math> </th>
<th bgcolor="#eeeeee" style="border-bottom: 1px solid black;"> P-Value </th>
</tr>
<tr align="center">
<th bgcolor="#eeeeee"> Treatments (Chemical agent) </th>
<td> 12.95 </td>
<td> 3 </td>
<td> 4.317 </td>
<td> 2.376 </td>
<td> 3.490 </td>
<td> 0.121 </td>
</tr>
<tr align="center">
<th bgcolor="#eeeeee"> Blocks (Bolts) </th>
<td> 157.00 </td>
<td> 4 </td>
<td> 39.250 </td>
</tr>
<tr align="center">
<th bgcolor="#eeeeee"> Error </th>
<td> 21.80 </td>
<td> 12 </td>
<td> 1.817 </td>
</tr>
<tr align="center">
<th bgcolor="#eeeeee" style="border-top: 1px solid black;"> Total </th>
<td style="border-top: 1px solid black;"> 191.75 </td>
<td style="border-top: 1px solid black;"> 19 </td>
<td style="border-top: 1px solid black;"></td>
<td style="border-top: 1px solid black;"></td>
<td style="border-top: 1px solid black;"></td>
<td style="border-top: 1px solid black;"></td>
</tr>
</table>

Problem 4-1, Randomized blocks - Revision history

Statadmin: Created page with '==Problem Statement== ''A 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 bo…'