Disaggregate and find the change point in time series using PyMC

[matrix_testing]

Thank you! That fixed the problem. Also I had to change

import pymc3 as pm

to

 import pymc as pm

otherwise it complained that it could not find MutableData .