Я встречаю странное поведение при вызове lm
внутри lapply
с помощью аргумента weights
.
Мой код состоит из списка формул, в котором я запускаю линейную модель, которую я вызываю в lapply
. Пока он работал:
dd <- data.frame(y = rnorm(100),
x1 = rnorm(100),
x2 = rnorm(100),
x3 = rnorm(100),
x4 = rnorm(100),
wg = runif(100,1,100))
ls.form <- list(
formula(y~x1+x2),
formula(y~x3+x4),
formula(y~x1|x2|x3),
formula(y~x1+x2+x3+x4)
)
res.no.wg <- lapply(ls.form, lm, data = dd)
Однако, когда я добавляю аргумент weights
, я получаю странную ошибку:
res.with.wg <- lapply(ls.form, lm, data = dd, weights = dd[,"wg"])
Error in eval(extras, data, env) :
..2 used in an incorrect context, no ... to look in
Как будто ...
из lapply
противоречил ...
вызова lm
, но только из-за аргумента weights
.
Любая идея была причиной этой проблемы и как ее исправить?
ПРИМЕЧАНИЕ: использование вызова без lapply
работает как ожидалось:
lm(ls.form[[1]], data = dd, weights = dd[,"wg"] )
Call:
lm(formula = ls.form[[1]], data = dd, weights = dd[, "wg"])
Coefficients:
(Intercept) x1 x2
-0.12020 0.06049 -0.01937
РЕДАКТИРОВАТЬ Последний вызов - lapply
в пределах function
типа:
f1 <- function(samp, dat, wgt){
res.with.wg2 <- lapply(ls.form, function(x) {lm(formula = x, data=dat[samp,], weights=dat[samp,wgt])})
}
f1(1:66, dat=dd, wgt = "wg")