Question: trying to do MACS!!!
0
gravatar for tali.mandel
4.0 years ago by
tali.mandel10
Hebrew University
tali.mandel10 wrote:

Hello,

I have a file in BAM format of chip-seq I did in galaxy

when i try to do MACS on it I get this error message that I have no idea what it means:

INFO  @ Wed, 19 Nov 2014 04:40:52: 
# ARGUMENTS LIST:
# name = MACS_in_Galaxy
# format = BAM
# ChIP-seq file = /galaxy-repl/main/files/009/715/dataset_9715201.dat
# control file = None
# effective genome size = 9.00e+07
# tag size = 25
# band width = 300
# model fold = 6
# pvalue cutoff = 1.00e-05
# Ranges for calculating regional lambda are : peak_region,1000,5000,10000 
INFO  @ Wed, 19 Nov 2014 04:40:52: #1 read tag files... 
INFO  @ Wed, 19 Nov 2014 04:40:52: #1 read treatment tags... 
INFO  @ Wed, 19 Nov 2014 04:41:02:  1000000 
INFO  @ Wed, 19 Nov 2014 04:41:12:  2000000 
INFO  @ Wed, 19 Nov 2014 04:41:22:  3000000 
INFO  @ Wed, 19 Nov 2014 04:41:35:  4000000 
INFO  @ Wed, 19 Nov 2014 04:41:45:  5000000 
INFO  @ Wed, 19 Nov 2014 04:41:55:  6000000 
INFO  @ Wed, 19 Nov 2014 04:42:05:  7000000 
INFO  @ Wed, 19 Nov 2014 04:42:15:  8000000 
INFO  @ Wed, 19 Nov 2014 04:42:25:  9000000 
INFO  @ Wed, 19 Nov 2014 04:42:35:  10000000 
INFO  @ Wed, 19 Nov 2014 04:42:44:  11000000 
INFO  @ Wed, 19 Nov 2014 04:42:54:  12000000 
INFO  @ Wed, 19 Nov 2014 04:43:04:  13000000 
INFO  @ Wed, 19 Nov 2014 04:43:14:  14000000 
INFO  @ Wed, 19 Nov 2014 04:43:22:  15000000 
INFO  @ Wed, 19 Nov 2014 04:43:30:  16000000 
INFO  @ Wed, 19 Nov 2014 04:43:38:  17000000 
INFO  @ Wed, 19 Nov 2014 04:43:46:  18000000 
INFO  @ Wed, 19 Nov 2014 04:43:54:  19000000 
INFO  @ Wed, 19 Nov 2014 04:44:01:  20000000 
INFO  @ Wed, 19 Nov 2014 04:44:09:  21000000 
INFO  @ Wed, 19 Nov 2014 04:44:17:  22000000 
INFO  @ Wed, 19 Nov 2014 04:44:32: #1  Background Redundant rate: 0.03 
INFO  @ Wed, 19 Nov 2014 04:44:32: #1 finished! 
INFO  @ Wed, 19 Nov 2014 04:44:32: #2 Build Peak Model... 
INFO  @ Wed, 19 Nov 2014 04:44:37: #2 number of paired peaks: 1 
WARNING @ Wed, 19 Nov 2014 04:44:37: Too few paired peaks (1) so I can not build the model! Lower your MFOLD parameter may erase this error. 
WARNING @ Wed, 19 Nov 2014 04:44:37: Process is terminated! 

Can I do MACS on one sample? (with no control?)

Help please

Tali 

galaxy chip-seq • 1.1k views
ADD COMMENTlink modified 4.0 years ago by Dannon Baker3.7k • written 4.0 years ago by tali.mandel10
0
gravatar for Dannon Baker
4.0 years ago by
Dannon Baker3.7k
United States
Dannon Baker3.7k wrote:

Hi Tali,

Have you tried adjusting the MFOLD parameter, lowering it?

See the FAQ here: http://liulab.dfci.harvard.edu/MACS/FAQ.html

ADD COMMENTlink written 4.0 years ago by Dannon Baker3.7k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 175 users visited in the last hour