Hi @Jev,
Using list comprehension to create pymc3 variables are highly inefficient, if it is necessary you can do newvar = tt.stack([list]) to compress the list into a theano RV.
It is a bit unclear to me what the observation would be in your model, is it [K, I]? Also could you please detail the shape of the variables on the 5 steps above?