cal_post: cal_post

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

Description

Calculate the posterior probability for one iteration in the EM algorithm.

Usage

cal_post(phi_atac_new,w_acc_new,qi_new,phi_rna_new,w_exp_new,pi_rna_new,phi_met_new,w_met_new,pi_met_new)

Arguments

phi_atac	the probability for scCAS data belonging to cluster k (i.e. z[i,k]=1).
w_acc	the probability for scCAS data that gene g is active in cluster k.
qi	the probability that scCAS data has high gene score when gene g is active in cell i (i.e. pi_i1).
phi_rna	the probability for scRNA-Seq data belonging to cluster k (i.e. z[l,k]=1).
w_exp	the probability for scRNA-Seq data that gene g is active in cluster k.
pi_rna	the probability for scRNA-Seq data that gene g is expressed in cell l.
phi_met	the probability for sc-methylation data belonging to cluster k (i.e. z[d,k]=1).
w_met	the probability for sc-methylation data that gene g is active in cluster k.
pi_met	the probability for sc-methylation data that gene g is methylated in cell d.

Value

pst	the posterior probability for one iteration.

Author(s)

Jiaxuan WANGWU

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