Duality-based calculations for transition probabilities in stochastic chemical reactions

An idea for evaluating transition probabilities in chemical reaction systems is proposed, which is efficient for repeated calculations with various rate constants. The idea is based on duality relations; instead of direct time-evolutions of the original reaction system, the dual process is dealt with. Usually, if one changes rate constants of the original reaction system, the direct time-evolutions should be performed again, using the new rate constants. On the other hands, only one solution of an extended dual process can be re-used to calculate the transition probabilities for various rate constant cases. The idea is demonstrated in a parameter estimation problem for the Lotka-Volterra system.