ENCODE-DCC/chip-seq-pipeline2

Can I run chip-seq-pipeline without input files?

Reemann opened this issue · 3 comments

Describe the bug

If I run without chip.ctl* set in .json, it will report:

* Error: SPP requires control inputs. Define control input files (chip.ctl_*) in an input JSON file.

I know tf mode runs with input files which is different from histone mode, but is there any possibility to run ChIP-seq-pipeline without input files? Because I don't have input datasets in my research.

OS/Platform

  • OS/Platform: Ubuntu
  • Conda version: conda 4.10.3
  • Pipeline version: [e.g. v1.6.0]
  • Caper version: [e.g. v1.2.0]

Input JSON file

{
    "chip.pipeline_type" : "tf",
    "chip.genome_tsv" : "./genome_ENCODE/hg38.tsv",
    "chip.fastqs_rep1_R1" : ["./0_raw_data/TF1_rep1_R1.fastq.gz"
    ],
    "chip.fastqs_rep1_R2" : ["./0_raw_data/TF1_rep1_R2.fastq.gz"
    ],
    "chip.paired_end" : true,
    "chip.title" : "TF1",
    "chip.description" : "TF1"
}

Troubleshooting result

Paste troubleshooting result.

* Error: SPP requires control inputs. Define control input files (chip.ctl_*) in an input JSON file.

I think you can try with macs2 peak caller. Please add the following to your input JSON and try again. This will change your peak caller from spp to macs2. spp always requires inputs but macs2 doesn't.

{
    "chip.pipeline_type": "tf",
    "chip.peak_caller": "macs2"
}

I'll try as you recommend. Thanks a lot!

Hello all,
I'm facing the same issue (also running a 'tf' chip). I was wondering whether adding macs2 (instead of spp) did help solve the issue.
Also, I wanted to know where exactly in the JSON file we should add "chip.peak_caller": "macs2" toward the end or just after "chip.pipeline_type": "tf".

Thanks.

Best,
Karim