quetions on sigularities in AI Matrix
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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/ 
************************************************************************