I'm trying to run CEAS on two different sets of genes that are present in my ChIP-seq, but it seems as though I'm doing something incorrectly. When I run CEAS without the --gn-groups parameter, I get great results. I see reasonable looking profiles in the wiggle section of the .pdf output.
However, when I run CEAS with the --gn-groups parameter and two lists of genes, the average profiles don't look right at all. Both gene groups are simply graphed as perfectly vertical lines. This is definitely not correct, they should have a broad looking profile.
Here's my command line:
ceas --name CEAS_TSGTDG --gn-groups=TSG_list.txt,TDG_list.txt --gn-group-names='TSG,TDG' -g ../../../hg19.refGene -b MACSpeaks.bed -w MACS2_treat_pileup.wig
I know that my bed and wig files are good, because I get good output without supplying --gn-groups. I'm wondering if I'm not using the proper format for my "TSG_list.txt" and "TDG_list.txt"? They're simply a list of gene names followed by newlines. Here's an excerpt of TDG_list.txt:
BZRAP1
PEX13
NEK4
RNASEH2B
C12orf57
Am I doing something wrong?
However, when I run CEAS with the --gn-groups parameter and two lists of genes, the average profiles don't look right at all. Both gene groups are simply graphed as perfectly vertical lines. This is definitely not correct, they should have a broad looking profile.
Here's my command line:
ceas --name CEAS_TSGTDG --gn-groups=TSG_list.txt,TDG_list.txt --gn-group-names='TSG,TDG' -g ../../../hg19.refGene -b MACSpeaks.bed -w MACS2_treat_pileup.wig
I know that my bed and wig files are good, because I get good output without supplying --gn-groups. I'm wondering if I'm not using the proper format for my "TSG_list.txt" and "TDG_list.txt"? They're simply a list of gene names followed by newlines. Here's an excerpt of TDG_list.txt:
BZRAP1
PEX13
NEK4
RNASEH2B
C12orf57
Am I doing something wrong?
Comment