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 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 §
§ 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