Dear Dave,
Hi
>
> Q1. Trying to fit a spatial power model to both G and R structures with
> heterogeneous variance. Data comprises repeated measures over 8 days
(unequally
> spaced). Program error is "Failed to parse G structure line". Any suggestions
?
Try removing 8 !I from the day line and
lnggt ~ fac(day) !r tag.fac(day)
1 2 1
anim !S2==1
8 0 POW .9 .05 .22 .38 .45 .65 .47 .50 .50
# assumes data is sorted dat within animal
tag.fac(day) 2
tag 0 AINV
8 fac(day) POW .9 1.04 1.16 1.27 1.24 1.05 1.16 1.00 .91
Check the .res file for encoding of day values.
>
> Q2. For the Cholesky (CHOL) model one needs the !B qualifier. What is the
number
> of bands - is this the rank ?
>
Two bands (default) means matrix looks like
A B 0 0 0 0
B C D 0 0 0
0 D E F 0 0
0 0 F G H 0
0 0 0 H I J
0 0 0 0 J K
Three bands looks like
A B M 0 0 0
B C D N 0 0
M D E F Q 0
0 N F G H P
0 0 Q H I J
0 0 0 P J K
I have implemented but am still testing a reduced rank model
which models S=LL' and L looks like
A 0 0
B C 0
D E F
G H I
J K L
M N O
Arthur
> Thanks, Dave.
> FE 92 J
> tag !P
> sire
> dam
> date 8 !A
> day 8 !I
> lnggt
> anim 154 !I
> jped1.dat
> jggt92d.txt !skip 1 !ASMV 8
> lnggt ~ day !r tag.day
> 1 2 1
> anim !S2==1
> day day POW .9 .05 .22 .38 .45 .65 .47 .50 .50
> tag.day 2
> tag 0 AINV
> day day POW .9 1.04 1.16 1.27 1.24 1.05 1.16 1.00 .91
>
>
> ___________________________________________________________________________
> THE FOLLOWING WARNINGS ATTACH TO THIS COMMUNICATION
> This material is intended for the named recipient only and has been
created
> solely for the purpose and function of Livestock Improvement Corporation
Limited
> and may contain information that is subject to the Privacy Act /
Confidentiality
> / Privilege / Copyright. Any person other than the named recipient is
prohibited
> from retaining / disclosing / copying / distributing / using all or part of
this
> material.
>
>
> --
> 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
<> <> <> <> <> <> <>
"Christ Jesus came into the world to save sinners" I Timothy 1:15.
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml