nLogLik {rxSeq}

Negative log likelihood for coefficients provided in results of the fit using process function

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)