sas to R
Thank all of you. Ben and Kevin, you were right! The main reason of my problem that I experienced was not specifying factors (i.e., trial and block)... Once I factorized them, I was able to run the models with both lme and lmer in the SAME R session. I think it was possible since I did not use functions in lmer such as lmList. Also, the results, especially df and F-stats, were same or very close to results from SAS. I greatly appreciate all of your kind/patient suggestions. Thank you again!!! Steve
On Tue, Jun 26, 2012 at 2:10 AM, Ben Bolker <bbolker at gmail.com> wrote:
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]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models