Одним из решений является создание формулы с использованием paste()
и преобразование ее в формулу:
> ## your example plus some dummy data
> var1 <- "V001"
> var2 <- "V002"
> var3 <- "V003"
> dat <- data.frame(V001 = runif(10), V002 = runif(10), V003 = runif(10))
> f <- formula(paste(var1, "~", var2, "+", var3))
Теперь мы можем посмотреть на f
> f
V001 ~ V002 + V003
> class(f)
[1] "formula"
на самом деле это формула.Теперь мы можем передать это в rlm()
в качестве первого аргумента:
> require(MASS)
> mod <- rlm(f, data = dat)
> mod
Call:
rlm(formula = f, data = dat)
Converged in 8 iterations
Coefficients:
(Intercept) V002 V003
0.2725538 -0.1281576 0.1617250
Degrees of freedom: 10 total; 7 residual
Scale estimate: 0.251
HTH