If I use a series of operations like
array = pm.Normal( ... )
array_reshaped = array[:, np.newaxis]
and then go on to use array_reshaped
later in the model, is that expected to cause problems for PyMC/PyTensor?
In other words, when the documentation for PyTensor says that it has tight integration with numpy, does that extend to these kinds of operations? Or do I have to use the PyTensor equivalent? And if so, what would it be?
Apologies if this isn’t clear, and thanks in advance for any guidance!