Dear Justo,
The negative covariances in the first column below
sugggest an ARH model will not do too well but to demonstrate the syntax;
If you currently have
5 0 US
0.43
-0.09 0.67
-0.08 0.03 0.72
-0.07 0.01 0.10 0.82
-0.04 -0.02 0.03 0.1 0.94
replace it with
5 0 AR1H .03
.43 .67 .72 .82 .94
Another option is to try the constrained US model (symbollically written)
A
B C
D E F
G H E I
J H H E K
by changing the US line above to
5 0 US !=0000a00ba00bba0
Arthur
>
> Hi ASREMLers!
>
> Fitting an R unstructured (co)variance matrix to my data, I get the estimates:
>
> 0.43
> -0.09 0.67
> -0.08 0.03 0.72
> -0.07 0.01 0.10 0.82
> -0.04 -0.02 0.03 0.1 0.94
>
> Comp/SE for non-diagonal elements are
>
> 9
> -2 15
> -2 1 15
> -2 0 2 17
> -1 0 0 3 17
>
> Based on this, I want to try ARH structure:
>
> a.a
> a.b.r b.b
> a.c.r2 b.c.r c.c
> a.d.r3 b.d.r2 c.d.r d.d
> a.e.r4 b.e.r3 c.d.r2 d.e.r e.e
>
> How could I make that and where find I estimates of (a,b,c,d,e) and r?
>
> Thanks.
>
>
> 3000
> 5 0 US
> §§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§
> § §
> § Justo Lorenzo Bermejo §
> § Institut für Tierzucht und Tierhaltung §
> § Olshausenstr. 40 §
> § 24118 Kiel §
> § Germany &_____ §
> § / ^..^ §
> § T.n. 0431 880 7315 \_____(oo) §
> § Fax 04340 402 761 WW WW §
> § jlorenzo@tierzucht.uni-kiel.de §
> § §
> §§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§
>
> --
> Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD mailto:Arthur.Gilmour@agric.nsw.gov.au
Principal Research Scientist (Biometrics) fax: <61> 2 6391 3899
NSW Agriculture <61> 2 6391 3922
Orange Agricultural Institute telephone work: <61> 2 6391 3815
Forest Rd, ORANGE, 2800, AUSTRALIA home: <61> 2 6362 0046
ASREML is still free by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
or point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/
To join the asreml discussion list, send the message
subscribe
mailto:asreml-request@chiswick.anprod.CSIRO.au
To send messages to the list, mailto:asreml@chiswick.anprod.CSIRO.au
Asreml list archive: http://www.chiswick.anprod.csiro.au/lists/asreml
<> <> <> <> <> <> <>
"Blessed are the peacemakers,
for they shall be called sons of God" Jesus; Matthew 5:9
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml