[R-sig-ME] main effect significance disappears when interaction is modeled

Francesco Sigona francesco.sigona at unisalento.it
Fri Oct 30 13:48:55 CET 2015


Hi all,

I've searching the list but I could't find an answer.

In my design, where all my 3 fixed effects are categorical (binary), I 
found that HEIGHT is significant (also anova tell me the same), when I 
model main effects only.

STUDYING MODEL: X ~ HEIGHT + LOWNESS + PLACE + (1|Soggetto)
               Estimate Std. Error    t value
(Intercept) -52.204545   1.277885 -40.852292
HEIGHThigh    3.681818   1.475575   2.495175
LOWNESSlow   -7.386364   1.952002  -3.783994
PLACEback    -0.500000   1.475575  -0.338851

also, I get that  :
LOWNESS affected X (Chi2(1)=13.61, p=0.00023) adding about -7.39 ± 2 
(standard error) to it
HEIGHT affected X (Chi2(1)=6.335, p=0.012) adding about 3.68 ± 1.5 
(standard error) to it

(despite the names, LOWNESS is not the opposite of HEIGHT)

Then, if try to model the interaction between PLACE and HEIGHT this 
comes out:
STUDYING MODEL WITH INTERACTIONS:  X~LOWNESS + HEIGHT*PLACE + (1|Soggetto)
                        Estimate Std. Error    t value
(Intercept)          -48.909091   1.162784 -42.062068
LOWNESSlow            -4.090909   1.520004  -2.691380
HEIGHThigh            -2.909091   1.520004  -1.913870
PLACEback             -7.090909   1.520004  -4.665058
HEIGHThigh:PLACEback  13.181818   2.149611   6.132188
************HEIGHT*PLACE affected X (Chi2(1)=29.17, p=6.6e-08)

(then I am supposed to report the post-hoc that I would get by means of 
multcomp package.)

This posits the following questions to me:

a) Does HEIGHT is still significant? If I make an anova() comparison 
between X~LOWNESS + HEIGHT + PLACE + HEIGHT:PLACE +  (1|Soggetto) and 
X~LOWNESS + PLACE + HEIGHT:PLACE +  (1|Soggetto),    anova() thinks that 
the two models are the same.

b) Should I report that I found a main effect for HEIGHT (with 3.681818 
estimates) and an interaction with PLACE, or should I report that I 
found just the interaction?

c) Since I find the interaction, in this model the main effect of LOW is 
still significant but Chi2square related stuff changes to Chisq=7.3232, 
p=0.0068), as well the slope (now is -4.09). So , what I am supposed to 
report? WHat exactly you would report in my case?

I really need some direction about this, thank you in advance.

Francesco


-- 
*Francesco SIGONA*
Electronics engineer 	
Piazza Filippo Muratore
73100 - Lecce - Italy 	<https://maps.google.com/maps?q=40.331002,18.156462>
tel.: +39 0832 335006
fax.: +39 0832 335007

============================================================
*Center for Interdisciplinary Research on Language (CRIL) 
<http://www.cril.unisalento.it> &
Cognitive Neuroscience of Language and Speech Sciences Lab (CNLSS) * 	
*Dipartimento di Studi umanistici
Università del Salento * 	

============================================================
*Laboratorio Diffuso di Ricerca Interdisciplinare Applicata alla Medicina
(DReAM) * 	



More information about the R-sig-mixed-models mailing list