library( MixtureFitting )
load( "outputs/smm_002.RData" )
set.seed( 42 )
x = rt( 100, 10 )
sinit = smm_init_vector( x, 3 )
sfjt = smm_fit_em( x, sinit, polyroot.solution = 'jenkins_taub' )
sfnr = smm_fit_em( x, sinit, polyroot.solution = 'newton_raphson' )
if( all.equal( sinit, output_init ) != TRUE ) {
stop( all.equal( sinit, output_init ) )
}
if( all.equal( sfjt$p, output_parameters_jt ) != TRUE ) {
stop( all.equal( sfjt$p, output_parameters_jt ) )
}
if( all.equal( sfjt$p, output_parameters_nr ) != TRUE ) {
stop( all.equal( sfjt$p, output_parameters_nr ) )
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.