PoonLab / Kaphi

Kernel-embedded ABC-SMC for phylodynamic inference
GNU Affero General Public License v3.0
4 stars 2 forks source link

Align is not working #129

Closed mavino closed 6 years ago

mavino commented 6 years ago

This is the command line I am using:

R
library('Kaphi')
Ptree1 <- read.tree("tree1.nwk")
Ptree2 <- read.tree("tree2.nwk")
Align(Ptree1, Ptree2)

and I get this error:

Error in if (cost[r, c] == 0 && rcov[r] == 0 && ccov[c] == 0) { : 
  missing value where TRUE/FALSE needed

These are my trees ( they have the same number of tips and same tip labels): tree1.nwk:

(((Cichlidae_P:128.0,Centrarchidae_P:128.0):101.90276096,Catostomidae_P:229.90276097):205.42131365,(Dicroglossidae_P:351.75848321999996,((Anatidae_P:98.04286929,(Psittacidae_P:83.10924465,(Gruidae_P:81.30000002,(Ciconiidae_P:81.3,Ardeidae_P:81.30000000000001):0.0):1.8092446500000001):14.93362464):213.86105206000002,(((Rhinolophidae_P:46.1,Hipposideridae_P:46.1):15.561207860000001,(Phyllostomidae_P:52.6460235,Vespertilionidae_P:52.646023500000005):9.01518435):34.80118132,(Sciuridae_P:89.82318741,(Atelidae_P:43.15129634,(Hylobatidae_P:20.18921354,Hominidae_P:20.18921354):22.962082799999997):46.671891079999995):6.63920175):215.44153219999998):39.85456185):83.5655914):361.23157320999997;

tree2.nwk:

((Catostomidae_P:229.995,(Cichlidae_P:128.783,Centrarchidae_P:128.783):101.213):206.238,(Dicroglossidae_P:354.816,((((Rhinolophidae_P:46.2441,Hipposideridae_P:46.2441):16.1834,(Phyllostomidae_P:53.1492,Vespertilionidae_P:53.1492):9.27837):34.0423,(Sciuridae_P:94.5938,(Atelidae_P:43.2172,(Hylobatidae_P:22.665,Hominidae_P:22.665):20.5522):51.3766):1.87607):215.521,(Anatidae_P:99.3081,(Ciconiidae_P:85.058,(Psittacidae_P:83.4908,(Gruidae_P:81.8608,Ardeidae_P:81.8608):1.63):1.56723):14.25):212.682):42.8261):81.4164);
gtng92 commented 6 years ago

Please let me know if you encounter any more problems with the distance metrics. Closing for now.

mavino commented 6 years ago

I confirm it's working fine. Thanks