数据分析之peakfinder

时间:2022-07-10 19:23:56 其他 我要投稿
  • 相关推荐

数据分析之peakfinder

ChIPSeq Peak Finder

程序下载地址

总体而言,因为程序都是一堆 python 脚本,写的很分散,所以感觉用

起来不是很好用,所以现在开始测试这个程序。

Peak finder 解压,数了数,一共有17 *.py 文件,也没作什么合并

所以几天都没有跑起来

I.程序文档的基本解读

1.

You will want to first convert Solexa output for the chip

and the control sample into bed files using one of the

following scripts:

maketrackfromeland.py

maketrackfromrealign.py

覆盖 Solexa 输出到 chip, 使用这两个脚本控制 示例 到 基准文件

2.

The following scripts are used to read the output from the

0.3 version of ELAND run with the --multi option:

maketrackfromeland2.py

maketrackmulti.py

下面的脚本用于读 ELAND 0.3版本的输出, 使用 --multi 选项

3.

You can also create a bed-formatted WIG file, for display

The following scripts are used to read the output from the

0.3 version of ELAND run with the --multi option:

maketrackfromeland2.py

maketrackmulti.py

你也能创建一个 基准 WIG 文件,以上的脚本用于读 ELAND 0.3 版本

的输出, 使用 --multi 选项

4.

You will want to first convert Solexa output for the chip

and the control sample into bed files using one of the

following scripts:

maketrackfromeland.py

maketrackfromrealign.py

Chip 到 Solexa 输出的转换,控制 示例 到 基准文件.

5.

on the UCSC browser:

makewiggle.py

USCE 浏览器, 这个脚本什么作用?

6.

The main script actually implements the peak finder:

findall.py

peak finder 实际执行的主脚本

7.

You will want to first convert Solexa output for the chip

and the control sample into bed files using one of the

following scripts:

maketrackfromeland.py

maketrackfromrealign.py

on the UCSC browser:

findallnocontrol.py

文件转换 和 示例 矫正 到 基准,作者推荐使用第一个脚本

8.

NEW FEATURE of findall.py : as of version 2.0, you can

/ should use the -normalize option to calculate

everything as Reads Per Million (RPM). While we have

kept the original behavior as default, we will switch

-normalize to be the default in the next release.

findall.py 脚本的新特征: version 2.0 可以使用-normalize

选项计算每个RPM(Reads Per Million). 我们默认保持原样,下

一个版本将会打开 -normalize

The philosophy of this peak finder is to define regions,

and then search for the motif. However, the findall

script can report the actual peaks in the region with

the -listpeak option.

peak finder 的哲学是定义区域, 搜索模体。尽管这样, findall

脚本报告实际的峰的区域,选项, -listpeak

9.

The rest of the analysis depends heavily on Cistematic

to run. The following scripts find associated genes and

anlyze their GO ontology enrichment, if any:

getallgenes.py

analyzego.py

基于 Cistematic 的其余分析,关联 基因 和 GO 富集

10.

The following scripts, also requiring Cistematic,

the sequence in the enriched regions, find motifs using

Meme and map motif sites in regions around the peaks:

getfasta.py

findMotifs.py

getallsites.py

其余脚本, 也要求 Cistematic, 恢复富集区域的序列,使用

MEME 寻找模体,比对peak附近的模体区域

11.

The output of findMotifs.py and input of getallsites.py

are motifs in the Cistematic .mot format. A modified

version of getallsites.py to output NRSEs that uses

multiple motifs is:

getallNRSE.py

NRSE2.mot

NRSE2left.mot

NRSE2right.mot

findMotif.py 的输出 以及 getallsites.py 的输入均是 Cistematic .mot格式。

一个修饰的版本是getallsites.py 到 NRSEs 使用 多个 模体。

12.

The remaining scripts are just helper scripts to allow

comparison between runs and/or move data into UCSC format.

bedtoregion.py

makesitetrack.py

regiontobed.py

regionintersects.py

siteintersects.py

剩余的脚本是一些帮助脚本,帮助比较运行或转换数据到UCSC格式

II. 程序测试实例.

[数据分析之peakfinder]相关文章:

1.数据分析之peakfinder

【数据分析之peakfinder】相关文章:

数据分析报告07-28

大数据分析07-20

销售数据的分析方法07-25

大数据分析07-25

多维数据分析方法04-07

数据分析工作职责07-28

数据分析常用方法07-26

数据分析主管的职责01-13

数据分析报告【推荐】03-07

数据分析简历模板02-17