Multivariate version of Rolling Regression Example

If you are using pymc>4, the time dimension of the MvGaussianRandomWalk (as well as the univariate) has to be the rightmost dimension