OWN instead of USxI
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

OWN instead of USxI



Dear ASREMLers,
I am trying to replace next code (please, see the complete code below) 

1 2
79 !S2==1
Trait 0 US !+6
3.062
-4.144872 6.088
1.339425 -2.0818153 0.7566

and use the facility OWN in order to model the product US x I (or I x US,
but not US alone and multiply afterwards). 

The code:

1 1
237 0 OWN6 3.062 -4.144872 6.088 1.339425 -2.0818153 0.7566 

produces the fault '2 Sorting data into field order' (.asr file below). How
could I circumvent this?   

Thanks in advance,
Justo

-------------------------------------------------------------
Code to be replaced

 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


----------------------------------------------------------------
New code
 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 1 1
237 0 OWN6 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.777 1.09
ident

----------------------------------------------------------------
.asr file of new code

  ASReml [13 Jun 2001]  Two step methods
 10 Aug 2001 17:00:24.829   8.00 Mbyte  Unix   olsq1
 Your ASReml license expires 31 Mar 2001
 * ASREML 2001 *********************************************
 * Jan Scale of spline and pol design matrices has changed.*
 *        This will change the components but not the fit. *
 * Feb predict syntax has changed.                         *
 ***************************************************** ARG *
  sex 2 !I
  b !*100
  c !*10000
 A-inverse retrieved from ainverse.bin
 PEDIGREE [pedigree ] has      128 identities,     353 Non zero elements
 QUALIFIERS: !ASUV !MAXIT 50
  Reading olsq.dat  FREE FORMAT skipping 0 lines
Multivariate analysis of a              b              c
 Using       79 records of      79 read
  Model term      Size Type    COL   Minimum    Mean      Maximum   #zero #miss
   1 ident         128 Direct    1   49.00     88.16      128.0         0     0
   2 sex             2 Factor    2      1     1.5063          2         0     0
   3 iniwt           1 Covariat  3 -10.00     0.2722      12.00         6     0
   4 a               1 Variate   4 -7.304     -1.143      6.069         0     0
   5 b               1 Variate   5 -6.078      3.423      13.99         0     0
   6 c               1 Variate   6 -5.046    -0.6556      2.624         0     0
   7 Trait           3 Traits/Variat
   8 Tr.sex          6 Interaction  7 Trait   :    3    2 sex            :    2
   9 Tr.iniwt        3 Interaction  7 Trait   :    3    3 iniwt          :    1
  10 sex.iniwt       2 Interaction  2 sex     :    2    3 iniwt          :    1
  11 Tr.sex.iniwt    6 Interaction  7 Trait   :    3   10 sex.iniwt      :    2
  12 Tr.ident      384 Interaction  7 Trait   :    3    1 ident          :  128
    237  OWNxx          3.0620   -4.1449    6.0880    1.3394   -2.0818    0.756\
6
 Warning: Spatial mapping information for side 1 of order 237
          ranges from       0.0 to       1.0
 Fault     2 Sorting data into field order
  Last line read was:  237 0 OWN6 3.062 -4.144872 6.088 1.339425 -2.0818153 0.7\
566 0 0 0 0 0
 olsq1 olsq.dat
  Model specification:  TERM LEVELS GAMMAS
 Tr.sex.iniwt                          6     0.000
 Tr.iniwt                              3     0.000
 Tr.sex                                6     0.000
 Trait                                 3     0.000
 Tr.ident                            384     0.100
SECTIONS     237       6       1
    TYPE       0       0       0       0       0
  STRUCT     237       6       7       7       0       1       0
   12 factors defined [max 500].
   12 variance parameters [max 600].   3 special structures
 Final parameter values                       0.10000    -.10000E-36 3.0620
  -4.1449     6.0880     1.3394    -2.0818    0.75660
  Last line read was:  237 0 OWN6 3.062 -4.144872 6.088 1.339425 -2.0818153 0.7\
566 0 0 0 0 0
    12     2     79   237   1000
 Finished: 10 Aug 2001 17:00:24.966   Sorting data into field order            \




 			                                                                     
  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