write_macs2_code: Write code for macs2 (only for auprc purpose).

In ZW-xjtlu/exomePeak2Test: Easy test of the exomePeak2 performance.

Description

Write code for macs2 (only for auprc purpose).

Usage

write_macs2_code(bam_ip, bam_input, name)

Arguments

bam_ip	a vector containing the BAM files for IP samples.
bam_input	a vector containing the BAM files for input samples.
name	a character for the name of the output sh file and peak calling file.

Value

A file named by name.

Examples

OneStep_PRC(bam_ip = c("bam_local/SRR1035214_sorted.bam",
                       "bam_local/SRR1035216_sorted.bam",
                       "bam_local/SRR1035222_sorted.bam",
                       "bam_local/SRR1035224_sorted.bam"),
            bam_input = c("bam_local/SRR1035213_sorted.bam",
                          "bam_local/SRR1035215_sorted.bam",
                          "bam_local/SRR1035221_sorted.bam",
                          "bam_local/SRR1035223_sorted.bam"),
            name = "hNPC_control")

ZW-xjtlu/exomePeak2Test documentation built on Jan. 6, 2020, 3:36 p.m.