Re: Output G and R matrices?

From: <brian.cullis_at_DPI.NSW.GOV.AU>
Date: Mon, 15 Jun 2009 16:51:18 +1000

Dear jason
ok, are you using ASReml-R or stand alone?
I have functions written in R for generating AR1 structures, which are
easy to use and I would be happy to provide. I caution that you may be
challenged in various ways
1. size of the problem
2. if you are fitting AR1xAR1 models then presumably you have estimated
these and therefore would need to incorporate such in to a fully
efficient one-stage modelling approach
3. the animal model variance matrix is returned from asreml.Ainverse() as
a sparse data frame (see manual for details), but again if the problem is
large you would then have to invert this and include this somehow with the
R package?

I have colleagues who are extending the WGAIM QTL approach (TAG Verbyla et
al, 2007?) which may be another approach to consider

warm regards

Brian Cullis
Research Leader, Biometrics &
Senior Principal Research Scientist
NSW Department of Primary Industries
Wagga Wagga Agricultural Institute

Visiting Professorial Fellow
School of Mathematics and Applied Statistics
Faculty of Informatics
University of Wollongong

Faculty of Agriculture, Food & Natural Resources
The University of Sydney

Adjunct Professor
School of Computing and Mathematics
Charles Sturt University

Phone: 61 2 6938 1855
Fax: 61 2 6938 1809
Mobile: 0439 448 591

jason <asremlforum_at_VSNI.CO.UK>
Sent by: ASReml users discussion group <ASREML-L_at_DPI.NSW.GOV.AU>
15/06/2009 04:24 PM
Please respond to


Re: Output G and R matrices?

Thank, Brian. I am trying some R packages for Bayesian association
analysis for SNPs which require design matrices as well as the full R
matrix as input. Since I am using an animal model with AR1xAR1 structure,
the R matrix seems a bit complicated. It would be helpful if ASReml can
provide an option to output the R matrix (maybe G matrix as well).


-------------------- m2f --------------------

Sent using Mail2Forum (

Read this topic online here:

-------------------- m2f --------------------
Received on Tue Jun 15 2009 - 16:51:18 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.)