The functions implement the different steps of the paper. The example.m file contains sample calculations on how to use the functions. The implementation follows the paper as much as possible. Whenever programming choices had to be made that are not described in the paper there are comments that explain the choice.
Does anyone know how to set the start values for the parameters?
Excellent implementation, one remark: I believe the restriction p >= -2 psi should also be included in the GA optimization algorithm, to ensure positive variances. Thank you for this nice implementation !
Made a typo, restriction: p >= -2 * psi.