getNullCount {mseq} | R Documentation |
Get the counts of reads for each position under the null hypothesis (all reads have the same sequencing preferences).
getNullCount(data)
data |
the data frame generated by expData.R |
a numeric vector stores the counts
# read and expand the data data(g1_part) # for real data, please use read.csv, like g1 <- read.csv("g1.csv") data <- expData(g1_part, 2, 3) null.count <- getNullCount(data) null.dev <- getDev(null.count, data$count)