需要Vif值而不是gvif

Question

如何在没有Gvif和GVIF^(1/(2*Df))的情况下获得VIF值。我已经尝试了命令vif(model)，只需要vif值，但我得到输出为gvif

Answer 1

如果您想要回归模型中的VIF值，最简单的解决方案是使用car包：

library(car)
vif(model)

返回vif值：

             gdp labour_participation           m_per1000f 
        1.100277             1.457567             1.667722 
     time_prison 
        1.247356 

如果你想手动计算VIF值（更难的方法但不使用任何库），你可以这样做来验证你从上面得到的结果是否确实正确：

vif_lp <- 1/(1-(summary(lm(labour_participation ~ gdp + m_per1000f + time_prison, crime))$r.squared))
vif_lp
# returns [1] 1.457567