[R] "eval" inside a function call in connection with updating the data slot in the call of lmer

[Duncan Murdoch]

On 12-09-17 6:26 PM, Søren Højsgaard wrote:
> Dear list,
> Given a linear mixed model (from lme4) I want to 1) first change the input dataset and then 2) change the model formula. I want this to happen in a function call;
> Please see below. Options 1) and 2) below work whereas 3) fails with the message

I get the failure in 1), not in 3).  I think it's a bug.  The problem 
appears to be that the update method for mer objects looks in the wrong 
place for its variables.  It is looking in parent.frame() (i.e. in the 
caller), but the caller isn't you.  The method should be looking in 
parent.frame(2).

You can fix this yourself if you have the lme4 source (it's line 1483 of 
lmer.R in version 0.999999-0) but you probably want to send your sample 
code to the maintainers.  Making that change might break something else.

Duncan Murdoch






>> foo()
> Error in is.data.frame(data) : object 'beets2' not found
>
> Question: What is it one must to in case 3) to have R look "inside" the function to figure out what "beets2" is?
>
> Best regards
> Søren
> ________________
>
> library(pbkrtest)
> data(beets)
> lgs    <- lmer(sugpct~block+sow+harvest+(1|block:harvest), data=beets, REML=F)
>
> foo <- function(){
> 	## 1)
> 	beets2 <- transform(beets, yy = sugpct * yield)
> 	ma1    <- lmer(yy~block+sow+harvest+(1|block:harvest), data=beets2, REML=F)
> 	ma0    <- update(ma1, yy~.)
> 	## 2)
> 	cl <- getCall(lgs)
> 	cl[["data"]] <- beets2
> 	mb1 <- eval(cl)
> 	mb0 <- update(mb1, yy~.)
> 	mb0
> 	## 3)
> 	cl <- getCall(lgs)
> 	cl[["data"]] <- as.name("beets2")
> 	mc1 <- eval(cl)
> 	mc0 <- update(mc1, yy~.)
> 	mc0
> }
> foo()
>
> ______________________________________________
> R-help at 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.
>