[R] Coeficients estimation in a repeated measures linear model

[Bert Gunter]

Sergio:

1. You do not have a "repeated measures linear model" .

2. This list is not designed to replace your own efforts to learn the
necessary R background, in this case, factor coding and contrasts in linear
models. I would suggest you spend some time with any of the many fine R
linear model tutorials that can be found on the web. Here is one place to
look for suggestions: https://www.rstudio.com/online-learning/#R  . But
just googling around you'll probably find something that may suit even
better.

3. This list is primarily for R programming help, not statistics help
(although they do sometimes intersect). For the latter, try a statistics
site like stats.stackexchange.com  .

4. Finally, as always, consulting with a local statistical resource, if
available, is always worth considering.

HTH.

Cheers,
Bert



Bert Gunter

"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )

On Wed, Dec 6, 2017 at 6:17 AM, Sergio PV <serpalma.v at gmail.com> wrote:

> Dear Users,
>
> I am trying to understand the inner workings of a repeated measures linear
> model. Take for example a situation with 6 individuals sampled twice for
> two conditions (control and treated).
>
> set.seed(12)
> ctrl <- rnorm(n = 6, mean = 2)
> ttd <- rnorm(n = 6, mean = 10)
> dat <- data.frame(vals = c(ctrl, ttd),
>                   group = c(rep("ctrl", 6), rep("ttd", 6)),
>                   ind = factor(rep(1:6, 2)))
>
> fit <- lm(vals ~ ind + group, data = dat)
> model.matrix(~ ind + group, data = dat)
>
> I am puzzled on how the coeficients are calculated. For example, according
> to the model matrix, I thought the intercept would be individual 1 control.
> But that is clearly not the case.
> For the last coeficient, I understand it as the mean of all differences
> between treated vs control at each individual.
>
> I would greatly appreciate if someone could clarify to me how the
> coefficients in this situation are estimated.
>
> Thanks
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/
> posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

	[[alternative HTML version deleted]]