methylMnM (version 1.10.0)

countMeDIPbin: Compute the total MeDIP-seq number of each bin.

Description

The function is used to compute the total MeDIP-seq number of each bin.

Usage

countMeDIPbin (file.Medipsite,file.blacklist=NULL,file.bin=NULL,
file.CNV=NULL, writefile=NULL, reportfile=NULL, binlength=500)

Arguments

file.Medipsite

The path of MeDIP-seq site file or sequence tag file.

file.blacklist

The path of blacklist file (If we do not use the file, there will be defaulted as NULL).

file.bin

The path of all bins file. For computing the number of sequence tag of each window, we use the file as a normalization window position. (If we do not use the file, there will be defaulted as NULL).

file.CNV

If need, we should input CNV file to normalize count of each bin.

writefile

The path of output results. (If writefile=NULL, there will return the results back to main program.)

reportfile

The path of output results of bin length, the number of bin, total reads before processing and total reads after processing.

binlength

The length of each window.(Defaulted length is 500).

Value

The MeDIP-seq site should include at least three columns "chromosome", "start position" and "end position". The output file is include four columns, that is "chromosome", "start position", "end position" and "MeDIP-seq count". Also, the function output a report for some parameters.

Examples

  datafile<-system.file("extdata",  package = "methylMnM") 
  filepath<-datafile[1]
  file.Medipsite<-paste(filepath,"/all_CpGsite_chr18.txt",sep="")
  f<-countMeDIPbin(file.Medipsite, binlength=5000)