Hello,
I was wondering whether there had been any progress on this bug yet?
I have checked for updates to the package and it doesn't seem that there
have been any, but I might be mistaken!
Similarly to Mark, I am getting the following error message when running a
Gls with two (factor) covariates (
Hi Mark,
I'm sorry I forgot about this. I found that Gls was out of date with
respect to gls. It's now fixed and you can
source('http://biostat.mc.vanderbilt.edu/tmp/Gls.s') to override the old
version with the fixed version until the next release of rms. You won't
need the workaround below.
Fr
Frank, sorry if I seem impatient, but do you know approximately when
this will be fixed?
In the meantime, if I have
library(nlme)
library(rms)
d <- data.frame(x = rnorm(50), y = rnorm(50))
and I'm interested in the 3 df test for x that I would get from
anova(Gls(y ~ rcs(x, 4), data=d, correlati
Appears to be a definite bug, probably caused by having more than one
correlation parameter. I hope to have this fixed within 3 days.
Frank
Mark Seeto wrote
>
> Dear R-help,
>
> I don't understand why Gls gives me an error when trying to fit a
> model with AR(2) errors, while gls (from nlme) do
Dear R-help,
I don't understand why Gls gives me an error when trying to fit a
model with AR(2) errors, while gls (from nlme) does not. For example:
library(nlme)
library(rms)
set.seed(1)
d <- data.frame(x = rnorm(50), y = rnorm(50))
gls(y ~ x, data=d, correlation = corARMA(p=2)) #This works
Gl
5 matches
Mail list logo
