[R] Contrasts for 2x4 interaction in mixed effects model

David Winsemius dwinsemius at comcast.net
Tue Sep 11 09:14:36 CEST 2012


On Sep 10, 2012, at 5:59 PM, semperparatus wrote:

> I want to change it because I don't want to compare in this instance between
> conditions, but I simply want to see the contrast t-statistic between
> patient and control at every level of condition (1, 2, and 3).
> 
>> From there I'd like to be able to plot the t-statistic for the contrast
> between patient and control at level 1 of conditon, level 2 of condition,
> and level 3 of condition, each with error bars.
> 
> In the post I responded to the output gave fixed effect output for
> SecA:Fir1, SecB:Fir1, SecC:Fir1, and SecD:Fir1.  I'm hoping to get the same
> sort of output but for mine it would be Cond1:Patient1,  Cond2:Patient1,
> Cond3:Patient1.

It does not appear that you have the same situation as was being discussed earlier: 

Yours was:

' *When I tried using the syntax you used with my model:
lmer(H.y. ~ patient*stance*cond  +(cond/patient) + (1|subj), data=H), I got this
result, which seems to be using condition 1 as a part of the baseline. Any
idea how to change that?*'

The other was:

test <- lmer(Latency ~ (Nuisance1*Nuisance2) + (Sec/Fir) + (1|Subject) +
(1|Item), datatotest)

He had separated his nuisance parameters from the 2 variables (Sec and Fir)  for which he was interested in examining contrasts.

PLEASE learn to include context.



-- 
David Winsemius, MD
Alameda, CA, USA




More information about the R-help mailing list