Unconstrained DeCovarT log-likelihood

Evaluates the conditional log-likelihood \(\ell_{\boldsymbol{y}\,|\,\boldsymbol{\zeta}}(\boldsymbol{p})\) of a bulk profile under the Gaussian convolution model of the article, $$ \boldsymbol{y}\,|\,(\boldsymbol{\zeta},\boldsymbol{p}) \sim\mathcal{N}_{G}\!\bigl(\boldsymbol{\mu}\boldsymbol{p},\, \boldsymbol{\Sigma}(\boldsymbol{p})\bigr), $$ with plug-in parameters \(\boldsymbol{\zeta}=(\boldsymbol{\mu},\{\boldsymbol{\Sigma}_j\}_{j=1}^{J})\) and mixture covariance \(\boldsymbol{\Sigma}(\boldsymbol{p})=\sum_{j}p_j^{2}\boldsymbol{\Sigma}_j\).

Usage

loglik_multivariate(p, y, mean_signature_matrix, Sigma)

Arguments

p

Numeric vector \(\boldsymbol{p}\in\mathbb{R}^{J}\).

y

Numeric vector (or one-column matrix) \(\boldsymbol{y}\in\mathbb{R}^{G}\).

mean_signature_matrix

Numeric matrix \(\boldsymbol{\mu}\in\mathcal{M}_{G\times J}\) (plug-in means).

Sigma

Array of cell-type covariances in \(\mathcal{M}_{G\times G\times J}\).

Value

Scalar log-likelihood value.

Details

Up to an additive constant independent of \(\boldsymbol{p}\), $$ \ell_{\boldsymbol{y}\,|\,\boldsymbol{\zeta}}(\boldsymbol{p}) = -\log\det\boldsymbol{\Sigma}(\boldsymbol{p}) -\tfrac{1}{2} (\boldsymbol{y}-\boldsymbol{\mu}\boldsymbol{p})^{\mathsf{T}} \boldsymbol{\Sigma}(\boldsymbol{p})^{-1} (\boldsymbol{y}-\boldsymbol{\mu}\boldsymbol{p}). $$ Argument mean_signature_matrix stores the plug-in mean signature \(\boldsymbol{\mu}\). Latent sample-specific profiles \(\boldsymbol{x}_{\cdot j}\) are not observed; the frequentist likelihood treats \(\boldsymbol{\mu}\) as a fixed proxy. Estimating those latents jointly with \(\boldsymbol{p}\) requires a Bayesian / MAP step (see .map_gaussian_convolution()).

See also

gradient_loglik_unconstrained(), additive_logistic()