docs: def loss_ratio_exceedance_matrix(self, loss_ratios):
          """
          Compute the LREM (Loss Ratio Exceedance Matrix).
          """
          # LREM has number of rows equal to the number of loss ratios
          # and number of columns equal to the number of imls
          lrem = numpy.empty((len(loss_ratios), len(self.imls)))
          for row, loss_ratio in enumerate(loss_ratios):
              for col, (mean_loss_ratio, stddev) in enumerate(
                      zip(self.mean_loss_ratios, self.stddevs)):
                  lrem[row, col] = self.distribution.survival(
                      loss_ratio, mean_loss_ratio, stddev)
          return lrem