On the multivariate spatial median for clustered data