viq854 / lichee

Multi-sample cancer phylogeny reconstruction
http://viq854.github.io/lichee/
Other
34 stars 15 forks source link

-minPrivateClusterSize error #18

Open Ithinky opened 1 year ago

Ithinky commented 1 year ago

Hi, Thanks for your great tool!

When I set -minPrivateClusterSize 120 -minClusterSize 20 -cp, only minClusterSize works and all the clusters are larger than 20. Howerver, the -minPrivateClusterSize seems fails, and some private clusters are smaller than 120 such as node1 cluster.

Would you give me some help? The log is as follows: --- PHYLOGENETIC CONSTRAINT GRAPH --- numNodes = 29, numEdges = 50 NODES: level = 8: Node 0: level = 7: level = 6: Node 3: group tag = 0111111, Size: 157 VAF Mean: [ 0.75 0.85 0.99 0.99 0.99 0.93 ] Stdev: [ 0.19 0.14 0 0 0.01 0.02 ] level = 5: Node 4: group tag = 0101111, Size: 125 VAF Mean: [ 0.74 0.98 0.98 0.98 0.93 ] Stdev: [ 0.18 0.07 0.07 0.08 0.05 ] Node 7: group tag = 0011111, Size: 614 VAF Mean: [ 0.8 0.98 0.98 0.99 0.93 ] Stdev: [ 0.19 0.02 0.05 0.03 0.02 ] level = 4: Node 8: group tag = 0011110, Size: 35 VAF Mean: [ 0.75 0.99 0.98 0.99 ] Stdev: [ 0.19 0 0.02 0.01 ] Node 11: group tag = 0001111, Size: 591 VAF Mean: [ 0.98 0.98 0.98 0.93 ] Stdev: [ 0.07 0.06 0.07 0.03 ] Node 26: group tag = 0011011, Size: 80 VAF Mean: [ 0.27 0.98 0.99 0.93 ] Stdev: [ 0.21 0 0.01 0.01 ] level = 3: Node 5: group tag = 0011100, Size: 21 VAF Mean: [ 0.78 0.94 0.95 ] Stdev: [ 0.21 0.17 0.16 ] Node 6: group tag = 0001110, Size: 66 VAF Mean: [ 0.95 0.97 0.96 ] Stdev: [ 0.12 0.1 0.1 ] Node 21: group tag = 0001011, Size: 130 VAF Mean: [ 0.98 0.99 0.93 ] Stdev: [ 0 0.01 0.01 ] level = 2: Node 9: group tag = 0010100, Size: 229 VAF Mean: [ 0.62 0.99 ] Stdev: [ 0.23 0.01 ] Node 10: group tag = 0010100, Size: 21 VAF Mean: [ 0.29 0.88 ] Stdev: [ 0.14 0.16 ] Node 22: group tag = 0001010, Size: 361 VAF Mean: [ 0.97 0.93 ] Stdev: [ 0.05 0.14 ] Node 27: group tag = 0010001, Size: 24 VAF Mean: [ 0.21 0.93 ] Stdev: [ 0.14 0.02 ] level = 1: Node 1: group tag = 0000010, Size: 40 VAF Mean: [ 0.89 ] Stdev: [ 0.11 ] Node 2: group tag = 0000010, Size: 56 VAF Mean: [ 0.31 ] Stdev: [ 0.06 ] Node 12: group tag = 0000100, Size: 379 VAF Mean: [ 0.98 ] Stdev: [ 0.02 ] Node 13: group tag = 0000100, Size: 114 VAF Mean: [ 0.23 ] Stdev: [ 0.01 ] Node 14: group tag = 0000100, Size: 38 VAF Mean: [ 0.54 ] Stdev: [ 0.2 ] Node 15: group tag = 0001000, Size: 232 VAF Mean: [ 0.18 ] Stdev: [ 0 ] Node 16: group tag = 0001000, Size: 84 VAF Mean: [ 0.97 ] Stdev: [ 0.03 ] Node 17: group tag = 0100000, Size: 327 VAF Mean: [ 0.41 ] Stdev: [ 0.06 ] Node 18: group tag = 0100000, Size: 38 VAF Mean: [ 0.98 ] Stdev: [ 0.02 ] Node 19: group tag = 0100000, Size: 6494 VAF Mean: [ 0.16 ] Stdev: [ 0.04 ] Node 20: group tag = 0100000, Size: 384 VAF Mean: [ 0.68 ] Stdev: [ 0.1 ] Node 23: group tag = 0010000, Size: 6091 VAF Mean: [ 0.14 ] Stdev: [ 0.03 ] Node 24: group tag = 0010000, Size: 269 VAF Mean: [ 0.38 ] Stdev: [ 0.07 ] Node 25: group tag = 0010000, Size: 147 VAF Mean: [ 0.67 ] Stdev: [ 0.11 ] Node 28: group tag = 0000001, Size: 307 VAF Mean: [ 0.93 ] Stdev: [ 0.01 ] level = 0: EDGES: 0 -> 3 0 -> 18 1 -> 2 3 -> 4 3 -> 7 4 -> 11 5 -> 9 5 -> 10 6 -> 22 7 -> 8 7 -> 11 7 -> 26 8 -> 5 8 -> 6 9 -> 10 9 -> 12 9 -> 13 9 -> 14 9 -> 23 9 -> 24 9 -> 25 10 -> 12 10 -> 13 10 -> 14 10 -> 23 10 -> 24 11 -> 6 11 -> 21 12 -> 13 12 -> 14 14 -> 13 16 -> 15 17 -> 19 18 -> 17 18 -> 19 18 -> 20 20 -> 17 20 -> 19 21 -> 22 22 -> 1 22 -> 2 22 -> 15 22 -> 16 24 -> 23 25 -> 23 25 -> 24 26 -> 21 26 -> 27 27 -> 23 27 -> 28

Found 0 valid tree(s) Adjusting the network... 1000000 2000000 3000000 4000000 5000000 6000000 7000000 8000000 9000000 10000000 Found 0 valid trees after network adjustments