Multivariate analysis using the new ASREML version
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Multivariate analysis using the new ASREML version



Dear all,
using an old ASREML version, I carried out a multivariate analysis (Code I).
Convergence was achieved after five iterations (please, see results below). 
Now, I would like to use the facility OWN of the new ASREML version. Running
'Code I' on the new version, convergence failed.
Could anyone tell me how to carry out these analyses? Should I modify the
code or should I verify the installation of ASREML?
Thanks in advance,
Justo

-------------------------------------------------------------
Code I

 ident !P
 sex 2 !I
 iniwt
 a
 b !*100
 c !*10000

pedigree
olsq.dat !asuv !MAXIT 50

a b c ~ Trait Tr.sex Tr.iniwt Tr.sex.iniwt,
        !r Tr.ident

1 2 1
79 !S2==1
Trait 0 US !+6
3.062
-4.144872 6.088
1.339425 -2.0818153 0.7566
Tr.ident 2
3 0 US !+6
2.396
-4.075 7.216
1.535254 -2.804 1.09
ident
----------------------------------------------------------------
Results old version

 LogL=-40.2310     S2=  1.0000        225 df   1.0000     2.4987    -3.0832
   4.0561    0.90063    -1.2675    0.42588     3.3127    -5.7797     10.389
   2.2482    -4.0703     1.6122
 Final parameter values                        1.0000     2.4967    -3.0791
   4.0478    0.89888    -1.2639    0.42430     3.3167    -5.7875     10.404
   2.2515    -4.0769     1.6150

 Source                 Model  terms     Gamma     Component  Compnt/StndErr %C
 Residual            US=UnStr     1   2.49668       2.49668       1.85   0 U
 Residual            US=UnStr     1  -3.07908      -3.07908      -1.40   0 U
 Residual            US=UnStr     2   4.04777       4.04777       1.09   0 U
 Residual            US=UnStr     1  0.898877      0.898877       1.08   0 U
 Residual            US=UnStr     2  -1.26389      -1.26389      -0.89   0 U
 Residual            US=UnStr     3  0.424301      0.424301       0.77   0 U
 Tr.ident            US=UnStr     1   3.31670       3.31670       1.54   0 U
 Tr.ident            US=UnStr     1  -5.78748      -5.78748      -1.58   0 U
 Tr.ident            US=UnStr     2   10.4044       10.4044       1.64   0 U
 Tr.ident            US=UnStr     1   2.25153       2.25153       1.60   0 U
 Tr.ident            US=UnStr     2  -4.07690      -4.07690      -1.66   0 U
 Tr.ident            US=UnStr     3   1.61501       1.61501       1.68   0 U
-------------------------------------------------------------------
Results new version

 LogL=-21.3653     S2=  1.0000        225 df
 LogL=-23.8260     S2=  1.0000        225 df
 LogL=-19.2292     S2=  1.0000        225 df
 LogL=-55.5319     S2=  1.0000        225 df
 Logliklihood decreased to     -55.53: Trying again with reduced updates
 LogL=-17.6629     S2=  1.0000        225 df
 LogL=-16.7251     S2=  1.0000        225 df
 LogL=-45.8796     S2=  1.0000        225 df
 Logliklihood decreased to     -45.88: Trying again with reduced updates
 LogL=-18.2237     S2=  1.0000        225 df
 LogL=-17.7640     S2=  1.0000        225 df
 LogL=-4.33345     S2=  1.0000        225 df
 LogL=-9.55048     S2=  1.0000        225 df
 LogL=-16.7605     S2=  1.0000        225 df
 LogL=-81.9788     S2=  1.0000        225 df
 Logliklihood decreased to     -81.98: Trying again with reduced updates
 LogL=-21.7926     S2=  1.0000        225 df
 LogL=-22.1793     S2=  1.0000        225 df
 LogL=-25.9518     S2=  1.0000        225 df
 LogL=0.642927     S2=  1.0000        225 df
 LogL=-8.58751     S2=  1.0000        225 df
 LogL=-19.4729     S2=  1.0000        225 df
 LogL=-45.8992     S2=  1.0000        225 df
 Fault     0 Convergence Failed
--------------------------------------------------------------------

 			                                                                     
  Justo Lorenzo Bermejo                                    
  Institut fr 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                           
  http://www.tierzucht.uni-kiel.de/j_lorenzo.htm           
                                                           


          La vida est hecha solo de momentos

--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml