tests/snmm_fit_em_001.R

In MixtureFitting: Fitting of Univariate Mixture Distributions to Data using Various Approaches

library( MixtureFitting )

x = read.csv("inputs/enzyme.csv", header = TRUE)[,1]

# Resulting values taken from Lin et al. (2007), Statistica Sinica 17, 909-92, Table 1
p = c( 0.6240, 1 - 0.6240, 0.0949, 0.7802, 0.1331, 0.7150, 3.2780, 6.6684 )

sf = snmm_fit_em( x, snmm_init_vector( x, 2 ) )

if( !all( abs( sf - p ) < 0.1 ) ) {
    stop( sf )
}