関数内からlm（）を使用するにはどうすればよいですか？

Question

機能内またはlapplyからlm()を呼び出すと、フィットに関連付けられた$callがねじ込まれているようです。最小の作業例：

> library(MASS) 
> dat <- data.frame(x = 1:100, y=1:100) 
> dat <- within(dat, z <- x + log(y) + rnorm(100)) 
> fits <- lapply(list(z ~ x + y, z ~ x + log(y)), lm, dat) 
> stepAIC(fits[[1]])    # <-- error when I try to use the fit in other functions 
Error in eval(expr, envir, enclos) : could not find function "FUN" 

> fits[[1]]$call 
FUN(formula = X[[i]], data = ..1) # Aha -- this must be why -- $call is screwed up 

この問題を解決し、上記のエラーを防ぐにはどうすればよいですか？

Answer 1

時にはそれがlapplyに無名関数を供給すると良いでしょう：

fits <- lapply(list(z ~ x + y, z ~ x + log(y)), 
        function(f) lm(f, data = dat)) 
stepAIC(fits[[1]]) 
#works 

注DFがstepAICによって発見されていないため、この動作しません（通常は私の好ましい方法は、明示的なスコープ作る）こと：

fits <- lapply(list(z ~ x + y, z ~ x + log(y)), 
        function(f, DF) lm(f, data = DF), DF = dat) 

ステップワイズ回帰is a bad method anywayに注意してください。

Answer 2

代替youldは直接lapply内stepAICを適用すること：

AICs <- lapply(list(z ~ x + y, z ~ x + log(y)), 
       function(x) stepAIC(lm(x,dat))) 

これはあなたのすべてのモデルのstepAICの出力のリストを与えます。

Answer 3

lapplyの関数としてこれを使用してみてください。ノー@coffeinjunky

> fits[[1]] 

Call: 
lm(formula = z ~ x + y, data = dat) 

Coefficients: 
(Intercept)   x   y 
     2.154  1.031   NA 


> stepAIC(fits[[1]]) 
Start: AIC=-3.34 
z ~ x + y 


Step: AIC=-3.34 
z ~ x 

     Df Sum of Sq RSS AIC 
<none>     93 -3.34 
- x  1  88600 88693 680.78 

Call: 
lm(formula = z ~ x, data = dat) 

Coefficients: 
(Intercept)   x 
     2.154  1.031