2017-03-03 32 views
0

我有一个调用对象,我想添加一个参数,并且我不想使用解析方式this answer does更改“调用”对象中的参数

所以,说我有一个lm对象,所以从lm

lma <- lm(mpg ~ cyl, data=mtcars) 
lma$call 
# lm(formula = mpg ~ cyl, data = mtcars) 

现在的口号,说我想补充的说法,weights=wt,使用呼叫。我意识到有一个非常简单的方法来创建一个新的通话,但我想知道我是否可以使用一个通话对象。也有编辑weights如果它已经在那里

lmb <- lm(mpg ~ cyl, data=mtcars, wei=wt) 
cl <- lmb$call 
wtpos <- which.max(pmatch(names(cl), "weights")) 
cl[[wtpos]] <- mtcars$qsec 
eval(cl) 

但因为在lma$call没有权重的说法,这将不会对lma$call工作的方式。

所以,感觉就像我应该能够通过添加另一个元素来简单地“增长”这个调用,但是我不知道该怎么做。例如,下面的失败:

cl <- lma$call 
cl <- c(cl, weights=quote(wt)) 
eval(cl) 
# [[1]] 
# lm(formula = mpg ~ cyl, data = mtcars) 
# 
# $weights 
# wt 

,我就希望得到的结果是一个新的“LM”对象等于lmb,不只是一个列表。

尽管有一种解决方法,那就是不使用解析(修改lm的副本使wt =权重为缺省值similar to in this solution),这同样不涉及编辑调用对象。

回答

1

我相信pryr包提供了一些有用的功能操纵呼叫:

lma <- lm(mpg ~ cyl, data=mtcars) 
lm_call <- lma$call 

library(pryr) 
modify_call(lm_call,list(weights = runif(32))) 

> lm_call2 <- modify_call(lm_call,list(weights = runif(32))) 
> eval(lm_call2) 

Call: 
lm(formula = mpg ~ cyl, data = mtcars, weights = c(0.934802365722135, 
0.983909613220021, 0.762353664264083, 0.23217184189707, 0.850970500381663, 
0.430563687346876, 0.962665138067678, 0.318865151610225, 0.697970792884007, 
0.389103061752394, 0.824285467388108, 0.676439745584503, 0.344414771301672, 
0.292265978176147, 0.925716639030725, 0.517001488478854, 0.726312294835225, 
0.842773627489805, 0.669753148220479, 0.618112818570808, 0.139365098671988, 
0.843711007386446, 0.851153723662719, 0.134744396666065, 0.92681276681833, 
0.00274682720191777, 0.732672147220001, 0.4184603120666, 0.0912447033915669, 
0.427389309043065, 0.721000595251098, 0.614837386412546)) 

Coefficients: 
(Intercept)   cyl 
    38.508  -2.945 

你可以看看里面pryr::modify_call,看看它在做什么,如果你想手工做的,我想。

+0

感谢您的提示。在Wickham的Advanced R下的Expressions:Calls:修改一个调用(在线地址为http://adv-r.had.co.nz/Expressions.html#calls),寻找这个结果将导致更简单的解决方案。 'cl < - lma $ call; cl $ weights < - quote(wt); eval(cl)' 给了我想要的结果。 –