機能は機能しますが（boot.stepAIC）、別の機能 - 環境問題の中でエラーを投げますか？

Question

私は今日Rコードで奇妙な動作を実現しました。 AICで段階的に回帰した結果のブートストラップ機能を含むパッケージ{boot.StepAIC}を試しました。しかし、私は統計的な背景がここに問題（私はそう望む）だとは思わない。
私はこれは私のコード例であるR.のトップレベルでの機能を使用することができます。

require(MASS) 
require(boot.StepAIC) 

n<-100 
x<-rnorm(n); y<-rnorm(n,sd=2); z<-rnorm(n,sd=3); res<-x+y+z+rnorm(n,sd=0.1) 
dat.test<-as.data.frame(cbind(x,y,z,res)) 
form.1<-as.formula(res~x+y+z) 
boot.stepAIC(lm(form.1, dat.test),dat.test) # should be OK - works at me 

しかし、私はそれを独自の関数でラップしたかったのです。私はデータとその式をその関数に渡します。

モデルの適合が strsplitで100ブートストラップ標本のエラーで失敗した（nam.vars、 "："）：しかし、私は言って）（boot.stepAIC内のエラーを取得する非文字引数

# custom function 
fun.boot.lm.stepAIC<-function(dat,form) { 
    if(!inherits(form, "formula")) stop("No formula given") 
    fit.lm<-lm(formula=form,data=dat) 
    return(boot.stepAIC(object=fit.lm,data=dat)) 
} 
fun.boot.lm.stepAIC(dat=dat.test,form=form.1) 
# results in an error 

だからここで間違いはありますか？私はそれが地方や地球環境と何か関係がなければならないと思われますね。

Answer 1

答えを提供anova test fails on lme fits created with pasted formulaのようにdo.callを使用します。

関数内で実行formへのアクセス権を持っていません。このような地球環境で再現することができます。我々はlm式としてform.1を使用して、そしてそれはが失敗します削除されていることを参照してください。

> form.1<-as.formula(res~x+y+z) 
> mm <- lm(form.1, dat.test) 
> mm$call 
lm(formula = form.1, data = dat.test) 
> rm(form.1) 
> boot.stepAIC(mm,dat.test) 
# same error as OP 

do.callを使用します。ここでも私はas.nameを使っています。そうでなければ、mmオブジェクトはその名前の代わりにデータセット全体を持ち歩きます。元の問題にこれを適用するに

> form.1<-as.formula(res~x+y+z) 
> mm <- do.call("lm", list(form.1, data=as.name("dat.test"))) 
> mm$call 
lm(formula = res ~ x + y + z, data = dat.test) 
> rm(form.1) 
> boot.stepAIC(mm,dat.test) 

、私はこれを行うだろう：

fun.boot.lm.stepAIC<-function(dat,form) { 
    if(!inherits(form, "formula")) stop("No formula given") 
    mm <- do.call("lm", list(form, data=as.name(dat))) 
    do.call("boot.stepAIC", list(mm,data=as.name(dat))) 
}  
form.1<-as.formula(res~x+y+z) 
fun.boot.lm.stepAIC(dat="dat.test",form=form1) 

これはあまりにも動作しますが、データセット全体が最終的な出力オブジェクト、およびコンソールへの最終的な出力に含まれますが、 同じように。

fun.boot.lm.stepAIC<-function(dat,form) { 
    if(!inherits(form, "formula")) stop("No formula given") 
    mm <- do.call("lm", list(form, data=dat)) 
    boot.stepAIC(mm,data=dat) 
}  
form.1<-as.formula(res~x+y+z) 
fun.boot.lm.stepAIC(dat=dat.test,form=form.1)