Dear Luis,
There seems to be several problems.
1 is confusion caused by have terms labelled age_fac and age
Change age_fac to Age_fac
However, you can just have age in the data and use dev(age)
to have a factor version of it.
2 There is a bug in ASREML when checking for interactions involving
Pedigree terms. Hence 'Define structure for spl(age) ' arises
becasue a bug meant it thought this involved and interaction of subject.
3 Do you really want a spline of the genetic variation of age [for each
subject]? In general this would be very computationally demanding.
I'll update the PC version of ASREML on BBSRC within the hour if I can.
The other versions will have to wait a little while.
Arthur
> Date: Tue, 26 May 1998 07:44:58 -0400 (EDT)
> From: cronopio@altavista.net
> Mime-Version: 1.0
> To: asreml@chiswick.anprod.CSIRO.AU
> Content-Transfer-Encoding: 7bit
> Subject: Re: orthog polynomials and splines
>
> Hi Arthur,
>
> Yes, the fixed spline was a typo. I was asking about the structures for
spline because whenever I try to run a program fitting splines I get the error:
>
> Fault 1 Define structure for spl(age)
>
> that is not in the manual.
>
>
> My program is:
>
> Analysis of FA8102, spline (model 1a)
> subject !P
> mother 49
> father 1
> block 8
> age_fac 10 !I
> age
> ht
> di
> c:\arauco\fa8102un.dat !MAKE !REPEAT
> c:\arauco\fa8102un.dat !maxit 30
> ht~mu lin(age) !r spl(age) subject subject.age subject.spl(age) dev(age)
>
> or just
> ht~mu age !r spl(age) subject subject.age subject.spl(age)
>
> with the same problem.
>
> Luis Apiolaza
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> Institute of Veterinary, Animal and Biomedical Sciences
> Massey University
> Palmerston North
> New Zealand
> L.A.Apiolaza@massey.ac.nz
>
> "Have little and you will get,
> have much and you will be confused"
> Lao Tzu
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>
>
>
>
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD email: Arthur.Gilmour@agric.nsw.gov.au
Senior 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 currently free by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
Point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/
in the IACR-Rothamsted information system http://www.res.bbsrc.ac.uk/
To join the asreml discussion list, send the message
subscribe
to asreml-request@chiswick.anprod.CSIRO.au
The address for messages to the list is asreml@chiswick.anprod.CSIRO.au
<> <> <> <> <> <> <>
"Seek first the kingdom of God and His righteousness" Jesus; Matthew 6: 33
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>