def setUp(self):
filename = os.path.join(BASE_DATA_PATH,'completeness_test_cat.csv')
parser0 = CsvCatalogueParser(filename)
self.catalogue = parser0.read_file()
self.config = {'b-value': 1.0,
'sigma-b': 0.05,
'input_mmin': 5.0,
'input_mmax': None,
'input_mmax_uncertainty': None,
'tolerance': 0.001,
'maximum_iterations': 1000}
self.model = KijkoSellevolBayes()
mmax, mmax_sigma = mmax_ks.get_mmax(catalogue, mmax_config)
print 'Mmax = %8.3f +/- %8.3f' % (mmax, mmax_sigma)
# In[ ]:
mmax_config = {
'b-value': 1.0,
'sigma-b': 0.05,
'input_mmin': 4.5,
'input_mmax': None,
'input_mmax_uncertainty': 0.5
}
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()
