A novel unbiased method links variability of co-expression between multiple proteins on single cells to a clinical phenotype