Re: R-structure
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: R-structure




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