TAD方法

1.DomainCaller 【任兵-densenty】

http://yuelab.org/hi-c/

2.insulationscore 【dekker- densenty 】

https://github.com/dekkerlab/crane-nature-2015

# 安裝:ok

perl scripts/matrix2insulation.pl

perl?/home/shen/biosoft/crane-nature-2015-master/scripts/matrix2insulation.pl -i xxx.matrix?-is 500000 -ids 200000 -im mean -bmoe 3 -nt 0.1 -v

-b????????[] size (bp) of the insulation square

-v????? ? ?[]? ? ? ? FLAG, verbose mode

-ids? ? ? ?[]? ? ? ? insulation delta span (size (bp) of insulation delta window)

-im? ? ? ? ?[]? ? ? ? insulation mode (how to aggregrate signal within insulation square), mean,sum,median

-nt? ? ? ? ?[0.1]? ? ? noise threshold, minimum depth of valley

-bmoe? ? ? [3]? ? ? ? boundary margin of error (specified in number of BINS), added to each side of the boundary

軟件給的例子:

perl scripts/matrix2insulation.pl -i 10kb_chrX.matrix.gz -is 500000 -ids 200000 -im mean -bmoe 3 -nt 0.1 -v

perl scripts/matrix2insulation.pl -i 10kb_chrX.matrix.gz -is 500000 -ids 200000 -im mean -bmoe 3 -nt 0.1 -v

3.juicer:[.hic]

ref:https://github.com/aidenlab/juicer/blob/master/README.md

4.hicexplorer: 【.h5】

ref:?https://hicexplorer.readthedocs.io/en/latest/content/tools/hicFindTADs.html

5.tadtool:[可以對(duì)接hicpro輸出]

ref:https://github.com/vaquerizaslab/tadtool

pip intall tadtools

tadtools -h

Commands:

? ? plot? ? ? ? ? ? ? ? Main interactive TADtool plotting window

? ? tads? ? ? ? ? ? ? ? Call TADs with pre-defined parameters

? ? subset? ? ? ? ? ? ? Reduce a matrix to a smaller region.

my:

tadtool tads iced/50000/TU-2_50000_iced.matrix00_abs.bed 50000 0.005 tadtools_output.bed

tadtool plot iced/50000/TU-2_50000_iced.matrix raw/50000/TU-2_50000_abs.bed chr5:34000000-37000000

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖棚亩,帶你破解...
    沈念sama閱讀 219,539評(píng)論 6 508
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異苟鸯,居然都是意外死亡面睛,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,594評(píng)論 3 396
  • 文/潘曉璐 我一進(jìn)店門该镣,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)冻璃,“玉大人响谓,你說(shuō)我怎么就攤上這事损合。” “怎么了娘纷?”我有些...
    開封第一講書人閱讀 165,871評(píng)論 0 356
  • 文/不壞的土叔 我叫張陵嫁审,是天一觀的道長(zhǎng)。 經(jīng)常有香客問我赖晶,道長(zhǎng)律适,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,963評(píng)論 1 295
  • 正文 為了忘掉前任遏插,我火速辦了婚禮捂贿,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘钧忽。我一直安慰自己茎刚,他們只是感情好怜校,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,984評(píng)論 6 393
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著颜屠,像睡著了一般。 火紅的嫁衣襯著肌膚如雪鹰祸。 梳的紋絲不亂的頭發(fā)上甫窟,一...
    開封第一講書人閱讀 51,763評(píng)論 1 307
  • 那天,我揣著相機(jī)與錄音蛙婴,去河邊找鬼粗井。 笑死,一個(gè)胖子當(dāng)著我的面吹牛街图,可吹牛的內(nèi)容都是我干的浇衬。 我是一名探鬼主播,決...
    沈念sama閱讀 40,468評(píng)論 3 420
  • 文/蒼蘭香墨 我猛地睜開眼台夺,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼径玖!你這毒婦竟也來(lái)了?” 一聲冷哼從身側(cè)響起颤介,我...
    開封第一講書人閱讀 39,357評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤梳星,失蹤者是張志新(化名)和其女友劉穎赞赖,沒想到半個(gè)月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體冤灾,經(jīng)...
    沈念sama閱讀 45,850評(píng)論 1 317
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡前域,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,002評(píng)論 3 338
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了韵吨。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片匿垄。...
    茶點(diǎn)故事閱讀 40,144評(píng)論 1 351
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖归粉,靈堂內(nèi)的尸體忽然破棺而出椿疗,到底是詐尸還是另有隱情,我是刑警寧澤糠悼,帶...
    沈念sama閱讀 35,823評(píng)論 5 346
  • 正文 年R本政府宣布届榄,位于F島的核電站,受9級(jí)特大地震影響倔喂,放射性物質(zhì)發(fā)生泄漏铝条。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,483評(píng)論 3 331
  • 文/蒙蒙 一席噩、第九天 我趴在偏房一處隱蔽的房頂上張望班缰。 院中可真熱鬧,春花似錦悼枢、人聲如沸埠忘。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,026評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)给梅。三九已至,卻和暖如春双揪,著一層夾襖步出監(jiān)牢的瞬間动羽,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,150評(píng)論 1 272
  • 我被黑心中介騙來(lái)泰國(guó)打工渔期, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留运吓,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 48,415評(píng)論 3 373
  • 正文 我出身青樓疯趟,卻偏偏與公主長(zhǎng)得像拘哨,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子信峻,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,092評(píng)論 2 355