[R] Summary: proc mixed vs. lme
f0z6305 at labs.tamu.edu
Wed Oct 9 22:25:46 CEST 2002
Sorry to bother you.
I want to know if there are some Logistical Regression
functions in R?
----- Original Message -----
<dominik.grathwohl at rdls.nestle.com>
To: <r-help at stat.math.ethz.ch>
Sent: Wednesday, October 09, 2002 1:19 PM
Subject: [R] Summary: proc mixed vs. lme
> Summary: proc mixed vs. lme
> The objective of this summary is to help people
> to get more familiar with the specification of
> random effects with proc mixed or lme.
> Very useful are the examples of Ramon Littell's book:
> "SAS System for Mixed Models (1996)"
> The same data set's are kindly made available
> by Douglas Bates in the library(SASmixed).
> In the help file are examples of the lme statements
> equivalent to the proc mixed ones.
> To explain the different estimates,
> Hein and Brian suppose to check whether both
> analyses with SAS and R uses ML estimates
> or REML estimates. However, this was not the problem,
> the default in proc mixed and lme is already REML.
> Douglas advise me to use his option:
> options( contrasts = c(unordered = "contr.SAS",
> ordered = contr.poly"))
> However, I already used this option,
> because I copied the code from the SASmixed help file.
> Peter gave me the first hint
> and this solves the problem:
> To change the model formula in lme,
> from: strength ~ Program * Time
> to: strength ~ factor(Program) * factor(Time)
> Now the option statement grasp!
> For people like me who try to get familiar
> with the specification of random effects
> would it helpful if the help file of SASmixed
> would be updated or the variables time and program
> would be already introduced as factors
> in the Weights data set.
> Thank you all for the useful advises,
> r-help mailing list -- Read
> Send "info", "help", or "[un]subscribe"
> (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
More information about the R-help