Before revision:   In this paper we introduce a novel algorithm for this problem that makes use of sampling in order to handle large datasets.
After revision:  Previous methods have addressed this issue heuristically or by focusing on acyclic networks and specific classes of regulation functions. In this paper we introduce a novel algorithm for this problem that makes use of sampling in order to handle large datasets.