Thank you very much for the reply. Yes, it is a bit similar to pm.deterministic(), in that keyvalue is completely dependent on the values of m and c through my external program. Syntax-wise I am a bit uncertain on how to incorporate keyvalue into the pymc workflow (e.g. so that I can include keyvalue in pm.traceplot) through the theano Op. Any advice on this matter would be greatly appreciated.