Dear Tore,
re:
Regarding computation of "real" breeding values when using an individual
(tree) model with genetic groups
Model: y ~ block !r genid
In the sln file, I identify
Block effects: 0, b2, b3, b4, b5 (blue for five blocks, a singularity
due to groups)
Group effects: g1, g2, g3, g4 (four genetic groups)
Individual breeding values (BVs): One for each tree in the data file (or
in the pedigree file)
ARG: I like to check this out on a simple example but my interpretation is
that the overall mean is included in the group effects and therefore in
the 'BV's. i.e. these values are relative to Block_1
So a 'least squares mean' overall would be the average of the block and
group effects.
You would get relative BV's centred on average of Groups by subtracting
the
average of the 4 group means.
Problem: How should I compare and/or combine the mean of (1) BVs, (2)
group effects, and (3) block effects?
So, you could predict group effects with
predict genid 1:4 !ave block
You can predict block means with
predict block !ave genid { 4*1 ???*0}/4
I assume that the group effects are included in the BVs, and that the
group effects represent the group means in some way.
YES: for BLOCK_1
The averages of BVs and group effects, respectively, do not agree very
well for me. Also, adjusting either of them with tha mean of the block
effects (zero included) makes it still worse.
Only in a very simple model would the average of tree effects be same as
group means. They are regressed to that mean.
Question: How should I best combine these effects to get values that are
adequate for computation of relative breeding values?
The breeding values should be directly comparable as they are
(differencing will cancel any offset added).
Presumably you are not wanting to subtract out the 'group' effects on the
assumption these are genetic.
Trust this helps.
Grace, mercy and peace to you from God our Saviour Jesus Christ
Arthur Gilmour PhD (see http://www.cargovale.com.au/)
mailto:Arthur.Gilmour_at_dpi.nsw.gov.au
Principal Research Scientist (Biometrics)
NSW Department of Primary Industries
Orange Agricultural Institute, Forest Rd, ORANGE, 2800, AUSTRALIA
fax: 02 6391 3899; 02 6391 3922 Australia +61
telephone work: 02 6391 3815; home: 02 6364 3288; mobile: 0438 251 426
http://www.dpi.nsw.gov.au/reader/17263
ASREML 2 is now available from http://www.VSNi.co.uk/products/asreml
The ASReml discussion group is at ASREML-L_at_dpi.nsw.gov.au
To join it, mailto:arthur.gilmour_at_dpi.nsw.gov.au
Cookbook: http://uncronopio.org/ASReml
Proposed travel:
Aust Genstat Conf 5-8 Dec Victor Harbour
Leave 18 Dec to 5 Jan
ASReml workshop (Forestry) 16-18 April Tasmania
ASReml workshop (Animal Breeding) May South Africa
AAABG 24-26 September Armidale
<><><><><><><><><><><><><><><><><><><><><><><><><>
This message is intended for the addressee named and may contain
confidential information. If you are not the intended recipient or
received it in error, please delete the message and notify sender. Views
expressed are those of the individual sender and are not necessarily the
views of their organisation.
Received on Wed Sep 22 2006 - 16:47:42 EST
This webpage is part of the ASReml-l discussion list archives 2004-2010. More information on ASReml can be found at the VSN website. This discussion list is now deprecated - please use the VSN forum for discussion on ASReml. (These online archives were generated using the hypermail package.)