Thanks Luc, I used the overlay feature yesterday and it reduced to a more
comfortable 4 hours per iteration. I have one question re omission of the !MGS
flag on the pedigree file. If I do that, how will the program then know how to
form the relationship matrix (inverse) for the sire factor on sire/mgs rules
rather than animal model rules?
Dave
"Janss, L.L.G." <l.l.g.janss@id.dlo.nl> on 14/08/2000 19:48:11
To: asreml@chiswick.anprod.csiro.au
cc: (bcc: Dave Johnson/ham/LIC)
Subject: RE: sire-mgs model
Dave,
it is much more efficient to fit a sire-mgs model using the 'overlay'
option, rather than using the ASREML build-in facility (which transforms it
into an animal model when I understood correctly). To overlay mgs with sire
your model would become:
prot% ~ mu season period sea.per jersey het cowage region,
!r sire 0.207 -mgs and(mgs,0.5) 0.05 key 1.003 locn 2.066
and omit the !MGS on the pedigree file name. That should work similarly fast
as the plain sire model.
Luc.
Dr. Luc Janss, Scientist Quantitative Genetics
ID-Lelystad Institute for Animal Science and Health
P.O. Box 65 Tel. +31.320.23.82.59
8200 AB Lelystad, The Netherlands Fax. +31.320.23.80.50
L.L.G.Janss@id.wag-ur.nl !! NEW EMAIL
> -----Original Message-----
> From: Dave Johnson [SMTP:djohnson@lic.co.nz]
> Sent: Friday, August 11, 2000 3:04 AM
> To: asreml@chiswick.anprod.csiro.au
> Subject:
>
> Thanks to earlier responses I have the sire model working and it took
> about 10
> minutes per iteration. I then added the maternal grandsire (same factor as
> sire,
> ~10800 levels) without fitting the covariance and it seems to be
> converging but
> taking about 1 day per iteration. I also had to increase memory at least
> 4-fold
> (-s5 to -s7). I'm puzzled at the great increase in memory and time - am I
> expecting too much of this model? The mgs factor has some missings and
> I'm not
> sure what the program does with those records? Thanks.
>
> Region comparison of concentrations
> key 106031 !I
> region 6
> jersey !/16
> het
> sire !P
> mgs !M0 !P
> locn 14105
> season 11 !I
> period 3
> cowage 4 !I
> dim
> logdim
> pro_fat
> fat%
> prot%
> ped.txt !MGS
> regions.txt
> prot% ~ mu season period sea.per jersey het cowage region,
> !r sire 0.207 mgs 0.05 key 1.003 locn 2.066
>
>
>
> __________________________________________________________________________
> _
> 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
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml