> Code I have used thus far without being able to replicate the
> data includes:
>
> Fm<-lmer(Score~(1|Line%in%Set)+Set+(1|Block))
> (I figured out how to get a p-value, but it didn't yield the
> same results as those obtained in SAS)
%in% doesn't generally mean 'nested in' in R. It is a set membership test and
will return TRUE for those labels in Line that are also in Set and FALSE
otherwise.
Did you mean Score~(1|Set/Line)...?
If you did, bear in mind that , combined with the fixed Set term, (1|Set/Line)
implies a random Set grouping effect as well as a fixed effect - not sure that
makes sense in your circumstance unless Set is a continuous predictor. May be
safer to define SetLine<-interaction(Set, Line) and do
Score~Set + (1|SetLine) + (1|Block)
S Ellison
*******************************************************************
This email and any attachments are confidential. Any use...{{dropped:8}}
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
