Doubly structured sparsity for grouped multivariate responses with application to functional outcome score modeling