On Jun 25, 2009, at 11:35 AM, Rebecca Sela wrote:
I am trying to estimate models with subsets using the NLME package.
However, I am getting an error in the case below (among others):
subset <- c(rep(TRUE, 107), FALSE)
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1,
subset=subset)
Error in xj[i] : invalid subscript type 'closure'
fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1,
subset=1:107)
The second estimation works.
Does anyone know if there is another work-around? (I have also e-
mailed the package maintainers, but one of the e-mails bounced, so I
am trying this list as well.)
Hard to tell. Closure is a type of function. It appears that R might
be mistaking its own function "subset" for the object that you
intended to pass. Perhaps if you stopped calling your dogs, "Dog"?
> library(fortunes)
> fortune("dog")
Thank you!
Rebecca
--
Rebecca Sela
IOMS/Statistics Group
Stern School of Business
New York University
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
David Winsemius, MD
Heritage Laboratories
West Hartford, CT
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.