search

morrislab / pairtree

Pairtree is a method for reconstructing cancer evolutionary history in individual patients, and analyzing intratumor genetic heterogeneity. Pairtree focuses on scaling to many more cancer samples and cancer cell subpopulations than other algorithms, and on producing concise and informative interactive characterizations of posterior uncertainty.
MIT License
37 stars 11 forks source link

ValueError when plotting the tree #36

Closed OumaimaHamza closed 1 year ago

OumaimaHamza commented 1 year ago

Hi, I got this error when running the /bin/plotree command,

Screen Shot 2023-02-24 at 8 30 05 AM I dont understand why that's happening, since the computing and the sampling of trees worked and i got my .results.npz file and all my var_read_prob are between [0,1]

Would really appreciate your feedback on this

Thank you !

ethanumn commented 1 year ago

Hi there --

This is quite unusual. Is there any chance that the var_read_prob for one of the mutations you used is 0? I don't believe there is a check in that function that avoids dividing by a var_read_prob that is equal to 0, and this could very well be the issue.

Ethan

OumaimaHamza commented 1 year ago

Thank you so much Ethan, it works now !