Fault 1009

From: Gustavo A. Gutierrez <ggr03_at_IASTATE.EDU>
Date: Wed, 7 Feb 2007 13:53:01 -0600

Dear ASREML listing members,

I want to fit a random regression model for days open in dairy cattle, each cow
had information from its first to fith calving (5 records per cow). The
covariate is age at calving in months (acm). The code that I wrote is:

 anim * !P
 line * !A
 parity * !I
 ys * !I
 acm
 do
newped.txt !MAKE #pedigree file
datdofp1.csv !CSV #data file
do ~mu line parity ys leg(acm,-1) !r leg(acm,1).anim ide(anim).leg(acm,1) !f mv
0 0 2
leg(acm,1).anim 2
2 acm US !GP !+3
.27
0.08 0.07
1644 0 AINV
ide(anim).leg(acm,1) 2
1644 0 IDEN
2 0 US !GP !+3
0.16
0.147 0.21

There is 1644 individuals in the pedigree file, and 766 with records (data
file). Some cows have missing values for the response variable and some fixed
effects. At this point I worked under the assumption of homogenous residual
variance.

The code was not converged, and I got the following message from *.asr file.

NOTICE: 3 (more) singularities,
   1 LogL=-8844.07 S2= 5572.3 1784 df
   2 LogL=-8843.12 S2= 5583.3 1784 df : 2 components constrained
   3 LogL=-8842.98 S2= 5571.3 1784 df : 1 components constrained
 Warning: EM updates for 1 positive definite US structure(s).
   4 LogL=-8843.05 S2= 5520.3 1784 df : 2 components constrained
   5 LogL=-8843.05 S2= 5520.3 1784 df

 Warning: Since fault 1009 occured during the last iteration
           results reported may be erroneous

Last line read was:
    13 1009 3830 3830 1000
 Finished: 06 Feb 2007 13:55:52.437 Unable invert R or G [US?] matrix

What is the meaning of fault 1009?, and How I can solve it?

Thanks in advance,

Gustavo
Gustavo A. Gutierrez
Graduate student
Animal Breeding and Genetics
Iowa State University
227 Kildee Hall, Ames
phone office:515-2942712
IA 50011-3150
http://www.ans.iastate.edu/stud/grad/directory/?id=ggr03
Received on Sun Feb 07 2007 - 13:53:01 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.)