Computational overflow problem

[zjm960521]

When I submit the final script of the script, there is a computational overflow problem.

I don't have NA in my data.
How do I modify and successfully run the script to produce the final result?
Thank you.

[yangymargaret]

Hi @zjm960521 ,
I really apologize for the delay in my reply to you. I am very sorry about this delay. I think there is an overflow problem in updating the parameters in the optimization step based on the screenshot you showed. Would you please first try using the same version of NumPy and the other packages as shown in the readme? If the error is the same, would you please trying printing the beta1, lambda1, and theta1 values after Line 518 in the phylo_hmgp.py? Please consider using the command python phylo_hmpy.py [options] > log.txt. Then the intermediate results will be kept in the log.txt file. Would you please then share the log file with me? I will look into what the problems may be. Thank you very much!

Best regards!
Yang Yang