Hi ASREMLers!
Any suggestions about these two problems?
1. A question about the structure of the residual (R).
I have individual data of feed intake of approx. 3000 pigs at the 25 days
3,4,5,6,7, 17-21, 31-35, 45-49 and 59-63. I would like to try heterogeneous
AR and heterogeneous ARMA for these data, so that R(ARH)=
a.a a.b.r a.c.r2 a.d.r3 a.e.r4 (~0) ...
b.b b.c.r b.d.r2 b.e.r3
c.c c.d.r c.e.r2
(sim) d.d d.e.r
e.e
f.f f.g.r f.h.r2 ...
g.g g.h.r ...
... ...
How could I do that? Perhaps using 'POW'? Gamma results using DIAG were:
1.18 0.91 1.00 0.98 1.09
1.21 1.22 1.20 1.22 1.35
1.44 1.45 1.46 1.31 1.35
1.66 1.72 1.67 1.52 1.63
1.96 1.95 1.80 1.64 1.76
2. Is the G matrix correct as follows or should I modify it according to the
data structure?
fi~ mu !r ident ident.x ident.x2 !f mv
0 0 1
ident 2
3 0 US !+6
0.03681
0.001 0.08332
0.001 0.001 0.08332
ident 0 AINV
§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§
§ §
§ 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 §
§ §
§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§§
La vida está hecha solo de momentos
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml