Learn R Programming

ashr (version 2.0.5)

estimate_mixprop: Estimate mixture proportions of a mixture g given noisy (error-prone) data from that mixture.

Description

Estimate mixture proportions of a mixture g given noisy (error-prone) data from that mixture.

Usage

estimate_mixprop(data, g, prior, optmethod = c("mixEM", "mixVBEM", "cxxMixSquarem", "mixIP"), control)

Arguments

data
list to be passed to log_comp_dens_conv; details depend on model
g
an object representing a mixture distribution (eg normalmix for mixture of normals; unimix for mixture of uniforms). The component parameters of g (eg the means and variances) specify the components whose mixture proportions are to be estimated. The mixture proportions of g are the parameters to be estimated; the values passed in may be used to initialize the optimization (depending on the optmethod used)
prior
numeric vector indicating parameters of "Dirichlet prior" on mixture proportions
optmethod
name of function to use to do optimization
control
list of control parameters to be passed to optmethod, typically affecting things like convergence tolerance

Value

list, including the final loglikelihood, the null loglikelihood, an n by k likelihood matrix with (j,k)th element equal to $f_k(x_j)$, the fit and results of optmethod

Details

This is used by the ash function. Most users won't need to call this directly, but is exported for use by some other related packages.