Mean and quantile loss curves computation

Bug #1168454 reported by matley
6
This bug affects 1 person
Affects Status Importance Assigned to Milestone
OpenQuake Engine
Fix Released
Critical
matley

Bug Description

Currently, the mean/quantile loss curve computation fails when there is at least one degenerate individual loss curve (e.g. a loss curve where all the losses are set to 0).

Once that risklib supports the new event based algorithm ( https://bugs.launchpad.net/oq-risklib/+bug/1168446 ) we will be ready to compute mean/quantile loss curves with the following strategy:

Given a set of loss curves (for the same asset)

0) Let L(i) the losses of the i-th curve, P(i) the poes of the i-th curve, M(i) = max(L(i))

1) let i* be the index of the loss curve having the maximum M(i) and L'(i) a subarray of L(i) where each loss <= M(i*).

2) For each i compute P'(i) by interpolating each L'(i) on the reference losses L(i*). If some L(i) is out of bounds, set the corresponding P'(i) to 0.

3) Compute the mean/quantile curves (as done previously) by considering the curves with losses L'(i) and poes P'(i)

matley (matley)
Changed in oq-engine:
status: New → In Progress
Revision history for this message
matley (matley) wrote :

The condition "if zero does not belong to L(i*) add 0 with P(0) = 1" never holds as of https://bugs.launchpad.net/oq-risklib/+bug/1168446

description: updated
matley (matley)
description: updated
Revision history for this message
Michele Simionato (michele-simionato) wrote :
Revision history for this message
matley (matley) wrote :
Revision history for this message
matley (matley) wrote :

Description updated. Now we pick up the loss curve with the maximum loss as the reference curve.

description: updated
Revision history for this message
matley (matley) wrote :
matley (matley)
Changed in oq-engine:
status: In Progress → Fix Committed
Changed in oq-engine:
status: Fix Committed → Fix Released
To post a comment you must log in.
This report contains Public information  
Everyone can see this information.

Other bug subscribers

Remote bug watches

Bug watches keep track of this bug in other bug trackers.