プログラムで構築された数式を使用するうえでの落とし穴はありますか？

Question

私は、潜在的な説明変数であるベクトル を順番に回帰したいと考えています。 を一緒に貼り付けるのではなく、reformulate()、 as demonstrated hereと考えています。

以下の機能fun()は、希望のモデルに合った仕事をしているようです。しかし、 は、の値ではなく、という名前の構成式オブジェクト のという名前を呼び出し要素に記録することに注意してください。

## (1) Function using programatically contructed formula 
fun <- function(XX) { 
    ff <- reformulate(response="mpg", termlabels=XX) 
    lm(ff, data=mtcars) 
} 
fun(XX=c("cyl", "disp")) 
# 
# Call: 
# lm(formula = ff, data = mtcars)     <<<--- Note recorded call 
# 
# Coefficients: 
# (Intercept)   cyl   disp 
# 34.66099  -1.58728  -0.02058 

## (2) Result of directly specified formula (just for purposes of comparison) 
lm(mpg ~ cyl + disp, data=mtcars) 
# 
# Call: 
# lm(formula = mpg ~ cyl + disp, data = mtcars) <<<--- Note recorded call 
# 
# Coefficients: 
# (Intercept)   cyl   disp 
# 34.66099  -1.58728  -0.02058 

私の質問：はこの内の任意の危険性はありますか？たとえば、後でupdate、またはpredictまたは をモデルフィットオブジェクトに適用したい場合（おそらく他の環境から）、これは 問題になることがありますか？

それにもかかわらず、記録されている のコール権を取得するのは少し面倒な選択肢です。eval(substitute())を使用します。どのようにしてもこれは一般により安全な構造ですか？

fun2 <- function(XX) { 
    ff <- reformulate(response="mpg", termlabels=XX) 
    eval(substitute(lm(FF, data=mtcars), list(FF=ff))) 
} 
fun2(XX=c("cyl", "disp"))$call 
## lm(formula = mpg ~ cyl + disp, data = mtcars) 

Answer 1

は、私はいつもここで何かR環境とかむかもしれないスコープを含む何状況が存在しないと主張することを躊躇んだけど、...もう少し探検した後、私の最初の使用は、上記安全に見えるん。

印刷された電話はちょっとした赤ちゃんです。

実際が他の機能（およびformula()とas.formula()により抽出されたもの）によって使用されますが、フィットオブジェクトのterms要素に格納されているものであることを式、及びそれ実際式権を取得します。 （terms要素は、ちょうど付属属性の束と"formula"あるクラス"terms"のオブジェクトが含まれています。）

はまで（私の質問の提案と関連するコメントがすべて同じ"formula"オブジェクトを格納することを確認するには関連する環境）、以下を実行します。

## First the three approaches in my post 
formula(fun(XX=c("cyl", "disp"))) 
# mpg ~ cyl + disp 
# <environment: 0x026d2b7c> 

formula(lm(mpg ~ cyl + disp, data=mtcars)) 
# mpg ~ cyl + disp 

formula(fun2(XX=c("cyl", "disp"))$call) 
# mpg ~ cyl + disp 
# <environment: 0x02c4ce2c> 

## Then Gabor Grothendieck's idea 
XX = c("cyl", "disp") 
ff <- reformulate(response="mpg", termlabels=XX) 
formula(do.call("lm", list(ff, quote(mtcars)))) 
## mpg ~ cyl + disp 

formula()が本当にフィットオブジェクトのterms要素から、その出力を導出していることを確認するために、stats:::formula.lmとstats:::formula.termsを見てみましょう。