Variance inflation factor, tolerance, eigenvalues and condition indices.

blr_coll_diag(model) blr_vif_tol(model) blr_eigen_cindex(model)

model | An object of class |
---|

`blr_coll_diag`

returns an object of class `"blr_coll_diag"`

.
An object of class `"blr_coll_diag"`

is a list containing the
following components:

tolerance and variance inflation factors

eigen values and condition index

Collinearity implies two variables are near perfect linear combinations of
one another. Multicollinearity involves more than two variables. In the
presence of multicollinearity, regression estimates are unstable and have
high standard errors.
*Tolerance*

Percent of variance in the predictor that cannot be accounted for by other predictors.
*Variance Inflation Factor*

Variance inflation factors measure the inflation in the variances of the parameter estimates due to
collinearities that exist among the predictors. It is a measure of how much the variance of the estimated
regression coefficient \(\beta_k\) is inflated by the existence of correlation among the predictor variables
in the model. A VIF of 1 means that there is no correlation among the kth predictor and the remaining predictor
variables, and hence the variance of \(\beta_k\) is not inflated at all. The general rule of thumb is that VIFs
exceeding 4 warrant further investigation, while VIFs exceeding 10 are signs of serious multicollinearity
requiring correction.
*Condition Index*

Most multivariate statistical approaches involve decomposing a correlation matrix into linear combinations of variables. The linear combinations are chosen so that the first combination has the largest possible variance (subject to some restrictions), the second combination has the next largest variance, subject to being uncorrelated with the first, the third has the largest possible variance, subject to being uncorrelated with the first and second, and so forth. The variance of each of these linear combinations is called an eigenvalue. Collinearity is spotted by finding 2 or more variables that have large proportions of variance (.50 or more) that correspond to large condition indices. A rule of thumb is to label as large those condition indices in the range of 30 or larger.

Belsley, D. A., Kuh, E., and Welsch, R. E. (1980). Regression Diagnostics: Identifying Influential Data and Sources of Collinearity. New York: John Wiley & Sons.

# model model <- glm(honcomp ~ female + read + science, data = hsb2, family = binomial(link = 'logit')) # vif and tolerance blr_vif_tol(model)#> # A tibble: 3 x 3 #> Variable Tolerance VIF #> <chr> <dbl> <dbl> #> 1 female1 0.982 1.02 #> 2 read 0.602 1.66 #> 3 science 0.594 1.68# eigenvalues and condition indices blr_eigen_cindex(model)#> Eigenvalue Condition Index intercept female1 read science #> 1 3.57391760 1.000000 0.002062502 0.02386799 0.001675904 0.001561977 #> 2 0.39409893 3.011408 0.003037479 0.91508462 0.003889972 0.004232601 #> 3 0.01888407 13.757025 0.961793702 0.04701680 0.291423043 0.090184221 #> 4 0.01309940 16.517583 0.033106318 0.01403058 0.703011081 0.904021201# collinearity diagnostics blr_coll_diag(model)#> Tolerance and Variance Inflation Factor #> --------------------------------------- #> # A tibble: 3 x 3 #> Variable Tolerance VIF #> <chr> <dbl> <dbl> #> 1 female1 0.982 1.02 #> 2 read 0.602 1.66 #> 3 science 0.594 1.68 #> #> #> Eigenvalue and Condition Index #> ------------------------------ #> Eigenvalue Condition Index intercept female1 read science #> 1 3.57391760 1.000000 0.002062502 0.02386799 0.001675904 0.001561977 #> 2 0.39409893 3.011408 0.003037479 0.91508462 0.003889972 0.004232601 #> 3 0.01888407 13.757025 0.961793702 0.04701680 0.291423043 0.090184221 #> 4 0.01309940 16.517583 0.033106318 0.01403058 0.703011081 0.904021201