Ошибка «Объект не найден» при передаче формулы модели другой функции

Этот вопрос был перенесен из Перекрестная проверка, потому что на нее можно ответить при переполнении стека. Миграция 8 лет назад .

У меня странная проблема с R, которую я не могу решить.

Я ' я пытался написать функцию, которая выполняет K-кратную перекрестную проверку для модели, выбранной пошаговой процедурой в R. (я знаю проблемы с пошаговыми процедурами, это чисто для целей сравнения) :)

Теперь проблема состоит в том, что если я определяю параметры функции (linmod, k, direction) и запускаю содержимое функции, она работает безупречно. НО, если я запускаю его как функцию, я получаю сообщение об ошибке, что объект datas.train не может быть найден.

Я пробовал выполнить функцию с помощью debug (), и объект явно существует, но R говорит, что нет, когда я действительно запускаю функцию. Если я просто подбираю модель с помощью lm (), она работает нормально, поэтому я считаю, что это проблема с пошаговой функцией в цикле, находящейся внутри функции. (попробуйте закомментировать пошаговую команду и установить прогнозы на основе обычной линейной модели.)

#CREATE A LINEAR MODEL TO TEST FUNCTION
lm.cars <- lm(mpg~.,data=mtcars,x=TRUE,y=TRUE)


#THE FUNCTION
cv.step <- function(linmod,k=10,direction="both"){
  response <- linmod$y
  dmatrix <- linmod$x
  n <- length(response)
  datas <- linmod$model
  form <- formula(linmod$call)

  # generate indices for cross validation
  rar <- n/k
  xval.idx <- list()
  s <- sample(1:n, n) # permutation of 1:n
  for (i in 1:k) {
    xval.idx[[i]] <- s[(ceiling(rar*(i-1))+1):(ceiling(rar*i))]
  }

  #error calculation
  errors <- R2 <- 0

  for (j in 1:k){
     datas.test <- datas[xval.idx[[j]],]
       datas.train <- datas[-xval.idx[[j]],]
       test.idx <- xval.idx[[j]]

       #THE MODELS+
       lm.1 <- lm(form,data= datas.train)
       lm.step <- step(lm.1,direction=direction,trace=0)

      step.pred <- predict(lm.step,newdata= datas.test)
        step.error <- sum((step.pred-response[test.idx])^2)
        errors[j] <- step.error/length(response[test.idx])

        SS.tot <- sum((response[test.idx] - mean(response[test.idx]))^2)
        R2[j] <- 1 - step.error/SS.tot
  }

  CVerror <- sum(errors)/k
  CV.R2 <-  sum(R2)/k

  res <- list()
  res$CV.error <- CVerror
  res$CV.R2 <- CV.R2

return(res)
}


#TESTING OUT THE FUNCTION
cv.step(lm.cars)

Есть мысли?

r formula

задан joran 22 November 2011 в 03:44

0 ответов

Другие вопросы по тегам:

r formula

Ошибка «Объект не найден» при передаче формулы модели другой функции

0 ответов

Похожие вопросы: