=====methylKit差异DMR分析=====

methylkit：用于高通量亚硫酸氢盐测序的DNA甲基化分析和注释的R包。该包的目的是处理RRBS以及WGBS的测序数据。此外，也可以处理从Tab-seq或oxBS-seq获得的5hmC的碱基对分辨率数据。

DSS：该软件包专为高通量测序数据的差异分析而设计，依赖于bsseq包。它提供了用于分析RNA-seq 差异表达和亚硫酸氢盐测序 (BS-seq) 差异甲基化的功能。DSS 的核心是一个基于分层贝叶斯模型的程序，用于估计和缩小基因或 CpG 位点特异性分散，然后进行 Wald 检验以检测差异表达/甲基化。

dmrseq：这个包建立在bsseq包的基础上，bsseq包提供了亚硫酸氢盐测序数据的有效存储和操作以及对差异甲基化cpg的推断。dmrseq的主要目的是检测差异甲基化区域或CpGs。

^功能	|methylkit	|DSS	|dmrseq|
^生物学重复	|不需要	|不需要	|需要|
^筛选DML/DMR	|fisher精确回归/逻辑回归	|wald检验	|GLS回归模型|
^DML	|✔	|✔	|×|
^DMR	|✔	|✔	|✔|
^注释	|✔	|×	|×|
^可视化	|✔	|✔	|✔|

使用methylKit分析过程如下：
  1、/TJPROJ6/RNA_SH/personal_dir/zhangxin/jiaoben/jiajihua/diff_DMR/methylKit/Diff_dmr_prepare.py  #拆分cx report数据
  usage: Diff_dmr_prepare.py [opthions] [value]
  
  This program is used to prepare methylKit or DSS infile
  
  optional arguments:
    -h, --help            show this help message and exit
    -CX , --cx_repot      the cx report file
    -F , --fai            the fai file
    -N SPLIT_NUMBER, --split_number SPLIT_NUMBER
                        the cx repot split number, default:20
    -MD MIN_DEPTH, --min_depth MIN_DEPTH
                        the min reads depth for C site, default:4
    -MMC MIN_MC, --min_mC MIN_MC
                        the min mC reads depth for C site, default:1
    -OF {DSS,methylKit}, --outformate {DSS,methylKit}
                        the out file formate, DSS prepare suggest to use
                        script in pipline
    -S , --sample         the sample name
    -O , --out            the out dir
  2、/TJPROJ6/RNA_SH/personal_dir/zhangxin/jiaoben/jiajihua/diff_DMR/methylKit/product_diff.py  #生成差异分析脚本
    optional arguments:
    -h, --help            show this help message and exit
    -SD SPLIT_CX_DIR, --split_cx_dir SPLIT_CX_DIR
                        the split cx file dir
    -C CONDITION, --condition CONDITION
                        the condition file
    -CP COMPARE, --compare COMPARE
                        the compare
    -D DIFFERENCE, --difference DIFFERENCE
                        the difference , default:25
    -PJ PADJ, --padj PADJ
                        the padj , default:0.05
    -PV PVALUE, --pvalue PVALUE
                        the pvalue
    -A ANALYSIS, --analysis ANALYSIS
                        dml or dmr, default:dmr
    -O OUTDIR, --outdir OUTDIR
                        the out dir
  
                        
 
DMC的脚本使用：/TJPROJ6/RNA_SH/personal_dir/fengjie/Personal_analysis/methylKit/methylKit.R  
                     
                     
  3、/TJPROJ6/RNA_SH/personal_dir/zhangxin/jiaoben/jiajihua/diff_DMR/methylKit/merge_methylKit.py  #合并拆分分析结果
  This program is used to merge methylKit result
  
  optional arguments:
    -h, --help            show this help message and exit
    -I INDIR, --indir INDIR
                        the methylKit.R analysis out dir
    -O OUT, --out OUT     the outdir