Re: ASReml-R bivariate example

Hi, The following works fine for me in Version 2 orange <- asreml.read.table("orange.dat",header=T) orange$gfw <- orange$gfw/10 orange$fdiam <- orange$fdiam/10 wth.asr <- asreml(cbind(gfw,fdiam) ~ trait+trait:Year, random = ~ us(trait,init=c(0.4,0.3,1.3)):Team + us(trait,init=c(0.2,0.2,2.0)):Tag, rcov = ~ units:us(trait,init=c(0.2,0.2,0.4)), data = orange, trace=T) using the data > orange[1:5,] Tag X1 X2 Team Year gfw X3 fdiam 1 101 3 21 1 1 5.6 743 18.5 2 101 3 21 1 2 6.0 712 19.6 3 101 3 21 1 3 8.0 757 21.5 4 102 3 21 1 1 5.3 709 20.8 5 102 3 21 1 2 5.7 661 20.9 > asreml(): 2.00 Library: 2.00bi Run: Tue Apr 21 14:26:19 2009 Equations: 1120 (8 dense) Initial update shrinkage factor: 0.316 2 singularities detected in design matrix LogLik S2 DF -886.5213 1.0000 2964 14:26:20 -818.5083 1.0000 2964 14:26:20 -755.9111 1.0000 2964 14:26:20 -725.3739 1.0000 2964 14:26:20 -723.4750 1.0000 2964 14:26:20 -723.4617 1.0000 2964 14:26:20 -723.4617 1.0000 2964 14:26:20 Finished on: Tue Apr 21 14:26:20 2009 Can you send me some more details of the problem? Dave. ________________________________ From: ASReml users discussion group [mailto:ASREML-L@dpi.nsw.gov.au] On Behalf Of petergs Sent: Tuesday, 21 April 2009 10:40 AM To: ASREML-L@dpi.nsw.gov.au Subject: ASReml-R bivariate example G'day, Everyone. I've been trying to run the example bivariate analysis on p. 34 of the ASReml-R Reference Manual (Feb., 2007) but the code generates a termination error. The data file I've used is wether.dat (with columns renamed appropriately) supplied with ASReml-W. The code given in the ASReml User Guide (Release 3, Jul., 2008) for the same analysis, on the other hand, runs happily in ASReml-W. I'd be very pleased if someone could show me how to get the ASReml-R code working! -------------------- m2f -------------------- Sent using Mail2Forum (http://www.mail2forum.com). Read this topic online here: http://www.vsni.co.uk/forum/viewtopic.php?p=492#492 -------------------- 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 Wed Apr 21 2009 - 14:36:14 EST