[R-sig-ME] Characterizing correlation between binomial effects
David Duffy
David.Duffy at qimr.edu.au
Mon Sep 10 05:53:59 CEST 2012
On Wed, 5 Sep 2012, Mike Lawrence wrote:
> I have data from multiple human participants who each completed two
> tasks. In each task there were two conditions, and in each condition
> there are multiple trials that yield a binomial (1/0) outcome. I would
> like to evaluate whether there is a correlation between Ss' condition
> effect in task 1 and their condition effect in task 2.
>
> One idea I had was to combine the two task data frames to a single
> data frame, code task as a variable, and fit a full glmm model:
>
> fit = glmer(
> family = binomial
> , data = both_tasks
> , formula = response ~ task*condition + (1+task*condition|participant)
> )
If I understand this correctly, which I probably don't, you want something
like a P-value for the likelihood ratio test comparing
(1+task*condition|participant) with (1 + condition + task|participant). Or
should that be (1|condition) + (1 + task|participant)?
More information about the R-sig-mixed-models
mailing list