[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