I realized that the issue might not be related to how PyTensor is used, but rather to the fact that the CPT of my combination node consists of 0 and 1, with a large number of 0. I found that in PyMC, categorical no\des defined with CPTs of 0 and 1 often face difficulties in sampling 
