ISOpureS1.model_optimize.kappa.kappa_deriv_loglikelihood {ISOpureR}R Documentation

Compute derivative of loglikelihood with respect to kappa for step 1

Description

Computes the derivative of the part of the loglikelihood function relevant to optimizing kappa for step 1. Instead of performing constrained optimization on kappa directly, we optimize the log of kappa in an unconstrained fashion. Thus, if y=log(kappa) and L is the loglikelihood function w.r.t. y, to optimize L w.r.t. y, dL/dy = dL/dkappa * dkappa/dy, where dkappa/dy = exp(y)= exp( log(kappa)). The input into the derivative function is log(kappa - model\$MIN_KAPPA).

Usage

ISOpureS1.model_optimize.kappa.kappa_deriv_loglikelihood(log_kappa, tumordata, model)

Arguments

log_kappa

the scalar log(kappa - model\$MIN_KAPPA)

tumordata

a GxD matrix representing gene expression profiles of tumour samples

model

list containing all the parameters to be optimized

Value

The negative derivative of the part of the loglikelihood function relevant to kappa with respect to log kappa (a scalar given that for step 1 of ISOpure kappa is a scalar)

Author(s)

Gerald Quon, Catalina Anghel, Francis Nguyen


[Package ISOpureR version 1.1.3 Index]