Open
Description
Description
It is currently possible to create dims without coords via pm.Data
:
import pymc as pm
with pm.Model() as m:
x = pm.MutableData("x", [0, 1, 2], dims=("dim_without_coord",))
y = pm.Normal("y", x.sum(), dims=("dim_without_coord",))
We used to allow this also from RVs in the past but reverted in #6112
I don't see a fundamental difference between the two cases, and in that PR I just missed this.
I would suggest we remove this functionality because it adds quite some complexity in the Model class, including fancy checks to see if a Data variable defined from dims without coords of another Data variable is being resized before the original one (but after is fine?):
Lines 923 to 937 in 02cbac6
(The unawareness) of this option has also led to bugs that were only patched recently: #6882