Я не смог найти пакет или функцию для этого, но я нашел способ сделать это вручную:
using RDatasets, DataFrames, CSV, GLM
airquality = rename(dataset("datasets", "airquality"), "Solar.R" => "Solar_R")
model = lm(@formula(Wind ~ Temp+ Solar_R), airquality)
print(1/(1-r2(model)))
Это возвращает 1.267492, что совпадает с VIF(lm(Wind ~ Temp+Solar.R, data=airquality))
в R