def setUp(self):
filename = os.path.join(BASE_DATA_PATH,'completeness_test_cat.csv')
parser0 = CsvCatalogueParser(filename)
self.catalogue = parser0.read_file()
self.config = {'maximum_iterations': 1000,
'number_earthquakes': 100,
'number_samples': 51,
'tolerance': 0.05}
self.model = KijkoNonParametricGaussian()
}
mmax_ksb = KijkoSellevolBayes()
mmax, mmax_sigma = mmax_ksb.get_mmax(catalogue, mmax_config)
print 'Mmax = %8.3f +/- %8.3f' % (mmax, mmax_sigma)
# In[ ]:
mmax_config = {
'number_earthquakes':
100, # Selects the N largest earthquakes in the catalogue for analysis
'input_mmax': None,
'input_mmax_uncertainty': 0.5
}
mmax_knpg = KijkoNonParametricGaussian()
mmax, mmax_sigma = mmax_knpg.get_mmax(catalogue, mmax_config)
print 'Mmax = %8.3f +/- %8.3f' % (mmax, mmax_sigma)
# In[ ]:
mmax_config = {
'number_bootstraps': 1000
} # Number of samples for the uncertainty analyis
mmax_cum_mo = CumulativeMoment()
mmax, mmax_sigma = mmax_cum_mo.get_mmax(catalogue, mmax_config)
print 'Mmax = %8.3f +/- %8.3f' % (mmax, mmax_sigma)