Re: Binomial glm in asreml-r

From: <arthur.gilmour_at_DPI.NSW.GOV.AU>
Date: Sat, 22 Dec 2007 13:26:33 +1100

Dear Kevin,

You need to use
family=asreml.binomial() instead of the default family=gaussian()

though I can't see where the binomial bit is mentioned in my copy of the
User Guide.

I hope thats enough. It should let you define r and n
if the data is not binary.

Our annual family update is at

May Jesus Christ be gracious to you in 2008,

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

Personal website:,
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
The ASReml discussion group is at
To join it,
Archives are at

Proposed travel:
UK 28 Jan
USA 11-15 Feb
Turkey 8-14 Apr
Argentina 16-19 September Biometrics 

Kevin Wright <kw.statr_at_GMAIL.COM>
Sent by: ASReml users discussion group <ASREML-L_at_AGRIC.NSW.GOV.AU>
22/12/2007 09:15 AM
Please respond to
ASReml users discussion group <ASREML-L_at_AGRIC.NSW.GOV.AU>


Binomial glm in asreml-r

In asreml-r, how do I specify the model for a binomial glm?

I think SAS uses something like "model x/n = fixedeffects" and I think
lme4 uses "cbind(x,n-x)~fixedeffects".

I tried variations on those syntaxes, but the only thing I could get
to run in asreml is to define a new variable "proportion = x/n" and
asreml(proportion ~ fixedeffects, random, data,

That doesn't use the information about "n" anywhere that I can see, so
I'm wondering if I have the right syntax.

Any pointers would be appreciated.

Kevin Wright
Received on Wed Dec 22 2007 - 13:26:33 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.)