Skip to content

lme problems

2 messages · Tommaso Pizzari, Mark Difford

#
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: