>
>
> Dear ASREMLers
>
> I have a VERY simple data set, with a univariate model and 40 provenances,
> with unequal number of families in each provenance. (see prog. below)
>
> Some of my family within provenance variance is coming out as Bound (B). I
> therefore want to change the !GP that is assumed to !US for each of the
> family within variance components. I have tried !US after the factor in
> the model, but this doesn't work for all the variance components.
>
For !US read !GU
> I also tried specifying starting values and !US as a G structure, but I
> couldn't seem to get it to work.
> PROV/FAM 2
> PROV 0 US !+40 !GU
> (40 starting values for 40 different family within provenance effects)
> FAM 0 0
>
> is this right?
>
> How can I specify !US for a nested structure for a univariate analysis?
>
> Any help will be really appreciated!!!
> cheers
> Heidi
>
Dear Heidi,
Rather than use PROV/FAM in the model, use PROV.FAM
and then the G structure
PROV.FAM 2
PROV 0 DIAG !GU
(40 starting values for 40 different family within provenance effects)
FAM 0 0
should do the trick.
Best wishes
Arthur
PROV/FAM is expanded into 40 terms like PROV_1.FAM
so the syntax does not allow specification of indivivdual starting values
and !GU qualifications.
> ******************************************************************
> OP PROV PROGENY DATA
> TREE
> MUM 100
> DAD 1
> REP 5
> IBLOCK 95
> FAM 444
> PROV 40 !A
> HT3
> DATA.dat
> HT3 ~ mu REP PROV !r REP/IBLOCK PROV/FAM
> 0
>
>
>
> Heidi Dungey
> Forest Research
> Private Bag 3020
> Rotorua
> New Zealand
>
> Phone 64-7-343 5629
> Fax 64-7-343 5330
>
> --
> 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
Cargo: <61> 2 6364 3288
Until July 2001, I will be mainly building a house at Cargo.
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
<> <> <> <> <> <> <>
"Blessed are the peacemakers,
for they shall be called sons of God" Jesus; Matthew 5:9
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml