Re: GIV matrix

From: <arthur.gilmour_at_DPI.NSW.GOV.AU>
Date: Thu, 7 Jun 2007 15:57:23 +1000

Dear Vivi,

There are no obvious syntactical problems with your code
except it should have an itial value for the genetic variance, but
it looks like you are fitting a HUGE model.

BLUE !WT wt ~ mu !r gen yea.occ.gen

1 1 1
0 0 0 !s2== 0.4430 # No point using so many decimal places; 36672875

gen 1
gen 0 GIV 0.1

Should be exactly equivalent to

BLUE !WT wt ~ mu !r giv(gen,1) yea.occ.gen

1 1 1
0 0 0 !s2== 0.4430 # No point using so many decimal places; 36672875

I say HUGE because yea.occ.gen has 25 x 208 x 686 = 3567200 levels
so presumably you have more than that records of data, depending on the

Furthur, if you want a GIV structure on the main effects, you presumably
want it on the interaction as well.
This does require seting up a G structure line [yea.occ.giv(gen,1) does
not actually apply GIV1 ton the gen dimension]
yea.occ.gen 2
5200 0 I
gen 0 GIV1 .1 !GP

If there are many empty cells in yea.occ, then recoding the interaction
as a single factor
will reduce the dimensionality and maybe allow the model to run.

I trust this helps

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

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,

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 Jun 07 2007 - 15:57:23 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.)