cosmo confusion matrix hdr

function [confusion_matrix, classes] = cosmo_confusion_matrix(ds, varargin)