We present a multiresponse estimation procedure for parameters of systems described by ⋅ = As. The procedure features a generalized Gauss–Newton algorithm for optimizing the determinant criterion, efficient evaluation of the expectation function and its derivatives directly from the reaction network or compartment diagram, automatic determination of starting values, and eficient computational procedures for handling linear constraints on the responses.