rrBLUP is actually still using R

[oddoking]

Hi,
I want to speak about my understanding. This code rrBUP.py is a wrapper to get data ready and call R functions to get the answer. So in essence, it is still R?

[zhejiangzhuque]

Thanks for your comment. You are right, this code is just a R function calling wrapper temporarily, to ensure the consistency.

As you can see, for instance, in line 413 to 416 within the code with a comment format, we originally had used scipy.optimize.fminbound function in Python to replace optimize function in R, but found a different result.

We will find a solution of this issue in the future.

[oddoking]

@zhejiangzhuque

I created a pull request in an effort to use pure python code. The key change is sort eigen values/vestors before optimsing lambda.. See if my change is useful