[R-sig-ME] Running repeated measures in lme4
Emma Stone
Emma.Stone at bristol.ac.uk
Wed Mar 3 13:45:59 CET 2010
Hi Ben,
Thanks for your reply, however I have re ran this and I still get zero
group variance. Also, forgive my ignorance but doesn't this code remove the
repeated measures element of the design? Also for note, I ran my code with
a poisson distribution and I did get variance outputs for the Random Group
level.
Emma
--On 03 March 2010 07:28 -0500 Ben Bolker <bolker at ufl.edu> wrote:
> Emma Stone wrote:
>> Dear all,
>>
>> I wander if you can help.
>>
>> I am running a repeated measures model with lmer.I have 8 sites, at each
>> I conducted an experiment with treatments which have 4 levels. Each
>> treatment was conducted only once per site, but each site has one of
>> each of the treatment levels so it is balanced. Obviously the
>> treatments within sites are the repeated measures component, and I want
>> to look at differences between treatments within sites, I am not really
>> interested in the site group effect but need to incorporate it. So I
>> have set up my model as follows:
>>
>> model1<-lmer~(minutes~treatment+(treatment|site), data = data)
>
> This model fits a site by treatment interaction, which you don't
> have enough information for since you don't have replication. I believe
> you want
>
> model1<-lmer~(minutes~treatment+(1|site), data = data)
>
>
> --
> Ben Bolker
> Associate professor, Biology Dep't, Univ. of Florida
> bolker at ufl.edu / people.biology.ufl.edu/bolker
> GPG key: people.biology.ufl.edu/bolker/benbolker-publickey.asc
----------------------
Emma Stone
Postgraduate Researcher
Bat Ecology and Bioacoustics Lab
& Mammal Research Unit
School of Biological Sciences,
University of Bristol, Woodland Road,
Bristol, BS8 1UG
Email: emma.stone at bristol.ac.uk
More information about the R-sig-mixed-models
mailing list