Re: asreml-r multivariate animal model

From: Butler, David <David.Butler_at_DPI.QLD.GOV.AU>
Date: Mon, 23 Mar 2009 08:31:00 +1000

MM, You don't need the "var=T" argument to ped() when it is part of an interaction and the other term is a variance model. Dave. ________________________________ From: ASReml users discussion group [mailto:ASREML-L@DPI.NSW.GOV.AU] On Behalf Of MM Sent: Sunday, 22 March 2009 2:31 AM To: ASREML-L@DPI.NSW.GOV.AU Subject: asreml-r multivariate animal model Hi, I have followed the asreml-R manual example on page 70 and I seem to have been able to make univariate animal models run for a small simulated dataset of my own. However, I am getting an error about a singularity in a bivariate analysis that I suspect may be coming from a problem in specifying the model rather than a true lack of information . The data are of 100 individuals measured for two traits, with overlapping maternal and paternal sibships. I've simulated no data for the parents. My model and error are copied below. But more generally, is there any place I can go for examples of different types of animal models coded for asreml-R? Thanks, MM <>< > test4<-asreml(cbind(P1,P2)~trait, random=~us(trait,init=c(0.1,0,0.1)):ped(ID, var=T), ginverse=list(ID=data2.ainv), rcov=~units:us(trait, init=c(2,0,2)), data=data5) asreml(): 2.00 Library: 2.00bi Run: Sat Mar 21 16:23:09 2009 Equations: 362 (2 dense) Initial update shrinkage factor: 0.316 LogLik S2 DF 1195.8666 1.0000 198 16:23:09 (6 component(s) constrained) Warning: 1 singularities inAverage Information Matrix Exit status: -1 - Singularity in Average Information Matrix -------------------- m2f -------------------- Sent using Mail2Forum ( Read this topic online here: -------------------- m2f -------------------- This message is intended for the addressee named and may contain confidential information. If you are not the intended recipient, please delete it and notify the sender. Views expressed in this message are those of the individual sender, and are not necessarily the views of their organisation. ********************************DISCLAIMER**************************** The information contained in the above e-mail message or messages (which includes any attachments) is confidential and may be legally privileged. It is intended only for the use of the person or entity to which it is addressed. If you are not the addressee any form of disclosure, copying, modification, distribution or any action taken or omitted in reliance on the information is unauthorised. Opinions contained in the message(s) do not necessarily reflect the opinions of the Queensland Government and its authorities. If you received this communication in error, please notify the sender immediately and delete it from your computer system network.
Received on Tue Mar 23 2009 - 08:31:00 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.)