(57)【要約】【課題】効率的な診断を行うための適切な情報を操作者に提示すること。【解決手段】実施形態に係る医用情報処理装置は、処理部と、出力部とを備える。処理部は、判定対象の画像データが疾患を有する確率と、その確率の信頼度を示す指標である信頼区間とを算出する処理を行う。出力部は、確率と、確率の信頼区間とに基づく、疾患に関する情報を出力する。【選択図】  図２
2020149684.tif 000002
【特許請求の範囲】
【請求項１】
判定対象の画像データが疾患を有する確率と、その確率の信頼度を示す指標である信頼区間とを算出する処理を行う処理部と、
前記確率と、前記確率の信頼区間とに基づく、疾患に関する情報を出力する出力部と、
を備える医用情報処理装置。
【請求項２】
前記処理部は、統計情報取得部と信頼区間算出部とを備え、
前記統計情報取得部は、複数の画像データの特徴に基づき集計された疾患の頻度に関する統計情報を取得し、
前記信頼区間算出部は、前記判定対象の画像データが疾患を有する確率と、その確率の信頼度を示す指標である信頼区間とを、前記統計情報が示す頻度分布を利用して判定対象の計測値を含む任意の範囲から求め、
前記出力部は、前記確率と、前記確率の信頼区間とに基づく、疾患に関する情報を出力する、
請求項１に記載の医用情報処理装置。
【請求項３】
前記信頼区間算出部は、前記判定対象の特徴が所属する特徴量の範囲を計数範囲として算出する、
請求項２に記載の医用情報処理装置。
【請求項４】
前記信頼区間算出部は、前記複数の画像データの中で、前記計数範囲に所属する症例数を群ごとに求め、少なくとも１つの群について、所属する確率の平均値と、信頼区間とを算出する、
請求項３に記載の医用情報処理装置。
【請求項５】
前記出力部は、前記疾患に関する情報として、前記確率と、前記確率の信頼区間との関係から前記判定対象の疾患に関する判定結果を出力する、
請求項１乃至４のいずれか１項に記載の医用情報処理装置。
【請求項６】
前記出力部は、前記疾患に関する情報として、前記確率と、前記確率の信頼区間とを出力する、
請求項１乃至４のいずれか１項に記載の医用情報処理装置。
【請求項７】
前記出力部は、前記信頼区間の判定条件に応じた判定結果を出力する、
請求項５又は６に記載の医用情報処理装置。
【請求項８】
複数の画像データの特徴に基づき疾患の頻度に関する統計情報を収集する統計情報収集部をさらに備える、
請求項１乃至７のいずれか１項に記載の医用情報処理装置。
【請求項９】
前記処理部は、尺度算出部と抽出部とを備え、
前記尺度算出部は、前記判定対象の画像データを複数の画像パターンに分類して、各画像パターンの疾患確率を示す前記信頼区間を識別妥当性尺度として算出し、
前記抽出部は、前記複数の画像パターンの各画像パターンの識別妥当性尺度に基づいて、前記複数の画像パターンから、疾患の識別が行いやすい一部の画像パターンを抽出し、
前記出力部は、前記一部の画像パターンを優先して表示部に表示させる、
請求項１乃至８のいずれか１項に記載の医用情報処理装置。
【請求項１０】
前記尺度算出部は、前記判定対象の画像データに基づいて識別妥当性尺度を生成する学習済みモデルに対して、前記判定対象の画像データを入力することで、前記識別妥当性尺度を生成する、
請求項９に記載の医用情報処理装置。
【請求項１１】
前記尺度算出部は、前記判定対象の画像データの関心領域を設定し、前記関心領域に基づいて、前記判定対象の画像データを複数の画像パターンに分類する、
請求項９又は１０に記載の医用情報処理装置。
【請求項１２】
前記出力部は、前記一部の画像パターンの名称を文字情報として前記表示部に表示させる、
請求項９乃至１１のいずれか１項に記載の医用情報処理装置。
【請求項１３】
前記尺度算出部は、前記識別妥当性尺度として、予測悪性率又は予測良性率を示す信頼区間の広さを用いる、
請求項９乃至１２のいずれか１項に記載の医用情報処理装置。
【請求項１４】
前記抽出部は、優先度に関係なく、標準で表示する特定の画像パターンを予め規定し、
前記出力部は、前記特定の画像パターンを前記表示部に表示させると共に、前記一部の画像パターンを優先して前記表示部に表示させる、
請求項９乃至１３のいずれか１項に記載の医用情報処理装置。
【請求項１５】
前記疾患を有する確率は、腫瘤の良性又は悪性の確率である、
請求項９乃至１４のいずれか１項に記載の医用情報処理装置。
【請求項１６】
判定対象の画像データが疾患を有する確率と、その確率の信頼度を示す指標である信頼区間とを算出する処理を行い、
前記確率と、前記確率の信頼区間とに基づく、疾患に関する情報を出力する、
医用情報処理方法。
【請求項１７】
コンピュータに、
判定対象の画像データが疾患を有する確率と、その確率の信頼度を示す指標である信頼区間とを算出する処理を行う機能と、
前記確率と、前記確率の信頼区間とに基づく、疾患に関する情報を出力する機能と、
を実現させる医用情報処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本明細書及び図面に開示の実施形態は、医用情報処理装置、医用情報処理方法、及び医用情報処理プログラムに関する。
【背景技術】
【０００２】
放射線画像の撮影を目的としたＸ線ＣＴ（Computed Tomography）装置やＭＲＩ（Magnetic Resonance Imaging）装置や、病理標本を取り込むバーチャルスライドスキャナ等の医療機器の発達により、デジタル化された高精細医用画像が大量に取得可能となっている。近年では、医用画像を保存するＰＡＣＳ（Picture Archiving and Communication Systems）の普及が進み、これらの画像が医師の診断所見とともに順次格納されている。このような、医用画像のデジタル化、症例の蓄積が進むにつれ、コンピュータを用いて病変検出や自動診断を行うＣＡＤ（Computer Aided Diagnosis）の重要性が増している。
【０００３】
コンピュータで病変検出や自動診断を行う場合、過去の症例を教師とした機械学習を用いるのが一般的である。医用情報処理における機械学習の問題点として、学習に用いる症例の確保の難しさが挙げられる。このような問題を解決するための手段として、知識の動的な更新が考えられる。
【０００４】
医用情報処理装置は、知識を動的に更新するために、オンライン学習によって、知識の動的な更新を実現する。また、正解率、ＳＶＭ（Support Vector Machine）のサポートベクトルの数を用いて学習の終了判定を行うことができる。
【０００５】
また、医用情報処理装置は、学習用の症例データから特徴生成パラメータを生成し、特徴生成パラメータを用いて生成された各画像データの特徴ベクトルを基にして識別パラメータを生成する。そして、医用情報処理装置は、特徴生成パラメータおよび識別パラメータを新たな診断知識として、診断知識データベースに格納する技術が知られている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１５−１１６３１９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、効率的な診断を行うための適切な情報を操作者に提示することである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成に対応する課題を他の課題として位置づけることもできる。
【課題を解決するための手段】
【０００８】
実施形態に係る医用情報処理装置は、処理部と、出力部とを備える。処理部は、判定対象の画像データが疾患を有する確率と、その確率の信頼度を示す指標である信頼区間とを算出する処理を行う。出力部は、確率と、確率の信頼区間とに基づく、疾患に関する情報を出力する。
【図面の簡単な説明】
【０００９】
【図１】図１は、第１及び第２の実施形態に係る医用情報処理装置の構成例を示す概略図。
【図２】図２は、第１の実施形態に係る医用情報処理装置の機能を示すブロック図。
【図３】図３は、第１の実施形態に係る医用情報処理装置の機能の説明を補助するための図。
【図４】図４は、第１の実施形態に係る医用情報処理装置の機能の概要を示す図。
【図５】図５は、第１の実施形態に係る医用情報処理装置の設計段階における動作をフローチャートとして示す図。
【図６】図６は、第１の実施形態に係る医用情報処理装置の識別参照情報の実装段階における動作をフローチャートとして示す図。
【図７】図７は、第１の実施形態に係る医用情報処理装置において、パターン設定表の一例を示す図。
【図８】図８は、第１の実施形態に係る医用情報処理装置において、擬似コードの一例を示す図。
【図９】図９は、第１の実施形態に係る医用情報処理装置において、表示フラグの算出過程を示す図。
【図１０】図１０は、第１の実施形態に係る医用情報処理装置において、抽出された一部のパターンの表示例を示す図。
【図１１】図１１は、第１の実施形態に係る医用情報処理装置において、パターン設定表の一例を示す図。
【図１２】図１２は、第２の実施形態に係る医用情報処理装置の機能を示すブロック図。
【図１３】図１３は、第２の実施形態に係る医用情報処理装置において、第１の機能及び第２の機能を示すデータフロー図。
【図１４Ａ】図１４Ａは、第２の実施形態に係る医用情報処理装置において、信頼区間の算出方法を説明するための図。
【図１４Ｂ】図１４Ｂは、第２の実施形態に係る医用情報処理装置において、信頼区間の算出方法を説明するための図。
【図１５】図１５は、第２の実施形態に係る医用情報処理装置において、信頼区間の算出方法を説明するための図。
【図１６】図１６は、第２の実施形態に係る医用情報処理装置において、第１Ａの機能及び第２Ａの機能を示すデータフロー図。
【図１７】図１７は、第２の実施形態に係る医用情報処理装置において、第１Ｂの機能及び第２Ｂの機能を示すデータフロー図。
【図１８】図１８は、第２の実施形態に係る医用情報処理装置において、信頼区間の算出方法を説明するための図。
【図１９】図１９は、第２の実施形態に係る医用情報処理装置において、第１Ｃの機能及び第２Ｃの機能を示すデータフロー図。
【図２０Ａ】図２０Ａは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２０Ｂ】図２０Ｂは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２０Ｃ】図２０Ｃは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２０Ｄ】図２０Ｄは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２１Ａ】図２１Ａは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２１Ｂ】図２１Ｂは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２１Ｃ】図２１Ｃは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２２】図２２は、第２の実施形態に係る医用情報処理装置において、第１Ｄの機能及び第２Ｄの機能を示すデータフロー図。
【図２３】図２３は、第２の実施形態に係る医用情報処理装置において、第１Ｅの機能及び第２Ｅの機能を示すデータフロー図。
【図２４Ａ】図２４Ａは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２４Ｂ】図２４Ｂは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２５】図２５は、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図２６Ａ】図２６Ａは、第２の実施形態に係る医用情報処理装置において、判定方法を説明するための図。
【図２６Ｂ】図２６Ｂは、第２の実施形態に係る医用情報処理装置において、判定方法を説明するための図。
【図２７】図２７は、第２の実施形態に係る医用情報処理装置において、判定方法を説明するための図。
【図２８】図２８は、第２の実施形態に係る医用情報処理装置において、判定方法を説明するための図。
【図２９】図２９は、第２の実施形態に係る医用情報処理装置において、第１Ｆの機能及び第２Ｆの機能を示すデータフロー図。
【図３０】図３０は、第２の実施形態に係る医用情報処理装置において、第１Ｇの機能及び第２Ｇの機能を示すデータフロー図。
【図３１】図３１は、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図３２】図３２は、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図３３】図３３は、敵対的生成ネットワークの構成を示す図。
【図３４】図３４は、第２の実施形態に係る医用情報処理装置において、信頼区間の算出方法を説明するための図。
【図３５Ａ】図３５Ａは、第２の実施形態に係る医用情報処理装置において、信頼区間の補正方法を説明するための図。
【図３５Ｂ】図３５Ｂは、第２の実施形態に係る医用情報処理装置において、信頼区間の補正方法を説明するための図。
【図３６】図３６は、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図３７】図３７は、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図３８】図３８は、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図３９Ａ】図３９Ａは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図３９Ｂ】図３９Ｂは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４０Ａ】図４０Ａは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４０Ｂ】図４０Ｂは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４１Ａ】図４１Ａは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４１Ｂ】図４１Ｂは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４１Ｃ】図４１Ｃは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４１Ｄ】図４１Ｄは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４１Ｅ】図４１Ｅは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４１Ｆ】図４１Ｆは、第２の実施形態に係る医用情報処理装置において、判定結果の例を示す図。
【図４２】図４２は、第２の実施形態に係る医用情報処理装置において、ユーザインターフェースの例を示す図。
【図４３】図４３は、第２の実施形態に係る医用情報処理装置において、ユーザインターフェースの例を示す図。
【図４４】図４４は、第２の実施形態に係る医用情報処理装置において、ユーザインターフェースの例を示す図。
【発明を実施するための形態】
【００１０】
以下、図面を参照しながら、医用情報処理装置、医用情報処理方法、及び医用情報処理プログラムの実施形態について詳細に説明する。
【００１１】
本実施形態は、判定対象である被検体の画像データが疾患を有する確率（以下、「疾患確率」という）と、疾患確率の信頼度を示す指標である信頼区間とに基づく、疾患に関する情報を出力することを目的とする。疾患に関する情報とは、疾患の有無、種別、予後、治療適正等を意味する。複数の画像パターン（以下、単に「パターン」という）のうち、「陰性的中率／陽性的中率」の信頼区間が比較的狭く信頼できる１又は複数の画像パターンに絞って疾患確率と信頼区間とを算出する場合について、第１の実施形態で説明する。陽性的中率（ＰＰＶ：Positive Predictive Value）とは、陽性と判定された例が実際に陽性であった比率を意味する。陰性的中率（ＮＰＶ：Negative Predictive Value）とは、陰性と判定された例が実際に陰性であった比率を意味する。また、判定対象である被検体に関する画像データの特徴量により、疾患に関する情報の確率の信頼区間を算出し、それに基づき疾患確率を算出する場合について、第２の実施形態で説明する。
【００１２】
（第１の実施形態）
図１は、第１の実施形態に係る医用情報処理装置の構成例を示す概略図である。
【００１３】
図１は、第１の実施形態に係る医用情報処理装置１を示す。医用情報処理装置１は、医用画像管理装置（画像サーバ）や、ワークステーションや、読影端末等であり、ネットワークを介して接続された医用画像システム上に設けられる。なお、医用情報処理装置１は、オフラインの装置であってもよい。
【００１４】
医用情報処理装置１は、処理回路１１と、記憶回路１２と、入力インターフェース１３と、ディスプレイ１４と、ネットワークインターフェース１５とを備える。
【００１５】
処理回路１１は、専用又は汎用のＣＰＵ（Central Processing Unit）又はＭＰＵ（Micro Processor Unit）の他、特定用途向け集積回路（ＡＳＩＣ：Application Specific Integrated Circuit）、及び、プログラマブル論理デバイス等の処理回路を意味する。プログラマブル論理デバイスとしては、例えば、単純プログラマブル論理デバイス（ＳＰＬＤ：Simple Programmable Logic Device）、複合プログラマブル論理デバイス（ＣＰＬＤ：Complex Programmable Logic Device）、及び、フィールドプログラマブルゲートアレイ（ＦＰＧＡ：Field Programmable Gate Array）等の回路が挙げられる。処理回路１１は、記憶回路１２に記憶された、又は、処理回路１１内に直接組み込まれたプログラムを読み出し実行することで後述する機能を実現する。
【００１６】
また、処理回路１１は、単一の回路によって構成されてもよいし、複数の独立した回路要素の組み合わせによって構成されてもよい。後者の場合、複数の記憶回路１２が複数の回路要素の機能に対応するプログラムをそれぞれ記憶するものであってもよいし、１個の記憶回路１２が複数の回路要素の機能に対応するプログラムを記憶するものであってもよい。なお、処理回路１１は、制御部の一例である。
【００１７】
記憶回路１２は、ＲＡＭ（Random Access Memory）、フラッシュメモリ（Flash Memory）等の半導体メモリ素子、ハードディスク、及び光ディスク等によって構成される。記憶回路１２は、ＵＳＢ（Universal Serial Bus）メモリ及びＤＶＤ（Digital Video Disk）等の可搬型メディアによって構成されてもよい。記憶回路１２は、処理回路１１において用いられる各種処理プログラム（アプリケーションプログラムの他、ＯＳ（Operating System）等も含まれる）や、プログラムの実行に必要なデータ等を記憶する。また、ＯＳに、読影者等の操作者に対するディスプレイ１４への情報の表示にグラフィックを多用し、基礎的な操作を入力インターフェース１３によって行うことができるＧＵＩ（Graphical User Interface）を含めることもできる。
【００１８】
また、記憶回路１２は、Ｘ線ＣＴ（Computed Tomography）装置及びＭＲＩ（Magnetic Resonance Imaging）装置等の医用画像診断装置で生成された撮像画像を記憶する。なお、記憶回路１２は、記憶部の一例である。
【００１９】
入力インターフェース１３は、操作者によって操作が可能な入力デバイスと、入力デバイスからの信号を入力する入力回路とを含む。入力デバイスは、トラックボール、スイッチ、マウス、キーボード、走査面に触れることで入力操作を行うタッチパッド、表示画面とタッチパッドとが一化されたタッチスクリーン、光学センサを用いた非接触入力デバイス、及び音声入力デバイス等によって実現される。操作者により入力デバイスが操作されると、入力回路はその操作に応じた信号を生成して処理回路１１に出力する。なお、医用情報処理装置１は、入力デバイスがディスプレイ１４と一体に構成されたタッチパネルを備えてもよい。なお、入力インターフェース１３は、入力部の一例である。
【００２０】
ディスプレイ１４は、液晶ディスプレイパネル、プラズマディスプレイパネル、及び有機ＥＬ（Electro Luminescence）パネル等の表示デバイスである。ディスプレイ１４は、処理回路１１の制御に従って生成された医用画像を表示する。なお、ディスプレイ１４は、表示部の一例である。
【００２１】
ネットワークインターフェース１５は、パラレル接続仕様やシリアル接続仕様に合わせたコネクタによって構成される。ネットワークインターフェース１５は、各規格に応じた通信制御を行い、電話回線を通じてネットワークに接続することができる機能を有しており、これにより、医用情報処理装置１をネットワークに接続させることができる。なお、ネットワークインターフェース１５は、通信部の一例である。
【００２２】
図２は、医用情報処理装置１の機能を示すブロック図である。
【００２３】
処理回路１１が記憶回路１２等の非一過性の記録媒体に記憶されたコンピュータプログラムを実行することによって、医用情報処理装置１は、処理機能２１と、抽出機能２２と、出力機能２３とを実現する。なお、機能２１〜２３の全部又は一部は、医用情報処理装置１にＡＳＩＣ等の回路として実現されるものであってもよい。
【００２４】
処理機能２１は、判定対象である被検体の画像データを複数のパターンに分類して、各パターンの疾患確率を示す指標値を識別妥当性尺度として算出する機能を含む。例えば、処理機能２１は、被検体の画像データに基づいて識別妥当性尺度を生成する学習済みモデルに対して、被検体の画像データを入力することで、識別妥当性尺度を生成する。また、例えば、疾患確率は、腫瘤の良性／悪性の確率である。なお、処理機能２１は、処理部の一例である。
【００２５】
抽出機能２２は、複数のパターンの各パターンの識別妥当性尺度に基づいて、複数のパターンから、疾患の識別が行いやすい一部のパターンを抽出する機能を含む。なお、抽出機能２２は、抽出部の一例である。
【００２６】
出力機能２３は、抽出機能２２によって抽出された一部のパターンを優先して出力する機能を含む。例えば、出力機能２３は、抽出機能２２によって抽出された一部のパターンを優先した表示画面を生成し、その表示画面をディスプレイ１４に表示させる機能を含む。なお、出力機能２３は、出力部の一例である。
【００２７】
ここで、処理機能２１は、関心領域設定機能２１１と、特徴算出機能２１２と、尺度取得機能２１３とを含む。
【００２８】
関心領域設定機能２１１は、設計段階において、複数症例を含む画像データセットの各画像データから関心領域を抽出する機能を含む。また、関心領域設定機能２１１は、実装段階において、新規症例の画像データから関心領域を抽出する機能を含む。つまり、関心領域設定機能２１１が被検体の画像データの関心領域を設定することで、関心領域に基づいて、被検体の画像データを複数のパターンに分類することが可能になる。なお、関心領域設定機能２１１は、関心領域設定部の一例である。
【００２９】
特徴算出機能２１２は、関心領域設定機能２１１によって抽出された設計段階における各画像データの関心領域に基づいて特徴量を特徴ベクトルとして算出し、当該特徴ベクトルに基づいて分類尺度ベクトルを算出する機能を含む。また、特徴算出機能２１２は、関心領域設定機能２１１によって抽出された実装段階における画像データの関心領域に基づいて特徴量を特徴ベクトルとして算出し、当該特徴ベクトルに基づいて分類尺度ベクトルを算出する機能を含む。なお、特徴算出機能２１２は、特徴算出部の一例である。
【００３０】
尺度取得機能２１３は、実装段階における各群の分類尺度ベクトルの分布情報と、新規症例の疾患確率（不明なら均等などと仮定）とを取得する機能を含む。尺度取得機能２１３は、特徴算出機能２１２から、群（クラス）に属する分類尺度ベクトルを直接取得する。本明細書及び図面において、疾患あり、疾患なし等の、利用する複数の画像データの属性を「群」と称する。尺度取得機能２１３は、特徴算出機能２１２の各群の分類尺度ベクトルの統計的分布のパラメータを取得する。なお、尺度取得機能２１３は、尺度取得部の一例である。
【００３１】
抽出機能２２は、尺度取得機能２１３によって取得された分類尺度ベクトルから各パターンに分類する分類基準とパターンの表示条件を取得し、パターン候補の各々について、各群の分類尺度ベクトルの分布（と関心領域の分類尺度ベクトル）から予測疾患確率等の参考情報を算出する機能を含む。抽出機能２２は、パターン候補の各々について、そのパターンで識別する妥当性の尺度を識別妥当性尺度として算出し、それに基づいて表示するパターンの組み合わせを決定する機能を含む。抽出機能２２は、識別ベクトルの値、パターン分類基準、群尺度により表示フラグを設定する（パターン分類基準は、より一般にはパターン重み関数）機能を含む。
【００３２】
出力機能２３は、抽出機能２２によって抽出され設定された表示フラグに基づいて、パターンの名称と、それぞれのパターンの参考情報（群尺度情報）をディスプレイ１４に表示させる機能を有する。
【００３３】
比較例において、識別目的と対象画像データとの関係を参考情報として提示している。例として病名が識別目的であり、その病名の典型画像に対して類似症例を提示して、対象画像と近いかどうかを目視で判断するための情報を提供する。また、対象画像がその病名に対し典型的な画像かを示す典型度を提示する。これらは診断に当たって参考となる情報である。しかし、判断が妥当であるかどうかを読影者が追確認する状況、あるいは、判断の過程や追確認の過程を読影者の判断の記録として言語により記載する状況（カルテへの記載や診断レポートの作成）を考えると、従来例が提示する参考情報は、読影者が追確認や判断過程の記録に有用な情報とはいえない。
【００３４】
追確認や判断過程の記録のためには、識別目的と異なる観点での参考情報の提示が必要になる。例えば、風邪かどうかを判断するため、また、判断過程を記録するためには、咳などの症状の有無や体温など、風邪かどうか以外の情報の提示が重要である。また、提示される参考情報は、類似画像のようなものではなく言語として表現できるものである必要がある。さらには、画像診断の目的に供するためには、提示される参考情報は症状の有無や体温といった直接計測できるものではなく、画像データの特徴を表す情報である必要がある。
【００３５】
そこで、医用情報処理装置１は、処理機能２１と、抽出機能２２と、出力機能２３とを備えるものとする。
【００３６】
なお、機能２１〜２３の具体的な動作については、図３〜図１１を用いて後述する。
【００３７】
図３は、医用情報処理装置１の機能の説明を補助するための図である。
【００３８】
医用情報処理装置１は、画像データの疾患領域を設定し、提示される情報と画像データの観察結果とに基づいて、読影者が疾患領域の疾患の状態、例えば、「陰性」又は「陽性」を識別する例について述べる。提示される情報は、読影者が「陰性」又は「陽性」を識別するための参考となる識別参考情報である。
【００３９】
識別方法を構成するための識別対象のデータとして、例えば、ＮＬ個の乳腺造影ＭＲＩ画像データが使用される。ＮＬ個の乳腺造影ＭＲＩ画像データのそれぞれは、確定診断の結果に基づき「良性」又は「悪性」が教師データとして与えられる。悪性が陽性に、良性が陰性にそれぞれ対応する。
【００４０】
医用情報処理装置１は、これらＮＬ個の症例からなるデータセットに基づいて分類尺度ベクトルＮｘを算出することで、パターン名称等からなるパターンを生成する。そして、ＮＬ個の症例からなるデータセットの識別結果を、識別対象の画像データに対して適用し、識別参考情報として読影者に提示する。読影者は、提示された識別参考情報から容易に「陰性」又は「陽性」の識別を行うことができる。さらに、読影者は、識別参考情報に基づいて、識別対象の画像データ等の画像データを観察して、その識別が妥当であるか確認し、必要なら識別を修正することもできる。
【００４１】
医用情報処理装置１は、識別参考情報として、複数のパターンを表示する。医用情報処理装置１は、パターンとして、画像データの特徴を表す名称や、識別対象の画像データがそのパターンに属するかどうかの参考情報や、そのパターンにおいて識別対象の画像データは良性か悪性かを判断するための参考情報を表示する。パターンは常に同じものを表示するのではなく、必要なパターンのみが表示される。パターンの表示または非表示については、そのパターンでの識別が妥当性であるかを表す値（識別妥当性尺度）に基づいて決定される。
【００４２】
このような識別方法を設計するにあたり、陽性的中率（ＰＰＶ）が「ＰＰＶ＜定数」となる識別能力が得られることを設計の条件として導入する。さらに「特異度が最大」となる識別能力が期待できるよう設計の目標を設定する。例えば、定数は「０．０５」であり、これは、「ＰＰＶ＜０．０５」となる識別能力を意味する。この条件により、検査で関心領域があるパターンに分類されると判断した場合、関心領域が悪性である確率が５％以下であること、すなわち、良性である可能性が高いことが要求される。
【００４３】
さらに、特異度が最大であることを目標としているため、良性の関心領域が悪性と判断される可能性が小さくなることが期待される。条件や目標は感度、陰性的中率（ＮＰＶ）等の他の基準を用いて設定することも可能である。
【００４４】
実施形態に係る医用情報処理装置は、次の処理の２段階（１），（２）からなる。
（１）設計段階（データセットの識別）
（２）実装段階（新規症例での処理）
【００４５】
上記（１）の設計段階は、教師データセットの確率分布に対する解析を行い、その結果を保存するものである。当該設計段階は、識別の対象となる新規症例が与えられる前に１度だけ実施すればよい。
【００４６】
上記（２）の実装段階は、新規症例ごとに実行する実装段階の処理である。この処理では新規症例の１つについて識別の参考となる情報を算出し提示するものとして説明する。
【００４７】
なお、上記（１），（２）の各段階における処理をどこで分割するかは任意であるが、以下、代表例に基づいてそれらの処理を記載する。
【００４８】
図４は、医用情報処理装置１の機能の概要を示す図である。
【００４９】
図４に示すように、処理機能２１は、被検体の画像データに基づいて識別妥当性尺度を生成する学習済みモデルに対して、被検体の画像データを入力することで、識別妥当性尺度を生成する。抽出機能２２は、複数のパターンの各パターンの識別妥当性尺度に基づいて、複数のパターンから、疾患の識別が行いやすい一部のパターンを抽出する。出力機能２３は、抽出機能２２によって抽出された一部のパターンを優先して出力する（図４の下側の表示画面）。
【００５０】
図５は、医用情報処理装置１の設計段階における動作をフローチャートとして示す図である。図６は、医用情報処理装置１の識別参照情報の実装段階における動作をフローチャートとして示す図である。図５及び図６において、「ＳＴ」に数字を付した符号はフローチャートの各ステップを示す。
【００５１】
処理機能２１は、複数症例を含む画像データセットを取得する（ステップＳＴ１）。関心領域設定機能２１１は、画像データセットの各画像データについて関心領域を設定する（ステップＳＴ２）。ステップＳＴ２において、例えば、関心領域設定機能２１１は、乳がんの造影ＭＲＩにおいて収集された乳腺造影ＭＲＩ画像データから、解析対象領域、例えば、腫瘤領域を抽出して関心領域として設定する。関心領域設定機能２１１は、複数症例を含むＮＬ個の乳腺造影ＭＲＩ画像データから、合計ＮＬ個の関心領域を設定する。
【００５２】
特徴算出機能２１２は、ステップＳＴ２によって設定されたＮＬ個の関心領域に対して、教師データをそれぞれ設定する（ステップＳＴ３）。具体例としては、各関心領域に対して、腫瘤の状態、例えば、「良性」又は「悪性」が割り当てられる。又は、各関心領域に対して、腫瘤の状態、例えば、陰性、陽性、正常、疾患、異常等を割り当てられてもよい。また、具体的な疾患名を用いて３つ以上の値が割り当てられてもよい。教師データとして確定診断の結果を用いる場合、教師データは確定診断が得られたときに初めて設定可能となる。
【００５３】
特徴算出機能２１２は、ステップＳＴ２によって設定された複数の関心領域にそれぞれ対応する複数の特徴量を特徴ベクトルとして算出する（ステップＳＴ４）。特徴算出機能２１２は、ステップＳＴ４において、関心領域ｉに対するｊ番目の特徴量をｆｉｊ（図３に図示）と表現する。また、関心領域ｉに対する複数の特徴量は、特徴ベクトルＦｉを構成する。
【００５４】
例えば、特徴量は、造影ダイナミック画像の初期造影強度の平均値である。初期造影強度は、造影後の第１相、又は、第２相の画素値から造影前の画素値を引いた画素値の差を造影前画素値で割った値を関心領域の全画素で平均をとった値である。その他、特徴量は、ピーク造影強度や後期造影強度等の様々な造影関連パラメータであってもよく、これらの平均であってもよい。さらに、特徴量は、平均値のみでなく各種の統計量であってもよい。造影関連パラメータとテキスチャ解析値の組み合わせで多種多様な特徴量が算出され得る。また、円形度等の形態的な画像の特徴量でもよい。通常、造影ダイナミック画像には位置合わせ処理や、動き補正処理が適用される。
【００５５】
尺度取得機能２１３は、分類尺度ベクトルの算出方法を決定し、その算出方法により分類尺度ベクトルを算出する（ステップＳＴ５）。尺度取得機能２１３は、教師データに対応する各画像データの特徴ベクトルを変換することで、分類尺度ベクトルを算出する。これは、例えば、特開２０１５−１１６３１９号公報において識別ベクトルの算出方法を決めることに相当する。関心領域ｉに関するｋ番目の分類尺度をｘｉｊ（図３に図示）と表現する。関心領域ｉに関する複数の分類尺度を分類尺度ベクトルＸｉで表す。特徴ベクトルＦｘは、分類尺度ベクトルＸｉとの関係をモデルにより次の式（１）で表される。
【００５６】
【数１】
2020149684.tif 000003
【００５７】
モデルとモデルパラメータを適切にとれば、分類尺度ベクトルの各々の要素が表す特徴に言語的な意味づけを与えることができる。そこで、尺度取得機能２１３は、ステップＳＴ５において、各々の要素が適切に意味づけできるようにモデルパラメータを決定する。
【００５８】
分類尺度の各々の要素（以降、要素番号をｋで表す）は、その値の大小に応じて名称が与えられる。例えば、分類尺度ベクトルｘｉの最初の要素（ｋ＝１）は辺縁や形状の複雑性を表す尺度である。この尺度が小さな値（ｘｉｊ≦Ｔ）には「ｓｍｏｏｔｈ」や「ｓｉｍｐｌｅ」、大きな値（ｘｉｊ＞Ｔ）には、「ｉｒｒｅｇｕｌａｒ」や「ｃｏｍｐｌｅｘ」といった名称が与えられる。分類尺度ベクトルｘｉの２番目の要素（ｋ＝２）は、造影曲線の後半が下降するとき小さな値（名称の例は、「ｗａｓｈｏｕｔ」、「ｄｏｗｎｓｌｏｐｅ」等）、後半が下降しないとき大きな値（名称の例は、「ｎｏｔ  ｗａｓｈｏｕｔ」、「ｎｏｔ  ｄｏｗｎｓｌｏｐｅ」等」となる。
【００５９】
なお、特徴ベクトルから分類尺度ベクトルへの変換パラメータは、全領域のデータと教師データを基に、種々の線形／非線形最適化手法を用いて決定する。人工ニューラルネットワークやランダムツリー等の機械学習の技術も適用可能である（後述する第２の変形例）。比較的単純な方法は、特徴ベクトルの要素のうち適切な要素のみを選び、それを分類尺度ベクトルとするものである。具体的には、２つの要素の組み合わせのうち、２つの要素で２変数識別を行ったとき、識別の妥当性についての尺度が最大となる組み合わせを選ぶといった構成が可能である。識別妥当性尺度については後述する。
【００６０】
尺度取得機能２１３は、分類尺度ベクトルの算出方法を決定した後、その決定方法に従って分類尺度ベクトルＸｉを算出する。
【００６１】
抽出機能２２は、ステップＴＳ５によって取得された分類尺度ベクトルから各パターンに分類する（ステップＳＴ６）。最も単純なパターン分類法として閾値による分類である。
【００６２】
図７は、パターン設定表を示す図である。図７は、閾値によるパターン分類を行うものである。
【００６３】
閾値により分類されるパターン候補の一例を表として示す図である。図７に示す表において、ｈ＝１，２，…の各々が１つのパターンに対応している。図７に示す表において、複数のパターン番号ｈ＝１，２，…，９，１０，…が表現される。実際には、多数のパターンが表現されるが、当該表はその一部のみを示したものである。なお、説明の便宜のため、正規分布モデルでいうｈ１〜ｈ９の類型を、実施形態において、「パターン」と呼称する。
【００６４】
当該表のパターン分類基準は、判別不等式として表現される。当該表では、パターン分類尺度のどの要素に関する閾値であるか、不等号の向き、閾値の値により構成される。当該表では、第１パターンの分類に辺縁や形状の複雑性を表す分類尺度であるｋ＝１の要素ｘｓ１を用い、複雑性が低いときにこのパターンに属すると判定する。そのため、パターン名称として「辺縁平滑」が設定される。
【００６５】
判定に用いる要素、不等号の向き、閾値を決めるために、種々の判別分析手法を適用することができる。これらの手法を用いれば各群の分類尺度ベクトルから判別式を決定することができる。
【００６６】
閾値を決める具体的な方法の１つは、「ＰＰＶ＜定数」の条件を満たす閾値のうち「特異度が最大」の閾値、およびそのときの不等号の向きを選ぶものである。ＰＰＶは陽性適中率であり、このパターンに判定される関心領域のうち、教師データが「腫瘤」である腫瘤群の割合である。ＰＰＶの代わりに予測陽性率の信頼区間（後述）を基に閾値を決定することも可能である。これは選んだパターンが陽性と判定されることを想定した場合であるが、選んだパターンが陰性であると捉えることもでき、このとき算出されるＰＰＶは、実際にはＰＰＶではなく、１−ＮＰＶ（陰性適中率）となる。
【００６７】
閾値による分類は、不等式だけでなく、パターン重み関数によっても表現することができる。当該表のパターン分類基準の欄の下段に記載されている式がパターン重み関数による表現である。
【００６８】
上記は、全関心領域の分類尺度ベクトルと教師データを基に領域の数をカウントして閾値を決める方法の場合について記載しているが、その場合に限定されるものではない。他の方法として、各群のパターン分類尺度ベクトルを基に、各群の確率分布を、例えば多次元正規分布等で表し、分布関数（又は、確率密度関数）をフィッティングにより求める方法を用いてもよい。この方法の場合も同様に各群の分布を基に、ＰＰＶや特異度等の条件を基に閾値および不等号の向きを決めることができる（後述する第１の変形例）。
【００６９】
また、当該表の表示基準は、パターンごとに表示基準と子パターンを有する。これらの詳細については次の実装段階で説明する。
【００７０】
図５の説明に戻って、医用情報処理装置１は、以下の（ａ）〜（ｄ）のデータを、次の実装段階で利用可能なように保存する。
【００７１】
（ａ）分類尺度ベクトルの分布
分類尺度ベクトルの分布は、ＮＬ個の算出された分類尺度ベクトル、分類尺度ベクトルの分布を表すパラメータ、関心領域の疾患確率を意味する。分布情報の表現方法として、上記のほか、各群の（１又は複数の閾値）による各レベルの発生数、各群の各レベルの発生確率、各群のヒストグラム等の表現方法で保存する構成もあり、これらは同等である。
（ｂ）パターン名称
（ｃ）パターン分類基準
（ｄ）表示基準（と子パターン）
【００７２】
図６の説明に進んで、医用情報処理装置１は、新規症例に対する識別参考情報の提示処理を行うために、識別対象の画像データ（例えば、乳腺造影ＭＲＩ画像データ）を取得する（ステップＳＴ１１）。そして、関心領域設定機能２１１は、ステップＳＴ１１によって取得された識別対象の画像データから腫瘤領域を関心領域として設定する（ステップＳＴ１２）。特徴算出機能２１２は、関心領域に対応する特徴量を特徴ベクトルとして算出する（ステップＳＴ１３）。ステップＳＴ１３における特徴ベクトルの算出は、図５に示すステップＳＴ４と同等の方法で行われる。
【００７３】
特徴算出機能２１２は、識別対象の画像データについて、分類尺度ベクトルの算出方法を決定し、その算出方法により分類尺度ベクトルを算出する（ステップＳＴ１４）。ステップＳＴ１４における分類尺度ベクトルの算出方法の決定と、分類尺度ベクトルの算出とは、図５に示すステップＳＴ５と同等の方法で行われる。単純に特徴ベクトルの２要素を選ぶ方法を用いる場合は、特徴ベクトルは最初からその２要素のみを計算し、それらをそのまま分類尺度とすればよい。
【００７４】
尺度取得機能２１３は、図５に示すステップＳＴ３によって設定された複数の教師データに対応する分類尺度ベクトルを取得する（ステップＳＴ１５）。又は、群ごとに決定した分布関数のパラメータ、例えば平均と標準偏差を取得する（後述する第１の変形例）。また、新規症例における各群の発生比率（関心領域の疾患確率）を読み込む。例えば、本方法での検査対象となる症例の７割が悪性、３割が良性であることが想定される場合、これらの割合を読み込む。７割という割合は「関心領域の疾患確率」であり、教師データセットでの各群の発生比率とは一般に異なる。関心領域の疾患確率は有病率と呼ばれることもある。
【００７５】
抽出機能２２は、ステップＳＴ１５によって取得された分類尺度ベクトルから各パターンに分類する（ステップＳＴ１６）。抽出機能２２は、群ごとに取得した分類尺度ベクトルの分布（平均・標準偏差・発生確率、又は教師データセット全部の分類尺度ベクトル）をもとに、そのパターンがどの群に属するかの指標になる値である群尺度情報を算出する（ステップＳＴ１７）。群尺度情報は、パターンごとに算出される。パターンの数値的な定義は「パターン分類基準」で指定されているので、群尺度情報は、図７の表の「パターン分類基準」に従って算出する。
【００７６】
「パターン分類基準」は、より一般には「パターン重み関数」で表現される。図７の表に示す例では、「ｘｓ１≦Ｔ１ｈ」は関数ｕ（Ｔ１ｈ−ｘｓ１）で表される。ｕ（ｘ）は、次の式（２）に示すように階段関数である。
【数２】
2020149684.tif 000004
【００７７】
パターン重み関数としては、新規症例の分類尺度ベクトルＸｓが、パターンｈの典型的な点μｈに近いときは「１」に近く、遠ざかるほど「０」に近い値となる釣鐘上の関数を用いることができる。
【数３】
2020149684.tif 000005
【００７８】
ｗｈ（Ｘｓ）は、分類尺度ベクトルＸｓがパターンｈに属するかどうかを判断する関数である。
【００７９】
重み関数の別の構成では（後述する第１の変形例）、典型点μｈが新規症例の分類尺度ベクトルＸｓに等しいとしたとき（μｈ＝Ｘｓ）、任意のベクトルＸがパターンｈに属するかを判断するよう以下の関数を用いる。
【数４】
2020149684.tif 000006
【００８０】
このとき、パターンｈのパターン分類基準は「釣鐘上、ｋ＝１，２」のような表記となり、この表記に従い、ｋ＝１とｋ＝２の２つの要素以外はＲｈ−１の要素を「０」に置き換えることもできる（後述する第１の変形例）。
【００８１】
このように、群尺度情報の算出に新規症例の分類尺度ベクトル が参照される構成も可能である。閾値による判定も同様に に依存して閾値が変化する構成が可能である。
【００８２】
続いて、群尺度情報を算出する具体例（Ａ）〜（Ｇ）について説明する。なお、「Ｎｐｂ」はパターンに含まれる領域のうち良性の領域の数であり、「Ｎｐｍ」は、パターンに含まれる領域のうち悪性の領域の数であり、Ｎｐｂ＝Ｎｂ−Ｎｐｂ、Ｎｎｍ＝Ｎｍ−Ｎｐｍとする。ここで、Ｎｂは、良性群の領域数、Ｎｍは、悪性群の領域数である。Ｐｍは、関心領域の疾患確率であり、Ｐｂ＝１−Ｐｍである。
【００８３】
（Ａ）精度＝（Ｎｐｍ＋Ｎｎｍ）／（Ｎｂ＋Ｎｍ）
（Ｂ）良性群中陽性率（期待値）Ｐｐｂ＝Ｎｐｂ／（Ｎｐｂ＋Ｎｎｂ）  −（１−特異度）
（Ｃ）良性群中陽性率の信頼区間（下限と上限）−省略
（Ｄ）悪性群中陽性率（期待値）Ｐｐｍ＝Ｎｐｍ／（Ｎｐｍ＋Ｎｎｍ）  −（感度）
（Ｅ）悪性群中陽性率の信頼区間（下限と上限）−省略
（Ｆ）予測陽性率（期待値）＝Ｎｐｍ／（Ｎｐｂ＋Ｎｐｍ）  −（陽性的中率ＰＰＶ）、又は、予測陽性率(期待値)＝Ｐｐｍ×Ｐｍ／（Ｐｐｂ×Ｐｂ＋Ｐｐｍ×Ｐｍ）
（Ｇ）予測陽性率の信頼区間（下限と上限）  −省略
【００８４】
なお、群ごとの領域の数が取得されるのではなく、群ごとに決定した密度関数Ｐｐｂ（Ｘ）、Ｐｐｍ（Ｘ）、疾患確率Ｐｂが取得されるとき（後述する第１の変形例）、ＰｐｂとＰｐｍは重みつき積分で求められる。
【００８５】
【数５】
2020149684.tif 000007
【００８６】
予測陽性率は、良性／悪性の事前確率Ｐｂ、Ｐｍ＝１−Ｐｂ、データの確率分布Ｐｐｂ（Ｘ）、Ｐｐｍ（Ｘ）、測定ノイズの確率分布Ａｗｈ（Ｘ,Ｘｓ）が与えられ（Ａ：積分を「１」にするための係数）、データＸｓが測定されたとき、データが悪性領域から得られたものである確率をベイズの定理により事後確率として求めたものである。
【００８７】
良性群中陽性率、悪性群中陽性率、予測陽性率の信頼区間を算出するさまざまな方法が提案されている。例えば、二項分布に基づく方法が正確であることが知られており、これを採用することが望ましい。
【００８８】
上記（Ａ）〜（Ｇ）の他、オッズ比などさまざまな指標値が考えられ、これらの期待値と信頼区間が群尺度情報として算出される。
【００８９】
図６の説明に戻って、抽出機能２２は、識別妥当性尺度を算出する（ステップＳＴ１８）。抽出機能２２は、そのパターンに従って陽性または陰性と判断することをどれだけ信頼できるかを表す尺度として群尺度の算出値を用いて識別妥当性尺度を算出する。精度や予測陽性率など群尺度として算出された期待値または信頼区間を直接に識別妥当性尺度として用いることもできるが、実施形態では、予測陽性率の信頼区間などから以下の擬似コードで識別妥当性尺度を求める。擬似コード内の「陽性判定向け」、「陽性判定時ＰＰＶ設計下限値」、「予測陽性率の信頼区間上限値」は、図７に示す表の「設計目標」の記載内容である。
【００９０】
図８は、擬似コードの一例を示す図である。
【００９１】
図８に示す擬似コードの１行目と２行目は設計範囲からの信頼区間のはみ出し量となっている。３行目を満たすときは陽性と判定することが妥当であり、この場合、１行目の値を識別妥当性尺度として使用する。ただし、値が大きい方の妥当性が高くなるように正負を反転している。
【００９２】
５行目は信頼区間が設計条件に収まっていることを示しており、この場合は感度の信頼区間下限を識別妥当性尺度として使用する。３行目を満たさないときは陰性と判断することが妥当なため２行目の値を正負反転して識別妥当性尺度として使用する。１０行目は信頼区間が設計条件に収まっていることを示しており、この場合は特異度の信頼区間下限を識別妥当性尺度として使用する。
【００９３】
図６の説明に戻って、抽出機能２２は、抽出された一部のパターンに対し表示フラグを設定する（ステップＳＴ１９）。表示フラグは、図７に示す表の表示基準に従って算出する。ｈ＝９のパターンは表示基準が「すべての上位パターン（親パターン）より識別妥当性尺度が大きい」となっている。この条件を満たすときだけそのパターンの表示フラグを「ｔｒｕｅ」に設定し、そうでないときは表示フラグを「ｆａｌｓｅ」に設定する。さらに追加の条件としてパターン重み関数が「０」でないという条件を付加する。このためパターン分類基準を満たさないパターンは表示フラグが「ｆａｌｓｅ」となり、さらにそれ以外で表示基準を満たす場合は表示フラグが「ｔｒｕｅ」、満たさない場合は表示フラグが「ｆａｌｓｅ」に設定される。
【００９４】
図７に示す表で、パターン１，４，９の３つのパターンがパターン分類基準を満たす場合を考える。表示フラグを「ｔｒｕｅ」にする候補はこの３つである。ｈ＝９のパターンの上位パターンは、子パターンに「９」が記載されているｈ＝４のパターン（パターン９の親パターン）と、子パターンに4が記載されているｈ＝１のパターン（パターン４の親パターン）の２つである。
【００９５】
図７に示す表では、これらの予測陽性率上限はそれぞれ、３４．４％（パターン１）、６．５％（パターン４）、１６．７％（パターン９）と表示されており、識別妥当性尺度は、−０．３４４％（パターン１）、−０．０６５％（パターン４）、−０．１６７％（パターン９）である。パターン９の識別妥当性尺度は、パターン４の識別妥当性尺度よりも小さいため、パターン９は表示基準を満たしていないことがわかる。そのため、パターン９の表示フラグは「ｆａｌｓｅ」に設定され、そのため、パターン９は実際には画面には表示されない。パターン４は、表示基準は満たすため「ｔｒｕｅ」となり、パターン１は常に「ｔｒｕｅ」という表示基準のため、表示フラグは「ｔｒｕｅ」となる。表示フラグの算出過程について、図９に示す。図９において、図７に示すパターン９の識別妥当性尺度は表示基準を満たしていないので、パターン９の表示フラグは「ｆａｌｓｅ」に設定される。一方で、パターン４の識別妥当性尺度は表示基準を満たすので、パターン４の表示フラグは「ｔｒｕｅ」に設定される。
【００９６】
医用情報処理装置１は、抽出された一部のパターンに関するパターン名称、群尺度の算出値を図７に示す表に従ってディスプレイ１４に表示させる（ステップＳＴ２０）。
【００９７】
図１０は、抽出された一部のパターンの表示例を示す図である。
【００９８】
図１０に示すように、最初の「辺縁平滑パターン」は、図７に示す表の先頭パターンのパターン名称を基に文字情報として表示されている。各パターンに表示している数値は群尺度の算出値である。「良性群中」の数値は良性群中陽性率（期待値）、「悪性群中」の数値は悪性群中陽性率（期待値）の算出値である。また、予測陽性率（又は、予測悪性率）の期待値や信頼区間の広さも識別妥当性尺度として表示される。
【００９９】
パターンは表示フラグが「ｔｒｕｅ」のものだけが表示される。従って、図１０の３番目のパターン（パターン９）は実際には表示されない。なお、抽出機能２２が、優先度に関係なく、標準で表示する特定のパターンを予め規定することで、出力機能２３は、抽出された一部のパターンを優先してディスプレイ１４に表示させる一方で、優先度に関係なく特定のパターンをディスプレイ１４に表示させることもできる。
【０１００】
このように、医用情報処理装置１によれば、識別妥当性尺度の値に応じて、各パターンを表示するかどうかが判断された結果、一部のパターンのみが表示されるので、効率的な診断を行うための適切な情報を操作者に提示することができる。
【０１０１】
なお、学習データ群を予め複数パターンに分類して、各パターンの良性確率／悪性確率（例えば、「良性＝２〜５０％」のように分布をもつ場合もある）を予め算出しておいてもよい。その場合、被検体の画像データに該当する複数のパターンを設定した時に、出力機能２３は、複数のパターンのうち、悪性、又は、良性の識別を行いやすい一部パターンを優先してディスプレイ１４に表示させることもできる。
【０１０２】
医用情報処理装置１によれば、目的とする判断と異なる観点の情報を、画像データから得られる複数のパターンとして表示することができるし、パターンの参考情報として、そのパターンを示す文字情報（具体的にはパターンの名称）を表示することもできる。一方で、医用情報処理装置１によれば、複数のパターンを網羅的に表示したいが、反面、提示するパターンが過多にならないようにすることができる。
【０１０３】
医用情報処理装置１によれば、特に、パターン候補ごとに参考情報として画像の特徴をあらわす名称を表示することで、乳がんＭＲＩの読影とその結果のレポート記載に有用な情報を提供するができる。その際、表示するパターンに不足があると正しい診断のために必要な情報を提示できず、過多であると効率的な診断の妨げになるが、医用情報処理装置１によれば、必要最小限のパターンの提示を実現できる。
【０１０４】
（第１の変形例）
図５に示す設計段階において、各群のパターン分類ベクトルの確率分布を多次元正規分布で表し、確率密度関数をフィッティングにより求める。これらの値は分類尺度ベクトルの分布情報として保存される。フィッティングの結果、各群の平均値と標準偏差が求められる。
【０１０５】
図１１は、パターン設定表の一例を示す図である。図１１は、正規分布モデルによる分類尺度ベクトルの分布の近似によりパターン分類を行うものである。
【０１０６】
パターン重み関数は、前述したｗｈ（Ｘ，Ｘｓ）の式で算出する。Ｒｈ−１の要素のうち、表の「パターン重み関数」に記載のない要素は、値を「０」に設定する。例えば、図１１に示す表のパターン６のパターン重み関数にはｘｓ１，ｘｓ２，ｘｓ３の３要素が記載されている。これは分類尺度ベクトルのこれら３要素についての正規分布用の重み関数を用いることを示しており、これら３要素以外のＲｈ−１の要素は「０」に置き換える。
【０１０７】
図６に示す実装段階において、保存された平均値と標準偏差を分類尺度ベクトルの分布として取得し、パターン重み関数は、設計段階と同様に算出する。良性群中陽性率（期待値）Ｐｐｂと悪性群中陽性率（期待値）Ｐｐｍとは、上記式（５）により算出される。
【０１０８】
なお、予測陽性率（期待値）は、ＰｐｂとＰｐｍとから、前述のように算出される。また、信頼区間についても前述のように算出される。
【０１０９】
医用情報処理装置１の第１の変形例によれば、識別妥当性尺度の値に応じて、各パターンを表示するかどうかが判断された結果一部のパターンのみが表示されるので、同様に、効率的な診断を行うための適切な情報を操作者に提示することができる。
【０１１０】
（第２の変形例）
前述の図３は、特徴ベクトルから分類尺度ベクトルへの変換のパラメータを決定するための構成を示す。分類尺度ベクトルへの変換のモデルパラメータを決めたとき、教師データでの識別成績を求めるステップを簡潔にすることもできる。関心領域ｉごとに次の式（６）からＸｉを算出し、パターンｈごとに次の式（７）からＰｐｂを算出する。
【数６】
2020149684.tif 000008
【数７】
2020149684.tif 000009
【０１１１】
ここで、パターンｈの識別妥当性尺度は、擬似コード（図８に図示）によりＰｐｂ、Ｐｐｍから算出する。関心領域ｉの最大識別妥当性尺度は、識別妥当性尺度（ｈ＝１，２，…）の最大値である。識別成績は、最大識別妥当性尺度（ｉ＝１，２，…）の信頼区間下限値である。
【０１１２】
変換パラメータを決定するためには、上記のステップを反復的に実行し、識別成績が最大となるモデルパラメータを求める。識別成績の最大化には公知の最適化技術が用いられる。また、人工ニューラルネットワークやランダムフォレストを用いた回帰分析や、線形アルゴリズムによるモデル推定も適用可能である。
【０１１３】
このように求めた変換パラメータは、新規症例での処理に適用される。
【０１１４】
医用情報処理装置１の第２の変形例によれば、識別妥当性尺度の値に応じて、各パターンを表示するかどうかが判断された結果一部のパターンのみが表示されるので、同様に、効率的な診断を行うための適切な情報を操作者に提示することができる。
【０１１５】
（第２の実施形態）
第２の実施形態に係る医用情報処理装置１Ａは、被検体に関するの画像の特徴量により、疾患に関する情報の確率の信頼区間を算出し、それに基づき疾患、あるいは疾患に関わる確率区分を判定するものである。また、医用情報処理装置１Ａは、必要に応じて、信頼区間により症例数（又は、検査数）の計数範囲を選択することができるものである。
【０１１６】
第１の比較例として、被検体のＣＴ画像やＭＲＩ画像の病変領域などで画像の特徴量を算出し、その特徴量を閾値で判定して疾患の有無、良性／悪性などを判定する技術がある。このような技術の評価のための多症例での臨床研究では、陽性的中率（ＰＰＶ）や、陰性的中率（ＮＰＶ）や、それらの信頼区間が評価に用いられることがある。
【０１１７】
第１の比較例では、陽性と判定される場合、陽性的中率を疾患確率あるいは悪性の確率とみることができ、陰性と判定される例では陰性的中率を疾患のない確率あるいは良性の確率などとみることができる。しかし、陽性と判断された場合に疾患のない確率や陰性と判断された場合に疾患がある確率を表すものではない。信頼できる結果を得るためには、陽性的中率や陰性的中率の信頼区間が小さくする必要があり、そのためには十分大きな症例数を用いて評価を行うことが必要である。また、被検体に関する画像の特徴量が閾値に近いとき、あるいは逆に閾値より離れ大多数の検査より偏った位置にあるとき、症例数を数える計数範囲の偏った位置に特徴量があることになる。そのため、閾値での判定で評価した陽性的中率や陰性的中率は、当該被検体に関する特徴量に対しては計数範囲の非対称性が大きいことになり、陽性的中率や陰性的中率は、その被検体の実際の陽性確率や陰性確率とは大きく隔たりのある値となりえる。陽性的中率や陰性的中率の信頼区間も同様に、その被検体の実際の確率の範囲とは異なる可能性がある。
【０１１８】
なお、陽性的中率（ＰＰＶ）の信頼区間と、陰性的中率（ＮＰＶ）の信頼区間とは、検査を多数実施した時の平均的な性能と幅を表わすものであり、新規検査に対する性質を直接には表していない。新規検査の実際の性能値は、陽性的中率、陰性的中率での値とは異なる。被検体に関する画像の特徴量の値が閾値の近傍であるとき、この差は特に大きい。陽性的中率の信頼区間と、陰性的中率の信頼区間との両方を表示しても参考にすべき値はどちらかのみである。新規検査が陽性判定だが、陰性の確率を見たいとき、どれも実際の正答率を表していない。なお、以下、説明の便宜上、疾患の判定、疾患の確率、悪性の判定、悪性の確率、疾患有無の判定、疾患なしの確率、それらの信頼区間等の語を使用する。これは説明のための語であり、本明細書及び図面で示すようにこれらは、「疾患有無」、「疾患種別」、「予後良好／不良」、「治療適性あり／なし」の判定、これらの確率、およびこれらの信頼区間を代表するものである。
【０１１９】
第２の比較例として、研究に用いた多数の症例に関する画像の特徴量の値で良性群のヒストグラムと悪性群のヒストグラムを求め、被検体に関する画像の特徴量の値での良性の頻度、悪性の頻度から、「悪性確率＝悪性頻度／（良性頻度＋悪性頻度）」により悪性確率を求める単純ベイズ法と呼ばれる技術がある。
【０１２０】
第２の比較例では、悪性確率や陽性確率を求めることができるが、信頼区間を算出することができないという問題がある。そのため、信頼できる結果を得るためには十分に症例数を多くした大規模な臨床研究が必要になる。
【０１２１】
第３の比較例として、統計解析による様々な量、例えば画像の特徴量の標準偏差や平均値の信頼区間を使用する例が、臨床研究に広く用いられている。また、第４の比較例として、被検体の画像で測定した海馬体積が正常症例群での海馬体積の何パーセンタイルにあたるか（パーセンタイルランク）を表示するソフトウェアがある。また、体積の５パーセンタイル値と、９５パーセンタイル値を表示し、被検体の海馬体積がこの間に含まれるかを判定するソフトウェアもある。
【０１２２】
第３の比較例では、統計量の信頼区間は、複数症例からなるある群での分布の特徴を表すものであり、１例の被検体の検査データ（１例）の信頼度を示さない。従って、１検査データの解析結果の信憑性を表す数値として用いることはできない。第４の比較例では、パーセンタイル表示は被検体の海馬体積が正常範囲であるかの指標であるが、正常である確率を表すものではなく、その信頼区間でもない。正常か異常かの判定の信憑性を示すものとして使うことはできない。
【０１２３】
そこで、第２の実施形態に係る医用情報処理装置１Ａは、被検体に関する画像の特徴量により、疾患に関する情報の確率の信頼区間を算出し、それに基づき疾患あるいは疾患に関わる確率区分を判定し、また、確率の信頼区間により症例数の計数範囲を選択するものである。
【０１２４】
第２の実施形態に係る医用情報処理装置１Ａの構成は、図１に示す第１の実施形態に係る医用情報処理装置１の構成と同等であるので、説明を省略する。
【０１２５】
図１２は、医用情報処理装置１Ａの機能を示すブロック図である。
【０１２６】
処理回路１１が記憶回路１２等の非一過性の記録媒体に記憶されたコンピュータプログラムを実行することによって、医用情報処理装置１は、出力機能２３と、統計情報収集機能２４と、処理機能２５とを実現する。なお、機能２３〜２５の全部又は一部は、医用情報処理装置１にＡＳＩＣ等の回路として実現されるものであってもよい。
【０１２７】
統計情報収集機能２４は、複数の画像データの特徴に基づき疾患の頻度に関する統計情報を収集する機能と、収集した統計情報を記憶回路１２に記憶させる機能とを含む。なお、統計情報収集機能２４は、統計情報収集部の一例である。
【０１２８】
処理機能２５は、判定対象である被検体の画像データが疾患を有する確率と、その確率の信頼度を示す指標である信頼区間とを算出する処理を行う機能を有する。具体的には、処理機能２５は、統計情報取得機能２５１と、信頼区間算出機能２５２とを有する。なお、処理機能２５は、処理部の一例である。
【０１２９】
統計情報取得機能２５１は、統計情報収集機能２４によって収集された複数の画像データの特徴に基づき集計された疾患の頻度に関する統計情報を、記憶回路１２から取得する機能を有する。あるいは統計情報取得部は複数の画像データまたは画像特徴値を記憶回路１２から取得し、統計情報を算出する。なお、統計情報取得機能２５１は、統計情報取得部の一例である。
【０１３０】
信頼区間算出機能２５２は、判定対象である被検体の画像データが疾患を有する確率と、その確率の信頼度を示す指標である信頼区間とを、前記統計情報が示す頻度分布を利用して判定対象である被検体の計測値を含む任意の範囲から求める機能を含む。なお、信頼区間算出機能２５２は、信頼区間算出部の一例である。
【０１３１】
出力機能２３は、処理機能２５による確率と、確率の信頼区間とに基づく、疾患に関する情報を出力する機能を含む。例えば、出力機能２３は、第１の機能ｆ１〜第１Ｇの機能ｆ１Ｇを用いて後述するように、疾患に関する情報として、確率と、確率の信頼区間とを出力する（例えば、ディスプレイ１４に表示する）。また、出力機能２３は、第２の機能ｆ２〜第２Ｇの機能ｆ２Ｇを用いて後述するように、疾患に関する情報として、確率と、確率の信頼区間との関係から被検体の疾患に関する、あるいは疾患の確率区分(確率カテゴリ)に関する判定結果を出力する（例えば、ディスプレイ１４に表示する）。
【０１３２】
なお、機能２３〜２６の具体的な動作については、図１３〜図４４を用いて後述する。
【０１３３】
（第１の機能、及び、第２の機能）
図１３は、処理機能２５による疾患などの確率の信頼区間の算出機能と、出力機能２３による信頼区間の出力機能とを含む第１の機能ｆ１を示す。例えば、信頼区間の出力機能は、信頼区間を示す表示画面をディスプレイ１４に表示する。また、図１３は、第１の機能ｆ１と、処理機能２５による確率の信頼区間に基づく疾患の確率区分等の判定機能と、出力機能２３による判定結果の出力機能とを含む第２の機能ｆ２を示す。例えば、判定結果の出力機能は、判定結果を示す表示画面をディスプレイ１４に表示する。図１３は、医用情報処理装置１Ａの第１及び第２の機能ｆ１，ｆ２をデータフローとして示す図である。
【０１３４】
図１３の第１の機能ｆ１について先に説明する。信頼区間算出機能２５２は、解析に使用するＮ個の検査の画像データの特徴に基づき集計された疾患の頻度に関する統計情報（例えば、頻度分布（ヒストグラム））に、被検体に関する画像の特徴量の値（計測値）をプロットし、計測値を中心として幅を設けて症例数を計数して信頼区間を求める。なお、群症例数とは、複数画像データ（計数範囲に含まれる検査）の群ごとの症例数を意味する。前述の第１の比較例では、既定の２つの区域のいずれかのみで陽性的中率やその信頼区間を求める。なお、陽性的中率は、陽性判定される検査の平均的な陽性率である。また、陰性的中率は、陰性判定される検査の平均的な陰性率である。一方で、前述の第２の比較例では、信頼区間を求めることができない。
【０１３５】
それに対して、第１の機能ｆ１は、図１４Ａの矢印で示すように、被検体に関する計測値に応じた計数範囲で信頼区間を求めるものである。図１４Ａは、統計情報（例えば、頻度分布（ヒストグラム））に配置された、被検体に関する計測値を示す。疾患確率の信頼区間は、被検体が疾患であるかの確率のばらつきの範囲を表すものであるため、疾患であるかの判定結果の信憑性を示すものとして利用できる。また、信頼区間算出機能２５２は、被検体に関する計測値に応じた適切な計数範囲で信頼区間を求める。そのため、信頼区間算出機能２５２は、平均的な確率ではなく、特定の被検体に関する計測値に対する確率を示すことができる。そのため、大規模な症例数の臨床データを用いなくても、小規模スタディでも被検体に関する計測値によっては信頼区間の幅を狭くすることができ、信頼性のある利用価値のある結果が得られるという効果がある。
【０１３６】
例えば、予め、解析に使用するＮ個の検査の画像とそれぞれの検査に対応した群の値が用意されている。群は、「疾患なし」の群と、「疾患あり」の群との２つの群とする。図１５に示すように、Ｎ個の検査のうち「疾患なし」のデータの特徴量に基づくヒストグラムと、「疾患あり」のデータの特徴量からなるヒストグラムとを重ねたグラフＵ１５が図示される。このヒストグラムのグラフＵ１５に対し、被検体に関する計測値（グラフＢ１５に図示）をプロットし、計測値周辺の一定幅あたりの「疾患なし」の症例数を算出した結果がグラフＭ１５の点「ａ」で示される。同様に、一定幅あたりの「疾患あり」の症例数を算出した結果が、点「ｂ」で示される。信頼区間算出機能２５２は、疾患確率を、「ｂ／（ａ＋ｂ）」により算出し、疾患確率の信頼区間（上限値及び下限値）を、次の式（８），（９）により算出する。
【数８】
2020149684.tif 000010
【０１３７】
【数９】
2020149684.tif 000011
【０１３８】
さらに、被検体に関する計測値によっては信頼区間の幅が広く得られた疾患確率に対するばらつきが大きいことが示唆される場合もあるが、第１の機能ｆ１によれば、そのことを出力機能２３により操作者に提示できるため、操作者は、信頼性の低い判定結果を利用しないという判断が可能となる。
【０１３９】
なお、信頼区間算出機能２５２は、計数範囲を、特徴量の個々の値に応じて適切な幅に最適化するように構成してもよい。これにより、被検体に関する計測値に応じた最適な計数範囲を使用するため、信頼区間をより狭くできる可能性がある。そのため、信頼性のある結果を得られる可能性が増大する。
【０１４０】
また、信頼区間算出機能２５２は、図１４Ｂの複数の矢印で示すように、計数範囲を、信頼区間の値を使って適切な幅を選択するように構成してもよい。その場合、まず、信頼区間算出機能２５２は、多数の計数範囲に対応する多数の信頼区間のうち、幅の狭い信頼区間を決める。そして、信頼区間算出機能２５２は、悪性確率とその信頼区間を表示し、良性／悪性／不明を判定する。これにより、計数範囲を固定ではなく計測値に応じた適切な幅とすることができ、計測値がどのような値であっても、計測値に応じた適切なばらつきを抑えた判定を行うことができる。計数範囲は、計数値を中心としたある幅を持った範囲とし、その幅は例えば信頼区間下限が最も大きくなるように決める。この場合、出力機能２３は、疾患の確率の平均値、疾患の確率の信頼区間(上限値と下限値)、あるいは信頼区間の判定条件に応じた判定結果を出力することができる。
【０１４１】
図１３の第２の機能ｆ２について説明する。信頼区間算出機能２５２は、第１の機能ｆ１によって算出された、疾患確率などの信頼区間の上限値または下限値に基づいて、確率区分（以下、「確率カテゴリ」あるいは単に「カテゴリ」という）を判定して表示する。例えば、信頼区間算出機能２５２は、疾患可能性大、中、小のいずれであるかを出力する。ここで、カテゴリとは、疾患確率の信頼区間に基づく分類を意味する。例えば、疾患の可能性低、疾患の可能性中、疾患の可能性大等の分類である。
【０１４２】
（第１Ａの機能、及び、第２Ａの機能）
図１６は、処理機能２５による疾患などの確率の信頼区間の算出機能と、出力機能２３による信頼区間の出力機能とを含む第１Ａの機能ｆ１Ａを示す。また、図１６は、第１Ａの機能ｆ１Ａと、処理機能２５による確率の信頼区間に基づく疾患の確率区分等の判定機能と、出力機能２３による判定結果の出力機能とを含む第２Ａの機能ｆ２Ａを示す。第１Ａの機能ｆ１Ａ（第２Ａの機能ｆ２Ａも同様）は、計数範囲を重み関数で表すものである。図１６は、医用情報処理装置１Ａの第１Ａ及び第２Ａの機能ｆ１Ａ，ｆ２Ａに係るデータフローを示す図である。
【０１４３】
前述の最小構成の説明、および以降の説明では、単に計数範囲と表現する。計数範囲は、特徴量空間にて単一または複数の閾値で囲まれた範囲であってもよいが、一般的には重み関数（計数重み関数）で表すことができる。好適には、被検体に関する画像の特徴量の値で「１」となり、そこから離れるにしたがって値が「０」に近づく釣鐘状の重み関数が用いられる。
【０１４４】
特徴量（例えば、特徴量ベクトル）がｘ（ｘ＝ｘ１，ｘ２，…）で表されるとき、計数重み関数は、次のように表される。
（１）単純閾値：計測値がｕ（ｕ＝ｕ１，ｕ２，…）の場合、重み関数「ｗ」は全てのｉについて、ui-ri<xi<ui+riのときは「１」で、その他のときは「０」
（２）多次元ガウス関数：w=exp(-0.5(x-u)TR-1 (x-u))
なお、「Ｒ」は正値対称行列である。
【０１４５】
計数重み関数を用いた場合、信頼区間算出機能２５２は、計数範囲ｗ（「重み関数w」と同義）内のデータ総数とその中で群ｉに属するデータ数を求める。各群のデータのヒストグラムをＨｃｉ（ｘ）とし、全群のヒストグラムの和をＨａ（ｘ）とすると、計数範囲ｗ（ｘ）にあるデータの総数は、次の式（１２）で求められる。ｗ（ｘ）は前述した重み関数、好適には釣鐘状の重み関数である。また、その中で群ｉに属するデータの数は、次の式（１３）で求められる。なお、「ｉ」は群の番号、「ｘ」は特徴量ベクトルを示す。
【数１０】
2020149684.tif 000012
【０１４６】
（第１Ｂの機能、及び、第２Ｂの機能）
図１７は、処理機能２５による疾患などの確率の信頼区間の算出機能と、出力機能２３による信頼区間の出力機能とを含む第１Ｂの機能ｆ１Ｂを示す。図１７は、第１Ｂの機能ｆ１Ｂと、処理機能２５による確率の信頼区間に基づく疾患の確率区分等の判定機能と、出力機能２３による判定結果の出力機能とを含む第２Ｂの機能ｆ２Ｂを示す。第１Ｂの機能ｆ１Ｂ（第２Ｂの機能ｆ２Ｂも同様）は、計数範囲を最適化（最大化又は最小化）するものである。図１７は、医用情報処理装置１Ａの第１Ｂ及び第２Ｂの機能ｆ１Ｂ，ｆ２Ｂに係るデータフローを示す図である。
【０１４７】
計数範囲は、特徴量空間で常に一定の幅や形状である必要はなく、特徴量の値に応じて適切な幅や形状を用いるのが良い（以降、適切な計数範囲を定めると記載する）。適切な計数範囲を定めるためには、ある群に着目し、その群に所属する確率の信頼区間の下限値がなるべく大きな値になるような計数範囲を、被検体に関する計測値で「１」となる計数重み関数の中から選ぶ。あるいは、その群に所属する確率の信頼区間の上限値がなるべく小さな値になるような計数範囲を、被検体に関する計測値で「１」となる計数重み関数の中から選ぶ。
【０１４８】
図１８に示すように、「疾患なし」の確率の信頼区間下限値を最大化する計数範囲を選ぶ。グラフＵ１８にて、曲線ＲＡで表される値ＬＡの計数範囲の場合、その計数範囲の「疾患なし」の群の症例数が「ａ」であり（グラフＭ１８に図示）、「疾患あり」の群の症例数が「ｂ」であり（グラフＭ１８に図示）、「疾患なし」の確率（＝１−疾患確率）の信頼区間上限が最大になる。グラフＵ１８にて、曲線ＲＢは、被検体に関する計測値（グラフＢ１８に図示）が別の値ＬＢの時の計数範囲を示している。このように、計数範囲は、特徴量の値により最適な範囲が異なっている。つまり、個々の被検体に関する計測値に応じて適切な計数範囲が選ばれており、単に２種の範囲から選ぶ単純閾値による方法に比べて信頼性の高い判定結果を得ることができる。なお、信頼区間算出機能２５２は、疾患確率を、「ｂ／（ａ＋ｂ）」により算出し、疾患確率の信頼区間（上限値及び下限値）を、上記式（８），（９）、又は、上記式（１０），（１１）により算出する。
【０１４９】
（第１Ｃの機能、及び、第２Ｃの機能）
図１９は、処理機能２５による疾患などの確率の信頼区間の算出機能と、出力機能２３による信頼区間の出力機能とを含む第１Ｃの機能ｆ１Ｃを示す。図１９は、第１Ｃの機能ｆ１Ｃと、処理機能２５による確率の信頼区間に基づく疾患の確率区分等の判定機能と、出力機能２３による判定結果の出力機能とを含む第２Ｃの機能ｆ２Ｃを示す。前述の第１Ｂの機能ｆ１Ｂで説明したように、複数の群に対して計数範囲の最適化を行うと、群の数分の計数範囲が得られることになる。第１Ｃの機能ｆ１Ｃ、及び、第２Ｃの機能ｆ２Ｃでは、各群の計数範囲で求められた信頼区間の上限値あるいは下限値を用いて１つの群を選択し、その群の計数範囲を他の群に適用して信頼区間を計算しなおすものである。図１９は、医用情報処理装置１Ａの第１Ｃ及び第２Ｃの機能ｆ１Ｃ，ｆ２Ｃに係るデータフローを示す図である。
【０１５０】
信頼区間算出機能２５２は、各群について、その群の信頼区間の下限を最大化する計数範囲を決める。これを繰り返し、全群の計数範囲の中で信頼区間下限最大値が最も大きい計数範囲を選び、他の群に適用して疾患確率と信頼区間を計算する。図１９において、図２０Ａ〜図２０Ｄに示すように４種類の脳腫瘍を群として用いる場合で説明する。
【０１５１】
図２０Ａ〜図２０Ｄに示すように、複数の検査の最も細かい群分類は「グリオーマＧｒ１，２」、「グリオーマＧｒ３，４」、「ＰＣＮＳＬ（中枢神経性原発悪性リンパ腫）」、「脳転移」、「髄膜腫」、「下垂体腺腫」、「その他」であるが、これらをグループ化し「グリオーマ」、「ＰＣＮＳＬ」、「脳転移」、「その他」の４つの群(クラス)を使用している。
【０１５２】
信頼区間算出機能２５２の特徴的な部分は、次の４つの概略手順（１）〜（４）において、（３）の手順である。
（１）１つの疾患群をとった時、その疾患群の信頼区間の下限を最大化する計数範囲を決める。
（２）これを繰り返し全ての群(クラス)の計数範囲の中で信頼区間下限最大値が最も大きい計数範囲を選び主クラス(主群)とする。
（３）主クラス(主群)の計数範囲を他の群(クラス)に再度適用し、各群(クラス)の疾患確率とその信頼区間を求める。
（４）主クラス(主群)の信頼区間下限と、他の群(クラス)の信頼区間上限に基づいて確率カテゴリを判定する。
【０１５３】
以下、信頼区間算出機能２５２の処理手順を具体的に説明する。
ｉ）研究データは確定診断されており、様々な脳腫瘍が含まれている。それぞれの脳腫瘍の種類ごとに特徴量の値のヒストグラムが作成されている。選択された４つの大分類を判定群とした場合、４つの統合したヒストグラムが作成される。
ii）被検体での特徴量の計測値の周囲に計数範囲を設定し、「グリオーマ」の確率の信頼区間下限を最大化する計数範囲を求める（図２０Ａの破線）。同様に、他の群の信頼区間下限を最大化する計数範囲をそれぞれ決定する（図２０Ｂ〜図２０Ｄの破線）。そして各々の群に属する確率（平均）と群に属する確率の信頼区間をそれぞれ算出する。
iii）図２０Ａ〜図２０Ｄに示す各群の信頼区間下限を示す丸印「●」を比較し、それが最も大きいグリオーマ（図２０Ａ）の計数範囲を選択する。グリオーマの信頼区間下限を示す位置の表示形態（例えば、背景色）を他の位置と異ならせることで区別する。このグリオーマが主クラス(主群)となる。操作者は別の群を主クラス(主群)として選択しなおすことが可能である。
iv）表示するのは選択された図２０Ａの計数範囲から得られる主クラスの情報のみである。
v）選択された主クラス「グリオーマ」の計数範囲で、各群に属する確率（平均）と信頼区間を計算し表示する。
vi）選択されたグリオーマの信頼区間下限と、その他の群の信頼区間上限の値を表示し、その値がどのカテゴリに属するかを判定して表示する（図２０Ａに示すカテゴリ「可能性大」）。
【０１５４】
図２１Ａ〜図２１Ｃは、判定結果の表示例を示す図である。図２１Ａ〜図２１Ｃはそれぞれ表示されてもよいし、同一の表示画面上に並列表示されてもよい。特徴量としての「ＣＢＶ（Cerebral Blood Volume）」は、脳血液量を示し、「Ｋｔｒａｎｓ」は、造影剤がＥＥＳ（Eextracellular Extravascular Space）に漏出する移行速度定数を示す。図２１Ｂに示すように、主群(主クラス)として選択された「グリオーマ」が視認可能なように、例えば、「グリオーマ」の部分に、他のクラスとは異なる色が付与されて表示される。また、「グリオーマ」の平均確率の値等も表示される。加えて、サブ群である他のカテゴリについても、平均確率の値等が表示される。
【０１５５】
また、図２１Ｃに示すように、病変部位にマークがされた画像（例えば、ＤＷＩ（diffusion weighted image）画像）等を表示することもできる。「ＦＬＡＩＲ（Fluid-Attenuated Inversion Recovery）」は、ＭＲＩにおいて水信号が抑制された画像である。「Ｔ１Ｗ＿ＣＥ（T1 Weighted contrast enhanced）」は、ＭＲＩにおいてコントラストが強調されたＴ１強調画像である。「ＣＢＶ（Cerebral Blood Volume）」は、脳血流量像である。また、図２１Ｃは、造影剤がＥＥＳ（Extracellular Extravascular Space）に漏出する移行速度定数（Ktrans）を示す画像と、血管内の血漿ボリューム（Vp）を示す画像と、ＥＥＳボリューム（Ve）を示す画像と、血管に戻る移行速度定数（Kep）を示す画像とを示す。
【０１５６】
なお、１つの群内での計腫範囲の上記手順ii）を行い、主群の選択とそれに続く処理iii）〜v）を行わず、手順vi）にて、手順ii）で求めた確率の平均値と信頼区間およびカテゴリを表示してもよい。この構成は、表示する確率の各群の平均値の総和が「１」にならないなど、統計上の専門的な分析には向かないが、主群を選ばないため各群の処理が同等であり、公平な選択が可能である。従って、ルーチン検査用の解析ソフトウェアに好適な構成である。
【０１５７】
また、１つの群内での計腫範囲の上記手順ii）に加えて、計数範囲の群間の選択iii）〜v）も行う構成も可能である。その後に、手順vi）を実施する。この構成の場合、十分な分析を行うためには操作者は別の群を主群にする場合と比較する必要があり、主群を切り替えるために煩雑な操作が必要になる。一方、各群の平均値の総和が「１」になり、統計上の厳密で専門的な分析に向いている。従って、この構成は研究用の専門的分析ソフトウェアに好適な構成である。
【０１５８】
（第１Ｄの機能、及び、第２Ｄの機能）
図２２は、処理機能２５による疾患などの確率の信頼区間と、処理機能２５によるグラフの算出機能と、出力機能２３による信頼区間とグラフとの出力機能とを含む第２Ｄの機能ｆ１Ｄを示す。図２２は、第１Ｄの機能ｆ１Ｄと、処理機能２５による確率の信頼区間に基づく疾患の確率区分等の判定機能と、出力機能２３による判定結果の出力機能とを含む第２Ｄの機能ｆ２Ｄを示す。第１Ｄの機能ｆ１Ｄ（第２Ｄの機能ｆ２Ｄも同様）は、ヒストグラムと計数範囲を表示するものである。図２２は、医用情報処理装置１Ａの第１Ｄ及び第２Ｄの機能ｆ１Ｄ，ｆ２Ｄに係るデータフローを示す図である。
【０１５９】
出力機能２３は、複数検査の群ごとのヒストグラムに、被検体に関する計測値に応じて決定された計数範囲と、被検体に関する計測値とを重畳したグラフを作成することで表示画面を生成する（例えば、図１４Ａ及び図１４Ｂに示す計数範囲及び計測値を含む表示画面）。
【０１６０】
（第１Ｅの機能、及び、第２Ｅの機能）
図２３は、処理機能２５による疾患などの確率の信頼区間に基づく群所属確率のカテゴリ判定と、カテゴリ区域図との算出機能と、出力機能２３によるカテゴリ判定結果とカテゴリ区域図との出力機能とを含む第１Ｅの機能ｆ１Ｅを示す。図２３は、第１Ｅの機能ｆ１Ｅと、出力機能２３による、カテゴリ判定結果とカテゴリ区域図との同時の出力機能とを含む第２Ｅの機能ｆ２Ｅを示す。第１Ｅの機能ｆ１Ｅ（第２Ｅの機能ｆ２Ｅも同様）は、カテゴリの区域図をグラフとして表示するものである。図２３は、医用情報処理装置１Ａの第１Ｅ及び第２Ｅの機能ｆ１Ｅ，ｆ２Ｅに係るデータフローを示す図である。信頼区間算出機能２５２は、群症例数の計数の際、畳み込み演算を行う。
【０１６１】
図２４Ａに示すように、各群の群所属確率の信頼区間上限値、および信頼区間下限値の値に応じてどのカテゴリに判定するかを表すカテゴリ区域図を表示する。群所属確率とは、ある特徴量の値を持つ被検体がある群に属するものである確率を意味する。例えば、良性確率の信頼区間下限値と、悪性確率の信頼区間下限値を用いてカテゴリ区域図を表示する。判定のための水準、例えば、「９９％」及び「９０％」等により、図２４Ａ及び図２４Ｂに示す区域（ａ）は「確実に悪性」に対応し、区域（ｂ）は「おそらく悪性」に対応し、区域（ｃ）は「不明（悪性寄り）に対応し、区域（ｄ）は「不明（良性寄り）」に対応し、区域（ｅ）は「おそらく良性」に対応し、区域（ｆ）は「確実に良性」に対応する。
【０１６２】
カテゴリ区域図内に被検体の信頼区間下限値および上限値をプロットし、判定結果の区域ごとに異なる色を付与して表示するすることにより、被検体のデータでは確率がどの程度であるかをわかりやすく提示できる。また、信頼区間の値によりどう判定されるか、あるいはどう判定されたかをわかりやすく示すことができる。
【０１６３】
図２５に示す、閾値範囲の幅に対する疾患確率を示すグラフＵ２５は、疾患確率の信頼区間下限の最大値を丸印「〇」で示す。当該最大点をカテゴリ区域図Ｒ２５の縦軸に充てる。一方、閾値範囲の幅に対する疾患確率を示すグラフＢ２５は、疾患確率の信頼区間上限の最小値を丸印「〇」で示す。当該最小点をカテゴリ区域図の横軸に充てる。なお、疾患確率の平均値および信頼区間上限値、下限値は、第１の実施形態における「群尺度」に相当する。また、信頼区間上限値および下限値を計数範囲の選択や最適化に利用する場合、疾患確率の平均値および信頼区間上限値、下限値は、第１の実施形態の「識別妥当性尺度」に相当する。
【０１６４】
カテゴリ区域図Ｒ２５で（ａ）から（ｆ）で示される区域は、それぞれ「確実に悪性」、「おそらく悪性」、「不明（悪性寄り）」、「不明（良性寄り）」、「おそらく良性」、「確実に良性」の６つのカテゴリに対応している。このカテゴリ区域図上のどこにプロットされるかで、カテゴリを判定する。
【０１６５】
なお、カテゴリ区域図Ｒ２５は、軸ごとに判定してもよいが、その場合に限定されるものではない。例えば、カテゴリ区域図は、図２６Ａに示すような斜線の判定曲線及び判定直線の両方を含むカテゴリ区域図であってもよく、また、図２６Ｂに示すような判定直線のみを含むカテゴリ区域図であってもよい。
【０１６６】
また、群として「良性」、「悪性」の２つの群を用いた例を図２７に示す。図２７に示すように、悪性らしさ、良性らしさに応じたカテゴリ区域が用いられる。一方で、群として「治療効果あり」、「治療効果なし」の２つの群を用いた例を図２８に示す。図２８に示すように、有効か無効かに応じたカテゴリ区域が用いられる。
【０１６７】
（第１Ｆの機能、及び、第２Ｆの機能）
図２９は、処理機能２５による疾患などの確率の信頼区間と、画像表示と、領域表示との算出機能と、出力機能２３による信頼区間と、画像表示と、領域表示との出力機能とを含む第１Ｆの機能ｆ１Ｆを示す。図２９は、第１Ｆの機能ｆ１Ｆと、処理機能２５による群所属確率のカテゴリ判定機能と、出力機能２３によるカテゴリ判定結果の出力機能とを含む第２Ｆの機能ｆ２Ｆ示す。第１Ｆの機能ｆ１Ｆ（第２Ｆの機能ｆ２Ｆも同様）は、画像及び特徴量算出領域を表示するものである。図２９は、医用情報処理装置１Ａの第１Ｆ及び第２Ｆの機能ｆ１Ｆ，ｆ２Ｆに係るデータフローを示す図である。
【０１６８】
出力機能２３は、解析の基になる画像と、特徴量を計算した対象領域（例えば、関心領域）を画面に表示する。操作者は、表示されている各特徴量の値が画像と照らして妥当であるか確認することができる。
【０１６９】
（第１Ｇの機能、及び、第２Ｇの機能）
処理機能２５による第１Ｇの機能ｆ１Ｇ（第２Ｇの機能ｆ２Ｇも同様）は、施設の有病率を設定するものである。図３０は、医用情報処理装置１Ａの第１Ｇ及び第２Ｇの機能ｆ１Ｇ，ｆ２Ｇに係るデータフローを示す図である。
【０１７０】
操作者またはシステム管理者は、各群の当施設での発生確率を設定（変更）可能である。この場合、上記手順ii）とv）の確率と信頼区間の値は設定した発生確率を用いて補正される。
【０１７１】
この場合、検査を行う前に被検体がどのような疾患であるかの疾患確率（クラス事前確率）は、解析に使用するデータベースに含まれる症例の各疾患の比率（クラス構成比）と同一とみなすことが多い。これは人種なども含めてデータベースの症例に登録する基準と、施設での検査対象であるかを判断する基準が同一とみなせることを前提としている。しかし、施設の特性（クリニックであるか専門病院であるか）や対象検査の前にどのような別の検査をしているかにより、疾患確率（群の事前確率）は変わり得るものである。この場合、両者を同一とみなして算出疾患確率や信頼区間は不正確なものとなりえる。症例データベースでの群の比率と、施設での発生確率あるいは、その被検体の経緯に依存する群の確率との違いを吸収するために、群の事前確率の補正が必要になる。
【０１７２】
【数１１】
2020149684.tif 000013
【０１７３】
オッズ比Ｏ１は、次の式（１５）で表される。
【数１２】
2020149684.tif 000014
変形すると、次の式（１６）のようになる。
【数１３】
2020149684.tif 000015
被検体の検査を実施した施設でも、右辺の「感度／（１−特異度）」は変わらないが、群が「Ａ」である疾患確率は、ａ２に変わると考える。
【数１４】
2020149684.tif 000016
【０１７４】
式（１８）をＰ２について解くと、次の式（１９）で表される。ここで、ｋを次の式（２０）のように定義した。
【数１５】
2020149684.tif 000017
【０１７５】
結果として、補正された確率値と確率の信頼区間は次のように計算できる。
事前確率が、「ａ２」の時の予測される疾患確率Ｐ２は、事前確率が「ａ１」であるであるときの疾患確率Ｐ１から、上記次の式（１９），（２０）から求める。事前確率が「ａ２」の時の信頼区間についても同様に求める。
【０１７６】
（判定項目のバリエーション）
信頼区間算出機能２５２が、疾患の有無の判定を行う例を中心に説明したが、その場合に限定されるものではない。例えば、判定項目は、群の取り方や、特徴量空間の次元数により様々な応用が可能である。判定項目を次に例示する。
（I）信頼区間を算出する量（確率の信頼区間）
（I−ｉ）１変数の場合：病変の存在確率、重症確率、悪性の確率
（I−ii）複数変数の場合：各疾患の存在確率、カテゴリ分類、ステージ分類、病変タイプの各々の確率、各組織型の確率
（I−iii）予測：薬剤の種類ごとの効果ありの確率、治療(処置)の効果がある確率
（II）群
（II−ｉ）被検体の状態（疾患の有無、種類）
（II−ii）病変の有無、病変の悪性／良性の判定、病変のカテゴリ判定（カテゴリ０〜５）、病変のタイプ分類（サブタイプなど）、組織型の判定
（III）予測
（III−ｉ）薬効有無、治療適性（効果有無）、予後（重症化する／しない）、重篤イベントの発生有無、完治する／しない
（IV）判定する属性の数
（IV−ｉ）２水準、１変数〜多水準、複数変数（多変数による複数群への分類として考えることが可能）
【０１７７】
なお、「多水準、複数変数」については、多変数による複数群への分類として考えることができる。図３１はその一例である。
【０１７８】
図３１に示すように、肝臓線維化の程度を、Ｆ０〜Ｆ４の４つの群（ステージ）で表し、これらの群を２つの変数（特徴量１と特徴量２）で判定する場合を考える。前述の実施例によれば被検体に関する計測値（特徴量１と特徴量２）に応じて計数範囲（グラフ中Ｆ０からＦ４で図示）が設定され、Ｆ０の確率（平均値）、Ｆ１の確率（平均値）、…、Ｆ４の確率（平均値）、およびそれらの信頼区間が算出される。これらの信頼区間は、図３２のグラフＬ３２に上限・下限・平均値がプロットされる。例えば、これら確率の上限値が閾値を超えている群のうち、線維化程度が最も大きい群が、判定結果として提示される。図３２のグラフＲ３２は他の例で、確率の下限値が閾値を超え、かつ下限値が最も大きい群を判定結果として提示する。
【０１７９】
（特徴量のバリエーション）
信頼区間算出機能２５２が、疾患の有無の判定を行う例を中心に説明したが、その場合に限定されるものではない。特徴量は、次のようなものでもよい。
（Ｖ）画像処理の組み合わせで算出する特徴量
（Ｖ−ｉ）円形度、球形度
（Ｖ−ii）平均値、energy, entropy, contrast etc. (Texture特徴量)
（Ｖ−iii）最大直径、体積
（VI）従来の手法による機能マップの統計値
（VI−ｉ）ADC(Apparent diffusion coefficient), FA (fractional anisotropy)
（VI−ii）rCBF (relative regional cerebral blood flow, rCBV (relative regional cerebral blood volume, MTT(mean transit time), TTP (time to peak)
（VI−iii）Ktrans, Kep (Tofts model)
（VII）機械学習により決定する特徴量
（VII−ｉ）Neural network
（VII−ii）Support vector machine
（VII−iii）Random forest, extremely random tree
（VII−iv）Logistic regression
（VIII）画像の特徴量のほかに、心拍のR-R間隔や血液検査の結果、症状有無など、非画像の情報を加えてもよい。
【０１８０】
機械学習による画像の特徴量の例として敵対的生成ネットワーク（Generative Adversarial Networks）を使用する例を挙げる。敵対的生成ネットワーク（図３３に図示）を用いると比較的少数のユニットで構成される中間層の各ユニットが画像の何らかの特徴を表すように学習されることが知られている。このようにして得られた中間ユニットの出力を特徴量として利用することができる。
【０１８１】
また、各種機械学習を回帰に用いることができる。図３４は、ＳＶＭ（Support vector machine）を回帰に用いた場合である。回帰されたパラメータを特徴量として利用することができる。回帰にはまた人工ニューラルネットワークやランダムツリーを応用することもできる。グラフＵ３４，Ｍ３４，Ｂ３４の説明については、図１５でしたグラフＵ１５，Ｍ１５，Ｂ１５の説明と同様であるので、説明を省略する。
【０１８２】
（計数範囲の最適化方法の例）
以下、良性／悪性を判定する場合において、良性群の計数範囲を最適化する場合を例にとって説明する。特徴量の次元が２以上の場合について、図３５Ａ及び図３５Ｂを用いて説明する。
【０１８３】
図３５Ａに示すように、第１に、個別の特徴量の上閾値／下閾値の１つを、変数の分布幅に対して、例えば１％拡げる。１つの変数に対して信頼区間上限を拡げる／信頼区間下限を拡げるの２つがあるから、拡げるやり方は、「変数の数×２」の可能性があり、その中で最も信頼区間の幅が小さいものを選ぶ。これを繰り返して徐々に閾値範囲を拡げゆく。「信頼区間の幅が小さい」とは、信頼区間上限の最小値を求める場合、信頼区間上限の最小値が小さいことを指す。上閾値／下閾値の計測値とのずれ（上側閾値幅と下側閾値幅）の差が小さいことも同時に加味する。
【０１８４】
第２に、上と同じだが、閾値範囲を一般化して釣鐘状積分重み関数で範囲を表す。釣鐘状積分重み関数のパラメータを少し変えた多数の中から最も信頼区間の幅が小さいものを選ぶ。
【０１８５】
図３５Ｂに示すように、第３に、計測値を含む閾値範囲や積分重み関数を乱数で多数生成して、横軸に範囲の体積、縦軸に信頼区間の幅(信頼区間の上限や下限)をプロットする。その包絡線の最小値を求める。言い換えれば乱数で生成した積分範囲での信頼区間の幅の最小値を求める。
【０１８６】
第４に、閾値範囲は一定幅で分割する。一変数の場合はＮ分割、ｍ変数の場合は、Ｎ＾ｍ個の範囲に分割する。１変数の分割数Ｎを変えて信頼区間の幅を計算する。
【０１８７】
第５に、これらは信頼区間上限値を最小化、または信頼区間下限値を最大化するように計数範囲のパラメータを最適化する問題であり、上記のほか様々な最適化手法を適用することができる。
【０１８８】
（カテゴリの判定のバリエーション）
以下の例で、「下限」は信頼区間下限または信頼区間下限の最大値を表す。また、「上限」は、信頼区間上限、又は、信頼区間上限の最小値を表す。
【０１８９】
図３６に示すように、疾患なし群と疾患あり群のカテゴリを判定する例について考える。「疾患なし」が確実な被検体を選びたい場合や、少しでも疾患の可能性のある被検体をもれなく選びたい場合や、少しでも疾患なしの可能性のある被検体を選びたい場合や、疾患有が確実な被検体を選びたい場合に有効である。
【０１９０】
図３６に示すような表示により、追加の診断なしに疾患を否定できる被検体を選べる。追加の診断なしに確定できる被検体を選べる。「下限＋上限＝１」なので、左右の２つの判定は同等である。
【０１９１】
図３７は、疾患Ａ，Ｂ，Ｃ，Ｄの４つの群の場合で、各群の疾患のカテゴリを判定する一例を示す。複数疾患の中から可能性の高いものを選び、疾患の鑑別を行う。
【０１９２】
図３８は、肝線維化のステージＦ０，Ｆ０以外，Ｆ１，Ｆ２，Ｆ３，Ｆ４の各群のカテゴリを判定する例を示す。複数疾患の中から可能性の高いものを選び、疾患の鑑別を行う。なお、図３８の下部は、図３７の下部と同等であるので、表記を省略する。
【０１９３】
図３９Ａ及び図３９Ｂは、薬剤Ａ，Ｂ，Ｃの効果をそれぞれカテゴリに判定する例を示す。複数の薬剤の中で被検体に効果のある薬剤を選び、薬効の予測を行う。なお、同様の判定方法を複数の治療方法の中で被検体に効果のある治療法を選ぶためのカテゴリの判定や、合併症の少ない治療法を選ぶためのカテゴリの判定に利用することもできる。
【０１９４】
図４０Ａ及び図４０Ｂは、重篤化イベントの各々の発生率をカテゴリとして判定する例、および、いくつかの観点で良好なＱＯＬ（Quality of Life）を回復・維持できるかをカテゴリとして判定する例を示す例を示す。冠動脈イベントの種類ごとに発生の可能性がどの程度あるか、被検体ごとに予測する。同時に、良好なＱＯＬを回復・維持する可能性を予測する。
【０１９５】
（階層的な判定の例）
図４１Ａ及び図４１Ｂは、疾患確率を３階層で判定する例である。図４１Ｃ及び図４１Ｄは、各治療法による治療の効果を３階層で判定する例である。図４１Ｅ及び図４１Ｆは、ＱＯＬ良好の確率及び重篤化の確率を３階層で判定する例である。このような判定には、画像データまたは特徴量のデータ、その他非画像検査データ、疾患の種類、選択した治療法とその効果有無、各種ＱＯＬの良好／不良および各種重篤化の発生有無について複合的に記録されている症例データを用いる。
【０１９６】
判定は、まず、全症例を４つの疾患の群（大分類群）に分割し、特徴量の計数範囲を設定する。疾患１〜４の各大分類群の症例数を求め、各疾患の確率の信頼区間により各疾患の可能性をカテゴリで判定するとともに、信頼区間の値を用いて大分類群（疾患の種類）を判定する（図４１Ａ及び図４１Ｂ）。
【０１９７】
次に、判定された群（仮に疾患１と判定されたとする）の症例に絞った症例を対象として、各治療法の効果の判定を行う。大分類群が判定された疾患１である症例を、「治療法１の効果あり」と、「治療法１の効果なし」と、「治療法２の効果あり」と、「治療法２の効果なし」と、「治療法３の効果あり」と、「治療法３の効果なし」と、「治療法４の効果あり」と、「治療法４の効果なし」との８つの群（中分類群）に分割する。特徴量の計数範囲を設定することで、これら中分類群の症例数を各々求める。治療法１の効果の有無の確率は「治療法１の効果あり」と「治療法１の効果なし」の症例数をもとに効果ありの確率およびその信頼区間を求める。治療法２〜４についても同様に、効果ありの確率と信頼区間を求める。これらをもとに、効果ありの可能性をカテゴリで判定するとともに、信頼区間の値を用いて中分類群を判定する（図４１Ｃ及び図４１Ｄ）。図４１Ｄに示す例では、治療法１群が選択されたものとする。
【０１９８】
最後に、大分類群が「疾患１」で中分類群が「治療法１の効果あり」である症例を対象として、ＱＯＬ良好の確率および重篤化の確率を判定する。あるいは、中分類群に「治療法２の効果なし」を加える構成もありうる。これら対象症例をさらに「身体機能良好」群と「身体機能不良」群（小分類群１）とに、また同じ対象症例を「活動量良好」群と「活動量不良」群（小分類群２）とに、さらに同じ対象症例を「突然死」群と「非突然死」群（小分類群３）に、さらに同じ対象症例を「心筋梗塞」群と「非心筋梗塞」群（小分類群４）とに分割する。特徴量の計数範囲を設定することでこれら小分類群の症例数を各々求める。「身体機能良好」の確率と確率の信頼区間は「身体機能良好」の症例数と「身体機能不良」の症例数から求める。「活動量良好」の確率と信頼区間、「突然死」の確率と信頼区間、「心筋梗塞」の確率と信頼区間も同様に求める。これら信頼区間の値を用いて小分類１〜４のカテゴリを判定する（図４１Ｅ及び図４１Ｆ）。ＱＯＬ良好の確率および重篤化の確率については小分類１から４の各々について、ＱＯＬ良好か不良化、あるいは重篤化が発生するかしないかを小分類群１〜４の各々で判定する。
【０１９９】
（計数範囲を設定及び変更するＧＵＩの例）
ヒストグラムの初期表示では、最適化で求めた計数範囲を、図４２に示すように表示する。図４２に示すグラフにて、丸印「〇」の特徴量の値と四角印「□」の特徴量の値とで挟まれた範囲が計数範囲である。操作者が丸印「〇」と四角印「□」を入力インターフェース１３としてのマウスで操作すると、当該範囲が拡げられるか、あるいは狭められる。その結果、群に所属する確率およびその信頼区間が再計算され、更新された値が表示される。
【０２００】
図４３に示すグラフには、計数範囲が拡げられた場合、狭められた場合の信頼区間上限と下限が予め表示されている。したがって、操作者が実際に計数範囲を変更する前に、再計算される信頼区間下限、上限がどのような値になるか知ることができ、試行錯誤せずに、どのように変更すればよいかを判断でき、操作時間の短縮につながる。
【０２０１】
図４４は、特徴量が３変数の場合のヒストグラム表示を示す。多変数の場合、複数次元のうち１次元のみのグラフを表示し、他の次元は計数範囲を固定する。
【０２０２】
図４２〜図４４に示すように、自動決定した閾値範囲をもとに、予測陽性率の信頼区間の上下閾値の値に対する変化を示すグラフを、特徴量ごとに表示する。そして、閾値を変更する操作を行うとグラフが更新される。
【０２０３】
以上説明した少なくとも１つの実施形態によれば、効率的な診断を行うための適切な情報を操作者に提示することができる。
【０２０４】
なお、図１２に示す処理回路１１は、機能２３〜２５に加え、図２に示す処理回路１１の機能２１〜２３の全部又は一部を含むこともできる。また、処理機能２１，２５は、処理部の一例である。抽出機能２２は、抽出部の一例である。出力機能２３は、出力部の一例である。関心領域設定機能２１１は、関心領域設定部の一例である。特徴算出機能２１２は、特徴算出部の一例である。尺度取得機能２１３は、尺度取得部の一例である。統計情報収集機能２４は、統計情報収集部の一例である。統計情報取得機能２５１は、統計情報取得部の一例である。信頼区間算出機能２５２は、信頼区間算出部の一例である。
【０２０５】
以上、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これらの実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【０２０６】
１  医用情報処理装置
１１  処理回路
１２  記憶回路
２１，２５  処理機能
２１１  関心領域設定機能
２１２  特徴算出機能
２１３  尺度取得機能
２２  抽出機能
２３  出力機能
２４  統計情報収集機能
２５１  統計情報取得機能
２５２  信頼区間算出機能
【図１】
2020149684.tif 000018
【図２】
2020149684.tif 000019
【図３】
2020149684.tif 000020
【図４】
2020149684.tif 000021
【図５】
2020149684.tif 000022
【図６】
2020149684.tif 000023
【図７】
2020149684.tif 000024
【図８】
2020149684.tif 000025
【図９】
2020149684.tif 000026
【図１０】
2020149684.tif 000027
【図１１】
2020149684.tif 000028
【図１２】
2020149684.tif 000029
【図１３】
2020149684.tif 000030
【図１４Ａ】
2020149684.tif 000031
【図１４Ｂ】
2020149684.tif 000032
【図１５】
2020149684.tif 000033
【図１６】
2020149684.tif 000034
【図１７】
2020149684.tif 000035
【図１８】
2020149684.tif 000036
【図１９】
2020149684.tif 000037
【図２０Ａ】
2020149684.tif 000038
【図２０Ｂ】
2020149684.tif 000039
【図２０Ｃ】
2020149684.tif 000040
【図２０Ｄ】
2020149684.tif 000041
【図２１Ａ】
2020149684.tif 000042
【図２１Ｂ】
2020149684.tif 000043
【図２１Ｃ】
2020149684.tif 000044
【図２２】
2020149684.tif 000045
【図２３】
2020149684.tif 000046
【図２４Ａ】
2020149684.tif 000047
【図２４Ｂ】
2020149684.tif 000048
【図２５】
2020149684.tif 000049
【図２６Ａ】
2020149684.tif 000050
【図２６Ｂ】
2020149684.tif 000051
【図２７】
2020149684.tif 000052
【図２８】
2020149684.tif 000053
【図２９】
2020149684.tif 000054
【図３０】
2020149684.tif 000055
【図３１】
2020149684.tif 000056
【図３２】
2020149684.tif 000057
【図３３】
2020149684.tif 000058
【図３４】
2020149684.tif 000059
【図３５Ａ】
2020149684.tif 000060
【図３５Ｂ】
2020149684.tif 000061
【図３６】
2020149684.tif 000062
【図３７】
2020149684.tif 000063
【図３８】
2020149684.tif 000064
【図３９Ａ】
2020149684.tif 000065
【図３９Ｂ】
2020149684.tif 000066
【図４０Ａ】
2020149684.tif 000067
【図４０Ｂ】
2020149684.tif 000068
【図４１Ａ】
2020149684.tif 000069
【図４１Ｂ】
2020149684.tif 000070
【図４１Ｃ】
2020149684.tif 000071
【図４１Ｄ】
2020149684.tif 000072
【図４１Ｅ】
2020149684.tif 000073
【図４１Ｆ】
2020149684.tif 000074
【図４２】
2020149684.tif 000075
【図４３】
2020149684.tif 000076
【図４４】
2020149684.tif 000077