(57)【特許請求の範囲】
【請求項１】
プロセッサと記憶装置とを含み、
グラフィカルモデルを用いて、確率変数の状態を予測する経年変化予測システムであって、
前記記憶装置は、前記グラフィカルモデルを保持し、
前記グラフィカルモデルは、第１時点におけるエビデンス群を入力として、前記確率変数の前記第１時点から第１期間経過後の状態を示す情報を出力とし、
前記エビデンス群のエビデンスは、前記確率変数が従う確率分布又は前記確率変数の確定状態であり、
前記プロセッサは、前記グラフィカルモデルを使用して予測処理を実行し、
前記予測処理において、
前記第１時点におけるエビデンス群を取得し、
前記取得したエビデンス群に含まれる確率分布について、当該確率分布から算出される特徴量及び／又は前記グラフィカルモデルにおいて定義された当該確率分布に対応する確率変数の特徴量、に基づいて、当該確率分布から当該確率分布に対応する確率変数の確定状態への変換を行うか否かを決定し、
前記変換の対象に決定した確率分布を、当該確率分布に対応する確率変数の確定状態に変換し、
前記変換後の前記第１時点におけるエビデンス群を前記グラフィカルモデルに入力し、  前記確率変数の前記第１時点から前記第１期間経過後の状態を示す情報を出力する、経年変化予測システム。
【請求項２】
請求項１に記載の経年変化予測システムであって、
前記プロセッサは、
前記予測処理を繰り返し、
２回目以降の予測処理において、前回の予測処理における前記第１期間経過後の状態を示す情報を前記第１時点におけるエビデンス群として取得する、経年変化予測システム。
【請求項３】
請求項１に記載の経年変化予測システムであって、
前記グラフィカルモデルにおいて定義された当該確率分布に対応する確率変数の特徴量は、当該確率変数の状態数である、経年変化予測システム。
【請求項４】
請求項１に記載の経年変化予測システムであって、
前記グラフィカルモデルにおいて定義された当該確率分布に対応する確率変数の特徴量は、当該確率変数の状態の順序性の有無である、経年変化予測システム。
【請求項５】
請求項１に記載の経年変化予測システムであって、
前記グラフィカルモデルにおいて定義された当該確率分布に対応する確率変数の特徴量は、当該確率変数が従う事前分布の偏りである、経年変化予測システム。
【請求項６】
請求項１に記載の経年変化予測システムであって、
当該確率分布から算出される特徴量は、当該確率分布に対応する確率変数の平均情報量である、経年変化予測システム。
【請求項７】
請求項１に記載の経年変化予測システムであって、
前記プロセッサは、
前記予測処理において、
前記確率分布から算出される特徴量及び／又は前記グラフィカルモデルにおいて定義された前記確率分布に対応する確率変数の特徴量、に基づいて、前記取得したエビデンス群から前記変換の対象となる確率分布の候補を示す変換候補確率分布を選択し、
前記変換候補確率分布に対して前記変換を行った場合において、前記変換後の前記第１時点におけるエビデンス群を前記グラフィカルモデルに入力し前記第１時点から前記第１期間経過後における前記確率変数の状態を示す情報を出力する、場合の計算コストを、
前記グラフィカルモデルにおいて定義された前記確率変数間の関係性と、前記取得したエビデンス群に含まれ、前記変換候補確率分布と異なる確率分布、に対応する確率変数の前記グラフィカルモデルにおいて定義された状態数と、に基づいて、算出し、
前記計算コストが閾値を超える場合、前記変換候補確率分布を再選択し、
前記計算コストが前記閾値以下である場合、前記変換候補確率分布を前記変換の対象として決定する、経年変化予測システム。
【請求項８】
請求項１に記載の経年変化予測システムであって、
前記プロセッサは、
前記予測処理において、
前記変換の対象に決定した確率分布に対応する確率変数における、前記グラフィカルモデルにおいて定義された状態の順序性の有無、に基づいて、前記変換の対象に決定した確率分布を、当該確率分布に対応する確率変数の確定状態に変換する方法を決定し、
前記変換の対象に決定した確率分布を、前記決定した方法で当該確率分布に対応する確率変数の確定状態に変換する、経年変化予測システム。
【請求項９】
請求項８に記載の経年変化予測システムであって、
前記プロセッサは、
前記予測処理において、
前記変換の対象に決定した確率分布であって、対応する確率変数の状態が順序性を有する確率分布に対して、
前記グラフィカルモデルにおいて定義された当該確率変数が取り得る状態に含まれ、前記順序性が示す順序関係において互いに隣接する状態、の間隔を算出し、
前記算出した間隔に基づいて、前記変換の対象に決定した確率分布を、当該確率分布に対応する確率変数の確定状態に変換する方法を決定する、経年変化予測システム。
【請求項１０】
請求項１に記載の経年変化予測システムであって、
前記グラフィカルモデルにおいて、前記確率変数間の関係性が定義され、
前記記憶装置は、前記定義された関係性と異なる前記確率変数間の関係性を示す関係性情報を保持し、
前記プロセッサは、
前記予測処理において、
前記関係性情報に含まれる確率変数の指定を受け付け、
前記取得したエビデンス群に含まれる確率分布について、当該確率分布に対応する確率変数と前記指定された確率変数との前記関係性情報が示す関係性に基づいて、当該確率分布から当該確率分布に対応する確率変数の確定状態への変換を行うか否かを決定する、経年変化予測システム。
【請求項１１】
経年変化予測システムが、グラフィカルモデルを用いて、確率変数の状態を予測する方法であって、
前記経年変化予測システムは、前記グラフィカルモデルを保持し、
前記グラフィカルモデルは、第１時点におけるエビデンス群を入力として、前記確率変数の前記第１時点から第１期間経過後の状態を示す情報を出力とし、
前記エビデンス群のエビデンスは、前記確率変数が従う確率分布又は前記確率変数の確定状態であり、
前記方法は、
前記経年変化予測システムが、前記グラフィカルモデルを使用して予測処理を実行し、
前記予測処理において、
前記第１時点におけるエビデンス群を取得し、
前記取得したエビデンス群に含まれる確率分布について、当該確率分布から算出される特徴量及び／又は前記グラフィカルモデルにおいて定義された当該確率分布に対応する確率変数の特徴量、に基づいて、当該確率分布から当該確率分布に対応する確率変数の確定状態への変換を行うか否かを決定し、
前記変換の対象に決定した確率分布を、当該確率分布に対応する確率変数の確定状態に変換し、
前記変換後の前記第１時点におけるエビデンス群を前記グラフィカルモデルに入力し、  前記確率変数の前記第１時点から前記第１期間経過後の状態を示す情報を出力する、方法。
【請求項１２】
経年変化予測システムに、グラフィカルモデルを用いた確率変数の状態を予測させるプログラムであって、
前記経年変化予測システムは、前記グラフィカルモデルを保持し、
前記グラフィカルモデルは、第１時点におけるエビデンス群を入力として、前記確率変数の前記第１時点から第１期間経過後の状態を示す情報を出力とし、
前記エビデンス群のエビデンスは、前記確率変数が従う確率分布又は前記確率変数の確定状態であり、
前記プログラムは、前記グラフィカルモデルを使用した予測処理を実行させ、
前記予測処理において、
前記第１時点におけるエビデンス群を取得する手順と、
前記取得したエビデンス群に含まれる確率分布について、当該確率分布から算出される特徴量及び／又は前記グラフィカルモデルにおいて定義された当該確率分布に対応する確率変数の特徴量、に基づいて、当該確率分布から当該確率分布に対応する確率変数の確定状態への変換を行うか否かを決定する手順と、
前記変換の対象に決定した確率分布を、当該確率分布に対応する確率変数の確定状態に変換する手順と、
前記変換後の前記第１時点におけるエビデンス群を前記グラフィカルモデルに入力する手順と、
前記確率変数の前記第１時点から前記第１期間経過後の状態を示す情報を出力する手順と、を含む手順を実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、経年変化予測システムに関する。
【背景技術】
【０００２】
保険事業の運営を支援する背景技術として、特開２０１４−２２５１７６号公報（特許文献１）がある。この公報には、「病態を表す確率変数に対応するノードと病態の変化に影響を与える因子を表す確率変数に対応するノードとの間の確率的依存性が有向辺又は無向辺によって定義されたグラフ構造を作成する因果・遷移構造算出部と、前記確率変数の事象空間を作成するノード生成部と、前記グラフ構造の条件付確率を算出する確率テーブル算出部と、指定された確率変数からなるグラフ構造、事象空間及び条件付確率によってモデルを再構成するモデル再構築部と、病態遷移確率及び医療費を予測する病態遷移・医療費予測部と、保健指導の対象者及び保健指導内容を選定する保健指導支援部と、を備える。」と記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１４−２２５１７６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１に記載の技術は、例えば、ある年とその翌年のデータを用いて、１年後の有病率や医療費を予測するベイジアンネットワークを構築している。なお、ベイジアンネットワーク等のグラフィカルモデルに入力されるエビデンスには、確率変数の確定状態をエビデンスとするハードエビデンスと、確率変数の状態が１つに確定していないが当該確率変数が従う確率分布が判明している場合における、当該確率分布をエビデンスとするソフトエビデンスと、がある。
【０００５】
特許文献１に記載の技術は、現在の状態をハードエビデンスとしてベイジアンネットワークに入力し、１年後の状態の予測を行う。従って、現在の状態に、（確率変数の確定状態が判明していない）確率分布が含まれている場合、特許文献１に記載の技術を用いて予測を行うためには、当該確率分布を確定状態に変換する必要がある。
【０００６】
確率分布から確定状態への変換方法として、例えば、確率分布上もっとも確率の高くなる状態、即ち最頻値を選択する方法や、確率分布をもとに期待値を計算して該当する状態を選択する方法等がある。しかし、確率分布から変換された確定状態が入力される場合、当該確率分布が入力される場合と比較して、予測精度が低下する。
【０００７】
特に、例えば、グラフィカルモデルにおける出力結果を当該グラフィカルモデルにおける次の入力とする予測、を繰り返す場合には、実情と反する結果が得られる可能性が高い。例えば、有病率や医療費は一般的に年齢とともに上昇していく傾向があるが、医療費や有病率の確率分布を確定状態に変換して経年変化予測を行うと、年々医療費や有病率が減少するという予測結果が得られてしまう。
【０００８】
これは、確率分布に偏りがあることが原因である。一般的に、ある病気にかかっている人は、かかっていない人よりも圧倒的に少ないため、当該病気にかかっているか否かについての確率分布を確定状態に変換すると、当該病気にかかってないという状態が選択されてしまう。その結果、翌年の予測を繰り返すごとにほぼ全員が病気でない健康な状態に収束してしまう。
【０００９】
また、グラフィカルモデルにソフトエビデンスが入力される場合、例えば、確率変数の状態数等に応じて計算コストが増大する。そこで本発明の一態様は、エビデンス群に含まれる確率分布から、当該確率分布から当該確率分布に対応する確率変数の確定状態への変換を行う確率分布を適切に選択することで、高精度かつ低計算コストで予測を行うことができる、経年変化予測システムを提供する。
【課題を解決するための手段】
【００１０】
上記課題を解決するため、本発明の一態様は、例えば、以下の構成を採用する。プロセッサと記憶装置とを含み、グラフィカルモデルを用いて、確率変数の状態を予測する経年変化予測システムであって、前記記憶装置は、前記グラフィカルモデルを保持し、前記グラフィカルモデルは、第１時点におけるエビデンス群を入力として、前記確率変数の前記第１時点から第１期間経過後の状態を示す情報を出力とし、前記エビデンス群のエビデンスは、前記確率変数が従う確率分布又は前記確率変数の確定状態であり、前記プロセッサは、前記グラフィカルモデルを使用して予測処理を実行し、前記予測処理において、前記第１時点におけるエビデンス群を取得し、前記取得したエビデンス群に含まれる確率分布について、当該確率分布から算出される特徴量及び／又は前記グラフィカルモデルにおいて定義された当該確率分布に対応する確率変数の特徴量、に基づいて、当該確率分布から当該確率分布に対応する確率変数の確定状態への変換を行うか否かを決定し、前記変換の対象に決定した確率分布を、当該確率分布に対応する確率変数の確定状態に変換し、前記変換後の前記第１時点におけるエビデンス群を前記グラフィカルモデルに入力し、前記確率変数の前記第１時点から前記第１期間経過後の状態を示す情報を出力する、経年変化予測システム。
【発明の効果】
【００１１】
本発明の一態様によれば、一定期間後の確率変数の状態を予測するグラフィカルモデルを用いて、高精度かつ計算コストを抑えた予測を行うことができる。
【００１２】
また、上記した以外の課題、構成及び効果は、以下の実施形態の説明により明らかにされる。
【図面の簡単な説明】
【００１３】
【図１】疾病発症予測装置の構成例を示すブロック図である。
【図２Ａ】確率変数テーブルのデータ構造の一例である。
【図２Ｂ】リンクテーブルのデータ構造の一例である。
【図２Ｃ】確率テーブルのデータ構造の一例である。
【図３】出力用変数の確率分布テーブルのデータ構造の一例である。
【図４】Ｓ／Ｈ選択テーブルのデータ構造の一例である。
【図５】混合エビデンステーブルのデータ構造の一例である。
【図６】疾病発症予測処理の一例を示すフローチャートである。
【図７】Ｓ／Ｈ選択決定処理の一例を示すフローチャートである。
【図８】Ｓ／Ｈ選択条件設定部の設定用インタフェースの例である。
【図９Ａ】確率的に従属な変数の集合に分割されたベイジアンネットワークの例である。
【図９Ｂ】経路が遮断されるパターンの例である。
【図９Ｃ】経路が遮断されないパターンの例である。
【図１０Ａ】将来の疾病発症予測と医療費予測のための入出力例である。
【図１０Ｂ】生活習慣からの将来の測定値予測のための入出力例である。
【図１０Ｃ】生活習慣の推定のための入出力例である。
【図１１Ａ】疾病発症予測装置の出力を表示するインタフェースの例である。
【図１１Ｂ】疾病発症予測装置の出力を表示するインタフェースの例である。
【発明を実施するための形態】
【００１４】
以下、添付図面を参照して本発明の実施形態を説明する。本実施形態は本発明を実現するための一例に過ぎず、本発明の技術的範囲を限定するものではないことに注意すべきである。各図において共通の構成については同一の参照符号が付されている。
【００１５】
本実施形態の経年予測システムは、確率変数が従う確率分布を入力として確率推論を行うソフトエビデンス推論を用いることで、翌年の予測モデルを使った経年変化予測を実現する。経年予測システムは、確率分布を予測モデルに入力した場合、確率変数が従う確率分布を確率変数の確定状態へ変換する必要がないため、高精度に予測を行うことができる。ひいては、経年変化予測システムは、出力された確率分布を再度予測モデルに入力する予測、を繰り返す場合における、確率変数がとる値が収束する問題を回避できる。
【００１６】
一方、ソフトエビデンスを用いた確率推論処理の特徴として、ソフトエビデンスが入力される全ての確率変数の状態数等に応じて計算コストが増大する、という問題点がある。従って、経年変化予測システムが、例えば、全ての入力用確率変数に対してソフトエビデンスを入力する場合、計算コストが肥大化してしまう。
【００１７】
そこで、本実施形態の経年変化予測システムは、各入力用確率変数に対して、例えば、各入力確率変数の特徴量及び／又は各入力用変数に対応する出力用確率変数が従う確率分布の特徴量、に従って、ソフトエビデンスかハードエビデンスのどちらを入力するかを選択する。さらに、本実施形態の経年変化予測システムは、ソフトエビデンスからハードエビデンスへの変換方法を、例えば、各入力用確率変数の特徴量に従って決定してもよい。本実施形態の経年変化予測システムは、上記処理により、予測精度を向上しつつ、計算量を抑えた予測を実現する。
【実施例１】
【００１８】
本実施例は、経年変化予測システムの一例として、健診結果、問診結果、病歴、及び診療記録などといった医療データを基に、分析対象者の疾病発生確率の所定期間ごとの経年変化を予測する疾病発症予測装置の例を説明する。本実施例では当該所定期間を１年として説明する。
【００１９】
医療データは、例えば、対象者毎の診療記録や検査値など、個人毎の医療及び健康に関する情報を含むデータである。また、例えば、医療データは、身長、体重、ＢＭＩ、血圧、コレステロール値、及び血糖値などの健康診断や問診時に測定される検査値などを含んでもよい。また、例えば、医療データは、喫煙の有無、日常の汗運動の有無、飲酒の有無、及び睡眠の状態などの生活習慣に関するデータを含んでもよい。さらに、例えば、医療データは、病歴、すなわち、医療機関を受診した際の傷病名の履歴などのデータを含んでもよい。また、例えば、医療データは、処方された医薬品、実施された診療行為、及び医療費などの診療記録に関するデータを含んでもよい。
【００２０】
図１は、本実施例の疾病発症予測装置の構成例を示す。疾病発症予測装置１００は、入力部１１１と、出力部１１２と、演算装置１１３と、メモリ１１４と、記憶媒体１１５と、を含む。入力部１１１は、例えば、マウス、キーボードなどのヒューマンインターフェースであり、ユーザからデータ等の入力を受け付ける。出力部１１２は、例えば、ディスプレイやプリンタなどであり、疾病発症予測装置１００による演算結果を出力する。記憶媒体１１５は、例えば、疾病発症予測装置１００による処理を実現する各種プログラム、及び処理結果等の各種データを格納する。
【００２１】
演算装置１１３は、例えば、ＣＰＵ、ＧＰＵなどである。演算装置１１３は、例えば、プログラムに従って動作するプロセッサ及び／又は論理回路を含み、データの入力／出力、読み込み／書き込みを行い、さらに各プログラムを実行する。以下に説明する処理及び演算は、演算装置１１３が実行する。
【００２２】
メモリ１１４は、例えば、演算装置１１３が実行する各種プログラム、及びデータを一時的にロードして記憶する。記憶媒体１１５が保持するプログラム及びデータの、一部又は全部は、メモリ１１４に格納されていてもよい。なお、各種プログラムは、例えば、ＣＤ−ＲＯＭ等のコンピュータ読み取り可能な可搬性の非一時的記憶媒体に格納されていてもよく、当該非一時的記憶媒体を扱う外部記憶装置を介して、必要に応じて、当該非一時的記憶媒体からメモリ１１４にロードされてもよい。
【００２３】
記憶媒体１１５は、プログラムである、単年予測実行部１０１、Ｓ／Ｈ選択決定部１０３、Ｓ／Ｈ変換部１０５、再帰予測実行部１０６、予測結果出力部１０７、及びエビデンス入力部１０９を含む。また、記憶媒体１１５は、データ等を格納する領域である、年次予測結果記憶部１０２、医療知識記憶部２０３、Ｓ／Ｈ選択条件記憶部２０８、Ｓ／Ｈ選択テーブル記憶部１０４、及び単年予測モデル記憶部１０８を含む。
【００２４】
プログラムは演算装置１１３によって実行されることで、定められた処理を記憶装置及び通信ポート（通信デバイス）を用いながら行う。従って、本実施形態においてプログラムを主語とする説明は、演算装置１１３を主語とした説明でもよい。若しくは、プログラムが実行する処理は、そのプログラムが動作する計算機及び計算機システムが行う処理である。
【００２５】
演算装置１１３は、プログラムに従って動作することによって、所定の機能を実現する機能部（手段）として動作する。例えば、演算装置１１３は、単年予測実行部１０１に従って動作することで単年予測実行部（単年予測実行手段）として機能し、Ｓ／Ｈ選択決定部１０３に従って動作することでＳ／Ｈ選択決定部（Ｓ／Ｈ選択決定手段）として機能する。他のプログラムについても同様である。さらに、演算装置１１３は、各プログラムが実行する複数の処理のそれぞれを実現する機能部（手段）としても動作する。計算機及び計算機システムは、これらの機能部（手段）を含む装置及びシステムである。
【００２６】
単年予測モデル記憶部１０８は、例えば、ユーザ等によって予め入力された、ある時点の状態を表す医療データを基に当該ある時点から１年後の状態を予測するベイジアンネットワークの情報、を保持する。以下、このベイジアンネットワークを単年予測モデルとも呼ぶ。単年予測モデル記憶部１０８は、後述する確率変数テーブル、リンクテーブル、及び確率テーブルを保持する。単年予測モデルは、過去の医療データから統計的に構築され、医療データ内の各項目同士の統計確率的な因果関係について記述している。
【００２７】
各項目同士の因果関係は、例えば、疾病と当該疾病に対して処方される医薬品との関係、検査値の異常と当該異常が示す疾病の兆候との関係、疾病と当該疾病発症によりもたらされる増加医療費との関係等である。
【００２８】
また、単年予測モデルは少なくとも２種類の確率変数を含み、当該少なくとも２種類の確率変数は、ある時点の状態を表す入力用確率変数と、１年後の当該状態を表す出力用確率変数と、を含む。例えば、単年予測モデルは、ある時点で疾病Ａが発症しているかという入力用確率変数を含む場合、当該ある時点から１年後に疾病Ａが発症しているかという当該入力用確率変数に対応する出力用確率変数を含む。
【００２９】
単年予測実行部１０１は、単年予測モデル記憶部１０８が保持する単年予測モデルを用いて、ある時点の状態を表す医療データ等を入力とし、当該ある時点から１年後の状態を出力とする予測処理を行う。ここで、入力となる、ある時点の状態を表す医療データ等をエビデンスと呼ぶ。エビデンスは入力用確率変数に入力される。
【００３０】
エビデンスにはハードエビデンスとソフトエビデンスの２種類が存在する。ハードエビデンスは、ある確率変数の状態が１つに確定しているときの、その確定状態を示す情報である。「疾病Ａが発症している」という状態を示す情報は、ハードエビデンスの一例である。また、例えば、「血糖値が８０ｍｇ/ｄｌ〜１００ｍｇ/ｄｌである」という状態を示す情報のように、測定値が、予め定められた複数段階の範囲のうちの特定の１つの範囲に含まれることが確定している状態を示す情報も、ハードエビデンスの一例である。
【００３１】
ソフトエビデンスは、ある確率変数の状態は１つに確定していないが、当該確率変数が従う確率分布が判明しているときの、当該確率分布を示す情報である。「疾病Ａが１５％の確率で発症している」、「血糖値が８０ｍｇ／ｄｌ以下の確率が１０％、８０ｍｇ／ｄｌ〜１００ｍｇ／ｄｌの確率が５０％、１００ｍｇ／ｄｌ以上の確率が４０％である」、等の情報はソフトエビデンスの例である。
【００３２】
年次予測結果記憶部１０２は、単年予測実行部１０１により出力された予測結果、即ち出力用確率変数それぞれの確率分布を示す確率分布テーブルを記憶する。確率分布テーブルのデータ構造については後述する。疾病発症予測装置１００が経年推移の予測を行う場合、年次予測結果記憶部１０２は、例えば年次ごとの予測結果を別々に記憶する。
【００３３】
Ｓ／Ｈ選択決定部１０３は、単年予測モデル記憶部１０８が保持する単年予測モデル内の入力用確率変数それぞれに、ソフトエビデンスを入力するかハードエビデンスを入力するかを選択し、当該選択の結果を示すＳ／Ｈ選択テーブルを作成する。Ｓ／Ｈ選択テーブルのデータ構造については後述する。さらに、Ｓ／Ｈ選択決定部１０３は、ハードエビデンスの入力を選択した場合、ソフトエビデンス或いは確率分布を、ハードエビデンスへ変換する変換方式を決定してもよい。
【００３４】
また、Ｓ／Ｈ選択決定部１０３は、予測モデル入力部２０１と、確率分布入力部２０２と、予測計算コスト評価部２０４と、Ｓ／Ｈ判定部２０５と、Ｓ／Ｈ選択テーブル出力部２０６と、Ｓ／Ｈ選択条件設定部２０７と、を含む。
【００３５】
予測モデル入力部２０１は、単年予測モデル記憶部１０８から、単年予測モデルの入力を受け付ける。確率分布入力部２０２は、年次予測結果記憶部１０２から、入力用確率変数の確率分布の入力を受け付ける。予測計算コスト評価部２０４は、単年予測モデルとＳ／Ｈ選択テーブルとに基づき、単年の予測に必要となる計算コストを評価する。Ｓ／Ｈ判定部２０５は、単年予測モデル内の入力用確率変数それぞれについてソフトエビデンスを入力するかハードエビデンスを入力するかを判定し、Ｓ／Ｈ選択テーブルを作成する。Ｓ／Ｈ選択テーブル出力部２０６は、例えば、作成されたＳ／Ｈ選択テーブルを、出力部１１２に出力する。Ｓ／Ｈ選択条件設定部２０７は、Ｓ／Ｈ判定部２０５が判定を行うための条件の設定を、例えば入力部１１１を介してユーザ等から受け付ける。
【００３６】
医療知識記憶部２０３は、例えば、予めユーザ等によって入力された、単年予測モデル内の確率変数同士の関係性を示す情報を保持する。例えば、当該関係性を示す情報は、単年予測モデルに含まれる関係性を示す情報と異なる情報を含んでもよい。また、例えば、確率変数間の相関係数等は、当該関係性の一例である。また、各確率変数間の関係の有無等も、当該関係性の一例である。医療知識記憶部２０３は、入力用確率変数同士の関係性、出力用変数同士の関係性を示す情報を保持していてもよい。
【００３７】
Ｓ／Ｈ選択条件記憶部２０８は、Ｓ／Ｈ選択条件設定部２０７が受け付けた選択条件を保持する。Ｓ／Ｈ選択テーブル記憶部１０４は、Ｓ／Ｈ選択決定部１０３により作成されたＳ／Ｈ選択テーブルを記憶する。
【００３８】
Ｓ／Ｈ変換部１０５は、Ｓ／Ｈ選択テーブル記憶部１０４が保持するＳ／Ｈ選択テーブルに基づき、年次予測結果記憶部１０２が保持する予測結果をハードエビデンスに変換する。変換されなかった予測結果は確率分布の形式のままソフトエビデンスとして扱われる。変換されたハードエビデンスとソフトエビデンスは、次に単年予測実行部１０１が予測処理を行う際の入力となる。
【００３９】
再帰予測実行部１０６は、例えば、単年予測実行部１０１、年次予測結果記憶部１０２、Ｓ／Ｈ選択決定部１０３、Ｓ／Ｈ選択テーブル記憶部１０４、Ｓ／Ｈ変換部１０５、予測結果出力部１０７、及び単年予測モデル記憶部１０８等を用いて再帰的な経年推移予測を行う。予測結果出力部１０７は、年次予測結果記憶部１０２に記録された年次ごとの予測結果を出力部１１２に出力する。エビデンス入力部１０９は、例えば、単年予測実行部１０１による初回予測のためのエビデンスの入力を、入力部１１１を介して、ユーザ等から受け付ける。
【００４０】
次に、疾病発症予測装置１００が保持する各テーブルの構造例を説明する。図２Ａ、図２Ｂ、図２Ｃは単年予測モデル、即ちベイジアンネットワークのデータ構造を示す。ベイジアンネットワークは、例えば、単年予測モデル記憶部１０８が保持する、確率変数テーブル９００、リンクテーブル９１０、及び確率テーブル９２０の３つの要素から特定される。
【００４１】
図２Ａは、確率変数テーブル９００のデータ構造の一例を示す。確率変数テーブル９００は、ベイジアンネットワークに含まれるノードである各確率変数の特徴を示す。確率変数テーブル９００は、例えば、確率変数９０１、状態数９０２、状態９０３、順序性９０４、及び入力／出力９０５、及び対応変数９０６を含む。
【００４２】
確率変数９０１は、ベイジアンネットワークに含まれるノードである確率変数を示す。状態数９０２は、確率変数９０１に格納された確率変数の取り得る値のパターン数を示す。状態９０３は当該確率変数が取り得る状態の内容を示す。順序性９０４は、当該確率変数の状態が順序性を持つか否かを示す。
【００４３】
入力／出力９０５は、当該確率変数が入力用確率であるか、出力用確率変数であるかを示す。対応変数９０６は、当該確率変数が入力用確率変数であれば当該確率変数に対応する出力用確率変数を、当該確率変数が出力用確率変数であれば当該確率変数に対応する入力用確率変数を示す。当該確率変数に対応する入力用／出力用変数が存在しない場合、対応変数９０６は、例えば、Ｎｕｌｌ値や、「なし」という値、等を保持する。
【００４４】
図２Ｂは、リンクテーブル９１０のデータ構造の一例を示す。リンクテーブル９１０は、ベイジアンネットワークに含まれる確率変数同士の因果関係を示す。リンクテーブル９１０は、例えば、リンク９１１、始点９１２、及び終点９１３を含む。リンク９１１は、始点９１２に格納された確率変数と終点９１３に格納された確率変数との間に因果関係があることを示す。なお、２つの確率変数間の関係性において、始点となる確率変数を親変数、終点となる確率変数を子変数と呼ぶ。
【００４５】
図２Ｃは、確率テーブル９２０のデータ構造の一例を示す。確率テーブル９２０は、ベイジアンネットワークに含まれる確率変数における、当該確率変数の親変数の状態に応じた、確率分布を示す。確率テーブル９２０は、例えば、テーブル９２１、子変数９２２、親変数９２３、及び確率分布９２４を含む。テーブル９２１は、確率テーブル９２０の各レコードを識別する情報を示す。子変数９２２は子変数を示す。親変数９２３は親変数を示す。
【００４６】
確率分布９２４は、親変数９２３に格納された親変数の状態それぞれに応じた、子変数９２２に格納された子変数が従う確率分布それぞれ、を行とする行列を示す。親変数９２３に「なし」が格納されている場合、確率分布９２４は子変数が従う事前分布を示す。なお、子変数に対応する親変数の数は、０個であってもよいし、２個以上存在してもよい。
【００４７】
図３は、予測結果出力部１０７が出力する、出力用変数の確率分布テーブルのデータ構造の一例を示す。確率分布テーブル１０００は、例えば、確率変数１００１、及び確率分布１００２を含む。確率変数１００１は、出力用確率変数を示す。確率分布１００２は、当該出力用確率変数が従う確率分布を示す。
【００４８】
図４は、Ｓ／Ｈ選択テーブルのデータ構造の一例を示す。Ｓ／Ｈ選択テーブル１１００は、例えば、確率変数１１０１、Ｓ／Ｈ選択１１０２、及び変換方法１１０３を含む。確率変数１１０１は、入力用確率変数を示す。Ｓ／Ｈ選択１１０２は、当該入力用確率変数に、ソフトエビデンスとハードエビデンスのどちらが入力されるかを示す。なお、「ソフト」はソフトエビデンスを、「ハード」はハードエビデンスを、それぞれ示す。変換方法１１０３は、Ｓ／Ｈ選択１１０２が「ハード」を格納する場合における、確率分布からハードエビデンスへの変換方法かを示す。
【００４９】
図５は、混合エビデンステーブルのデータ構造の一例を示す。混合エビデンステーブル１２００は、例えば、確率変数１２０１、Ｓ／Ｈ選択１２０２、状態１２０３、及び確率分布１２０４を含む。確率変数１２０１は、入力用確率変数を示す。Ｓ／Ｈ選択１２０２は、当該入力用確率変数にソフトエビデンスとハードエビデンスのどちらが入力されるかを示す。
【００５０】
状態１２０３は、Ｓ／Ｈ選択１２０２が「ハード」を格納する場合における、当該入力用確率変数の確定状態を示す。Ｓ／Ｈ選択１２０２が「ソフト」を保持する場合、状態１２０３は「なし」を保持する。確率分布１２０４は、Ｓ／Ｈ選択１２０２が「ソフト」を保持する場合、当該入力用確率変数が従う確率分布を示す。Ｓ／Ｈ選択１２０２が「ハード」を保持する場合、確率分布１２０４は「なし」を保持する。
【００５１】
次に、疾病発症予測装置１００の動作例を説明する。図６は、疾病発症予測装置１００によるＮ年後までの疾病発生率の経年変化の予測処理の一例を示す。疾病発症予測装置１００は、単年予測モデルを繰り返し使うことでＮ年後の予測を行う。また、疾病発症予測装置１００は、単年予測モデル内の確率変数の特徴量に基づき、ソフトエビデンスとハードエビデンスを使い分けることで計算量の増大を抑えながら予測精度を向上する。ひいては、疾病発症予測装置１００は、予測処理を繰り返すことによる、予測結果の間違った値への収束、を防ぐ。
【００５２】
まずステップＳ３０１において、予測モデル入力部２０１が、単年予測モデルの情報、即ち確率変数テーブル９００、リンクテーブル９１０、及び確率テーブル９２０の情報の入力を受け付ける。次に、ステップＳ３０２において、エビデンス入力部１０９が、予測の初期入力である、予測対象者０年目の、１以上の確率変数の確定状態からなる集合Ｈ（０）の入力を受け付ける。Ｈ（０）の各要素はハードエビデンスであるとする。
【００５３】
次に、ステップＳ３０３において、単年予測実行部１０１が単年予測モデルを用いて１回目の予測を行う。１回目の予測における入力はＨ（０）であり、得られる出力は１年後の予測結果を示す１以上の確率分布からなる集合Ｓ（１）である。Ｓ（１）に含まれる確率分布それぞれは確率分布テーブル１０００に格納される。また、ステップＳ３０４において、単年予測実行部１０１は、単年予測モデルの入力用確率変数に対して、入力がない状態で予測を行い、１以上の事前分布からなる集合Ｓｐを出力として得る。
【００５４】
ステップＳ３０５において、再帰予測実行部１０６が、経過年数ｙをｙ＝１と設定し、再帰的な予測を開始する。ステップＳ３０６において、Ｓ／Ｈ選択決定部１０３が単年予測モデル、Ｓ（ｙ）、Ｓｐに基づきＳ／Ｈ選択テーブル１１００を作成する。Ｓ／Ｈ選択テーブル１１００の具体的な作成処理については、後述する。
【００５５】
ステップＳ３０７において、Ｓ／Ｈ変換部１０５がＳ／Ｈ選択テーブル１１００に基づき、ｙ年目の予測結果であるＳ（ｙ）を、ｙ年目の混合エビデンスＭ（ｙ）へと変換する。混合エビデンスＭ（ｙ）は入力用確率変数それぞれに対して入力される、ソフトエビデンスまたはハードエビデンスの集合であり、Ｍ（ｙ）の情報は混合エビデンステーブル１２００に格納される。
【００５６】
具体的には、ステップＳ３０７において、Ｓ／Ｈ変換部１０５は、Ｓ／Ｈ選択テーブル１１００の確率変数１１０１及びＳ／Ｈ選択１１０２を、混合エビデンステーブル１２００の確率変数１２０１、及びＳ／Ｈ選択１２０２にそれぞれ格納する。
【００５７】
Ｓ／Ｈ変換部１０５は、Ｓ／Ｈ選択１２０２が「ソフト」である入力用確率変数それぞれについて、混合エビデンステーブル１２００の、状態１２０３に「なし」を、確率分布１２０４にＳ（ｙ）の当該入力用確率数が従う確率分布を、それぞれ格納する。
【００５８】
Ｓ／Ｈ変換部１０５は、Ｓ／Ｈ選択１２０２が「ハード」である入力用確率変数それぞれについて、Ｓ（ｙ）の当該入力用確率変数が従う確率分布を、変換方法１１０３で指定された方法によって確定状態に変換する。Ｓ／Ｈ変換部１０５は、Ｓ／Ｈ選択１２０２が「ハード」である入力用確率変数それぞれについて、混合エビデンステーブル１２００の、状態１２０３に当該確定状態を、確率分布１２０４に「なし」を、それぞれ格納する。
【００５９】
ステップＳ３０８において、単年予測実行部１０１は、混合エビデンスＭ（ｙ）を入力とし、ｙ＋１年目の予測を行う。ステップＳ３０９において、再帰予測実行部１０６は、ｙ＋１＝Ｎであるか否かを判定する。再帰予測実行部１０６がｙ＋１＜Ｎであると判定した場合（Ｓ３０９：ＮＯ）、ステップＳ３１０に進む。ステップＳ３１０において、再帰予測実行部１０６は、ｙに１を加えて再度ステップＳ３０６に戻って処理を継続する。
【００６０】
ステップＳ３０９において、再帰予測実行部１０６がｙ＋１＝Ｎであると判定した場合（Ｓ３０９：ＹＥＳ）、即ち目的のＮ年目の予測結果が得られた場合、ステップＳ３１１に進む。ステップＳ３１１において、予測結果出力部１０７が予測結果Ｓ（１）、Ｓ（２）、…、Ｓ（Ｎ）を出力し、疾病発症予測装置１００による疾病発生予測処理は終了する。  なお、ステップＳ３０２において、Ｈ（０）の各要素は全てハードエビデンスであるとしたが、確率変数が従う確率分布、即ちソフトエビデンスがＨ（０）に含まれていてもよい。このとき、Ｓ／Ｈ選択決定部１０３は、Ｈ（０）に含まれるソフトエビデンスそれぞれを、ハードエビデンスに変換するか否かを判定し、当該判定結果に基づいてＳ／Ｈ変換を行ってもよい。
【００６１】
次に、Ｓ／Ｈ選択決定部１０３の動作例を説明する。図７は、Ｓ／Ｈ選択決定部１０３の動作例を示す。ステップＳ４０１において、予測モデル入力部２０１が、単年予測モデル記憶部１０８から単年予測モデル、即ち確率変数テーブル９００、リンクテーブル９１０、及び確率テーブル９２０の情報を受け取る。ステップＳ４０２において、確率分布入力部２０２が、年次予測結果記憶部１０２からＳ（ｙ）に含まれる確率分布、即ち確率分布テーブル１０００の情報を受け取る。
【００６２】
ステップＳ４０３において、Ｓ／Ｈ判定部２０５は、Ｓ／Ｈ選択条件記憶部２０８が保持するＳ／Ｈ選択条件を読み込む。Ｓ／Ｈ選択条件は、例えば、予めユーザにより設定されたＳ／Ｈ選択テーブル１１００作成のための条件である。Ｓ／Ｈ選択条件は、例えば、１回の予測にかかる計算コストの上限、特に精度よく予測したい出力用確率変数等の、ソフトエビデンスを選択する優先度を定めるためのルールを含む。
【００６３】
図８は、Ｓ／Ｈ選択条件設定部２０７がＳ／Ｈ選択条件の設定を受け付けるためのユーザインタフェースの例である。設定用ユーザインタフェース６００は、例えば、入力受付領域６０１〜６０２、チェックボックス６０３、及び確定ボタン６０４を含む。入力受付領域６０１は、推論速度の入力を受け付ける。当該推論速度は、ステップＳ３０６〜ステップＳ３１０の１回のループ処理における推論速度を示してもよいし、ステップＳ３０１〜ステップＳ３１１の処理における推論速度を示してもよい。
【００６４】
入力受付領域６０２は、出力用確率変数のうち、特に精度よく予測したい項目の入力を受け付ける。チェックボックス６０３は、検査値を全てハードエビデンスに変換するためのチェックボックスである。確定ボタン６０４は、入力受付領域６０１〜６０２、及びチェックボックス６０３に入力されたＳ／Ｈ選択条件を確定するためのボタンである。
【００６５】
図７の説明に戻る。ステップＳ４０４において、Ｓ／Ｈ判定部２０５は、確率変数テーブル９００に含まれる入力用確率変数それぞれをＳ／Ｈ選択テーブル１１００の確率変数１１０１に格納し、例えば、Ｓ／Ｈ選択１１０２の全ての値を「ハード」に設定する。また、Ｓ／Ｈ判定部２０５は、確率変数テーブル９００を参照し、入力用確率変数に対応する出力用確率分布が従う確率分布をＳ（ｙ）から特定する。
【００６６】
Ｓ／Ｈ判定部２０５は、確率変数テーブル９００を参照し、確率変数１１０１に格納された確率変数それぞれの特徴に従って、変換方法１１０３それぞれを決定する。例えば、当該確率変数の状態が順序性を有する場合、Ｓ／Ｈ判定部２０５は、変換方法１１０３を期待値に決定する。なお、変換方法１１０３が期待値である場合、Ｓ／Ｈ変換部１０５は、確率分布を、当該確率分布の期待値に最も近い状態に変換する。
【００６７】
また、当該確率変数の状態が順序性を有する場合、例えば、Ｓ／Ｈ判定部２０５は、当該順序性が示す順序関係において、隣接する値同士の差を全て算出してもよい。このとき、当該差の分散が所定以上である、又は最大の当該差から最小の当該差を引いた差が所定以上であれば、当該確率変数の状態が順序性を有する場合であっても、Ｓ／Ｈ判定部２０５は、変換方法１１０３を最頻値に決定してもよい。
【００６８】
また、例えば、当該確率変数の状態が順序性を有しない場合、Ｓ／Ｈ判定部２０５は、変換方法１１０３を最頻値に決定する。Ｓ／Ｈ判定部２０５が上述のように変換方法を決定することにより、予測精度が向上する。なお、入力用確率変数のハードエビデンスへの変換方法は、例えば、確率変数テーブル９００上に予め定められていてもよく、Ｓ／Ｈ判定部２０５は当該予め定められた変換方法を変換方法１１０３に格納してもよい。
【００６９】
ステップＳ４０５において、Ｓ／Ｈ判定部２０５は、例えば、選択済み確率変数のリストとして、空のリストを作成する。ステップＳ４０６において、Ｓ／Ｈ判定部２０５は、例えば、単年予測モデル、Ｓ（ｙ）に含まれる確率分布、及びＳ／Ｈ選択条件に基づき、選択済み確率変数のリストに含まれない入力用確率変数から、ソフトエビデンスを入力する優先度の高い入力用確率変数を選択する。Ｓ／Ｈ判定部２０５は、Ｓ／Ｈ選択テーブル１１００の、当該選択した入力用確率変数に対応するＳ／Ｈ選択１１０２の値を「ソフト」に設定する。この時点でＳ／Ｈ選択１１０２の値が「ハード」である確率変数が従う確率分布は、変換候補確率分布である。なお、優先度の最も高い複数の入力用確率変数が存在する場合、Ｓ／Ｈ判定部２０５は、例えば、当該複数の入力用確率変数からランダムに１つ入力用確率変数を選択すればよい。
【００７０】
なお、Ｓ／Ｈ判定部２０５は、例えば、入力用確率変数の状態数、入力用確率変数が従う事前分布の偏り、及び入力用確率変数の状態の順序性等の、入力用確率変数の特徴量に基づいて、当該優先度を定める。また、Ｓ／Ｈ判定部２０５は、入力用確率変数の平均情報量等の入力用確率変数に対応する確率分布の特徴量に基づいて、当該優先度を定めてもよい。Ｓ／Ｈ判定部２０５は、例えば、入力用確率変数の状態数が小さいほど、当該優先度を高くすることで、計算コストの増加量が小さいソフトエビデンスを選択することができる。
【００７１】
また、Ｓ／Ｈ判定部２０５は、例えば、入力用確率変数の平均情報量、又は当該入力用確率変数が従う事前分布の偏りが大きいほど、当該入力用確率変数の状態に順序性がない場合に、当該優先度を高く定めることにより、予測精度向上への寄与率が高いソフトエビデンスを選択することができる。
【００７２】
なお、Ｓ／Ｈ判定部２０５は、例えば、入力用確率変数の状態数、及び入力用確率変数の状態の順序性等の特徴量を確率変数テーブル９００から取得する。また、Ｓ／Ｈ判定部２０５は、例えば、確率テーブル９２０から取得した確率分布から、入力用確率変数の事前分布の偏りを算出する。また、Ｓ／Ｈ判定部２０５は、確率分布テーブル１０００から出力用確率変数の確率分布を取得し、当該確率分布の平均情報量を算出し、当該算出した平均情報量を、当該出力用確率変数に対応する入力用確率用変数の平均情報量とすればよい。
【００７３】
例えば、薬Ａの投薬の有無という確率変数Ｘと、検査値Ｂの値という確率変数Ｙと、があった場合、Ｓ／Ｈ判定部２０５は、確率変数Ｘの優先度を高くする。確率変数Ｘの状態数は「あり」、「なし」の２つと少なく、さらに「あり」の確率が極端に小さいため平均情報量が大きい。
【００７４】
これに対して、確率変数Ｙは、数値データであるため状態数が多く、さらに確率分布の偏りもそれほど極端ではなく平均情報量が小さい。さらに確率変数Ｙの状態には順序性があるため、期待値を計算することができる。従って、Ｓ／Ｈ判定部２０５は、確率変数Ｘの優先度を確率変数Ｙの優先度より高くすることにより、予測精度低下を抑制することができ、かつ計算コストを低下させることができると見込めるからである。
【００７５】
なお、Ｓ／Ｈ判定部２０５は、複数の特徴量を用いて優先度を定める場合、例えば、当該複数の特徴量それぞれを変数とする所定の関数によって算出された値が大きい順に優先度を高く設定すればよい。当該複数の特徴量が状態数、平均情報量、及び順序性の有無からなる場合、当該所定の関数は、例えば、状態数、及び順序性の有無（例えば順序性ありを１、順序性なしを０とする）の減少関数、かつ平均情報量の増加関数であればよい。
【００７６】
また、Ｓ／Ｈ判定部２０５は、複数の特徴量を用いて優先度を定める場合、例えば、当該複数の特徴量間における所定の優先順位に基づき、当該複数の特徴量をキーとして入力用確率変数をソートすることにより、優先度を定めてもよい。当該複数の特徴量が状態数、平均情報量、及び順序性の有無からなり、特徴量間の優先順位が高い順に状態数、平均情報量、及び順序性の有無である場合を考える。このとき、Ｓ／Ｈ判定部２０５は、例えば、状態数を第１キーとする昇順、平均情報量を第２キーとする降順、順序性の有無を第３キーとする昇順で、入力用確率変数をソートし、当該ソート後の順番を優先度とする。
【００７７】
また、Ｓ／Ｈ選択条件に特に精度よく予測したい項目が指定されている場合、Ｓ／Ｈ判定部２０５は、例えば、当該項目に対応する確率変数、及び／又は医療知識記憶部２０３に格納された情報が示す当該項目に関連する項目に対応する確率変数の優先度を上げてもよい。具体的には、Ｓ／Ｈ判定部２０５は、例えば、優先度を上げる対象として選択された確率変数の優先度を最も高くしてもよいし、また上述の所定の関数の値によって優先度を算出する場合、当該確率変数においては関数から算出された値に所定値を加えてもよい。Ｓ／Ｈ判定部２０５は、医療知識記憶部２０３に格納された情報を使用することにより、ベイジアンネットワークが示す関係性と異なる関係性を用いることができるため、予測精度が向上する。
【００７８】
ステップＳ４０７において、予測計算コスト評価部２０４は、単年予測モデルとその時点でのＳ／Ｈ選択テーブル１１００に基づき、予測処理を行うための計算コストを評価する。ステップＳ４０８において、予測計算コスト評価部２０４は、評価した計算コストがＳ／Ｈ選択条件で定められた計算コストの上限以内か否かを判定する。
【００７９】
予測計算コスト評価部２０４が、評価した計算コストが当該上限以内でないと判定した場合（Ｓ４０８：ＮＯ）、ステップＳ４０９において、Ｓ／Ｈ判定部２０５は、Ｓ／Ｈ選択テーブル１１００においてステップＳ４０６で選択した入力用確率変数のＳ／Ｈ選択１１０２を「ハード」に戻す処理を行い、ステップＳ４１０に遷移する。また、予測計算コスト評価部２０４が、評価した計算コストが当該上限以内でないと判定した場合（Ｓ４０８：ＹＥＳ）、ステップＳ４１１に遷移する。
【００８０】
次にステップＳ４１０において、Ｓ／Ｈ判定部２０５は、ステップＳ４０６で選択した入力用確率変数を選択済み確率変数のリストに加える。ステップＳ４１１において、Ｓ／Ｈ判定部２０５は、選択済み確率変数のリストに含まれない入力用確率変数があるか否かを判定する。
【００８１】
Ｓ／Ｈ判定部２０５が、選択済み確率変数のリストに含まれない入力用確率変数があると判定した場合（Ｓ４１１：ＹＥＳ）は、ステップＳ４０６に戻り処理を継続する。Ｓ／Ｈ判定部２０５が、選択済み確率変数のリストに含まれない入力用確率変数がないと判定した場合（Ｓ４１１：ＮＯ）、ステップＳ４１３において、Ｓ／Ｈ選択テーブル出力部２０６は、Ｓ／Ｈ選択テーブル１１００をＳ／Ｈ選択テーブル記憶部１０４に出力し、Ｓ／Ｈ選択決定部１０３の処理を終了する。
【００８２】
なお、Ｓ／Ｈ選択決定部１０３は、全ての予測年次について、Ｓ／Ｈ選択テーブル１１００を作成しなくてもよい。例えば、Ｓ／Ｈ選択決定部１０３は、Ｍ年目（Ｍは１≦Ｍ＜Ｎを満たす自然数）までのＳ／Ｈ選択テーブル１１００を作成し、Ｓ／Ｈ変換部１０５はＭ年目のＳ／Ｈ選択テーブル１１００を用いてＭ＋１年目以降のＳ／Ｈ変換を行ってもよい。
【００８３】
また、特に状態数、事前確率分布の偏り、順序性の有無は、どの年次においても不変である。従って、Ｓ／Ｈ判定部２０５は、これらの特徴量から選択した１以上の特徴量に基づいて、Ｓ／Ｈ選択テーブル１１００を作成する場合、１年目のみＳ／Ｈ選択テーブル１１００を作成してもよい。また、Ｓ／Ｈ選択決定部１０３は、予測年次ごとに、Ｓ／Ｈ選択テーブル１１００のＳ／Ｈ選択１１０２及び変換方法１１０３の決定方法を変更してもよい。
【００８４】
以下、予測処理を行う場合の計算コストを評価するステップＳ４０８の処理の一例について説明する。図９Ａは、確率的に従属な変数の集合に分割されたベイジアンネットワークの例を示す。なお、図９Ａ、並びに後述する図９Ｂ及び図９Ｃにおいて、各円はノードを、各矢印は各ノード間の因果関係を、各円がドットで塗りつぶされているか否かは各ノードに入力されるエビデンスの種類を示す。ドットで塗りつぶされていないノードは、ソフトエビデンスが入力されるノード、又はエビデンスが入力されない第１種ノードを示し、ドットで塗りつぶされたノードはハードエビデンスが入力される第２種ノードを示す。
【００８５】
ベイジアンネットワーク５０１は、リンクテーブル９１０で定義された構造を持ち、現時点でのＳ／Ｈ選択テーブル１１００に従って、各ノードに入力されるエビデンスの種類が特定されている、ベイジアンネットワークの一例である。予測計算コスト評価部２０４は、各ノードの種類、及び各ノードを結ぶ経路に従って、ベイジアンネットワーク５０１に含まれる確率変数（ハードエビデンス以外が入力されるノード）であって、確率的に従属な確率変数、からなる集合である確率変数集合５０２それぞれを特定する。確率変数集合５０２の特定方法は後述する。
【００８６】
予測計算コスト評価部２０４は、例えば、確率変数集合５０２それぞれに対して、ソフトエビデンスが入力される変数の状態数の積を算出し、当該積を当該確率変数集合５０２の計算コストとする。予測計算コスト評価部２０４は、例えば、算出した確率変数集合５０２の計算コストの総和を、全ての確率変数を推論するための計算コスト、即ちステップＳ４０８における計算コストとする。
【００８７】
以下、確率変数集合を特定する方法を説明する。ベイジアンネットワーク５０３は、リンクの経路で結ばれた確率変数５０４と確率変数５０５とを含む。ベイジアンネットワークは、上述した経路のみを含むため、確率変数５０４と確率変数５０５は、その他の確率変数に対して互いに独立である。従って、予測計算コスト評価部２０４は、ベイジアンネットワーク５０３を、確率変数５０４と確率変数５０５からなる集合を確率変数集合として特定する。
【００８８】
つまり、予測計算コスト評価部２０４は、経路で結ばれた確率変数からなる集合を確率変数集合として特定する作業を繰り返す。但し、確率変数間の経路が存在する場合であっても、経路中のハードエビデンスが入力されるノードの有無と、グラフの構造と、によっては、当該経路が遮断される場合がある。確率変数間の経路が遮断されている場合は、当該確率変数は互いに独立である。
【００８９】
図９Ｂは、経路が遮断されるパターンの例を示す。パターン５１１〜５１３は、いずれも破線部分で経路が遮断される。パターン５１１はｔａｉｌ−ｔｏ−ｔａｉｌと呼ばれる、親ノードから複数の子ノードへの経路を含むパターンである。パターン５１１における親ノードは第２種ノード、子ノードは第１種ノードである。
【００９０】
パターン５１２は、ｈｅａｄ−ｔｏ−ｔａｉｌと呼ばれる、親ノードから子ノード、及び子ノードから孫ノードへの経路を含むパターンである。パターン５１２における親ノードは第１種ノード、子ノードは第２種ノード、孫ノードは第１種ノードである。パターン５１３は、ｈｅａｄ−ｔｏ−ｈｅａｄと呼ばれる、複数の親ノードから子ノードへの経路を含むパターンである。パターン５１３における親ノード及び子ノードは、第１種ノードである。
【００９１】
図９Ｃは、経路が遮断されないパターンの例を示す。ｔａｉｌ−ｔｏ−ｔａｉｌのパターン５２１における親ノード及び子ノードは第１種ノードである。ｈｅａｄ−ｔｏ−ｔａｉｌのパターン５２２における、親ノード、子ノード、及び孫ノードは第１種ノードである。ｈｅａｄ−ｔｏ−ｈｅａｄのパターン５２３及びパターン５２４は、第１種ノードである親ノード、子孫ノードに第２種ノードを含む。パターン５２３及びパターン５２４のように、ｈｅａｄ−ｔｏ−ｈｅａｄのパターンであって、子ノード又は子ノードの子孫ノードのいずれかが第２種ノードであるパターンにおいては、遮断が解かれる。
【００９２】
次に、本実施例の疾病発症予測装置１００における入出力例を説明する。図１０Ａは、将来の疾病発症予測と医療費予測に本実施例の疾病発症予測装置１００を適用した場合の入出力例を示す。疾病発症予測装置１００は、例えば、体重、身長等の現在の測定値、生活習慣、及び病歴等を入力として、各疾病の将来の発症確率などの確率、及び医療費の期待値などを出力してもよい。
【００９３】
図１０Ｂは、生活習慣からの将来の測定値予測に本実施例の疾病発症予測装置１００を適用した場合の入出力例を示す。疾病発症予測装置１００は、例えば、体重、身長等の現在の測定値、及び生活習慣を入力として、体重、血圧等の将来の予測値を出力してもよい。なお、疾病発症予測装置１００は、体重、血圧などの測定値の予測値を、具体的な数値でなく、複数段階で分割された測定値の範囲として出力してもよい。
【００９４】
図１０Ｃは、生活習慣の推定に本実施例の疾病発症予測装置１００を適用した場合の入出力例を示す。疾病発症予測装置１００は、例えば、体重、身長等の現在の測定値、及び病歴を入力として、現在の生活習慣を出力してもよい。
【００９５】
図１１Ａは、予測結果出力部１０７によって出力部１１２に表示させた個人における予測結果出力画面の例である。予測結果出力画面７００は、例えば、タブ７０１、及び表示領域７０２〜７０８を含む。表示領域７０２〜７０８それぞれは、タブ７０１で指定された年次の値を表示する。表示領域７０２は、当該個人の年間医療費の期待値を表示する。表示領域７０３〜７０８それぞれは、当該個人の各疾病の発生率を表示する。
【００９６】
図１１Ｂは、複数人の団体の予測結果出力画面の例である。予測結果出力画面８００は、例えば、タブ８０１、及び表示領域８０２〜８０８を含む。表示領域８０２〜８０８それぞれは、タブ８０１で指定された年次の値を表示する。表示領域７０２は、当該団体の年間医療費の期待値を表示する。表示領域７０３〜７０８それぞれは、当該団体の各疾病の発生予測人数を表示する。
【００９７】
また、出力内容は、確率推論による予測の情報に限定されず、Ｓ／Ｈ選択テーブル１１００の情報、ソフトエビデンスを入力する優先度が高いが計算コストの制限によりソフトエビデンスが入力されなかった確率変数の情報なども出力部１１２に表示してもよい。
【００９８】
上述したように、本実施例の疾病発症予測装置１００は、分析対象者に関する既知の医療データを入力し、ベイジアンネットワークである単年予測モデル上での確率推論によって１年後の医療データの状態を予測する。疾病発症予測装置１００は、予測された１年後の状態を入力として再度単年予測モデルにより２年後の医療データの状態を予測し、さらに再帰的にＮ回の予測を行うことでＮ年後までの医療データの年次推移を予測することができる。
【００９９】
その際、疾病発症予測装置１００は、確率変数の特徴量に基づいて、入力用確率変数それぞれにソフトエビデンス又はハードエビデンスのいずれか入力するかを選択することで、小さい計算コストで、高精度に予測を行うことができる。特に、疾病発症予測装置１００は、確率分布に偏りがある医療データの値が予測処理の繰り返しによって収束してしまうことを回避することができる。
【０１００】
なお、本実施例の疾病発症予測装置１００は、ベイジアンネットワークを用いた予測処理を行っているが、ベイジアンネットワークの代わりに、例えば、マルコフネットワーク等の他のグラフィカルモデルを用いて同様の予測処理を行ってもよい。
【０１０１】
なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることも可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
【０１０２】
また、上記の各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリや、ハードディスク、ＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）等の記録装置、または、ＩＣカード、ＳＤカード、ＤＶＤ等の記録媒体に置くことができる。
【０１０３】
また、制御線や情報線は説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。
【符号の説明】
【０１０４】
１００  疾病発症予測装置、１０１  単年予測実行部、１０２  年次予測結果記憶部、１０３  Ｓ／Ｈ選択決定部、１０４  Ｓ／Ｈ選択テーブル記憶部、１０５  Ｓ／Ｈ変換部、１０８  単年予測モデル記憶部、１１３  演算装置、１１４  メモリ、１１５  記憶媒体、２０１  予測モデル入力部、２０２  確率分布入力部、２０３  医療知識記憶部、２０４  予測計算コスト評価部、２０５  Ｓ／Ｈ判定部、２０６  Ｓ／Ｈ選択テーブル出力部
【図１】
0006301853.tif 000002
【図２Ａ】
0006301853.tif 000003
【図２Ｂ】
0006301853.tif 000004
【図２Ｃ】
0006301853.tif 000005
【図３】
0006301853.tif 000006
【図４】
0006301853.tif 000007
【図５】
0006301853.tif 000008
【図６】
0006301853.tif 000009
【図７】
0006301853.tif 000010
【図８】
0006301853.tif 000011
【図９Ａ】
0006301853.tif 000012
【図９Ｂ】
0006301853.tif 000013
【図９Ｃ】
0006301853.tif 000014
【図１０Ａ】
0006301853.tif 000015
【図１０Ｂ】
0006301853.tif 000016
【図１０Ｃ】
0006301853.tif 000017
【図１１Ａ】
0006301853.tif 000018
【図１１Ｂ】
0006301853.tif 000019