[R] update.default evaluating in wrong environment?
Martin C. Martin
martin at martincmartin.com
Tue Oct 10 19:32:29 CEST 2006
Gabor Grothendieck wrote:
> As a workaround use evaluate=FALSE argument to update and
> evaluate it yourself fetching the environment from the innards
> of the lm structure:
>
> f <- function() {
> DF <- data.frame(y = 1:12, x1 = gl(2, 1, 12), x2 = gl(2,6))
> lm(y ~ x1, DF)
> }
>
> f.lm <- f()
> e <- attr(terms(f.lm), ".Environment")
> eval(update(f.lm, formula = y ~ x2, evaluate = FALSE), e)
Thanks, or even just:
e <- environment(formula(f.lm))
But this was more of a bug report. Is update.default wrong? Should it
be changed? I don't see how evaluating in update's parent environment
would ever be better default behavior than the formula's environment.
- Martin
>
>
> On 10/10/06, Martin C. Martin <martin at martincmartin.com> wrote:
>> Hi all,
>>
>> update.default, which is the method used to update "lm" objects (among
>> others), extracts the "call" element from it's first argument, updates
>> it, then evaluates it in the parent.frame(). Shouldn't it be evaluated
>> in environment(formula(object)), if that's non-NULL?
>>
>> I ask because I call "lm" from within a function, and the data argument
>> is a local variable of that function. After that, I can't update the
>> model any more, since the new lm() call (the one evaled in
>> parent.frame()) can't find the data.
>>
>> Best,
>> Martin
>>
>> ______________________________________________
>> R-help at 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.
>>
More information about the R-help
mailing list