How to calculate VIF in python 3.6 with scipy 1.0.0?
First, ss module is in scipy.stats not scipy.stats.stats as of scipy 0.17. Also, it is deprecated, but you seem to have figured that out.

Second, in the current linear_model.py on github, the ss package is no longer imported:

from scipy.linalg import toeplitz
from scipy import stats
from scipy import optimize


What version is your statsmodels? Sounds like you may need to update it.