methylMnM (version 1.10.0)

calculatecount: Call C programs to R.

Description

Call C programs to R for calculate MeDIP-seq or CpG count of each bin.

Usage

calculatecount(data2, data3, cpg2, cpg3, datalength, cpglength, 
count=rep(0,cpglength))

Arguments

data2

Start position of each tag.

data3

End position of each tag.

cpg2

Start position of each bin.

cpg3

End position of each bin.

datalength

The number of tags

cpglength

The number of bins

count

Read count of each bin.

Value

Read count of each bin.

Examples

  data<-matrix( 1:800, nrow=400 )
  data[,2]<-data[,1]+37
  cpg<-matrix( 1:20, nrow=10)
  cpg[,1]<-seq(0,360,length=10)
  cpg[,2]<-seq(40,400,length=10)
  f <- calculatecount(data[,1], data[,2], cpg[,1], cpg[,2], length(data[,1]),
 length(cpg[,2]), count=rep(0,length(cpg[,2])))