Dear TOre,
First let me apologise for not replying on Friday.
I simply ran out of time.
Both of your proposals for fitting PARENT effects
in nthe diallel are possible in ASREML. In the first case, the elements
of the job would be
My DIALLEL
father 50
mother 50 # coded identically with father
response
diallel.dat
response ~ mu !r father and(mother)
The and() function lays the 'mother' design matrix over the
'father' matrix. Any 'selfed' plants would have a 2.
All others will have a one in the columns relating to father and mother.
If you wanted to add maternal effects or family effects
just add the terms.
The alternative where you actually supply the design matrix
is also possible. The elements of that job would be
My DIALLEL
parent 50 !F
response
diallel.des
response ~ mu !r parent
In this case, the !F qualifier specifies that the factor 'parent'
is already in the form of a design matrix (50 columns in the
data file).
Obviously, the first method is much easier to set up
but both are possible.
Arthur
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD mailto:Arthur.Gilmour@agric.nsw.gov.au
Senior Research Scientist (Biometrics) fax: <61> 2 6391 3899
NSW Agriculture <61> 2 6391 3922
Orange Agricultural Institute telephone work: <61> 2 6391 3815
Forest Rd, ORANGE, 2800, AUSTRALIA home: <61> 2 6362 0046
ASREML is still free by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
or point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/
To join the asreml discussion list, send the message
subscribe
mailto:asreml-request@chiswick.anprod.CSIRO.au
To send messages to the list, mailto:asreml@chiswick.anprod.CSIRO.au
Asreml list archive: http://www.chiswick.anprod.csiro.au/lists/asreml
<> <> <> <> <> <> <>
"Christ Jesus came into the world to save sinners" I Timothy 1:15.
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml