nLogLik: Negative log likelihood for coefficients provided in results...

In rxSeq: Combined Total and Allele Specific Reads Sequencing Study

Description

Calculates negative log(likelihood) of an X chromosome joint TReC and ASE counts model at a given set of parameters

Usage

nLogLik(res, rc, genei, hessian=FALSE)

Arguments

res	result object from process function
rc	Read count data object created by readCounts function
genei	get results for i'th gene
hessian	a logical option whether to calculate a Hessian matrix, the default values is set to FALSE.

Value

output - list(nll=-log.likelihood,hessian=hessian matrix)

Author(s)

Vasyl Zhabotynsky [e-mail]

See Also

process, rcA, readCounts.

Examples

## Not run: 
# get negative-log likelihood at the given point:
nLogLik(res=trecase.X.out, rc=rcX, genei=1, hessian=TRUE)

## End(Not run)

rxSeq documentation built on May 1, 2019, 11:31 p.m.