Анкова с помощью переменной, выбранной из цикла (R)

Question

У меня есть 1000 зависимых переменных (у) для использования в Ancova. Независимые переменные (x) одинаковы во всех моделях.

Yvar1, Yvar2, …, Yvar1000, x1, x2
1       2          5       11  16
2       3          6       18  23

Мне нужен код R, который заменяет example$Y.to.use значения для каждого Ивара.

Я попробовал это:

for (i in 1:1000){ 
example$Y.to.use<- example$paste("Yvar",i, sep="")

# Error: attempt to apply non-function

paste("fitted.model",i)<-lme(log(Y.to.use+1)  ~ x1*x2, data=example, random=~1| Exp/Person)
} 

# Error in paste("fitted.model", i, sep = "") <- lme(log(Y.to.use+1)  ~ x1*x2  : target of assignment expands to non-language object

Я создам таблицу с коэффициентами из каждой подобранной модели каждого Ивара.