Oops, accidentally overwrote this script.
But I used category=1 for hyper and category=2 for hypo.
Data came from /ddn/gs1/home/grimmsa/ruifang/WGBS_analysis/FIX/final_DMR_set/*.bed.