load( "inputs/2-2-cryptand-dihedrals-trimmed.RData" )
load( "outputs/2-2-cryptand-dihedrals.RData" )
library( MixtureFitting )
init = list()
parameters = list()
ll = list()
for( i in 1:10 ) {
vinit = vmm_init_vector( i )
vf = vmm_fit_em( dihedrals, vinit )
init[[i]] = vinit
parameters[[i]] = vf$p
ll[[i]] = llvmm( dihedrals, vf$p )
}
if( all.equal( init, output_init ) != TRUE ) {
stop( all.equal( init, output_init ) )
}
if( all.equal( parameters, output_parameters ) != TRUE ) {
stop( all.equal( parameters, output_parameters ) )
}
if( all.equal( ll, output_ll ) != TRUE ) {
stop( all.equal( ll, output_ll ) )
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.