Bivariate binomial

From: Macbeth, Michael <Michael.Macbeth_at_DPI.QLD.GOV.AU>
Date: Thu, 31 May 2007 12:54:13 +1000

I need some assistance to estimate the genetic correlations of two
binomial traits. The univariate analysis converged with expected

The data file for the bivariate binomial analysis has been set up as
follows (one field defines the trait and another field has the binary
score of that trait):
Binrait score
1 1
1 0
1 0
2 1
2 1
2 0

The model I have is as follows:
ASREML Linear threshold model
# anim !P
 sire !P
 dam !P
 bintrait # trait 1 or 2
 binscore # binary score 0 or 1 with trait defined in field
bi.ped # pedigree file
bi.dat !ASUV !MAXIT 100 # data file
binscore !BINOMIAL !PROBIT ~ Trait bintrait !r Trait.bintrait.sire
#ran but with singularity
#binsc !BINOMIAL !PROBIT ~ Trait !r Trait.bintrait.sire
#ran unable to invert R or G matrix
1 2 1
0 0 0
Trait 0 US 1 !GP
Trait.bintrait.sire 2
2 0 US !GP !+3
.05 .06
 0 0 AINV

I was unsuccessful in getting convergence. Is the model above correct
to determine genetic correlations of binomial traits?
From what I have seen in the literature a sire model is preferred in
binomial analyses as the animal model does not behave well. Has anyone
had experiences with the stability of bivariate binomial estimates?

Michael Macbeth
The information contained in the above e-mail message or messages
(which includes any attachments) is confidential and may be legally
privileged. It is intended only for the use of the person or entity
to which it is addressed. If you are not the addressee any form of
disclosure, copying, modification, distribution or any action taken
or omitted in reliance on the information is unauthorised. Opinions
contained in the message(s) do not necessarily reflect the opinions
of the Queensland Government and its authorities. If you received
this communication in error, please notify the sender immediately

delete it from your computer system network.
Received on Mon May 31 2007 - 12:54:13 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.)