em_iter: EM Iteration

Description Usage Arguments Value

View source: R/RcppExports.R

Description

Complete one iteration of the EM-algorithm (combines E-step and M-step). Calculate posterior expectations given current estimate for π,the proportion of observations (e.g. SNPs) coming from each association pattern (i.e. E-Step). Re-estimate the π vector that maximizes the posterior expectations (i.e. M-step).

Usage

1
em_iter(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 estimating the proportion of observations coming from each association pattern.


kjgleason/primo documentation built on Sept. 7, 2021, 5:21 p.m.