Как устранить переменную в модели в R? - PullRequest
Винтажный Клуб для гитаристов
0 голосов
/ 23 декабря

Я хочу написать функцию обратного исключения в R.
Я пробовал это:

B.Elimination <- function(model, alpha.drop){
 n.of.var <- length(model$effect.sources)
 model.after.e <- model
 for(i in 1:n.of.var){
   if(max(summary(model.after.e)$coefficients[,4]) > alpha.drop){
     index <- which.max(summary(model.after.e)$coefficients[,4])
     model.after.e <- update(model.after.e, .~. - get(model.after.e$effect.sources[index]))
     }else{
       break
     }
   }
  return(model.after.e)
}

Но как-то переменная, которую я хочу исключить, это p-значение> alpha.dropне может быть устранено.
Я всегда получал оригинальную модель в конце.

...