# To detect enriched high-scoring clusters, get highest modscore that is not uniformly distributed
Rscript /data/images/proton/DKlab/mr/parclip/paralyzer/get_modescore_hist1.r
modscore_cutoff for IFN 0.957 0.1134793327234
modscore_cutoff for IL4 0.969 0.216298175051088
modscore_cutoff for 0h_all_reps 0.959 0.100461110973894
modscore_cutoff for 2h_all_reps 0.96 0.207257556532147
modscore_cutoff for 6h_all_reps 0.955 0.179819664237969



# per sample
Rscript /data/images/proton/DKlab/mr/parclip/paralyzer/get_modescore_hist1.r
modscore_cutoff for 0hrep1 0.962 0.0560802059206641
modscore_cutoff for 0hrep2 0.963 0.0526080540848064
modscore_cutoff for 0hrep3 0.961 0.120913954333128
modscore_cutoff for 2hrep1 0.965 0.139925455405083
modscore_cutoff for 2hrep2 0.964 0.0715323650379684
modscore_cutoff for 2hrep3 0.962 0.149587618299848
modscore_cutoff for 6hrep1 0.966 0.0579317957861902
modscore_cutoff for 6hrep2 0.967 0.0800039136567543
modscore_cutoff for 6hrep3 0.962 0.0930348364107952