更改“调用”对象中的参数

pdb*_*pdb 3 r

我有一个调用对象,我想添加一个参数,并且我不想像这个答案那样使用解析。

所以,假设我有一个lm对象,那么来自lm

lma <- lm(mpg ~ cyl, data=mtcars)
lma$call
# lm(formula = mpg ~ cyl, data = mtcars)
Run Code Online (Sandbox Code Playgroud)

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)
Run Code Online (Sandbox Code Playgroud)

但这不起作用,lma$call因为 中没有权重参数lma$call

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

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

所以,我希望结果是一个等于 的新“lm”对象lmb,而不仅仅是一个列表。

虽然有一个解决方法,但它不使用解析(修改 lm 的副本以将 wt=weights 作为默认值,类似于此解决方案),同样,也不涉及编辑调用对象。

jor*_*ran 5

我相信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  
Run Code Online (Sandbox Code Playgroud)

pryr::modify_call我想,如果您想手动执行此操作,您可以查看内部以了解它在做什么。

  • 谢谢你的提示。在 Wickham 的 Advanced R 中的 Expressions:Calls:Modifying a call(在线地址:http://adv-r.had.co.nz/Expressions.html#calls)下搜索此问题可以找到更简单的解决方案。`cl &lt;- lma$call; cl$权重 &lt;- 报价(wt); eval(cl)` 给了我我想要的结果。 (2认同)