Skip to content
Prev 156841 / 398503 Next

lme problems

Hi Tommaso,
To begin with, why try to compare things that are obviously quite different?
Surely you can see that the error structure of the two models are different?

aov(Mean1~treatment*layingday+Error(male.pair/treatment/layingday))
lme(Mean1 ~ treatment*layingday, random = ~1|male.pair)

If you want to compare them then at least make them "equal," otherwise what
is the point? (And one might ask, What is the point, anyway?)

## This would be a reasonable comparison
aov(Mean1~treatment*layingday+Error(male.pair))
lme(Mean1 ~ treatment*layingday, random = ~1|male.pair)

Regards, Mark.
Tommaso Pizzari wrote: