TMD！！！

前言

最近在整理WGBS分析的流程，下游需要找 Differentially Methylated Loci (DML) / Region (DMR)，类似普通转录组中的差异分析。之前看的一篇文章提到一个R package ： DSS，看Bioconductor 的教程：近期更新的，写的还不错，所以我就试了一下，结果就是 196 cores的服务器直接卡死

196 cores 的服务器直接卡死

DSS github 地址：https://github.com/haowulab/DSS

??? 代码这么写？196 个 core，直接用193个？然后R中并行，直接内存爆炸！！正常人的思维难道不是默认1，根据使用者需求而提升么？还有12个star ？

再看看另一个工具：methylKit

这才是正常思维好吧！！

反思

当然，事情的发生也怪我，在使用之前没有认认真真的、仔仔细细的、好好的看他的文档和代码，TMD！！！