Implements variational empirical Bayes methods for linear regression. Coordinate-wise optimization algorithms for model fitting are implemented in C++ for efficient handling of large-scale data sets.