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

From: Olori, Victor <volori_at_AVIAGEN.COM>
Date: Tue, 8 Dec 2009 09:47:22 -0000

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.
*****************************************************************************
Received on Wed Dec 08 2009 - 09:47:22 EST

This webpage is part of the ASReml-l discussion list archives 2004-2010. More information on ASReml can be found at the VSN website. This discussion list is now deprecated - please use the VSN forum for discussion on ASReml. (These online archives were generated using the hypermail package.)