You don't get the mu term because brrank is place after the random
effects in the model. This means it is fitted before other effects (mu
in particular). Put brrank before the random effects and you will get a
mu term.
Your variable list contains one more than the number of data points on
your example line.
> -----Original Message-----
> From: Yuandan Zhang [SMTP:yzhang@metz.une.edu.au]
> Sent: Tuesday, 9 June 1998 14:30
> To: asreml@chiswick.anprod.CSIRO.AU
> Subject: quetions on sigularities in AI Matrix
>
> Hi, every one,
> I have not getten any response to my question that I posted last
> night. I
> guess I did not identify my questions well.
>
> Here are my questions regarding the output in asr file:
> 1. the Gamma term peculiar to animal (random effect) increased
> dramatically, from 0.1 to 7846.1 after 7 iterations. why?
> 2. There is a warning: 1 singularities in AI matrix. What does this
> mean?
> 3. How to identify which fixed effect is significant?
> 4. I included mu in the model, but the mu was zero in sln file.
>
> Thanks a lot
> Yuandan
>
> PS. Following are the as file and asr file for reference.
> the as file:
> Lamb fec and weight data from mating97 FEC -6/3, wt birth-6/3
> animal !P
> sire
> dam
> sex 2
> brrank 3
> sire_gd
> ewegd
> sirebd
> genotype 5
> lambplot
> bage
> wage
> bwt
> wwt
> wt4298
> wt6398
> FEC181297
> Nem181297
> FEC6398
> nem6398
> #97002 94087 93193 2 2 0 0 1 4 11 94 2.6 19.5 21.7 20.8 4
> 3 3 0
> wt_fec.txt
> wt_fec.txt
> wwt ~mu sex genotype !r animal !f brrank
>
> part of asr file
>
> Forming 410 equations: 9 dense
> NOTICE: 3 (more) singularities,
> LogL=-331.254 S2= 8.292 199 df 0.10000 0.00000
> 1.00000
> LogL=-331.012 S2= 7.681 199 df 0.20413 0.00000
> 1.00000
> LogL=-329.910 S2= 5.400 199 df 0.81331 0.00000
> 1.00000
> LogL=-327.708 S2= 1.568 199 df 5.89486 0.00000
> 1.00000
> LogL=-326.870 S2= 0.1704 199 df 64.84345 0.00000
> 1.00000
> LogL=-326.777 S2= 0.1577E-01 199 df 713.27793 0.00000
> 1.00000
> LogL=-326.768 S2= 0.1436E-02 199 df7846.05719 0.00000
> 1.00000
>
> Final parameter values 7846.12527 0.00000
> 1.00000
>
> WARNING: 1 singularities in AI matrix.
> Source Model terms Gamma Component Compnt/StndErr
> animal 398 398 7846.13 11.2667 9.97
> Variance 207 199 1.00000 0.143596E-02 0.00
> Solution Standard Error T-value
> T-prev
> 9 genotype 4 5.16 5.16 [ DF F_inc
> F_all ]
> 2 -1.96544 1.87758 -1.05
> 3 -3.62408 2.28853 -1.58
> -0.72
> 4 4.42244 1.71869 2.57
> 3.65
> 5 4.42008 1.74527 2.53
> 0.00
> 4 sex 1 3.83 3.99 [ DF F_inc
> F_all ]
> 7 -0.834937 0.417911 -2.00
> 1 animal 398 effects fitted
> 5 brrank 3 effects fitted
> Finished: 22:37:03.03 LogL Converged
>
>
>
>
>
> **********************************************************************
> **
> Yuandan Zhang
>
> Department of Animal Science
> University of New England
> Armidale, NSW 2351, Australia
>
> Phone +61 + 2 + 6773 2756(Lab)
>
>
> Fax +61 + 2 + 6773 3275
> E-mail yzhang@metz.une.edu.au
> homepage http://metz.une.edu.au/~yzhang/
> **********************************************************************
> **