Hello,
1. Illustrating example:
sessionInfo at the end of this email
library(nlme)
fm1 <- gls(weight ~ Time * Diet, BodyWeight)
Variogram(fm1, form = ~ Time | Rat)[1:10,]
Error in FUN(X[[1L]], ...) : unused argument(s) (method = "euclidean")
2. It appears that the error is caused by the sta
Hello,
I would like to propose to make a minor change in method print.VarCov
by inserting one line
stdevs <- FALSE # Inserted line
as shown later in this email.
By adding this line we would modify the output from
print(corrmtx) as shown in the illustratio
This is a multi-part message in MIME format.
--040502080208060001050807
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
Hello,
1. It appears that "coef<-.varIdent" method does not work properly in
some instances.
Execution error:
Error in
This is a multi-part message in MIME format.
--000706080709050507050106
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
Hello,
1. The following code (below) and output
illustrate that
argument fixed= in varIdent function is ignored during
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
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.