Implementing multilabel confusion matrix - Issue #3452

The marginals don't need to be shown, so it's possible to implement with a numpy array as long as that additional dimension signified there by `/` (it seems) is represented explicitly as an axis of size 2 (or via returning a pair of n_classes x n_classes matrices).