Logic tree branch selector

This feature of the Risk Modeller's Toolkit allows users to analyse sets of hazard curves derived using the Classical PSHA-based approach, and identify the branch of the logic tree that generated the hazard curve closest to the mean hazard curve, or a specific hazard curve fractile. This IPython notebook allows the user to sort a list of hazard curves in descending order of the root mean square relative error between the branch hazard curves and mean hazard curve. An average of the RMSE distance across all locations is used for the sorting.


In [1]:
import closest_curves as cc

Please specify below the path of the directory containing the hazard curve NRML files:


In [2]:
input_dir = '../../../../rmtk_data/hazard_curves/'

In [3]:
cc.list_closest_curves(input_dir)


gsimlt_path=CauzziFaccioli2008_Toro2002SHARE_*_*_*_*_*,smlt_path=b2 0.0847919470655
gsimlt_path=ZhaoEtAl2006Asc_Toro2002SHARE_*_*_*_*_*,smlt_path=b2 0.0848318418032
gsimlt_path=AkkarBommer2010_Toro2002SHARE_*_*_*_*_*,smlt_path=b2 0.0848420322196
gsimlt_path=ChiouYoungs2008_Toro2002SHARE_*_*_*_*_*,smlt_path=b2 0.0849064677496
gsimlt_path=ChiouYoungs2008_Toro2002SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.0850259901133
gsimlt_path=ChiouYoungs2008_Toro2002SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.0850260944879
gsimlt_path=AkkarBommer2010_Toro2002SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.0850653783651
gsimlt_path=AkkarBommer2010_Toro2002SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.0850654826913
gsimlt_path=ZhaoEtAl2006Asc_Toro2002SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.0850668072792
gsimlt_path=ZhaoEtAl2006Asc_Toro2002SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.0850669116037
gsimlt_path=CauzziFaccioli2008_Toro2002SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.0854499431266
gsimlt_path=CauzziFaccioli2008_Toro2002SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.0854500469833
gsimlt_path=ChiouYoungs2008_CauzziFaccioli2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.130789343431
gsimlt_path=ChiouYoungs2008_CauzziFaccioli2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.13078963624
gsimlt_path=AkkarBommer2010_CauzziFaccioli2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.130873927961
gsimlt_path=AkkarBommer2010_CauzziFaccioli2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.13087422058
gsimlt_path=CauzziFaccioli2008_CauzziFaccioli2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.130929956138
gsimlt_path=CauzziFaccioli2008_CauzziFaccioli2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.130930248633
gsimlt_path=ZhaoEtAl2006Asc_CauzziFaccioli2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.133892636771
gsimlt_path=ZhaoEtAl2006Asc_CauzziFaccioli2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.133892922793
gsimlt_path=CauzziFaccioli2008_Campbell2003SHARE_*_*_*_*_*,smlt_path=b2 0.135733885976
gsimlt_path=CauzziFaccioli2008_Campbell2003SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.135764812219
gsimlt_path=CauzziFaccioli2008_Campbell2003SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.135764993002
gsimlt_path=ZhaoEtAl2006Asc_Campbell2003SHARE_*_*_*_*_*,smlt_path=b2 0.135773047011
gsimlt_path=AkkarBommer2010_Campbell2003SHARE_*_*_*_*_*,smlt_path=b2 0.135779376957
gsimlt_path=ChiouYoungs2008_Campbell2003SHARE_*_*_*_*_*,smlt_path=b2 0.135800339924
gsimlt_path=ChiouYoungs2008_Campbell2003SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.135921928347
gsimlt_path=ChiouYoungs2008_Campbell2003SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.135922108921
gsimlt_path=ZhaoEtAl2006Asc_Campbell2003SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.135983794726
gsimlt_path=ZhaoEtAl2006Asc_Campbell2003SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.135983975218
gsimlt_path=AkkarBommer2010_Campbell2003SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.135987024029
gsimlt_path=AkkarBommer2010_Campbell2003SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.135987204516
gsimlt_path=CauzziFaccioli2008_AkkarBommer2010_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.136894940324
gsimlt_path=CauzziFaccioli2008_AkkarBommer2010_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.136894956804
gsimlt_path=ZhaoEtAl2006Asc_AkkarBommer2010_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.137067562408
gsimlt_path=ZhaoEtAl2006Asc_AkkarBommer2010_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.137067578867
gsimlt_path=ChiouYoungs2008_AkkarBommer2010_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.137070243489
gsimlt_path=ChiouYoungs2008_AkkarBommer2010_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.137070259948
gsimlt_path=AkkarBommer2010_AkkarBommer2010_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.137070995268
gsimlt_path=AkkarBommer2010_AkkarBommer2010_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.137071011727
gsimlt_path=CauzziFaccioli2008_AkkarBommer2010_*_*_*_*_*,smlt_path=b2 0.137218098757
gsimlt_path=ZhaoEtAl2006Asc_AkkarBommer2010_*_*_*_*_*,smlt_path=b2 0.13723069386
gsimlt_path=AkkarBommer2010_AkkarBommer2010_*_*_*_*_*,smlt_path=b2 0.137241747195
gsimlt_path=ChiouYoungs2008_AkkarBommer2010_*_*_*_*_*,smlt_path=b2 0.137294748373
gsimlt_path=CauzziFaccioli2008_ChiouYoungs2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.140200756381
gsimlt_path=CauzziFaccioli2008_ChiouYoungs2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.140200770312
gsimlt_path=ZhaoEtAl2006Asc_ChiouYoungs2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.140367452207
gsimlt_path=ZhaoEtAl2006Asc_ChiouYoungs2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.140367466121
gsimlt_path=AkkarBommer2010_ChiouYoungs2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.140372046773
gsimlt_path=AkkarBommer2010_ChiouYoungs2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.140372060687
gsimlt_path=ChiouYoungs2008_ChiouYoungs2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.140384961604
gsimlt_path=ChiouYoungs2008_ChiouYoungs2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.140384975516
gsimlt_path=CauzziFaccioli2008_ChiouYoungs2008_*_*_*_*_*,smlt_path=b2 0.140579439424
gsimlt_path=ZhaoEtAl2006Asc_ChiouYoungs2008_*_*_*_*_*,smlt_path=b2 0.140586132361
gsimlt_path=AkkarBommer2010_ChiouYoungs2008_*_*_*_*_*,smlt_path=b2 0.140597554658
gsimlt_path=ChiouYoungs2008_ChiouYoungs2008_*_*_*_*_*,smlt_path=b2 0.140656249795
gsimlt_path=ChiouYoungs2008_AkkarBommer2010_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.140794652897
gsimlt_path=ChiouYoungs2008_AkkarBommer2010_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.14079490728
gsimlt_path=AkkarBommer2010_AkkarBommer2010_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.140903568143
gsimlt_path=AkkarBommer2010_AkkarBommer2010_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.140903822329
gsimlt_path=CauzziFaccioli2008_AkkarBommer2010_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.140947809394
gsimlt_path=CauzziFaccioli2008_AkkarBommer2010_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.1409480635
gsimlt_path=ZhaoEtAl2006Asc_AkkarBommer2010_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141073944924
gsimlt_path=ZhaoEtAl2006Asc_AkkarBommer2010_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141074198803
gsimlt_path=ChiouYoungs2008_ChiouYoungs2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141484129241
gsimlt_path=ChiouYoungs2008_ChiouYoungs2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141484368305
gsimlt_path=AkkarBommer2010_ChiouYoungs2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141613977877
gsimlt_path=AkkarBommer2010_ChiouYoungs2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141614216722
gsimlt_path=CauzziFaccioli2008_ChiouYoungs2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141646548848
gsimlt_path=CauzziFaccioli2008_ChiouYoungs2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141646787638
gsimlt_path=ZhaoEtAl2006Asc_ChiouYoungs2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.141738048883
gsimlt_path=ZhaoEtAl2006Asc_ChiouYoungs2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.14173828752
gsimlt_path=ChiouYoungs2008_Campbell2003SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.145216604883
gsimlt_path=ChiouYoungs2008_Campbell2003SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.145216654481
gsimlt_path=CauzziFaccioli2008_Campbell2003SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.145362515121
gsimlt_path=CauzziFaccioli2008_Campbell2003SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.145362564668
gsimlt_path=AkkarBommer2010_Campbell2003SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.145407138941
gsimlt_path=AkkarBommer2010_Campbell2003SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.145407188474
gsimlt_path=ZhaoEtAl2006Asc_Campbell2003SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.145522047807
gsimlt_path=ZhaoEtAl2006Asc_Campbell2003SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.145522097301
gsimlt_path=AkkarBommer2010_CauzziFaccioli2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.272436870971
gsimlt_path=AkkarBommer2010_CauzziFaccioli2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.272436878179
gsimlt_path=ZhaoEtAl2006Asc_CauzziFaccioli2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.27243847263
gsimlt_path=ZhaoEtAl2006Asc_CauzziFaccioli2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.272438479838
gsimlt_path=ChiouYoungs2008_CauzziFaccioli2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.272445204902
gsimlt_path=ChiouYoungs2008_CauzziFaccioli2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.27244521211
gsimlt_path=CauzziFaccioli2008_CauzziFaccioli2008_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.273085939673
gsimlt_path=CauzziFaccioli2008_CauzziFaccioli2008_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b1 0.273085946864
gsimlt_path=ZhaoEtAl2006Asc_CauzziFaccioli2008_*_*_*_*_*,smlt_path=b2 0.273209868842
gsimlt_path=CauzziFaccioli2008_CauzziFaccioli2008_*_*_*_*_*,smlt_path=b2 0.273212080723
gsimlt_path=AkkarBommer2010_CauzziFaccioli2008_*_*_*_*_*,smlt_path=b2 0.273213499085
gsimlt_path=ChiouYoungs2008_CauzziFaccioli2008_*_*_*_*_*,smlt_path=b2 0.273243965883
gsimlt_path=ChiouYoungs2008_Toro2002SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.999535378154
gsimlt_path=ChiouYoungs2008_Toro2002SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.999535398355
gsimlt_path=AkkarBommer2010_Toro2002SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.999558146551
gsimlt_path=AkkarBommer2010_Toro2002SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.999558166751
gsimlt_path=CauzziFaccioli2008_Toro2002SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.999567009306
gsimlt_path=CauzziFaccioli2008_Toro2002SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.999567029506
gsimlt_path=ZhaoEtAl2006Asc_Toro2002SHARE_Campbell2003SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 0.999999979808
gsimlt_path=ZhaoEtAl2006Asc_Toro2002SHARE_Toro2002SHARE_*_*_FaccioliEtAl2010_*,smlt_path=b3 1.0