R Under development (unstable) (2025-05-13 r88200 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> 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 ) )
+ }
> 
> proc.time()
   user  system elapsed 
  13.70    0.60   14.31