Dear Kim,
The values in the data are used as indexes to pick values from the SET list.
Thus you need to supply 21 numbers in the SET list even though the
first one (in your case) is never actually used.
[ unless you subtract 1 from the data first !-1 !SET 1 2 3 ... ]
If there are not enough values in the SET list as in yourb first form,
zeros are inserted i.e. 21 becomes zero.
ASREML processes the list record by record so never realises that
some positions are never referenced.
Arthur
> X-Sender: kbunter@metz.une.edu.au
> Date: Tue, 25 Jan 2000 11:59:05 +1100
> To: Arthur Gilmour <gilmoua@apollo.agric.nsw.gov.au>
> From: Kim Bunter <kbunter@metz.une.edu.au>
> Subject: Re: third time lucky!
> Mime-Version: 1.0
>
> Hi Arthur,
>
> thanks for your reply - it works of course! However, I am a little bit
> unsure exactly what happens with missing values. Say 1 is missing in my
> data (ie data range is 2-21). I assumed that giving a sequence list (length
> = 20) would be the correct option as the 2 value would be the first index
> value. So my order was:
>
> data coded 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
> I want 1 2 3 4 5 6 7 7 8 8 9 9 10 10 11 11 11 11 11 11 (ie 20 in list)
>
> asreml told me this gives the incorrect no. of levels if I want to describe
> a factor with 11 levels. (ie there were 12). Thus the correct option
> appears to be:
>
> data coded - 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
> 1 1 2 3 4 5 6 7 7 8 8 9 9 10 10 11 11 11 11 11 11 (ie 21
> in list)
>
> Is this option correct (and giving me the grouping that I want)? Do you
> always have to provide the same number of entries in the !SET list as the
> top data value (including codes missing in the data?).
>
> Thanks for your help!
>
> Cheers
>
> Kim
>
> At 10:21 25/01/00 +1100, you wrote:
> >Dear Kim,
> >
> >If the data is coded 1 2 3 4 5
> >and you want 1 2 2 3 3
> >use
> > agegrp !SET 1 2 2 3 3
> >
> > The brackets in the !SUB expression would have confused ASREML because
> it is
> >looking for a list.
> >
> > The !SUB expression works quite differently and is designed for when
> > the data is say 2.3 4.5 99.2 101
> > and you want to recode to 1 2 3 4 respectively.
> > recode !SUB 2.3 4.5 99.2 101
> > would do that.
> >
> > Arthur
> >> X-Authentication-Warning: lamb.chiswick.anprod.csiro.au: petidomo set
> sender
> >to asreml-owner@lamb.chiswick.anprod.csiro.au using -f
> >> Date: Tue, 25 Jan 2000 09:58:03 +1100
> >> To: asreml@chiswick.anprod.csiro.au
> >> From: Kim Bunter <kbunter@metz.une.edu.au>
> >> Subject: third time lucky!
> >> Mime-Version: 1.0
> >>
> >> I am obviously not good with addresses. Third time lucky in sending this -
> >> hope someone can help!
> >>
> >> Cheers
> >>
> >> Kim
> >>
> >> >Hi all,
> >> >
> >> >I'm sure Arthur has already told me how to do this , but would like to
> >> >group factors together according to value - and the syntax I've tried is
> >> >not working. Any ideas?
> >> >
> >> >eg. group 1yolds separately and group 2 with 3yr olds, and 4 with 5
> year olds
> >> >
> >> >agegrp !SUB 1 (2,3) (4,5) #does not work..
> >> >
> >> >I have put commas in, taken them out, but never seem to get the right
> >> >number of groups even if asreml manages not to crash on a line like this.
> >> >(I am using the October 1999 version...).
> >> >
> >> >I would appreciate any help.
> >>
> >> Kim Bunter (M.Rur.Sc)
> >> PhD Student
> >> Animal Genetics and Breeding Unit
> >> University of New England
> >> Armidale, NSW, 2351
> >> AUSTRALIA
> >>
> >> Ph: (02) 6773 3788
> >> Fax: (02) 6773 3266
> >> email: kbunter@metz.une.edu.au
> >> --
> >> 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.
> ><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
> >
> >
> >
> >
> >
> >
> >
> Kim Bunter (M.Rur.Sc)
> PhD Student
> Animal Genetics and Breeding Unit
> University of New England
> Armidale, NSW, 2351
> AUSTRALIA
>
> Ph: (02) 6773 3788
> Fax: (02) 6773 3266
> email: kbunter@metz.une.edu.au
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
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