I’m facing a very similar problem atm where I have a masked array in a pm.Minibatch in a theano.shared var. This could explain wrong estimates I am facing. Is there an issue for this open somewhere?
Where and how is the missing observed treatment for numpy masked arrays implemented? I would like to look into it to understand if it is my problem.