danselechnik at gmail.com
Sat Jul 15 03:26:42 CEST 2017
My name is Dan and I'm a PhD student in Australia. I was hoping that I
could ask you for some help with using lme4. I have a dataset in which I
have PC1 as a response variable. Population, treatment, RBC, and
population*treatment are my explanatory variables. ID is a random factor.
(I have attached the CSV file here)...
I am trying to run a power analysis, and first to fit my data using lmer.
First I read my data into R:
Then I attempt to fit:
fm1=lmer(pc1$PC1 ~ pc1$RBC + pc1$Population + pc1$Treatment +
pc1$Population*pc1$Treatment + (1|pc1$ID), data=pc1, REML=FALSE)
However, this fails, returning the message:
Error: number of levels of each grouping factor must be < number of
My number of populations and treatments is much less than my number of
observations, so I am not sure why I am getting this error... Also, when I
run this, it works fine:
fm1=lm(pc1$PC1 ~ pc1$RBC + pc1$Population + pc1$Treatment +
pc1$Population*pc1$Treatment + (1|pc1$ID), data=pc1)
I was hoping I could ask for your assistance in figuring out what may be
the problem. Thank you very much.
More information about the R-sig-mixed-models