Displaying 2 results from an estimated 2 matches for "pretot".
Did you mean:
oretot
2002 Mar 31
1
lme degrees of freedoms: SAS and R
...-1 where n is total # of observations, v is the # of
levels for the grouping factor. From SAS df = v -1. Am I wrong about this
or can somebody explain which is correct and why?
Thanks a lot!
Kaiya Liu
-------------------------------------
Here are the codes:
For R:
> formula (SIMS)
Gain ~ Pretot | Class
> data(SIMS)
> fm1SIMS <- lme(Gain ~ Pretot, data = SIMS, random = ~ Pretot | Class,
control = list(msVerbose = TRUE))
> summary (fm1SIMS)
_________________________________
For SAS:
proc mixed data=sims;
class class;
model gain = pretot / solution;
random interc...
1998 Jun 22
0
lme_2.9-2.tar.gz available soon on CRAN
...used in the "SAS Sytems for Mixed Models" book.
> library( lme )
> data( SIMS )
> ### Analysis of the Second International Mathematics Study (SIMS)
> ### described in section 7.2 of "SAS System for Mixed Models"
> unix.time(assign("fm1RSIMS", lme(Gain ~ Pretot, SIMS, ~ Pretot | Class,
+ REML = TRUE)))
[1] 14.25 0.10 15.00 0.00 0.00
> ### Timing was done on a Pentium Pro 200 MHz system running Linux
> ### Your mileage may vary.
> summary(fm1RSIMS) # compare to output 7.4, p. 262, "SAS System ..."
Linear mixed-ef...