[R-sig-ME] sas to R
Ben Bolker
bbolker at gmail.com
Tue Jun 26 09:10:21 CEST 2012
Steve Hong <emptican at ...> writes:
>
> Thank you, all.
>
> I restarted R and checked package list using 'search()'. And then,
> loaded 'lme4' and ran the model. However, the results are same...
>
> Below is what I did.
>
> Thank you much again!!!
>
> Steve
As Kevin Wright said:
Make sure that trial, block, trt are factors,
[snip snip snip]
> > bcwrear1 <- read.table("C:/bcwrear1.txt", header=T)
> > df <- bcwrear1
> > df=bcwrear1
> > df=transform(df, y=day10, trt=turf)
df <- transform (trial=factor(trial), block=factor(block), trt=factor(trt))
> > fm.lmer <- lmer(log10(y) ~ trt + (1|trial/block/trt), data=df,
na.action=na.omit)
[snip snip snip]
More information about the R-sig-mixed-models
mailing list