Signal flanking rDNA

This function allows you to pull out the ChIP signal flanking the rDNA region on chromosome 12. The function takes as input the wiggle data as a list of 16 chromosomes. (output of readall_tab).

signal_at_rDNA(inputData, saveFile = FALSE)

Arguments

inputData

As a list of the 16 chr wiggle data (output of readall_tab). No default.
saveFile

Boolean indicating whether output should be written to a .txt file (in current working directory). If saveFile = FALSE, output is returned to screen or an R object (if assigned). Defaults to FALSE.

Value

A local data frame with four columns:

  1. chr Chromosome number

  2. position Nucleotide coordinate (in normalized total length of 1 kb)

  3. signal ChIP-seq signal at each position (1 to 1000)

  4. gene Systematic gene name

Examples

# NOT RUN {
signal_at_rDNA(WT)

signal_at_rDNA(WT, saveFile = TRUE)
# }