(57)【要約】【課題】目に係る診断情報と組となっていない情報も有効に活用しつつ、機械学習処理を行い、当該機械学習の効率を向上しつつ、診断等の処理に役立てることのできる、情報処理装置、機械学習データの製造方法、及びプログラムを提供する。【解決手段】目の症状に関する情報のパターンを表す情報を取得し、目に係る診断情報を入力としたときのニューラルネットワークの出力と、当該入力した診断情報に対応する目の症状に関する情報との相違に関する情報と、受け入れた目の症状に関する情報のパターンに係る情報とを用いてニューラルネットワークのパラメータを制御して機械学習データを生成する。【選択図】図１
2019152977.tif 000002
【特許請求の範囲】
【請求項１】
出力情報のパターンを表す情報を取得する取得手段と、
入力情報と、対応する出力情報とを関連付ける機械学習データを生成する手段であって、前記入力情報を入力したときに学習器が出力する情報と、当該入力情報に対応する前記出力情報との相違に関する情報と、前記取得した前記出力情報のパターンに係る情報とを用いて前記学習器のパラメータを制御して機械学習データを生成する学習処理手段と、
を含み、
前記学習処理手段により生成された機械学習データを、処理の対象となる前記入力情報に対応する、前記出力情報を推定する処理に供する情報処理装置。
【請求項２】
目の症状に関する情報のパターンを表す情報を取得する取得手段と、
目に係る診断情報と、対応する目の症状に関する情報とを関連付ける機械学習データを生成する手段であって、目に係る診断情報を入力としたときの学習器の出力と、当該入力した診断情報に対応する目の症状に関する情報との相違に関する情報と、前記受け入れた目の症状に関する情報のパターンに係る情報とを用いて前記学習器のパラメータを制御して機械学習データを生成する学習処理手段と、
を含み、
前記学習処理手段により生成された機械学習データを、処理の対象となる目に係る診断情報に対応する、目の症状に関する情報を推定する処理に供する情報処理装置。
【請求項３】
請求項２に記載の情報処理装置であって、
前記取得手段は、目の症状に関する情報を複数受け入れて、当該複数の目の症状に関する情報に基づく教師なし学習により得られる目の症状に関する情報のパターンを表す情報を生成して取得する情報処理装置。
【請求項４】
請求項２または３に記載の情報処理装置であって、
前記目に係る診断情報は、網膜構造を表す情報であり、前記目の症状に関する情報は、視野に関する情報である情報処理装置。
【請求項５】
請求項３に記載の情報処理装置であって、
前記取得手段が受け入れる目の症状に関する情報は少なくとも視野を表す情報を含む情報である情報処理装置。
【請求項６】
目の症状に関する情報を複数受け入れて、当該複数の目の症状に関する情報に基づく教師なし学習により得られる診断情報のパターンを表す情報を取得する工程と、
目に係る診断情報と、対応する目の症状に関する情報とを関連付ける機械学習データを生成する工程であって、目に係る診断情報を入力としたときの学習器の出力と、当該入力した診断情報に対応する目の症状に関する情報との相違に関する情報と、前記受け入れた目の症状に関する情報のパターンに係る情報とを用いて前記学習器のパラメータを制御する工程と、
を含む、
機械学習データの製造方法。
【請求項７】
コンピュータを、
目の症状に関する情報のパターンを表す情報を取得する取得手段と、
目に係る診断情報と、対応する目の症状に関する情報とを関連付ける機械学習データを生成する手段であって、目に係る診断情報を入力としたときの学習器の出力と、当該入力した診断情報に対応する目の症状に関する情報との相違に関する情報と、前記受け入れた目の症状に関する情報のパターンに係る情報とを用いて前記学習器のパラメータを制御して機械学習データを生成する学習処理手段と、
として機能させ、
前記学習処理手段により生成された機械学習データが、処理の対象となる目に係る診断情報に対応する、目の症状に関する情報を推定する処理に供されるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、眼科医療用の情報等を処理する情報処理装置、機械学習データの製造方法、及びプログラムに関する。
【背景技術】
【０００２】
近年では種々の情報処理に機械学習が用いられるようになっており、医療用の情報についても、機械学習の結果を用いた情報処理を行うものが提案されている（例えば特許文献１）。
【０００３】
一方、緑内障など、非可逆的な視機能の喪失を伴う疾病については、早期の発見が求められる。しかしながら、緑内障では、その確定診断のための検査が一般に時間のかかるものであり、負担の大きいものであるため、予めスクリーニングによって緑内障の可能性の有無を簡便に検出する方法が求められている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１５−１１６３１９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記従来の機械学習を用いる装置等においては、入力されるデータと出力するデータとが組となっているデータを用いていわゆる教師あり学習処理を行っているが、医療用の情報等では、入力するデータと関わりなく、出力されるデータが検査結果の情報として蓄積されている例が多く、必ずしも入力するデータと出力されるデータとが組となっていない場合がある。
【０００６】
例えば、緑内障の診断において、視野検査の結果については過去から数多くの蓄積があるが、網膜構造などの情報は、必ずしも視野検査の結果と組となって保持されているものではない。
【０００７】
本発明は上記実情に鑑みて為されたもので、組となっていない情報も有効に活用しつつ、機械学習処理を行い、当該機械学習の効率を向上しつつ、診断等の処理に役立てることのできる、情報処理装置、機械学習データの製造方法、及びプログラムを提供することを、その目的の一つとする。
【課題を解決するための手段】
【０００８】
上記従来例の問題点を解決する本発明は、情報処理装置であって、目の症状などの患者の症状に関する情報等、出力情報のパターンを表す情報を取得する取得手段と、目の診断情報など、患者の診断情報等、入力情報と、対応する出力情報とを関連付ける機械学習データを生成する手段であって、入力情報を入力としたときの学習器の出力と、当該入力情報に対応する出力情報との相違に関する情報と、前記受け入れた入力情報のパターンに係る情報とを用いて前記学習器のパラメータを制御して機械学習データを生成する学習処理手段と、を含み、前記学習処理手段により生成された機械学習データを、処理の対象となる診断情報等の入力情報に対応する、症状に関する情報等の出力情報を推定する処理に供することとしたものである。
【発明の効果】
【０００９】
本発明によると、入力される情報と出力する情報とが組となっていない情報も有効に活用しつつ、機械学習処理を行い、当該機械学習の効率を向上でき、また、当該機械学習の結果を診断等の処理に利用できる情報処理装置が提供される。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施の形態に係る情報処理装置の構成ブロック図である。
【図２】本発明の実施の形態に係る情報処理装置の一例に係る機能ブロック図である。
【図３】本発明の実施の形態に係る情報処理装置の動作例を表すフローチャート図である。
【図４】本発明の実施の形態に係る情報処理装置の出力する情報の例を表す説明図である。
【発明を実施するための形態】
【００１１】
本発明の実施の形態について図面を参照しながら説明する。本発明の実施の形態に係る情報処理装置１は、図１に例示するように、制御部１１と、記憶部１２と、操作部１３と、表示部１４と、入出力部１５とを含んで構成される。
【００１２】
ここで制御部１１は、ＣＰＵ等のプログラム制御デバイスであり、記憶部１２に格納されたプログラムに従って動作する。本実施の形態の例では、この制御部１１は、出力情報のパターンを表す情報を取得し、入力情報と、対応する出力情報とを関連付ける機械学習データを生成する際に、入力情報を入力したときに学習器が出力する情報と、当該入力情報に対応する出力情報との相違に関する情報と、取得した、出力情報のパターンに係る情報とを用いて学習器のパラメータを制御して機械学習データを生成する。そして制御部１１は、当該生成された機械学習データを、処理の対象となる入力情報に対応する、出力情報を推定する処理に供する。
【００１３】
具体的に、本実施の形態の一例では、制御部１１は、出力情報としての、目の症状に関する情報のパターンを表す情報を取得する。またこのときには、制御部１１は、入力情報としての、目に係る診断情報と、対応する目の症状に関する情報とを関連付ける機械学習データを生成する。具体的に制御部１１は、目に係る診断情報を入力情報として入力したときのニューラルネットワーク等、非線形関数を用いる学習器の出力（出力情報）と、当該入力した入力情報である診断情報に対応する出力情報である、目の症状に関する情報との相違に関する情報と、取得した、目の症状に関する情報のパターンに係る情報とを用いて学習器のパラメータを制御して機械学習データを生成する。
【００１４】
そしてこの制御部１１は、学習処理により生成された機械学習データを、処理の対象となる目に係る診断情報に対応する、目の症状に関する情報を推定する処理を行ってもよい。この制御部１１の詳しい動作の内容については後に説明する。
【００１５】
なお、本実施の形態の以下の説明では、入力情報と出力情報とは上記の例のものとして説明するが、これは一例であり、他の種類の入力情報や出力情報であっても同様の処理を適用できる。例えば入力情報は、目以外の部位に関する患者の診断情報等であってもよい。また出力情報は、目以外の部位に関する症状に関する情報等であってもよい。このとき、入力情報が関係する患者の部位と、出力情報が関係する患者の部位とが互いに異なっていてもよい。さらに、本実施の形態の入力情報と出力情報との例は、人体に関係する情報に限られず、非線形関数を用いる学習器が適用される種々の分野、例えば言語解析、経営情報解析、気象情報の解析、等であっても、出力情報についてパターンの情報を抽出できるものであれば、いかなる情報であっても構わない。またここで出力情報は多次元の情報であってよい。
【００１６】
記憶部１２は、ディスクデバイスやメモリデバイス等であり、制御部１１によって実行されるプログラムを保持する。このプログラムは、コンピュータ可読かつ非一時的な記録媒体に格納されて提供され、この記憶部１２に格納されたものであってもよい。
【００１７】
操作部１３は、キーボードやマウス等であり、利用者の指示を受け入れて、制御部１１に出力する。表示部１４は、ディスプレイ等であり、制御部１１から入力される指示に従い、情報を表示出力する。
【００１８】
入出力部１５は、例えばＵＳＢ（Universal Serial Bus）等の入出力インタフェースであり、本実施の形態の一例では、処理の対象となる種々の情報を外部の装置（例えば種々の検査装置やカードリーダー等）から受け入れて制御部１１に出力する。
【００１９】
次に制御部１１の動作について説明する。本実施の形態では、この制御部１１は、図２に例示するように、学習処理部２０と、推定処理部３０とを含む。また学習処理部２０は、パターン取得部２１と、学習データ生成部２２とを含んで構成される。
【００２０】
本実施の形態の一例において、このパターン取得部２１は、目の症状に関する情報を複数受け入れて、当該複数の目の症状に関する情報に基づく教師なし学習により得られる目の症状に関する情報のパターンを表す情報を生成して取得する。
【００２１】
具体的にこのパターン取得部２１は、目の症状に関する情報として、視野に関する情報を複数受け入れる。ここで視野に関する情報は、視野検査の結果の情報、すなわち視野中心から所定の角度（例えば１０度または３０度）以内に予め設定された、複数の測定点Ｐ１，Ｐ２…における視野感度の情報でよい。なお、以下の例において各測定点は固有の識別情報（連番など）によって識別可能とされているものとする。
【００２２】
パターン取得部２１は、受け入れた視野に関する情報に基づいて、教師なし学習により症状の進行のパターンを抽出する。例えば緑内障であれば、視野欠損が上側から進行する症例や、下側から進行する症例等、症例がいくつかのパターンに分類できることが知られている。また、いずれの症例においても、比較的近い位置にある視野の測定点間では比較的似た視野感度の情報が得られること（共起性を有すること）が知られている。
【００２３】
そこでこのパターン取得部２１では、教師なし学習として、例えば非負値行列因子分解（ＮＭＦ）や、主成分分析（ＰＣＡ）、オートエンコーダ、混合分布を用いる方法等の方法を用いて、受け入れた視野に関する情報から症状の進行のパターンを抽出する。一例として非負値行列因子分解を用いる場合、視野に関する情報のそれぞれに含まれる各測定点での視野感度の情報を、予め定めた順序としておく。つまりパターン取得部２１は、複数（以下ではＮ個（Ｎは２以上の整数）とする）の視野に関する情報に含まれる視野感度の情報を、対応する測定点Ｐi（ｉ＝１，２，…）の順で配列した、ベクトル情報Ｙj（ｊ＝１，２，…，Ｎ）としておく。
【００２４】
パターン取得部２１は、このベクトル情報Ｙj（ｊ＝１，２，…，Ｎ）の列ベクトルを連接して配列した行列
【数１】
2019152977.tif 000003
を求め、この行列Ｙを、一対の非負値の行列ＷYとΘYとの積に分解する。すなわち、
【数２】
2019152977.tif 000004
とする。ここで、ΘYは、Ｋ個の列ベクトルθYを連接して配列した行列を転置した基底行列
【数３】
2019152977.tif 000005
である。この基底行列ΘYに含まれる各行（転置前の各列θYi（ｉ＝１，２，…Ｋ））が進行パターンを表す。
【００２５】
この非負値行列因子分解の具体的方法については、D. D. Lee and H. S. Seung, Algorithms for Non-negative Matrix Factorization, In NIPS, Vol. 13, pp. 556-562, 2000等に記載されているので、ここでの詳しい説明を省略する。
【００２６】
本実施の形態のこの例に係るパターン取得部２１は、ここで得た基底行列ΘYを、目の症状に関する情報のパターンを表す情報として出力する。
【００２７】
学習データ生成部２２は、学習器としての所定のニューラルネットワークの学習を行う。ここでのニューラルネットワークは例えば、ＣＮＮ（畳み込みニューラルネットワーク）であり、目に係る診断情報と、それに対応する目の症状に関する情報とを関連付ける機械学習データである。本実施の形態のここでの例では、目の症状に関する情報は、上述の通り、視野検査の結果をベクトルの形に表現したものでよい。
【００２８】
また目に係る診断情報は、例えば、光干渉断層計（ＯＣＴ）による、眼底計測結果（網膜の構造に関する情報）または、眼底写真の画像データ等のいずれかであり、以下の例では光干渉断層計（ＯＣＴ）による、眼底計測結果であるものとする。また以下の説明では、ニューラルネットワークとしてＣＮＮを用いた場合を例として説明するが、本実施の形態のニューラルネットワークはＣＮＮに限られず、残差（Residual）ネットワーク（ResNet：Kaiming He, et.al., Deep Residual Learning for Image Recognition, https://arxiv.org/pdf/1512.03385v1.pdf）など、他の構造を有するニューラルネットワークであってもよい。ただし、ニューラルネットワークの出力は、視野検査の結果を表すベクトルの次元と同じ次元のベクトル量であるようにしておく。
【００２９】
学習データ生成部２２は、目に係る診断情報と当該診断情報が得られたときの目の症状に関する情報との組の入力を受け入れる。この情報の組は、ＯＣＴによる検査と、視野検査とを実質的に同時に（例えば数日の期間内に）行った特定の対象者から得られたものを用いる。このような情報の組の数は、パターン取得部２１が用いる、目の症状に関する情報（診断情報と組となっていないもの）の総数よりもその数は少ないものの、学習データの生成に十分な数だけあるものとする。
【００３０】
学習データ生成部２２は、受け入れた目に係る診断情報と当該診断情報が得られたときの、正解情報である目の症状に関する情報との組を順次選択し、選択した組に含まれる診断情報を、ニューラルネットワークの入力としたとき、当該ニューラルネットワークが出力する情報（目の症状に関する情報との差が演算可能な情報）ｙpredと、選択した組に含まれる目の症状に関する情報（正解情報）ｙtrueとの二乗残差平均（Ｌ２罰則項）ｌ2（ｙpred，ｙtrue）を求める。
【００３１】
またこの学習データ生成部２２は、パターン正則化項Ｒ（ｙpred）を次のように求める。すなわち、このパターン正則化項を
【数４】
2019152977.tif 000006
と定義する。ここでｙpredNMFは、
【数５】
2019152977.tif 000007
で規定される。ただし、
【数６】
2019152977.tif 000008
とする。
【００３２】
これは、ニューラルネットの予測結果であるｙpredを、症状の進行のパターン（基底）がつくる部分空間に正射影したものに相当し、パターン正則化項全体としては、この正射影の結果と、ニューラルネットの予測結果であるｙpredとのユークリッド距離の二乗とする。
【００３３】
これにより、予測結果が症状の進行のパターンから外れるほど、当該予測結果の部分空間への正射影と、予測結果自体とのユークリッド距離は大きくなり、より大きい正則化が行われることとなる。
【００３４】
学習データ生成部２２は、診断情報と目の症状に関する情報との組に基づくＬ２罰則項ｌ2（ｙpred，ｙtrue）とパターン正則化項Ｒ（ｙpred）との和（重みつきの和ｌ2（ｙpred，ｙtrue）＋λ・Ｒ（ｙpred）、でよい（ここでλは実数であり、重みの一例である））を小さくする方向に、ニューラルネットワークのパラメータを制御して設定し直す。
【００３５】
このような設定方法は、一般的なニューラルネットワークにおける勾配降下法として種々の知られた方法を適宜選択して採用できるので、ここでの詳しい説明を省略する。
【００３６】
学習データ生成部２２は、受け入れた目に係る診断情報と当該診断情報が得られたときの、正解情報である目の症状に関する情報との組のそれぞれについて、上述の処理を繰り返して行い、ニューラルネットワークのパラメータを制御して、機械学習データを生成する。
【００３７】
また、制御部１１の推定処理部３０は、学習処理部２０が生成した機械学習データを用い、学習処理部２０と同じニューラルネットワークを用いて、目の症状に関する情報を推定する。具体的にこの例では、目に係る診断情報として光干渉断層計（ＯＣＴ）による、眼底計測結果（網膜の構造に関する情報）を用いて、機械学習データを得ているので、この推定処理部３０も、診断の対象となる対象者の光干渉断層計（ＯＣＴ）による、眼底計測結果を、機械学習データによりパラメータを設定したニューラルネットワークに入力する。そして当該ニューラルネットワークが出力する情報（ここではそれぞれの要素が視野内の各測定点における視野感度を表すベクトルの情報）を、視野内の各測定点における視野感度の推定結果（つまり目の症状に関する情報の推定結果）として出力する。
【００３８】
なお、本実施の形態の例では、制御部１１が学習処理部２０と推定処理部３０との双方を備えることとしたが、本実施の形態の制御部１１は、学習処理部２０のみを備えて、機械学習データを出力するものであってもよい。
【００３９】
この場合、他のコンピュータ装置において、当該学習処理部２０が出力する機械学習データを用いた推定処理部３０としての処理を実行させる。
【００４０】
またここで推定処理部３０と学習処理部２０とが用いるニューラルネットワークは同じものであるとしたが、この例に限られず、例えば学習処理部２０がＣＮＮ、推定処理部３０が２層の全結合ネットワークを用いてもよい。この場合、学習処理部２０が出力する機械学習データと、学習処理部２０とが用いるニューラルネットワークとを用いた推定処理の結果を用いて、当該推定処理部３０が用いるニューラルネットワークを学習させる、いわゆる蒸留を行って、推定処理部３０が用いるニューラルネットワークのパラメータである機械学習データ（学習処理部２０が出力する機械学習データとの区別のため蒸留データと呼ぶ）を得ておく。そして推定処理部３０は、この蒸留データと、蒸留時に用いたニューラルネットワークとを用い、診断の対象となる対象者の光干渉断層計（ＯＣＴ）による、眼底計測結果をこのニューラルネットワークに入力して、視野内の各測定点における視野感度の推定結果（つまり目の症状に関する情報の推定結果）を得て出力する。
【００４１】
［動作］
本実施の形態は以上の構成を備えており、次のように動作する。本発明の実施の形態に係る情報処理装置１は、学習処理を行うよう指示を受けると、図３に例示するように、光干渉断層計（ＯＣＴ）による検査と、視野検査とを実質的に同時に（例えば数日の期間内に）行った特定の対象者から得られた、光干渉断層計（ＯＣＴ）による、眼底計測結果（網膜の構造に関する情報）と、視野検査の結果との組（以下の説明ではペアデータと呼ぶ）と、光干渉断層計（ＯＣＴ）による検査の結果と組となっていない視野検査の結果の情報（以下、ノンペアデータと呼ぶ）との入力をそれぞれ複数受け入れる（Ｓ１）。
【００４２】
そして情報処理装置１は、処理Ｓ１で受け入れたノンペアデータに基づいて非負値行列因子分解（ＮＭＦ）等の教師なし学習処理を行い、目の症状に関する情報のパターンを表す情報を生成して取得する（Ｓ２）。一例としてこの処理Ｓ２では、非負値行列因子分解（ＮＭＦ）により、目の症状に関する情報である視野検査の結果の情報のパターンを表す情報として、いくつかの互いに異なる代表的なパターンを基底として得る。
【００４３】
情報処理装置１は、次に、ＣＮＮのニューラルネットワークを設定する（Ｓ３）。この設定では、過去に学習して得た機械学習データがあれば、ＣＮＮのニューラルネットワークに当該機械学習データのパラメータを設定する（学習済みモデルをロードする）こととし、過去に学習して得た機械学習データがなければ、ＣＮＮのニューラルネットワークを初期化する。初期化の方法は、広く知られた方法を採用してよい。
【００４４】
そして情報処理装置１は、処理Ｓ１で受け入れた複数組のペアデータを一つずつ選択して次のループを繰り返し実行する（Ｓ４）。すなわち情報処理装置１は、選択したペアデータに含まれる眼底計測結果（網膜の構造に関する情報）をＣＮＮのニューラルネットワークに入力したときの出力を得る（Ｓ５）。
【００４５】
情報処理装置１は、当該処理Ｓ５で得た出力と、選択したペアデータに含まれる正解情報である視野検査の結果の情報との二乗残差平均（ｌ2（ｙpred，ｙtrue））を求め、また、処理Ｓ２で得たパターンを表す情報と処理Ｓ５で得た出力で得た情報との適合性を評価する値（Ｒ（ｙpred））を求める。この値は例えばニューラルネットの出力を、処理Ｓ２で得たパターンを表す情報（基底）がつくる部分空間に正射影し、この正射影の結果と、ニューラルネットの出力とのユークリッド距離の二乗等、ニューラルネットの出力が、処理Ｓ２で得た情報が表すパターンから外れるほど、大きい値となるような値を用いる。
【００４６】
情報処理装置１は、ここで求めた二乗残差平均ｌ2（ｙpred，ｙtrue）と、パターンへの適合性を表す値Ｒ（ｙpred）とを用いた損失関数を演算し、当該損失関数の演算結果に基づいてニューラルネットワークのパラメータ（機械学習データ）を更新する（Ｓ６）。ここで損失関数は、上記二乗残差平均とパターンへの適合性を表す値との重みづけ和でよい。なお、損失関数に基づく更新の処理自体は、一般的なニューラルネットワークにおける勾配降下法として種々の知られた方法を適宜選択して採用できる。
【００４７】
情報処理装置１は、処理Ｓ５，Ｓ６を、処理Ｓ１で受け入れた複数組のペアデータの一つ一つについて繰り返し実行し、選択していないペアデータがなくなると、その時点での機械学習データ（学習済みモデル）を記憶部１２に格納して処理を終了する。
【００４８】
この機械学習データは、外部のコンピュータに出力されて処理の対象となる目に係る診断情報に対応する、目の症状に関する情報を推定する処理に用いられてもよいし、この情報処理装置１自身が利用して、処理の対象となる目に係る診断情報に対応する、目の症状に関する情報を推定する処理（推定処理）を行ってもよい。
【００４９】
情報処理装置１は、推定処理を行うよう指示を受けると、記憶部１２に格納されている機械学習データ（学習済みモデル）を読み出して、学習処理において用いたものと同じＣＮＮのニューラルネットワークに当該機械学習データのパラメータを設定する（学習済みモデルをロードする）。
【００５０】
そして情報処理装置１は、光干渉断層計（ＯＣＴ）による検査の結果の情報の入力を受けて、当該光干渉断層計（ＯＣＴ）による検査の結果の情報を設定したニューラルネットワークに入力する。そして情報処理装置１は、当該入力に対するニューラルネットワークの出力する情報（推定された視野検査の結果に相当する情報）を得て、当該情報を表示部１４に表示するなどして出力する。
【００５１】
このように本実施の形態によると、過去に蓄積されているノンペアデータを有効に活用して症状の推定に資するニューラルネットワークを学習処理することが可能となり、機械学習の効率を向上でき、また、当該機械学習の結果を診断等の処理に利用できる。
【００５２】
［パターン抽出の他の例］
また、教師なし学習により症状の進行のパターンを抽出する方法は、非負値行列因子分解（ＮＭＦ）の方法に限られない。
【００５３】
情報処理装置１の制御部１１が用いる教師なし学習の方法は、主成分分析（ＰＣＡ）やオートエンコーダであってもよい。
【００５４】
具体的に主成分分析を用いる場合、情報処理装置１の制御部１１は、主成分ベクトルΘYを、視野検査の結果から学習することとなる。そして制御部１１は、この主成分ベクトルを用い、ニューラルネットワークが出力する推定結果ｙpredと、主成分分析により得られたパターンとの一致度を表す値Ｒ（ｙpred）を、
【数７】
2019152977.tif 000009
ただし、
【数８】
2019152977.tif 000010
として規定する。なお、
【数９】
2019152977.tif 000011
は、視野検査結果から得た視野検査の情報の要素ごとの平均（各測定点での視野感度の情報の平均を配列したベクトル）である。
【００５５】
またオートエンコーダを用いる場合、情報処理装置１の制御部１１は、ニューラルネットワークが出力する推定結果ｙpredと、オートエンコーダにより得られたパターンとの一致度を表す値Ｒ（ｙpred）を、
【数１０】
2019152977.tif 000012
として規定する。
【００５６】
ここで関数ｆY(Ｘ）は、視野検査の結果に相当するベクトルをＸとして、学習したオートエンコーダに入力したときの、オートエンコーダの出力を表す。ここでのオートエンコーダで用いるニューラルネットワークは例えば隠れ層が１層のみで、入力層よりも次元の低い層となっている全結合層とし、入力層での活性化関数を、例えばシグモイド関数とし、隠れ層での活性化関数をtanh関数としておく。オートエンコーダであるため、出力層は入力層と同じ次元とする。
【００５７】
この例のオートエンコーダにより得られたパターンとの一致度を表す値Ｒ（ｙpred）はつまり、ペアデータを学習しようとするニューラルネットワークの出力を、ノンペアデータで学習したオートエンコーダのニューラルネットワークに入力し、当該入力した情報と、オートエンコーダのニューラルネットワークの出力（パターンに一致させたもの、つまりオートエンコーダが視野検査結果から学習した部分空間に、ペアデータを学習しようとするニューラルネットワークの出力を射影したものに相当する）した情報との差が大きいほど、パターンに一致していないことを表すものとするものである。
【００５８】
さらに、本実施の形態のある例では、教師なし学習の方法としてガウスモデル（ＧＭ）を用いてもよい。この例では、情報処理装置１の制御部１１は、ノンペアデータに含まれる複数の視野検査結果のベクトル情報から、当該ベクトルの要素ごとの平均ベクトルμYを求め、分散共分散行列ΣYを推定する。この分散共分散行列の推定方法は広く知られた方法を採用できるので、ここでの詳しい説明は省略する。
【００５９】
制御部１１は、視野検査結果がガウス分布に従って発生していると仮定して、ペアデータを学習しようとするニューラルネットワークの出力の、上記平均ベクトルμYと、分散共分散行列ΣYとで特徴づけられるガウスモデルに対する一致度を表す値Ｒ（ｙpred）を、マハラノビス距離を用いて、
【数１１】
2019152977.tif 000013
として規定する。
【００６０】
このように、非負値行列因子分解（ＮＭＦ），主成分分析（ＰＣＡ），オートエンコーダ、あるいはガウスモデル（ＧＭ）のいずれを用いる方法でも、制御部１１は、いずれかの方法で得られたパターンと、ペアデータを学習しようとするニューラルネットワークの出力との一致度を表す値であるパターン正則化項の値Ｒ（ｙpred）を求め、診断情報と目の症状に関する情報との組に基づくＬ２罰則項ｌ2（ｙpred，ｙtrue）とこのパターン正則化項Ｒ（ｙpred）との和（重みつきの和ｌ2（ｙpred，ｙtrue）＋λ・Ｒ（ｙpred）、でよい）を小さくする方向に、ニューラルネットワークのパラメータを制御して設定し直すことで、ニューラルネットワークの学習処理を行う。
【００６１】
［パターンを表す情報を予め指定しておく例］
また、ここまでの説明では、ノンペアデータに基づいてノンペアデータが表すニューラルネットワークの出力に関するパターンを推定し、当該パターンとニューラルネットワークの出力とが一致するように、ニューラルネットワークの学習処理を行うようにしていたが、本実施の形態はこれに限られない。
【００６２】
すなわちパターンを表す情報（パターン正規化項Ｒ（ｙpred））は、ノンペアデータではなく、事前の知識によってユーザが任意に設定してもよい。
【００６３】
一例として、視野検査結果において、互いに隣接する測定点同士での測定結果は相関する傾向にあることを考慮し、隣接する測定点の値の差を小さくする方向にニューラルネットワークのパラメータを制御して設定し直すことで、ニューラルネットワークの学習処理を行うこととしてもよい。この例では、例えば各測定点について、当該測定点の値の推定結果（ニューラルネットワークの出力）と、隣接する測定点の値の推定結果との差を求め、各測定点で求めた差の二乗和をＲ（ｙpred）とすればよい。
【００６４】
［可視化処理］
上述の機械学習の結果を用いると、診断処理のほか、次のように、患者の状態に係る入力情報（上述の例であれば、光干渉断層計（ＯＣＴ）による、眼底計測結果（網膜の構造に関する情報）や眼底写真の画像データなどの目に係る診断情報）と、患者の状態に係る出力情報（上述の例であれば視野検査の結果に相当する情報）との関係を可視化する処理を行うこともできる。
【００６５】
すなわち、本実施の形態の情報処理装置１は、上述の方法で機械学習したニューラルネットワーク等の非線形の学習器による入出力関係の微分を用い、患者の状態に係る入力情報についての教師なし学習で得た複数のパターン（入力情報のパターン）Ｐinｉ（ｉ＝１，２…）と、患者の状態に係る出力情報についての教師なし学習で得た複数のパターン（出力情報のパターン）Ｐoutｊ（ｊ＝１，２，…）との間の重みＤij（Ｐinｉと、Ｐoutｊとの間の重み）を求める。
【００６６】
そして情報処理装置１は、入力情報におけるどのようなパターンが、出力情報におけるどのパターンにより大きく寄与するのかを示す情報として、この重みの情報を提示することとしてもよい。
【００６７】
具体的に、本実施の形態の一例では、上述の方法で機械学習したニューラルネットワーク等の非線形の学習器を用いる。以下の例では、この学習器にｘを入力したときの当該学習器による推定結果θの出力をｙpred（ｘ；θ）とする。
【００６８】
また以下では、入力される情報ｘの変化量Δｘと、この変化量Δｘに対する、出力情報ｙpredの変化量Δｙpredとが、Δｙpred＝ＡΔｘという線形の関係にあると仮定する。
【００６９】
このとき、Ａは、（Δｘ，Δｙがベクトルであるので）複数のベクトルａ1からａMを行方向に配列して転置したヤコビ行列（そのｍ行ｌ列目の成分をａmlとする）であり、
【数１２】
2019152977.tif 000014
で与えられる。情報処理装置１は、上述の方法で機械学習した学習器を用い、このヤコビ行列Ａを求めておく。なお、ｙpred_mは、推定結果θの出力のうち、ｍ番目の成分を表す。また∂ｙpred_m／∂ｘlは、ｘlを変化させたときのｙpred_mの変化量（勾配）を表す。
【００７０】
また、情報処理装置１は、光干渉断層計（ＯＣＴ）による検査の結果と組となっていない、複数のノンペアの視野検査の結果の情報ＸVFと、視野検査の結果の情報と組となっていない、複数のノンペアの光干渉断層計（ＯＣＴ）による検査の結果の情報ＸRTとを得て、次の処理を行う。
【００７１】
すなわち、情報処理装置１は、複数のノンペアの視野検査の結果の情報ＸVFと、複数のノンペアの光干渉断層計（ＯＣＴ）による検査の結果の情報ＸRTとのそれぞれについて、非負値行列因子分解（ＮＭＦ）等の教師なし学習の方法を用いて、それぞれの基底行列ＢVF＝ΘＸVFと、ＢRT＝ΘＸRTとを得る。
【００７２】
このとき、視野検査の結果の情報ＸVFの微小変化ΔＸVFは、
【数１３】
2019152977.tif 000015
と記述でき、また、光干渉断層計（ＯＣＴ）による検査の結果の情報ＸRTの微小変化ΔＸRTもまた、
【数１４】
2019152977.tif 000016
と記述できる。
【００７３】
この関係を、上記に仮定したΔｙpred＝ＡΔｘの関係に当てはめて、ΔＸVF＝ＡΔＸRTと記述できるとし、ｚVFとｚRTとの関係を、一般化逆行列を考慮して求めると、
【数１５】
2019152977.tif 000017
となる。ただし、
【数１６】
2019152977.tif 000018
である。
【００７４】
情報処理装置１は、先に求めた基底行列ＢVF，ＢRTを用いて、（１６）式によりＤを求める。
【００７５】
なお、情報処理装置１は、基底行列の演算に非負値行列因子分解を用いる場合など、Ｄのどの成分も負になるはずがない場合には、（１６）式で求めたＤの成分のうち、負のものについては無視する（使用しない）ものとする。
【００７６】
ここで各基底行列ＢVF，ＢRTに含まれる各行（転置前の各列）のベクトルは、それぞれ視野検査の結果のパターンＰVF1，ＰVF2，…ＰVFｉ（ｉ＝１，２，…）と、光干渉断層計（ＯＣＴ）による検査の結果のパターンＰRT1，ＰRT2，…ＰRTｊ（ｊ＝１，２，…）とを表す（対応する視野の位置でのそれぞれの検査結果の値に対応する）ので、ここで求めたＤの各成分Ｄijは、各パターンＰVFiとＰRTjとの間の関係の重みを表すものとなる。
【００７７】
そこで情報処理装置１は、図４に例示するように、各基底行列ＢVF，ＢRTに含まれる各行が表すパターンの画像（Ｐ）を表示し、その間の関係重み（無視されたものを除く）とを、二部グラフのように表示する。このとき、グラフのリンクに相当する線分を、Ｄの大きさが大きいほど太い線として表示する。
【００７８】
この可視化処理によると、患者の状態に係る入力情報の代表的なパターンと、患者の状態に係る出力情報の代表的なパターンとの間でのどのような寄与の関係があるかを、学習器の微分情報を用いて構成し、視覚的に示すことができ、症状の進行等の理解に供することができる。
【符号の説明】
【００７９】
１  情報処理装置、１１  制御部、１２  記憶部、１３  操作部、１４  表示部、１５  入出力部、２０  学習処理部、２１  パターン取得部、２２  学習データ生成部、３０  推定処理部。
【図１】
2019152977.tif 000019
【図２】
2019152977.tif 000020
【図３】
2019152977.tif 000021
【図４】
2019152977.tif 000022