Nice, this worked out exactly. I would never have guessed this; from what I recall from linear algebra, column vectors cannot be subtracted from row vectors (and vice versa.)
Could you link a resource on “broadcasting”? Seems like something to master early on for good PyMC3 modeling.