[R-sig-ME] Priors for MCMCglmm with haplodiploid relatedness matrix

[Jarrod Hadfield]

Hi Mike,

Can you tell us how haplo.mat.inverse is formed or provide the code. The 
error suggests that the matrix is (near) singular.

Cheers,

Jarrod



On 04/10/2016 17:00, Michael Sheehan wrote:
> Dear all,
>
>
> I am attempting to examine the heritability of a number of traits in a haplodiplod species using MCMCglmm. This has presented a number of challenges since the standard relatedness matrix calculated from a pedigree as part of the program assumes diploid relatedness, which underestimates relatedness among sisters. I have been able to make a haplodiploid appropriate relatedness matrix but when I run the following:
>
>
> model1=MCMCglmm(trait~1, random=~animal, ginverse=list(animal=haplo.mat.inverse), data=data, family="ordinal", prior=prior1, nitt=105000, thin=100, burnin=5000, verbose=T)
>
>
> I immediately get the following error:
>
>
> "G-structure 1 is ill-conditioned (possibly because of ginverse): use proper priors if you haven't, or rescale data if you have"
>
>
> If I run the same code with a diploid relatedness matrix (calling it through pedigree or ginverse) the program runs without error. I have tried a wide range of priors and even run without priors and get the same error for the haplodiploid matrix while the diploid matrix runs.
>
>
> I am curious if others have received a similar message and how it was dealt with.
>
>
> Thanks in advance for any help!
>
> Mike Sheehan
>
>
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models


-- 
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.