Re: How to write this code when age and it's quadratic fitted as fixed covariates in model?

Dear Dr. Victor：
thanks for your detailed reply.
 I found that we can apply pow () function to replace of forming the data
field. Then the model was:
Wt ~ mu Age pow(Age, 2) !r AnimalID.

I found that there are differences between Wt ~ mu Age Age2 !r AnimalID
model and Wt ~ mu pol(Age, -2) !r AnimalID model for the solution of
regression coefficients. What's the reason?

luan sheng
----------------------------------------
  Yellow Sea Fisheries Research Institute
  Chinese Academy of Fisheries Sciences
  Nanjing Road 106
  Qingdao 266071
  China
  luansheng_at_gmail.com
  luansheng_at_ysfri.ac.cn
----------------------------------------

Sent from Qingdao, Shandong, China

2009/12/8 Olori, Victor <volori_at_aviagen.com>

> Luan,
> You do not necessarily need to derive age2
>
> Lin(age) will fit age as a linear covariate while Pol(age,2) will fit
> age with a quadratic regression (Orthogonal polynomial) If you want to
> use Legendary polynomials then leg(age,2) will do the same thing.
> Regards
>
> Dr. Victor E. Olori
> ----------------------------
> Aviagen Limited
> Newbridge Midlothian
> EH28 8SZ
> Scotland, UK.
> Tel +44 (0) 131 3331056
>
>
> -----Original Message-----
> From: ASReml users discussion group [mailto:ASREML-L_at_DPI.NSW.GOV.AU] On
> Behalf Of luansheng
> Sent: 08 December 2009 00:24
> To: ASREML-L_at_DPI.NSW.GOV.AU
> Subject: How to write this code when age and it's quadratic fitted as
> fixed covariates in model?
>
> Dear friends:
> I am estimating the heritability of body weight in catfish. In this
> program, because there was difference between the birth date of
> families, fish's age and its quadratic were therefore fitted as fixed
> covariates in the analysis models of body weight. I didn't find the
> right function to finish this analysis, then I have to form a new
> covariate Age2 , which raised the Age field to the power 2. In the
> model, Age and Age2 were all fitted as the fixed covariates. My
> question is: Is this analysis process right? Can anyone give me better
> suggestions?
>
> !RENAME !ARGS 1
> AnimalID !P
> SireID !P
> DamID !P
> FamilyID * !ALPHA
> Birthday * !ALPHA
> Wt
> Lg
> Age
> Age2 !=Age !^2
> G0.dat !ALPHA !SKIP 1 !MAKE
> !DOPATH $1
> !PATH 1
> G0.dat !SKIP 1 !MAXIT 150 !EXTRA 5
> Wt ~ mu Age Age2 !r AnimalID
>
>
> Luan sheng
> ----------------------------------------
> Yellow Sea Fisheries Research Institute
> Chinese Academy of Fisheries Sciences
> Nanjing Road 106
> Qingdao 266071
> China
> luansheng_at_gmail.com
> luansheng_at_ysfri.ac.cn
> ----------------------------------------
>
> This message is intended for the addressee named and may contain
> confidential information. If you are not the intended recipient, please
> delete it and notify the sender. Views expressed in this message are
> those of the individual sender, and are not necessarily the views of
> their organisation.
>
>
>
>
> *****************************************************************************
> This e-mail contains confidential information and is intended solely for
> use by the individual named or entity to whom it is addressed. Please
> notify the sender and postmaster_at_aviagen.com immediately by e-mail if you
> have received this e-mail by mistake and delete this e-mail from your
> system. If you are not the named addressee, you are notified that
> disclosing, copying, distributing or taking any action in reliance on the
> contents of this information is strictly prohibited. Aviagen accepts no
> liability for the content of this e-mail and any consequences of actions
> taken in reliance on the contents, unless that information is subsequently
> confirmed in writing. Any views, opinions or actions presented in this
> e-mail are solely the author's and have not been approved by Aviagen. Any
> defamatory statements or infringing communication is contrary to Aviagen
> policy and outside the scope of the employment of the individual concerned.
> No employee or agent is authorized to conclude any binding agreement on
> behalf of Aviagen with another party by e-mail. Aviagen has taken reasonable
> precautions to ensure no viruses are present in this e-mail, but cannot
> accept responsibility for any loss or damage arising from the use of this
> e-mail or attachments.
>
> *****************************************************************************
>
> This message is intended for the addressee named and may contain
> confidential information. If you are not the intended recipient, please
> delete it and notify the sender. Views expressed in this message are those
> of the individual sender, and are not necessarily the views of their
> organisation.
>
>
Received on Thu Dec 09 2009 - 01:31:26 EST