def get_site_residuals(self, database, component="Geometric"):
"""
Calculates the total, inter-event and within-event residuals for
each site
"""
imt_dict = dict([(imtx, {}) for imtx in self.imts])
for site_id in self.site_ids:
print site_id
selector = SMRecordSelector(database)
site_db = selector.select_from_site_id(site_id, as_db=True)
resid = Residuals(self.input_gmpe_list, self.imts)
resid.get_residuals(site_db, normalise=False, component=component)
setattr(resid, "site_analysis", self._set_empty_dict())
setattr(resid, "site_expected", self._set_empty_dict())
self.site_residuals.append(resid)
def get_site_residuals(self, database):
"""
Calculates the total, inter-event and within-event residuals for
each site
"""
imt_dict = dict([(imtx, {}) for imtx in self.imts])
for site_id in self.site_ids:
print site_id
selector = SMRecordSelector(database)
site_db = selector.select_from_site_id(site_id, as_db=True)
resid = Residuals(self.input_gmpe_list, self.imts)
resid.get_residuals(site_db, normalise=False)
setattr(
resid,
"site_analysis",
self._set_empty_dict())
setattr(
resid,
"site_expected",
self._set_empty_dict())
self.site_residuals.append(resid)
def get_site_residuals(self, database):
"""
Calculates the total, inter-event and within-event residuals for
each site
"""
imt_dict = dict([(imtx, {}) for imtx in self.imts])
for site_id in self.site_ids:
print site_id
selector = SMRecordSelector(database)
site_db = selector.select_from_site_id(site_id, as_db=True)
resid = Residuals(self.gmpe_list, self.imts)
resid.get_residuals(site_db, normalise=False)
setattr(
resid,
"site_analysis",
OrderedDict([(gmpe, imt_dict) for gmpe in self.gmpe_list]))
setattr(
resid,
"site_expected",
OrderedDict([(gmpe, imt_dict) for gmpe in self.gmpe_list]))
self.site_residuals.append(resid)
