Re: ASreml-R additve genetic-maternal genetic covariance

From: Butler, David <David.Butler_at_DPI.QLD.GOV.AU>
Date: Fri, 17 Apr 2009 15:44:17 +1000

Craig, The syntax of link() is limited and virtually restricted to the random regression type case. General covariance structures that parallel the capabilities in asreml standalone are implemented in Release 3 with the "str()" function. Release 3 will be available in the coming weeks as soon as the documentation is updated. Dave. ________________________________ From: ASReml users discussion group [mailto:ASREML-L@DPI.NSW.GOV.AU] On Behalf Of craig walling Sent: Wednesday, 15 April 2009 9:24 PM To: ASREML-L@DPI.NSW.GOV.AU Subject: ASreml-R additve genetic-maternal genetic covariance Hi, I am trying to use ASreml-R to estimate the covariance between additive and maternal genetic variances for birth weight. The model runs fine without trying to specify the covariance, coded as: model4<-asreml(fixed= wt ~ 1, random = ~ped(animal, var=T) + ped(mother, var=T) + ide(mother, var=T), ginverse=list(animal=Dped, mother=Dped), data=juvtraits) However, I can't get the model to run linking animal and mother together. I have tried specifying this model as both: model5<-asreml(fixed= wt ~ 1, random = ~ us(link(ped(animal, var=T))):ped(mother, var=T) + ide(mother, var=T), ginverse=list(animal=Dped, mother=Dped), data=juvtraits) model6<-asreml(fixed= wt ~ 1, random = ~ us(link(~ped(animal, var=T))):ped(mother, var=T) + ide(mother, var=T), ginverse=list(animal=Dped, mother=Dped), data=juvtraits) neither of which will run. The first attempt give an error message "Error in switch(mode(x), "NULL" = structure(NULL, class = "formula"), : invalid formula". The second seems to make more of an attempt to run (it takes longer to return an error message), but then fails and gives "Error in data[[obj]] : object is not subsettable". Has anyone tried to do this before? Am I doing something obviously wrong? Is it possible to fit these covariances in ASreml-R, given that the information in the manual is very limited? Many thanks for any help you can offer Craig -------------------- 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 Sat Apr 17 2009 - 15:44:17 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.)