christiaanjs / treeflow

GNU General Public License v3.0
13 stars 4 forks source link

Numerical stability of very short branches with scaled rate approximation #27

Open christiaanjs opened 3 years ago

christiaanjs commented 3 years ago

Fitting a variational approximation on sequences of length 10000 simulated from this tree create NaNs when using the scaled rate approximation (mean field succeeds): ((((((((((T10_9.473684210526315:0.01339,(T11_10.526315789473683:0.00028,(T12_11.578947368421051:0.02809,(T13_12.631578947368421:0.01531,(T14_13.68421052631579:0.00846,(T15_14.736842105263158:0.01458,(T16_15.789473684210526:0.00462,(T17_16.842105263157894:0.01904,(T18_17.894736842105264:0.02114,(T19_18.94736842105263:0.00396,T20_20.0:1.05659):1.06981):1.05054):1.03821):1.06260):1.04651):1.05948):1.06541):1.02482):1.06574):1.04154,T9_8.421052631578947:0.00230):1.05765,T8_7.368421052631579:0.00731):1.04649,T7_6.315789473684211:0.00117):1.05244,T6_5.263157894736842:0.00097):1.06488,T5_4.2105263157894735:0.01322):1.08433,T4_3.1578947368421053:0.04492):1.09101,T3_2.1052631578947367:0.08330):0.99583,T2_1.0526315789473684:0.02650):1.06187,T1_0.0:0.03574):0.00000;

Full simulation replicate: error-replicate.tar.gz