Re: [R] Updating lmer - object is not subsettable?

From: Jarrett Byrnes <jebyrnes_at_ucdavis.edu>
Date: Wed 13 Sep 2006 - 04:10:49 GMT

Of course. The traceback is as follows. If you wish, I can privately email you the data, as well as the function I'm working on.

  1. all.names(x)
  2. inherits(x, "factor")
  3. is.factor(table)
  4. match(x, table, nomatch = 0)
  5. "/" %in% all.names(x)
  6. slashTerms(x[[3]])
  7. FUN(X[[1]], ...)
  8. lapply(bb, function(x) { if (is.list(trms <- slashTerms(x[[3]]))) return(lapply(unlist(makeInteraction(trms)), function (trm) substitute(foo | bar, list(foo = x[[2]], bar = trm)))) x })
  9. unlist(lapply(bb, function(x) { if (is.list(trms <- slashTerms(x[[3]]))) return(lapply(unlist(makeInteraction(trms)), function (trm) substitute(foo | bar, list(foo = x[[2]], bar = trm)))) x }))
  10. expandSlash(findbars(formula[[3]]))
  11. lmer(formula = sim.data ~ (1 | Trial), data = exp.data, method = "ML")
  12. lmer(formula = sim.data ~ (1 | Trial), data = exp.data, method = "ML")
  13. eval(expr, envir, enclos)
  14. eval(call, parent.frame())
  15. .local(object, ...)
  16. update(model.2, sim.data ~ .)
  17. update(model.2, sim.data ~ .)

On Sep 12, 2006, at 8:16 PM, Douglas Bates wrote:

> Can you show a traceback on this example? It may be related to a
> problem that I just fixed in the development version of the lme4
> package.
>
> Alternatively if you can make the data available I can generate a
> traceback myself.
>
>
> On 9/12/06, Jarrett Byrnes <jebyrnes@ucdavis.edu> wrote:
>> I'm attempting to write a general function to implement Faraway's
>> bootstrapping algorithm for mixed models with lmer, but have run into
>> a curious problem. I'm comparing two models
>>
>> model.1<-lmer(Response ~ Treatment + (1|Trial), data=exp.data,
>> method="ML")
>> model.2<-lmer(Response ~ 1 + (1|Trial), data=exp.data, method="ML")
>>
>>
>> When I attempt to update model.2 with simulated data, however, I get
>> the following error:
>>
>> sim.data<-unlist(simulate(model.1))
>> sim.model.2<-update(model.2, sim.data~.)
>>
>> Error in x[[3]] : object is not subsettable
>>
>>
>> Now, the following
>> sim.model.1<-update(model.1, sim.data~.)
>>
>> appears to work just fine. Does anyone know why update won't work,
>> and is there something I can do about this?
>>
>> -Jarrett
>>
>> ______________________________________________
>> R-help@stat.math.ethz.ch 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.
>>



R-help@stat.math.ethz.ch 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. Received on Wed Sep 13 14:14:07 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 13 Sep 2006 - 04:30:07 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.