evalLogLikImpulse: Cost function for impulse model

View source: R/srcImpulseDE2_CostFunctionsFit.R

evalLogLikImpulseR Documentation

Cost function for impulse model

Description

Log likelihood cost function for numerical optimisation of impulse model. Implements log linker function for the amplitude parameters and the batch correction factors. Implements upper and lower sensitivity bound of likelihood with respect to batch correction factors and lower bound for amplitude paramters.

Usage

evalLogLikImpulse(vecTheta, vecCounts, scaDisp, vecSizeFactors,
  vecTimepointsUnique, vecidxTimepoint, lsvecidxBatch, vecboolObserved)

Arguments

vecTheta

(numeric vector number of parameters to be estimated) Impulse model parameter and batch correction factor estimates.

vecCounts

(numeric vector number of samples) Read count data.

scaDisp

(scalar) Gene-wise negative binomial dispersion hyper-parameter.

vecSizeFactors

(numeric vector number of samples) Model scaling factors for each sample which take sequencing depth into account (size factors).

vecTimepointsUnique

(numeric vector length number of unique time points) Unique time points of set of time points of given samples.

vecidxTimepoint

(index vector length number of samples) Index of of time point assigned to each sample in vector vecTimepointsUnique.

lsvecidxBatch

(list length number of confounding variables) List of index vectors. One vector per confounding variable. Each vector has one entry per sample with the index batch within the given confounding variable of the given sample. Batches are enumerated from 1 to number of batches.

vecboolObserved

(bool vector number of samples) Whether sample is observed (finite and not NA).

Value

scaLogLik (scalar) Value of cost function (loglikelihood) for given gene.

Author(s)

David Sebastian Fischer

See Also

Compiled version: evalLogLikImpulse_comp


YosefLab/ImpulseDE2 documentation built on Sept. 17, 2022, 2:45 a.m.