[S] LMM

Shyh-Forng Guo (guo@calshp.cals.wisc.edu)
Thu, 19 Feb 1998 09:27:52 -0600


Dear S user,

Some days ago Laszlo Tothfalusi asked a question about LME.
I have a similar queation about LME. How to specify a covariance
structure of random effects using LME ?
For example, SAS code would do the job:

*--------------------Known G and R------------------*
| Data represent a genetic example from Henderson | | (1984, p.53). The
data are artificial and consist | | of measurements of two traits on three
animals. | | Both G and R are known. |
*---------------------------------------------------*;

data h;
input trait animal y;
datalines;
1 1 6
1 2 8
1 3 7
2 1 9
2 2 5
2 3 .
;

/* specify g matrix */
data g;
input row col1-col6;
datalines;
1 2 1 1 2 1 1
2 1 2 .5 1 2 .5
3 1 .5 2 1 .5 2
4 2 1 1 3 1.5 1.5
5 1 2 .5 1.5 3 .75
6 1 .5 2 1.5 .75 3
;

proc mixed data=h mmeq mmeqsol;
class trait animal;
model y = trait / noint s p;
random trait*animal / type=un gdata=g g gi s; repeated / type=un sub=animal
r ri;
parms (4) (1) (5) / noiter;
run;

How can I do this type of analysis in Splus ?

Shyh-Forng Guo
-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu. To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message: unsubscribe s-news