chức năng hoạt động (boot.stepAIC) nhưng ném một lỗi bên trong một vấn đề môi trường chức năng khác?

Question

Tôi nhận ra một hành vi lạ hôm nay với mã R của tôi. Tôi đã thử một gói {boot.StepAIC} bao gồm một hàm bootstrap cho kết quả hồi quy từng bước với AIC. Tuy nhiên tôi không nghĩ rằng nền thống kê là ở đây vấn đề (tôi hy vọng như vậy).
Tôi có thể sử dụng chức năng ở cấp cao nhất của R. Đây là mã ví dụ của tôi.

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 

Tuy nhiên, tôi muốn kết hợp với chức năng riêng. Tôi chuyển dữ liệu và công thức đến hàm đó. Nhưng tôi nhận được một lỗi trong vòng boot.stepAIC() nói:

mô hình phù hợp thất bại trong 100 mẫu bootstrap Lỗi trong strsplit (nam.vars, ":"): đối số phi vật

# 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 

Vậy đâu là sai lầm? Tôi cho rằng nó phải có một cái gì đó để làm với môi trường địa phương và toàn cầu, phải không?

Answer 1

Sử dụng do.call như trong anova test fails on lme fits created with pasted formula cung cấp câu trả lời.

không có quyền truy cập vào form khi chạy trong một hàm; có thể được tái tạo trong môi trường toàn cầu như thế này; chúng ta thấy rằng lm đang sử dụng form.1 làm công thức và xóa nó làm cho không thành công.

> 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 

Sử dụng do.call không hoạt động. Ở đây tôi cũng sử dụng as.name; nếu không, đối tượng mm sẽ mang theo toàn bộ tập dữ liệu thay vì chỉ tên của nó.

> 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) 

Để áp dụng điều này cho vấn đề ban đầu, tôi muốn làm điều này:

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) 

này hoạt động quá nhưng toàn bộ tập dữ liệu được bao gồm trong đối tượng đầu ra cuối cùng, và kết quả cuối cùng để an ủi , cũng.

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)