Description on the average comparison sub topic:
- Hypothesis Test
- Error Type
- Level of significant
- Analysis of variance (Anova)
- Average comparison
- The number of possible combinations of pairs
- Problems in Multiple Comparison
Introduction
Hypothesis Test
Hypothesis Test Procedure
- Determining H0 and HA
- Determining its statistical value
- Confidence level
- Degrees of freedom
- Number of samples
- Determines the value of the calculated statistic, this value depends on the parametric method used..
- Making decisions, this is determined by comparing the calculated statistical value with the statistical value of the table or its critical value.
Error Type
Conclusion | The Significant State of Affairs | |
True Hypothesis | False Hypothesis | |
Accept Hypothesis | True | Wrong (type II error) |
Reject Hypothesis | Wrong (type I error) | True |
- If we reject the hypothesis (H0 is rejected), while the true hypothesis says we are committing a type I error
- Conversely if we accept the hypothesis (H0 is accepted) while the hypothesis is wrong, it says we committed a type II error
Level of Significant
- The probability of making a Type I error is known as the significant level or level of significance and is usually expressed by α.
- Commonly used values are α = 0.05 and α = 0.01.
- Example: α = 0.05
- α = 0.05 is often also called the significant rate of 5% which means roughly 5 out of 100 conclusions that we would reject a hypothesis that should be accepted.
- This means we may be wrong with a 5% chance or in other words: roughly 95% are sure that we have made the correct conclusion
Analysis of Variance (Anova)
Table of Analysis of Variance for Nitrogen Content
Variance Sources | DF | SS | MS | F-stat | F .01 |
Treatment (P) | 5 | 847.046667 | 169.409333 | 14.37 ** | 3.895 |
Error | 24 | 282.928 | 11.7886667 | - |
|
Total | 29 | 1129.97467 |
|
|
|
cv = 17.27 %; |
- Test F (Anova) was used to determine whether there were differences between the treatment means
- H0: m1 = m2 = m3 = . . . = mk
- HA: one or more of the mean treatments differ from others
- When H0 is accepted, we conclude that all treatment means do not differ
- If H0 is rejected (HA is accepted), we conclude that there are one or more different treatment means from others, but we cannot conclude that all treatment means are different.
- If H0 is rejected, which of the mean treatments is different from others??
- Post hoc test => to find out the mean of which treatments are different and how they are in order
Treatment (P) | Mean N Content |
3Dok1 | 28.82 |
3Dok5 | 23.98 |
3Dok4 | 14.64 |
3Dok7 | 19.92 |
3Dok13 | 13.26 |
Combined | 18.70 |
Mean comparison
- Planned comparison
- It has been planned before the data (experimental results) are obtained
- Linear Contrasts (Complex Comparisons)
- Scheffé's Test
- Comparison with controls:
- Dunnet
- Bonferroni
- Sidak
- Unplanned comparison – post hoc test
- Pair-wise comparisons or Multiple Comparisons:
- LSD (not recommended)
- Tukey HSD (recommended)
- Scheffé
- Bonferroni
- Sidak
- Gabriel
- Hochberg
- Multistage test (Multiple Range Test): comparison of all combinations of mean pairs determined after the results of the experiment were obtained (Post-Hoc Comparison)
- SNK (Student Newman Keul)
- Duncan
- Tukey HSD
- Tukey B
- Scheffé
- Gabriel
- REGWQ (Ryan, Einot, Gabriel and Welsh. Q = the studentized range statistic) recommended if supporting software is available)
- Pair-wise comparisons or Multiple Comparisons:
The number of possible combinations of pairs
In general, the many possible combinations of comparisons for the level of treatment are:
$$\frac{k(k-1)}{2}$$
In the example above, k=3, so that the number of possible combinations:
$$\frac{3(3-1)}{2}=\frac{3(2)}{2}=3$$
Problems in Multiple Comparison
- Any pairwise comparison at a certain significant level (e.g., a = 5%) will result in a type I error a
- Example: a= 0.05
- a = 0.05 is often also called the significant rate of 5% which means roughly 5 out of 100 conclusions that we would reject a hypothesis that should be accepted (type I error);
- This means we may be wrong with a 5% chance or in other words: roughly 95% are sure that we have made the correct conclusion
- When we compare 20 treatments at a significant level of 0.05, where all the null hypotheses are actually significant:
- 0.9520 = 0.3585 and
- Type I error = 1-0.3585 = 0.6415
- In general, the odds of making a type I error on the comparison j are:
- $\alpha_{EW}=1-(1-\alpha_{PC})^j$
- It is known as Experiment wise (EW), or Familywise type I error rate.
- Suppose we want to do 3 treatment comparisons at a significant level: aPC = 0.05. The probability to make Type I Mistakes are:
$\begin{matrix}\alpha_{EW}=1-(1-\alpha_{PC})^j\\=1-(1-0.05)^3\\=1-(0.95)^3\\=1-0.8574\\=0.1426\\\end{matrix}$
General Procedure Comparison of Treatments
Suppose the mean values of t: ${\bar{x}}_1,{\bar{x}}_2,\cdots,\bar{x}t$
- F-test showed significant differences between treatment means
- Perform an analysis (post hoc test) to determine exactly where there is a difference
- Two means are expressed differently when the difference is greater than the critical value
Tips in Comparison means
- Select a pairwise comparison when no treatment will be compared to the control.
- Use comparison with controls when we are going to compare treatment with controls (Dunnet).
Presentation of Treatment Mean Table
In addition to the methods already mentioned above, there are several ways of presenting differences between treatment means, namely in the form of crosstabulation of mean differences or odds, giving the same line to the mean of different treatments, and in the form of graphs.
In the difference in the difference in the mean value of the treatment:
Treatment | Mean | 28.82 | 23.98 | 14.64 | 19.92 | 13.26 | 18.7 |
3Dok1 | 28.82 | 4.84 | 14.18 | 8.9 | 15.56 | 10.12 | |
3Dok5 | 23.98 | 4.84 | 9.34 | 4.06 | 10.72 | 5.28 | |
3Dok4 | 14.64 | 14.18 | 9.34 | 5.28 | 1.38 | 4.06 | |
3Dok7 | 19.92 | 8.9 | 4.06 | 5.28 | 6.66 | 1.22 | |
3Dok13 | 13.26 | 15.56 | 10.72 | 1.38 | 6.66 | 5.44 | |
Combined | 18.7 | 10.12 | 5.28 | 4.06 | 1.22 | 5.44 |
Description: if the mean difference is >4.482 mg, then the two treatments are significantly different
In the probability form:
Treatment | Mean | {1} | {2} | {3} | {4} | {5} | {6} |
3Dok1 | 28.82 |
| 0.035 | 0.000 | 0.000 | 0.000 | 0.000 |
3Dok5 | 23.98 | 0.035 |
| 0.000 | 0.074 | 0.000 | 0.023 |
3Dok4 | 14.64 | 0.000 | 0.000 |
| 0.023 | 0.531 | 0.074 |
3Dok7 | 19.92 | 0.000 | 0.074 | 0.023 |
| 0.005 | 0.579 |
3Dok13 | 13.26 | 0.000 | 0.000 | 0.531 | 0.005 |
| 0.019 |
Combined | 18.7 | 0.000 | 0.023 | 0.074 | 0.579 | 0.019 |
|
Description: if the probability value < α=0.05, then the two treatments are significantly different
In the notation form:
No. | Treatment | Mean | 1 | 2 | 3 | 4 | 5 | Notation |
5 | 3Dok13 | 13.26 | 13.26 | a | ||||
3 | 3Dok4 | 14.64 | 14.64 | 14.64 | ab | |||
6 | Combined | 18.70 | 18.70 | 18.70 | bc | |||
4 | 3Dok7 | 19.92 | 19.92 | 19.92 | cd | |||
2 | 3Dok5 | 23.98 | 23.98 | d | ||||
1 | 3Dok1 | 28.82 | 28.82 | e |
After the order is returned based on the order of treatment:
No. | Treatment | Mean | |
1 | 3Dok1 | 28.82 | e |
2 | 3Dok5 | 23.98 | d |
3 | 3Dok4 | 14.64 | ab |
4 | 3Dok7 | 19.92 | cd |
5 | 3Dok13 | 13.26 | a |
6 | Combined | 18.70 | bc |
Description: the number followed by the same letter in one column does not differ markedly at a significant level of 5% according to the LSD test