Bian Li

Thanks for the reply. My PyMC and aesara version are ``` pymc 4.1.3 aesara 2.7.7 ```

Changing ```python logp_fn = m.logp logp_censored_fn = m_censored.logp ``` to ```python logp_fn = m.compile_logp() logp_censored_fn = m_censored.compile_logp() ``` did resolve the error. Thank you so much! I typed up everything...

Is it likely that it is because the two pipelines are using different default idr thresholds? The WDL-based uses 0.05 whereas this one uses 0.1?