cpg<-matrix( 1:800, nrow=400 )
cpg[,2]<-cpg[,1]+2
data<-cpg[3:100,]
data[,1]<-data[,1]-1
data[,2]<-data[,2]+1
f <- cpgcount(data[,1], data[,2], cpg[,1], cpg[,2], length(data[,1]),
length(cpg[,2]), count=rep(0,length(cpg[,2])))
Run the code above in your browser using DataLab