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.