Defining random for a CustomDist time series

[ricardoV94]

The biggest change is that automatic imputation of multivariate distributions (which all timeseries are) is not supported in pymc>=4.

Everything else should work the same.