Estep: The E-step in the EM algorithm

Description Usage Arguments Value

View source: R/Estep.R

Description

A function to compute the expectation of the logarithm of complete likelihood w.r.t the conditional distribution of hidden status given the observed data and parameters.

Usage

1
Estep(obs, par)

Arguments

obs

A S by 5 dataframe returned by data.trans function where S is total number of SNPs.

par

A list contains the current estimates of parameters.

Value

A list contains the following elements:

e.loglik

The expectation of log-likelihood.

P.ks

a S by 9 matrix contains P(S_s=k|O,θ)

P.kls

a S by 81 matrix contains P(S_s=k,S_(s+1)=l|O,θ)


JingXieMIZZOU/hmmASE documentation built on Dec. 31, 2019, 12:52 a.m.