[R-sig-ME] Mixed linear model with nested and interaction term

[Lin, Heng-An]

Hi all,

I am analyzing my data with following model,

model1 <- lmer(Yield~Treatment+(1|Location)+(1|Location:Treatment)+(1|Location:Block), data=A)

in here, I want to set an random interaction term (Location*treatment) and an random nested term (block nested within location).

But I couldn't get similar ANOVA results when I compare the output with SAS porc mixed output.

So, I think i might make some mistake in the model in R...

Can anyone give me some suggestion?

Thanks in advance!

Heng-An

	[[alternative HTML version deleted]]