Linear Gaussian Models in R