Some R-code, I am developing for around half a year:
https://bitcointalk.org/index.php?topic=4586924.0and
https://github.com/trantuteThe bitcointalk-link is in german and I am too lazy to translate it in the moment. However, if there is real interest for a translation, I would invest the time. The code should speak for itself.
Tell me what you think or even improve it.
Have fun.