I have been trying to use a GP as an input for a function which is used for the likelihood. Some pseudo code would be:
f=gp.prior('f',X=X) func=func_call(f,...) likelihood = pm.Normal('obs', mu=func, sigma=sd, observed=obs)
But the function only take in a 4-D array so i have reshape f since gp.prior only takes in arrays of the form (length(array),1). Luckily that is possible to do inside of the model
f = f.reshape((1,1,len(X),1))
and the sampler does seem to accept this. My question is does pymc “know” what is going on with the gp prior after i reshape? Will it still sample f (f_rotated) after i have reshaped it?