cal_E_rna: cal_E_rna

In cuhklinlab/scAMACE: scAMACE (integrative Analysis of single-cell Methylation, chromatin ACcessibility, and gene Expression)

Description

Perform E-step (i.e. calculate the expectations of missing data) for one iteration of scRNA-Seq or sc-methylation data in the EM algorithm.

Usage

cal_E_rna(rlg,replaceg0,phi_rna, w_exp, pi_rna)

Arguments

rlg	ratio g1/(g1+g0) for scRNA-Seq data or ratio h1/h0 for sc-methylation data.
replaceg0	(1-rlg) for scRNA-Seq data or a matrix of constant 1 for sc-methylation data.
phi_rna	the probability for scRNA-Seq or sc-methylation data belonging to cluster k (i.e. z[l,k]=1).
w_exp	the probability for scRNA-Seq or sc-methylation data that gene g is active in cluster k.
pi_rna	the probability for scRNA-Seq or sc-methylation data that gene g is expressed in cell l.

Value

prob_z_rna	Expectation of z.
prob_z_u_rna	Expectation of z*u.
prob_z_1_u_rna	Expectation of z*(1-u).
prob_z_u_v_rna	Expectation of z*u*v for scRNA-Seq data or expectation of z*u*m for sc-methylation data.
prob_z_1_u_v_rna	Expectation of z*(1-u)*v for scRNA-Seq data or expectation of z*(1-u)*m for sc-methylation data.
prob_z_v_rna	Expectation of z*v for scRNA-Seq data or expectation of z*m for sc-methylation data.

Author(s)

Jiaxuan WANGWU

cuhklinlab/scAMACE documentation built on Dec. 19, 2021, 7:03 p.m.