Many biomolecular systems can be described by chemical reaction networks. Determining which chemical reaction network models are inconsistent with observed data can be done via model invalidation. In this work, we formulate and solve a robust version of the model invalidation problem for the case where only measurements from the stationary distribution are available. This problem corresponds to determining if an observed distribution could have been generated by the given chemical reaction network for some value of the parameters, plus a perturbation of bounded size with respect to total variation distance. The main technical tool we introduce to solve the problem is a set of generalized moments that make the problem amenable to an algorithmic solution.
