Description Usage Arguments Details Value
Given the individual component likelihoods for a mixture model, estimates the mixture proportions.
1 |
matrix_lik, |
a n by k matrix with (j,k)th element equal to f_k(x_j). |
prior, |
a k vector of the parameters of the Dirichlet prior on π. Recommended to be rep(1,k) |
pi_init, |
the initial value of π to use. If not specified defaults to (1/k,...,1/k). |
control |
A list of control parameters for the SQUAREM algorithm, default value is set to be control.default=list(K = 1, method=3, square=TRUE, step.min0=1, step.max0=1, mstep=4, kr=1, objfn.inc=1,tol=1.e-07, maxiter=5000, trace=FALSE). |
weights |
weights to be assigned to the observations (an n vector) |
Optimizes
L(pi)= sum_j w_j log(sum_k pi_k f_{jk}) + h(pi)
subject to pi_k non-negative and sum_k pi_k = 1. Here
h(pi)
is a penalty function h(pi) = sum_k (prior_k-1) log pi_k. Calls REBayes::KWDual in the REBayes package, which is in turn a wrapper to the mosek convex optimization software. So REBayes must be installed to use this. Used by the ash main function; there is no need for a user to call this function separately, but it is exported for convenience.
A list, including the estimates (pihat), the log likelihood for each interation (B) and a flag to indicate convergence
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.