Theano error with element wise multiplication

This line works:
radon_est = (b[county_idx] * data_set).sum(axis=1)

So, using a theano function is probably just not the right way.