Re: ASReml-R and using weights with the ped() function

From: Butler, David <David.Butler_at_DPI.QLD.GOV.AU>
Date: Thu, 11 Dec 2008 10:26:41 +1000

Dear kiwi, Weights and ginv structures reside internally in contiguous memory locations so from your description it is likely there is a bug in the code. I'm on leave and will take a look as soon as I can. Perhaps you can forward me your script and data file for testing. dave. ________________________________ From: ASReml users discussion group [mailto:ASREML-L@DPI.NSW.GOV.AU] On Behalf Of kiwi Sent: Wednesday, 10 December 2008 9:18 AM To: ASREML-L@DPI.NSW.GOV.AU Subject: Re: ASReml-R and using weights with the ped() function I've not received any comments about this one yet. There's a further development: I mentioned the 'singularity' error I got when I ran this: Code: > discard.asr <- asreml(fixed = Discard_Total ~ 1 + Season, random = ~ped(Vine, var = T) + Block, data = model.df, na.method.X = "include", ginverse = list(Vine = G.inv$ginv)) However, it didn't prevent convergence. I mentioned that adding a weights argument causes an error to do with a Code: Error : Abnormal termination Singularity in Average Information Matrix result. At that time, I considered it possible that the singularity had a cabability of creating a problem with weighting, but otherwise it was benign. I've found an instance where even without a singularity problem, it's still a problem. The first one works fine; the second falls over: Code: resp.asreml <- try(asreml(fixed = form, random = ~ped(Vine, var = T) + Block, data = model.df, na.method.X = "include", ginverse = list(Vine = G.inv$ginv))) model.df$Wt <- rep(1, nrow(model.df)) respW.asreml <- try(asreml(fixed = form, random = ~ped(Vine, var = T) + Block, weights = Wt, data = model.df, na.method.X = "include", ginverse = list(Vine = G.inv$ginv))) The 'singularity' problem returns even though it would appear to the ignorant that the same model was being tried. That leads me to conclude that it can't just be me. So my question is this: Has anyone ever got weights to work when using the ped() function? -------------------- 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 Dec 11 2008 - 10:26:41 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.)