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)