Тестирование на мультиколлинеарность ковариат для моделей занятости в R

Question

Я пытаюсь проверить мультиколлинеарность среди моих ковариат в модели занятости в R. Я использую пакет unmarked для этих моделей занятости. Я пытался использовать VIF, но он не поддерживается в этом пакете.

Кто-нибудь может знать, как я могу проверить мультиколлинеарность для моделей занятости?

У меня 26 ковариат, поэтому было бы полезно узнать, смогу ли я устранить некоторые из них.

Ниже мой ум для этого:

siteCovs <- data.frame(
  SEff=SEff, Undis=Undis, Dis=Dis, Undis_=Undis_, 
  Dis_=Dis_, NDVI=NDVI, WS=WS, WS_=WS_, HS=HS, 
  HS_=HS_, Rdist=Rdist, Tdist=Tdist, N17dist=N17dist, 
  R550dist=R550dist, R29dist=R29dist, R548dist=R548dist, 
  LZKGdist=LZKGdist, Fdist=Fdist, R=R, T=T, N17=N17, 
  R550=R550, R29=R29, R548=R548, LZKG=LZKG, F=F
  )
umf <- unmarkedFrameOccu(
  y=ydat, 
  siteCovs= siteCovs
  )

Спасибо!

Answer 1

Это может быть проблема именования: в пакете есть функция unmarkedтакже как и vifфункция в упаковке car.

car::vifэто то, что вы использовали бы в линейной модели формы lm(), тогда как unmarked::vifэто то, что вы использовали бы в своей модели занятости. Чтобы рассчитать VIFна занятость, вы бы вошли

      unmarked::vif(mod = yourmodelhere, type = "state")

Answer 2

Вы можете проверить мультиколлинеарность с помощью пакета olsrr. Пример ниже.

library(olsrr)
model <- lm(mpg ~ disp + hp + wt + qsec, data = mtcars)
ols_coll_diag(model)