Пользовательская функция, основанная на step(), дает: Ошибка в eval (expr, envir, enclos): объект «X» не найден

Я пытаюсь создать функцию, которая охватывает функцию ordiR2step() из пакета vegan. Эта функция основана на функции step().

Это код, который отлично работает вне функции:

    install.packages("vegan")
    require(vegan)
    data(mite) 
    data(mite.env) #explanatory variables
    mite.hel = decostand(mite, "hel") #response variable
    mod0 <- rda(mite.hel ~ 1, mite.env)  # Model with intercept only 
    mod1 <- rda(mite.hel ~ ., mite.env)  # Model with all explanatory variables
    step.res <- ordiR2step(mod0, scope = formula(mod1), direction="forward") 
    step.res$anova 

Однако, если я попытаюсь обернуть это в функцию:

    forward <- function(X, Y) {
    intercept <- rda(X ~ 1, data = Y)  # Model with intercept only
    model <- rda(X ~ ., data = Y)  # Model with all explanatory variables
    # this is where debugging is stuck
    forStep <- ordiR2step(object=intercept, scope = formula(model), 
                    direction = "forward", trace = FALSE)
    list(forStep$anova)
    }
    forward(mite.hel, mite)

Я получаю следующую ошибку: Error in eval(expr, envir, enclos) : object 'X' not found

Я знаю, что при попытке обернуть эти типы функций могут возникнуть проблемы, и эти проблемы много раз обсуждались в stackoverflow и других местах, таких как здесь, здесь и здесь. Однако ни одно из этих решений, похоже, не сработало должным образом.

Насколько я понимаю это поведение, функция step() и, соответственно, ordiR2step() не могут читать из среды, определенной внутри функции, а только из среды рабочей области, потому что, если я определяю X перед вызовом функции, все работает хорошо. Однако это противоречит цели, поэтому я попробовал некоторые предлагаемые решения, такие как:

    forward2 <- function(X, Y) {
    intercept <- do.call("rda",list(X ~ 1, data = Y))  
    model <- do.call("rda", list(X ~ ., data = Y))  
    forStep <- ordiR2step(object=intercept, scope = formula(model), 
                    direction = "forward", trace = FALSE)
    list(forStep$anova)
    }
    forward2(mite.hel, mite.env)

Та же ошибка... Нет игральных костей... Есть предложения? Спасибо за ваше время!