This function takes a snap object and peak list as input and creates cell-by-peak matrix
createPmat(obj, peaks, ncell.chunk = 20, do.par = TRUE, num.cores = 1)
A snap object.
A GRanges object contains the peaks
A numeric class that indicates the number of cells to calculate per processing core [20].
A logical variable indicates if run this using multiple processors [TRUE].
Number of processers to use [1].
Returns a snap objects contains cell-by-peak matrix