Multicollinearity in WLS regression

I run a weighted-least squares regression to account for the heteroscedasticity in my data. When I examine the Pearson correlations between all predictor variables, I can't detect high collinearity. However, when I apply WLS with the inverse of the variance as my weights and then screen the variance inflation factors of the regression, I find VIFs far above 10.

My question is: Does WLS always suffer from multicollinearity in my case as the weighting with the inverse of the variance necessarily leads to higher dependencies of the predictors?

• Just to clarify: why should weighing with the inverse of squared residuals generate dependence among the regressors? The off-diagonal elements are all still 0. Also, is there any good reason you are using WLS as opposed to OLS with robust standard errors? – ChinG Oct 30 '15 at 19:10
• How many variables do you have? What are the VIFs if you fit an OLS model ignoring the heteroscedasticity? – gung - Reinstate Monica Oct 30 '15 at 20:15
• @ChinG Off-diagonal elements of what? Dependence among the regressors would naturally be evaluated by means of their covariance matrix (or, in the case of WLS, their weighted covariance matrix). The only times those off-diagonal elements are zero is when all regressors are orthogonal (which is rare except for designed experiments with no missing values). – whuber Oct 30 '15 at 20:27
• Off diagonal elements of the variance-covariance matrix. For WLS, we set the off diagonal elements to 0. If one wants dependence (autocorrelation, clustering etc) across error terms, one could implement a more general GLS estimator. – ChinG Nov 1 '15 at 14:17

It's a great question, because it concerns an important and fundamental issue about multiple regression. The answer is that weighting can change the VIF (a measure of collinearity among the regressors) by arbitrary amounts up or down (with a lower limit of $1$).
To see why that is, consider a model with two regressors $X_1$ and $X_2$. (We do not have to consider the response variable: it plays no role in computing the VIF.) Here is a scatterplot of 50 observations, with the symbol areas proportional to the weights:
If those weights were not applied, clearly the correlation would be strongly positive, because the points line up closely along a positively sloping line. Indeed, their (unweighted) VIF is $13.54$: large enough to drive one to investigate this model carefully for effects of collinearity. (Rules of thumb assert that VIFs above $5$ or $10$ begin to be of concern.)
The weights, though, throw the correlation in altogether a different direction. The points with the heaviest weights tend to be negatively correlated: some towards the upper left, others toward the lower right. These cause the weighted correlation to be nearly zero. Indeed, the VIF for these weighted data is merely $1.44$: low and benign.