E_step: Computes E step weights

In DavidKLim/FSCseq: Feature Selection and Clustering of RNA-seq Count Data

Description

Takes current parameter estimates and likelihoods, and computes E step weights. Also updates CEM temperature Tau, if applicable

Usage

E_step(wts, l, pi, CEM, Tau, PP_filt)

Arguments

wts	matrix of dimension k by n, current E step weights. Used as template matrix to update
l	matrix of dimension k by n, log-likelihood summed over all genes, for each k clusters and n subjects
pi	vector of k mixture proportions
CEM	logical, FALSE for EM and TRUE for CEM
Tau	numeric, current temperature of CEM. Tau=1 for EM

Value

list containing updated wts: matrix of E step weights, keep: logical matrix of whether a sample is included in cluster baseline calculation, Tau: numeric: resultant CEM temperature (after annealing), CEM: logical: whether CEM is continued to be run

DavidKLim/FSCseq documentation built on Dec. 12, 2021, 3:46 a.m.