Re: Singularity In Average Information Matrix Using Ainverse

Thank you for your response. I understand what you're saying and would now like to use a sire or sire+MGS model but am unsure how this should be setup.
This is what I have so far:

geno = read.table(file='geno.txt', header=TRUE)
SNP = read.table(file='SNP.txt', header=TRUE)
geno$Animal<-as.factor(geno$Animal)
test.asr = asreml(Fert ~ SNP, random = ∼ giv(Sire, var=T), ginverse=list(Sire=sire.giv), data=geno)

The MGS-3.ped.txt needs to be altered to make the sire.giv file but I am not sure how. The example in the ASREML-R manual setsup sire.giv as:

       row column value
Sire 1 1 1 0.5
Sire 2 2 2 0.5
Sire 3 3 3 0.5
Sire 4 4 4 0.5
Sire 5 5 5 0.5
Sire 6 6 6 0.5
Sire 7 7 7 0.5
Sire 8 8 8 0.5
Sire 9 9 9 0.5

I'm unsure how this would relate to my pedigree file and how I could set mine up similarly.

Any help would be much appreciated- I'm very thankful for the guidance you have provided.

Michelle

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

Sent using Mail2Forum (http://www.mail2forum.com).

Read this topic online here:
http://www.vsni.co.uk/forum/viewtopic.php?p=869#869

-------------------- m2f --------------------
Received on Thu Sep 23 2009 - 20:06:50 EST