I am finding this as well. It seems that pymc struggles with missing values in multidimensional distribution arrays.