Missing Values and the Multidimensional Indexing FutureWarning

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