Re: Binomial glm in asreml-r

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
http://www.cargovale.com.au/Christmas07r.htm

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
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
Archives are at
https://gatekeeper.dpi.nsw.gov.au/Listserv/archives/asreml-l.html
Cookbook: http://uncronopio.org/ASReml

Proposed travel:
UK 28 Jan
USA 11-15 Feb
Turkey 8-14 Apr www.icapss.com
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>

To
ASREML-L_at_AGRIC.NSW.GOV.AU
cc

Subject
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
then
asreml(proportion ~ fixedeffects, random, data,
family=asreml.binomial(link="logit"))

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