Optimum_KernelC {DeMixT}R Documentation

Kernel function for optimizing parameters and hidden variables in DeMixT

Description

This function is invoked by DeMixT_S1 and DeMixT_S2 to finish parameter estimation and expression deconvolution.

Usage

Optimum_KernelC(inputdata, groupid, nhavepi, givenpi, givenpiT, 
niter, ninteg, tol, 
sg0 = 0.5^2, mu0 = 0.0, nthread = 1)

Arguments

inputdata

A matrix of expression data (e.g gene expressions) from reference (e.g. normal) and mixed samples (e.g. mixed tumor samples). It is a GxS matrix where G is the number of genes and S is the number of samples including reference and mixed samples. Samples with the same tissue type should be placed together in columns (e.g. cbind(normal samples, mixed tumor samples)

groupid

A vector of indicators to denote if the corresponding samples are reference samples or mixed tumor samples. DeMixT is able to deconvolve mixed tumor samples with at most three components. We use 1 and 2 to denote the samples referencing the first and the second known component in mixed tumor samples. We use 3 to indicate mixed tumor samples prepared to be deconvolved. For example, in two-component deconvolution, we have c(1,1,...,3,3) and in three-component deconvolution, we have c(1,1,...,2,2,...,3,3).

nhavepi

If it is set to 0, then deconvolution is performed without any given proportions; if set to 1, deconvolution with given proportions for the first and the second known component is run; if set to 2, deconvolution is run with given tumor proportions. This option helps to do deconvolution in different settings. Because in estimation of component-specific proportions, we just use a subset of genes ; so when it is required to deconvolve another subset of genes, we just easily plug back our estimated proportions by setting this option to 1. In our two-step estimation strategy in a three-component setting, this option is set to 2 to implement the second step.

givenpi

S_{T_N}-Vector of proportions. Given the number of mixed tumor samples is S_T(S_T<S), S_{T_N} is set to 2*S_T in a three-component setting and S_T in a two-component setting. When nhavepi is 1, it is fixed with the given proportions for the first and the second known component of mixed tumor samples, or just for one known component when there is just one type of reference tissues. It has the form of Vector π^1_{N_1},pi^2_{N_1},...,π^{S_T}_{N_1}, π^1_{N_2},π^2_{N_2},...,π^{S_T}_{N_2}.

givenpiT

S_T-Vector of proportions. When nhavepi is set to 2, givenpiT is fixed with given proportions for unknown component of mixed tumor samples. This option is used when we adopt a two-step estimation strategy in deconvolution. It has the form of Vector π^1_T, π^2_T, \cdots, π^{S_T}_T. If option is not 2, this vector can be given with any element.

niter

The number of iterations used in the algorithm of iterated conditional modes. A larger value can better guarantee the convergence in estimation but increase the computation time.

ninteg

The number of bins used in numerical integration for computing complete likelihood. A larger value can increase accuracy in estimation but also increase the running time. Especially in three-component deconvolution, the increase of number of bins can greatly lengthen the running time.

tol

The convergence criterion. The default is 10^(-5).

nthread

The number of threads used for deconvolution when OpenMP is availble in the system. The default is the number of whole threads minus one. In our no-OpenMP version, it is set to 1.

sg0

Initial value for σ. The default is 0.5^2.

mu0

Initial value for μ. The default is 0.

Value

pi

Matrix of estimated proportions for each known component. The first row corresponds to the proportion estimate of each sample for the first known component (groupid = 1) and the second row corresponds to that for the second known component (groupid = 2)

decovExpr

A matrix of deconvolved expression profiles corresponding to unknown (e.g tumor) component in mixed samples for a given subset of genes. Each row corresponds to one gene and each column corresponds to one sample.

decovMu

Estimated μ of log2-normal distribution for tumor component.

decovSigma

Estimated σ of log2-normal distribution for tumor component

pi1

An S_T x I Matrix of estimated proportions for each iteration i \in \{1, \cdots, I\} for the first known component

pi2

An S_T x I Matrix of estimated proportions for each iteration i \in \{1, \cdots, I\} for the second known component

Author(s)

Zeya Wang, Wenyi Wang

See Also

http://bioinformatics.mdanderson.org/main/DeMixT

Examples

# Example 1: simulated two-component data
data.comp1 <- SummarizedExperiment::assays(test.data1.comp1)[[1]]
data.Y <- SummarizedExperiment::assays(test.data1.y)[[1]]
inputdata <- cbind(data.comp1, data.Y)
groupid <- c(rep(1, ncol(data.comp1)), rep(3, ncol(data.Y)))
Optimum_KernelC(inputdata, groupid, nhavepi = 0, 
givenpi = rep(0, 2 * ncol(data.y)), 
givenpiT = rep(0, ncol(data.y)), 
niter = 10, ninteg = 30, tol = 10^(-4))

[Package DeMixT version 1.0.2 Index]