[R] Problem with 2-ways ANOVA interactions

Frodo Jedi frodo.jedi at yahoo.com
Thu Jan 6 00:10:24 CET 2011


Dear All, 
I have a problem in understanding how the interactions of 2 ways ANOVA work, 
because I get conflicting results
from a t-test and an anova. For most of you my problem is very simple I am sure.

I need an help with an example, looking at one table I am analyzing. The table 
is in attachment
and can be imported in R by means of this command: 
scrd<- 
read.table('/Users/luca/Documents/Analisi_passi/Codice_R/Statistics_results_bump_hole_Audio_Haptic/tables_for_R/table_realism_wood.txt',
 header=TRUE, colClasse=c('numeric','factor','factor','numeric')) 


This table is the result of a simple experiment. Subjects where exposed to some 
stimuli and they where asked to evaluate the degree of realism
of the stimuli on a 7 point scale (i.e., data in column "response").
Each stimulus was presented in two conditions, "A" and "AH", where AH is the 
condition A plus another thing (let´s call it "H").

Now, what means exactly in my table the interaction stimulus:condition?

I think that if I do the analysis anova(response ~ stimulus*condition) I will 
get the comparison between 

the same stimulus in condition A and in condition AH. Am I wrong?

For instance the comparison of stimulus flat_550_W_realism presented in 
condition A with the same stimulus, flat_550_W_realism, 

presented in condition AH.

The problem is that if I do a t-test between the values of this stimulus in the 
A and AH condition I get significative difference,
while if I do the test with 2-ways ANOVA I don´t get any difference. 
How is this possible?

Here I put the results analysis


#Here the result of ANOVA:
> fit1<- lm(response ~ stimulus + condition + stimulus:condition, data=scrd) 
>#EQUIVALE A lm(response ~ stimulus*condition, data=scrd) 
>
> anova(fit1)
Analysis of Variance Table

Response: response
                    Df Sum Sq Mean Sq F value   Pr(>F)    
stimulus             6  15.05   2.509  1.1000   0.3647    
condition            1  36.51  36.515 16.0089 9.64e-05 ***
stimulus:condition   6   1.47   0.244  0.1071   0.9955    
Residuals          159 362.67   2.281                     
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 


#As you can see the p-value for stimulus:condition is high.


#Now I do the t-test with the same values of the table concerning the stimulus 
presented in A and AH conditions:

flat_550_W_realism                         =c(3,3,5,3,3,3,3,5,3,3,5,7,5,2,3) 
flat_550_W_realism_AH                   =c(7,4,5,3,6,5,3,5,5,7,2,7,5,   5) 

> t.test(flat_550_W_realism,flat_550_W_realism_AH, var.equal=TRUE) 

Two Sample t-test 

data:  flat_550_W_realism and flat_550_W_realism_AH 
t = -2.2361, df = 27, p-value = 0.03381 
alternative hypothesis: true difference in means is not equal to 0 
95 percent confidence interval: 
-2.29198603 -0.09849016 
sample estimates: 
mean of x mean of y 
3.733333  4.928571 


#Now we have a significative difference between these two stimuli (p-value = 
0.03381) 



Why I get this beheaviour?


Moreover, how by means of ANOVA I could track the  significative differences 
between the stimuli presented in A and AH  condition
whitout doing the t-test? 

Please help!

Thanks in advance


      
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: table_realism_wood.txt
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110105/8e3524de/attachment.txt>


More information about the R-help mailing list