Mixture with multiple observations

handling shape in pm.NormalMixture might be tricky, maybe try reshaping the output data into a vector (also the mu and sd for the mixture)