hmm, in pymc3 we could subtract tensor variables from the observed to isolate the cause of x on y given x has some relationship to an unobserved Z the mu is supposed to control for this . I have ran with just the data in y but im unsure if this is mathematically sound for controlling the confounder.