Question: trying to do MACS!!!
0
gravatar for tali.mandel
4.0 years ago by
tali.mandel • 10
Hebrew University
tali.mandel • 10 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 COMMENT • link • modified 4.0 years ago by Dannon Baker ♦ 3.7k • written 4.0 years ago by tali.mandel • 10
0
gravatar for Dannon Baker
4.0 years ago by
Dannon Baker ♦ 3.7k
United States
Dannon Baker ♦ 3.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 COMMENT • link written 4.0 years ago by Dannon Baker ♦ 3.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