SVC and code to reproduce

This is the working Q-MGP code for univariate modeling of spacetime data with spatially-varying intercept. y = XB + w + eps