getNullCount {mseq}R Documentation

counts under null model

Description

Get the counts of reads for each position under the null hypothesis (all reads have the same sequencing preferences).

Usage

getNullCount(data)

Arguments

data the data frame generated by expData.R

Value

a numeric vector stores the counts

Examples

 # 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)

[Package mseq version 1.1 Index]