e_step_withColSums: E-step, with column sums of posterior probabilities

In kjgleason/Primo: Package in R for Integrative Multi-Omics association analysis

Description

Calculate posterior expectations given maximized estimate for π (the proportion of observations coming from each association pattern). Return column sums of the posterior expectations. This allows the M-step to be completed in chunks (dividing a running sum by the total number of rows at the end to obtain the mean), alleviating potential memory issues.

Usage

e_step_withColSums(old_pi, Dmat)

Arguments

old_pi	vector of association pattern proportions, fit through maximization.
Dmat	matrix of conditional joint densities under each association pattern.

Value

A row vector of column sums of posterior expectations (allows M-step to be performed in chunks).

kjgleason/Primo documentation built on Sept. 7, 2021, 3:58 a.m.