Transcript POST HOC TESTS - Imperial College London
Week 2 – PART III
POST-HOC TESTS
POST HOC TESTS
• When we get a significant F test result in an ANOVA test for a main effect of a factor with more than two levels, this tells us we can reject H o • i.e. the samples are not all from populations with the same mean. • We can use post hoc tests to tell us which groups differ from the rest.
POST HOC TESTS
• There are a number of tests which can be used. SPSS has them in the ONEWAY and General Linear Model procedures • SPSS does post hoc tests on repeated measures factors, within the Options menu
Sample data
Group 1 12 14 15 13 12 2 25 22 19 18 3 13 14 17 14 23 4 24 25 23 16 34 22
Post Hoc test button
Select desired test
T e s t s o f B
7
e t w
T y
e
p f
e
e
n S u
D I I I e p S a 2 6 6 7 4 7 5 .
.
0 7
b
e
je
o r r e e r c R r r O o r o t o a r r l e S q c e c
c
t p t r
t
e e t e
s
d d e d
E
M T
f
V a o o =
f
d t
e
r ia e a .
c
l 5 l 7
t s
b 6 le ( : A S d C O t R e E d R S q ANOVA Table r e d = .
4 9 7 )
M u l t
9 5 %
i p l e
L 1 -
C o
D L e S p D
a
e
r
n
i s
d
o
e
n
n I a n c t e r
s
v a U P l 3 r ia b le : S C O R E 3 4 2 * * B a h s e e d m o e n a o n b d s e if f e r r v e e n d c m e e is a n s s ig .
n if ic Post Hoc Tests a n t a t t h e . 0 5 le v e l.
Choice of post-hoc test
• There are many different post hoc tests, making different assumptions about equality of variance, group sizes etc. • The simplest is the Bonferroni procedure
Bonferroni Test
• first decide which pairwise comparisons you will wish to test (with reasonable justification) • get SPSS to calculate t-tests for each comparison • set your significance criterion alpha to be .05 divided by the total number of tests made
Bonferroni test
• repeated measures factors are best handled this way • ask SPSS to do related t-tests between all possible pairs of means • only accept results that are significant below .05/k as being reliable (where k is the number of comparisons made)
PLANNED COMPARISONS/ CONTRASTS
• It may happen that there are specific hypotheses which you plan to test in advance, beyond the general rejection of the set of null hypotheses
PLANNED COMPARISONS
•
For example:
– a) you may wish to compare each of three patient groups with a control group – b) you may have a specific hypothesis that for some subgroup of your design – c) you may predict that the means of the four groups of your design will be in a particular order
PLANNED COMPARISONS
• Each of these can be tested by specifying them beforehand - hence planned comparisons. • The hypotheses should be orthogonal that is independent of each other
PLANNED COMPARISONS
• To compute the comparisons, calculate a t-test, taking the difference in means and dividing by the standard error as estimated from MS within from the ANOVA table
TEST OF LINEAR TREND – planned contrast
• for more than 2 levels, we might predict a constantly increasing change across levels of a factor • In this case we can try fitting a model to the data with the constraint that the means of each condition are in a particular rank order, and that they are equidistant apart.
TEST OF LINEAR TREND
• The Between Group Sum of Squares is then partitioned into two components. – the best fitting straight line model through the group means – the deviation of the observed group means from this model
TEST OF LINEAR TREND
• The linear trend component will have one degree of freedom corresponding to the slope of the line. • Deviation from linearity will have (k-2) df.
• Each of these components can be tested, using the Within SS, to see whether it is significant.
TEST OF LINEAR TREND
• If there is a significant linear trend, and non-significant deviation from linearity, then the linear model is a good one.
• For k>3, The same process can be done for a quadratic trend - a parabola is fit to the means. For example, you may be testing a hypothesis that as dosage level increases, the measure initially rises and then falls (or vice versa).
1 1 1 1 S t d
R
.
TEST OF LINEAR TREND
e p o
S
r
C 4 1 2 3 T
t
O .
.
.
0 0 0 o 0 0 0 t a R l E P
TEST OF LINEAR TREND
16 14 12 22 20 18 10 8 6 4 2 0 Group 1 Group 2 GROUP Group 3
TEST OF LINEAR TREND
A
1 1 2
N
M 0 7 2
O
.
S 2
V A
u m
T a
it
b
h
l e
e in y G r f ) r R p O s U L P in it y