tupdate_modek: Update the mode k variational density.

Description Usage Arguments Value Author(s)

Description

Update the mode k variational density.

Usage

1
2
tupdate_modek(Y, ex_list, ex2_vec, esig, k, mixcompdist = "normal",
  which_na = NULL, nullweight = 10)

Arguments

Y

An array of numerics. The data array.

ex_list

A list of vectors of numerics. The expected values of the components of each mode.

ex2_vec

A vector of positive numerics. The expected value of x'x for each mode.

esig

A positive numeric. The expected value of the precision.

k

A positive integer. The current mode to update.

mixcompdist

The mixing distribution to assume. Defaults to normal. Options are those available in the ashr package.

which_na

Either NULL (when complete data) or an array of logicals the same dimension as Y, indicating if the observation is missing (TRUE) or observed (FALSE).

nullweight

A numeric greater than or equal to 1. The penalty term on the probability of a component being zero.

Value

ex_list A list of vectors of the starting expected values of each component.

ex2_vec A vector of starting values of E[x'x]

Author(s)

David Gerard


kkdey/flashr documentation built on May 20, 2019, 10:36 a.m.