Control of initialization bias in multivariate simulation response