Re: [Rd] Bug in getVarCov.gls method (PR#9752)

2007-06-26 Thread ripley
Thank you, will be fixed in the next release of nlme (not yet scheduled, as everythign is code frozen for the release of R 2.5.1 tomorrow). On Tue, 26 Jun 2007, Andrzej Galecki wrote: > > Two attachments: > > 1. getVarCovBugReport.R - Rcode with an example illustrating the problem and > how to

Re: [Rd] Bug in getVarCov.gls method (PR#9752)

2007-06-26 Thread agalecki
This is a multi-part message in MIME format. --000206090008050103050209 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit Two attachments: 1. getVarCovBugReport.R - Rcode with an example illustrating the problem and how to fix it 2. getVarC

Re: [Rd] Bug in getVarCov.gls method (PR#9752)

2007-06-26 Thread ripley
On Mon, 25 Jun 2007, [EMAIL PROTECTED] wrote: > I am using R2.5 under Windows. I presume you mean 2.5.0 (there is no R2.5: see the posting guide). But which version of nlme, which is the relevant fact here? The R posting guide suggests showing the output of sessionInfo() to establish the env

[Rd] Bug in getVarCov.gls method (PR#9752)

2007-06-26 Thread agalecki
Hello, I am using R2.5 under Windows. Looks like the following statement vars <- (obj$sigma^2)*vw in getVarCov.gls method (nlme package) needs to be replaced with: vars <- (obj$sigma*vw)^2 With best regards Andrzej Galecki Douglas Bates wrote: >I'm not sure when the getVarCov.