Re: OWN structures

Dear Ian,

Best wishes for Christmas.

re:
There seems to be a max size to the order of OWN structures. Does
anybody know what it is?

With the structure specified in sparse format, is it feasible to create
a diagonal OWN structure of order approx 10000?

As you anticipate, a dense form of the matrix would be too large
but a diagonal form of this size, indeed 100000 would be reasonable.

While it does not appear in the User Guide, the program does appear to
have a facility
of handling a DIAGonal OWN matrix.

The key is to write a 2 as a fourth field (column 20) of the header line
in the .gdg file
and then just to provide the diagonal elements rather than the whole
matrices.

Let me know if it works.
There are other options in this fourth field but I'm unsure how well they
have been tested.

 Check out http://www.cargovale.com.au/christmas06
for my personal greetings and update (I don't know what year it is!).

Thanking God for His mercy through Jesus Christ.

Arthur Gilmour PhD
mailto:Arthur.Gilmour_at_dpi.nsw.gov.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

http://www.dpi.nsw.gov.au/reader/17263
http://www.asreml.co.uk/
The ASReml discussion group is at ASREML-L_at_dpi.nsw.gov.au
To join it, mailto:arthur.gilmour_at_dpi.nsw.gov.au
Cookbook: http://uncronopio.org/ASReml

Proposed travel:
REML Mixed model RR 19Jan
Wool meets meat 22-23 Feb Orange
IBC Quebec 16-21 July 2006
WCGALP8 13-18 Aug Brazil
      <><><><><><><><><><><><><><><><><><><><><><><><><>

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 Tue Dec 21 2005 - 12:54:13 EST