Several lines of code were improved. The function now works with geiger 3.
The Version 8 now have options to tune the value of the critical difference between proposals of qij that will be accepted of denied.
New version of the package.
This is the current version of the function