Re: Singularity In Average Information Matrix Using Ainverse

From: Michelle <asremlforum_at_VSNI.CO.UK>
Date: Wed, 23 Sep 2009 20:06:50 +0100

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

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.)