Re: Bivariate binomial

Dear Michgael,

I'm not convinced the model will do the right thing.

Assuming the data is in N sets of 2 traits, try

1 2 1
N 0 0 # N is number of records
2 0 US 1 .5 1 !GFPF !S2==1
Trait.bintrait.sire 2
2 0 US !GP
.06
.05 .06
 0 0 AINV

To check that it is probably working, fix the two covariances to ZERO
and compare the result with separate univariate analyses.

May you know Jesus Christ and His blessing in 2007,

Arthur Gilmour, His servant .

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

Profile: http://www.dpi.nsw.gov.au/reader/17263
Personal website: http://www.cargovale.com.au/

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
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:
Leave early September
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 Mon May 31 2007 - 15:03:08 EST