Thanks, Bert.

2012/8/17 Bert Gunter <gunter.berton@gene.com>

> 1. Post on R-sig-mixed-models, not here.
>
> 2. Models 2 and 3 make no sense (to me, anyway). What do you think
> they mean? (Don't answer here -- explain on the mixed models list).
>
> -- Bert
>
> On Thu, Aug 16, 2012 at 9:13 PM, li li <hannah.hlx@gmail.com> wrote:
> > Dear all,
> >   I am trying to use R to fit mixed models.
> >   Take the following example, where ind is a random effect and
> > sample is fixed. I wanted to fit
> >    Model 1: values = ind + sample
> >    Model 2: values =ind * sample
> >    Model 3:  values=ind(sample) + sample
> >
> > Tried to use the below for mod1, but it did not work.
>  Can anyone give
> > some help on this. Thanks so much.
> > mod1 <- lme(values ~ sample +
> >                      random=~ 1|ind, data=y)
>
> Your syntax is wrong. Re-read ?lme
> (random is a separate argument, not part of the formula).
>  >
> >
> >
> >
> >
> >
> >
> >
> >        values ind sample
> > 1  0.03325   1      1
> > 2  0.03305   1      1
> > 3  0.03185   1      1
> > 4  0.03515   1      1
> > 5  0.03375   1      1
> > 6  0.01180   1      2
> > 7  0.01850   1      3
> > 8  0.02915   1      4
> > 9  0.06200   1      5
> > 10 0.03230   2      1
> > 11 0.03345   2      1
> > 12 0.03385   2      1
> > 13 0.03605   2      1
> > 14 0.03225   2      1
> > 15 0.01145   2      2
> > 16 0.01805   2      3
> > 17 0.02950   2      4
> > 18 0.05995   2      5
> > 19 0.03425   3      1
> > 20 0.03575   3      1
> > 21 0.03535   3      1
> > 22 0.03380   3      1
> > 23 0.03410   3      1
> > 24 0.01050   3      2
> > 25 0.01735   3      3
> > 26 0.03140   3      4
> > 27 0.06170   3      5
> > 28 0.03430   4      1
> > 29 0.03185   4      1
> > 30 0.03165   4      1
> > 31 0.03380   4      1
> > 32 0.03235   4      1
> > 33 0.01100   4      2
> > 34 0.01825   4      3
> > 35 0.03045   4      4
> > 36 0.06060   4      5
> >
> >         [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@r-project.org mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html>
> > and provide commented, minimal, self-contained, reproducible code.
>
>
>
> --
>
> Bert Gunter
> Genentech Nonclinical Biostatistics
>
> Internal Contact Info:
> Phone: 467-7374
> Website:
>
> http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm
>

	[[alternative HTML version deleted]]