Thanks for your quick help Junpeng! I’ve now raised an issue: Multiple independent random walks with differing innovation sds · Issue #4409 · pymc-devs/pymc3 · GitHub . That diagonal normal version should be a good workaround. I didn’t realise we can just use theano to manipulate tensors like that, that’s really great.