オペレーショナルリスクの計量化装置、方法、およびプログラム
【課題】少ない損失事例から巨額損失のオペレーショナルリスクを合理的に計量化できるオペレーショナルリスクの計量化装置、方法、およびプログラムを提供すること。
【解決手段】取引金額を読み込む手段と、損失率密度を読み込む手段と、前記取引金額と前記損失率密度とに基づいて巨額損失密度を算出する巨額損失密度算出手段と、前記巨額損失密度算出手段により算出された巨額損失密度からリスク量を算出するリスク量算出手段と、を具備する。
【解決手段】取引金額を読み込む手段と、損失率密度を読み込む手段と、前記取引金額と前記損失率密度とに基づいて巨額損失密度を算出する巨額損失密度算出手段と、前記巨額損失密度算出手段により算出された巨額損失密度からリスク量を算出するリスク量算出手段と、を具備する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、統計学的手法に基づいてオペレーショナルリスクを計量化する装置、方法、およびプログラムに関する。
【背景技術】
【０００２】
新BIS(Bank for International Settlements；国際決済銀行)規制では新しいオペレーショナルリスクが導入される。これまでのBIS規制では信用リスクと市場リスクに対して８％以上の自己資本を準備しておくことが要求されていた。これに対して新BIS規制では、信用リスク、市場リスク、およびオペレーショナルリスクの合計に対して８％の自己資本を準備しておくことが要求される。より正確には、９９．９％水準のVaR（バリューアットリスク）に対応するリスク量に対して８％以上の自己資本を用意しておくことが要求される。
【０００３】
リスク量の評価には、基礎的手法、標準的手法、先進的手法のいずれかを採用してよいことになっているが、前者になるほど精度の低い評価であることから、余裕を持った評価が求められ、リスク量としては大きくならざるを得ないと考えられる。このとき、リスク相応の自己資本を調達しておかなければならない。一般に、自己資本の調達にはコストが掛かるため望ましいことではないし、対外的にはリスクが小さくて自己資本比率が高いということは良い評価に結びつくため、銀行にとって先進的手法を導入しようというインセンティブが存在する。
【０００４】
しかし、先進的手法を採用した場合、オペレーショナルリスクを合理的に計量化する必要が生じるが、過去の損失事故の件数が十分にないために正確な計量化が難しいという問題がある。特に、リスク評価で重要な巨額損失のデータが少ないことが問題である。この課題を解決するために、各種方法が提案されているが、決定的な方法がまだ確定していないのが現状である。オペレーショナルリスクの計量化については、例えば非特許文献１に記載されている。
【非特許文献１】「オペレーショナル・リスクのすべて」，三菱信託銀行オペレーショナル・リスク研究会編，東洋経済新聞社
【発明の開示】
【発明が解決しようとする課題】
【０００５】
本発明は、少ない損失事例から巨額損失のオペレーショナルリスクを合理的に計量化できるオペレーショナルリスクの計量化装置、方法、およびプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の一観点に係るオペレーショナルリスクの計量化装置は、取引金額を入力する手段と、損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を入力する手段と、損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を前記取引金額と前記損失率密度とに基づいて算出する巨額損失密度算出手段と、前記巨額損失密度算出手段により算出された巨額損失密度からオペレーショナルリスクを算出するリスク量算出手段と、を具備する。
【発明の効果】
【０００７】
本発明によれば、少ない損失事例から巨額損失のオペレーショナルリスクを合理的に計量化できるオペレーショナルリスクの計量化装置、方法、およびプログラムを提供できる。
【発明を実施するための最良の形態】
【０００８】
以下、図面を参照しながら本発明の実施形態を説明する。
【０００９】
先ずは、幾つかの仮定の下にオペレーショナルリスク評価のためのモデルを導出する過程を第１の実施形態により説明したのち、同モデルを用いた具体的な定量化に係る実施形態（第２〜第５の実施形態）を説明する。第２及び第３の実施形態では巨額損失のみならず少額損失をも分析可能である。第４及び第５の実施形態では少額損失の分析を行わず処理手順を簡単化したものである。
【００１０】
（第１の実施形態）
次のような確率変数、すなわち取引金額の確率変数をX（ただし、Xはゼロ以上の実数）、損失額の確率変数をY（ただし、Yはゼロ以上の実数）、損失率の確率変数をΘ（ただし、Θゼロ以上１以下の実数）、損失発生の有無：L（０：損失無し、１：損失有り）を考えるとき、損失率Θは次の式(1)のように表すことができる。
Θ=Y/X…(1)
なお、取引金額、損失額は実際にはゼロ以上の整数であるが、説明を簡単にするために実数であるものとして説明する。
【００１１】
始めに、P(Y=y,L=l,X=x)について、以下の関係が成立する。
【００１２】
P(Y=y,L=l,X=x)=P(Y=y|L=l,X=x)P(L=l|X=x)P(X=x)…(2)
また、金額xの取引において損失事故が発生する確率P(L=1|X=x)は、取引金額xの広義単調減少関数であると考えられる。ただし、人間の注意力には限界があるため、ある程度大きな金額になると一定値になると考える。そこで、次の仮定１を置く。
【００１３】
（仮定１）x≧xth1であればP(L=1|X=x)=Pl
このとき、
P(L=1|X=x)=Pl…(3)であり、
P(L=0|X=x)=1-Pl…(4)である。さらに、P(Y=y|L=0,X=x)では、次の仮定２を置く。
【００１４】
（仮定２）損失事故が発生しない場合の損失金額は０である。このとき、
P(Y-y|L=0,X=x)=δ(y)…(5)
である。ここで、δ(y)はディラックのデルタ関数である。また、
【数１】
である。ここで、Θは回収率を表し、P(Θ=θ|L=1,X=x)は、x円の取引で事故が発生したときに回収できる割合がθである確率密度を表している。ここでは、一定以上の金額における回収率は取引金額xに依らない密度関数に従うと考えて、次の仮定３を置く。
【００１５】
（仮定３）x≧xth2であればP(Θ=θ|L=1,X=x)=P∞(Θ=θ)
次に、P(Y=y,L=l,X=x)の周辺密度関数P(Y=y,L=l)において、特にl=1の場合を考える。
【００１６】
【数２】
【００１７】
通常、取引金額xの取引でxより大きな損失が出ることは無いので、次の仮定４を置く。
【００１８】
（仮定４）P(Y=y|L=l,X=x)=0 for {(x,y)|x<y}
特に、y≧xth=max(xth1,xth2)であれば、仮定１と仮定３を考慮すると、式(7)より、
【数３】
である。一方、P(Y=y,L=1)=P(Y=y|L=1)P(L=1)…(9)
であるから、
【数４】
である。大数の法則よりサンプル数が十分多いとき、
【数５】
である。ここで、Ωは全サンプルを表す集合である。さらに、yについて積分することで、損失金額が0からyの区間にある確率は、
【数６】
である。ここで、
【数７】
である。また、
【数８】
であるから、サンプル数が十分多いとき、
【数９】
と算出できる。ここで、δj,kはクロネッカのデルタであり、ｊとｋが等しければ１、そうでなければ０である。また、Plは、Ωth={i|xi≧xth}として、仮定１のもとで、
【数１０】
となる。特に、Pl=P(L=1)であれば、
【数１１】
である。
【００１９】
以上の仮定１，２，４の下で巨額損失の分布は式(11)のように表すことができる。これを利用してオペレーショナルリスクの計量化を次の手順により行うことができる。
【００２０】
先ずP∞(Θ=θ)のモデルを作成する（ステップＳ１）。このモデルは、一定以上の取引金額における損失率の確率密度に相当し、後述する実施形態における損失率密度の読込及び記憶の処理手順に対応する。より具体的に説明すると、損失率密度とは、損失率（０〜１の値）を確率変数としたときの確率密度のことである。例えば、１００件の損失事故のなかで、０．５〜０．６（５０％〜６０％）の損失率の事故が１０件だとすると、０．５〜０．６に事象が発生する確率は１０／１００＝０．１となる。さらに、確率密度は区間の幅で割って、０．１／（０．６−０．５）＝１となる。これが損失率密度である。
【００２１】
P∞(Θ=θ)は、Ωthに属するサンプルにより、θについての密度関数をモデル化すればよい。ノンパラメトリック手法（参考文献：「平滑化とノンパラメトリック回帰への招待」J・S・シモノフ 農林統計協会）が有力な方法であるが、ここでは、最も単純なヒストグラムを使った例を説明する。
【００２２】
P∞(Θ=θ)は、0≦θ≦1の範囲でゼロ以上の値を取り得るが、これ以外の領域では０となる。そこで、モデル化は0≦θ≦1におけるP∞(Θ=θ)の関数を決めることである。
【００２３】
以下のようにGk(k=1,2,…,K)を定義する。
【数１２】
である。また、
【数１３】
である。このとき、
【数１４】
として、P∞(Θ=θ)=Pk（θ∈Rkの場合）となる。この場合、P∞(Θ=θ)=Pkの例を図１に示す。この例では、閾値すなわち、k/K(k=0,…,K)の間隔は一定であったが、閾値を任意の値に変化させても良い。最適な閾値を算出する方法については、例えば「情報量統計学」坂本慶行他、共立出版株式会社等に記載されている。また、前出の「平滑化とノンパラメトリック回帰への招待」には、ヒストグラムに代えて頻度ポリゴンを使用したり、カーネル関数を使う方法が記載されている。
【００２４】
次に、
【数１５】
を算出する（ステップＳ２）。ここでは、全取引サンプルの数をNとし、このうちで、事故の発生した数をN1とし、P(L=1)=N1/Nとして算出する。算出される値は、後述する実施形態における、損失事故率の読込及び記憶に対応する。
【００２５】
次に、
【数１６】
を算出する（ステップＳ３）。
【００２６】
xthを外部からうけとり、サンプルの中で取引金額がxth以上のサンプルの数を数えてnとする。さらに、このうちで事故が発生した数をn1とし、Pl=n1/nとして算出する。算出される値は、後述する実施形態における、巨額取引損失事故率の読込及び記憶に対応する。
【００２７】
次に、
【数１７】
を計算する（ステップＳ４）。
【００２８】
ただし、ymaxは取引金額の最大値とし、ymin=xthとする。
【００２９】
ymin≦y≦ymaxの範囲で複数の代表点を決めておく。第１の実施形態では外部から与える数値Mを用いて、
【数１８】
としておく。そして全てのｍについて、以下のステップＳ４−１〜Ｓ４−３の処理を繰り返す。
【００３０】
まず、y=ymをセットする（ステップＳ４−１）。
【００３１】
次に、
【数１９】
を算出する。そして、全ての取引サンプルiについて、θi=y/xiを算出し、上記ステップＳ１で求めた関数P∞(Θ=θ)を使って、
【数２０】
を算出する。算出された値の全てのサンプルについての平均を求める（ステップＳ４−２）。
【００３２】
次に、上記ステップＳ２，Ｓ３，およびＳ４−２の結果を用い、
【数２１】
に従って、pmmaxを算出する（ステップＳ４−３）。
【００３３】
次に、P(Y=y|L=1)をヒストグラムで表す。
【数２２】
では、P(Y=y|L=1)=pmmaxとする。これにより、図２のようなヒストグラムが得られる。ここで、η=(y0+y1)/2である（ステップＳ５）。得られるヒストグラムは、後述する実施形態において算出される巨額損失密度に対応する。ここで、損失額を確率変数としたときの確率密度を「損失密度」と呼ぶ。特に、損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する損失密度のことを「巨額損失密度」と呼ぶ。
【００３４】
次に、ｙが小さいときのモデルを作成する。このモデルは、後述する実施形態における少額条件付損失密度の読込及び記憶に対応する。
【００３５】
上記ステップＳ１と同様に、ヒストグラムを使ってｙがyminよりも小さい範囲のヒストグラムを描く（図３）。
【００３６】
全サンプルから、事故が発生しており、かつ、損失額yがηより小さいサンプルを集めて、Ωminとする。これをU個に分ける以下のような集合を考える。
【数２３】
【００３７】
このとき、事故が発生しており、かつ、損失額yがηより小さいサンプルの数をNminとし、Huに属するサンプルの数をNuminとして、
【数２４】
を算出する（ステップＳ６）。
【００３８】
次に、P(Y=y|L=1)の分布を求める（ステップＳ７）。これは、後述する実施形態における損失密度の算出に対応する。
【００３９】
上記ステップＳ６で求めたpuminは事故が発生しており、かつ、損失額yがηより小さいという条件の下での分布関数であるため、P(Y=y|L=1)を算出するには、事故が発生しており損失額yがηより小さい確率qを乗じる必要がある。これは以下のように求める。
【００４０】
すなわち、
【数２５】
であることから、
【数２６】
とする。
【００４１】
最後に、オペレーショナルリスクの計量化を行う（ステップＳ８）。計量化の手法については別の実施形態において説明するが、公知の手法を適用可能である。
【００４２】
以上説明した第１の実施形態によれば、上述した４つの仮定の下にオペレーショナルリスク評価のためのモデルを導出することができ、同モデルを用い、ステップＳ８においては合理的にオペレーショナルリスクの計量化を行うことが可能になる。尚、ステップＳ８におけるオペレーショナルリスクの計量化には、公知の手法を適用することができる。
【００４３】
（第２の実施形態）
本発明の第２の実施形態では巨額損失のみならず少額損失をも分析可能である。
【００４４】
図４は本発明の第２の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図である。本装置はコンピュータを利用して実現できる。また、本発明に係る機能を実現するためのプログラムが図示しないプログラム記憶装置に格納される。プログラム記憶装置は、例えば、磁気ディスク装置や光ディスク装置等の外部記憶装置からなる。上記プログラムがＣＰＵからの制御でランダムアクセスメモリ（ＲＡＭ）等からなる主記憶装置に読み込まれ、同ＣＰＵにより実行されることにより、コンピュータを本発明に係るオペレーショナルリスク定量化装置として機能させることができる。なお、オペレーショナルリスク定量化装置には、各種コンピュータ資源を管理し、ファイルシステムやネットワーク通信機能、ならびにグラフィカルユーザインタフェース（ＧＵＩ）等を提供するオペレーティングシステムも導入されている。なお、後述する第３乃至第５の実施形態に係るオペレーショナルリスク定量化装置についても、コンピュータを利用して実現できる。
【００４５】
先ず、取引金額読込部４０１により、過去１年間の銀行の各取引における取引金額を読み込む。次の表２７のような形式のデータから、指定された期間（例えば2003年４月１日から2004年３月３１日）の取引を抜き出して、取引金額記憶部４０２に蓄積する。
【００４６】
【表２７】
【００４７】
取引金額記憶部４０２には主記憶装置が利用されるが、データ量が多い場合には外部記憶装置が用いられる。ここでは、取引金額の集合を{xi|i∈Ω}と書くが、Ωは指定された期間に含まれる取引のレコード番号の集合であり、xiはレコード番号iの取引の取引金額である。
【００４８】
次に、損失率密度読込部４０３により損失率密度関数P∞(Θ=θ)を読み込む。損失率密度関数は少なくとも０以上１以下の範囲で定義されており、出力は０以上の値を取る。また、０から１までの区間で積分すると１になる。本実施形態では図５のような損失率密度関数のヒストグラムを考えるが、具体的には次の表２８のような形式でデータを損失率密度記憶部４０４に記憶する。
【００４９】
【表２８】
【００５０】
以下の実施形態の説明において、X〜Yと書いたときは、X以上Y未満を表すものとする。ただし、Y=1の場合だけ、X以上Y以下を意味するものとする。０．２１という値がくれば、この表２８の０．２以上０．３未満の欄をみて、０．３という値を返す。同様に考えて、値０．２が入力されると、０．３を出力し、値０．３の場合は０．２を出力する。また、特に１．０の場合は３．５であるとする。
【００５１】
次に、損失事故率読込部４０８により、全取引に対する損失事故率P(L=1)を読み込で損失事故率記憶部４０９に記憶する。なお、乱数に基づく事故発生モデルを含むオペレーション品質評価モデルから損失事故率を算出してもよい。
【００５２】
次に、巨額取引損失事故率読込部４１０により巨額取引損失事故率（本実施形態では、１００万円以上の取引に対する損失事故率とする）Plを読み込んで巨額取引損失事故率記憶部４１１に記憶する。
【００５３】
次に、少額条件付損失密度読込部４１２は、１００万円未満の損失が発生したという条件付の損失密度を読み込んで、少額条件付損失密度記憶部４１３に記憶する。損失密度は次の表２９のような形で記録される。
【００５４】
【表２９】
【００５５】
損失分布は０円から999,999円までの間で和をとると１となる（損失額０円の確率＋損失額１円の確率＋・・・損失額999,999の確率＝１）。
【００５６】
次に、巨額損失密度算出部４０５は、取引金額記憶部４０２から取引金額{xi|i∈Ω}、損失率密度記憶部４０４から損失密度関数P∞(Θ=θ)の情報を受け取り、損失事故率記憶部４０９から損失事故率P(L=1)、巨額取引損失事故率記憶部４１１から巨額取引損失事故率Plを受取り、
【数３０】
を計算する。この計算はｙを変化させながら繰り返し実行する。次の表３１は１００万円毎にこの計算を実行した結果をまとめた例である。
【００５７】
【表３１】
【００５８】
ｙを１円毎変化させて計算を実行してもかまわないが、その場合は計算量が膨大になるので、本実施形態では１００万円ごとに計算を行うことする(例えばy=1.5百万円の場合の数値を１〜２の欄に記載する)。
【００５９】
次に、損失密度算出部４１４は、巨額損失密度算出部４０５で算出した巨額損失密度と少額条件付損失密度記憶部４１３に記憶された少額条件付損失密度から損失密度関数を算出する。具体的には、１００万円以上の損失密度は巨額損失密度算出部４０５で算出した巨額損失密度とする。巨額損失密度算出部４０５で算出した巨額損失密度に基づき、１００万円以上の損失が発生する確率は以下の式に従って算出できる。
【数３２】
【００６０】
１００万円未満の損失密度については、少額条件付損失密度記憶部４１３に記憶された少額条件付損失密度に1-ηを乗じた値とする。η=0.01であったとすると、次の表３３のような１００万円未満の損失密度を得る。
【００６１】
【表３３】
【００６２】
これを、１００万円以上の巨額損失密度とあわせると、次の表３４を得る。
【００６３】
【表３４】
【００６４】
さらに、損失累積確率として次の表３５を得る。
【００６５】
【表３５】
【００６６】
次に、損失事故率記憶部４０９に記憶した損失事故率と損失密度算出部４１４で算出した損失累積確率とを用いてリスク量算出部４０６は、VaR等のリスクおよび損失金額の期待値をオペレーショナルリスクとして算出する。VaRについては、例えば参考文献（「金融リスクの計量化 上 バリューアットリスク」，木島正明編著）に詳しい。ここで、リスク量算出部４０６により算出されるオペレーショナルリスクは、将来の取引に対するリスク量のことをいう。損失額が大方の予想に反して非常に大きくなることを「リスク」と呼ぶ。オペレーショナルリスクは、所定期間の損失額合計を横軸とし、個々の損失額合計値に対する確率密度を縦軸とする損失額合計密度関数において、所定の上側のパーセンタイル値として求めることができる。あるいは、所定期間の損失額合計を横軸とし、個々の損失額合計値に対する確率密度を縦軸とする損失額合計密度関数において、所定の上側のパーセンタイル値以上の損失額合計値の平均値として求めてもよい。
【００６７】
より具体的に、リスク量算出のための処理手順の一例を説明する。はじめに、λ=P(L=1)とし、パラメータλのポアッソン分布を考える。
【数３６】
【００６８】
λ＝１０の場合のpkの例を次の表３７に示す。
【００６９】
【表３７】
【００７０】
表３７の累積確率は、
【数３８】
であり、たとえば、ｋ＝２であれば、事故件数が２件以下である確率を意味する。
【００７１】
パラメータλにしたがって、予めこのようなテーブルを用意しておく。
次に、一様乱数発生器により０から１までの実数をランダムに発生する。
次に、乱数の値をｘとしたときに、上記テーブルをもとに、事故件数を決定する。表３７では、4.53999E-05未満であれば事故は０件、4.53999E-05以上0.000499399未満であれば事故は１件、・・・とする。たとえば、ｘ＝０．８であれば、0.791556476以上、0.864464423未満であるから、事故件数は１３件とする。
【００７２】
次に、上記損失事故のそれぞれについての損失金額をセットする。事故件数回一様乱数を発生させ、表３５を用いて損失額を確定する。たとえば、乱数の値が０．５であれば、0.4851以上0.5544未満であるから、損失金額は７円となる。１３回事故が発生すれば、それぞれについて同様の処理により損失金額を定める。
【００７３】
次に、事故回数分の損失金額を合算する。たとえば、１３回分の損失金額を合計して、例えば15,568を得る。
【００７４】
以上の処理を指定された回数（N回）繰り返す。リスク量出力部４０７は、それぞれのラウンドで算出される損失金額合計を大きい順に並べて、上位０．１％の値を９９．９％水準のVaRとして出力する。また、損失金額の平均も求めて期待損失額として出力する。たとえば、N＝１００万回であるとすると、上位０．１％は上位１０００番目の損失金額に対応している。また、損失額の上位０．１％の計算結果について抜き出し、抜き出したサンプルの損失額の期待値を計算してこれをCVaR（コンディショナルバリューアットリスク）として出力する。
【００７５】
（第３の実施形態）
図６は本発明の第３の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図である。第３の実施形態が第２の実施形態と相違する点は、損失の確率密度を読み込むのではなく、損失の累積確率を読み込む点である。つまり、図６における損失率累積確率読込部６０３および少額条件付損失累積確率読込部６１２が第２の実施形態（図４）との違いである。この結果、引き続き行われる巨額損失累積確率算出部６０５、損失累積確率算出部６１４の処理手順も第２の実施形態とは異なる。
【００７６】
先ず、取引金額読込部６０１により、過去１年間の銀行の各取引における取引金額を読み込む。次の表３９のような形式のデータから、指定された期間（例えば2003年４月１日から2004年３月３１日）の取引を抜き出して、取引金額記憶部６０２に蓄積する。
【００７７】
【表３９】
【００７８】
取引金額記憶部６０２には主記憶装置が利用されるが、データ量が多い場合には外部記憶装置が用いられる。ここでは、取引金額の集合を{xi|i∈Ω}と書くが、Ωは指定された期間に含まれる取引のレコード番号の集合であり、xiはレコード番号iの取引の取引金額である。
【００７９】
次に、損失率累積確率読込部６０３により、損失率累積確率すなわち、
【数４０】
を読み込む。損失率累積確率は少なくとも０以上１以下の範囲で定義されており、出力も０以上１以下の値を取る。本実施形態では次の表４１のような形式でデータを損失率累積確率関数記憶部６０４に記憶する。
【００８０】
【表４１】
【００８１】
次に、損失事故率読込部６０８により、全取引に対する損失事故率P(L=1)を読み込んで損失事故率記憶部６０９に記憶する。
【００８２】
次に、巨額取引損失事故率読込部６１０により巨額取引損失事故率（本実施形態では、１００万円以上の取引に対する損失事故率とする）Plを読み込んで巨額取引損失事故記憶部６１１に記憶する。
【００８３】
次に、少額条件付損失累積確率読込部６１２は、１００万円未満の損失密度を読み込み、少額条件付損失累積確率記憶部６１３に記憶する。損失累積確率は次の表４２のような形で記録される。
【００８４】
【表４２】
【００８５】
次に、巨額損失累積確率算出部６０５は、取引金額記憶部６０２から取引金額{xi|i∈Ω}、損失率累積確率関数記憶部６０４から損失率累積確率関数F∞(Θ=θ)の情報を受け取り、損失事故率記憶部６０９から損失事故率P(L=1)、巨額取引損失事故記憶部６１１から巨額取引損失事故率Plを受取り、
【数４３】
を計算する。この計算はｙを変化させながら繰り返し実行する。次の表４４は１００万円毎にこの計算を実行した結果をまとめた例である。
【００８６】
【表４４】
【００８７】
ｙを１円毎変化させて計算を実行してもかまわないが、その場合は計算量が膨大になるので、本実施形態では１００万円ごとに計算を行うこととし、その間の値は直線で補完するものとする（もちろん、二次曲線等の曲線で近似しても良い）。
【００８８】
次に、損失累積確率算出部６１４では、巨額損失累積確率算出部６０５で算出した巨額損失累積確率と少額条件付損失累積確率記憶部６１３に記憶された少額条件付損失累積確率から損失累積確率を算出する。具体的は、１００万円以上の損失累積確率は巨額損失累積確率算出部６０５で算出した巨額損失累積確率とする。巨額損失累積確率算出部６０５で算出した巨額損失累積確率をつかい、１００万円未満の損失が発生する確率は以下の式に従って算出できる。
【数４５】
【００８９】
１００万円未満の損失累積確率については、少額条件付損失累積確率記憶部６１３に記憶された少額条件付損失密度に1-ηを乗じた値とする。
【００９０】
η=0.01であったとすると、次の表４６のような１００万円未満の損失累積確率を得る。
【００９１】
【表４６】
【００９２】
１００万円以上とあわせて、次の表４７を得る。
【００９３】
【表４７】
【００９４】
次に、損失事故率記憶部６０９に記憶した損失事故率と損失累積確率算出部６１４で算出した損失累積確率を用いて、リスク量算出部６０６はVaR等のリスクおよび損失金額の期待値を算出する。リスク量算出のための処理手順の一例は次の通りである。
【００９５】
はじめに、λ=P(L=1)とし、パラメータλのポアッソン分布を考える。
【数４８】
【００９６】
λ＝１０の場合のpkの例を次の表４９に示す。
【００９７】
【表４９】
【００９８】
表４９の累積確率は、
【数５０】
であり、たとえば、ｋ＝２であれば、事故件数が２件以下である確率を意味する。
【００９９】
パラメータλにしたがって、予めこのようなテーブルを用意しておく。
【０１００】
次に、一様乱数発生器により０から１までの実数をランダムに発生する。
【０１０１】
次に、乱数の値をｘとしたときに、上記テーブルをもとに、事故件数を決定する。表４９では、4.53999E-05未満であれば事故は０件、4.53999E-05以上0.000499399未満であれば事故は１件、・・・とする。たとえば、ｘ＝０．８であれば、0.791556476以上、0.864464423未満であるから、事故件数は１３件とする。
【０１０２】
次に、上記損失事故のそれぞれについての損失金額をセットする。事故件数回一様乱数を発生させ、表４７を用いて損失額を確定する。たとえば、乱数の値が０．５であれば、0.693以下であるから、損失金額は０円から１０円の間にあることが分かる。（０，０）→（10,0.693）の間を直線で補間すると、
【数５１】
となり、損失額は７円とする。１３回事故が発生すれば、それぞれについて同様の処理により損失金額を定める。
【０１０３】
次に、事故回数分の損失金額を合算する。たとえば、１３回分の損失金額を合計して、15,568を得る。
【０１０４】
以上の処理を指定された回数（N回）繰り返す。リスク量出力部６０７は、それぞれのラウンドで算出される損失金額合計を大きい順に並べて、上位０．１％の値を９９．９％水準のVaRとして出力する。また、損失金額の平均も求めて期待損失額として出力する。たとえば、N＝１００万回であるとすると、上位０．１％は上位１０００番目の損失金額に対応している。また、損失額の上位０．１％の計算結果について抜き出し、抜き出したサンプルの損失額の期待値を計算してこれをCVaRとして出力する。
【０１０５】
（第４の実施形態）
第４の実施形態は少額損失の分析を行わず処理手順を簡単化したものである。図７は本発明の第４の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図である。取引金額読込部７０１により、過去１年間の銀行の各取引における取引金額を読み込む。次の表５２のような形式のデータから、指定された期間（例えば２００３年４月１日から２００４年３月３１日）の取引を抜き出して、取引金額記憶部７０２に蓄積する。
【０１０６】
【表５２】
【０１０７】
取引金額記憶部７０２には主記憶装置が利用されるが、データ量が多い場合には外部記憶装置が用いられる。取引金額の集合を{xi|i∈Ω}と書くが、Ωは指定された期間に含まれる取引のレコード番号の集合であり、xiはレコード番号iの取引の取引金額である。
【０１０８】
次に、損失率密度読込部７０３により、損失率密度関数P∞(Θ=θ)を読み込む。損失率密度関数は少なくとも０以上１以下の範囲で定義されており、出力は０以上の値を取る。また、０から１までの区間で積分すると１になる。第４の実施形態では図８のようなヒストグラムを考えるが、具体的には次の表５３のような形式でデータを損失率密度記憶部７０４に記憶する。
【０１０９】
【表５３】
【０１１０】
次に、損失事故率読込部７０８により、全取引に対する損失事故率P(L=1)を読み込み、損失事故率記憶部７０９に記憶する。
【０１１１】
次に、巨額取引損失事故率読込部６１０により巨額取引損失事故率（第４の実施形態では、１００万円以上の取引に対する損失事故率とする）Plを読み込み、巨額取引損失事故率記憶部７１１に記憶する。
【０１１２】
次に、巨額損失密度算出部７０５では、取引金額記憶部７０２から取引金額{xi|i∈Ω}、損失率密度記憶部７０４から損失密度関数P∞(Θ=θ)の情報を受け取り、損失事故率記憶部７０９から損失事故率P(L=1)、巨額取引損失事故率記憶部７１１から巨額取引損失事故率Plを受取り、
【数５４】
を計算する。この計算はｙを変化させながら繰り返し実行する。次の表５５は１００万円毎にこの計算を実行した結果をまとめた例である。
【０１１３】
【表５５】
【０１１４】
ｙを１円毎変化させて計算を実行してもかまわないが、その場合は計算量が膨大になるので、第４の実施形態では１００万円ごとに計算を行うことする(例えばy=1.5百万円の場合の数値を１〜２の欄に記載する)。ここではさらに、表５５をもとに次の表５６を作成しておく。
【０１１５】
【表５６】
【０１１６】
表５５によると10,000超の損失発生確率が０であるから、表５６の10,000超の累積確率は１（＝１−０）である。さらに、表５５で、9,999〜10,000の損失発生確率が2.39348E-08であるから、9,999までの累積確率は0.999999976（＝1−2.39348E-08）となる。以下同様の手続きを繰り返して表５６を完成させることが出来る。
【０１１７】
次に、損失事故率記憶部７０９に記憶した損失事故率と、算出された損失累積確率とを用いてリスク量算出部７０６は、VaR等のリスクおよび損失金額の期待値を算出する。
【０１１８】
リスク量算出のための処理手順の一例は次の通りである。
【０１１９】
はじめに、λ=事故発生回数の期待値とし、パラメータλのポアッソン分布を考える。
【数５７】
【０１２０】
λ＝１０の場合のpkの例を次の表５８に示す。
【０１２１】
【表５８】
【０１２２】
表５８の累積確率は、
【数５９】
で表され、例えば、ｋ＝２であれば、事故件数が２件以下である確率を意味する。
【０１２３】
パラメータλにしたがって、予めこのようなテーブルを用意しておく。
【０１２４】
次に、一様乱数発生器により０から１までの実数をランダムに発生する。
【０１２５】
乱数の値をｘとしたときに、上記テーブルをもとに、事故件数を決定する。表５８では、4.53999E-05未満であれば事故は０件、4.53999E-05以上0.000499399未満であれば事故は１件、・・・とする。例えば、ｘ＝０．８であれば、0.791556476以上、0.864464423未満であるから、事故件数は１３件とする。
【０１２６】
次に、上記損失事故のそれぞれについての損失金額をセットする。事故件数回一様乱数を発生させ、表５６を用いて損失額を確定する。第４の実施形態では巨額損失のみを対象に考察しているため、表５６で百万円未満の累積確率が定まっていないが、表５６から損失が百万円未満である確率は９９％であることは読み取れる。したがって、０〜１の間で一様乱数を発生させて値が０．９９以下であれば百万円未満の損失であることがわかるが、このときの金額は０円以上百万円以下の値uをセットする。あらかじめこの値を設定しておく必要がある。ここではu=千円とする。
【０１２７】
具体的には、乱数の値が０．５であれば、0.99未満であるから、損失金額は千円となる。また、乱数の値が、0.9930 であれば、損失額が５百万円と6百万円の間であることがわかるが、ここでは最悪のケースを想定して６百万円とする。１３回事故が発生すれば、それぞれについて同様の処理により損失金額を定める。
【０１２８】
次に、事故回数分の損失金額を合算する。例えば、１３回分の損失金額を合計して、例えば15,568千円を得る。
【０１２９】
以上の処理を指定された回数（N回）繰り返す。リスク量出力部７０７は、それぞれのラウンドで算出される損失金額合計を大きい順に並べて、上位０．１％の値を９９．９％水準のVaRとして出力する。また、損失金額の平均も求めて期待損失額として出力する。例えば、N＝１００万回であるとすると、上位０．１％は上位１０００番目の損失金額に対応している。また、損失額の上位０．１％の計算結果について抜き出し、抜き出したサンプルの損失額の期待値を計算してこれをCVaRとして出力する。
【０１３０】
（第５の実施形態）
図９は本発明の第５の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図である。第５の実施形態は第４の実施形態と同様に、少額損失の分析を行わず処理手順を簡単化したものである。第５の実施形態が第４の実施形態と相違する点は、確率密度を読み込むのではなく、累積確率を読み込む点である。
【０１３１】
つまり、図９において損失率累積確率読込部９０３が第４の実施形態（図６）との違いである。この結果、引き続き行われる巨額損失累積確率算出部９０５による処理も第４の実施形態とは異なる。
【０１３２】
先ず、取引金額読込部９０１により、過去１年間の銀行の各取引における取引金額を読み込む。次の表６０のような形式のデータから、指定された期間（例えば２００３年４月１日から２００４年３月３１日）の取引を抜き出して、取引金額記憶部９０２に蓄積する。
【０１３３】
【表６０】
【０１３４】
取引金額記憶部９０２には主記憶装置が利用されるが、データ量が多い場合には外部記憶装置が用いられる。ここでは、取引金額の集合を{xi|i∈Ω}と書くが、Ωは指定された期間に含まれる取引のレコード番号の集合であり、xiはレコード番号iの取引の取引金額である。
【０１３５】
次に、損失率累積確率読込部９０３により、損失率累積確率すなわち、
【数６１】
を読み込む。損失率累積確率は少なくとも０以上１以下の範囲で定義されており、出力も０以上１以下の値を取る。第５の実施形態では次の表６２のような形式でデータを損失率累積確率関数記憶部９０４に記憶する。
【０１３６】
【表６２】
【０１３７】
次に、損失事故率読込部９０８により、全取引に対する損失事故率P(L=1)を読み込んで損失事故率記憶部９０９に記憶する。
【０１３８】
次に、巨額取引損失事故率読込部９１０により、巨額取引損失事故率（第５の実施形態では、１００万円以上の取引に対する損失事故率とする）Plを読み込み、巨額取引損失事故記憶部９１１に記憶する。
【０１３９】
次に、巨額損失累積確率算出部９０５は、取引金額記憶部９０２から取引金額{xi|i∈Ω}、損失率累積確率関数記憶部９０４から損失率累積確率関数F∞(Θ=θ)の情報を受け取り、損失事故率記憶部９０９から損失事故率P(L=1)、巨額取引損失事故記憶部９１１から巨額取引損失事故率Plを受取り、
【数６３】
を計算する。この計算はｙを変化させながら繰り返し実行する。次の表６４は１００万円毎にこの計算を実行した結果をまとめた例である。
【０１４０】
【表６４】
【０１４１】
ｙを１円毎変化させて計算を実行してもかまわないが、その場合は計算量が膨大になるので、第５の実施形態では１００万円ごとに計算を行うこととし、その間の値は直線で補完するものとする（もちろん、二次曲線等の曲線で近似しても良い）。
【０１４２】
次に、損失事故率記憶部９０９に記憶した損失事故率と、算出された損失累積確率とを用いて、リスク量算出部９０６はVaR等のリスクおよび損失金額の期待値を算出する。
【０１４３】
リスク量算出のための処理手順の一例は以下の通りである。
【０１４４】
はじめに、λ=事故発生回数の期待値とし、パラメータλのポアッソン分布を考える。
【数６５】
【０１４５】
λ＝１０の場合のpkの例を次の表６６に示す。
【０１４６】
【表６６】
【０１４７】
表６６の累積確率は、
【数６７】
であり、例えば、ｋ＝２であれば、事故件数が２件以下である確率を意味する。
【０１４８】
パラメータλにしたがって、予めこのようなテーブルを用意しておく。
【０１４９】
次に、一様乱数発生器により０から１までの実数をランダムに発生する。
【０１５０】
次に、乱数の値をｘとしたときに、上記テーブルをもとに、事故件数を決定する。表６６では、4.53999E-05未満であれば事故は０件、4.53999E-05以上0.000499399未満であれば事故は１件、・・・とする。例えば、ｘ＝０．８であれば、0.791556476以上、0.864464423未満であるから、事故件数は１３件とする。
【０１５１】
上記損失事故のそれぞれについての損失金額をセットする。事故件数回、一様乱数を発生させ、表６４を用いて損失額を確定する。第５の実施形態では巨額損失のみを対象に考察しているため、表６４で百万円未満の累積確率が定まっていないが、表６４から損失が百万円未満である確率は９９％であることは読み取れる。したがって、０〜１の間で一様乱数を発生させて値が０．９９以下であれば百万円未満の損失であることがわかるが、このときの金額は０円以上百万円以下の値uをセットする。あらかじめこの値を設定しておく必要がある。ここではu=千円とする。
【０１５２】
具体的には、乱数の値が０．５であれば、0.99未満であるから、損失金額は千円となる。また、乱数の値が、0.9930 であれば、損失額が５百万円と6百万円の間であることがわかるが、ここでは最悪のケースを想定して６百万円とする。これ以外に、直線で補完する方法も考えられる。その場合は、
【数６８】
により計算できる。１３回事故が発生すれば、それぞれについて同様の処理により損失金額を定める。
【０１５３】
次に、事故回数分の損失金額を合算する。例えば、１３回分の損失金額を合計して、15,568千円を得る。
【０１５４】
上記の処理を指定された回数（N回）繰り返す。リスク量出力部９０７は、それぞれのラウンドで算出される損失金額合計を大きい順に並べて、上位０．１％の値を９９．９％水準のVaRとして出力する。また、損失金額の平均も求めて期待損失額として出力する。例えば、N＝１００万回であるとすると、上位０．１％は上位１０００番目の損失金額に対応している。また、損失額の上位０．１％の計算結果について抜き出し、抜き出したサンプルの損失額の期待値を計算してこれをCVaRとして出力する。
【０１５５】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【図面の簡単な説明】
【０１５６】
【図１】本発明の第１の実施形態に係り、密度関数をモデル化したヒストグラムの一例を示す図
【図２】本発明の第１の実施形態に係り、損失額のヒストグラムの一例を示す図
【図３】本発明の第１の実施形態に係り、損失額のヒストグラムの別の例を示す図
【図４】本発明の第２の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図
【図５】本発明の第２の実施形態に係り、損失率密度関数のヒストグラムを示す図
【図６】本発明の第３の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図
【図７】本発明の第４の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図
【図８】本発明の第４の実施形態に係り、損失率密度関数のヒストグラムを示す図
【図９】本発明の第５の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図
【符号の説明】
【０１５７】
４０１…取引金額読込部
４０２…取引金額記憶部
４０３…損失率密度読込部
４０４…損失率密度記憶部
４０５…巨額損失密度算出部
４０６…リスク量算出部
４０７…リスク量出力部
４０８…損失事故率読込部
４０９…損失事故記憶部
４１０…巨額取引損失事故率読込部
４１１…巨額取引損失事故率記憶部
４１２…少額条件付損失密度読込部
４１３…少額条件付損失密度記憶部
４１４…損失密度算出部
【技術分野】
【０００１】
本発明は、統計学的手法に基づいてオペレーショナルリスクを計量化する装置、方法、およびプログラムに関する。
【背景技術】
【０００２】
新BIS(Bank for International Settlements；国際決済銀行)規制では新しいオペレーショナルリスクが導入される。これまでのBIS規制では信用リスクと市場リスクに対して８％以上の自己資本を準備しておくことが要求されていた。これに対して新BIS規制では、信用リスク、市場リスク、およびオペレーショナルリスクの合計に対して８％の自己資本を準備しておくことが要求される。より正確には、９９．９％水準のVaR（バリューアットリスク）に対応するリスク量に対して８％以上の自己資本を用意しておくことが要求される。
【０００３】
リスク量の評価には、基礎的手法、標準的手法、先進的手法のいずれかを採用してよいことになっているが、前者になるほど精度の低い評価であることから、余裕を持った評価が求められ、リスク量としては大きくならざるを得ないと考えられる。このとき、リスク相応の自己資本を調達しておかなければならない。一般に、自己資本の調達にはコストが掛かるため望ましいことではないし、対外的にはリスクが小さくて自己資本比率が高いということは良い評価に結びつくため、銀行にとって先進的手法を導入しようというインセンティブが存在する。
【０００４】
しかし、先進的手法を採用した場合、オペレーショナルリスクを合理的に計量化する必要が生じるが、過去の損失事故の件数が十分にないために正確な計量化が難しいという問題がある。特に、リスク評価で重要な巨額損失のデータが少ないことが問題である。この課題を解決するために、各種方法が提案されているが、決定的な方法がまだ確定していないのが現状である。オペレーショナルリスクの計量化については、例えば非特許文献１に記載されている。
【非特許文献１】「オペレーショナル・リスクのすべて」，三菱信託銀行オペレーショナル・リスク研究会編，東洋経済新聞社
【発明の開示】
【発明が解決しようとする課題】
【０００５】
本発明は、少ない損失事例から巨額損失のオペレーショナルリスクを合理的に計量化できるオペレーショナルリスクの計量化装置、方法、およびプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の一観点に係るオペレーショナルリスクの計量化装置は、取引金額を入力する手段と、損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を入力する手段と、損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を前記取引金額と前記損失率密度とに基づいて算出する巨額損失密度算出手段と、前記巨額損失密度算出手段により算出された巨額損失密度からオペレーショナルリスクを算出するリスク量算出手段と、を具備する。
【発明の効果】
【０００７】
本発明によれば、少ない損失事例から巨額損失のオペレーショナルリスクを合理的に計量化できるオペレーショナルリスクの計量化装置、方法、およびプログラムを提供できる。
【発明を実施するための最良の形態】
【０００８】
以下、図面を参照しながら本発明の実施形態を説明する。
【０００９】
先ずは、幾つかの仮定の下にオペレーショナルリスク評価のためのモデルを導出する過程を第１の実施形態により説明したのち、同モデルを用いた具体的な定量化に係る実施形態（第２〜第５の実施形態）を説明する。第２及び第３の実施形態では巨額損失のみならず少額損失をも分析可能である。第４及び第５の実施形態では少額損失の分析を行わず処理手順を簡単化したものである。
【００１０】
（第１の実施形態）
次のような確率変数、すなわち取引金額の確率変数をX（ただし、Xはゼロ以上の実数）、損失額の確率変数をY（ただし、Yはゼロ以上の実数）、損失率の確率変数をΘ（ただし、Θゼロ以上１以下の実数）、損失発生の有無：L（０：損失無し、１：損失有り）を考えるとき、損失率Θは次の式(1)のように表すことができる。
Θ=Y/X…(1)
なお、取引金額、損失額は実際にはゼロ以上の整数であるが、説明を簡単にするために実数であるものとして説明する。
【００１１】
始めに、P(Y=y,L=l,X=x)について、以下の関係が成立する。
【００１２】
P(Y=y,L=l,X=x)=P(Y=y|L=l,X=x)P(L=l|X=x)P(X=x)…(2)
また、金額xの取引において損失事故が発生する確率P(L=1|X=x)は、取引金額xの広義単調減少関数であると考えられる。ただし、人間の注意力には限界があるため、ある程度大きな金額になると一定値になると考える。そこで、次の仮定１を置く。
【００１３】
（仮定１）x≧xth1であればP(L=1|X=x)=Pl
このとき、
P(L=1|X=x)=Pl…(3)であり、
P(L=0|X=x)=1-Pl…(4)である。さらに、P(Y=y|L=0,X=x)では、次の仮定２を置く。
【００１４】
（仮定２）損失事故が発生しない場合の損失金額は０である。このとき、
P(Y-y|L=0,X=x)=δ(y)…(5)
である。ここで、δ(y)はディラックのデルタ関数である。また、
【数１】
である。ここで、Θは回収率を表し、P(Θ=θ|L=1,X=x)は、x円の取引で事故が発生したときに回収できる割合がθである確率密度を表している。ここでは、一定以上の金額における回収率は取引金額xに依らない密度関数に従うと考えて、次の仮定３を置く。
【００１５】
（仮定３）x≧xth2であればP(Θ=θ|L=1,X=x)=P∞(Θ=θ)
次に、P(Y=y,L=l,X=x)の周辺密度関数P(Y=y,L=l)において、特にl=1の場合を考える。
【００１６】
【数２】
【００１７】
通常、取引金額xの取引でxより大きな損失が出ることは無いので、次の仮定４を置く。
【００１８】
（仮定４）P(Y=y|L=l,X=x)=0 for {(x,y)|x<y}
特に、y≧xth=max(xth1,xth2)であれば、仮定１と仮定３を考慮すると、式(7)より、
【数３】
である。一方、P(Y=y,L=1)=P(Y=y|L=1)P(L=1)…(9)
であるから、
【数４】
である。大数の法則よりサンプル数が十分多いとき、
【数５】
である。ここで、Ωは全サンプルを表す集合である。さらに、yについて積分することで、損失金額が0からyの区間にある確率は、
【数６】
である。ここで、
【数７】
である。また、
【数８】
であるから、サンプル数が十分多いとき、
【数９】
と算出できる。ここで、δj,kはクロネッカのデルタであり、ｊとｋが等しければ１、そうでなければ０である。また、Plは、Ωth={i|xi≧xth}として、仮定１のもとで、
【数１０】
となる。特に、Pl=P(L=1)であれば、
【数１１】
である。
【００１９】
以上の仮定１，２，４の下で巨額損失の分布は式(11)のように表すことができる。これを利用してオペレーショナルリスクの計量化を次の手順により行うことができる。
【００２０】
先ずP∞(Θ=θ)のモデルを作成する（ステップＳ１）。このモデルは、一定以上の取引金額における損失率の確率密度に相当し、後述する実施形態における損失率密度の読込及び記憶の処理手順に対応する。より具体的に説明すると、損失率密度とは、損失率（０〜１の値）を確率変数としたときの確率密度のことである。例えば、１００件の損失事故のなかで、０．５〜０．６（５０％〜６０％）の損失率の事故が１０件だとすると、０．５〜０．６に事象が発生する確率は１０／１００＝０．１となる。さらに、確率密度は区間の幅で割って、０．１／（０．６−０．５）＝１となる。これが損失率密度である。
【００２１】
P∞(Θ=θ)は、Ωthに属するサンプルにより、θについての密度関数をモデル化すればよい。ノンパラメトリック手法（参考文献：「平滑化とノンパラメトリック回帰への招待」J・S・シモノフ 農林統計協会）が有力な方法であるが、ここでは、最も単純なヒストグラムを使った例を説明する。
【００２２】
P∞(Θ=θ)は、0≦θ≦1の範囲でゼロ以上の値を取り得るが、これ以外の領域では０となる。そこで、モデル化は0≦θ≦1におけるP∞(Θ=θ)の関数を決めることである。
【００２３】
以下のようにGk(k=1,2,…,K)を定義する。
【数１２】
である。また、
【数１３】
である。このとき、
【数１４】
として、P∞(Θ=θ)=Pk（θ∈Rkの場合）となる。この場合、P∞(Θ=θ)=Pkの例を図１に示す。この例では、閾値すなわち、k/K(k=0,…,K)の間隔は一定であったが、閾値を任意の値に変化させても良い。最適な閾値を算出する方法については、例えば「情報量統計学」坂本慶行他、共立出版株式会社等に記載されている。また、前出の「平滑化とノンパラメトリック回帰への招待」には、ヒストグラムに代えて頻度ポリゴンを使用したり、カーネル関数を使う方法が記載されている。
【００２４】
次に、
【数１５】
を算出する（ステップＳ２）。ここでは、全取引サンプルの数をNとし、このうちで、事故の発生した数をN1とし、P(L=1)=N1/Nとして算出する。算出される値は、後述する実施形態における、損失事故率の読込及び記憶に対応する。
【００２５】
次に、
【数１６】
を算出する（ステップＳ３）。
【００２６】
xthを外部からうけとり、サンプルの中で取引金額がxth以上のサンプルの数を数えてnとする。さらに、このうちで事故が発生した数をn1とし、Pl=n1/nとして算出する。算出される値は、後述する実施形態における、巨額取引損失事故率の読込及び記憶に対応する。
【００２７】
次に、
【数１７】
を計算する（ステップＳ４）。
【００２８】
ただし、ymaxは取引金額の最大値とし、ymin=xthとする。
【００２９】
ymin≦y≦ymaxの範囲で複数の代表点を決めておく。第１の実施形態では外部から与える数値Mを用いて、
【数１８】
としておく。そして全てのｍについて、以下のステップＳ４−１〜Ｓ４−３の処理を繰り返す。
【００３０】
まず、y=ymをセットする（ステップＳ４−１）。
【００３１】
次に、
【数１９】
を算出する。そして、全ての取引サンプルiについて、θi=y/xiを算出し、上記ステップＳ１で求めた関数P∞(Θ=θ)を使って、
【数２０】
を算出する。算出された値の全てのサンプルについての平均を求める（ステップＳ４−２）。
【００３２】
次に、上記ステップＳ２，Ｓ３，およびＳ４−２の結果を用い、
【数２１】
に従って、pmmaxを算出する（ステップＳ４−３）。
【００３３】
次に、P(Y=y|L=1)をヒストグラムで表す。
【数２２】
では、P(Y=y|L=1)=pmmaxとする。これにより、図２のようなヒストグラムが得られる。ここで、η=(y0+y1)/2である（ステップＳ５）。得られるヒストグラムは、後述する実施形態において算出される巨額損失密度に対応する。ここで、損失額を確率変数としたときの確率密度を「損失密度」と呼ぶ。特に、損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する損失密度のことを「巨額損失密度」と呼ぶ。
【００３４】
次に、ｙが小さいときのモデルを作成する。このモデルは、後述する実施形態における少額条件付損失密度の読込及び記憶に対応する。
【００３５】
上記ステップＳ１と同様に、ヒストグラムを使ってｙがyminよりも小さい範囲のヒストグラムを描く（図３）。
【００３６】
全サンプルから、事故が発生しており、かつ、損失額yがηより小さいサンプルを集めて、Ωminとする。これをU個に分ける以下のような集合を考える。
【数２３】
【００３７】
このとき、事故が発生しており、かつ、損失額yがηより小さいサンプルの数をNminとし、Huに属するサンプルの数をNuminとして、
【数２４】
を算出する（ステップＳ６）。
【００３８】
次に、P(Y=y|L=1)の分布を求める（ステップＳ７）。これは、後述する実施形態における損失密度の算出に対応する。
【００３９】
上記ステップＳ６で求めたpuminは事故が発生しており、かつ、損失額yがηより小さいという条件の下での分布関数であるため、P(Y=y|L=1)を算出するには、事故が発生しており損失額yがηより小さい確率qを乗じる必要がある。これは以下のように求める。
【００４０】
すなわち、
【数２５】
であることから、
【数２６】
とする。
【００４１】
最後に、オペレーショナルリスクの計量化を行う（ステップＳ８）。計量化の手法については別の実施形態において説明するが、公知の手法を適用可能である。
【００４２】
以上説明した第１の実施形態によれば、上述した４つの仮定の下にオペレーショナルリスク評価のためのモデルを導出することができ、同モデルを用い、ステップＳ８においては合理的にオペレーショナルリスクの計量化を行うことが可能になる。尚、ステップＳ８におけるオペレーショナルリスクの計量化には、公知の手法を適用することができる。
【００４３】
（第２の実施形態）
本発明の第２の実施形態では巨額損失のみならず少額損失をも分析可能である。
【００４４】
図４は本発明の第２の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図である。本装置はコンピュータを利用して実現できる。また、本発明に係る機能を実現するためのプログラムが図示しないプログラム記憶装置に格納される。プログラム記憶装置は、例えば、磁気ディスク装置や光ディスク装置等の外部記憶装置からなる。上記プログラムがＣＰＵからの制御でランダムアクセスメモリ（ＲＡＭ）等からなる主記憶装置に読み込まれ、同ＣＰＵにより実行されることにより、コンピュータを本発明に係るオペレーショナルリスク定量化装置として機能させることができる。なお、オペレーショナルリスク定量化装置には、各種コンピュータ資源を管理し、ファイルシステムやネットワーク通信機能、ならびにグラフィカルユーザインタフェース（ＧＵＩ）等を提供するオペレーティングシステムも導入されている。なお、後述する第３乃至第５の実施形態に係るオペレーショナルリスク定量化装置についても、コンピュータを利用して実現できる。
【００４５】
先ず、取引金額読込部４０１により、過去１年間の銀行の各取引における取引金額を読み込む。次の表２７のような形式のデータから、指定された期間（例えば2003年４月１日から2004年３月３１日）の取引を抜き出して、取引金額記憶部４０２に蓄積する。
【００４６】
【表２７】
【００４７】
取引金額記憶部４０２には主記憶装置が利用されるが、データ量が多い場合には外部記憶装置が用いられる。ここでは、取引金額の集合を{xi|i∈Ω}と書くが、Ωは指定された期間に含まれる取引のレコード番号の集合であり、xiはレコード番号iの取引の取引金額である。
【００４８】
次に、損失率密度読込部４０３により損失率密度関数P∞(Θ=θ)を読み込む。損失率密度関数は少なくとも０以上１以下の範囲で定義されており、出力は０以上の値を取る。また、０から１までの区間で積分すると１になる。本実施形態では図５のような損失率密度関数のヒストグラムを考えるが、具体的には次の表２８のような形式でデータを損失率密度記憶部４０４に記憶する。
【００４９】
【表２８】
【００５０】
以下の実施形態の説明において、X〜Yと書いたときは、X以上Y未満を表すものとする。ただし、Y=1の場合だけ、X以上Y以下を意味するものとする。０．２１という値がくれば、この表２８の０．２以上０．３未満の欄をみて、０．３という値を返す。同様に考えて、値０．２が入力されると、０．３を出力し、値０．３の場合は０．２を出力する。また、特に１．０の場合は３．５であるとする。
【００５１】
次に、損失事故率読込部４０８により、全取引に対する損失事故率P(L=1)を読み込で損失事故率記憶部４０９に記憶する。なお、乱数に基づく事故発生モデルを含むオペレーション品質評価モデルから損失事故率を算出してもよい。
【００５２】
次に、巨額取引損失事故率読込部４１０により巨額取引損失事故率（本実施形態では、１００万円以上の取引に対する損失事故率とする）Plを読み込んで巨額取引損失事故率記憶部４１１に記憶する。
【００５３】
次に、少額条件付損失密度読込部４１２は、１００万円未満の損失が発生したという条件付の損失密度を読み込んで、少額条件付損失密度記憶部４１３に記憶する。損失密度は次の表２９のような形で記録される。
【００５４】
【表２９】
【００５５】
損失分布は０円から999,999円までの間で和をとると１となる（損失額０円の確率＋損失額１円の確率＋・・・損失額999,999の確率＝１）。
【００５６】
次に、巨額損失密度算出部４０５は、取引金額記憶部４０２から取引金額{xi|i∈Ω}、損失率密度記憶部４０４から損失密度関数P∞(Θ=θ)の情報を受け取り、損失事故率記憶部４０９から損失事故率P(L=1)、巨額取引損失事故率記憶部４１１から巨額取引損失事故率Plを受取り、
【数３０】
を計算する。この計算はｙを変化させながら繰り返し実行する。次の表３１は１００万円毎にこの計算を実行した結果をまとめた例である。
【００５７】
【表３１】
【００５８】
ｙを１円毎変化させて計算を実行してもかまわないが、その場合は計算量が膨大になるので、本実施形態では１００万円ごとに計算を行うことする(例えばy=1.5百万円の場合の数値を１〜２の欄に記載する)。
【００５９】
次に、損失密度算出部４１４は、巨額損失密度算出部４０５で算出した巨額損失密度と少額条件付損失密度記憶部４１３に記憶された少額条件付損失密度から損失密度関数を算出する。具体的には、１００万円以上の損失密度は巨額損失密度算出部４０５で算出した巨額損失密度とする。巨額損失密度算出部４０５で算出した巨額損失密度に基づき、１００万円以上の損失が発生する確率は以下の式に従って算出できる。
【数３２】
【００６０】
１００万円未満の損失密度については、少額条件付損失密度記憶部４１３に記憶された少額条件付損失密度に1-ηを乗じた値とする。η=0.01であったとすると、次の表３３のような１００万円未満の損失密度を得る。
【００６１】
【表３３】
【００６２】
これを、１００万円以上の巨額損失密度とあわせると、次の表３４を得る。
【００６３】
【表３４】
【００６４】
さらに、損失累積確率として次の表３５を得る。
【００６５】
【表３５】
【００６６】
次に、損失事故率記憶部４０９に記憶した損失事故率と損失密度算出部４１４で算出した損失累積確率とを用いてリスク量算出部４０６は、VaR等のリスクおよび損失金額の期待値をオペレーショナルリスクとして算出する。VaRについては、例えば参考文献（「金融リスクの計量化 上 バリューアットリスク」，木島正明編著）に詳しい。ここで、リスク量算出部４０６により算出されるオペレーショナルリスクは、将来の取引に対するリスク量のことをいう。損失額が大方の予想に反して非常に大きくなることを「リスク」と呼ぶ。オペレーショナルリスクは、所定期間の損失額合計を横軸とし、個々の損失額合計値に対する確率密度を縦軸とする損失額合計密度関数において、所定の上側のパーセンタイル値として求めることができる。あるいは、所定期間の損失額合計を横軸とし、個々の損失額合計値に対する確率密度を縦軸とする損失額合計密度関数において、所定の上側のパーセンタイル値以上の損失額合計値の平均値として求めてもよい。
【００６７】
より具体的に、リスク量算出のための処理手順の一例を説明する。はじめに、λ=P(L=1)とし、パラメータλのポアッソン分布を考える。
【数３６】
【００６８】
λ＝１０の場合のpkの例を次の表３７に示す。
【００６９】
【表３７】
【００７０】
表３７の累積確率は、
【数３８】
であり、たとえば、ｋ＝２であれば、事故件数が２件以下である確率を意味する。
【００７１】
パラメータλにしたがって、予めこのようなテーブルを用意しておく。
次に、一様乱数発生器により０から１までの実数をランダムに発生する。
次に、乱数の値をｘとしたときに、上記テーブルをもとに、事故件数を決定する。表３７では、4.53999E-05未満であれば事故は０件、4.53999E-05以上0.000499399未満であれば事故は１件、・・・とする。たとえば、ｘ＝０．８であれば、0.791556476以上、0.864464423未満であるから、事故件数は１３件とする。
【００７２】
次に、上記損失事故のそれぞれについての損失金額をセットする。事故件数回一様乱数を発生させ、表３５を用いて損失額を確定する。たとえば、乱数の値が０．５であれば、0.4851以上0.5544未満であるから、損失金額は７円となる。１３回事故が発生すれば、それぞれについて同様の処理により損失金額を定める。
【００７３】
次に、事故回数分の損失金額を合算する。たとえば、１３回分の損失金額を合計して、例えば15,568を得る。
【００７４】
以上の処理を指定された回数（N回）繰り返す。リスク量出力部４０７は、それぞれのラウンドで算出される損失金額合計を大きい順に並べて、上位０．１％の値を９９．９％水準のVaRとして出力する。また、損失金額の平均も求めて期待損失額として出力する。たとえば、N＝１００万回であるとすると、上位０．１％は上位１０００番目の損失金額に対応している。また、損失額の上位０．１％の計算結果について抜き出し、抜き出したサンプルの損失額の期待値を計算してこれをCVaR（コンディショナルバリューアットリスク）として出力する。
【００７５】
（第３の実施形態）
図６は本発明の第３の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図である。第３の実施形態が第２の実施形態と相違する点は、損失の確率密度を読み込むのではなく、損失の累積確率を読み込む点である。つまり、図６における損失率累積確率読込部６０３および少額条件付損失累積確率読込部６１２が第２の実施形態（図４）との違いである。この結果、引き続き行われる巨額損失累積確率算出部６０５、損失累積確率算出部６１４の処理手順も第２の実施形態とは異なる。
【００７６】
先ず、取引金額読込部６０１により、過去１年間の銀行の各取引における取引金額を読み込む。次の表３９のような形式のデータから、指定された期間（例えば2003年４月１日から2004年３月３１日）の取引を抜き出して、取引金額記憶部６０２に蓄積する。
【００７７】
【表３９】
【００７８】
取引金額記憶部６０２には主記憶装置が利用されるが、データ量が多い場合には外部記憶装置が用いられる。ここでは、取引金額の集合を{xi|i∈Ω}と書くが、Ωは指定された期間に含まれる取引のレコード番号の集合であり、xiはレコード番号iの取引の取引金額である。
【００７９】
次に、損失率累積確率読込部６０３により、損失率累積確率すなわち、
【数４０】
を読み込む。損失率累積確率は少なくとも０以上１以下の範囲で定義されており、出力も０以上１以下の値を取る。本実施形態では次の表４１のような形式でデータを損失率累積確率関数記憶部６０４に記憶する。
【００８０】
【表４１】
【００８１】
次に、損失事故率読込部６０８により、全取引に対する損失事故率P(L=1)を読み込んで損失事故率記憶部６０９に記憶する。
【００８２】
次に、巨額取引損失事故率読込部６１０により巨額取引損失事故率（本実施形態では、１００万円以上の取引に対する損失事故率とする）Plを読み込んで巨額取引損失事故記憶部６１１に記憶する。
【００８３】
次に、少額条件付損失累積確率読込部６１２は、１００万円未満の損失密度を読み込み、少額条件付損失累積確率記憶部６１３に記憶する。損失累積確率は次の表４２のような形で記録される。
【００８４】
【表４２】
【００８５】
次に、巨額損失累積確率算出部６０５は、取引金額記憶部６０２から取引金額{xi|i∈Ω}、損失率累積確率関数記憶部６０４から損失率累積確率関数F∞(Θ=θ)の情報を受け取り、損失事故率記憶部６０９から損失事故率P(L=1)、巨額取引損失事故記憶部６１１から巨額取引損失事故率Plを受取り、
【数４３】
を計算する。この計算はｙを変化させながら繰り返し実行する。次の表４４は１００万円毎にこの計算を実行した結果をまとめた例である。
【００８６】
【表４４】
【００８７】
ｙを１円毎変化させて計算を実行してもかまわないが、その場合は計算量が膨大になるので、本実施形態では１００万円ごとに計算を行うこととし、その間の値は直線で補完するものとする（もちろん、二次曲線等の曲線で近似しても良い）。
【００８８】
次に、損失累積確率算出部６１４では、巨額損失累積確率算出部６０５で算出した巨額損失累積確率と少額条件付損失累積確率記憶部６１３に記憶された少額条件付損失累積確率から損失累積確率を算出する。具体的は、１００万円以上の損失累積確率は巨額損失累積確率算出部６０５で算出した巨額損失累積確率とする。巨額損失累積確率算出部６０５で算出した巨額損失累積確率をつかい、１００万円未満の損失が発生する確率は以下の式に従って算出できる。
【数４５】
【００８９】
１００万円未満の損失累積確率については、少額条件付損失累積確率記憶部６１３に記憶された少額条件付損失密度に1-ηを乗じた値とする。
【００９０】
η=0.01であったとすると、次の表４６のような１００万円未満の損失累積確率を得る。
【００９１】
【表４６】
【００９２】
１００万円以上とあわせて、次の表４７を得る。
【００９３】
【表４７】
【００９４】
次に、損失事故率記憶部６０９に記憶した損失事故率と損失累積確率算出部６１４で算出した損失累積確率を用いて、リスク量算出部６０６はVaR等のリスクおよび損失金額の期待値を算出する。リスク量算出のための処理手順の一例は次の通りである。
【００９５】
はじめに、λ=P(L=1)とし、パラメータλのポアッソン分布を考える。
【数４８】
【００９６】
λ＝１０の場合のpkの例を次の表４９に示す。
【００９７】
【表４９】
【００９８】
表４９の累積確率は、
【数５０】
であり、たとえば、ｋ＝２であれば、事故件数が２件以下である確率を意味する。
【００９９】
パラメータλにしたがって、予めこのようなテーブルを用意しておく。
【０１００】
次に、一様乱数発生器により０から１までの実数をランダムに発生する。
【０１０１】
次に、乱数の値をｘとしたときに、上記テーブルをもとに、事故件数を決定する。表４９では、4.53999E-05未満であれば事故は０件、4.53999E-05以上0.000499399未満であれば事故は１件、・・・とする。たとえば、ｘ＝０．８であれば、0.791556476以上、0.864464423未満であるから、事故件数は１３件とする。
【０１０２】
次に、上記損失事故のそれぞれについての損失金額をセットする。事故件数回一様乱数を発生させ、表４７を用いて損失額を確定する。たとえば、乱数の値が０．５であれば、0.693以下であるから、損失金額は０円から１０円の間にあることが分かる。（０，０）→（10,0.693）の間を直線で補間すると、
【数５１】
となり、損失額は７円とする。１３回事故が発生すれば、それぞれについて同様の処理により損失金額を定める。
【０１０３】
次に、事故回数分の損失金額を合算する。たとえば、１３回分の損失金額を合計して、15,568を得る。
【０１０４】
以上の処理を指定された回数（N回）繰り返す。リスク量出力部６０７は、それぞれのラウンドで算出される損失金額合計を大きい順に並べて、上位０．１％の値を９９．９％水準のVaRとして出力する。また、損失金額の平均も求めて期待損失額として出力する。たとえば、N＝１００万回であるとすると、上位０．１％は上位１０００番目の損失金額に対応している。また、損失額の上位０．１％の計算結果について抜き出し、抜き出したサンプルの損失額の期待値を計算してこれをCVaRとして出力する。
【０１０５】
（第４の実施形態）
第４の実施形態は少額損失の分析を行わず処理手順を簡単化したものである。図７は本発明の第４の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図である。取引金額読込部７０１により、過去１年間の銀行の各取引における取引金額を読み込む。次の表５２のような形式のデータから、指定された期間（例えば２００３年４月１日から２００４年３月３１日）の取引を抜き出して、取引金額記憶部７０２に蓄積する。
【０１０６】
【表５２】
【０１０７】
取引金額記憶部７０２には主記憶装置が利用されるが、データ量が多い場合には外部記憶装置が用いられる。取引金額の集合を{xi|i∈Ω}と書くが、Ωは指定された期間に含まれる取引のレコード番号の集合であり、xiはレコード番号iの取引の取引金額である。
【０１０８】
次に、損失率密度読込部７０３により、損失率密度関数P∞(Θ=θ)を読み込む。損失率密度関数は少なくとも０以上１以下の範囲で定義されており、出力は０以上の値を取る。また、０から１までの区間で積分すると１になる。第４の実施形態では図８のようなヒストグラムを考えるが、具体的には次の表５３のような形式でデータを損失率密度記憶部７０４に記憶する。
【０１０９】
【表５３】
【０１１０】
次に、損失事故率読込部７０８により、全取引に対する損失事故率P(L=1)を読み込み、損失事故率記憶部７０９に記憶する。
【０１１１】
次に、巨額取引損失事故率読込部６１０により巨額取引損失事故率（第４の実施形態では、１００万円以上の取引に対する損失事故率とする）Plを読み込み、巨額取引損失事故率記憶部７１１に記憶する。
【０１１２】
次に、巨額損失密度算出部７０５では、取引金額記憶部７０２から取引金額{xi|i∈Ω}、損失率密度記憶部７０４から損失密度関数P∞(Θ=θ)の情報を受け取り、損失事故率記憶部７０９から損失事故率P(L=1)、巨額取引損失事故率記憶部７１１から巨額取引損失事故率Plを受取り、
【数５４】
を計算する。この計算はｙを変化させながら繰り返し実行する。次の表５５は１００万円毎にこの計算を実行した結果をまとめた例である。
【０１１３】
【表５５】
【０１１４】
ｙを１円毎変化させて計算を実行してもかまわないが、その場合は計算量が膨大になるので、第４の実施形態では１００万円ごとに計算を行うことする(例えばy=1.5百万円の場合の数値を１〜２の欄に記載する)。ここではさらに、表５５をもとに次の表５６を作成しておく。
【０１１５】
【表５６】
【０１１６】
表５５によると10,000超の損失発生確率が０であるから、表５６の10,000超の累積確率は１（＝１−０）である。さらに、表５５で、9,999〜10,000の損失発生確率が2.39348E-08であるから、9,999までの累積確率は0.999999976（＝1−2.39348E-08）となる。以下同様の手続きを繰り返して表５６を完成させることが出来る。
【０１１７】
次に、損失事故率記憶部７０９に記憶した損失事故率と、算出された損失累積確率とを用いてリスク量算出部７０６は、VaR等のリスクおよび損失金額の期待値を算出する。
【０１１８】
リスク量算出のための処理手順の一例は次の通りである。
【０１１９】
はじめに、λ=事故発生回数の期待値とし、パラメータλのポアッソン分布を考える。
【数５７】
【０１２０】
λ＝１０の場合のpkの例を次の表５８に示す。
【０１２１】
【表５８】
【０１２２】
表５８の累積確率は、
【数５９】
で表され、例えば、ｋ＝２であれば、事故件数が２件以下である確率を意味する。
【０１２３】
パラメータλにしたがって、予めこのようなテーブルを用意しておく。
【０１２４】
次に、一様乱数発生器により０から１までの実数をランダムに発生する。
【０１２５】
乱数の値をｘとしたときに、上記テーブルをもとに、事故件数を決定する。表５８では、4.53999E-05未満であれば事故は０件、4.53999E-05以上0.000499399未満であれば事故は１件、・・・とする。例えば、ｘ＝０．８であれば、0.791556476以上、0.864464423未満であるから、事故件数は１３件とする。
【０１２６】
次に、上記損失事故のそれぞれについての損失金額をセットする。事故件数回一様乱数を発生させ、表５６を用いて損失額を確定する。第４の実施形態では巨額損失のみを対象に考察しているため、表５６で百万円未満の累積確率が定まっていないが、表５６から損失が百万円未満である確率は９９％であることは読み取れる。したがって、０〜１の間で一様乱数を発生させて値が０．９９以下であれば百万円未満の損失であることがわかるが、このときの金額は０円以上百万円以下の値uをセットする。あらかじめこの値を設定しておく必要がある。ここではu=千円とする。
【０１２７】
具体的には、乱数の値が０．５であれば、0.99未満であるから、損失金額は千円となる。また、乱数の値が、0.9930 であれば、損失額が５百万円と6百万円の間であることがわかるが、ここでは最悪のケースを想定して６百万円とする。１３回事故が発生すれば、それぞれについて同様の処理により損失金額を定める。
【０１２８】
次に、事故回数分の損失金額を合算する。例えば、１３回分の損失金額を合計して、例えば15,568千円を得る。
【０１２９】
以上の処理を指定された回数（N回）繰り返す。リスク量出力部７０７は、それぞれのラウンドで算出される損失金額合計を大きい順に並べて、上位０．１％の値を９９．９％水準のVaRとして出力する。また、損失金額の平均も求めて期待損失額として出力する。例えば、N＝１００万回であるとすると、上位０．１％は上位１０００番目の損失金額に対応している。また、損失額の上位０．１％の計算結果について抜き出し、抜き出したサンプルの損失額の期待値を計算してこれをCVaRとして出力する。
【０１３０】
（第５の実施形態）
図９は本発明の第５の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図である。第５の実施形態は第４の実施形態と同様に、少額損失の分析を行わず処理手順を簡単化したものである。第５の実施形態が第４の実施形態と相違する点は、確率密度を読み込むのではなく、累積確率を読み込む点である。
【０１３１】
つまり、図９において損失率累積確率読込部９０３が第４の実施形態（図６）との違いである。この結果、引き続き行われる巨額損失累積確率算出部９０５による処理も第４の実施形態とは異なる。
【０１３２】
先ず、取引金額読込部９０１により、過去１年間の銀行の各取引における取引金額を読み込む。次の表６０のような形式のデータから、指定された期間（例えば２００３年４月１日から２００４年３月３１日）の取引を抜き出して、取引金額記憶部９０２に蓄積する。
【０１３３】
【表６０】
【０１３４】
取引金額記憶部９０２には主記憶装置が利用されるが、データ量が多い場合には外部記憶装置が用いられる。ここでは、取引金額の集合を{xi|i∈Ω}と書くが、Ωは指定された期間に含まれる取引のレコード番号の集合であり、xiはレコード番号iの取引の取引金額である。
【０１３５】
次に、損失率累積確率読込部９０３により、損失率累積確率すなわち、
【数６１】
を読み込む。損失率累積確率は少なくとも０以上１以下の範囲で定義されており、出力も０以上１以下の値を取る。第５の実施形態では次の表６２のような形式でデータを損失率累積確率関数記憶部９０４に記憶する。
【０１３６】
【表６２】
【０１３７】
次に、損失事故率読込部９０８により、全取引に対する損失事故率P(L=1)を読み込んで損失事故率記憶部９０９に記憶する。
【０１３８】
次に、巨額取引損失事故率読込部９１０により、巨額取引損失事故率（第５の実施形態では、１００万円以上の取引に対する損失事故率とする）Plを読み込み、巨額取引損失事故記憶部９１１に記憶する。
【０１３９】
次に、巨額損失累積確率算出部９０５は、取引金額記憶部９０２から取引金額{xi|i∈Ω}、損失率累積確率関数記憶部９０４から損失率累積確率関数F∞(Θ=θ)の情報を受け取り、損失事故率記憶部９０９から損失事故率P(L=1)、巨額取引損失事故記憶部９１１から巨額取引損失事故率Plを受取り、
【数６３】
を計算する。この計算はｙを変化させながら繰り返し実行する。次の表６４は１００万円毎にこの計算を実行した結果をまとめた例である。
【０１４０】
【表６４】
【０１４１】
ｙを１円毎変化させて計算を実行してもかまわないが、その場合は計算量が膨大になるので、第５の実施形態では１００万円ごとに計算を行うこととし、その間の値は直線で補完するものとする（もちろん、二次曲線等の曲線で近似しても良い）。
【０１４２】
次に、損失事故率記憶部９０９に記憶した損失事故率と、算出された損失累積確率とを用いて、リスク量算出部９０６はVaR等のリスクおよび損失金額の期待値を算出する。
【０１４３】
リスク量算出のための処理手順の一例は以下の通りである。
【０１４４】
はじめに、λ=事故発生回数の期待値とし、パラメータλのポアッソン分布を考える。
【数６５】
【０１４５】
λ＝１０の場合のpkの例を次の表６６に示す。
【０１４６】
【表６６】
【０１４７】
表６６の累積確率は、
【数６７】
であり、例えば、ｋ＝２であれば、事故件数が２件以下である確率を意味する。
【０１４８】
パラメータλにしたがって、予めこのようなテーブルを用意しておく。
【０１４９】
次に、一様乱数発生器により０から１までの実数をランダムに発生する。
【０１５０】
次に、乱数の値をｘとしたときに、上記テーブルをもとに、事故件数を決定する。表６６では、4.53999E-05未満であれば事故は０件、4.53999E-05以上0.000499399未満であれば事故は１件、・・・とする。例えば、ｘ＝０．８であれば、0.791556476以上、0.864464423未満であるから、事故件数は１３件とする。
【０１５１】
上記損失事故のそれぞれについての損失金額をセットする。事故件数回、一様乱数を発生させ、表６４を用いて損失額を確定する。第５の実施形態では巨額損失のみを対象に考察しているため、表６４で百万円未満の累積確率が定まっていないが、表６４から損失が百万円未満である確率は９９％であることは読み取れる。したがって、０〜１の間で一様乱数を発生させて値が０．９９以下であれば百万円未満の損失であることがわかるが、このときの金額は０円以上百万円以下の値uをセットする。あらかじめこの値を設定しておく必要がある。ここではu=千円とする。
【０１５２】
具体的には、乱数の値が０．５であれば、0.99未満であるから、損失金額は千円となる。また、乱数の値が、0.9930 であれば、損失額が５百万円と6百万円の間であることがわかるが、ここでは最悪のケースを想定して６百万円とする。これ以外に、直線で補完する方法も考えられる。その場合は、
【数６８】
により計算できる。１３回事故が発生すれば、それぞれについて同様の処理により損失金額を定める。
【０１５３】
次に、事故回数分の損失金額を合算する。例えば、１３回分の損失金額を合計して、15,568千円を得る。
【０１５４】
上記の処理を指定された回数（N回）繰り返す。リスク量出力部９０７は、それぞれのラウンドで算出される損失金額合計を大きい順に並べて、上位０．１％の値を９９．９％水準のVaRとして出力する。また、損失金額の平均も求めて期待損失額として出力する。例えば、N＝１００万回であるとすると、上位０．１％は上位１０００番目の損失金額に対応している。また、損失額の上位０．１％の計算結果について抜き出し、抜き出したサンプルの損失額の期待値を計算してこれをCVaRとして出力する。
【０１５５】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【図面の簡単な説明】
【０１５６】
【図１】本発明の第１の実施形態に係り、密度関数をモデル化したヒストグラムの一例を示す図
【図２】本発明の第１の実施形態に係り、損失額のヒストグラムの一例を示す図
【図３】本発明の第１の実施形態に係り、損失額のヒストグラムの別の例を示す図
【図４】本発明の第２の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図
【図５】本発明の第２の実施形態に係り、損失率密度関数のヒストグラムを示す図
【図６】本発明の第３の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図
【図７】本発明の第４の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図
【図８】本発明の第４の実施形態に係り、損失率密度関数のヒストグラムを示す図
【図９】本発明の第５の実施形態に係るオペレーショナルリスク定量化装置を示すブロック図
【符号の説明】
【０１５７】
４０１…取引金額読込部
４０２…取引金額記憶部
４０３…損失率密度読込部
４０４…損失率密度記憶部
４０５…巨額損失密度算出部
４０６…リスク量算出部
４０７…リスク量出力部
４０８…損失事故率読込部
４０９…損失事故記憶部
４１０…巨額取引損失事故率読込部
４１１…巨額取引損失事故率記憶部
４１２…少額条件付損失密度読込部
４１３…少額条件付損失密度記憶部
４１４…損失密度算出部
【特許請求の範囲】
【請求項１】
取引金額を入力する手段と、
損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を入力する手段と、
損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を前記取引金額と前記損失率密度とに基づいて算出する巨額損失密度算出手段と、
前記巨額損失密度算出手段により算出された巨額損失密度からオペレーショナルリスクを算出するリスク量算出手段と、を具備するオペレーショナルリスクの計量化装置。
【請求項２】
前記オペレーショナルリスクは、バリューアットリスク及びコンディショナルバリューアットリスクのいずれかを含む請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項３】
前記巨額損失密度に基づいて、損失額の累積確率を求める手段と、
取引金額によらず、損失事故が発生する確率に相当する損失事故率を入力する手段と、
乱数と、前記損失事故率とに基づいて、複数件の損失事故を発生する手段と、
発生した損失事故の各々の損失額を前記損失額の累積確率に基づいて算出する手段と
を具備し、
算出された損失額を値の大きい順に並べ、所定の基準で選択された値を前記オペレーショナルリスクとする請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項４】
損失額が前記閾値未満であるという少額条件の下に、該損失額を確率変数としたときの条件付き確率密度に相当する少額損失密度を入力する手段と、
前記巨額損失密度に基づく補正量を乗じることにより前記少額損失密度を補正する補正手段と、
前記巨額損失密度と、前記少額損失密度とを組み合わせて前記損失密度を求める手段と、
を具備し、
前記損失額の累積確率を前記損失密度から求める請求項１乃至３のいずれかに記載のオペレーショナルリスクの計量化装置。
【請求項５】
損失額が前記閾値未満である確率を前記巨額損失密度を用いて算出し、該確率を前記補正量とする請求項４に記載のオペレーショナルリスクの計量化装置。
【請求項６】
取引金額によらず、損失の事故が発生する確率に相当する損失事故率と、
巨額損失の事故が発生する条件付き確率に相当する巨額損失事故率との比を前記巨額損失密度算出手段により算出された巨額損失密度に乗じる請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項７】
乱数に基づく事故発生モデルを含むオペレーション品質評価モデルから、前記損失事故率を算出する手段を具備する請求項６に記載のオペレーショナルリスクの計量化装置。
【請求項８】
前記巨額損失密度算出手段は、
一の損失額を一の取引金額で除すことにより損失率を算出し、該損失率に対応する損失率密度を前記一の取引金額で除してなる値を求め、
前記値を複数の取引金額の全てについて求めるとともに、該複数の取引金額の全てについて求めた値の平均値を前記巨額損失密度として算出する請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項９】
取引金額を入力する手段と、
入力された取引金額を記憶する取引金額記憶手段と、
損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を入力する手段と、
入力された損失率密度を記憶する損失率密度記憶手段と、
損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を、前記取引金額記憶手段に記憶された取引金額と、前記損失率密度記憶手段に記憶された損失率密度とに基づいて算出する巨額損失密度算出手段と、
算出された巨額損失密度を記憶する巨額損失密度記憶手段と、
前記巨額損失密度記憶手段に記憶された巨額損失密度からオペレーショナルリスクを算出するリスク量算出手段と、
算出されたリスク量を出力する出力手段と、
を具備するオペレーショナルリスクの計量化装置。
【請求項１０】
所定期間の損失額合計を横軸とし、個々の損失額合計値に対する確率密度を縦軸とする損失額合計密度関数において、所定の上側のパーセンタイル値を前記オペレーショナルリスクの値とする請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項１１】
所定期間の損失額合計を横軸とし、個々の損失額合計値に対する確率密度を縦軸とする損失額合計密度関数において、所定の上側のパーセンタイル値以上の損失額合計値の平均値を前記オペレーショナルリスクの値とする請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項１２】
取引金額入力手段により取引金額を入力するステップと、
損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を損失率密度入力手段により入力するステップと、
損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を、巨額損失密度算出手段が前記取引金額と前記損失率密度とに基づいて算出する巨額損失密度算出ステップと、
前記巨額損失密度算出ステップにより算出された巨額損失密度から、損失額に対するリスク量をリスク量算出手段により算出するリスク量算出ステップと、を具備するオペレーショナルリスクの計量化方法。
【請求項１３】
前記オペレーショナルリスクは、バリューアットリスク及びコンディショナルバリューアットリスクのいずれかを含む請求項１２に記載のオペレーショナルリスクの計量化方法。
【請求項１４】
前記巨額損失密度に基づいて、損失額の累積確率を求めるステップと、
取引金額によらず、損失事故が発生する確率に相当する損失事故率を入力するステップと、
乱数と、前記損失事故率とに基づいて、複数件の損失事故を発生するステップと、
発生した損失事故の各々の損失額を前記損失額の累積確率に基づいて算出するステップと
を具備し、
算出された損失額を値の大きい順に並べ、所定の基準で選択された値を前記オペレーショナルリスクとする請求項１２に記載のオペレーショナルリスクの計量化方法。
【請求項１５】
損失額が前記閾値未満であるという少額条件の下に、該損失額を確率変数としたときの条件付き確率密度に相当する少額損失密度を入力するステップと、
前記巨額損失密度に基づく補正量を乗じることにより前記少額損失密度を補正する補正ステップと、
前記巨額損失密度と、前記少額損失密度とを組み合わせて前記損失密度を求めるステップと、をさらに具備し、
前記損失額の累積確率を前記損失密度から求める請求項１４に記載のオペレーショナルリスクの計量化方法。
【請求項１６】
損失額が前記閾値未満である確率を前記巨額損失密度を用いて算出し、該確率を前記補正量とする請求項１５に記載のオペレーショナルリスクの計量化方法。
【請求項１７】
取引金額によらず、損失の事故が発生する確率に相当する損失事故率と、
巨額損失の事故が発生する条件付き確率に相当する巨額損失事故率との比を前記巨額損失密度算出ステップにより算出された巨額損失密度に乗じる請求項１５に記載のオペレーショナルリスクの計量化方法。
【請求項１８】
乱数に基づく事故発生モデルを含むオペレーション品質評価モデルから、前記損失事故率を算出するステップを具備する請求項１７に記載のオペレーショナルリスクの計量化方法。
【請求項１９】
前記巨額損失密度算出ステップにおいて、
一の損失額を一の取引金額で除すことにより損失率を算出し、該損失率に対応する損失率密度を前記一の取引金額で除してなる値を求め、
前記値を複数の取引金額の全てについて求めるとともに、該複数の取引金額の全てについて求めた値の平均値を前記巨額損失密度として算出する請求項１２に記載のオペレーショナルリスクの計量化方法。
【請求項２０】
コンピュータに、
取引金額を入力する手順と、
損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を入力する手順と、
損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を前記取引金額と前記損失率密度とに基づいて算出する巨額損失密度算出手順と、
前記巨額損失密度算出手順により算出された巨額損失密度からオペレーショナルリスクを算出するリスク量算出手順と、を実行させるためのオペレーショナルリスクの計量化プログラム。
【請求項１】
取引金額を入力する手段と、
損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を入力する手段と、
損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を前記取引金額と前記損失率密度とに基づいて算出する巨額損失密度算出手段と、
前記巨額損失密度算出手段により算出された巨額損失密度からオペレーショナルリスクを算出するリスク量算出手段と、を具備するオペレーショナルリスクの計量化装置。
【請求項２】
前記オペレーショナルリスクは、バリューアットリスク及びコンディショナルバリューアットリスクのいずれかを含む請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項３】
前記巨額損失密度に基づいて、損失額の累積確率を求める手段と、
取引金額によらず、損失事故が発生する確率に相当する損失事故率を入力する手段と、
乱数と、前記損失事故率とに基づいて、複数件の損失事故を発生する手段と、
発生した損失事故の各々の損失額を前記損失額の累積確率に基づいて算出する手段と
を具備し、
算出された損失額を値の大きい順に並べ、所定の基準で選択された値を前記オペレーショナルリスクとする請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項４】
損失額が前記閾値未満であるという少額条件の下に、該損失額を確率変数としたときの条件付き確率密度に相当する少額損失密度を入力する手段と、
前記巨額損失密度に基づく補正量を乗じることにより前記少額損失密度を補正する補正手段と、
前記巨額損失密度と、前記少額損失密度とを組み合わせて前記損失密度を求める手段と、
を具備し、
前記損失額の累積確率を前記損失密度から求める請求項１乃至３のいずれかに記載のオペレーショナルリスクの計量化装置。
【請求項５】
損失額が前記閾値未満である確率を前記巨額損失密度を用いて算出し、該確率を前記補正量とする請求項４に記載のオペレーショナルリスクの計量化装置。
【請求項６】
取引金額によらず、損失の事故が発生する確率に相当する損失事故率と、
巨額損失の事故が発生する条件付き確率に相当する巨額損失事故率との比を前記巨額損失密度算出手段により算出された巨額損失密度に乗じる請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項７】
乱数に基づく事故発生モデルを含むオペレーション品質評価モデルから、前記損失事故率を算出する手段を具備する請求項６に記載のオペレーショナルリスクの計量化装置。
【請求項８】
前記巨額損失密度算出手段は、
一の損失額を一の取引金額で除すことにより損失率を算出し、該損失率に対応する損失率密度を前記一の取引金額で除してなる値を求め、
前記値を複数の取引金額の全てについて求めるとともに、該複数の取引金額の全てについて求めた値の平均値を前記巨額損失密度として算出する請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項９】
取引金額を入力する手段と、
入力された取引金額を記憶する取引金額記憶手段と、
損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を入力する手段と、
入力された損失率密度を記憶する損失率密度記憶手段と、
損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を、前記取引金額記憶手段に記憶された取引金額と、前記損失率密度記憶手段に記憶された損失率密度とに基づいて算出する巨額損失密度算出手段と、
算出された巨額損失密度を記憶する巨額損失密度記憶手段と、
前記巨額損失密度記憶手段に記憶された巨額損失密度からオペレーショナルリスクを算出するリスク量算出手段と、
算出されたリスク量を出力する出力手段と、
を具備するオペレーショナルリスクの計量化装置。
【請求項１０】
所定期間の損失額合計を横軸とし、個々の損失額合計値に対する確率密度を縦軸とする損失額合計密度関数において、所定の上側のパーセンタイル値を前記オペレーショナルリスクの値とする請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項１１】
所定期間の損失額合計を横軸とし、個々の損失額合計値に対する確率密度を縦軸とする損失額合計密度関数において、所定の上側のパーセンタイル値以上の損失額合計値の平均値を前記オペレーショナルリスクの値とする請求項１に記載のオペレーショナルリスクの計量化装置。
【請求項１２】
取引金額入力手段により取引金額を入力するステップと、
損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を損失率密度入力手段により入力するステップと、
損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を、巨額損失密度算出手段が前記取引金額と前記損失率密度とに基づいて算出する巨額損失密度算出ステップと、
前記巨額損失密度算出ステップにより算出された巨額損失密度から、損失額に対するリスク量をリスク量算出手段により算出するリスク量算出ステップと、を具備するオペレーショナルリスクの計量化方法。
【請求項１３】
前記オペレーショナルリスクは、バリューアットリスク及びコンディショナルバリューアットリスクのいずれかを含む請求項１２に記載のオペレーショナルリスクの計量化方法。
【請求項１４】
前記巨額損失密度に基づいて、損失額の累積確率を求めるステップと、
取引金額によらず、損失事故が発生する確率に相当する損失事故率を入力するステップと、
乱数と、前記損失事故率とに基づいて、複数件の損失事故を発生するステップと、
発生した損失事故の各々の損失額を前記損失額の累積確率に基づいて算出するステップと
を具備し、
算出された損失額を値の大きい順に並べ、所定の基準で選択された値を前記オペレーショナルリスクとする請求項１２に記載のオペレーショナルリスクの計量化方法。
【請求項１５】
損失額が前記閾値未満であるという少額条件の下に、該損失額を確率変数としたときの条件付き確率密度に相当する少額損失密度を入力するステップと、
前記巨額損失密度に基づく補正量を乗じることにより前記少額損失密度を補正する補正ステップと、
前記巨額損失密度と、前記少額損失密度とを組み合わせて前記損失密度を求めるステップと、をさらに具備し、
前記損失額の累積確率を前記損失密度から求める請求項１４に記載のオペレーショナルリスクの計量化方法。
【請求項１６】
損失額が前記閾値未満である確率を前記巨額損失密度を用いて算出し、該確率を前記補正量とする請求項１５に記載のオペレーショナルリスクの計量化方法。
【請求項１７】
取引金額によらず、損失の事故が発生する確率に相当する損失事故率と、
巨額損失の事故が発生する条件付き確率に相当する巨額損失事故率との比を前記巨額損失密度算出ステップにより算出された巨額損失密度に乗じる請求項１５に記載のオペレーショナルリスクの計量化方法。
【請求項１８】
乱数に基づく事故発生モデルを含むオペレーション品質評価モデルから、前記損失事故率を算出するステップを具備する請求項１７に記載のオペレーショナルリスクの計量化方法。
【請求項１９】
前記巨額損失密度算出ステップにおいて、
一の損失額を一の取引金額で除すことにより損失率を算出し、該損失率に対応する損失率密度を前記一の取引金額で除してなる値を求め、
前記値を複数の取引金額の全てについて求めるとともに、該複数の取引金額の全てについて求めた値の平均値を前記巨額損失密度として算出する請求項１２に記載のオペレーショナルリスクの計量化方法。
【請求項２０】
コンピュータに、
取引金額を入力する手順と、
損失事故における損失率を確率変数としたときの確率密度に相当する損失率密度を入力する手順と、
損失額を確率変数としたときの確率密度に相当する損失密度のうち、閾値以上の損失額に対応する巨額損失密度を前記取引金額と前記損失率密度とに基づいて算出する巨額損失密度算出手順と、
前記巨額損失密度算出手順により算出された巨額損失密度からオペレーショナルリスクを算出するリスク量算出手順と、を実行させるためのオペレーショナルリスクの計量化プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【公開番号】特開２００６−１５５４２７（Ｐ２００６−１５５４２７Ａ）
【公開日】平成１８年６月１５日（２００６．６．１５）
【国際特許分類】
【出願番号】特願２００４−３４７８０５（Ｐ２００４−３４７８０５）
【出願日】平成１６年１１月３０日（２００４．１１．３０）
【出願人】（０００００３０７８）株式会社東芝 (54,554)
【公開日】平成１８年６月１５日（２００６．６．１５）
【国際特許分類】
【出願日】平成１６年１１月３０日（２００４．１１．３０）
【出願人】（０００００３０７８）株式会社東芝 (54,554)
[ Back to top ]