mv function

Virtual factor for multivariate responses