Variational Inference for GARCH-family Models