Re: [R] Equivalent of intervals() in lmer

2008-04-23 Thread Michael Kubovy
Dear Friends, (1) There may be a solution for those (e.g., experimental psychologists) who are *not at all* interested in generalizing the absolute level of the response variable (say, reaction time, rt) to other subjects, but *only* to generalize the effect of the manipulated variables wi

Re: [R] Equivalent of intervals() in lmer

2008-04-21 Thread Michael Kubovy
Sorry, I meant to say: "For the moment I wonder if the solution is not to use CIs based on the two low SEs produced by the ~ time model, and to treat them as least-significant difference intervals." _ Professor Michael Kubovy University of Virginia Department of Ps

Re: [R] Equivalent of intervals() in lmer

2008-04-21 Thread Michael Kubovy
Thanks Doug, You write: "If you want to examine the three means then you should fit the model as lmer(rcl ~ time - 1 + (1 | subj), fr)" I do just that (which is what Dieter just sent). But the CIs are much too big compared to the CIs for differences between means (which should be bigger tha

Re: [R] Equivalent of intervals() in lmer

2008-04-21 Thread Dieter Menne
Douglas Bates stat.wisc.edu> writes: > If you want to examine the three means then you should fit the model as > lmer(rcl ~ time - 1 + (1 | subj), fr) > True, but for the notorious "error bars" in plots that reviewers always request the 0.35 is probable more relevant than the 1.87. Which I thin

Re: [R] Equivalent of intervals() in lmer

2008-04-21 Thread Douglas Bates
On 4/21/08, Michael Kubovy <[EMAIL PROTECTED]> wrote: > To help Kedar a bit: > > Here is one way: > > recall <- c(10, 13, 13, 6, 8, 8, 11, 14, 14, 22, 23, 25, 16, 18, 20, > 15, 17, 17, 1, 1, 4, 12, 15, 17, 9, 12, 12, 8, 9, 12) > fr <- data.frame(rcl = recall, time = factor(rep(c(1, 2, 5), 10)),

Re: [R] Equivalent of intervals() in lmer

2008-04-21 Thread Michael Kubovy
To help Kedar a bit: Here is one way: recall <- c(10, 13, 13, 6, 8, 8, 11, 14, 14, 22, 23, 25, 16, 18, 20, 15, 17, 17, 1, 1, 4, 12, 15, 17, 9, 12, 12, 8, 9, 12) fr <- data.frame(rcl = recall, time = factor(rep(c(1, 2, 5), 10)), subj = factor(rep(1:10, each = 3))) (fr.lmer <- lmer(rcl ~ time +

Re: [R] Equivalent of intervals() in lmer

2008-04-20 Thread Dieter Menne
kedar nadkarni gmail.com> writes: > I have been trying to obtain confidence intervals for the fit after having > used lmer by using intervals(), but this does not work. intervals() is > associated with lme but not with lmer(). What is the equivalent for > intervals() in lmer()? ci in Gregory

[R] Equivalent of intervals() in lmer

2008-04-20 Thread kedar nadkarni
Hi all, I have been trying to obtain confidence intervals for the fit after having used lmer by using intervals(), but this does not work. intervals() is associated with lme but not with lmer(). What is the equivalent for intervals() in lmer()? I could not get this information from the documentat