[R-sig-ME] [R] understanding I() in lmer formula
Ben Bolker
bbolker at gmail.com
Sat Jun 17 20:23:26 CEST 2017
On Sat, Jun 17, 2017 at 12:29 PM, Don Cohen <don-r-help at isis.cs3-inc.com> wrote:
> Ben Bolker writes:
> > For the level of detail you're getting into, it would be a really good
> > idea to read the paper that accompanies the lme4 package:
> > vignette("lmer",package="lme4") . This goes into a lot of detail
> > about the theory and data structures ...
>
> vignette("lmer",package="lme4")
> gives me
> Warning message:
> vignette 'lmer' not found
That's surprising ... what's packageVersion("lme4") ?
>
> Is this the same as
> https://cran.r-project.org/web/packages/lme4/vignettes/lmer.pdf ?
Yes.
>
> I think that's the same one that I was having trouble with before
> and gave up around eqn 15.
> Although, I had the impression that it (looks like eqn 14) was
> describing what I expected and asked about in a previous message,
> namely paying only once for each group and then once for each data
> point within the group.
>
> Are you saying that the vignette link above actually answers the
> questions in this last message about how to compute the loglik of
> the model? It doesn't look to me like it will.
> I view my current line of questions (and I have many more, but am
> trying to resist bombarding you with all at once) as a way to get
> the background I'll need to get through that paper.
I would also recommend checking out Pinheiro and Bates (2000),
which is a full (book-length) treatment of the same topic, so is a
tiny bit more discursive/friendlier ...
>
> In fact, one of my questions when I read that paper was what
> correlated vs uncorrelated intercept and slope meant - I didn't
> see any explanation. I think that Emmanuel Curis has now explained
> that, but I'm still trying to check my understanding.
>
> Since I'm writing anyway, let me indulge in one more question about
> the formulas. Since (x|g) means correlated intercept and slope for
> x, does (x+y|g) include a separate correlation between x slope and
> y slope? That is, the cost of specifying a group would involve a
> 3 dimensional normal distribution over intercept,x,y ?
Yes. (So in particular this model are 3*(3+1)/2=6 parameters
(equivalent to s^2{1}, s^2(x),s^2(y), cov(1,x), cov(1,y), cov(x,y))
More information about the R-sig-mixed-models
mailing list