Brug 'alias'-funktionen i R for at se, hvilke variabler der er lineært afhængige. Fjern de afhængige variable, og vif-funktionen skulle fungere korrekt.
formula <- as.formula(Spring_Autumn ~ Oct + Nov + Dec + Jan + Feb + Mar + Apr + May + Jun + Jul + Aug + Sep + X1min + X3min + X7min + X30min + X90min + X1max + X3max + X7max + X30max + X90max + BF + Dmin + Dmax+ LP + LPD + HP + HPD + RR + FR + Rev, data = IHA_stats)
fit <-lm(formula)
#the linearly dependent variables
ld.vars <- attributes(alias(fit)$Complete)$dimnames[[1]]
#remove the linearly dependent variables variables
formula.new <- as.formula(
paste(
paste(deparse(formula), collapse=""),
paste(ld.vars, collapse="-"),
sep="-"
)
)
#run model again
fit.new <-lm(formula.new)
vif(fit.new)
BEMÆRK:Dette vil ikke fungere i tilfælde af, at du har automatisk genererede dummy-variabler, der er identiske med andre variabler. Variabelnavnene bliver rodet. Du kan oprette dit eget hack for at komme uden om det.