Re: Estimation of additive effect?

Dear Khandker

You wrote
This is a follow up mail of my previous email entitled "How to get
additive and dominance effect value".

Thanks Arthur for your reply and by using the commands you mentioned I
can get allele substitution effect as well as dominance effect.
(resp ~ mu SNP at(SNP,1) !r animal),
This command gives me both measurement at one single run.

However, still I can't figure out the additive effect which is
estimated as the difference between the two homozygous means divided
by two (Falconer, D. S., and T. F. C. Mackay. 1996. Introduction to
Quantitative Genetics. 4th ed. Longman Scientific and Technical, New
York, NY.). Here, I am using SNP with three genotypes
(0=homozygous,1=heterozygous,2=homozygous), so it will be the estimate
of difference between 0 and 2 divided by 2.

1) I have not reverenced Falconer, so I may have missed something.

2)
Assuming SNP is a variable with values of 0,1 or 2
and is declared as a simple variate in ASReml (i.e. no factor
qualification appended)

the the design matrix for the fixed effects will be

mu SNP at(SNP,1)

1 0 0
1 1 1
1 2 0

and let us say the 3 fitted effects are c (for mu), a for SNP and d for
at(SNP,1)

So the SNP=0 class BLUE is c
      SNP=1 class, BLUE is c+a+d
   SNP=2 class, BLUE is c+2a

So additive effect is (c+2a -c)/2 = a

So, a is the additive effect and the average allele substitution effect.

Unless I have missed something.

May Jesus Christ continue to be gracious to you,

Arthur Gilmour, His servant .
 
Mixed model regression mapping for QTL detection in experimental crosses.
Computational Statistics and Data Analysis 51:3749-3764 at
http://dx.doi.org/10.1016/j.csda.2006.12.031

Profile:http://www.dpi.nsw.gov.au/research/staff/arthur-gilmour
      http://www.dpi.nsw.gov.au/reader/17263
Personal website: http://www.cargovale.com.au/

Skype: arthur.gilmour
mailto:Arthur.Gilmour_at_dpi.nsw.gov.au, arthur_at_cargovale.com.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

ASREML 2 is now available from http://www.VSNi.co.uk/products/asreml
Register on the ASReml forum at http://www.vsni.co.uk/forum
ASReml-L Archives are at
https://listserv.dpi.nsw.gov.au/cgi-bin/wa.exe
Cookbook: http://uncronopio.org/ASReml

The ASReml3 alpha download registration page is:
http://www.austatgen.org/asreml3alpha.php
There are asreml3alpha forums at: http://www.austatgen.org/asreml3forum

Proposed travel:
Argentina 9-19 September ASReml WS + Biometrics
www.gab2008.com.ar/actividades.php
Coffs Harbour WS 14-17 October (Forestry)
><><><><><><><><><><><><><><><><><><><><><><><>
Received on Sat Aug 28 2008 - 10:27:02 EST