[R-sig-ME] RE : Interactions between fixed and random effects
Julien Beguin
julien.beguin.1 at ulaval.ca
Mon Apr 5 13:23:13 CEST 2010
If you have 22 populations, and assuming that each indivudual in each population is randomly distributed within you water treatment (you did not specify your experimenta design... so I guess), you might have enough power to test the interaction between plant population and water treatment (21 df) in a similar way than in a classic ANOVA. So why to specify plant population as a random variable, rather than a fixed factor? are you (not) intersted to know which plant population respond differently regarding watering regimes?
proposition:
> lm(DryMass~Water*Population)
just a thought...
Julien Beguin
________________________________________
De : r-sig-mixed-models-bounces at r-project.org [r-sig-mixed-models-bounces at r-project.org] de la part de Christopher Eckert [chris.eckert at queensu.ca]
Date d'envoi : 4 avril 2010 23:39
À : r-sig-mixed-models at r-project.org
Objet : [R-sig-ME] Interactions between fixed and random effects
I apologize if the answer to this query is somewhere totally obvious, but i couldn't find it.
I am trying to analyze an experiment where a set of 22 populations of a dune plant species (populations were randomly chosen from across the species' geographic range) were grown in a glasshouse under two different watering regimes (Water = Control vs. Drought). DryMass is the response variable. There was about 20 individuals from each population grown in each Water treatment.
Population is a random effect, but I would like to test for an interaction between Population and Water -to ask the question: do different populations respond differently to drought?
>From what I can gather this is analogous the random intercepts and slopes model discussed in Zuur et al and elsewhere, except that I am examining a categorical predictor (Water) rather than a continuous predictor.
Am I right in thinking that the basic syntax using lme is:
lme(DryMass~Water,random=~Water|Population)
and the syntax using lmer is:
lmer(DryMass~Water+(Water|Population))
Thanks very much for any clarification you can send my way.
Chris Eckert
Department of Biology
Queen's University
Kingston, Ontario K7L 3N6 Canada
chris.eckert at queensu.ca
_______________________________________________
R-sig-mixed-models at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
More information about the R-sig-mixed-models
mailing list