[R-sig-ME] using lmer instead of 2-way repeated measure ANOVA for categorical response data

[David Duffy]

On Mon, 4 Jul 2011, Stefano Guidi wrote:

> I am trying to analyze with lmer some data I collected in a fully factorial
> 2x2 experiment, with only within-subject factors and a categorical response
> variable (binomial).
>
> model <- lmer(scelta =="antizolner" ~ phase * contrast  + (1 +
> phase*contrast| subject), data=ck1bis, family="binomial")
>
> At this point I would like to be able to plot a barchart of the
> probabilities for the 4 conditions with error bar representing 95% CI, just
> I would do with the marginal means after an ANOVA, but I don't know how to
> do it.

It seems mcmcsamp() and HPDinterval() are working OK, and I presume making 
a new 4 level variable phase_contrast would simplify the output.

Cheers, David Duffy.