R を用いると以下のようになる。
> ## トレランスとVIF > r <- cor(X) > VIF <- diag(solve(r)) > tolerance = 1/VIF > data.frame(tolerance, VIF) tolerance VIF X1 0.983576 1.016698 X2 0.983576 1.016698