rxSeq (version 0.99.3)

get.tausB: Produce Xce estimates for mice with allele specific reads

Description

Xce estimation for mice with allele specific reads.

Usage

get.tausB(n, n0B, geneids, min.cnt=50, exclude.prop=.05, Xist.ID="ENSMUSG00000086503")

Arguments

n

vector of allele specific counts for each mouse

n0B

vector of allele specific counts for allele B

geneids

gene IDs

min.cnt

minimum number of allele specific counts

exclude.prop

minimum proportion of allele specific counts for each allele

Xist.ID

and ID of Xist, to exclude it from estimating Xce, since Xce would 1-tausB

Value

output - matrix of 4 rows:each column represent respective mouse.

See Also

process,nLogLik, data.X, rcX.

Examples

# Estimating XCE effect for each mouse for X chromosome
get.tausB(n=data.X$n, n0B=data.X$n0B, geneids=data.X$geneids)