F-tests for lme results

Regula Suter (suter@inw.agrl.ethz.ch)
Wed, 14 Jan 1998 14:27:57 +0100


I posted this message a week ago and it seems it didn't appear in the
newsgroup. So I decided to send it again:

I am using S-Plus, Version 3.4 Release 1 for Sun SPARC, SunOS 5.3

I am analyzing an unbalanced repeated measures design using "lme". I got the
fixed effects estimates (different groups and regressors) listed below.
Then I tried to get Type III Sums of
Squares, Mean Squares, the corresponding F-Values and Probabilities as one
gets them by using "proc mixed" in SAS. I looked through the S-Archive and
didn't find much help.

What do I have to do in order to get all those values? Can anyone out there
help?

Thanks a lot in advance.
Regula Suter

> summary(lme(fixed = TS ~ C(race, sum) + C(group,sum) + BCS + LG + SKL,
random= ~ -1, cluster= ~ LaufNummer, re.structure = "identity", data =
co.na, est.method = "RML"))

... snip...

Fixed Effects Estimate(s):
Value Approx. Std.Error z ratio(C)
(Intercept) 48.57359531 6.624871631 7.3320055
C(race, sum)1 1.01175749 0.560594675 1.8047933
C(race, sum)2 0.23621308 0.269235597 0.8773471
C(group, sum)1 -2.55432600 0.305444022 -8.3626649
C(group, sum)2 0.93351974 0.288713437 3.2333782
C(group, sum)3 -0.92666219 0.290456390 -3.1903660
BCS -3.43810275 0.439285416 -7.8265807
LG 0.03051065 0.005684205 5.3676197
SKL -0.21553492 0.049827914 -4.3255858

...snip...

_____________________________________

Regula Suter
Swiss Federal Institute of Technology
Institute of Animal Science
Animal Breeding
Research-Station Chamau
CH-6331 Huenenberg
phone: +41 41 780 10 47
fax: +41 41 780 43 83
e-mail: suter@inw.agrl.ethz.ch
http://www.tz.inw.agrl.ethz.ch/~suter
_____________________________________