Re: regression analysis(2)

From: Olori, Victor <volori_at_AVIAGEN.COM>
Date: Wed, 17 Feb 2010 10:06:40 -0000

Hi Mariam,

I think you probably have too many elements in you starting covariance
matrix for a quadratic polynomial. i.e. for pol(assess,2) you need only
6 elements. This is the error message coming from your .asr file at
present. I have modified you code below by putting a hash at the
beginning of the last line of you starting values. Personally instead of
using ide(anim), I will create another column in the data with the hen
ids exactly the same as the animal id in each row but call this new
column HEN I will use this as the second random term to take care of the
permanent environmental term but lets try your approach. . I have
included HEN it in your as file below but commented it out for now. I
don't know how you laid out your data and pedigree files. Can I assume
you have a pedigree file indexped.txt with just animal sire and dam? Do
you have repeated animal ids in this pedigree file? I would ensure I
have unique animals. Now you need to sort your data by ASSESS so that
(assuming this trait ranges from 1 to 13) all animals with a value at
'assess' class 1 will be at the top followed by 2 etc to 13. This is
very important if you assume heterogeneous error variance which right
now you don't so perhaps it will not matter but I will do it to be on
the safe side. Another thing is that it is sometimes better to scale
you age trait (assess) especially if it does not range neatly from 1 to
13 but those are final refinements.

So Try this ;

 RECORD !M -99
# HEN 2518 # where 2518 is the number (levels) of animals in
the file.
alllnew.txt !MAXIT 20
pol(ASSESS,2).ide(ANIM) !f mv
0 0 2
pol(ASSESS,2).ANIM 2
3 0 US !GP +6
0.01 0.1
0.01 0.01 0.1
#0.01 0.01 0.01 0.1
ANIM 0 ainv
pol(ASSESS,2).ide(ANIM) 2
3 0 US !GP +6
0.01 0.1
0.01 0.01 0.1
#0.01 0.01 0.01 0.1
ide(ANIM) 0 0




Good luck.


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 maryam
Sent: 17 February 2010 08:42
Subject: Re: regression analysis(2)


Dear Ruhallah,
Thanks so much for your reply,
but again when i run the job i cant get results!
please if you have time help me more,
all the best,


-------------------- m2f --------------------

Sent using Mail2Forum (

Read this topic online here:

-------------------- m2f --------------------

31.asr <>
31.AS <>

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 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 Feb 17 2010 - 10:06:40 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.)