Re: why "Singularity in Average Information Matrix"? how to do?

Dear Luan Sheng,

Further to my response last night,
What if genetic is just the dam component?
I am not very familiar with fish experiments but I understand that
sometimes families are raised in different tubs. Since in this case,
families are represented by sires, if these families were raised in
different tubs, then the sire variance is tub variance (in the sire dam
model), and the dam variance is primarily genetic. Under this scenario,
sire is analgous to 'maternal environment' effects in animal experiments.
I have therefore fitted sire + dam (allowing for 6 outliers and on the log
scale with the model, see below)
 log(wt) ~ mu age out(53) out(301) out(302) out(342) out(1004) out(631) !r
sire dam

This gives
 Source Model terms Gamma Component Comp/SE %
C
 dam 13 13 0.747327 0.104939 1.09 0
P
 sire 26 26 1.40185 0.196846 2.50 0
P
 Variance 1300 1292 1.00000 0.140419 25.18 0
P

The dam variance is still too large to provide a plausible heritability
though. That is, the residual still goes to the boundary in the 'sire +
animal' model because 0.104939 > 0.140419/3
The bottom line is that there appears to be family effects over and above
simple genetic effects. Maybe you need to replicate at the family level so
that you can partition the variance better.

May Jesus Christ continue to be gracious to you,

Arthur Gilmour, His servant .
 
Mixed model regression mapping for QTL detection in experimental crosses.
Computational Statistics and Data Analysis 51:3749-3764 at
http://dx.doi.org/10.1016/j.csda.2006.12.031

Profile:http://www.dpi.nsw.gov.au/research/staff/arthur-gilmour
      http://www.dpi.nsw.gov.au/reader/17263
Personal website: http://www.cargovale.com.au/

Skype: arthur.gilmour
mailto:Arthur.Gilmour_at_dpi.nsw.gov.au, arthur_at_cargovale.com.au
Principal Research Scientist (Biometrics)
NSW Department of Primary Industries
Orange Agricultural Institute, Forest Rd, ORANGE, 2800, AUSTRALIA
 
fax: 02 6391 3899; 02 6391 3922 Australia +61
telephone work: 02 6391 3815; home: 02 6364 3288; mobile: 0438 251 426

ASREML 2 is available from http://www.VSNi.co.uk/products/asreml
Register on the ASReml forum at http://www.vsni.co.uk/forum
ASReml-L Archives are at
https://listserv.dpi.nsw.gov.au/cgi-bin/wa.exe
Cookbook: http://uncronopio.org/ASReml

The ASReml3 alpha download registration page is:
http://www.austatgen.org/asreml3alpha.php
There are asreml3alpha forums at: http://www.austatgen.org/asreml3forum

Proposed travel:
MFX WS 24-28 November
Leave 22Dec to 9 Jan
><><><><><><><><><><><><><><><><><><><><><><><>