Estimation of multiresponse simulation metamodels using control variates