JP6192432B2 - Risk weighing system - Google Patents

Risk weighing system Download PDF

Info

Publication number
JP6192432B2
JP6192432B2 JP2013172162A JP2013172162A JP6192432B2 JP 6192432 B2 JP6192432 B2 JP 6192432B2 JP 2013172162 A JP2013172162 A JP 2013172162A JP 2013172162 A JP2013172162 A JP 2013172162A JP 6192432 B2 JP6192432 B2 JP 6192432B2
Authority
JP
Japan
Prior art keywords
distribution
estimated
loss
accident
estimation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2013172162A
Other languages
Japanese (ja)
Other versions
JP2015041259A (en
Inventor
前川 拓也
拓也 前川
Original Assignee
株式会社日立ソリューションズ西日本
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社日立ソリューションズ西日本 filed Critical 株式会社日立ソリューションズ西日本
Priority to JP2013172162A priority Critical patent/JP6192432B2/en
Publication of JP2015041259A publication Critical patent/JP2015041259A/en
Application granted granted Critical
Publication of JP6192432B2 publication Critical patent/JP6192432B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Financial Or Insurance-Related Operations Such As Payment And Settlement (AREA)

Description

本発明は、リスク計量化技術に関し、特に事故発生頻度と損失額に対する確率分布とそのパラメータを推定し、推定した事故件数分布と損失額分布の確率分布関数を利用してリスクを計量するリスク計量システムにおいて、推定した分布と事故データとの当てはまり具合の評価のために表示を行うシステムに関する。   The present invention relates to a risk quantification technique, and in particular, a risk metric that estimates a probability distribution and its parameters with respect to an accident occurrence frequency and a loss amount, and measures the risk using a probability distribution function of the estimated number of accidents and loss distribution. The present invention relates to a system that performs display for evaluating the degree of fit between an estimated distribution and accident data.

一般に企業の業務は、システム障害、事務的なミス、詐欺など様々なリスクに遭遇する可能性がある。このため、事故データを収集してリスクを計量し、リスクに対する対策を講じることが求められている。一般的なリスク計量化は、例えば特許文献1及び非特許文献1のようにして、VaRと呼ばれるリスク量を算出する。   In general, business operations can encounter various risks such as system failures, clerical errors, and fraud. For this reason, it is required to collect accident data, measure the risk, and take measures against the risk. In general risk quantification, for example, as in Patent Document 1 and Non-Patent Document 1, a risk amount called VaR is calculated.

具体的には、まず、収集された事故データの一定期間当たりの事故件数などから事故件数分布を推定し、発生した事故の損失額などから事故1件当たりの損失額分布を推定する。そして、モンテカルロ・シミュレーションにより、上記の事故件数分布を用いて発生させた損失件数分の損失額を、上記の損失額分布から発生させて合計し、一定期間当たりの損失総額を算出する処理を何万、何十万回と繰り返して損失総額の分布を生成する。そして、この生成された損失総額の分布から、所定の分位点の損失総額をVaR(バリュー・アット・リスク)とする。例えば、1年間に発生し得る損失総額の分布をシミュレーションによって生成し、そこから百年に一度あるいは千年に一度発生し得る巨額損失を計量する。   Specifically, first, the accident number distribution is estimated from the number of accidents per certain period of the collected accident data, and the loss distribution per accident is estimated from the loss amount of the accident that occurred. Then, by Monte Carlo simulation, what is the processing to calculate the total loss per fixed period by summing up the loss amount generated from the above loss amount distribution using the above accident number distribution? Generate a distribution of the total loss by repeating the number of times. Then, from the generated distribution of the total loss, the total loss at a predetermined quantile is defined as VaR (Value at Risk). For example, a distribution of the total loss that can occur in one year is generated by simulation, and a huge loss that can occur once in a hundred years or once in a thousand years is measured therefrom.

ここで、推定した確率分布と事故データとの当てはまり具合の評価方法としては、非特許文献1のようにコルモゴロフ・スミルノフ検定やアンダーソン・ダーリング検定といった定量的な検定手法のほか、p-pプロットやq-qプロットなどを用いて、それらの当てはまり具合をグラフ表示して判断させる手法がある。   Here, as a method of evaluating the degree of fit between the estimated probability distribution and the accident data, in addition to quantitative test methods such as Kolmogorov-Smirnov test and Anderson-Darling test as in Non-Patent Document 1, pp plot and qq plot For example, there is a technique for determining the degree of fit using a graph.

またリスク計量法として、複数のリスクシナリオごとに分布を推定して、それらの損失総額をシミュレーション内で合計することで、高精度のリスク計量が行えるようにする手法も用いられている。   In addition, as a risk measurement method, a method is used in which a highly accurate risk measurement can be performed by estimating a distribution for each of a plurality of risk scenarios and totaling the total loss in a simulation.

特許第4241083号公報Japanese Patent No. 4241083

金融機関のオペレーショナル・リスク管理態勢(金融検査マニュアルハンドブックシリーズ)p64-65、p114-119Operational risk management system of financial institutions (Financial Inspection Manual Handbook Series) p64-65, p114-119

一般的に、損失額分布の推定には対数正規分布や一般化パレート分布などが用いられることが多いが、事故データとの当てはまり具合がよくないことが知られており、損失額分布を表す有効な理論分布関数は未だ発見されていないのが現状である。   In general, logarithmic normal distribution and generalized Pareto distribution are often used to estimate the loss distribution, but it is known that it does not fit well with accident data. The current theoretical distribution function has not yet been discovered.

そのため、検定による定量的な評価ではあまりよい結果は得られないことが多い。また検定は、データが特定の確率分布に従っているか否かを評価するものであり、複数の推定分布候補に対しての当てはまり具合の比較に用いることを想定したものではない。したがって、推定分布の当てはまり具合はq-qプロットなどのグラフを見ながら人間が経験的に評価し、妥当と思われる推定分布候補に対するVaRについて検討する必要がある。このとき、推定分布の当てはまり具合を評価するための表示方法として、累積分布関数のプロットのほかにp-pプロットやq-qプロットが用いられている。   For this reason, quantitative evaluation by testing often does not give very good results. The test evaluates whether or not the data follows a specific probability distribution, and is not intended to be used for comparison of the degree of fit to a plurality of estimated distribution candidates. Therefore, it is necessary for humans to empirically evaluate the fit of the estimated distribution while looking at graphs such as the q-q plot, and to examine VaR for the estimated distribution candidate that seems to be appropriate. At this time, as a display method for evaluating the fit of the estimated distribution, a p-p plot and a q-q plot are used in addition to the cumulative distribution function plot.

累積分布関数の表示方法は、推定分布の累積分布関数と、事故データの損失額とそれ以下の事故件数の全事故件数との比率を用いたグラフを重ねて表示するものである。p-pプロット及びq-qプロットは、上記累積分布を使ったグラフに対して、p-pプロットでは横軸の、q-qプロットでは縦軸の変数を媒介変数として、事故データのグラフの値と、推定分布のグラフの値の組を座標とした2次元のグラフに表示するものである。   The display method of the cumulative distribution function is a method of displaying a cumulative distribution function of an estimated distribution and a graph using a ratio between the loss amount of accident data and the total number of accidents less than that. The pp plot and qq plot are graphs of the accident data graph and the estimated distribution graph with the horizontal axis in the pp plot and the vertical axis variable in the qq plot as parameters. It is displayed on a two-dimensional graph with a set of values as coordinates.

事故データは一般的に、損失額が高額になるほど急激に発生頻度が下がる傾向があるため、グラフ表示によって高頻度低額損失から低頻度高額損失までの全体的な当てはまり具合を評価するのが困難である。
また、損失額の対数値で表示することが広く行われているが、この場合では損失額のずれの比がグラフ上でのずれとして表示されるため、千円と2千円の損失のずれと、百万円と2百万円の損失のずれが、同一のずれとして表されてしまい、当てはまり具合のずれを損失額ベースで評価しにくくなるという問題がある。
また、推定分布の累積分布は損失額が無限大のときに確率1となることが多いため、当てはまり具合の評価に事故データの最高損失額を考慮することが困難であるという問題がある。
Accident data generally tends to decrease in frequency more rapidly as the amount of loss increases, so it is difficult to evaluate the overall fit from high frequency low loss to low frequency high loss by graph display. is there.
In addition, the logarithm value of the loss amount is widely used. In this case, the ratio of the difference in loss amount is displayed as a difference on the graph. The difference between the loss of 1 million yen and the loss of 2 million yen is expressed as the same difference, and there is a problem that it is difficult to evaluate the difference in fit on a loss basis.
Further, since the cumulative distribution of the estimated distribution often has a probability of 1 when the amount of loss is infinite, there is a problem that it is difficult to consider the maximum amount of loss of accident data in the evaluation of the fit.

上記の問題点に鑑み本発明は、リスク計量システムにおいて事故データと推定分布のグラフ表示における当てはまり具合を評価しやすくすることを目的とする。   In view of the above-described problems, an object of the present invention is to make it easy to evaluate the degree of fit in the graph display of accident data and estimated distribution in a risk weighing system.

上記目的を達成するため、本発明は、以下の構成を提供する。
本発明の態様は、事故発生日時と損失金額を含む収集された事故データから所定の推定条件により選択された分布推定用事故データに対して、指定された確率分布関数及び指定されたパラメータ推定法を適用して事故件数分布と損失額分布を推定する分布推定手段と、
推定された前記事故件数分布と前記損失額分布に従う乱数を発生させて一定期間内に発生する損失総額を算出するシミュレーションを複数回実行するシミュレーション手段と、
前記シミュレーション手段によって得られた損失総額データにおける指定された分位点の損失総額を算出するリスク計量手段と、を備えた処理装置を有するリスク計量システムにおいて、前記処理装置は、
前記分布推定手段において事故件数分布と損失額分布を表す複数の確率分布関数及び複数のパラメータ推定法により推定された各推定分布の候補について、
当該推定分布と前記分布推定用事故データとの当てはまり具合の評価のために、当該推定分布の累積分布逆関数の積分関数と、前記分布推定用事故データに基づいた損失額と一定損失額以下の累積事故件数の全事故件数に対する比率と一定損失額以下の累積損失額とにより表されるデータとをグラフに表示する推定分布評価手段を有することを特徴とする。
In order to achieve the above object, the present invention provides the following configurations.
An aspect of the present invention provides a specified probability distribution function and a specified parameter estimation method for accident data for distribution estimation selected according to a predetermined estimation condition from collected accident data including an accident occurrence date and time and a loss amount. A distribution estimation means that estimates the accident number distribution and loss distribution by applying
Simulation means for generating a random number according to the estimated accident number distribution and the loss amount distribution to calculate a total loss occurring within a predetermined period, a simulation means for executing a plurality of times;
In a risk weighing system comprising: a risk weighing unit that calculates a total loss of a specified quantile in the total loss data obtained by the simulation unit, the processing device includes:
For each estimated distribution candidate estimated by a plurality of probability distribution functions and a plurality of parameter estimation methods representing the accident number distribution and loss distribution in the distribution estimation means,
In order to evaluate the degree of fit between the estimated distribution and the distribution estimation accident data, the integral function of the cumulative distribution inverse function of the estimated distribution, the loss amount based on the distribution estimation accident data, and a certain loss amount or less It has an estimated distribution evaluation means for displaying on a graph the data represented by the ratio of the cumulative number of accidents to the total number of accidents and the cumulative loss amount below a certain loss amount.

上記態様における前記推定分布評価手段において、
前記分布推定用事故データについては、横軸を一定損失額以下の累積事故件数の全事故件数に対する比率とし、縦軸を累積損失額若しくはその対数値、又は累積損失額の損失総額に対する比率としたグラフを表示し、
前記推定分布の累積分布逆関数の積分関数については、横軸を確率とし、縦軸を前記積分関数の値若しくはその対数値、又は確率1のときの前記積分関数の値に対する比率としたグラフを表示し、
これらのグラフを並べて表示又は重ね合わせて表示することが、好適である。
In the estimated distribution evaluation means in the above aspect,
For the accident data for distribution estimation, the horizontal axis is the ratio of the total number of accidents with a certain loss amount or less to the total number of accidents, and the vertical axis is the cumulative loss amount or its logarithm, or the ratio of the cumulative loss amount to the total loss Display the graph,
For the integral function of the cumulative distribution inverse function of the estimated distribution, a graph in which the horizontal axis is the probability and the vertical axis is the value of the integral function or its logarithmic value, or the ratio to the value of the integral function when the probability is 1. Display
It is preferable to display these graphs side by side or superimposed.

上記態様における前記推定分布評価手段において、
事故データに対するグラフの横軸と、推定分布の累積分布逆関数の積分関数に対するグラフの横軸を1件当たりの事故の損失額として表示することが、好適である。
In the estimated distribution evaluation means in the above aspect,
It is preferable to display the horizontal axis of the graph with respect to the accident data and the horizontal axis of the graph with respect to the integral function of the cumulative distribution inverse function of the estimated distribution as an accident loss amount per case.

上記態様における前記推定分布評価手段において、
事故データと推定分布のそれぞれのグラフを重ね合わせたときの、横軸の変数を媒介変数として、事故データと推定分布のそれぞれのグラフ上での値の組を縦軸と横軸にとったグラフを表示することが、好適である。
In the estimated distribution evaluation means in the above aspect,
A graph in which the horizontal axis variable is the parameter when the accident data and estimated distribution graphs are superimposed, and the vertical and horizontal axes are pairs of values on the accident data and estimated distribution graphs. Is preferably displayed.

本発明では、推定分布の累積分布逆関数に対する積分関数と、事故データでの損失額の最低額からの累積損失を利用し、これらを比較できるようなグラフを表示する。具体的には、横軸に確率をとり、事故データと推定分布に対してその確率で起こる損失額の最低損失からの累積損失額を算出し、その値を縦軸にとってグラフを表示する。また、横軸を損失額とし、その金額以下の事故に対する損失の累積額を縦軸にとって表示する。これら指定分布と事故データのグラフは重ね合わせて表示することができる。また、これらのグラフの横軸の変数を媒介変数として、推定分布と事故データのそれぞれの累積損失額の組を座標とし、縦軸と横軸にとったグラフを表示することで、当てはまり具合を損失額ベースで評価できるようにする。   In the present invention, the integral function for the cumulative distribution inverse function of the estimated distribution and the cumulative loss from the minimum loss amount in the accident data are used, and a graph that can compare these is displayed. Specifically, taking the probability on the horizontal axis, the cumulative loss amount from the lowest loss that occurs at that probability is calculated for the accident data and the estimated distribution, and the graph is displayed with the value as the vertical axis. In addition, the horizontal axis represents the amount of loss, and the vertical axis represents the cumulative amount of loss for accidents below that amount. These specified distributions and accident data graphs can be displayed superimposed. In addition, the horizontal axis of these graphs is used as a parameter, and the cumulative loss amounts of each of the estimated distribution and accident data are used as coordinates, and the vertical and horizontal axes are displayed to show the fit. Make it possible to evaluate on a loss basis.

これらのプロット方法によるグラフ表示により、既存のp-pプロットやq-qプロット及び検定などによる定量的な評価と、シミュレーションによるVaRの値を、推定分布の各候補に対して任意に組み合わせて表示し、推定分布候補から適切と思われるVaRを選択するための支援を行うことが可能となる。   By displaying the graphs using these plotting methods, quantitative evaluations using existing pp plots, qq plots, and tests, and VaR values from simulations are displayed in any combination for each estimated distribution candidate. It is possible to provide support for selecting VaR that seems appropriate from the candidates.

本発明により、リスク計量システムで事故データとそれに対する推定分布との当てはまり具合を、高頻度低額損失から低頻度高額損失まで全体的にバランスよく表示することができ、推定分布の平均値が有限な場合は事故データでの最大損失額の点がプロット可能となり、推定分布の当てはまり具合を損失額ベースで評価することができる。   According to the present invention, it is possible to display the fit between the accident data and the estimated distribution corresponding thereto in the risk weighing system in a well-balanced manner from the high frequency low loss to the low frequency high loss, and the average value of the estimated distribution is finite. In this case, the maximum loss point in the accident data can be plotted, and the fit of the estimated distribution can be evaluated on a loss basis.

図1は、本実施形態の一例のシステム構成図である。FIG. 1 is a system configuration diagram of an example of the present embodiment. 図2は、本実施形態のデータ構造の一例を表す図である。FIG. 2 is a diagram illustrating an example of the data structure of the present embodiment. 図3は、本実施形態の全体処理の一例を表すフローチャート図である。FIG. 3 is a flowchart showing an example of the overall processing of this embodiment. 図4は、図3の推定用事故データテーブル生成処理の一例を表すフローチャート図である。FIG. 4 is a flowchart showing an example of the estimation accident data table generation process of FIG. 図5は、図3の分布推定処理の一例を表すフローチャート図である。FIG. 5 is a flowchart showing an example of the distribution estimation process of FIG. 図6は、図3のリスク計量処理の一例を表すフローチャート図である。FIG. 6 is a flowchart showing an example of the risk measurement process of FIG. 図7は、図3の出力処理におけるグラフ表示例を表す図である。FIG. 7 is a diagram illustrating a graph display example in the output processing of FIG.

以下、本発明のリスク計量システムの実施形態を、一例として示した図面に従って説明する。リスク計量システムの本実施形態において、事故データの収集については、本システムにリスク管理手段を付加するか、あるいは別途リスク管理システムが行うものとし、本システムはそこから事故データを参照し、それに基づいて分布推定及びリスク計量を行うものとする。よって、事故データの収集の局面については、本システムに含まれないものとする。   Hereinafter, embodiments of the risk measurement system of the present invention will be described with reference to the drawings shown as an example. In this embodiment of the risk measurement system, accident data is collected by adding risk management means to this system or by a separate risk management system. This system refers to accident data from there and based on it. Distribution estimation and risk measurement. Therefore, accident data collection aspects are not included in this system.

また、一般的なリスク計量では、事故の属するリスクシナリオごとに分布を推定することが行われるが、リスク計量では各シナリオごとの損失額の乱数を合計して損失総額とすればよいので、本実施形態の一例ではリスクシナリオごとの推定及びリスク計量の合計などの処理は省略し、シナリオが1つであるものとして説明する。   In general risk measurement, the distribution is estimated for each risk scenario to which the accident belongs, but in risk measurement, it is only necessary to add the random number of losses for each scenario to make the total loss. In the example of the embodiment, processing such as estimation for each risk scenario and total risk metric is omitted, and description will be made assuming that there is one scenario.

図1は、本システムの実施形態の一例を示したシステム構成図である。本システムは、入力装置101、プリンタなどの出力装置102、ディスプレイなどの表示装置103、処理装置104、及びデータベース(DB)110を備える。   FIG. 1 is a system configuration diagram illustrating an example of an embodiment of the present system. The system includes an input device 101, an output device 102 such as a printer, a display device 103 such as a display, a processing device 104, and a database (DB) 110.

処理装置104は、推定用事故データ生成部106、分布推定部107、リスク計量部108、及び出力処理部109を含む一連の機能処理部105を備えている。これらの機能処理部105は、コンピュータである処理装置104に導入された所定のプログラムをメモリに読み込みCPUが当該プログラムを実行することにより実現される。   The processing device 104 includes a series of function processing units 105 including an estimation accident data generation unit 106, a distribution estimation unit 107, a risk measurement unit 108, and an output processing unit 109. These function processing units 105 are realized by reading a predetermined program installed in the processing device 104, which is a computer, into a memory and executing the program by the CPU.

DB110には、収集された事故データを格納した事故データテーブル111、分布推定に用いる分布推定用事故データテーブル112、推定された分布とそのパラメータなどを格納する推定分布テーブル113、分布推定用事故データに対する各種統計量を格納する事故データ統計量テーブル114、推定された分布に対する各種統計量を格納する推定分布統計量テーブル115、シミュレーションで算出された損失総額を格納するシミュレーション結果テーブル116を格納する。なお、事故データテーブル111は、事故データ収集処理を行う別のシステムが保有していてもよく、その場合は本システムからデータを参照可能であればよい。   The DB 110 includes an accident data table 111 storing collected accident data, a distribution estimation accident data table 112 used for distribution estimation, an estimated distribution table 113 storing estimated distributions and parameters thereof, and distribution estimation accident data. Accident data statistic table 114 that stores various statistics for the estimated distribution, an estimated distribution statistic table 115 that stores various statistics for the estimated distribution, and a simulation result table 116 that stores the total loss calculated in the simulation. The accident data table 111 may be held by another system that performs the accident data collection process. In this case, it is only necessary that the data can be referred to from this system.

図2(a)〜(f)は、本実施形態のデータ構造の一例を表す図である。
図2(a)の事故データテーブル111は、発生した事故に関する情報を格納し、事故データを特定するためのIDとしての「項番」201、事故の「発生日時」202、事故によって発生した「損失額」203、「事故コード」204などを含む。事故コードには、例えば、発生した事故の業務種別やその原因などを識別する情報や、担当従業員のIDなど、リスクシナリオごとの推定を行う際に必要となる情報が含まれている。
2A to 2F are diagrams illustrating an example of the data structure of the present embodiment.
The accident data table 111 in FIG. 2A stores information on the accident that occurred, “item number” 201 as the ID for identifying the accident data, “occurrence date and time” 202 of the accident, “ Loss amount "203," accident code "204, and the like. The accident code includes, for example, information necessary for estimation for each risk scenario, such as information for identifying the business type of the accident that has occurred and its cause, and the ID of the employee in charge.

図2(b)の分布推定用事故データテーブル112には、事故件数分布や損失額分布を推定するために用いられる事故データを、事故データテーブル111から選択して取得し、分布推定用事故データとして格納する。推定単位を識別するための「データセットID」205、推定に用いる事故データの期間やリスクシナリオなどの条件を格納する「推定条件」206、推定条件206を満たす項番201を格納するための「事故データ項番」207、「事故件数比率」208、「累積損失額」209、「累積損失額比率」210を格納する。   In the accident estimation data table 112 for distribution estimation in FIG. 2 (b), accident data used for estimating the number of accidents distribution and the loss amount distribution is selected and acquired from the accident data table 111, and the accident data for distribution estimation is acquired. Store as. “Data set ID” 205 for identifying the estimation unit, “estimation condition” 206 for storing conditions such as the period and risk scenario of the accident data used for estimation, and “item number 201 for storing the estimation condition 206” Accident data item number "207," accident number ratio "208," cumulative loss amount "209, and" accumulated loss amount ratio "210 are stored.

分布推定用事故データは、同一データセットIDに対して、損失額203の昇順にソートされている。事故件数比率208、累積損失額209及び累積損失額比率210は、最低損失額のデータからの事故件数又は損失額を順番に累積することにより、それぞれ算出される。事故件数比率208は、同一データセットIDの全事故件数に対する累積事故件数の比率である。累積損失額比率210は、同一データセットIDの損失総額に対する累積損失額209の比率である。また、これらは必要に応じて最高損失額からの累積として算出することもあるが、本実施例では省略する。   The accident data for distribution estimation is sorted in ascending order of the loss amount 203 with respect to the same data set ID. The accident number ratio 208, the cumulative loss amount 209, and the cumulative loss ratio 210 are calculated by sequentially accumulating the number of accidents or the loss amount from the minimum loss data. The accident number ratio 208 is a ratio of the cumulative number of accidents to the total number of accidents of the same data set ID. The cumulative loss amount ratio 210 is a ratio of the cumulative loss amount 209 to the total loss amount of the same data set ID. Further, these may be calculated as accumulation from the maximum loss amount as necessary, but are omitted in this embodiment.

図2(c)の推定分布テーブル113は、分布推定用事故データテーブル112に当てはまる理論分布を推定した結果を格納する。分布を特定するためのIDとしての「項番」211、推定に用いた「データセットID」212、「分布種別ID」213、「推定法ID」214、「推定パラメータ」215からなる。分布種別ID213は、例えばポアソン分布や対数正規分布などといった理論分布を識別するものである。推定法ID214は、分布種別ID213に対して、例えば最小二乗法や最尤法などといったパラメータ推定法を識別するものである。これらによって推定されたパラメータが推定パラメータ215に格納される。   The estimated distribution table 113 in FIG. 2C stores a result of estimating a theoretical distribution that applies to the distribution estimation accident data table 112. It consists of “item number” 211 as an ID for specifying the distribution, “data set ID” 212 used for estimation, “distribution type ID” 213, “estimation method ID” 214, and “estimation parameter” 215. The distribution type ID 213 identifies a theoretical distribution such as a Poisson distribution or a lognormal distribution. The estimation method ID 214 identifies a parameter estimation method such as a least square method or a maximum likelihood method with respect to the distribution type ID 213. The parameters estimated by these are stored in the estimation parameter 215.

ここで、分布種別が複数の分布の混合分布を表している場合は、推定パラメータ215には各分布関数とその混合率なども格納されるものとする。分布種別ID213は、事故件数分布と損失額分布のどちらを推定したものかも識別できるようになっている。   Here, when the distribution type represents a mixed distribution of a plurality of distributions, the estimation parameter 215 stores each distribution function and its mixing rate. The distribution type ID 213 can identify which of the accident number distribution and the loss amount distribution is estimated.

図2(d)の事故データ統計量テーブル114は、「データセットID」216、「 統計量ID」217、「統計量」218からなる。図2(e)の推定分布統計量テーブル115は、「推定分布項番」219、「統計量ID」220、「統計量」221からなる。これらは前者の事故データを特定するデータセットID216と、後者の推定分布を特定する推定分布項番219の違い以外は、構造が同一である。統計量ID217、220は、平均や標準偏差、メディアンなどの統計量の種類を識別するIDで、分布推定に用いるものや、推定分布の候補として棄却するかどうかを判断するための検定統計量、推定分布候補の比較などに用いるものが予め決められているものとする。そして、統計量218、221には、統計量ID217、220でそれぞれ指定された統計量を、データセットID216が示す事故データ、又は、推定分布項番219が示す推定分布に対して算出したものを格納する。   The accident data statistic table 114 in FIG. 2D includes a “data set ID” 216, a “statistic ID” 217, and a “statistic” 218. The estimated distribution statistic table 115 in FIG. 2E includes an “estimated distribution item number” 219, a “statistic ID” 220, and a “statistic” 221. These have the same structure except for the difference between the data set ID 216 that specifies the former accident data and the estimated distribution item number 219 that specifies the estimated distribution of the latter. Statistics IDs 217 and 220 are IDs for identifying types of statistics such as average, standard deviation, and median, and are used for distribution estimation, and test statistics for determining whether to reject as a candidate for estimated distribution, It is assumed that what is used for comparison of estimated distribution candidates is determined in advance. The statistics 218 and 221 are obtained by calculating the statistics specified by the statistics IDs 217 and 220 with respect to the accident data indicated by the data set ID 216 or the estimated distribution indicated by the estimated distribution item number 219, respectively. Store.

図2(f)のシミュレーション結果テーブル116は、1回のシミュレーション結果を識別するための「結果項番」222、シミュレーションに用いた事故データの「データセットID」223、シミュレーションでの試行回数、乱数シードなどを表す「シミュレーション条件」224、シミュレーションに用いた「事故件数分布項番」225、「損失額分布項番」226、シミュレーションで得られた試行回数分の損失総額に対する、上位からの比率を表す「確率」227、確率227が示す分位点の損失総額を表す「VaR」228からなる。   The simulation result table 116 in FIG. 2F includes a “result item number” 222 for identifying one simulation result, a “data set ID” 223 of accident data used for the simulation, the number of trials in the simulation, and a random number. “Simulation condition” 224 representing seeds, “accident number distribution item number” 225, “loss amount distribution item number” 226 used in the simulation, and the ratio from the top to the total loss amount of trials obtained in the simulation “Probability” 227 representing the probability, and “VaR” 228 representing the total loss of the quantile indicated by the probability 227.

事故件数分布項番225と損失額分布項番226には、それぞれ推定分布テーブル113の推定分布項番211を格納する。複数のリスクシナリオでシミュレーションを実行する場合には、これらの項番をシナリオ数のリストに拡張すればよい。   In the accident number distribution item number 225 and the loss amount distribution item number 226, the estimated distribution item number 211 of the estimated distribution table 113 is stored. When a simulation is executed with a plurality of risk scenarios, these item numbers may be expanded to a list of scenario numbers.

一般的に、確率227には、90%、99%、99.9%などが用いられることが多いが、これらの値は予め1つ又は複数設定されているものとする。確率227の代わりにシミュレーションの試行番号を格納し、試行回数分の損失総額を全て格納してもよい。この場合、必要に応じて指定された分位点のVaRなどをその都度算出すればよい。例えば、損失総額の平均やメディアンなど、リスクの計量結果としてはあまり重要ではない統計量などを算出することも可能となる。   Generally, 90%, 99%, 99.9%, etc. are often used for the probability 227, but it is assumed that one or more of these values are set in advance. Instead of the probability 227, a simulation trial number may be stored, and the total loss total for the number of trials may be stored. In this case, VaR of the quantile specified as needed may be calculated each time. For example, it is possible to calculate a statistic that is not so important as a risk measurement result, such as an average of the total loss or a median.

図3は、本実施形態の全体処理の一例を表すフロー図である。
まず、事故データテーブル111から分布推定用事故データテーブル112を生成する(ステップ301)。詳細は図4で説明する。
次に、ステップ301で生成された分布推定用事故データテーブル112に対して、各確率分布候補のパラメータを推定し、推定分布テーブル113を作成する(ステップ302)。詳細は図5で説明する。
次に、ステップ302で推定された推定分布テーブル113を用いて、リスク計量処理を実行する(ステップ303)。詳細は図6で説明する。
次に、出力処理を実行し(ステップ304)、処理を終了する。
FIG. 3 is a flowchart showing an example of the overall processing of this embodiment.
First, a distribution estimation accident data table 112 is generated from the accident data table 111 (step 301). Details will be described with reference to FIG.
Next, parameters of each probability distribution candidate are estimated with respect to the distribution estimation accident data table 112 generated in step 301, and an estimated distribution table 113 is created (step 302). Details will be described with reference to FIG.
Next, risk metric processing is executed using the estimated distribution table 113 estimated in step 302 (step 303). Details will be described with reference to FIG.
Next, an output process is executed (step 304), and the process ends.

ステップ304の出力処理では、ステップ301で生成された分布推定用事故データテーブル112とステップ302で推定された各推定分布との当てはまり具合を評価するためのグラフや、後述のステップ408、ステップ503で算出した統計量、ステップ303で得られたシミュレーション結果によるVaRなどを一覧表示したり、必要に応じて印刷する。グラフの表示方法の例は、図7で説明する。   In the output process of step 304, a graph for evaluating the degree of fit between the distribution estimation accident data table 112 generated in step 301 and each estimated distribution estimated in step 302, and in steps 408 and 503 described later, The calculated statistics, VaR based on the simulation result obtained in step 303, etc. are displayed as a list or printed as necessary. An example of a graph display method will be described with reference to FIG.

これら一連の処理は、ステップ301で新たに生成された分布推定用事故データのデータセットに対して行われる。また、過去に生成、推定又は計量されたデータを用いる場合は、対応するステップを省略し、データセットIDや結果項番などを入力して出力処理を行うことが可能であるものとする。なお、推定分布の候補に対して、出力処理によって分布の当てはまり具合を評価した後、採用された分布のみリスク計量処理を行うように処理の順序を変更することは容易である。   A series of these processes is performed on the data set of accident data for distribution estimation newly generated in Step 301. In addition, when using data generated, estimated, or measured in the past, it is possible to omit the corresponding step and input the data set ID, result item number, and the like to perform output processing. Note that it is easy to change the processing order so that the risk metric processing is performed only on the adopted distribution after evaluating the distribution fit to the estimated distribution candidate by the output processing.

図4は、図3のステップ301の分布推定用事故データテーブル生成処理の一例を表すフロー図である。
まず、推定条件を入力する(ステップ401)。このとき、推定条件が予め設定されていればその条件を利用することができる。例えば、現在蓄積されている事故データから、過去5年分のデータをリスク計量に利用することが決まっていれば、条件の入力を省略することができる。
FIG. 4 is a flowchart showing an example of the distribution estimation accident data table generation process in step 301 of FIG.
First, an estimation condition is input (step 401). At this time, if an estimation condition is set in advance, the condition can be used. For example, if it is determined from the currently accumulated accident data that data for the past five years will be used for risk measurement, the input of conditions can be omitted.

次に、事故データテーブル111から、推定条件を満たす事故データ項番201のリストを取得する(ステップ402)。次に、リストを損失額の昇順にソートする(ステップ403)。次に、事故データ項番リストを分布推定用事故データテーブル112に格納する(404)。このとき、データセットID205には新規に生成されたIDを格納し、推定条件206にはステップ401での推定条件を格納する。   Next, a list of accident data item numbers 201 satisfying the estimation condition is acquired from the accident data table 111 (step 402). Next, the list is sorted in ascending order of the loss amount (step 403). Next, the accident data item number list is stored in the accident data table 112 for distribution estimation (404). At this time, the newly generated ID is stored in the data set ID 205, and the estimation condition in step 401 is stored in the estimation condition 206.

次に、各事故データの事故件数比率を算出し、事故件数比率208に格納する(ステップ405)。
次に、各事故データの累積損失額を算出し、累積損失額209に格納する(ステップ406)。
次に、累積損失額209と事故データの損失総額との比率を算出し、累積損失比率210に格納する(ステップ407)。
Next, the accident number ratio of each accident data is calculated and stored in the accident number ratio 208 (step 405).
Next, the cumulative loss amount of each accident data is calculated and stored in the cumulative loss amount 209 (step 406).
Next, the ratio between the cumulative loss 209 and the total loss of accident data is calculated and stored in the cumulative loss ratio 210 (step 407).

次に、分布推定用事故データテーブル112に対して、予め設定されている各種統計量を算出し、事故データ統計量テーブル114に格納し(ステップ408)、処理を終了する。
ステップ408で算出する統計量には、検定統計量の有位水準ごとの値などを格納しておいてもよい。例えばコルモゴロフ・スミルノフ検定では、データ数によって基準が変化するため、このテーブルに有意水準1%や5%の検定統計量の値を格納することができる。
Next, various preset statistics are calculated for the accident estimation data table 112 for distribution estimation, stored in the accident data statistics table 114 (step 408), and the process ends.
The statistic calculated in step 408 may store a value for each significant level of the test statistic. For example, in the Kolmogorov-Smirnov test, the standard changes depending on the number of data, and therefore, the value of the test statistic at the significance level of 1% or 5% can be stored in this table.

図5は、図3のステップ302の分布推定処理の一例を表すフロー図である。
まず、事故件数分布の推定を行う(ステップ501)。ここでは、分布推定用事故データテーブル112の単位期間当たりの事故件数をもとに、指定された各理論分布と推定法に従って、分布のパラメータを推定する。推定結果は推定分布テーブル113に格納される。
FIG. 5 is a flowchart showing an example of the distribution estimation process in step 302 of FIG.
First, the accident number distribution is estimated (step 501). Here, based on the number of accidents per unit period in the accident data table 112 for distribution estimation, the parameters of the distribution are estimated according to each specified theoretical distribution and estimation method. The estimation result is stored in the estimated distribution table 113.

次に、損失額分布の推定を行う(ステップ502)。ここでは、分布推定用事故データテーブル112の損失額をもとに、指定された各理論分布と推定法に従って、分布のパラメータを推定する。推定結果は推定分布テーブル113に格納される。   Next, the loss amount distribution is estimated (step 502). Here, based on the loss amount in the distribution estimation accident data table 112, the parameters of the distribution are estimated according to each specified theoretical distribution and estimation method. The estimation result is stored in the estimated distribution table 113.

次に、ステップ501と502でそれぞれ推定された各分布に対して、各種統計量を算出し、推定分布統計量テーブル115に格納し(503)、処理を終了する。これは、事故データ統計量テーブル114と対応し、出力処理304で出力することで、推定分布と事故データとの当てはまり具合を評価するのに利用する。   Next, various statistics are calculated for each distribution estimated in steps 501 and 502, stored in the estimated distribution statistics table 115 (503), and the process is terminated. This corresponds to the accident data statistic table 114 and is output by the output process 304, and is used to evaluate the degree of fit between the estimated distribution and the accident data.

図6は、図3のステップ303のリスク計量処理の一例を表すフロー図である。
ここでは、分布推定処理302で推定された各推定分布によるモンテカルロ・シミュレーションを実行し、リスク量を計量する。
FIG. 6 is a flowchart showing an example of the risk measurement process in step 303 of FIG.
Here, a Monte Carlo simulation based on each estimated distribution estimated in the distribution estimation process 302 is executed, and the risk amount is measured.

まず、シミュレーション条件を入力する(ステップ601)。シミュレーション条件には、試行回数や乱数シードなどがある。その他に、出力処理304による分布当てはまり具合の評価などに基づいて、シミュレーションを行う特定の事故件数分布と損失額分布の組み合わせを候補から除外するような処理を前もって実行した場合には、当該除外する分布の組み合わせリストもシミュレーション条件として入力し、指定された推定分布の組み合わせのみシミュレーションを実行することができる。   First, simulation conditions are input (step 601). Simulation conditions include the number of trials and random number seed. In addition, if processing that excludes a combination of a specific accident number distribution and loss distribution to be simulated from candidates based on an evaluation of distribution fit by the output process 304 is excluded in advance A distribution combination list is also input as a simulation condition, and only a specified combination of estimated distributions can be simulated.

次に、推定された事故件数分布と損失額分布の推定分布項番とシミュレーション条件を組み合わせた情報に対するリストPを生成する(ステップ602)。ここでは、シミュレーションを行うデータセットIDに対して、推定された事故件数分布と損失額分布の全ての組み合わせを生成する。ただし、ステップ601において、それらの組み合わせも含めて入力された場合には、この処理は省略可能である。   Next, a list P is generated for information obtained by combining the estimated accident number distribution, the estimated distribution item number of the loss amount distribution, and the simulation conditions (step 602). Here, all combinations of the estimated accident number distribution and loss distribution are generated for the data set ID to be simulated. However, this processing can be omitted if the combination including those combinations is input in step 601.

次に、リストPの要素を表す変数sに"1"を代入する(ステップ603)。次に、sがリストPの要素数以下である間(ステップ604)、以下の処理をループする。   Next, “1” is substituted into a variable s representing an element of the list P (step 603). Next, while s is less than or equal to the number of elements in the list P (step 604), the following processing is looped.

まず、シミュレーション条件P[s]に従って初期化する(ステップ605)。ここでは、シミュレーションの試行回数や乱数シードなどの初期化とともに、試行回数分の各損失総額を格納する配列Lを用意し、各要素を0で初期化する。   First, initialization is performed according to the simulation condition P [s] (step 605). Here, along with initialization of the number of simulation trials, random number seeds, etc., an array L for storing the total loss amount for the number of trials is prepared, and each element is initialized with zero.

次に現在の試行回数を表す変数iに"1"を代入する(ステップ606)。次に、iがP[s]で指定された試行回数以下の間(ステップ607)、以下の処理をループする。   Next, “1” is substituted into a variable i representing the current number of trials (step 606). Next, while i is equal to or smaller than the number of trials specified by P [s] (step 607), the following processing is looped.

まず、P[s]の事故件数分布に従う乱数を1つ発生させる(ステップ608)。
次に、P[s]の損失額分布に従う乱数を、事故件数の個数だけ発生させ、L[i]に加算する(ステップ609)。次に、iに1を加えて(ステップ610)、ステップ607に戻る。
ステップ607で条件がNoなら、まず、配列Lを昇順にソートし(ステップ611)、VaRを算出する(ステップ612)。予め指定された分位点のVaRを算出する。補間してもよく、又は、配列Lをそのまま格納してもいい。
First, one random number according to the accident number distribution of P [s] is generated (step 608).
Next, random numbers according to the loss distribution of P [s] are generated by the number of accidents and added to L [i] (step 609). Next, 1 is added to i (step 610), and the process returns to step 607.
If the condition is No in step 607, first, the array L is sorted in ascending order (step 611), and VaR is calculated (step 612). Calculate VaR of quantile specified in advance. Interpolation may be performed, or the array L may be stored as it is.

次に、条件P[s]のシミュレーション結果をシミュレーション結果テーブル116に格納する(ステップ613)。次に、sに1を加えて(ステップ614)、ステップ604に戻る。
ステップ604での条件がNoなら、処理を終了する。
Next, the simulation result of the condition P [s] is stored in the simulation result table 116 (step 613). Next, 1 is added to s (step 614), and the process returns to step 604.
If the condition in step 604 is No, the process ends.

図7は、本実施形態の出力処理で表示するグラフの表示方法の一例を表す図である。
推定分布と事故データとの当てはまり具合を評価するためのグラフとしては、事故データについては昇順にソートされた損失額とその事故件数比率とをプロットし、そこに推定分布の累積分布関数グラフを対応させて表示するのが、一般的に利用されている方法である。推定分布に基づくグラフと事故データに基づくグラフは、これらを並べて又は重ねて表示することができる。
FIG. 7 is a diagram illustrating an example of a display method of a graph displayed in the output process of the present embodiment.
As a graph to evaluate the fit between the estimated distribution and the accident data, the accident data is plotted with the amount of loss sorted in ascending order and the number of accidents, and the cumulative distribution function graph of the estimated distribution corresponds to it. It is a commonly used method to display the image. The graph based on the estimated distribution and the graph based on the accident data can be displayed side by side or superimposed.

図7(a)では縦横の軸を入れ替えているので、事故データは、損失額の小さい順に事故件数比率とその損失額でプロットしたもので、推定分布は累積分布逆関数となっている。   Since the vertical and horizontal axes are interchanged in FIG. 7A, the accident data is plotted with the accident number ratio and the loss amount in ascending order of the loss amount, and the estimated distribution is an inverse cumulative distribution function.

また、図7(a)のグラフの横軸の変数を媒介変数として、事故データと推定分布それぞれの2点の組を使ってプロットしたものが、図7(b)のq−qプロットである。q−qプロットを使うと、プロットされた点が斜め45度の点線に近いほど当てはまりがよいことがわかる。しかしq−qプロットでは、事故データの最高損失額に対しては、推定分布が無限大になるためプロットすることはできない。このため、最高損失への折れ線部分は点線にしてある。   Also, the plot of the q-q plot of FIG. 7B is a plot using two sets of the accident data and the estimated distribution using the variable on the horizontal axis of the graph of FIG. 7A as a parameter. . Using the qq plot, it can be seen that the closer the plotted point is to the dotted line of 45 degrees, the better the fit. However, the q-q plot cannot be plotted because the estimated distribution becomes infinite with respect to the maximum loss amount of accident data. For this reason, the broken line portion to the maximum loss is a dotted line.

本実施形態では、上記のプロット方法を拡張し、事故データに対して昇順にソートされた損失額に対する事故件数比率と、累積損失額とをプロットし、そこに推定分布の累積分布逆関数の積分関数グラフを対応させて表示する。すなわち、図7(a)のグラフでの変数pの値に対して、確率"0"から該当確率までの面積をプロットしたものが図7(c)のグラフである。累積分布関数をF(x)とすると、累積分布逆関数の積分関数L(p)は、確率"0"から"p"までの積分として、   In the present embodiment, the above plotting method is extended to plot the ratio of the number of accidents against the amount of loss sorted in ascending order with respect to the accident data, and the cumulative loss amount, and integrate the inverse distribution function of the cumulative distribution of the estimated distribution there. Display function graphs correspondingly. That is, the graph of FIG. 7C is obtained by plotting the area from the probability “0” to the corresponding probability with respect to the value of the variable p in the graph of FIG. If the cumulative distribution function is F (x), the integral function L (p) of the inverse cumulative distribution function is the integral from the probability “0” to “p”,

Figure 0006192432
Figure 0006192432

となる。累積分布逆関数の積分は、初等関数の四則演算などに変形できない場合は、既存の数値積分法を用いて近似値を算出すればよい。 It becomes. If the integration of the inverse cumulative distribution function cannot be transformed into four basic arithmetic operations, an approximate value may be calculated using an existing numerical integration method.

さらに、図7(c)のグラフに対して、事故データと推定分布の横軸の変数を媒介変数として、上記のq−qプロットのように事故データと推定分布の値の組を縦軸と横軸にプロットしたものが図7(d)のグラフである。   Furthermore, with respect to the graph of FIG. 7 (c), the variable of the horizontal axis of the accident data and the estimated distribution is a parameter, and the set of the value of the accident data and the estimated distribution is the vertical axis as in the above q-q plot. What is plotted on the horizontal axis is the graph of FIG.

また、推定分布の累積分布逆関数の積分関数を、損失額xの関数として   Also, the integral function of the cumulative distribution inverse function of the estimated distribution is used as a function of the loss amount x.

Figure 0006192432
Figure 0006192432

としてもよい。さらに、縦軸は、積分関数の対数値、又は、確率1のときの積分関数の値との比率とすることも可能である。 It is good. Further, the vertical axis can be a logarithmic value of the integral function or a ratio with the value of the integral function when the probability is 1.

事故データからは、損失額(横軸)と累積損失額(縦軸)のグラフをプロットする。さらに、縦軸の累積損失額を対数値にしてプロットすることも可能である。さらに、縦軸を累積損失額の損失総額に対する比率としてプロットすることも可能である。
事故データに対するグラフの横軸と、推定分布の累積分布逆関数の積分関数に対するグラフの横軸を、1件当たりの事故の損失額として表示してもよい。
さらに、事故データと推定分布の横軸を媒介変数として、事故データと推定分布の値の組を縦軸と横軸にとったグラフをプロットすることが可能である。
From the accident data, a graph of loss (horizontal axis) and cumulative loss (vertical axis) is plotted. Furthermore, it is possible to plot the cumulative loss amount on the vertical axis as a logarithmic value. Further, it is possible to plot the vertical axis as a ratio of the accumulated loss amount to the total loss amount.
The horizontal axis of the graph for the accident data and the horizontal axis of the graph for the integral function of the cumulative distribution inverse function of the estimated distribution may be displayed as the amount of accident loss per case.
Furthermore, it is possible to plot a graph in which the horizontal axis of the accident data and the estimated distribution is a parameter, and the set of the values of the accident data and the estimated distribution is the vertical axis and the horizontal axis.

また、これらのプロット法は、基準となる推定分布候補を1つ選び、その他の推定分布を複数重ねて表示することも可能である。   In addition, in these plotting methods, it is also possible to select one candidate estimated distribution and display a plurality of other estimated distributions in a superimposed manner.

101:入力装置
102:出力装置
103:表示装置
104:処理装置
105:プログラム
106:推定用事故データ生成部
107:分布推定部
108:リスク計量部
109:出力処理部
110:DB
111:事故データテーブル
112:分布推定用事故データテーブル
113:推定分布テーブル
114:事故データ統計量テーブル
115:推定分布統計量テーブル
116:シミュレーション結果テーブル
101: input device 102: output device 103: display device 104: processing device 105: program 106: estimation accident data generation unit 107: distribution estimation unit 108: risk weighing unit 109: output processing unit 110: DB
111: Accident data table 112: Accident data table for distribution estimation 113: Estimated distribution table 114: Accident data statistic table 115: Estimated distribution statistic table 116: Simulation result table

Claims (4)

事故発生日時と損失金額を含む収集された事故データから所定の推定条件により選択された分布推定用事故データに対して、指定された確率分布関数及び指定されたパラメータ推定法を適用して事故件数分布と損失額分布を推定する分布推定手段と、
推定された前記事故件数分布と前記損失額分布に従う乱数を発生させて一定期間内に発生する損失総額を算出するシミュレーションを複数回実行するシミュレーション手段と、
前記シミュレーション手段によって得られた損失総額データにおける指定された分位点の損失総額を算出するリスク計量手段と、を備えた処理装置を有するリスク計量システムにおいて、前記処理装置は、
前記分布推定手段において事故件数分布と損失額分布を表す複数の確率分布関数及び複数のパラメータ推定法により推定された各推定分布の候補について、
当該推定分布と前記分布推定用事故データとの当てはまり具合の評価のために、当該推定分布の累積分布逆関数の積分関数と、前記分布推定用事故データに基づいた損失額と一定損失額以下の累積事故件数の全事故件数に対する比率と一定損失額以下の累積損失額とにより表されるデータとをグラフに表示する推定分布評価手段を有することを特徴とする
リスク計量システム。
The number of accidents by applying the specified probability distribution function and the specified parameter estimation method to the accident data for distribution estimation selected from the collected accident data including the date and time of the accident and the amount of loss according to the specified estimation conditions Distribution estimation means for estimating distribution and loss distribution;
Simulation means for generating a random number according to the estimated accident number distribution and the loss amount distribution to calculate a total loss occurring within a predetermined period, a simulation means for executing a plurality of times;
In a risk weighing system comprising: a risk weighing unit that calculates a total loss of a specified quantile in the total loss data obtained by the simulation unit, the processing device includes:
For each estimated distribution candidate estimated by a plurality of probability distribution functions and a plurality of parameter estimation methods representing the accident number distribution and loss distribution in the distribution estimation means,
In order to evaluate the degree of fit between the estimated distribution and the distribution estimation accident data, the integral function of the cumulative distribution inverse function of the estimated distribution, the loss amount based on the distribution estimation accident data, and a certain loss amount or less A risk weighing system comprising: an estimated distribution evaluation unit that displays a ratio of the total number of accidents to the total number of accidents and data represented by a cumulative loss amount equal to or less than a certain loss amount on a graph.
前記推定分布評価手段において、
前記分布推定用事故データについては、横軸を一定損失額以下の累積事故件数の全事故件数に対する比率とし、縦軸を累積損失額若しくはその対数値、又は累積損失額の損失総額に対する比率としたグラフを表示し、
前記推定分布の累積分布逆関数の積分関数については、横軸を確率とし、縦軸を前記積分関数の値若しくはその対数値、又は確率1のときの前記積分関数の値に対する比率としたグラフを表示し、
これらのグラフを並べて表示又は重ね合わせて表示することを特徴とする
請求項1に記載のリスク計量システム。
In the estimated distribution evaluation means,
For the accident data for distribution estimation, the horizontal axis is the ratio of the total number of accidents with a certain loss amount or less to the total number of accidents, and the vertical axis is the cumulative loss amount or its logarithm, or the ratio of the cumulative loss amount to the total loss Display the graph,
For the integral function of the cumulative distribution inverse function of the estimated distribution, a graph in which the horizontal axis is the probability and the vertical axis is the value of the integral function or its logarithmic value, or the ratio to the value of the integral function when the probability is 1. Display
The risk weighing system according to claim 1, wherein these graphs are displayed side by side or displayed in a superimposed manner.
前記推定分布評価手段において、
事故データに対するグラフの横軸と、推定分布の累積分布逆関数の積分関数に対するグラフの横軸を1件当たりの事故の損失額として表示することを特徴とする
請求項2に記載のリスク計量システム。
In the estimated distribution evaluation means,
The risk measurement system according to claim 2, wherein the horizontal axis of the graph for accident data and the horizontal axis of the graph for the integral function of the cumulative distribution inverse function of the estimated distribution are displayed as an accident loss amount per case. .
前記推定分布評価手段において、
事故データと推定分布のそれぞれのグラフを重ね合わせたときの、横軸の変数を媒介変数として、事故データと推定分布のそれぞれのグラフ上での値の組を縦軸と横軸にとったグラフを表示することを特徴とする
請求項2又は3に記載のリスク計量システム。
In the estimated distribution evaluation means,
A graph in which the horizontal axis variable is the parameter when the accident data and estimated distribution graphs are superimposed, and the vertical and horizontal axes are pairs of values on the accident data and estimated distribution graphs. The risk weighing system according to claim 2, wherein the risk weighing system is displayed.
JP2013172162A 2013-08-22 2013-08-22 Risk weighing system Active JP6192432B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013172162A JP6192432B2 (en) 2013-08-22 2013-08-22 Risk weighing system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013172162A JP6192432B2 (en) 2013-08-22 2013-08-22 Risk weighing system

Publications (2)

Publication Number Publication Date
JP2015041259A JP2015041259A (en) 2015-03-02
JP6192432B2 true JP6192432B2 (en) 2017-09-06

Family

ID=52695382

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013172162A Active JP6192432B2 (en) 2013-08-22 2013-08-22 Risk weighing system

Country Status (1)

Country Link
JP (1) JP6192432B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7207762B2 (en) 2017-06-20 2023-01-18 ヴァーダーグ パイプ テクノロジー リミテッド How to design pipe fittings for use in subsea pipelines

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6712921B2 (en) * 2016-07-21 2020-06-24 株式会社竹中工務店 Crime damage amount estimation device and program
CN107862470A (en) * 2017-11-27 2018-03-30 国网北京市电力公司 Determination method and device, the storage medium of accident risk grade
KR102234490B1 (en) * 2018-12-20 2021-03-31 대한민국 Estimation of insect developmental time using skewed model algorithm

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07230449A (en) * 1994-02-17 1995-08-29 Nippon Telegr & Teleph Corp <Ntt> Data prediction device
JP2003036346A (en) * 2001-07-23 2003-02-07 Mitsubishi Trust & Banking Corp Method for evaluating operational risk and its system
JP2005056024A (en) * 2003-07-31 2005-03-03 Hitachi East Japan Solutions Ltd Support system and support method for creating supply plan, and supporting computer program
JP5725547B2 (en) * 2011-03-29 2015-05-27 日本電気株式会社 Risk management device
JP5867349B2 (en) * 2011-09-21 2016-02-24 新日鐵住金株式会社 Quality prediction apparatus, operation condition determination method, quality prediction method, computer program, and computer-readable storage medium

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7207762B2 (en) 2017-06-20 2023-01-18 ヴァーダーグ パイプ テクノロジー リミテッド How to design pipe fittings for use in subsea pipelines

Also Published As

Publication number Publication date
JP2015041259A (en) 2015-03-02

Similar Documents

Publication Publication Date Title
Erdik Earthquake risk assessment
US9772923B2 (en) Fast OLAP for real user measurement of website performance
US9367382B2 (en) Apparatus, method, and program product for calculating abnormality based on degree of correlation destruction
US20150067648A1 (en) Preparing an optimized test suite for testing an application under test in single or multiple environments
JP6192432B2 (en) Risk weighing system
US10451416B1 (en) Optimizing sensor placement for structural health monitoring based on information entropy or total modal energy
US20160132798A1 (en) Service-level agreement analysis
US20190377653A1 (en) Systems and methods for modeling computer resource metrics
WO2019234716A1 (en) Methods and systems for blockchain testing
US9317387B2 (en) Methods and systems for reducing metrics used to monitor resources
WO2021021271A9 (en) Indiagnostics framework for large scale hierarchical time-series forecasting models
Wang An imperfect software debugging model considering irregular fluctuation of fault introduction rate
CN115841046A (en) Accelerated degradation test data processing method and device based on wiener process
CN112613983B (en) Feature screening method and device in machine modeling process and electronic equipment
KR101471797B1 (en) Risk-management device
JP2013036901A (en) Destructive inspection system, destructive inspection method, data processing program, and program recording medium
KR101566601B1 (en) Risk management device
JP5348351B2 (en) Risk profile generator
JP5650290B1 (en) Operational risk measurement method and apparatus
Agarwal Markovian software reliability model for two types of failures with imperfect debugging rate and generation of errors
CN112184415A (en) Data processing method and device, electronic equipment and storage medium
JP5804492B2 (en) Risk management device
US20240005259A1 (en) Index modeling
US20160379147A1 (en) Energy intensity variability analysis
Aranha et al. Model simulation for test execution capacity estimation

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20150223

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160809

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170615

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170725

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170808

R150 Certificate of patent or registration of utility model

Ref document number: 6192432

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250