BedGC.chr: Compute GC-content of genomic regions defined by a BED file...
In gustaveroussy/EaCoN: EaCoN : Easy Copy Number !

Description Usage Arguments Details Value Author(s)

This computes the GC-content of a BED file, and possibly extensions of these regions, using (samtools-indexed) FASTA files for each chromosome.

1
2
3

EaCoN.BedGC.chr(binned.bed.file = NULL, genome = "hg19", fasta.dir = NULL, na.to0 = TRUE,
nt.add = c(0,50,100, 250, 500, 1000, 2500, 5000), out.dir = getwd(), return.data = FALSE,
nthread = 2)

`binned.bed.file`	BED file with binned region (ie, neither the WES capture BED, nor its cleaned version. See `Details`).
`genome`	Genome build (one supported by datasets in the `chromosomes` package).
`fasta.dir`	A directory containing (samtools-indexed) FASTA files corresponding the nucleotidic sequence of the given `genome` chromosomes.
`na.to0`	Replace putative NA values of the computed GC-content by zero.
`nt.add`	A numeric vector corresponding to Window(s) of extension of the BED regions (see `Details`).
`out.dir`	Output directory.
`return.data`	Return the computed GC-content as a data.frame (and still write it as a BED-like file).
`nthread`	Number of threads for multithreading (see `Details`).

binned.bed.file : This BED file corresponds to the BED generated by EaCoN::EaCoN.BAMConv(), thus corresponding to binning of the capture BED after being processed through EaCoN::EaCoN.BedCheck().
nt.add : To simulate different states of degradation of DNA, GC-content can be computed on several tracks corresponding to different sizes of extension of the regions defined by te binned.bed.file. The nt.add parameter corresponds to a numeric vector containing the different extension size(s) desired. Please note that these sizes are expressed in nucleotides unit, and that the extension will be performed on both sides of the regions. Example : with an extension size of 50, the genomic region chr1:1000-1100 will get extended to chr1:950-1150. Default value for this parameter is recommended for WES data.
nthread : When nthread > 1, multithreading is performed using one chromosome (one extension) per thread, so please do not use a value greater than the length of chromosomes in binned.bed.file.