JP7737099B2 - 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム - Google Patents

脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム

Info

Publication number
JP7737099B2
JP7737099B2 JP2022536448A JP2022536448A JP7737099B2 JP 7737099 B2 JP7737099 B2 JP 7737099B2 JP 2022536448 A JP2022536448 A JP 2022536448A JP 2022536448 A JP2022536448 A JP 2022536448A JP 7737099 B2 JP7737099 B2 JP 7737099B2
Authority
JP
Japan
Prior art keywords
clustering
classifier
treatment
subject
brain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2022536448A
Other languages
English (en)
Other versions
JPWO2022014682A1 (ja
Inventor
雄人 柏木
智磯 徳田
雄史 高原
光男 川人
歩 山下
宙人 山下
雄希 酒井
潤一郎 吉本
剛 岡田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hiroshima University NUC
ATR Advanced Telecommunications Research Institute International
Shionogi and Co Ltd
Original Assignee
Hiroshima University NUC
ATR Advanced Telecommunications Research Institute International
Shionogi and Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hiroshima University NUC, ATR Advanced Telecommunications Research Institute International, Shionogi and Co Ltd filed Critical Hiroshima University NUC
Publication of JPWO2022014682A1 publication Critical patent/JPWO2022014682A1/ja
Application granted granted Critical
Publication of JP7737099B2 publication Critical patent/JP7737099B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/70ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mental therapies, e.g. psychological therapy or autogenous training
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room
    • A61B5/004Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0042Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/16Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
    • A61B5/165Evaluating the state of mind, e.g. depression, anxiety
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30016Brain
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Biomedical Technology (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Psychiatry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Medicinal Chemistry (AREA)
  • Data Mining & Analysis (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • General Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Computing Systems (AREA)
  • Fuzzy Systems (AREA)
  • Hospice & Palliative Care (AREA)
  • Psychology (AREA)
  • Social Psychology (AREA)
  • Developmental Disabilities (AREA)
  • Quality & Reliability (AREA)
  • Radiology & Medical Imaging (AREA)

Description

この発明は、複数の装置において脳機能画像法により計測される脳機能結合相関値のパターンをクラスタリングする技術、より特定的には、脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システムに関する。
(データ駆動型のクラスタリング方法)
近年の人工知能技術、特に、データ駆動型人工知能技術の進展により、音声認識、翻訳、画像認識などの分野で、一部、人間の能力に匹敵し、あるいは、一部の領域では、人間の能力を凌駕するような応用も実現されている(たとえば、特許文献1)。
医療技術の分野においても、画像診断などで、深層学習などの機械学習が利用されることが増えてきた。深層学習は、多層のニューラルネットワークを用いた機械学習であり、画像認識の分野では、深層学習の1つである畳込みニューラルネットワーク(Convolutional Neural Network、以下、CNNという)を用いた学習方法が従来の方法に比べて非常に高い性能を示すことが知られている(たとえば、特許文献2)。
たとえば、大腸がんの内視鏡による画像診断などでは、診断の正確度が人間を上回るような診断機器が実用化されるに至っている(非特許文献1)。
ただし、これらの人工知能技術は、機械学習の分類上は、正解データと入力データ(たとえば、画像データ)との組を大量に準備して、これを入力として、人工知能に学習処理をさせるという、いわゆる「教師あり学習」の範疇に入っているものがほとんどである。
一方で、データ駆動型人工知能の応用用途としては、与えられたデータを、その特徴量に基づいて、いくつかのクラスタに分類するというタスクの実行もある。この場合は、正解データの存在しない、いわゆる「教師なし学習」や、少量の「正解ラベル付きの学習データ」による学習と大量の「正解ラベルなし学習データ」による学習とを組み合わせた「「半教師あり学習」などが知られている(たとえば、特許文献3)。
たとえば、特許文献3では、「半教師あり学習は、比較的少ないラベル付きデータとラベル無しデータとに基づいて学習を行う学習方法であり、例えばラベル付きデータ(状態データSと判定データLを含む教師データT)を用いて分類を行う学習モデルを生成し、該学習モデルとラベル無しデータ(状態データS)とを用いて、該学習モデルに対して追加で学習を行うことで学習の精度を向上させるブートストラップ法や、ラベル付きデータ及びとラベル無しデータのデータ分布に基づいてグループ分けすることで分類器としての学習モデルを生成するグラフベースアルゴリズム等を含む。」とされている。ただし、この例にもあるとおり、「半教師あり学習」においては、教師データが存在するのは、少数の学習データであり、これで、まずは分類器を生成した後、大量の「正解ラベル無し学習データ」用いることで、この分類器自身に対して再学習することなどが前提となっていることになる。
(バイオマーカー)
以下では、人工知能技術による判別やクラスタリングを適用する分野として、医学分野を例にとる。
生体内の生物学的変化を定量的に把握するため、生体情報を数値化・定量化した指標のことを「バイオマーカー」と呼ぶ。
FDA(米国食品医薬品局)はバイオマーカーの位置づけを、「正常なプロセスや病的プロセス、あるいは治療に対する薬理学的な反応の指標として客観的に測定・評価される項目」としている。また疾患の状態や変化、治癒の程度を特徴づけるバイオマーカーは、新薬の臨床試験での有効性を確認するためのサロゲートマーカー(代用マーカー)として使われる。血糖値やコレステロール値などは、生活習慣病の指標として代表的なバイオマーカーである。尿や血液中に含まれる生体由来の物質だけでなく、心電図、血圧、PET画像、骨密度、肺機能なども含まれる。またゲノム解析やプロテオーム解析が進んできたことによって、DNAやRNA、生体蛋白等に関連したさまざまなバイオマーカーが見出されている。
バイオマーカーは、疾患にかかった後の治療効果の測定だけでなく、疾患を未然に防ぐための日常的な指標として疾患の予防に、さらに副作用を回避した有効な治療法を選択する個別化医療への応用が期待されている。
たとえば、肺疾患に対して、遺伝子情報を用いて、疾患への罹患の可能性を判断するためのバイオマーカーについて開示がある(特許文献4)。特許文献4では、「バイオマーカー」または「マーカー」とは、「その生体システムの生理学的な状態の特徴を示すものとして、客観的に測定することが可能な生物学的な分子である」とされている。そして、この特許文献4では、「通常、バイオマーカー測定値は典型的にはタンパク質またはポリペプチドである発現産物の定量的な測定に関する情報である。本発明は、バイオマーカー測定値をRNA(翻訳前)レベルでまたはタンパク質レベル(翻訳後修飾も含まれ得る)で決定することを想定している」とされている。そして、特許文献4では、このようなバイオマーカー測定値に対する「分類システム」として使用する分類器としては、決定木、ベイジアン分類器、ベイジアン信念ネットワーク、k-最近傍法、事例ベース推論、およびサポートベクトルマシンなどが例示されている。
一方で、神経・精神疾患の場合、現状の診断は、DSM-5(精神障害の診断・統計マニュアル(Diagnostic and Statistical Manual of Mental Disorders)第5版)などに基づいて、いわゆる症候に基づく診断となっていることもあり、生化学的もしくは分子遺伝学的観点から客観的な指標として利用可能な分子マーカーなども研究されているものの、検討段階というべき状況である。
もっとも、NIRS(Near-infraRed Spectroscopy)技術を用いて、生体光計測により計測されたヘモグロビン信号から特徴量に応じて、統合失調症、うつ病などの精神疾患について分類を行う疾患判定システムなどについては報告がある(特許文献5)。
(脳活動に基づくバイオマーカー)
一方で、いわゆる画像診断の分野においては、上述したような「生物学的な分子」というバイオマーカーの概念とは異なり、「画像バイオマーカー」と呼ばれるものも存在する。たとえば、脳神経領域での分子イメージングに、PET(positron emission tomography)を用いて、神経伝達機能や受容体機能解析を行おうとする試みもある。
さらに、核磁気共鳴映像法(MRI:Magnetic Resonance Imaging)では、血流量の変化に応じて、検出される信号に変化が現れることを用いて、外部刺激等に対する脳の活動部位を視覚化することも可能である。このような核磁気共鳴映像法を、特に、fMRI(functional MRI)と呼ぶ。
fMRIでは、装置としては通常のMRI装置に、さらに、fMRI計測に必要なハードおよびソフトを装備したものが使用される。
ここで、血流量の変化がNMR信号強度に変化をもたらすのは、血液中の酸素化および脱酸素化ヘモグロビンは磁気的な性質が異なることを利用している。酸素化ヘモグロビンは反磁性体の性質があり、周りに存在する水の水素原子の緩和時間に影響を与えないのに対し、脱酸素化ヘモグロビンは常磁性体であり、周囲の磁場を変化させる。したがって、脳が刺激を受け、局部血流が増大し、脱酸素化ヘモグロビンが変化すると、その変化分をMRI信号として検出する事ができる。このような被験者への刺激は、たとえば、視覚による刺激や聴覚による刺激、あるいは所定の課題(タスク)の実行等が用いられることが一般的である。
そして、脳機能研究においては、微小静脈や毛細血管における赤血球中の脱酸素化ヘモグロビンの濃度が減少する現象(BOLD効果)に対応した水素原子の核磁気共鳴信号(MRI信号)の上昇を測定することによって脳の活動の測定が行われている。
このように、fMRI装置により計測される脳活動を反映した血中酸素濃度依存信号をBOLD信号(Blood Oxygen Level Dependent Signal)と呼ぶ。
特に、人の運動機能に関する研究では、被験者に何らかの運動を行わせつつ、上記fMRI計測によって脳の活動を測定することが行われている。
ところで、ヒトの場合、非侵襲的な脳活動の計測が必要であり、この場合、fMRIデータから、より詳細な情報を抽出できるデコーディング技術が発達してきている。特に、fMRIが脳におけるボクセル単位(volumetric pixel : voxel)で脳活動を解析することで、脳活動の空間的パターンから、刺激入力や認識状態を推定することが可能となっている。
さらに、このようなデコーディング技術を発展させた技術として、特許文献6には、神経・精神疾患に対して、脳機能画像法による「診断バイオマーカー」を実現するための脳活動解析方法が開示されている。この方法では、健常群、患者群において測定された安静時機能結合的MRIのデータから、それぞれの被験者について、所定の脳領域間の活動度の相関行列(脳機能結合パラメータ)を導出する。被験者の疾患/健常ラベルを含む被験者の属性と相関行列とについて正則化正準相関解析により特徴抽出が行われる。正則化正準相関解析の結果に基づいて、スパースロジスティック回帰(Sparse Logistic Regression: SLR)による判別分析によりバイオマーカーとして機能する判別器が生成される。このような機械学習の技術により、安静時のfMRI データから導き出される脳領野間の結合にもとづいて神経疾患の診断結果を予測可能であることが示された。しかも、その予測性能の検証は1つの施設において計測された脳活動のみではなく、他の施設で計測された脳活動に対しても、ある程度の汎化が可能であることが示された。
さらに、このような「診断バイオマーカー」について、汎化性能を一層向上させるための技術改良もなされている(特許文献7)。
また、最近では、米国のヒューマンコネクトームプロジェクトのように、大規模な脳画像データを得て共有することは、基礎的な神経科学研究と、精神疾患の診断および治療のような臨床応用の間のギャップを埋めることにとって重要な意味を持つ、と認識されている (非特許文献2)。
2013年には、日本の国立研究開発法人 日本医療研究開発機構が、8つの研究所が2,239のサンプルおよび5つの疾患を含む複数サイトの安静時の機能的磁気共鳴(安静時の機能MRI)データを収集し、公にSRPBS(Strategic Research Program for Brain Sciences,https://www.amed.go.jp/program/list/01/04/001_nopro.html)の複数サイトの複数疾患のデータベース(https://bicr-resource.atr.jp/decnefpro/)を通じてそれらを共有するという、デコーディッドニューロフィードバック(DecNef)プロジェクトを組織した。このプロジェクトは、完全に独立したコホートに汎化することができる、いくつかの精神疾患の、安静時の機能的なコネクティビティ(安静時の機能結合MRI)に基づくバイオマーカーを同定している。
このように健常群と疾患群の診断については一定の成果が得られつつある。ところで、疾患群の中において、たとえば、一般に「うつ病」と診断されている患者群も、実は、複数のサブタイプに分かれていることが知られている。たとえば、通常の「抗うつ薬」の投与で寛解に向かう患者群がある一方で、寛解しにくい「治療抵抗性」の患者群があることなどが知られている。
このような「うつ病」の患者に対して、上述したような「脳機能結合パラメータ」に対して、データ駆動型人工知能によるクラスタリングを適用して、分類しようとする試みも存在し、一定の傾向が存在することを示した文献も存在する(非特許文献3、4)。
ただし、このような疾患群のサブタイプを分類する手法を実用化するためには、当該疾患群について、大規模データが必要となる。ところが、大規模に脳画像データを収集することは、健常者についても、そして、特に、患者に対しては容易ではない。
そのため、大規模なデータ収集をするために、複数のサイトでの計測を実施すると、各計測サイトでの計測データのサイト間差が問題となる。上述した非特許文献4でも、多施設の大量の計測データに対するクラスタリングの「汎化(generalization)」は、将来の課題であると言及されている。
たとえば、上述した非特許文献3では、うつ患者が、4つのサブタイプに層別化され、TMS(transcranial magnetic stimulation)に対する治療反応性の違いがあることなどを指摘された。ところが、脳機能結合指標を発見する過程で、うつ症状データを2度使用しており過学習のために、うつ症状との関連性に統計学的有意性が確認できず層別化の安定性も悪いことが、別の文献で指摘されている(非特許文献5)。
したがって、たとえば、うつ病に関しては、独立した検証データにおける層別化の精度確認は未実施というのが現状である。
一方で、たとえば、複数の計測サイトで、MRI計測をした場合の計測データのサイト間差を評価するために、複数のサイトへ多数の参加者が移動しつつ測定を受けるという、いわゆる「トラベリングサブジェクト(旅行被験者)」を採用することにより、安静時の機能的なコネクティビティに対する測定バイアスの効果を調査するという試みも実施されている(非特許文献6、非特許文献7)。
いずれにしても、fMRIデータから被験者の属性の分類を行う場合には、機械学習では過学習の問題を避けるために、一人の被験者を除いて検証用として用いる交差検証法:leave-one-subject-out cross validationやデータを10分割し、10分の9で学習し、残りの10分の1で検証を行う10-fold cross validationを用いて分類器の評価を行うことが多い。しかし、単一施設から得られた少数のサンプルに対して機械学習を適用すると、予測のインフレーションを起こす危険が、精神医学分野でも近年認識されるようになってきている。
少数のデータに対する機械学習では、学習用データにおける特定の施設のfMRI装置や測定方法、実験者、参加者群などに存在する特定の傾向、あるいはノイズに対して過学習してしまう可能性が高い。
例えば、脳の解剖画像から自閉スペクトラム症を判別する分類器は、開発に使われた英国の学習用データには感度も特異度も9割以上の高性能を示すが、日本人のデータでは5割になってしまうことが報告されている例もある。このことから、学習用データとは全く異なる施設と被験者群からなる独立検証コホートで検証していない分類器は、科学的にも実用的にも意義が少ないといえる。
上述したような計測サイト間のサイト間差を補償するための「ハーモナイゼーション法」については、本件の出願人による報告もある(非特許文献8)。
再表2018/147193号公報(国際公開WO2018/147193) 特開2019-198376号公報 特開2020-024139号公報 特表2019-516950号公報(国際公開WO2017/162773) 再表2005/025421号公報(国際公開WO2005/025421) 特開2015-62817号公報 特開2017-196523号公報
国立研究開発法人 日本医療研究開発機構 平成30年12月10日プレスリリース 「AIを搭載した内視鏡診断支援プログラムが承認―医師の診断補助に活用へ―」 https://www.amed.go.jp/news/release_20181210.html Glasser MF, et al. The Human Connectome Project's neuroimaging approach. Nat Neurosci 19, 1175-1187 (2016). Andrew T Drysdale, Logan Grosenick, Jonathan Downar, Katharine Dunlop, Farrokh Mansouri, Yue Meng1, Robert N Fetcho, Benjamin Zebley, Desmond J Oathes, Amit Etkin, Alan F Schatzberg, Keith Sudheimer, Jennifer Keller, Helen S Mayberg, Faith M Gunning, George S Alexopoulos, Michael D Fox, Alvaro Pascual-Leone, Henning U Voss, BJ Casey, Marc J Dubin & Conor Liston, "Resting-state connectivity biomarkers define neurophysiological subtypes of depression", nature medicine, VOLUME 23, NUMBER 1, JANUARY 2017 Tomoki Tokuda, Junichiro Yoshimoto,, Yu Shimizu, Go Okada, Masahiro Takamura, Yasumasa Okamoto, Shigeto Yamawaki, Kenji Doya, "Identification of depression subtypes and relevant brain regions using a data-driven approach", SCIENTIFIC REPORTS | (2018) 8:14082 | DOI:10.1038/s41598-018-32521-z Richard Dinga, Lianne Schmaal, Brenda W.J.H.Penninx, Marie Josevan Tol, Dick J.Veltman, Laura van Velzen, Maarten Mennes, Nic J.A.van der Wee, Andre F.Marquand, "Evaluating the evidence for biotypes of depression: Methodological replication and extension of Drysdale et al. (2017)", NeuroImage: Clinical 22 (2019) 101796 Noble S, et al. Multisite reliability of MR-based functional connectivity. Neuroimage 146, 959-970 (2017). Pearlson G. Multisite collaborations and large databases in psychiatric neuroimaging advantages, problems, and challenges. Schizophr Bull 35, 1-2 (2009). Ayumu Yamashita, Noriaki Yahata, Takashi Itahashi, Giuseppe Lisi, Takashi Yamada, Naho Ichikawa, Masahiro Takamura, Yujiro Yoshihara, Akira Kunimatsu, Naohiro Okada, Hirotaka Yamagata, Koji Matsuo, Ryuichiro Hashimoto, Go Okada, Yuki Sakai, Jun Morimoto, Jin Narumoto, Yasuhiro Shimada, Kiyoto Kasai, Nobumasa Kato, Hidehiko Takahashi, Yasumasa Okamoto, Saori C Tanaka, Mitsuo Kawato, Okito Yamashita, and Hiroshi Imamizu,"Harmonization of resting-state functional MRI data across multiple imaging sites via the separation of site differences into sampling bias and measurement bias.",PLOS Biology. DOI: 10.1371/journal.pbio.3000042, http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.3000042
以上のように、機能的磁気共鳴画像法などの脳機能画像法による脳活動の解析の神経・精神疾患の治療への応用を考えた場合、たとえば、上述したようなバイオマーカーとして、脳機能画像法による脳活動の解析は、非侵襲的な機能マーカーとして、診断法の開発、根本治療を実現するための創薬に向けた標的分子の探索・同定などへの応用も期待される。
たとえば、これまで、精神疾患に対しては、遺伝子を用いた実用的なバイオマーカーの完成には至っておらず、それゆえに、薬物の効果判定などが困難であるため、治療薬の開発も困難であった。
この発明は、上記のような問題点を解決するためになされたものであって、その目的は、脳活動の計測データに基づいて、機械学習により診断マーカーとしての判別器(識別器)や層別化マーカーとしての分類器を生成し、これをバイオマーカーとして使用し、うつ症状を呈する被験者の脳活動の計測結果に基づいて、前記被験者に対する治療法の選択に関する情報を提供するための治療法選択支援システム、治療法選択支援装置、治療法選択支援方法、治療法選択支援プログラムを提供する。
この発明の他の目的は、うつ症状に対する治療薬候補物質の臨床試験において、被験者の脳活動の計測結果に基づいて、被験者に対するスクリーニングを支援するためのスクリーニング支援システム、スクリーニング支援装置、スクリーニング支援方法、スクリーニング支援プログラムを提供する。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援システムに関する。治療法選択支援システムは、複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング装置を備え、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む。前記クラスタリング装置は、前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含む。前記演算装置は、クラスタリング分類器の生成処理において、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成する。前記治療法選択支援システムは、さらに、前記クラスタリング分類器による層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置と、前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果に応じて、対応する治療法情報を出力する支援情報提供装置とを、備える。
好ましくは、前記演算装置は、前記識別器モデルを生成する機械学習において、前記第1の群と前記第2の群から、アンダーサンプリングおよびサブサンプリングを実行して、複数の学習用サブサンプルを生成し、前記学習用サブサンプルのそれぞれについて、機械学習による識別器の生成において使用される特徴量の和集合から、前記和集合に属する特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、選択された前記クラスタリングのための特徴量に基づいて、前記多重共クラスタリング法により、前記クラスタリング分類器を生成する。
好ましくは、前記支援情報提供装置は、クラスタリング演算装置とインタフェース装置とを備え、前記クラスタリング演算装置は、前記クラスタリング分類器により前記第1の被験者が各前記クラスタに属する確率を算出し、前記データベース装置から、前記確率に応じて選択される少なくとも2つの前記治療法情報を読み出し、前記インタフェース装置は、前記選択されたクラスタとそれぞれ対応する前記治療法情報とを関連付けて表示するためのデータを出力する。
好ましくは、前記治療法情報は、特定の治療薬に対する応答性を示す情報である。
好ましくは、前記治療法情報は、特定の物理的治療法に対する応答性を示す情報である。
好ましくは、前記機械学習による識別器の生成の処理は、前記複数の学習用サブサンプルに対してそれぞれ複数の識別器サブモデルを生成し、前記複数の識別器サブモデルを統合して前記識別器モデルを生成するアンサンブル学習である。
好ましくは、前記クラスタリング装置は、複数の計測サイトにそれぞれ設けられた複数の脳活動計測装置から、各前記複数の第2の被験者の所定の複数の脳領野ペア間の脳活動の時間相関を表現する情報を受ける。
好ましくは、前記演算装置は、各前記複数の第2の被験者について前記複数の脳機能結合相関値について、前記計測サイトの測定バイアスを除去するように補正することで、補正された調整値を前記特徴量として前記記憶装置に格納するハーモナイゼーション算出手段を含む。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援装置に関する。治療法選択支援装置は、複数の第2の被験者のうち、うつ病の診断ラベルを有する被験者に対する層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置と、前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に基づく層別化の結果に応じて、対応する治療法情報を出力する支援情報提供装置とを備える。前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む。前記層別化の結果のクラスタは、クラスタリング装置による脳機能結合相関値の計測結果に対するクラスタリング処理により得られたクラスタリング分類器により得られるものである。前記クラスタリング装置は、前記第1の群について、前記クラスタリングの処理を実行するための演算装置と記憶装置とを含む。前記クラスタリング分類器の生成処理において、前記演算装置は、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングすることにより、前記クラスタリング分類器を生成する。
好ましくは、前記演算装置は、前記識別器モデルを生成する機械学習において、前記第1の対照者群と前記第2の対照者群から、アンダーサンプリングおよびサブサンプリングを実行して、複数の学習用サブサンプルを生成し、前記学習用サブサンプルのそれぞれについて、機械学習による識別器の生成において使用される特徴量の和集合から、前記和集合に属する特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、選択された前記クラスタリングのための特徴量に基づいて、前記多重共クラスタリング法により、前記クラスタリング分類器を生成する。
好ましくは、前記支援情報提供装置は、クラスタリング演算装置とインタフェース装置とを備える。前記クラスタリング演算装置は、前記クラスタリング分類器により前記第1の被験者が各前記クラスタに属する確率を算出し、前記データベース装置から、前記確率に応じて選択される少なくとも2つの前記治療法情報を読み出す。前記インタフェース装置は、前記選択されたクラスタとそれぞれ対応する前記治療法情報とを関連付けて表示するためのデータを出力する。
好ましくは、前記治療法情報は、特定の治療薬に対する応答性を示す情報である。
好ましくは、前記治療法情報は、特定の物理的治療法に対する応答性を示す情報である。
好ましくは、前記機械学習による識別器の生成の処理は、前記複数の学習用サブサンプルに対してそれぞれ複数の識別器サブモデルを生成し、前記複数の識別器サブモデルを統合して前記識別器モデルを生成するアンサンブル学習である。
好ましくは、前記クラスタリング装置は、複数の計測サイトにそれぞれ設けられた複数の脳活動計測装置から、各前記複数の被験者の所定の複数の脳領野ペア間の脳活動の時間相関を表現する情報を受け、前記演算装置は、各前記複数の被験者について前記複数の脳機能結合相関値について、前記計測サイトの測定バイアスを除去するように補正することで、補正された調整値を前記特徴量として前記記憶装置に格納するハーモナイゼーション算出手段を含む。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援方法に関する。治療法選択支援方法は、複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング分類器を生成して準備するステップを備え、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む。前記準備するステップは、前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算ステップを含む。前記演算ステップは、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成するステップとを有する。前記治療法選択支援方法は、さらに、前記第1の被験者の脳活動の計測結果に対する前記クラスタリング分類器による分類結果に応じて、前記クラスタリング分類器による層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベースから、対応する治療法情報を取得し出力する支援情報提供ステップを、備える。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援方法に関する。治療法選択支援方法は、前記第1の被験者の脳活動の計測結果に基づく層別化の結果のクラスタに応じて、複数の第2の被験者のうち、うつ病の診断ラベルを有する被験者に対する層別化の結果と対応する所定の治療法情報とを関連付けて格納するためのデータベースから、対応する治療法情報を取得し出力する支援情報提供ステップを備える。前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む。前記層別化の結果のクラスタは、脳機能結合相関値の計測結果に対するクラスタリング処理により得られたクラスタリング分類器により得られるものである。前記クラスタリング分類器は、前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算ステップにより生成されたものである。前記演算ステップは、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成する。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援プログラムに関する。前記治療法選択支援プログラムは、コンピュータに実行させたときに、コンピュータに、複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング分類器を生成するステップと、前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果に応じて、前記クラスタリング分類器による層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置から、対応する治療法情報を取得し出力するステップとを実行させる。前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む。前記クラスタリング処理は、前記複数の第2の被験者について、前記クラスタリングの処理を実行するための演算ステップを含む。記演算ステップは、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を記憶装置に格納するステップと、ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成するステップとを含む。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援プログラムに関する。前記治療法選択支援プログラムは、コンピュータに実行させたときに、コンピュータに、前記第1の被験者の脳活動の計測結果に基づく層別化の結果のクラスタに応じて、複数の第2の被験者のうち、うつ病の診断ラベルを有する被験者に対する層別化の結果と対応する所定の治療法情報とを関連付けて格納するためのデータベースから、対応する治療法情報を取得し出力する支援情報提供ステップを実行させる。前記層別化の結果のクラスタは、脳機能結合相関値の計測結果に対するクラスタリング処理により得られたクラスタリング分類器により得られるものである。前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む。前記クラスタリング分類器は、前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算ステップにより生成されたものである。前記演算ステップは、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成するステップとを含む。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援システムに関する。スクリーニング支援システムは、複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング装置を備える。前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む。前記クラスタリング装置は、前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含む。前記演算装置は、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成する。前記スクリーニング支援システムは、さらに、前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力する支援情報提供装置とを、備える。
好ましくは、前記演算装置は、前記識別器モデルを生成する機械学習において、前記第1の群と前記第2の群から、アンダーサンプリングおよびサブサンプリングを実行して、複数の学習用サブサンプルを生成し、前記学習用サブサンプルのそれぞれについて、機械学習による識別器の生成において使用される特徴量の和集合から、前記和集合に属する特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、選択された前記クラスタリングのための特徴量に基づいて、前記多重共クラスタリング法により、前記クラスタリング分類器を生成する。
好ましくは、前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援装置に関する。スクリーニング支援装置は、クラスタリング分類器を特定する情報を格納するための記憶装置を有する支援情報提供装置であって、前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器に基づく分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力する支援情報提供装置を備える。前記層別化の結果のクラスタは、クラスタリング装置による脳機能結合相関値の計測結果に対するクラスタリング処理により得られた前記クラスタリング分類器により得られる。前記クラスタリング装置は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む複数の第2の被験者について、前記クラスタリングの処理を実行するための演算装置と記憶装置とを含む。前記クラスタリング分類器の生成処理において、前記演算装置は、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成する。
好ましくは、前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援方法に関する。スクリーニング支援方法は、記憶装置に格納された情報で特定されるクラスタリング分類器に基づいて、前記脳活動の計測結果により、演算装置が前記第1の被験者の分類を実行するステップと、前記分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力するステップと、を備える。前記クラスタリング分類器の生成のための処理は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む複数の第2の被験者について、前記クラスタリング処理を実行するための演算ステップを含む。前記演算ステップは、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成するステップとを含む。
好ましくは、前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である。
この発明の1つの局面に従うと、本発明のある実施形態は、うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援プログラムに関する。前記スクリーニング支援プログラムは、コンピュータに実行させたときに、コンピュータに、記憶装置に格納された情報で特定されるクラスタリング分類器に基づいて、前記脳活動の計測結果により、演算装置が前記第1の被験者の分類を実行するステップと、前記分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力するステップと、を実行させる。前記クラスタリング分類器の生成のための処理は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む複数の第2の被験者について、前記クラスタリングの処理を実行するための演算ステップを含み、前記演算ステップにおいて、i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を記憶装置に格納するステップと、ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成するステップを含む。
好ましくは、前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である。
機械学習により生成された、診断マーカーとしての判別器や層別化マーカーとしての分類器を用いて、治療を必要とするうつ病の被験者に対する治療方法の選択を支援することができる。
また、診断マーカーとしての判別器や層別化マーカーとしての分類器を用いてうつ症状に対する治療薬候補物質の臨床試験に参加する被験者のスクリーニングを支援することが可能となる。
複数の計測サイトに設置されたMRI計測システムにより計測されたデータについて、ハーモナイゼーション処理を説明するための概念図である。 被験者の脳の関心領域(ROI:Region of Interest)について、安静時の機能結合の相関を示す相関行列を抽出する手続きを示す概念図である。 「測定パラメータ」および「被験者属性データ」の内容の例を示す概念図である。 各計測サイトに設置されるMRI装置100.i(1≦i≦Ns)の全体構成を示す模式図である。 データ処理部32のハードウェアブロック図である。 相関行列から、診断マーカーとなる判別器を生成する過程およびクラスタリング処理について説明する概念図である。 計算処理システム300の構成を説明するための機能ブロック図である。 計算処理システム300の構成を説明するための機能ブロック図である。 アンサンブル学習による疾患識別器を生成するための機械学習の手続きを説明するためのフローチャートである。 学習用のデータセット(データセット1)のデモグラフィック特性を示す図である。 独立検証のデータセット(データセット2)のデモグラフィック特性を示す図である。 全撮像サイトについて、学習用データセットに対するMDDの予測性能(出力の確率分布)を示す図である。 各撮像サイトについて、学習用データセットに対するMDDの予測性能(識別器の出力の確率分布)を示す図である。 独立した検証データセットにおけるMDDの識別器の出力の確率分布を示す図である。 各撮像サイトについて、独立検証データセットに対するMDDの識別器の出力の確率分布を示す図である。 特徴量を選択して、教師なし学習によりクラスタリングを行う処理を説明するためのフローチャートである。 「特徴量選択を伴う学習処理」により、複数個(たとえば、Nch個)の特徴量が存在する場合に、特徴量の選択が実施される概念を示す図である。 特徴量選択を伴う学習処理により1つの識別器を生成する際に最終的に選択される特徴量を示す概念図である。 アンダーサンプリングおよびサブサンプリングの処理を複数回実施して、識別器を生成される際に、特徴量が選択される様子を示す概念図である。 特徴量に依存して複数のクラスタの分け方が存在する場合を説明するための概念図である。 複数の特徴量により、複数の対象が特徴づけられる場合のクラスタリングの概念を説明するための概念図である。 多重クラスタリングおよび多重共クラスタリングを説明するための概念図である。 「多重共クラスタリング」において、1つのビューの中に、異なる種類の確率分布の確率モデルが想定されている場合を示す概念図である。 多重共クラスタリングの学習方法の概要を説明するためのフローチャートである。 多重共クラスタリングの学習方法におけるベイズ推定のグラフ表現を示す図である。 2つに分割されたデータセット1とデータセット2とを示す図である。 各データセットにおいてクラスタリングを実施する概念を説明する概念図である。 被験者データに対する多重共クラスタリングの例を示す概念図である。 データセット1およびデータセット2に対して、実際に多重共クラスタリング処理を実施した結果を示す図である。 データセット1とデータセット2において、それぞれのビューに割り当てられた脳機能結合(FC)の個数を示す表である。 クラスタリングの類似度(層別化の汎化性能)の評価手法を説明するための概念図である。 ARIについて説明するための概念図である。 クラスタリング1とクラスタリング1´の間、および、クラスタリング2とクラスタリング2´の間の類似度の評価結果を示す図である。 クラスタリング1とクラスタリング1´において、それぞれビュー1の各クラスタに割り当てられる被験者数の分布を示す表である。 サイト間を移動しながら計測を受ける移動被験者によるサイト間差の評価手法を説明するための概念図である。 被験者aのb番目の機能的なコネクティビティの表現を説明するための概念図である。 ハーモナイゼーションのために、測定バイアスを算出する処理を説明するためのフローチャートである。 データ収集、推定処理および対象者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。 治療法選択支援システム1000の機能構成を示す図である。 治療法情報データベース5100の例を示す図である。 被験者の治療法選択の支援処理のフローを説明するフローチャートである。 一般的な創薬のプロセスを示す図である。 被験者のスクリーニング支援処理のフローを説明するフローチャートである。 スクリーニング支援装置1000´の構成を説明する図である。 データセット1について多重共クラスタリングによるクラスタリングの結果のビューを示す図である。 データセット2について多重共クラスタリングによるクラスタリングの結果のビューを示す図である。 データセット1とデータセット2との間のクラスタリング安定性を示す図である。 全うつ病患者データ(データセット1+2)に対するクラスタリング結果のビューを示す図である。 データセット1+2のクラスタリングで生成されたビュー3について、全うつ病患者の数と、臨床データが存在するうつ病患者の数とを、サブタイプごとに示す図である。 クラスタリングに使用された脳機能結合(ビュー3)を示す図である。 ビュー3の各サブタイプついて、うつ病の重症度及び重症度の改善率との関係を示す図である。
I.学習フェーズ
以下、本発明の「脳機能結合相関値のクラスタリング装置」、「脳機能結合相関値のクラスタリング方法」等を説明するために、複数の脳活動計測装置からなる計測システムで計測した被験者(精神疾患の患者を含む)の脳機能結合画像データに対する人工知能技術による「クラスタリング」を例として説明することとする。
そこで、以下、本発明の実施の形態の計測システム、より具体的にはMRI計測システムの構成について、図に従って説明する。なお、以下の実施の形態において、同じ符号を付した構成要素および処理工程は、同一または相当するものであり、必要でない場合は、その説明は繰り返さない。
また、本実施の形態では、複数の施設に設置された「脳活動計測装置」、より具体的には、「MRI装置」によって、脳の複数の領野間の脳活動を時系列で計測し、これらの領野間の時間相関(「脳機能結合」と呼ぶ)のパターンに基づいて、複数施設に汎化可能に、特定の疾患を有する被験者を、さらに複数のグループ(サブグループ)に分類するものとして本発明を説明する。
なお、これも特に限定されないが、「特定の疾患」については、「大うつ病」を例として説明することとする。ただし、以下の説明にも示されるように、本発明は、被験者の「脳機能結合相関値」をデータ駆動による分類する技術に関するものであって、被験者の疾患は、「大うつ病」に限定されず、他の疾患であってもよい。また、被験者の「脳機能結合相関値」のパターンにより分類される被験者の属性であれば、必ずしも、疾患である必要もなく、他の属性であってもよい。
そして、このような「MRI計測システム」では、複数の「MRI装置」が、複数の異なる施設に設置されており、後述するように、これらの計測施設(計測サイト)間での計測に対するサイト間差を、測定機器を原因とする測定バイアスと、計測サイトにおける被験者のポピュレーションによる差(標本バイアス)とを独立に評価する。その上で、各計測サイトにおける計測値について、測定バイアスの効果を除いてサイト間差を補正する処理を実施することで、計測サイト間での測定結果に対する調和処理(ハーモナイゼーション)を実現する。そして、ハーモナイゼーション後の脳機能結合値について、特定疾患の診断ラベルを教師データとするアンサンブル学習を用いて、「特徴量選択」を行った後に、教師なし学習によるクラスタリングを実行することで、被験者属性(たとえば、精神疾患のサブタイプ)の分類を実行するものとして、説明する。
[実施の形態1]
図1は、複数の計測サイトに設置されたMRI計測システムにより計測されたデータについて、クラスタリング(層別化)処理を説明するための概念図である。
図1を参照して、計測サイトMS.1~MS.Ns(Ns:サイト数)には、それぞれ、MRI装置100.1~100.Nsが設置されているものとする。
また、計測サイトMS.1~MS.Nsにおいては、それぞれ、被験者群PA.1~PA.Nsの計測が行われる。被験者群PA.1~PA.Nは、本明細書において、第2の被験者群とも呼ばれる。また、第2の被験者群に属する各人は第2の被験者とも呼ばれる。すなわち、第2の被験者群は、複数の第2の被験者を含む。被験者群PA.1~PA.Nsの各々には、少なくとも2つ以上に分類される群、たとえば、患者群と健常者群とが含まれるものとする。患者群を、患者群ではない健常者群と特に区別して呼ぶ必要がある場合、本明細書において、第1の群と呼ぶ。第1の群に属する被験者は、たとえば、DSM-5などの診断法による診断に基づき、うつ病の診断ラベルを有する被験者である。健常者群を、患者群と特に区別して呼ぶ必要がある場合、本明細書において、第2の群と呼ぶ。第2の群に属する被験者は、うつ病の診断ラベルを有さない被験者である。また、患者群としては、特に限定されないが、たとえば、精神疾患の患者、より特定的には、「大うつ病の患者」の群が相当するものとする。なお、「診断ラベル」の決定方法としては、上述したような従来の「DSM-5などの症候に基づく診断法」に限定されるものではなく、たとえば、後に例示するように、脳活動の計測データの分析結果を補助情報として決定される方法であってもよい。
そして、各計測サイトにおいては、原則として、MRI装置の仕様上、可能な範囲で、統一した計測プロトコルで各被験者の計測が実施されるものとする。
ここで、特に限定されないが、計測プロトコルとしては、たとえば、以下のような内容が規定されているものとする。
1)頭部のスキャンを実行する方向 たとえば、頭部の後側(posterior:以下、“P”と略記する)から前側(anterior:以下、“A”と略記する)に向かう方向(以下、「P→A方向」と呼ぶ)と、逆の方向、すなわち、前側から後側に向かう方向(以下、「A→P方向」と呼ぶ)にスキャンを行うのかを規定する必要がある。状況によっては、両方の場合のスキャンを行うことを規定することもあり得る。
MRI装置によって、デフォルトの方向が違っていたり、あるいは、両方を任意に設定できない場合もあり得る。
スキャンの方向は、たとえば、画像としての「ひずみ方」を規定する可能性があり、プロトコルとして、条件を設定する。
2)脳構造画像の撮像条件
いわゆるスピンエコー法により、「T1強調画像」や「T2強調画像」のいずれか、または両方を撮像する条件を設定する。
3)脳機能画像の撮像条件
fMRI(functional Magnetic Resonance Imaging)法により、「安静時」の被験者の脳機能画像を撮像する条件を設定する。
4)拡散強調画像の撮像条件
拡散強調画像(DWI:diffusion (weighted) image)を撮像するか否か、そしてその条件について設定する。
拡散強調画像とは、MRI撮像のシーケンスの一種で、水分子の拡散運動を画像化したものである。通常用いられるスピンエコー法のパルス系列では拡散による信号の減衰は無視できるが、大きな傾斜磁場が長時間にわたって印加されると、その間の各磁化ベクトルの移動によって生じる位相のずれが無視できなくなり、拡散が活発な領域ほど低信号として現れる、ことを利用するものである。
5)EPI歪みを画像処理で補正するための撮像
たとえば、EPI歪みを画像処理で補正するための1つの方法として、「フィールドマップ法」が知られており、空間歪みの補正に対する撮像の条件を設定する。
フィールドマップ法は、多重エコー時間によってEPI画像を収集し、これらのEPI画像に基づいてEPI歪みの量を算出する。フィールドマップ法を適用して、新規の画像に含まれるEPI歪みを補正することができる。異なるエコー時間による同一の解剖学的構造の一組の画像を前提として、EPI歪みを計算して、画像の歪みを補正することが可能である。
たとえば、以下の公知文献に、「フィールドマップ法」が開示されている。
公知文献1 : 特開2015-112474号明細書
なお、計測プロトコルには、上記の条件の中から適宜、必要となるシーケンス部分を抜き出してもよいし、必要に応じて、他のシーケンスやその条件を追加してもよい。
再び、図1を参照して、各計測サイトMS.1~MS.Nsにおいて、計測を行う被験者として採用することを、「被験者をサンプリングする」と呼び、各計測サイトでのサンプリングの偏りにより生じる、計測値のサイト間差の原因を、「標本バイアス」と呼ぶ。
たとえば、上述した例では、従来の診断基準で、「大うつ病」と診断される患者には、実は、いくつかのサブタイプが含まれることが知られている。
典型的には、「メランコリー型」、「非定型」、「季節型」、「産後」うつ病などがある。また、「2種類以上の作用機序の異なる抗うつ薬による十分な治療により十分に改善せず、中等症以上の症状が続く場合」は、「治療抵抗性うつ病」と呼ばれ、うつ病の10~20%と推定されるとの報告もある。すなわち、一般に、「大うつ病」と診断される患者群は、決して、均質なものではないことが知られている。ただし、このようなサブタイプを、客観的な計測データに基づいて、分類する手法は、これまでは、必ずしも実用化に至っていない。
各計測サイトでは、当該計測サイトの病院に来院する患者の地域性に由来する性向の偏りや、当該病院での診断の傾向等々のさまざまの要因により、一口に、「大うつ病」の患者といったとしても、その患者群に含まれるサブタイプの分布が均一であるとは、必ずしも言えない。この結果、各計測サイトの患者群には、そのサブタイプの分布に、偏りが生じることが、むしろ一般的であり、その結果、上述した「標本バイアス」が生じると考えられる。
また、「健常者群」と呼ばれる被験者の群であっても、その中には、複数のサブタイプが存在することが一般であり、その点では、「健常者群」においても、「標本バイアス」が存在することになる。
また、MRI装置100.1~100.Nsも、各計測サイトで決して同一計測特性を有するMRI装置が使用されているとは言えない。
たとえば、MRI装置の製造メーカー、MRI装置の型番、MRI装置の静磁場強度、MRI装置における(送)受信コイルのコイル数(チャンネル数)、等々、MRI装置の条件、MRI装置の計測条件などの測定条件に応じて、計測データのサイト間差異が生じうる。このような測定条件により生じるサイト間差異のことを「測定バイアス」と呼ぶ。
仮に、同一のMRI装置の製造メーカーの同一の型番のMRI装置であっても、装置の固有な個性により、完全に同一の計測特性を実現するとも、必ずしも言えない。
ここで、(送)受信コイルは、計測される信号のSN比を向上させることを目的として、「マルチアレイコイル」が採用されることが一般的である。「受信コイルのコイル数」とは、マルチアレイコイルを構成する「エレメントコイル」の個数を指す。個々のエレメントコイルの感度を高め,この出力を束ねることで受信感度の向上を図るものである。
そして、特に限定されないが、後に説明するようなハーモナイゼーション手法により、本実施の形態では、「標本バイアス」と「測定バイアス」とを独立に評価することを可能とする。
再び、図1を参照して、各計測サイトMS.1~MS.Nsのそれぞれからの計測関連データDA100.1~DA100.Nsは、データセンター200内の記憶装置210に集積され、格納される。
ここで、「計測関連データ」には、各計測サイトでの「測定パラメータ」と、各計測サイトにおいて計測された「患者群データ」および「健常者群データ」とが含まれる。
さらに、「患者群データ」および「健常者群データ」には、それぞれ、各々の被験者に対応して「患者のMRI計測データ」および「健常者のMRI計測データ」とが含まれる。
以下では、このような「計測関連データ」について説明する。
図2は、被験者の脳の関心領域(ROI:Region of Interest)について、安静時の機能結合の相関を示す相関行列を抽出する手続きを示す概念図である。
ここで、図1において、「患者群データ」および「健常者群データ」における、「患者のMRI計測データ」および「健常者のMRI計測データ」には、少なくとも、以下のようなデータが含まれる。
i)相関行列のデータを算出するための時系列の「脳機能画像データ」、および/または、相関行列のデータ自身
すなわち、図1における計算処理システム300が、記憶装置210内に格納されたデータに基づいて、後述するような脳活動バイオマーカーを算出する際の基礎となるデータとして、これらのデータが使用される。
ここで、相関行列のデータについては、時系列の「脳機能画像データ」に基づいて、各計測サイトで計算された後に、記憶装置210内に格納され、計算処理システム300が、この記憶装置210内の相関行列のデータに基づいて、脳活動バイオマーカを算出する、という構成とすることができる。
あるいは、時系列の「脳機能画像データ」が記憶装置210内に格納され、計算処理システム300が、この記憶装置210内の「脳機能画像データ」に基づいて、相関行例のデータを算出し、さらに、脳活動バイオマーカを算出する、という構成とすることも可能である。
したがって、「患者のMRI計測データ」および「健常者のMRI計測データ」の各々については、少なくとも、相関行列のデータを算出するための時系列の「脳機能画像データ」、および、相関行列のデータ自身のいずれか一方が含まれることになる。
ii)被験者の構造画像データ、および、拡散強調画像データ
なお、特に限定されないが、EPI歪みを画像処理で補正する処理については、各計測サイトで演算処理がされた後に、記憶装置210内にデータが格納される構成とすることができる。
また、特に限定されないが、個人情報保護の観点から、記憶装置210内へデータが格納される前に、各計測サイトにおいて、匿名加工処理を実行する構成とすることが可能である。ただし、匿名加工処理については、計算処理システム300を運営する主体が、法規制上、個人情報を扱うことが認められている場合などは、計算処理システム300において匿名加工処理を実行する構成としてもよい。
図2に戻って、図2(a)に示すように、リアルタイムで測定した安静時のfMRIのn個(n:自然数)の時刻分のfMRIデータから、各関心領域の平均的な「活動度」を算出し、図2(b)に示すように、脳領域間(関心領域間)の機能的なコネクティビティ(「活動度の相関値」)の相関行列を算出する。
(脳領域の区画化(パーセレーション:Parcellation))
機能的なコネクティビティは、各参加者に対して2つの脳領域間の安静時の機能MRIの血中酸素濃度に依存する(BOLD)信号の時間的な相関性として計算される。
ここで、関心領域としては、上述したようにNr領域を考えているので、相関行列における独立な非対角成分は、対称性を考慮すると、
Nr×(Nr-1)/2(個)
ということになる。
関心領域の設定の方法としては、以下のような方法が想定される。
方法1)“解剖学的な脳領野に基づいて、関心領域を定義する。”
ここでは、脳活動バイオマーカーのために、たとえば、関心領域として140の領域を採用する。
つまり、この方法では、ROIについては、Brain Sulci Atlas (BAL)に含まれる137のROIに加え、 Automated Anatomical Labeling Atlas の小脳(左右)及び虫部のROIを用いる。これら合計140のROI間の機能的結合FCを特徴量として用いる。
ここで、Brain Sulci Atlas (BAL)およびAutomated Anatomical Labeling Atlasについては、以下に開示がある。
公知文献2:Perrot et al., Med Image Anal, 15(4), 2011
公知文献3:Tzourio-Mazoyer et al., Neuroimage, 15(1), 2002
このような関心領域としては、たとえば、以下のような領域である。
背内側前頭前皮質(DMPFC)
前頭前野腹内側部(VMPFC)
前帯状皮質(ACC)
小脳虫部、
左視床、
右下頭頂葉、
右尾状核、
右中後頭葉、
右中帯状皮質
ただし、採用する脳の領野は、このような領域に限定されるものではない。
たとえば、対象とする神経・精神疾患に応じて、選択する領域を変更してもよい。
方法2) “機能的なコネクティビティを、全体の脳をカバーする機能的な脳地図の脳領域に基づいて、定義する。”
ここで、このような機能的な脳地図の脳領域については、以下のような文献にも開示がされており、特に限定されないが、たとえば、268のノード(脳領域)から成る、という構成とすることができる。
公知文献4:Noble S, et al. Multisite reliability of MR-based functional connectivity. Neuroimage 146,959-970 (2017).
公知文献5:Finn ES, et al. Functional connectome fingerprinting: identifying individuals using patterns of brain connectivity. Nat Neurosci 18, 1664-1671 (2015).
公知文献6:Rosenberg MD, et al. A neuromarker of sustained attention from whole-brain functional connectivity. Nat Neurosci 19, 165-171 (2016).
公知文献7:Shen X, Tokoglu F, Papademetris X, Constable RT. Groupwise whole-brain parcellation from resting-state fMRI data for network node identification. Neuroimage 82, 403-415 (2013).
方法3)表面ベースの方法
脳領域のパーセレーションについては、ヒューマンコネクトームプロジェクト(HCP)スタイルのマルチモダリティなイメージング(ミエリン・タスク・ファンクショナル)を用いることにより、脳を脳溝に沿ってシート状に変換して作成された脳地図に基づく「表面ベースの方法」で、データを分析することも可能である。
このようなパーセレーション法のためには、以下のサイトに開示されるようなツールボックスを使用することができる(ciftifyツールボックスバージョン2.0.2)。
https://edickie.github.io/ciftify/#/
このツールボックスは、使用するデータを分析することを、(たとえば、HCPパイプラインに必要なT2強調画像を欠いている場合でも)HCP類似の表面ベースのパイプラインの中で可能にするものである。
そして、方法3の解析では、関心領域(ROI)として、以下の公知文献に開示される379の表面ベースの区画(皮質の360の区画+皮質下の19の区画)を使用する。
公知文献8:Glasser, M.F., Coalson, T.S., Robinson, E.C., Hacker, C.D., Harwell, J., Yacoub, E., et al. (2016). A multi-modal parcellation of human cerebral cortex. Nature 536(7615), 171-178. doi: 10.1038/nature18933.
したがって、BOLD信号の時間的変化は、これらの379個の関心領域(ROI)から抽出される。
さらに、以下の文献に開示されるような解剖学的な自動標識付(AAL)およびNeurosynth(http://neurosynth.org/locations/)の使用により、重要なROIの解剖学的名称、およびROIを含む内在性の脳ネットワークの名称を特定することができる。
公知文献9:Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15(1), 273-289. doi: 10.1006/nimg.2001.0978.
方法4)データドリブンに脳の領域を定めた方法
下記の文献にも開示があるように、事前情報(脳地図)なしに位相のそろったボクセルからネットワークを新規に同定する手法であり、「Canonical ICA」や「辞書学習」などと呼ばれる手法である。
公知文献:Kamalaker Dadi,Mehdi Rahim,Alexandre Abraham,Darya Chyzhyk,Michael Milham,Bertrand Thirion,Gael Varoquaux,“Benchmarking functional connectome-based predictive models for resting-state fMRI”,Preprint submitted to NeuroImage,October 31,2018.
以下では、基本的に「方法3」の表面ベースの脳地図の脳領域に基づいて機能的なコネクティビティを定義する方法を用いるものとして説明する。
また、相関値の計算としては、たとえば、タンジェント法や部分相関法のような機能的なコネクティビティの計測に対するいくつかの候補がある。
しかしながら、以下では、特に限定されないが、ピアソンの相関係数を使用するものとして説明する。
前処理されたBOLD信号の可能なノード組の各々の時間経過の間についてのフィッシャーのz変換後のピアソンの相関係数が計算され、要素がそれぞれ2つのノード間のコネクティビティの強度を表わす、379×379の対照的な脳機能結合マトリックスを構築するために使用される。
図3は、「測定パラメータ」および「被験者属性データ」の内容の例を示す概念図である。
「被験者属性データ」は、図1における「患者群データ」または「健常者群データ」において、「患者のMRI計測データ」および「健常者のMRI計測データ」に、それぞれ関連づけて格納されているものとする。
図3(a)に示すように、計測サイトを識別するためのサイトIDと、サイト名と、測定パラメータを識別するための条件IDと、計測装置に関する情報と、計測条件に関する情報とを含む。
「測定パラメータ」は、「計測装置に関する情報」と「計測条件に関する情報」とを含む。
「計測装置に関する情報」には、各計測サイトにおいて被験者の脳活動を計測するためのMRI装置のメーカー名と、型番と、(送)受信コイル数とを含む。
なお、「計測装置に関する情報」は、これらのものに限られず、たとえば、静磁場強度、シム調整後の磁場の均一性、等々の他の計測装置の性能を表す指標などを含んでいてもよい。
「計測条件に関する情報」には、画像再構成時の位相エンコードの方向(P→AまたはA→P)、画像種類(T1強調、T2強調、拡散強調など)、撮像シーケンス(スピンエコーなど)、撮像中の被験者の開眼/閉眼、等々の情報が含まれる。
「計測条件に関する情報」も、これらのものに限られるものではない。
図3(b)に示されるように、「被験者属性データ」は、被験者の識別ができないように仮名制御した被験者仮ID,その被験者を計測した際の測定条件を表す条件ID、被験者の属性情報とを含む。
そして、「被験者の属性情報」としては、被験者の性別、年齢、健常または疾患のいずれかを示すラベル、医師による被験者の診断病名、被験者への投薬履歴、診断履歴などの情報が含まれる。
なお、「被験者の属性情報」は、必要に応じて、たとえば、計測サイトにおいて、匿名加工処理が行われているものとする。
たとえば、年齢や性別などについては、準識別子(同一属性)のデータがk件以上になるようにデータを変換するなどの処理を行うことで、個人が特定される確率をk分の1以下に低減し特定を困難にする「k-匿名性」を維持するように加工することができる。ここで、「準識別子」とは、「年齢」「性別」「居住地」など単体では特定できないが、組み合わせることで個人の特定が可能になる属性をいう。
また、投薬履歴や診断履歴は、日付のランダム化やシフト化(相対化)など、必要に応じて、匿名化のための処理を行う。
そして、以下では、「患者のMRI計測データ」および「健常者のMRI計測データ」で、上述した方法で、各被験者の脳領野間の各々について経時的な活動の相関として算出される機能的結合のことを、各領野間について総称して「機能的なコネクティビティ(Connectivity)」(省略する際には、“FC”と表記)と呼ぶことにする。脳領野ごとに、機能的なコネクティビティを区別する必要がある場合は、後述するように、下添え字を付加して区別するものとする。
(MRI装置の構成)
図4は、各計測サイトに設置されるMRI装置100.i(1≦i≦Ns)の全体構成を示す模式図である。
図4においては、1番目の計測サイトのMRI装置100.1を、例示的に、詳しく説明している。他のMRI装置100.2~100.Nsについても、基本的な構成は同様である。
図4に示すように、MRI装置100.1は、被験者2(第1の被験者、又は第2の被験者であり得る)の関心領域に制御された磁場を付与してRF波を照射する磁場印加機構11と、この被験者2からの応答波(NMR信号)を受信してアナログ信号を出力する受信コイル20と、この被験者2に付与される磁場を制御するとともにRF波の送受信を制御する駆動部21と、この駆動部21の制御シーケンスを設定するとともに各種データ信号を処理して画像を生成するデータ処理部32とを備える。
なお、ここで、被験者2が載置される円筒形状のボアの中心軸をZ軸にとりZ軸と直交する水平方向にX軸及び鉛直方向にY軸を定義する。
MRI装置100.1は、このような構成であるので、磁場印加機構11により印加される静磁場により、被験者2を構成する原子核の核スピンは、磁場方向(Z軸)に配向するとともに、この原子核に固有のラーモア周波数でこの磁場方向を軸とする歳差運動を行う。
そして、このラーモア周波数と同じRFパルスを照射すると、原子は共鳴しエネルギーを吸収して励起され、核磁気共鳴現象(NMR現象;Nuclear Magnetic Resonance)が生じる。この共鳴の後に、RFパルス照射を停止すると、原子はエネルギーを放出して元の定常状態に戻る緩和過程で、ラーモア周波数と同じ周波数の電磁波(NMR信号)を出力する。
この出力されたNMR信号を被験者2からの応答波として受信コイル20で受信し、データ処理部32において、被験者2の関心領域が画像化される。
磁場印加機構11は、静磁場発生コイル12と、傾斜磁場発生コイル14と、RF照射部16と、被験者2をボア中に載置する寝台18とを備える。
被験者2は、寝台18に、たとえば、仰臥する。被験者2は、特に限定されないが、たとえば、プリズムメガネ4により、Z軸に対して垂直に設置されたディスプレイ6に表示される画面を見ることができる。このディスプレイ6の画像により、必要に応じて、被験者2に視覚刺激が与えることもできる。なお、被験者2への視覚刺激は、被験者2の目前にプロジェクタにより画像が投影される構成であってもよい。
このような視覚刺激は、被験者にニューロフィードバックを行う場合は、フィードバック情報の提示に相当する。
駆動部21は、静磁場電源22と、傾斜磁場電源24と、信号送信部26と、信号受信部28と、寝台18をZ軸方向の任意位置に移動させる寝台駆動部30とを備える。
データ処理部32は、操作者(図示略)から各種操作や情報入力を受け付ける入力部40と、被験者2の関心領域に関する各種画像及び各種情報を画面表示する表示部38と、各種処理を実行させるプログラム・制御パラメータ・画像データ(構造画像等)及びその他の電子データを記憶する記憶部36と、駆動部21を駆動させる制御シーケンスを発生させるなどの各機能部の動作を制御する制御部42と、駆動部21との間で各種信号の送受信を実行するインタフェース部44と、関心領域に由来する一群のNMR信号からなるデータを収集するデータ収集部46と、このNMR信号のデータに基づいて画像を形成する画像処理部48と、ネットワークとの間で通信を実行するためのネットワークインタフェース50を備える。
また、データ処理部32は、専用コンピュータである場合の他、各機能部を動作させる機能を実行する汎用コンピュータであって、記憶部36にインストールされたプログラムに基づいて、指定された演算やデータ処理や制御シーケンスの発生をさせるものである場合も含まれる。以下では、データ処理部32は、汎用コンピュータであるものとして説明する。
静磁場発生コイル12は、Z軸周りに巻回される螺旋コイルに静磁場電源22から供給される電流を流して誘導磁場を発生させ、ボアにZ軸方向の静磁場を発生させるものである。このボアに形成される静磁場の均一性の高い領域に被験者2の関心領域を設定することになる。ここで、静磁場発生コイル12は、より詳しくは、たとえば、4個の空芯コイルから構成され、その組み合わせで内部に均一な磁界を作り、被験者2の体内の所定の原子核、より特定的には水素原子核のスピンに配向性を与える。
傾斜磁場発生コイル14は、Xコイル、Yコイル及びZコイル(図示省略)から構成され、円筒形状を示す静磁場発生コイル12の内周面に設けられる。
なお、傾斜磁場の均一性を向上させるためにシムコイル(図示せず)が設けられ、「シム調整」が実施される。
これらXコイル、Yコイル及びZコイルは、それぞれX軸方向、Y軸方向及びZ軸方向を順番に切り替えながら、ボア内の均一磁場に対し傾斜磁場を重畳させ、静磁場に強度勾配を付与する。Zコイルは励起時に、磁界強度をZ方向に傾斜させて共鳴面を限定し、Yコイルは、Z方向の磁界印加の直後に短時間の傾斜を加えて検出信号にY座標に比例した位相変調を加え(位相エンコーディング)、Xコイルは、続いてデータ採取時に傾斜を加えて、検出信号にX座標に比例した周波数変調を与える(周波数エンコーディング)。
この重畳される傾斜磁場の切り替えは、制御シーケンスに従って、Xコイル、Yコイル及びZコイルにそれぞれ異なるパルス信号が送信部24から出力されることにより実現される。これにより、NMR現象が発現する被験者2の位置を特定することができ、被験者2の画像を形成するのに必要な三次元座標上の位置情報が与えられる。
ここで、上述のように、3組の直交する傾斜磁場を用いて、それぞれにスライス方向、位相エンコード方向、および周波数エンコード方向を割り当ててその組み合わせにより様々な角度から撮影を行える。たとえば、一般にX線CT装置で撮像されるものと同じ方向のトランスバーススライスに加えて、それと直交するサジタルスライスやコロナルスライス、更には面と垂直な方向が3組の直交する傾斜磁場の軸と平行でないオブリークスライス等について撮像することができる。
RF照射部16は、制御シーケンスに従って信号送信部26から送信される高周波信号に基づいて、被験者2の関心領域にRF(Radio Frequency)パルスを照射するものである。
なお、RF照射部16は、図1において、磁場印加機構11に内蔵されているが、寝台18に設けられたり、あるいは、受信コイル20と一体化されて、送受信コイル20として構成されていてもよい。
受信コイル20は、被験者2からの応答波(NMR信号)を検出するものであって、このNMR信号を高感度で検出するために、被験者2に近接して配置されている。
ここで、受信コイル20には、NMR信号の電磁波がそのコイル素線を切ると電磁誘導に基づき微弱電流が生じる。この微弱電流は、信号受信部28において増幅され、さらにアナログ信号からデジタル信号に変換されデータ処理部32に送られる。
(送)受信コイル20については、上述のとおり、SN比の向上のためにマルチアレイコイルが使用される。
すなわち、静磁界にZ軸傾斜磁界を加えた状態にある被験者2に、共鳴周波数の高周波電磁界を、RF照射部16を通じて印加すると、磁界の強さが共鳴条件になっている部分の所定の原子核、たとえば、水素原子核が、選択的に励起されて共鳴し始める。共鳴条件に合致した部分(たとえば、被験者2の所定の厚さの断層)にある所定の原子核が励起され、(古典的な描像では)スピンがいっせいに回転する。励起パルスを止めると、受信コイル20には、今度は、回転しているスピンが放射する電磁波が信号を誘起し、しばらくの間、この信号が検出される。この信号によって、被験者2の体内の、所定の原子を含んだ組織を観察する。そして、信号の発信位置を知るために、XとYの傾斜磁界を加えて信号を検知する、という構成になっている。
画像処理部48は、記憶部36に構築されているデータに基づき、励起信号を繰り返し与えつつ検出信号を測定し、1回目のフーリエ変換計算により、共鳴の周波数をX座標に還元し、2回目のフーリエ変換でY座標を復元して画像を得て、表示部38に対応する画像を表示する。
たとえば、このようなMRIシステムにより、上述したBOLD信号をリアルタイムで撮像し、制御部42により、時系列に撮像される画像について、後に説明するような解析処理を行うことで、安静時機能結合的MRI(rs-fcMRI)の撮像を行うことが可能となる。
図4においては、MRI装置100.1および他の計測サイトにおけるMRI装置100.2~100.Nsからの計測データ、測定パラメータ、被験者属性データが、データセンター200内の通信インタフェース202を介して記憶装置210に集積され、格納される。さらに、計算処理システム300は、通信インタフェース204を介して、記憶装置210内のデータにアクセスする構成となっている。
図5は、データ処理部32のハードウェアブロック図である。
データ処理部32のハードウェアとしては、上述のとおり、特に限定されないが、汎用コンピュータを使用することが可能である。
図5において、データ処理部32のコンピュータ本体2010は、メモリドライブ2020、ディスクドライブ2030に加えて、演算装置(CPU)2040と、ディスクドライブ2030及びメモリドライブ2020に接続されたバス2050と、ブートアッププログラム等のプログラムを記憶するためのROM2060とに接続され、アプリケーションプログラムの命令を一時的に記憶するとともに一時記憶空間を提供するためのRAM2070と、アプリケーションプログラム、システムプログラム、およびデータを記憶するための不揮発性記憶装置2080と、通信インタフェース2090とを含む。通信インタフェース2090は、駆動部21等と信号の授受を行うためのインタフェース部44および図示しないネットワークを介して他のコンピュータと通信するためのネットワークインタフェース50に相当する。なお、不揮発性記憶装置2080としては、ハードディスク(HDD)やソリッドステートドライブ(SSD:Solid State Drive)などを使用することが可能である。不揮発性記憶装置2080が、記憶部36に相当する。
CPU2040が、プログラムに基づいて実行する演算処理により、データ処理部32の各機能、たとえば、制御部42、データ収集部46、画像処理部48の各機能が実現される。
データ処理部32に、上述した実施の形態の機能を実行させるプログラムは、DVD-ROM2200、またはメモリ媒体2210に記憶されて、ディスクドライブ2030またはメモリドライブ2020に挿入され、さらに不揮発性記憶装置2080に転送されても良い。プログラムは実行の際にRAM2070にロードされる。
データ処理部32は、さらに、入力装置としてのキーボード2100およびマウス2110と、出力装置としてのディスプレイ2120とを備える。キーボード2100およびマウス2110が入力部40に相当し、ディスプレイ2120が表示部38に相当する。
上述したようなデータ処理部32として機能するためのプログラムは、コンピュータ本体2010に、情報処理装置等の機能を実行させるオペレーティングシステム(OS)は、必ずしも含まなくても良い。プログラムは、制御された態様で適切な機能(モジュール)を呼び出し、所望の結果が得られるようにする命令の部分のみを含んでいれば良い。データ処理部32がどのように動作するかは周知であり、詳細な説明は省略する。
また、上記プログラムを実行するコンピュータは、単数であってもよく、複数であってもよい。すなわち、集中処理を行ってもよく、あるいは分散処理を行ってもよい。
また、計算処理システム300内のハードウェアも、演算処理装置が並列化されていたり、GPGPU(General-purpose computing on graphics processing units)が使用される等の構成上の相違はあり得るものの、基本的な構成は、図5に示したものと同様である。
(脳機能結合に基づく、疾患/健常の判別器の生成処理、およびクラスタリング処理)
図6は、図2で説明したような相関行列から、診断マーカーとなる判別器を生成する過程およびクラスタリング処理について説明する概念図である。
機械学習の処理としては、判別器の生成は、いわゆる「教師あり学習」の処理が実行され、クラスタリング処理については、「教師なし学習」の処理が実行される。
そして、クラスタリング処理自体は、「教師なし学習」であって、医師の診断などの情報を用いているわけではないので、この結果得られる個々のクラスタは、データ駆動によって得られる患者のグループであり、患者がサブタイプに分かれる場合は、脳機能結合を特徴量とする「患者の層別化」の基礎となる。
図6に示すように、まず、複数のMRI装置において、健常者群、患者群において安静時のfMRI画像データが撮像され、計算処理システム300が、そのようなfMRI画像データに対して、後述するような「前処理」を実施する。続いて、測定された安静時機能結合的MRIのデータから、それぞれの被験者について、計算処理システム300が、脳領野のパーセレーション処理を行って、脳領域間(関心領域間)の活動度の相関行列を導出する。
続いて、相関行列の非対角成分について、後述するようにして、対応する測定バイアスを導出しておき、計算処理システム300は、その測定バイアスを相関行例の要素の値から減算することで、ハーモナイゼーション処理を実行する。
さらに、ハーモナイゼーション処理が行われた相関行列の要素値と、各被験者についての疾患ラベル(疾患あるいは健常を示すラベル)との間で、計算処理システム300は、過学習を抑制し、特徴選択を伴う識別器の生成を、後述するような「アンサンブル学習による識別器生成処理」として実施して、被験者の疾患あるいは健常を予測することが可能な疾患識別器(診断マーカ)として生成する。
一方で、計算処理システム300は、アンサンブル学習中に、疾患ラベルに対する識別器の生成処理の中で特定された特徴量(脳機能結合)のうちから、後述するようなクラスタリングのための特徴量選択処理を実施した上で、「教師なし学習」により、多重共クラスタリング処理を実行する。
以下、図6中の各処理について、さらに詳しく説明する。
[前処理から疾患識別器の生成およびクラスタリング処理までの概要]
(安静時機能的結合FCマトリックスの前処理および計算)
計測されたfMRIデータの、たとえば最初の10秒間はT1平衡を考慮に入れるために廃棄される。
前処理のステップでは、計算処理システム300は、スライスタイミングの校正、頭部に見られる体動アーチファクトの補正を行うためのリアライン処理、脳機能画像(EPI画像)と形態画像の共登録(co-registration)、歪補正、T1強調構造画像の分割、モントリオール神経学研究所(MNI)空間への正規化および、たとえば、6mmの半値幅の等方性のガウスカーネルによる空間の平滑化などの処理を実施する。
このような前処理のパイプライン処理については、たとえば、以下のサイトに開示がある。
http://fmriprep.readthedocs.io/en/latest/workflows.html
(脳領域の区画化(パーセレーション:Parcellation))
脳領域のパーセレーションについては、特に限定されないが、上述した「方法3」に従って、「表面ベースの方法」で実施することが可能である。
(生理的な雑音回帰)
生理的な雑音回帰は、以下の文献に開示されるCompCorを適用して実行される。
公知文献:Behzadi, Y., Restom, K., Liau, J., and Liu, T.T. (2007). A component based noise correction method (CompCor) for BOLD and perfusion based fMRI. Neuroimage 37(1), 90-101. doi: 10.1016/j.neuroimage.2007.04.042.
いくつかのスプーリアスの源(余計な信号源)を除去するために、6つの運動パラメータ、全脳、など回帰母数を備えた直線回帰が使用される。
(時間的なフィルタリング)
計算処理システム300は、時間的なバンドパスフィルターとして、たとえば、0.01Hzと0.08Hzの間の通過帯域を備えたバターワース・フィルタを使用して、時系列データに適用し、BOLD活動の特性である低周波の変動に分析を限定する。
(頭部運動)
フレームの位置ずれ(FD:Frame-wise displacement)が、個々の機能的なセッションにおいて計算され、頭部運動による機能的結合FCのスプーリアスの変化を低減するために、たとえば、FD>0.5mmであるボリュームが除去される。
FDは、スカラ量(つまり並進運動と回転における絶対的な変位の加算)として時間的に連続する2つのボリュームの間の頭部運動を表わす。
たとえば、後述する具体例においては、上述したような具体的なデータセットでは、スクラビングの後に除去されたボリュームの比率が、(平均±3×標準偏差)を超過した場合、その参加者のデータは分析から除外される。その結果、データセット全体では、35人の参加者が除去された。したがって、学習用データセットの中で683人の参加者(545人のHC、138人のMDD)を使用し、独立した検証データセットの中で444人の参加者(263人のHC、181人のMDDを患者)のデータが、以下の分析に対して使用される。
(機能的結合(FC)マトリックスの計算)
本実施の形態の具体例においては、機能的結合FCは、上述したようなパーセレーション法により領域の分割がされた後に、各参加者につき379個の関心領域(ROI)にわたってBOLD信号の時間的な相関として計算される。
機能的結合の計算では、特に限定されないが、ここでは、ピアソンの相関係数を使用する。
個々の可能な組のROIの前処理されたBOLD信号の時間的経過間のフィッシャーのz変換されたピアソンの相関係数が計算され、要素がそれぞれ2つのROIの間の結合の強度を表わす、379行×379列の対称的な結合マトリックスが構築される。
さらに、分析のために、結合マトリックスの下三角行列の71,631(=(379×378)/2)の機能的結合FCの値が使用される。
(脳活動バイオマーカーのためのハーモナイゼーション処理)
精神疾患と関連したビッグデータを収集するときには、上述した通り、1つのサイトが大規模な脳画像データ(人間の疾患と関係するコネクトーム)を集めることが、ほとんど不可能であるので、複数のサイトから画像データを収集することが必要である。
MRI装置(スキャナー)の型、プロトコルおよび患者層を完全に制御するのは難しい。したがって、収集したデータを解析するには、異種条件下で撮像された脳画像データが用いられる。
特に、疾患要因はサイト要因と交絡する傾向があるので、そのような異種条件下のデータに機械学習技術を適用することで疾患要因を抽出する場合、サイト間差は最も大きな障壁となる。
1つのサイト(あるいは病院)は、ほんの少数のタイプの精神疾患(例えば、サイトAからは主として統合失調症、サイトBからは主として自閉症、サイトCからは主として大うつ病など)をサンプリングしがちであるので、交絡が生じることになる。
適切にそのような異種条件下のデータを管理するために、サイト間のデータをハーモナイズ(調和)させることが必要である。
サイト間差は、本質的に2つのタイプのバイアスを含んでいる。
技術的なバイアス(つまり測定バイアス)および生物学的なバイアス(つまり標本バイアス)である。
測定バイアスは、撮像パラメータ、電界の強さ、MRI装置メーカーおよびスキャナの型のようなMRIスキャナの特性の違いを含んでおり、標本バイアスは、サイト間の被験者のグループの相違に関連している。
そこで、このようなサイト間差を補償するための「ハーモナイゼーション処理」が必要となる。このハーモナイゼーション処理の詳細については、上述した非特許文献8(Ayumu Yamashita他)に記載されており、その内容については後述する。
(アンサンブル学習による疾患識別器)
本明細書においては、「アンサンブル学習」という用語は、原学習データから復元抽出してK組の学習データを作成し、それぞれの学習データについて、独立に、機械学習処理によりK個の識別器を生成し、これらのK個の識別器を統合して、判別器を生成する処理のことをいうものとする。
特に、ここでは、特定疾患について、ある被験者の脳機能結合パターンにより、その被験者が疾患あるいは健常であるかを判別することを目的とするので、個々の識別器は、2クラス識別問題に対する識別器ということになる。
そして、原学習データから復元抽出してK組の学習データを作成する際には、後述するように、「アンダーサンプリング」および「サブサンプリング」が実行される。
ここで、L1正則化による識別器(LASSO(Least Absolute Shrinkage and Selection Operator)法)のような「特徴量選択を伴う学習処理による識別器」や、リッジ正則化法(L2正則化)などを使用する正則化学習法を用いることも可能である。
ここでは、「正則化学習法」とは、元の訓練データの特徴量全体を学習の対象とするものの、学習アルゴリズムにおいて、モデルを学習する際に、複雑さが増すことに対するペナルティを設け、このペナルティを訓練誤差に加えた量が最も小さくなる学習モデルを求めるようにすることにより汎化性能を高めようとする学習方法のことを意味する。そして、L1正則化は、ペナルティとして学習モデルの(特徴量に対応する)パラメータの絶対値の総和を用いるものであり、L2正則化は、ペナルティとして学習モデルのパラメータの二乗の総和を用いるものである。なお、ペナルティとしては、モデルに使用される特徴量の個数そのものを用いるL0正則化もあり得る。
また、LASSO法(L1正則化)は、いわゆるスパース推定が可能な手法であり、その派生形として、Elastic Net法、Group Lasso法、Fused Lasso法、Adaptive Lasso法、Graphical Lasso法などがある。
一方で、「特徴量選択を伴う学習処理による識別器」としては、「ランダムフォレスト法」のように、識別器の生成において、特徴量の選択とともに、特徴量の重要度が、併せて、得られる手法を利用することも可能である。
なお、このような「アンサンブル学習による識別器の生成手法」としては、以後は、LASSO法を中心として説明するものの、その手法は、上述のようなものに、限定されるものではない。たとえば、パーセレーション法として、辞書学習(DictionaryLearning)の手法を用いてデータ依存的に解析対象脳領域を設定し、機能的脳結合の値としてタンジェント法(tangent-space covariance)を用い、データセット内の脳機能結合FCを後述するComBat法を用いて施設間補正を行い、リッジ正則化で識別器を生成するアンサンブル学習法などであってもよい。パーセレーション法や、脳機能結合の計算方法、ハーモナイゼーション法、識別器の生成方法などは、別の組み合わせであってもよい。
後述するように、本実施の形態では、このようなアンサンブル学習において、識別器を学習させる際に、各特徴量の識別という機能を達成するための「重要度」が特定される。
(クラスタリングのための特徴量選択とクラスタリング)
K組の学習データに対して、「アンサンブル学習」として、K個の識別器を生成する過程で、それぞれの識別器生成において使用された「第1の特徴量」の和集合から、さらに、「教師なし学習」によるクラスタリングを実行する際の第2の特徴量の集合を特定するための処理が実行される。
特に、限定されないが、たとえば、以下のようにして「重要度」が決定される。
i)アンサンブル学習で、K個の識別器を生成する学習手法が、「特徴量選択を伴う学習処理」の場合は、識別器の生成において選択された「第1の特徴量」の和集合において、K個の識別器を生成する際に使用される頻度に応じて、特徴量をランキングする。
ii)アンサンブル学習で、K個の識別器を生成する学習手法が、「ランダムフォレスト法」のように、識別器の生成において、特徴量の重要度が得られる手法の場合は、特徴量のランキングリストは、このような重要度に従って生成される構成とすることができる。
iii)アンサンブル学習で、K個の識別器を生成する学習手法が、(必ずしも特徴量選択を伴わない)「リッジ正則化法(L2正則化)」であって、特徴量の重み付き和を引数とする識別器の生成手法である場合は、個々の識別器での各特徴量の重み係数の絶対値を、K個の識別器にわたって集計した中央値を重要度として特徴量のランキングリストが生成される。なお、重要度としては、このような「中央値」に必ずしも限定されるものではなく、たとえば、「K個の識別器にわたって積算した積算値」など他の統計的な代表値を用いることとしてもよい。
i)~iii)のようにして生成されたランキングのリストの中で、上位の所定の個数の特徴量を、「第2の特徴量」とする、という構成とすることが可能である。
なお、「第2の特徴量」を特定する条件としては、ランキングリストの上位の所定の個数ということに限られず、たとえば、ランキングリストにおいて、所定の頻度以上であること(K個の識別器の生成で、ある割合以上選択されるという頻度そのものを用いる)を条件としてもよい。
以上のようにして、選択された特徴量に基づいて、後述するように、教師なし学習である「多重共クラスタリング法」により、クラスタリング処理(患者層別化)が実行される。
[2クラス分類のための識別器生成処理]
以下では、図6で説明した処理のうち、アンサンブル学習による識別器生成について、より詳しく説明する。
すなわち、2クラス分類のための分類器生成処理、より特定的には、疾患識別器(「健常」または「疾患」についての2クラス分類器)の訓練データとして、学習用データセットを使用して、MDDのためのバイオマーカーを構築する処理を例として説明する。
ここでは、精神疾患のうち、大うつ病性障害を例として、すなわち、従来の症候に基づく診断手法により医師によって大うつ病性障害と診断された患者群を例として、分類器を生成する処理を説明する。そして、図8に示した疾患識別器生成部3008が、患者群と健常者群とを判別するための診断の補助情報を出力するような分類器を生成するために実行する処理の例を説明する。
そこで、以下では、機能的結合FCに基づいて、健常者(HC)およびMDD患者を識別するMDD識別器を構築する手続きについて説明する。
以下では、疾患識別器(MDD識別器)を作成するための「特徴量選択を伴う学習処理による識別器」としては、例示として、L1正則化による識別器の学習方法(LASSO法)を用いるものとして説明する。
そして、後述するように、MDD診断と関係する機能的結合FCを特定するために、各機能的結合FCの疾患識別器の構築に対する「重要度」に応じて、クラスタリングに使用する特徴量が選択される。
図7および図8は、データセンター200の記憶装置210に格納されたデータに基づき、ハーモナイゼーション処理、疾患識別器生成処理、クラスタリング分類器生成処理、判別処理を実行する計算処理システム300の構成を説明するための機能ブロック図である。
なお、ここでは、「判別処理」には、疾患の判別(疾患あるいは健常の判別)と、検査対象の被験者がいずれの「クラスタ」(サブタイプ)に属するのかを判断する分類処理を含むものとする。
図7を参照して、計算処理システム300は、記憶装置210からのデータおよび計算途中で生成されるデータを格納するための記憶装置2080と、記憶装置2080内のデータに対して、演算処理を実行する演算装置2040とを含む。演算装置2040としては、たとえば、CPUが相当する。
演算装置2040は、プログラムの実行により、患者群および健常者群のMRI計測データ3102について相関行列の要素を算出し、相関行列データ3106として、記憶装置2080に格納するための相関行列算出部3002と、ハーモナイゼーション処理を実行するハーモナイゼーション算出部3020と、ハーモナイゼーション処理の結果に基づいて、疾患識別器の生成処理、クラスタリング分類器の生成処理および生成された疾患識別器またはクラスタリング分類器による判別処理を実行する学習及び判別処理部3000とを含む。
図8は、図7の構成をより詳しく説明する機能ブロック図である。
また、図9は、アンサンブル学習による疾患識別器を生成するための機械学習の手続きを説明するためのフローチャートである。
そこで、まずは、図6に示したように、ハーモナイゼーション処理、アンサンブル学習による識別器生成(疾患識別器)の生成までの処理について、図8および図9を参照して説明する。
まず、データセンター200の記憶装置210には、各計測サイトから被験者(健常者および患者)のfMRI計測データ、被験者の属性データ、測定パラメータが、「学習用データセット」として収集されていることを前提とする。
図8および図9を参照して、疾患識別器(「健常」または「疾患」についての2クラス分類器)の訓練データとして、このような学習用データセットを使用して、MDD判別のためのバイオマーカーを構築する。それは、つまり、71,631の機能的結合FCの値に基づいて、健常者群(健常(HC)との診断ラベルの個体群)およびMDD患者群(大うつ病性障害との診断ラベルの個体群)を識別するものである。
以下に説明するように、MDDに対する識別器(以下、「MDD識別器」と呼ぶ)の生成のための学習処理において、L1正則化(LASSO法)によるロジスティク回帰分析(スパースモデリング法の1種)を使用して、71,631個の機能的結合FCの中からの機能的結合FCの最適な部分集合を選択する。
一般に、L1正則化を使用すると、いくつかのパラメータ(以下の説明では、ウェイトの要素)を0にすることができる。つまり、特徴選択を行っていることになり、スパースモデルになる。
ただし、スパースモデリングの手法としては、LASSO法には限定されず、後述するように、変分ベイズ法をロジスティック回帰に応用したスパースロジスティック回帰(SLR:Sparse Logistic Regression)を用いるなど、他の手法を用いることも可能である。
図9を参照して、MDD識別器に対する学習処理が開始されると(S100)、予め準備され(記憶装置2080に格納され)た学習データセットを使用して(S102)、相関行列算出部3002が、結合マトリックスの成分を算出する。
続いて、ハーモナイゼーション算出部3020が、算出した測定バイアスを使用して、ハーモナイゼーション処理を実施する(S104)。
後述するとおり、ハーモナイゼーション処理は、トラベリングサブジェクトを用いた方法が望ましいが、他の方法であってもよい。
たとえば、発見データセットおよび独立検証データセット間で、後述するようなcombat法を用いて、データセット間のハーモナイゼーションを行う構成とすることも可能である。
続いて、疾患識別器生成部3008は、学習用データに対して、いわゆる「アンサンブル学習の手法」であって、「入れ子構造の交差検証(Nested Cross Validation)」の手法を修正した方法により、MDD識別器を生成する。
まず、疾患識別器生成部3008は、学習用データに対して、「K分割交差検証」(K:自然数)(外側の交差検証)を使用して学習処理を実行するために、たとえば、K=10として、学習用データを10分割する。
すなわち、疾患識別器生成部3008は、K分割(10分割)されたうちの1つの部分データセットを検証用の「テストデータセット」とし、残りの(K-1)分割(9分割)分のデータを訓練データセット(トレーニングデータセット)に設定する(S108,S110)。
続いて、疾患識別器生成部3008は、トレーニングデータセットに対して、「アンダーサンプリング処理」および「サブサンプリング処理」を実行する(S112)。
ここで、「アンダーサンプリング処理」とは、トレーニングデータセットにおいて、分類の対象となる特定の属性データ(2種以上)にそれぞれ対応するデータの数がそろっていない場合に、この数をそろえるために数が多い方の属性のデータを除いて同数となるように行う処理を意味する。
ここでは、トレーニングデータセットにおいて、MDD患者群の被験者数と、健常者群の被験者数が等しくないために、これをそろえるための処理を行うことに相当する。
さらに、「サブサンプリング処理」とは、トレーニングデータセットから所定数のサンプルをランダムに抽出する処理を意味する。
すなわち、ステップS108~S118、S122を経由して、K回繰り返される交差検証において、各交差検証では、トレーニングデータセットがMDD患者および健常者HCの数に関して不均衡であるので、分類器を構築するためのアンダーサンプリング方法を行うこととし、かつ、サブサンプリング処理として、所定数、たとえば、130人のMDD患者および同数の130人の健常者が、トレーニングデータセットからランダムにサンプリングされる。
なお、130人という数値は、このような値に限定されるものではなく、学習データセット中のデータ数(後述する「データセット1」では683人)、分割数K(ここでは、たとえば、K=10)、並びに、分類の対象となる特定の属性に含まれるデータ数の不均衡の程度に応じて、適宜、上述したようなアンダーサンプリングが可能となるように決定される。
このようなサブサンプリング処理を実行するのは、アンダーサンプリングは、除外されたデータを使用して識別器が学習することができなくなるという点で不利であるため、この不利を除去するために、無作為抽出手続き(つまりサブサンプリング)を、M回(M:自然数、たとえば、M=10)繰り返す処理とするためである。
なお、後述するように「層別化」のための「分類器」を生成する際の「特徴量選択」のためにも、このようなアンダーサンプリング処理およびサブサンプリング処理を実行することには、技術的な意義が存在するので、その点については、後述する。
続いて、疾患識別器生成部3008は、サブサンプリングされたサブサンプル1~10の各々に対して、ハイパーパラメータの調整処理を実行する(S114.1~S114.10)。
ここで、各サブサンプルにおいて、以下のようなロジスティック関数を使用することで、識別器サブモデルが生成される。このようなロジスティック関数は、サブサンプル内において、MDDクラスに属する参加者の可能性を以下のように定義するために使用される。
ここで、ysubは、参加者のクラス・ラベル(MDD, y=1; HC, y=0)を表わし、csubは、与えられた参加者に対するFCベクトルを表わし、wはウェイトベクトルを表わす。
ウェイトベクトルwは、以下の評価関数(コスト関数)を最小化するように決定される(LASSO計算)。
LASSO計算では、コスト関数において、ウェイトベクトルの各要素の絶対値(1次)の総和(L1ノルム)が第2項として存在する。
ここで、λはハイパーパラメータを表わし、評価に適用される収縮量をコントロールする。
疾患識別器生成部3008は、各サブサンプルにおいて、これも特に限定されないが、所定数のデータをハイパ―パラメータ調整用データとし、残りのデータ(たとえば、n=250または248人のデータ)を用いてウェイトベクトルwを決定する。このとき、疾患識別器生成部3008は、特に限定されないが、たとえば、ハイパーパラメータλが0<λ≦1.0であるものとして、この区間を、P等分(P:自然数)、たとえば25等分した各値のλを使用して、上記のようなLASSO計算により、ウェイトベクトルwを決定する。
このとき、上述の通り、「入れ子構造の交差検証」として、ハイパーパラメータの調整については、「内側の交差検証」として実行する。内側の交差検証では、外側の交差検証の「テストデータセット」は、使用されない。
その上で、疾患識別器生成部3008は、生成された各λの値に対応するロジスティック関数によりハイパ―パラメータ調整用データについて判別性能(たとえば、正確度)を比較し、最も判別性能の高いλに対応するロジスティック関数を決定する(ハイパーパラメータの調整処理)。
続いて、疾患識別器生成部3008は、現在の交差検証のループの中で生成された各サブサンプルに対応するロジスティック関数の出力値の平均を出力するものとして「識別器サブモデル」を設定する(S116)。各サブサンプルにおいて算出された識別器の出力値の平均により識別性能を判定しているという点でも、これは、「アンサンブル学習」の一種ということができる。
疾患識別器生成部3008は、ステップS110で準備されたテストデータセットを入力として、現在の交差検証のループの中で生成された識別器サブモデルの検証を実行する(S118)。
なお、アンダーサンプリングとサブサンプリングにより、サブサンプルを生成して、各サブサンプルにおいて、特徴選択を実行して識別器サブモデルを生成する方法としては、上記のようなLASSO法とハイパーパラメータの調整とを実行する方法以外に、他のスパースモデリングの手法を利用してもよい。
疾患識別器生成部3008は、交差検証のループをK回(ここでは10回)分終了していないと判断すると(S122でN)、K分割されたデータにおいて、それまでのループで使用したのとは異なる別の部分データセットをテストデータセットに設定し、残りの部分データセットをトレーニングデータセットに設定して(S108,S110)、処理を繰り返す。
一方で、疾患識別器生成部3008は、交差検証のループをK回(10回)分終了している場合(S122でY)、入力データに対して、K×M個(この場合は、10×10=100個)のロジスティック関数(識別器)の出力の平均を出力するものとして、MDDに対する識別器モデル(MDD識別器)を生成する(S120)。
結局のところ、MDD識別器は、K×M個の識別器の出力の平均を、その識別出力とするという意味で、「アンサンブル学習」の結果得られる「識別器」であるといえる。
MDD識別器の出力(診断の確率値)が、0.5を超えるときは、MDD患者を示す指標と見なすことができる。
さらに、本実施の形態では、このようにして生成されたMDD識別器の性能の評価指標として、マシューズ相関係数(MCC:Matthews correlation coefficients)、ROC曲線(Receiver Operatorating Characteristic curve、受信者動作特性曲線)についてのROC曲線下面積(AUC: area under the curve)、正確度(Accuracy)、感度(Sensitivity)、特異度(Specificity)を使用する。
なお、それぞれのサブサンプルにおいて特徴選択された特徴量(この場合、測定バイアスに対するハーモナイゼーション処理後の相関行列の要素)を使用して、対象疾患(たとえば、MDD)の識別器を生成する手法は、このような複数個の識別器サブモデルの出力の平均処理によるものに限られず、多数決による処理としてもよいし、あるいは、特徴選択された特徴量に対して、他のモデリング手法、特に、他のスパースモデリング手法を使用して識別器を生成する構成としてもよい。
(MDD識別器に用いるデータの例と性能)
既述したように、機械学習アルゴリズムを使用する信頼できる分類器および回帰モデルの構築については、多数の撮像サイトから集められた大規模なサンプルサイズのデータを用いることが必要である。
そこで、以下では、4つの異なる撮像サイトから集められたMDD患者を含む、約700人の参加者の学習用の安静時fMRIデータセットを使用して検討する。
図10は、このような学習用のデータセット(データセット1)のデモグラフィック特性を示す図である。
データセット1は、上述したSRPBS中のデータである。
図11は、独立検証のデータセット(データセット2)のデモグラフィック特性を示す図である。
データセット2も、基本的には、上述したSRPBS中のデータである。
すなわち、以下の分析では、2つの安静時機能的MRI(rs-fMRI)データセットを使用する。
(1)図10に示すように、データセット1は、713人の参加者に対するデータを含んでいる(4つのサイトから564人の健常者群HC、3つのサイトからの149人のMDD患者群)。
(2)図11に示すように、データセット2は、449人の参加者に対するデータを含んでいる(4つのサイトからの264人の健常者群HC 、4つのサイトからの185人のMDD患者群)。
また、併せて「抑うつの症候」は、各データセットの大部分の参加者から得られたベックうつ病自己評価尺度(BDI:Beck Depression Inventory)IIを使用して評価されている。
データセット1は、「学習用データセット」であり、MDDの識別器およびクラスタリングの分類器を構築するのに使用される。
参加者は、それぞれ10分間の単一の安静時機能的MRI(rs-fMRI)セッションにて計測が実行されている。
ここでも、統一された撮像プロトコルの下で安静時機能的MRI(rs-fMRI)データが取得されている(http://www.cns.atr.jp/rs-fmri-protocol-2/)。
ただし、画像撮影が、すべてのサイトで同じパラメータを使用して行なわれたことを保証することは実際には困難であり、2つの位相変調方式の方向(P→AとA→P)、2社のMRI装置メーカー(シーメンスとGE)、コイルの3つの異なる数(12,24,32)および3つの型番のスキャナが計測に使用されている。
安静時機能的MRI(rs-fMRI)のスキャン中に、参加者は、原則として、以下のような教示がされている。
「リラックスしてください。眠らないでください。中央の十字線マーク上を注視して、特定のことに関して考えないでください。」
データセットにおける「デモグラフィック特性」は、いわゆる「人口統計学」において使用される特性であって、年齢、性別、などの他、診断名のような表中の属性を含む。
なお、図10および図11において、括弧中の人数は、BDIスコアのデータを有する参加者の数を示す。
デモグラフィック分布は、すべて学習用データセットの中でMDDとHCの個体群間で統計的に有意な差はない(p>0.05)。
データセット2は、 「独立した検証データセット」であり、MDDの分類器およびクラスタリングの分類器をテストするために使用される。
データセット2に対する撮像を行ったサイトは、データセット1に含まれていない。
年齢のデモグラフィックな分布は独立した検証データセットの中でMDDとHCの個体群間で一致しているものの(p>0.05)、性比のデモグラフィックな分布は独立した検証データセットの中でMDDとHCの個体群間で一致していない(p<0.05)。
(サイト効果の制御)
また、以下では、機能的結合FCの上のサイト効果を制御するために、後述するような、学習用データセットに対するトラベリングサブジェクトのハーモナイゼーション方法を使用するものとして説明する。
ただし、ハーモナイゼーション方法としては、この方法に限定されるものではなく、たとえば、これもComBat法など他の方法を使用してもよい。
なお、ComBat法については、たとえば、以下の文献に開示がある。
公知文献:Johnson WE, Li C, Rabinovic A. “Adjusting batch effects in microarray expression data using empirical Bayes methods.” Biostatistics 8, 118-127 (2007).
トラベリングサブジェクトのハーモナイゼーション方法を使用することで、純粋なサイト間差(測定バイアス)を除去することが可能となる。
なお、独立した検証データセットに含まれたサイトに対しては、トラベリングサブジェクトのデータセットが存在しなかったので、独立した検証データセットの中でサイト効果のコントロールのためには、ComBat法によるハーモナイゼーション方法を使用している。
図12は、全撮像サイトについて、学習用データセットに対するMDDの予測性能(出力の確率分布)を示す図である。
学習用データセットに対して、0.5のしきい値によって、識別器モデルからの出力において、MDD患者と健常者の個体群に対応する2つの診断の確率分布が、明白に、右(MDD)と左(HC)へ分離されている。
識別器モデルは、MDD患者を66%の正確度でHC個体群から分離する。
対応するAUCは、0.77であり、高い識別力を示した。
また、MCCは、約0.33である。
図13は、各撮像サイトについて、学習用データセットに対するMDDの予測性能(識別器の出力の確率分布)を示す図である。
図13から、全データセットだけでなく3つの撮像サイト(サイト1、サイト2、サイト4)の個々のデータセットに対して、ほとんど同じ程度に高い分類精度が達成されることがわかる。
なお、サイト3(SWA)のデータセットには、単に健常者群があるだけであるが、その確率分布は、他のサイトの健常者群のものに相当している。
(識別器の汎化性能)
図14は、独立した検証データセットにおけるMDDの識別器の出力の確率分布を示す図である。
すなわち、独立した検証データセットを使用して、識別器モデルの汎化性能がテストされる。
MDDに対して、図12の処理において100個(10分割×10サブサンプリング)のロジスティック関数の識別器が機械学習により生成されており、独立した検証データセットを、生成された100個の識別器のすべて(識別器の集合としての識別器モデル)に入力する。
そして、各参加者につき100個の識別器の出力の平均(診断の確率)をとり、平均された診断の確率値が>0.5だった場合、その参加者の診断ラベルとして、大うつ病性障害が該当するものとする。
独立した検証データセットにおいて、約70%の精度で、生成された識別器モデルは、MDD個体群をHC個体群から分離している。
対応するAUCは、0.75となり、高い識別能力となる能力を示した(並べかえ検定p<0.01)。
独立検証データセットに対して、0.5のしきい値によって、識別器モデルからの出力において、MDD患者と健常者の個体群に対応する2つの診断の確率分布が、明白に、右(MDD)と左(HC)へ分離されている。
感度は68%であり、特異度は71%である。これは、0.38という高いMCC値となっている(並べかえ検定 p<0.01)。
図15は、各撮像サイトについて、独立検証データセットに対するMDDの識別器の出力の確率分布を示す図である。
4つの撮像サイトの全データセットだけでなく個々のデータセットに対して高い分類精度が達成されることがわかる。
[被験者データに対するクラスタリング処理]
以下では、図6で説明した処理のうち、クラスタリングのための特徴量選択と、選択された特徴量によるクラスタリング処理について、より詳しく説明する。
すなわち、図6において、「特徴量選択」と「クラスタリング処理」として説明した処理を、図8における疾患識別器生成部3008とクラスタリング分類器生成部3010とが実行する処理として説明をする。
図16は、特徴量を選択して、教師なし学習によりクラスタリングを行う処理を説明するためのフローチャートである。
以下では、図9において説明したような「2クラス識別器」の学習処理において、各識別器サブモデルの生成において使用された特徴量(脳機能結合)を、ランク付けして、上位の所定数の特徴量を使用して、教師なし学習によりクラスタリングを行う処理について説明する。
上述したように、脳機能結合は、脳に対するパーセレーションの手法に応じて、70000次元を超えるような高次元であって、通常の方法では、教師なし学習でのクラスタリング処理を実行することは一般には困難である。本実施の形態の方法では、このようなクラスタリング問題に対して、「教師あり学習での識別器生成」において、特徴量の重要度に応じたランキングを実行し、このランキングに基づいて選択された特徴量による「教師なし学習によるクラスタリング」を組み合わせることで、このようなクラスタリング処理を可能としている。
なお、以下では、便宜上、疾患識別器の生成処理とは別のものとして、クラスタリング処理を説明するが、図16において、ステップS200~S210は、図9におけるステップS100~S120と同等の処理であり、疾患識別器の生成処理とクラスタリング処理とは、一連の処理として実行することが可能である。
図16を参照して、クラスタリングの学習処理が開始されると、疾患識別器生成部3008は、健常者群Nh人、うつ病群Nm人の被験者データを準備し(S202)、疾患識別器生成部3008は、被験者の機能的脳活動データにおいて、脳領野の分割(パーセレーション)処理、脳機能結合値の算出、および、ハーモナイゼーション処理とを実行する(S204)。
続いて、疾患識別器生成部3008は、Ncv分割交差検証(Ncv:自然数、かつ、Ncv≧2)のためのデータ分割を実行して、各分割されたデータについて、トレーニングデータセットと、テストデータセットとを準備し、各トレーニングデータセットについては、Ns個のテストデータのサブセットを生成するように、アンダーサンプリングおよびサブサンプリングを実行する(S206)。
さらに、疾患識別器生成部3008は、サブサンプルされた各サブサンプルに対して、
特徴量選択を伴う学習方法により識別器を生成する(S208)。
なお、ここでは、図9と同様にして、L1正則化(LASSO)で特徴量選択が実行されるものとする。
このようなステップS206~S208の処理を、Ncv個に分割された学習データセットについて、トレーニングデータセット(分割されたデータセットの(Ncv-1)個分)とテストデータセット(分割されたデータセットの1個分)との組み換えを順次実行して、交差検証がNcv回実施されるまで、繰り返す。
このようにして生成された(Ns×Ncv)個の識別器の平均を出力とする統合された識別器を、疾患識別器(診断マーカ)として生成する(S210)。
ここまでの処理は、上述したとおり、図9におけるステップS100~S120と同等の処理である。
一方で、クラスタリング分類器生成部3010は、Ncv回繰り返されるステップS206~S208において、特徴量選択を伴う学習方法により識別器を生成する際に選択された特徴量(脳機能結合)の和集合において、特に限定されないが、この和集合中の特徴量について、選択された回数でランキングを実行する(S220)。
ここでは、「特徴量として選択された回数」を、このランキングにおける、その特徴量の重要度と呼ぶことにする。
言い換えると、LASSO法で、たとえば、図9に示した例によれば、100(=10×10)個の識別器が生成されるが、各識別器において重みがゼロでない脳機能結合について、+1となるように、選択回数をカウントする。カウント回数の大きい順で、重要な結合としてランキングする。
続いて、クラスタリング分類器生成部3010は、うつ病患者群に対してクラスタリングを教師なし学習で実施するために、上記和集合中から重要度に応じて、たとえば所定数の特徴量を選択する(S222)。
さらに、クラスタリング分類器生成部3010は、教師なし学習の方法として、後述するような多重共クラスタリング法を用いて、クラスタリング処理を実施する(S224)。
以上の処理により、クラスタリング分類器生成部3010は、うつ病患者群に対するクラスタリング分類器を生成する(S226)。
II.モデルの使用フェーズ
すなわち、以上の処理により、クラスタリング分類器生成部3010により、各クラスタについて、観測データから、そのような観測データを生成するような確率分布のモデルがそれぞれ特定され、記憶装置2080に各モデルの情報が疾患識別器データ内に格納される。そして、判別値算出部3012は、クラスタリング分類器として、学習データ以外の入力データに対しては、このよう各確率分布のモデルに基づいて、それぞれのクラスタに当該入力データが属する事後確率を算出し、最大の事後確率となるクラスタに、当該入力データが属するとの分類結果を出力することになる(MAP推定法(Maximum A posteriori Probability estimation method))。
[本手法による学習方法および学習済みモデルについての追加的説明]
なお、以上の説明では、クラスタリング処理は、疾患識別器生成部3008により、健常者群Nh人、うつ病群Nm人の被験者データに対して実施された識別器の生成処理に基づくものとして説明した。しかし、本実施の形態のクラスタリング方法は、このような場合に限定されることなく、「うつ病患者群」以外の疾患群、たとえば、「統合失調症の患者群」「自閉症の患者群」「強迫神経症の患者群」など、他の精神疾患の患者群のクラスタリングに使用することも可能である。
あるいは、より一般的に、人間が経験的に分類してきた属性ラベル(たとえば、その人の性格、その人の得意分野など)と、脳活動の時間変化の領野間の相関のパターンとが一定の関係性を有することが判明した属性については、当該属性ラベルにより分類される被験者について、データ駆動により「当該属性に属する被験者群のクラスタリング」(サブタイプへの分類)を実施することに使用することも可能である。
(アンダーサンプリングおよびサブサンプリング処理)
以上説明した処理では、「アンダーサンプリングおよびサブサンプリング処理」を実行しているので、その技術的意義について、簡単に説明する。
まず、「アンダーサンプリング処理」の効果については、識別器における識別の境界が適切に設定されることが挙げられる。
たとえば、2クラス分類のタスクを考えた場合、学習データ中で、各クラスに属するデータの数に偏りが少ないほど、処理フロー中で、識別器の性能(たとえば、正確度)の評価の精度が高くなる。
図9におけるステップS114.1~114.10において、ハイパーパラメータの設定には、「最も判別性能の高いλに対応するロジスティック関数を決定する」との処理を実行するので、「判別性能」の評価が正確に行われることが必要になる。
極端な例として、クラス1に属するデータ数が100個で、クラス2に属するデータ数が1個というような学習データについて識別器の学習をする場合、識別器として、すべてのデータがクラス1であると判別したとしても、正確度等には、大きな影響が出ない、という事態が発生してしまう。この点で、2つのクラスにそれぞれ属するデータ数を、ランダムなサンプリングによりそろえることには、意義がある。
また、サブサンプリングについては、そもそも、これは、以下の理由により、アンダーサンプリングと併せて、複数回実施することが前提である。
第1には、アンダーサンプリングおよびサブサンプリング処理は、仮に、ランダムサンプリングにより実施されるとしても、1回の処理だけでは、データに偏りが生じる可能性がある。
第2に、以下に説明するように、図9のステップS108~S122で繰り返される「識別器」の生成は、上述したように、「アンサンブル学習」により実施される。
このとき、各識別器の生成において、識別に対する特徴量の重要度が決定される。
重要度は、「特徴量選択を伴う学習処理」においては、当該特徴量が選択されること、あるいは、「特徴量選択を伴わない学習処理」においては、算出された識別に対する当該特徴量の重みというように、識別に対する各特徴量の貢献度に応じて、決定される。
以下では、このような重要度の決定における「アンダーサンプリングおよびサブサンプリング処理」の意義を、「特徴量選択を伴う学習処理」を例として説明する。なお、L2正則化のような「特徴量選択を伴わない学習処理」であっても、「当該特徴量の重みが大きくなる」という事象は、「特徴量として選択される」という事象と、基本的には、同様の技術的理由から生じるものと考えることができる。
ここで、「特徴量選択を伴う学習処理」としては、たとえば、上述したLASSO法のような、いわゆる「スパースモデリング」の手法が挙げられる。
スパースモデリングでは、特徴量がスパースに選択される、すなわち、特定の特徴量については、重みが非ゼロとなるのに対して、他の特徴量については、重みが0となることで、特徴量が選択される。このような特徴量のスパースな選択が実現される理由の1としては、「判別(識別)処理」に対して「類似の寄与をする特徴量のグループ」が存在する場合に、当該グループのうちの1つの特徴が選択され、当該グループの他の特徴量の重みが0となるように、学習処理が進行させるための「特徴量の数に対応するペナルティ項」が存在していることが、挙げられる。LASSO法については、特に、この傾向が顕著である。
つまり、「判別処理」に対して、特徴量Aと特徴量Bが、同様に関与する、たとえば、特徴量Aと特徴量Bの相関が高い、という場合には、特徴量として、特徴量Aのみを選択しても、判別性能を落とすことなく、判別処理を行うことができることになる。
しかしながら、たとえば、クラスタリングにおいては、特徴量Aと特徴量Bとを双方考慮する必要がある場合が想定される。ところが、このような「スパース化」を実行して、1回の識別器の生成過程における判別処理に対する寄与度だけで、特徴量を選択してしまうと、クラスタリングに対する「特徴量選択」としては、十分でない可能性が生じてしまう。
図17は、このような「特徴量選択を伴う学習処理」により、複数個(たとえば、Nch個)の特徴量が存在する場合に、特徴量の選択が実施される概念を示す図である。
図17を参照して、学習処理の対象となる被験者群には、健常者群と患者群とが含まれるものとする。
健常者群の被験者には、ラベルHが対応付けられており、健常者群には、サブタイプh1とサブタイプh2とが含まれるものとする。
患者群の被験者には、ラベルMが対応付けられており、患者群には、サブタイプm1、サブタイプm2、サブタイプm3とが含まれるものとする。
ここで、観測量では、健常者群、患者群が、いくつのサブタイプに分かれるのかは未知であり、サブタイプの識別ラベルは明示的には被験者に対応付けられておらず潜在的なラベルであるものとする。
そして「クラスタリング」の目的は、観測量からこれらのサブタイプへのデータ駆動によるクラスタリングを実行することである。
上記のような健常者群と患者群とから、被験者を「アンダーサンプリング」および「サブサンプリング」をランダムに実施するので、図17の「被験者群」に示すように、健常者群および患者群のそれぞれから、たとえば、点線で囲んだ部分の被験者が選択されるものとする。
さて、ラベルMとラベルHを識別するために使用しうる特徴量(脳機能結合の相関値)は、各被験者を特徴づける特徴量(全体でNch個)の脳機能結合のうち、図17において一点鎖線で示した範囲の特徴量(図16のステップS220における「脳機能結合の和集合」)であるものとする。
そして、特徴量選択を伴う学習処理(ここでは、LASSO法による処理)により、ラベルMとラベルHを識別するための識別器を学習した結果として、図17における一点鎖線内のさらに、黒丸で示した特徴量が選択されるものとする。
図18は、アンダーサンプリングおよびサブサンプリングの処理後に、特徴量選択を伴う学習処理により1つの識別器を生成する際に最終的に選択される特徴量を示す概念図である。
図18に示すように、一点鎖線内のラベルMとラベルHを識別するために使用しうる特徴量は、さらに点線の枠で示すように、互いに相関が強い特徴量のグループに分かれているものとする。
LASSO法では、このような点線の枠内のグループごとに、1つの特徴量が選択されることで、スパース化が実現していることになる。
図19は、アンダーサンプリングおよびサブサンプリングの処理を複数回実施して、識別器を生成される際に、特徴量が選択される様子を示す概念図である。
図19に示すように、サブサンプリングを、たとえば、Ns回実施しするものとすると、各回において、健常者群と患者群から、それぞれ異なる被験者がサブサンプルされる。
そして、各サブサンプルに対して、特徴量選択を伴う学習処理により、ラベルMとラベルHを識別するための識別器を学習した結果として、各サブサンプル対応して、上述した和集合である一点鎖線内のうちで、相関の強い各グループからは、黒丸で示したように、それぞれの識別器で異なる特徴量が選択されることになる。
結果として、アンダーサンプリングおよびサブサンプリングの処理を複数回実施することで、ラベルMとラベルHを識別するために使用しうる特徴量の和集合が選択されることになる。
本実施の形態では、各サブサンプルに対して、LASSO法により、特徴量選択を伴う識別器の学習処理により選択される特徴量を、選択される頻度に応じて、ランク付けする。
そして、ランク付けのうちの上位から所定の個数、たとえば、100個の特徴量を用いて、図16のステップS224において、後述するような「多重共クラスタリング」により、教師なし学習によるクラスタリングを実施する。
なお、以上の説明では、「特徴量選択を伴う識別器の学習処理」として、LASSO法を例として説明し、選択された頻度のランキングでクラスタリングのための特徴量の選択を実施する構成としている。
ただし、上述のとおり、本実施の形態のクラスタリング処理において、「特徴量選択を伴う識別器の学習処理」は、このような方法に限定されず、たとえば、ランダムフォレスト法のような方法でもよく、所定の重要度に従って、クラスタリングのための特徴量の選択を実施することも可能である。
たとえば、上述したように、ランダムフォレスト法では、識別器の学習処理において、ジニ不純度(Gini impurity)をもとに特徴量の重要度(importance)が併せて算出されるので、この重要度に基づいて、特徴量をランク付けし、ランク付けのうちの上位から所定の個数の特徴量を用いて、図16のステップS224における「多重共クラスタリング」により、教師なし学習によるクラスタリングを実施するとの構成とすることも可能である。
また、「アンサンブル学習として識別器を学習する処理」としては、リッジ正則化法などを使用し、上述したように、重み係数の絶対値を集計した中央値に対応する重要度に従って、特徴量をランキングして、ランク付けのうちの上位から所定の個数の特徴量を用いて、クラスタリングのための特徴量の選択を実施することも可能である。
[多重共クラスタリング処理]
以下では、図16のステップS224における「多重共クラスタリング」について、その概念を説明し、「多重共クラスタリング」との用語の定義を行う。
(通常のクラスタリング手法)
前提として、「クラスタリング」とは、コンピュータにより実行される教師なし学習によるデータ分類手法をいい、より特定的には、与えられたデータを外的基準なしに自動的に分類する手法をいうものとする。これに対して、「クラス分類」は、一般に、「教師あり学習」による分類手法のことをいう。また、「クラスタ」とは、内的結合と外的分離の性質を持つデータの部分集合である、と定義される。ここで、外的分離とは、違うクラスタにある対象は類似していないという性質をいい、内的結合とは、同じクラスタ内の対象は互いに類似しているという性質をいう。さらに、「類似」の尺度として、集合の要素間の距離が定義される。一般には、いわゆる「距離の公理」を満たすように距離が定義され、距離としては、ユークリッド距離、マハラビノス距離、シティブロック距離、ミンコフスキー距離などが使われることがある。
また、一般に、教師なし学習によりクラスタリングを行う手法としては、非階層型手法である「k-means法」などのように、クラスタの良さを定義する目的関数を最適にする分割を探索する方法である「分割最適化クラスタリング」や、階層的クラスタリング手法である「凝集型階層的クラスタリング」や「分割型階層的クラスタリング」などと呼ばれる手法が知られている。
ただし、このような従来のクラスタリング手法は、特徴量をすべて使って対象をクラスタ(グループ)に分け、得られるクラスタの分け方は一通りとなる、という特徴がある。
従って、特徴量に依存して複数のクラスタの分け方が存在する場合、うまく対処できないという問題が存在することになる。一般に、特徴量の数が多くなるほど、このような複数のクラスタ構造が存在する可能性が高いと考えられる。
なお、クラスタリングの手法としては、上記のような手法だけではなく、むしろ、クラスタリングを実施する複数の対象が、それぞれのクラスタにおいて、ある確率分布に従って生起したものと想定して、このような「確率分布」を推定するという方向で、クラスタリングを実行するというアルゴリズムもあり、このようなクラスタリング手法としては、たとえば、「混合ガウス分布によるクラスタリング法」などが知られており、より柔軟なクラスタリングを実行することが可能であることが知られている。
(特徴量の選択に応じた異なるクラスタリング)
以下では、まずは、特徴量に依存して複数のクラスタの分け方が存在する場合を説明するために、クラスタリングを実行する複数の対象を含んだ対象群について、各対象は、複数の特徴量で特徴づけられるものとする。
図20は、特徴量に依存して複数のクラスタの分け方が存在する場合を説明するための概念図である。
図20に示すように、クラスタリングの対象となるデータ(以下、単に「対象」と呼ぶ)が、6つの文字「A」「B」「C」「D」「E」「F」であるものとする。
そして、これらの文字は、異なる背景パターンや、異なる書体(文字スタイル)を有しているものとする。
そこで、これらの文字を特徴づける特徴量としては、「背景パターン」「文字スタイル」「文字に含まれる穴の数(線で完全に囲まれた領域の数)」を考えることができる。
したがって、同一の文字の集合を考えた場合でも、いずれの特徴量に基づいて、クラスタリングするかによって、異なるクラスタに分けられることになる。
図20では、たとえば、「背景パターン」に基づく場合、{A,D}{B,E}{C,F}の3つのクラスタに分けられ、「文字スタイル」に基づく場合、{A,B,C}{D,E,F}の2つのクラスタに分けられ、「穴の数」に基づく場合、0個、1個、2個にそれぞれ対応して、{C,E,F}{A,D}{B}の3つのクラスタに分けられる。
図20では、1つのクラスタは、1つの特徴量により特徴づけられる場合を例としているが、一般には、1つのクラスタは、複数個の特徴量により特徴づけられる。
図21は、複数の特徴量により、複数の対象が特徴づけられる場合のクラスタリングの概念を説明するための概念図である。
まず、図21(i)に示すように、クラスタリングする対象が行方向に配置され、これらの対象を特徴づける特徴量が列方向に配置されている「データ行列」を考える。
図21(a)に示すように、対象をクラスタリングする(対象を複数の対象クラスタに分割する)と同時に、各対象クラスタと関連するように特徴量もクラスタリングする手法は、「共クラスタリング」と呼ばれ、例えば、以下の公知文献にその手法が開示されている。
公知文献:Madeira SC, Oliveira AL. Biclustering algorithms for biological data analysis: a survey. IEEE/ACM Transactions on Computational Biology and Bioinformatics (TCBB). 2004; 1(1):24±45. https://doi.org/10.1109/TCBB.2004.2
「共クラスタリング」では、図21(a)に示すように、データ行列の行、あるいは列を、入れ替えることで、すなわち、対象及び特徴量を類似度によって再配列することで、たとえば、(i,j)(i=1,2:j=1,2,3)で示されるクラスタブロックに分割する。
このとき、各クラスタに含まれる対象ついての生成モデル(確率モデル)を仮定して、観測データに対して、その尤度が高くなるように各確率モデルのパラメータを決定する。
このように、各クラスタについて、確率モデルが推定されると、特定の観測データ(テストデータ)について、当該データが、いずれのクラスタに属するかが判別できる(分類できる)ことになる。
図22は、多重クラスタリングおよび多重共クラスタリングを説明するための概念図である。
図21(a)に示した「共クラスタリング」では、「データ行列」の行と列を入れ替えることにより、ブロック構造のクラスタを生成したので、特徴量が、複数の特徴量クラスタに分割された場合に、各対象は、この複数の特徴量クラスタに共通に並んでいるものとしてクラスタリングがされていることになる。
ただし、特徴量を複数の特徴量クラスタに分割し、特徴量クラスタごとに対象も対象クラスタに分割することを想定する場合、特徴量クラスタごとに、対象クラスタ中の対象の並び方(ある対象クラスタ中に含まれる対象の並び方)も異なるものとした方が、より尤度の高い確率モデルが推定できることが想定される。
このような場合、特徴量クラスタごとに、対象の分割方法(対象のクラスタリング)が異なることに対応して、特徴量クラスタを、特に、「ビュー(視点)」と呼ぶ。
以上のようにして、特徴量の視点ごとに、異なる対象のクラスタリングを実行することを、図22(b)に示すように「多重クラスタリング」と呼ぶ。
さらに、各視点において、特徴量の列と対象の行を入れ替えてクラスタリングすることで、さらに、観測データに対する尤度の高い確率モデルが推定できる場合は、図22(c)に示すように、「多重共クラスタリング」と呼ぶ。
ここで、ビューが1つしかない場合や、複数のビューが存在する場合に、少なくとも1つのビューにおいて特徴量クラスタが1種類しかないときも含めて、「多重共クラスタリング」と呼ぶこととし、「共クラスタリング」および「多重クラスタリング」は、「多重共クラスタリング」の下位概念であるものとする。
なお、この実施の形態においては、単に「クラスタリング」とい用語は、1つのビューにおけるクラスタの組を生成することを言うものとし、たとえば、図22(b)のように特徴量のビューへの分割と対象のクラスタリングを行う場合はを「多重クラスタリング」と呼び、図22(c)のようにビューへの分割と共クラスタリングとを同時に実行する場合は、「多重共クラスタリング」と呼ぶことで区別する。
図23は、「多重共クラスタリング」において、1つのビューの中に、異なる種類の確率分布の確率モデルが想定されている場合を示す概念図である。
図23(d)においては、白のブロックと、網掛のブロックでは、異なる種類の確率分布の確率モデルに従っている場合を示している。
たとえば、白のブロックは、連続的な確率変数の確率分布であるのに対して、網掛の部分は、離散的な確率変数の確率分布を想定している場合などを示す。
以後に説明するように、本実施の形態の「多重共クラスタリングの学習方法」では、このように、異なる分布を含む分布族に対して、クラスタリング処理を実行することが可能である。
図24は、多重共クラスタリングの学習方法の概要を説明するためのフローチャートである。
多重共クラスタリングの学習方法の処理が開始されると(S300)、クラスタリング分類器生成部3010は、データ行列について、ランダムに特徴量を部分群に分割して、特徴量のビューとビュー内の特徴量クラスタを生成する(S302:後述するYの生成(Yの初期化)に対応)。
続いて、クラスタリング分類器生成部3010は、ステップS302で生成された特徴量のビューおよび特徴量クラスタに対応して、対象クラスタの分割を生成し最適化する(ステップS304:後述するZの生成に対応)。
さらに、クラスタリング分類器生成部3010は、得られた対象クラスタに対して、特徴量の分割を最適化する(S306:後述するYの生成処理に対応し、生成されたZを使ってYを最適化する)。
続いて、クラスタリング分類器生成部3010は、目的関数が所定の条件を満たして収束しているかどうかを判定する(S308)。なお、この目的関数は、後述するような関数L(q(φ))に相当する。関数L(q(φ))は、これも後述するYおよびZをアップデートするに従って単調に増加する性質があり、その増え方が十分小さくなったと判断された時に収束したものと判定することになる。クラスタリング分類器生成部3010は、収束していなければ(S308でN)、処理をステップS304に復帰させ、収束していれば(S308でY)、処理を次のステップに進める。
そして、クラスタリング分類器生成部3010は、目的関数の大きさを記憶装置2080に格納する(S310)。
次に、クラスタリング分類器生成部3010は、ステップS302~S310までの処理を、所定回数実施したかを判定する。クラスタリング分類器生成部3010は、所定回数実施していない場合は(S312でN)、処理をステップS302に復帰させ、所定回数実施している場合は(S312でY)、処理を次のステップに進める。
クラスタリング分類器生成部3010は、目的関数が最大となる特徴量分割、クラスタの分け方を最終結果として(S314)、多重共クラスタリングについての学習の処理を終了させ、クラスタリング分類器を生成する。
(多重共クラスタリングの処理の詳細)
以下では、図24で説明した多重共クラスタリングの学習方法について、さらに詳しく説明する。
なお、多重共クラスタリングの処理の詳細については、下記の文献に開示があるので、以下では、その概要について説明する。
公知文献:Tomoki Tokuda, Junichiro Yoshimoto, Yu Shimizu, Go Okada, Masahiro Takamura, Yasumasa Okamoto, Shigeto Yamawaki, Kenji Doya,“Multiple co-clustering based on nonparametric mixture models with heterogeneous marginal distributions”, PLOS ONE | https://doi.org/10.1371/journal.pone.0186566 October 19, 2017
図25は、図24の多重共クラスタリングの学習方法におけるベイズ推定のグラフ表現を示す図である。
多重共クラスタリングモデルは、図25のグラフィカルモデルにまとめられており、このグラフィカルモデルは、関連するパラメータとデータ行列間の因果関係のリンクを明確にしている。
(多重共クラスタリングモデル)
特徴量(脳機能結合値)と被験者(ここでは、患者群の被験者)とは、図21(i)に示したようなデータ行列として表現されるものとする。
そして、データ行列Xは、事前に知られているM個の分布からなる分布族で構成されていると仮定する。
分布族に属する確率分布としては、ガウス分布、ポアソン分布、およびカテゴリ分布/多項分布などが含まれうるものとする。
クラスタリング分類器生成部3010は、X(m)に対して、各データサイズが n× d(m)であるように、以下のように分割する。
X= {X(1),…,X(m),…,X(M) }
ここで、mは、分布族(m=1,…,M)を示す指標である。さらに、ビュー(視点)の数を V(すべての分布族に共通)とし、ビューvおよび分布族mの特徴量クラスタの数をGν (m)とし、ビューvのオブジェクトクラスタの数をKv(すべての分布族に共通) で示す。
さらに、表記の簡単化のために、特徴の数およびクラスタの数を示すために、空のクラスタが存在することを許して、 G(m)=maxvv (m) および K=maxvvと表記する。
この表記では、分布族mに対する独立同分布(i.i.d.)のd(m)次元ランダムベクトル X1 (m), …, Xn (m)について、d (m) × V × G(m)の(3階の)特徴量分割テンソルY(m)を考えることとし、分布族mの特徴量jがビューvの特徴量クラスタgに属するときは、Yj,v,g (m)=1(それ以外は0)であるものとする。
異なる分布族に対してこれを組み合わせて、Y= {Y(m)}mとする。
同様にして、n×V×Kの オブジェクト分割 (3階) テンソルZを考え、オブジェクトiがビューvのオブジェクトクラスタkに属する場合は Zi、v、k = 1であるとする。
特徴量jはビューのうちの1つ( Σv,gj, v, g (m) =1)に属し、オブジェクトiは各ビューに属している (すなわち、Σki,v,k (m)=1)。さらに、Zはすべての分布族に共通しており、これは、推定される確率モデルは、すべての分布族の情報を使用して被験者クラスタリング解を推定していることを意味する。
まず、図25に示されるように、Yの事前生成モデルに対しては、ビューと特徴量クラスタの階層構造を考え、最初にビューが生成され、次に、特徴量クラスタが生成されるものとする。したがって、特徴量はビューと特徴量クラスタのペアのメンバーシップで分割されることとなり、特徴量の分割の割り当ては、ビューと特徴量クラスタで共同的に決定される。
一方、図25に示されるように、オブジェクト(対象)は各ビューのオブジェクトクラスタに分割されるため、Zに対して、オブジェクトクラスタの1つの構造だけを考えている。これらの生成モデルは、以下に説明するように、すべて「棒折り過程」(SBP:Stick Breaking Process,スティックブレイキング過程)に基づいていると仮定する。
(特徴量クラスタ Yの生成モデル)
j.. (m)が、階層棒折りプロセスによって生成される、分布族mの特徴量jのビュー/特徴量クラスタメンバーシップベクトルを示すものとすると、以下の式が成り立つ。
ここで、τ(m) は、1 × GVベクトル (τ1 ,1 (m) ,…,τG,V (m))T を示す。 (上付きTは転置行列を示す)。
Mul(・|π)は確率パラメータπを持つ1つのサンプルサイズの多項分布である。
Beta(・|a,b)は、事前サンプルサイズ(a,b)を有するベータ分布である。
j.. (m)は1×GVベクトル(Yj,1,1 (m), … ,Y j,V,G (m))Tを示す。
ここでは、所定の条件に従って、十分な大きなVのビューの数とGの特徴量クラスタの数を切り捨てる。 この処理については、たとえば、以下の文献を参照。
公知文献:Blei DM, Jordan MI, et al. Variational inference for Dirichlet process mixtures. Bayesian analysis. 2006; 1(1) : 121-143,
https://doi.org/10.1214/06-BA104
j,v,g (m)=1の場合、特徴量jはビューvの特徴量クラスタgに属する。デフォルトでは、ハイパーパラメータである集中度パラメータα1とα2は、1に設定される。
(オブジェクトクラスタ Zの生成モデル)
ビュー vのオブジェクトiの被験者クラスタメンバーシップベクトルであって、Zi, v.と表記されるベクトルが以下の式によって生成される。
ここでZi, v.は、Zi,v= (Zi, v, 1,…, Zi, v, KTによって与えられる1×K(Kは十分に大きい値を取る)ベクトルである。集中度パラメータβは、1に設定される。
(尤度と事前分布)
各インスタンスXi,j (m) は、YとZに条件付きで、独立に特定の分布に従うと仮定する。ビューv、特徴量クラスタg、およびオブジェクトクラスタkのクラスタブロックにおける分布族mのパラメータを θv,g,k (m) と表す。
さらに、Θ={θv,g,k (m)} v,g,k,mと表記して、Xの尤度の対数は以下の式に従う:
ここで、I(x) は指標関数で、xが真の場合は1を返し、それ以外の場合は0を返す。
尤度は、w={wvv、w´={w´g,v (m)g,vおよび u={uk,vk,vに直接関連付けられていない。
未知変数の結合事前分布φ={Y,Z,w,w´,u,Θ}(すなわち、クラスのメンバーシップ変数とモデルパラメータ)は、次のように与えられる。
(変分推定)
YとZのMAP(最大事後)推定には変分ベイズEMアルゴリズムを使用する。
このような変分ベイズEMアルゴリズムについては、以下の文献に開示がある。
公知文献:Guan Y, Dy JG, Niu D, Ghahramani Z. Variational inference for nonparametric multiple clustering. In:MultiClust Workshop, KDD-2010; 2010.
対数周辺尤度p(X)はジェンセンの不等式を使用して以下のように近似される。
なお、ヤンセンの不等式については、以下の文献に開示がある。
公知文献:Jensen V. Sur les fonctions convexes et les inegalites entre les valeurs moyennes. Acta Mathematica. 1906; 30(1):175-193.
https://doi.org/10.1007/BF02418571
ここで、q(φ)はパラメータφの任意の分布である。q(φ)とp(φ)の間のカルバック・ライブラー・ダイバージェンス、すなわち、KL(q(φ), p(φ|X))によって左辺と右辺の差が与えられることが証明される。したがって、q(φ)を選択するアプローチは、KL(q(φ), p(φ|X))を最小化することであり、これは評価するのは、一般には難しい。
ここでは、異なるパラメータ(平均場近似)に対して因数分解されるq(φ)を選ぶ。
ここで、各 q(・)は、パラメータのサブセットwv, w´g,v (m),Yj.. (m), uk, v, Zi,v.およびθv,g,k (m)に対して、さらに因数分解される。
一般に、KL(Πl=1 L qll), p (φ|x))を最小化する分布 q ii)は、以下の式によって与えられる。
この性質については、以下の文献に開示がある。
公知文献:Murphy K. Machine Learning: A Probabilistic Perspective. Cambridge, Massachusetts: MIT Press;2012.
この性質を現在検討しているモデルに適用すると、以下を示すことができる。
ここで、以下の形で表される関数を考える。
上記の式を除くハイパーパラメータは、以下の式で表される。
ψ(・)は、ガンマ関数の対数の1次導関数として定義されるダイガンマ関数を示す。
τj,g,v (m) は、各ペア (j,m)についての複数のペア(g,v) にわたって正規化される。一方、ηi, v, kは 各(i, v) のペアでk について正規化される。
観測モデルとパラメータΘの事前分布は、後述する。
(観測モデル)
観測モデルでは、ガウス分布、ポアソン分布、およびカテゴリ分布/多項分布を考慮する。クラスタブロックごとに、クラスタブロック内の特徴が独立していることを前提に、これらのファミリーの1変量分布をフィッティングさせる。これらの分布族のパラメータは、共役事前分布を仮定する。
(最適化アルゴリズム)
ハイパーパラメータの更新方程式では、変分ベイズEMアルゴリズムを用いて、次のように計算が実行する。
まず、ランダムに{τ(m)}mと{ηv}vを初期化し、式(1)の下限L(q(φ)) が収束するまでハイパーパラメータを更新する。これは、L(q(φ))の観点から局所的に最適な分布q(φ)を生成する。この手順を何度も繰り返し、近似事後分布q*(φ)として最大の下限を持つ最良の解を選択する。
YとZのMAP推定値は、それぞれ argmaxY q* Y(Y)およびargmaxZ q* Z(Z)として評価される。
下限L(q(φ))は以下の式で与えられる。
右辺の両方の項は、閉じた形式で導出できる。q(φ)が最適化されるにつれて、この値は、単調に増加することが示される。すなわち、上述した通り、関数L(q(φ))は、YおよびZをアップデートするに従って単調に増加する性質があり、その増え方が十分小さくなったと判断された時(特に限定されないが、たとえば、増分が所定値以下となるなどの条件が成り立つ時)に収束したものと判定することになる。
まず、各特徴の分布族が特定され、対応する分布族のデータ行列を生成する。次に、データ行列のセットに対して、さらに、YとZのMAP推定値を生成し、YおよびZの推定値を使用して、各ビューにおけるオブジェクト/特徴量クラスタが分析される。
(モデル表現)
多重共クラスタモデルは、ビューの数と特徴/オブジェクトクラスタの数がデータ駆動のアプローチで導き出されるため、さまざまなクラスタリング モデルを表現するのに十分な柔軟性を備えている。たとえば、ビューの数が1の場合、モデルは、共クラスタモデルと一致する。特徴量クラスタの数がすべてのビューに対して1つである場合、多重クラスタリングモデルと一致する。さらに、ビューの数が1であり、特徴量クラスタの数が特徴の数と同じである場合、独立した特徴を持つ従来の混合モデルと一致する。さらに、このモデルは、オブジェクトクラスタを区別しない無情報特徴を検出することができる。このような場合、モデルはオブジェクトクラスタの数が 1つのビューを生成する。モデルの利点は、このような基になるデータ構造を自動的に検出することである。
以上のような「多重共クラスタリング法」により、以下のことが可能となる。
1)データの背後にある複数のクラスターの分け方(対象の分け方だけでなく、特徴量の分け方も含める)とそれに対応する特徴量群を、データ駆動により同定することができる。
2)この手法によって他の手法では見出せなかったクラスターを同定することが可能となる。
3)さらに、それぞれのクラスターの分け方を特徴量によって意味づけることができ、個々のクラスターの解釈を容易にできる。
[データセットに対するクラスタリングの結果の評価]
以下では、上述したようなSRPBSとして公開されている多数の被検者から収集された多施設大規模fMRIデータを2つに分割して、それぞれを上述したような多重共クラスタリング法で使用することで、クラスタリングについての汎化性能を検証する。
図26は、このように2つに分割されたデータセット1とデータセット2とを示す図である。
データセット1は、施設1~4までで取得された健常者545人、うつ患者138人のデータからなり、データセット2は、施設5~8までで取得された健常者263人、うつ患者181人のデータからなる。基本的には、図10および図11に示したデータセットに対応する。
図27は、各データセットにおいてクラスタリングを実施する概念を説明する概念図である。
図27に示すように、データセット1に対しては、それぞれ、独立に、図24に示したフローに従って多重共クラスタリング法により、クラスタリングが実行される。
ここで、関心の対象となるのは、このようにして、データセット1と、データセット2に対して、それぞれ個別に実行されたデータ駆動によるクラスタリング法により得られた各クラスタが、相互にどの程度、類似しているか(一致の程度は、どの程度か)という問題である。
データセット1と、データセット2におけるクラスタリングが、同一または類似する特性を有するクラスタ(被験者グループ)へ分類(グループ分け)できているのであれば、このようなデータ駆動によるクラスタリングは、施設や計測装置等に依存することなく、汎化性の高い状態で実行されていることになる。そこで、問題となるのは、「同一または類似する特性を有するクラスタへ分類」されていることを、どのように定量的に評価するのかということになる。
図28は、被験者データに対する多重共クラスタリングの例を示す概念図である。
図28(a)に示すように、入力となるデータ行列として、行方向に被験者が並び、列方向に特徴量が並んでいるものとする。
この入力のデータ行列に対して、多重共クラスタリングを実行すると、たとえば、図28(b)に示すように、2つのビューに特徴量が分割され、各ビューにおいて、被験者がクラスタリングされることになる。
図29は、データセット1およびデータセット2に対して、実際に多重共クラスタリング処理を実施した結果を示す図である。
図29においては、クラスタリングのための特徴量としては、データセット1およびデータセット2とも、それぞれ99個が選択されている。
その上で、データセット1については、138人のうつ患者に対して、データセット2については、181人のうつ患者に対して、多重共クラスタリング処理が実行されている。
データセット1については、ビュー1とビュー2の2つのビューに特徴量が分割される。ビュー1については、さらに2つの特徴量クラスタに共クラスタリングされ、被験者は5つの被験者クラスタに分割され、ビュー2についても、被験者は5つのクラスタに分割される。
データセット2についても、ビュー1とビュー2の2つのビューに特徴量が分割される。ビュー1については、さらに2つの特徴量クラスタに共クラスタリングされて、被験者は4つの被験者クラスタに分割され、ビュー2については、被験者は5つのクラスタに分割される。
図30は、データセット1とデータセット2において、それぞれのビューに割り当てられた脳機能結合(FC)の個数を示す表である。
データセット1では、ビュー1に特徴量としてFCが92個、ビュー2にFCが7個割り当てられる。
データセット2では、ビュー1に特徴量としてFCが93個、ビュー2にFCが6個割り当てられる。
また、この表では、データセット1とデータセット2で、割り当てられた脳機能結合のうち、一致する個数を表の対角上に記載している。データセット1とデータセット2で、ビュー1とビュー2に割り当てられた脳機能結合は、ほぼ一致していることがわかる。
(クラスタリング(層別化)の汎化性(データセット間類似性)検証方法)
以下では、データセット1と、データセット2に対して、それぞれ個別に実行されたデータ駆動によるクラスタリングにより得られた各クラスタが、相互にどの程度、類似しているか(一致の程度は、どの程度か)を定量的に評価する。
図31は、このようなクラスタリングの類似度(層別化の汎化性能)の評価手法を説明するための概念図である。
まず、図31(a)に示すように、データセット1とデータセット2を独立に、上述したような多重共クラスタリング法により、クラスタに分割した場合、被験者は、それぞれのデータセットのクラスタにおいて、相互に独立であるために、クラスタリングの類似度を比較する、ということは困難である。
ここで、データセット1で作成された分類器1で、データセット1の被験者の分類を実行した結果をクラスタリング1とする。一方で、データセット2で作成された分類器2で、データセット2の被験者の分類を実行した結果をクラスタリング2とする。
これに対して、図31(b)に示すように、データセット1で作成された分類器1で、データセット2の被験者の分類を実行した結果をクラスタリング1´とする。一方で、データセット2で作成された分類器2で、データセット1の被験者の分類を実行した結果をクラスタリング2´とする。
この場合、クラスタリング1とクラスタリング1´の間、および、クラスタリング2とクラスタリング2´の間では、それぞれ、共通な被験者について分類を行っているので、それぞれの類似度を評価することができる。
(クラスタリング間の類似度(再現度)を測る評価尺度)
ここで、データ駆動によりクラスタリング処理を実行しており、図29におけるクラスタインデックスの値自体(インデックスの値の順序)には、意義はないので、同じデータ集合に対して、異なるクラスタリングが実行された場合に、その類似度をどのように評価するかが問題になる。
たとえば、同じデータ集合Xに対する2つのクラスタリング結果πとρとがある場合において、この2つのクラスタリング結果の類似性(外的妥当性尺度)を評価する尺度として、ランド指数(Rand index)が知られている。
データ集合中のすべてのデータの対(ペア){x1,x2}∈X(M=N(N-1)/2個)については、対の種類として、以下のものがあり、それぞれの種類に属するペアの個数を以下のように定義する。
このとき、2つのクラスタリングにより分類されるクラスタが同じクラスタかどうかの判定の正解率として、ランド指数が以下の式により定義される。
ただし、たとえば、データ集合の各クラスタの要素数に偏りがある場合などでは、「ランダムにクラスタリングした場合でも、ランド指数が高い値となってしまう」場合があることが知られている。このため、より厳密には、以下のような調整ランド指数(ARI:Adjusted Rand Index)が使用される。
図32は、ARIについて説明するための概念図である。
なお、ARIについては、たとえば、以下の文献に開示がある。
公知文献:Jorge M. Santos and Mark Embrechts,“On the Use of the Adjusted Rand Index as a Metric for Evaluating Supervised Classification”,ICANN 2009, Part II, LNCS 5769, pp. 175-184, 2009.
上述のとおり、同じデータ集合に対して、2つのクラスタリング結果をそれぞれ適用する場合、図32(a)に示すように、2回とも同じクラスタに分類される場合と、2回とも別のクラスター分類される場合がある一方、図32(b)に示すように、1回は同じクラスタに分類されるが、もう1回は異なるクラスタに分類される場合も存在する。
ARIは、2つのクラスタリングが互いに独立だとした場合に、データ対がともに2つのクラスタリングで同一のクラスタに分類され、または、ともに異なるクラスタに分類される場合の期待値を計算して、その期待値をランド指数の分子と分母のそれぞれから減算して、計算される。したがって、ARIでは、クラスタリングに相関がない場合に、その値が0となるように調整されていることになる。
ここでは、Aは、2つのクラスタリングに対して、[(2回とも同じクラスタに分類)+(2回とも別のクラスタに分類)]された被検者ペアの数を表し、max(A)は、全ペアの個数を表し、Eは、2回のクラスタリングが独立にもかかわらず、割り当て結果が一致する被検者ペアの数を表す。
図33は、クラスタリング1とクラスタリング1´の間、および、クラスタリング2とクラスタリング2´の間の類似度の評価結果を示す図である。
図33(a)は、データセット1とデータセット2について、それぞれのビューに対するARIを算出した表である。
データセット1とデータセット2について、ビュー1については、ARI=0.47、ビュー2については、ARI=0.51となっており、有意な類似性があるといえる。
図33(b)は、図33(a)に対応するパーミュテーションテスト(Permutation Test)の結果を示す。ビュー1およびビュー2については、要素の交換を行ったとき(ヒストグラムで表示)と比較して、ビュー1とビュー2では、ARI値(実線で示す)が統計的に有意に高い値となっていることがわかる。
なお、「パーミュテーションテスト」とは、被験者のクラスタ属性ラベルを被験者間でランダムに交換した場合においてARI値を算出した結果であり、図中では、このような交換を所定回数実施したときのその分布がヒストグラムとして表示されている。もしも、クラスタリング間の類似が統計的に有意であるならば、要素をランダムに交換した場合に比べて、比較対象となるクラスタリング間のARI値が有意に高い値となることになる。
以上より、データセット間でのクラスタリング(層別化)について、有意な類似性が認められた、すなわち、汎化したクラスタリングが実現されたと判断できる。
以上説明した通り、データセット1に対する多重共クラスタリングも、データセット2に対する多重共クラスタリングも、ともに、データ駆動により実現されたものであるから、脳機能結合を特徴量とする「患者の層別化」の基礎となるといえる。
図34は、クラスタリング1とクラスタリング1´において、それぞれビュー1の各クラスタに割り当てられる被験者数の分布を示す表である。
被験者クラスタインデックスを適切に並べ替えることにより、ほとんどの被験者が、表の対角の近傍に分布するようにでき、2つのクラスタリングが、相互に類似していることが視覚的にも確認できる。
[ハーモナイゼーション処理]
以下では、以下の文献に開示され、図6において、ハーモナイゼーション処理と呼んでいる処理の内容について、説明する。
[トラベリングサブジェクト法のハーモナイゼーション]
以下では、以上説明してきた「疾患識別器」の生成や、層別化のための「クラスタリング処理」のために使用される手法であって、標本バイアスとは独立に、測定バイアスを評価して計測データをハーモナイズするための手法について説明する。
図35は、本実施の形態のrs-fcMRI法において、サイト間を移動しながら計測を受ける移動被験者(以下、「トラベリングサブジェクト:旅行被験者」)によるサイト間差の評価手法を説明するための概念図である。
以下に説明するように、本実施の形態では、トラベリングサブジェクトのデータセットを使用して、測定バイアスだけを除くことができるようなハーモナイゼーション法について説明する。
図35を参照して、計測サイトMS.1~Ms.Nsのサイトを横断して測定バイアスを評価するために、トラベリングサブジェクトTS1(人数:Nts人)のデータセットを取得する。
健康なNts人の参加者の安静時脳活動が、Ns個のサイトの各々で撮像されるものとして、Ns個のサイトは、患者データを撮像したすべてのサイトを含んでいるものとする。
取得されたトラベリングサブジェクトのデータセットは、データセンター200における記憶装置210に移動被験者データとして格納される。
そして、後述するように、計算処理システム300において、「脳活動バイオマーカーのハーモナイズ方法」のための処理が実行される。
トラベリングサブジェクトのデータセットは健常者群だけを含む。また、参加者はすべてのサイトにわたって同じとする。したがって、トラベリングサブジェクトについては、サイト間差は「測定バイアス」のみから成る。
以下に説明する本実施の形態のハーモナイゼーション法では、「脳活動バイオマーカーのハーモナイズ方法」として、各計測サイトにおける計測データについて、「測定バイアス」の影響を除いて補正する処理を行う。
つまり、以下では、「測定バイアス」および「標本バイアス」を、「統計モデリング」手法のなかの「一般化線形混合モデル(GLMM:Generalized Linear Mixed Model)」を用いて評価する。
ここで、通常は、GLM(Generalized Linear Model:一般化線形モデル)は、「応答変数」の確率分布を説明する「説明変数」を組み込んだモデルである。GLMには、「確率分布」「リンク関数」「線形予測子」という主に3つの部品があり、この部品の組み合わせかたを指定することによって、さまざまなタイプのデータを表現できる。
さらに、GLMM(一般化線形混合モデル)は、GLMでは説明できない「人間が測定できない・測定しなかった個体差など」を組み込むことができる統計モデルである。GLMMでは、たとえば、対象がいくつかの部分集合(たとえば、計測場所が異なる部分集合)からなるような場合も、その場所差をモデルに組み込むことができる。別の言い方をすれば、複数の確率分布を部品とする(混合させた)モデルということになる。
たとえば、GLMMについては、以下の文献に開示がある。
公知文献8:久保 拓弥著、「データ解析のための統計モデリング入門」、岩波書店、2012年 第1刷、2017年 第14刷
ただし、以下に説明する本実施の形態の統計モデルでは、通常、「効果」と呼ばれるものについて、「バイアス」と「要因(ファクター)」という用語を、「バイアス」に関しては「測定バイアス」と「標本バイアス」に使用し、「要因(ファクター)」に関しては、その他の要因(被験者要因や疾患要因)に関して使用する。
そして、以下の解析は、単純なGLMMの流れとは異なり、「固定効果」と「ランダム効果」という分け方をせずに要因の解析を行う。これは、通常、GLMMを使用すると、ランダム効果に関しては分散のみが推定され、各要因の効果の大きさがわからなくなるからである。そこで、以下では、各要因の効果の大きさを評価するために、各要因について、平均0の固定効果となるように、以下のように変数を変換して、推定を行う。
i) 各々の機能的なコネクティビティに対する相関値のすべてのサイトにわたる平均からの偏差として、各サイトの測定バイアスを定義する。
ii) 健常者および精神疾患の患者の標本バイアスは、互いに異なると仮定する。したがって、健常者群及び各疾患を有する患者群に対して個別に、各サイトの標本バイアスを計算する。
iii) 疾患要因は、健常者群の値からの偏差として定義する。
すなわち、以下では、患者が含まれるデータセットと、トラベリングサブジェクトのデータセットに対して、一般化線形混合効果モデルを以下のように適用する。
トラベリングサブジェクトが、Nts人であり、Ns個の計測サイトのうち、健常者の計測が行われたサイト数をNsh個とし、ある疾患(ここでは、添え字“dis”で表現する)の患者の計測が行われたサイトの個数をNsd個とする。
参加者要因(p)、測定バイアス(m)、標本バイアス(Shc, Sdis)および精神疾患要因(d)は、患者に対する計測結果のデータセットおよびトラベリングサブジェクトのデータセットからのすべての参加者の機能的なコネクティビティの相関値に、回帰モデルを適合させることにより評価される。
以下では、ベクトルは小文字(例えば、m)によって表示され、ベクトルはすべて列ベクトルであると仮定される。
ベクトルの要素は、mkのように添字によって表示される。
脳領野間のn個の相関値からなる機能的コネクティビティベクトル(列ベクトルとする)の回帰モデルは、以下の式のように表される。
参加者の特性を表わすために、1-of-Kのバイナリーコード体系を使用し、サイトkに属する測定バイアスmに対するターゲットベクトル(たとえば、xm)は、1に等しい要素k以外は、全てゼロに等しい。
参加者がどんなクラス(健常者、患者、トラベリングサブジェクト)にも属さなければ、ターゲットベクトルは、すべての要素が0に等しいベクトルである。
上付きの添字Tは、マトリックスかベクトルの置換を表示し、xTは行ベクトルを表わす。
ここで、mは、測定バイアス(Ns×1の列ベクトル)を表し、shcは、健常者群の標本バイアス(Nsh×1の列ベクトル)を表し、sdisは、患者の標本バイアス(Nsd×1の列ベクトル)を表し、dは、疾患要因(2×1の列ベクトル、健常および疾患を要素とする)、pは、参加者要因(Nts×1の列ベクトル)を表し、constは、すべての計測サイトからの全ての参加者(健常者、患者、トラベリングサブジェクトを含む)にわたる機能的コネクティビティの平均を表し、e~N(0,γ-1)は、ノイズを表す。
なお、ここでは、説明の簡略化のための疾患の種類は1種類であるものとして説明している。疾患の種類が複数ある場合については、後述する。
各機能的なコネクティビティの相関値に対しては、回帰モデルの計画行列がランク不足であるので、L2正規化による最小二乗回帰を使用して、それぞれのパラメータを評価する。なお、L2正規化による最小二乗回帰法以外にも、たとえば、ベイズ推定法など別の評価方法を用いることも可能である。
上記のような回帰計算の後には、被験者aのb番目のコネクティビティは、以下のように記述できる:
図36は、被験者aのb番目の機能的なコネクティビティの表現を説明するための概念図である。
図36においては、1項目と2項目のターゲットベクトルの意義と、測定バイアスベクトルおよび健常者の標本バイアスベクトルについて示している。
第3項目以降についても同様である。
(ハーモナイゼーション処理のフロー)
図37は、ハーモナイゼーションのために、測定バイアスを算出する処理を説明するためのフローチャートである。
まず、データセンター200の記憶装置210には、各計測サイトから被験者(健常者および患者)のfMRI計測データ、被験者の属性データ、測定パラメータが収集される(図37 S402)。
続いて、トラベリングサブジェクトTS1の脳活動が、特に限定されないが、たとえば、所定の周期(たとえば、1年周期)で、各計測サイトを巡回して計測され、データセンター200の記憶装置210には、各計測サイトからトラベリングサブジェクトのfMRI計測データ、被験者の属性データ、測定パラメータが収集される(図37 S404)。
ハーモナイゼーション算出部3020は、上述したようなGLMM(一般化線形混合モデル)を用いることで、機能的コネクティビティについて、各計測サイトの測定バイアスを評価する(図37 S406)。
このようにして算出された各計測サイトの測定バイアスを、ハーモナイゼーション算出部3020は、記憶装置2080に測定バイアスデータ3108として格納する(図37 S408)。
(判別器生成処理におけるハーモナイゼーション)
判別処理部3000が、被験者に対する疾患あるいは健常ラベルについての疾患識別器を生成する処理における脳機能結合値に対するハーモナイゼーションについて簡単に説明する。
このような疾患識別器は、被験者を診断するための補助情報(支援情報)を提供する。
相関値補正処理部3004は、記憶装置2080に格納されている計測サイトごとの測定バイアスデータ3108を読出し、疾患識別器生成の機械学習のための訓練対象となる各被験者の相関行列の非対角成分について、以下の式のようにして、ハーモナイゼーション処理を実行する。
ここで、機能的コネクティビティConnectivityは、ハーモナイゼーション前の機能的コネクティビティベクトルを表し、Csubは、ハーモナイゼーション後の機能的コネクティビティベクトルを表す。また、m(ハット)(以下、文字xの頭部に^が付されたものを“x(ハット)”と表記する)は、上述したようなL2正規化による最小二乗回帰により評価された計測サイトにおける測定バイアスを表す。これにより、機能的コネクティビティConnectivityが計測された計測サイトに対応する測定バイアスが、機能的コネクティビティConnectivityから減算され、ハーモナイゼーション処理を受けることになる。
補正処理が実行された後のデータは、補正後相関値データ3110として、記憶装置2080に格納される。
計測サイト間のバイアスを除去する方法としては、上述したようなトラベリングサブジェクトによる方法に限られるものではない。たとえば、これもComBat法など他の方法を使用してもよい。
[実施の形態2]
実施の形態1では、脳活動計測装置(fMRI装置)で、複数の測定場所で計測された脳活動データを計測し、この脳活動データに基づいて、バイオマーカーの生成およびバイオマーカーによる診断ラベルの推定(予測)を行う構成として、分散処理により行う一例の構成を説明した。
ただし、i)バイオマーカーを機械学習により訓練するための脳活動データの計測(データ収集)、ii)バイオマーカーの機械学習による生成処理および特定の被験者(推定の対象であり、以下「第1の被験者」とも呼ぶ)についてのバイオマーカーによる診断ラベルの推定(予測)の処理(推定処理)、iii)上記特定の被験者についての脳活動データの計測(第1の被験者の脳活動計測)を、それぞれ、異なる施設で分散して実行する構成とすることも可能である。
図38は、データ収集、推定処理および対象者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。
図38を参照して、サイト100.1~100.Nは、患者群、及び健常者群(すなわち第2の被験者群)のデータを脳活動計測装置により計測する施設であり、管理サーバ200は、サイト100.1~100.Nsからの計測データを管理する。
計算処理システム300は、サーバ200に格納されたデータから識別器を生成する。
また、計算処理システム300のハーモナイゼーション算出部3020は、サイト100.1~100.NsおよびMRI計測装置410のサイトを含めて、ハーモナイゼーション処理を実行するものとする。
MRI装置410は、計算処理システム300上の識別器の結果を利用する別サイトに設けられており、特定の被験者について脳活動のデータを計測する。
コンピュータ400は、MRI装置410が設けられる別サイトに設置され、MRI装置410の測定データから特定の被験者の脳の機能結合の相関データを算出し、機能結合の相関データを計算処理システム300に送信して、返信されてくる識別器の結果を利用する。
サーバ200は、サイト100.1~100.Nsから送信されてくる患者群および健常者群のMRI測定データ3102と、MRI測定データ3102に関連付けられた被験者の人属性情報3104とを格納し、計算処理システム300からのアクセスに従って、これらのデータを計算処理システム300に送信する。
計算処理システム300は、通信インタフェース2090を介して、サーバ200からのMRI測定データ3102および被験者の人属性情報3104を受信する。
なお、サーバ200、計算処理システム300、コンピュータ400のハードウェアの構成は、基本的に、図5で説明した「データ処理部32」の構成と同様であるので、その説明は繰り返さない。
図38に戻って、相関行列算出部3002、相関値補正処理部3004、疾患識別器生成部3008、クラスタリング分類器生成部3010および判別値算出部3012、ならびに、機能結合の相関行列のデータ3106、測定バイアスデータ3108、補正後相関値データ3110および識別器データ3112については、実施の形態1で説明したのと同様であるので、その説明は、繰り返さない。
MRI装置410は、診断ラベルの推定対象となる被験者の脳活動データを計測し、コンピュータ400の処理装置4040は、計測されたMRI測定データ4102を不揮発性記憶装置4100に格納する。
さらに、コンピュータ400の処理装置4040は、MRI測定データ4102に基づいて、相関行列算出部3002と同様にして、機能結合の相関行列のデータ4106を算出し、不揮発性記憶装置4100に格納する。
コンピュータ400のユーザから診断の対象となる疾患が指定され、当該ユーザの送信の指示に従い、コンピュータ400が、機能結合の相関行列のデータ4106を計算処理システム300に送信する。これに応じて、計算処理システム300は、MRI装置410の設置されるサイトに対応したハーモナイゼーション処理を実行して、判別値算出部3010は、指定された診断ラベルについての判別結果やサブタイプについての評価結果を算出し、計算処理システム300は、通信インタフェース2090を介して、コンピュータ400に送信する。
コンピュータ400では、図示しない表示装置などを介して、ユーザに対して、判別結果を知らせる。
このような構成とすることで、より多くの被験者について収集したデータに基づいて、識別器による診断ラベルの推定結果を提供することが可能となる。
また、サーバ200と計算処理システム300とを別個の管理者が管理する形態とすることも可能で、その場合、サーバ200にアクセスできるコンピュータを制限することで、サーバ200に格納される被験者の情報のセキュリティを向上させることも可能となる。
さらに、計算処理システム300の運営主体からみると、「識別器による判別のサービスを受ける側(コンピュータ400)」に対して、識別器についての情報や「測定バイアス」に関する情報は、一切提供しなくても、「判別結果を提供するサービス」を行うことが可能となる。
なお、以上の実施の形態1および実施の形態2の説明では、脳機能画像法により脳活動を時系列に計測するための脳活動検出装置としては、リアルタイムfMRIを用いるものとして説明した。ただし、脳活動検出装置として、上述したfMRI、脳磁計、近赤外光計測装置(NIRS)、脳波計、またはこれらの組み合わせを使用することができる。たとえば、これらの組合せを用いる場合、fMRIとNIRSとは、脳内の血流変化に関連する信号を検出するものであり、高空間分解能である。一方で、脳磁計や脳波計は、脳活動に伴う電磁場の変化を検出するための高時間分解能であるという特徴をもつ。したがって、たとえば、fMRIと脳磁計とを組み合わせれば、空間的にも時間的にも高分解能で脳活動を計測することができる。あるいは、NIRSと脳波計とを組み合わせても、同様に空間的にも時間的にも高分解能で脳活動を計測するシステムを小型で携帯可能な大きさで構成することも可能である。
以上のような構成により、神経・精神疾患に対して、脳機能画像法によるバイオマーカーとして機能する脳活動解析装置および脳活動解析方法を実現することが可能となる。
また、以上の説明では、被験者の属性として「診断ラベル」を含んだ場合について、機械学習による識別器の生成により、当該識別器をバイオマーカーとして機能させる例について説明したが、本発明は、必ずしもこのような場合に限定されず、事前に機械学習の対象となる測定結果を得る対象の被験者群が、客観的な方法により、複数のクラスに分けられており、被験者の脳領域間(関心領域間)の活動度の相関(結合)を測定し、測定結果に対する機械学習により、クラスに対する識別器が生成できるものであれば、他の判別のために使用されるものであってもよい。
また、上述のとおり、このような判別は、ある属性に属する可能性を確率として表示するものであってもよい。
したがって、たとえば、ある「訓練」や「行動パターン」をとることが、被験者にとっての健康増進に役立つかを客観的に評価することができる。また、実際には、疾患に至っていない状態(「未病」)であっても、ある「食物」「飲料」等のような摂取物や、ある活動などが、より健康状態に近づくために効果があるのか、ということを客観的に評価することも可能である。
また、未病の状態においても、上述したように、たとえば、「健常である確率は、○○%」といような表示が出力されれば、ユーザに対して、健康状態について客観的な数値としての表示を行うことができる。このとき、出力されるのは、必ずしも確率でなくとも、「健康度合いの連続値、たとえば、健常である確率」をスコアに変換したものを表示することとしてもよい。このような表示を行うことで、本実施の形態の装置を、診断の支援以外にも、ユーザの健康管理のための装置として使用することも可能である。
[実施の形態3]
[治療法選択支援システム]
(治療法選択支援システムの構成)
本発明のある実施形態は、うつ症状を呈する被験者の脳活動の計測結果に基づいて、この被験者に対する治療法の選択に関する情報を提供するための治療法選択支援システム1000に関する。
図39は、治療法選択支援システム1000の機能構成を示す図である。
治療法選択支援システム1000は、医療機関4において、診断ラベルの推定対象となる被験者の脳活動データを計測するためのMRI装置410からのデータを受信可能なように接続されたコンピュータ400と、支援情報提供装置300aと、クラスタリング分類器の生成処理を実施する分類器生成装置300bと、データサーバ200´と、治療法情報提供サーバ500とを備える。
図39を参照して、データサーバ200´中の記憶装置210の健常者・患者データベースは、図示しないMRI装置100.1~100.Nが設定されているサイトから収集された、MRI測定データ3102および被験者の人属性情報3104を管理する。
支援情報提供装置300aは、i)医療機関4のコンピュータ400から伝送される被験者の脳活動の計測結果を入力として受け、演算装置2040aが、計測結果に対して判別値算出部3012として実行するクラスタリング分類器による分類結果に応じて、対応する治療法情報を出力する治療法情報生成部3200の処理を実行する点と、ii)疾患識別器生成部3008およびクラスタリング分類器生成部3010の処理が、支援情報提供装置300aとは異なる計算システムである分類器生成装置300b内の演算装置2040bにより実行される構成となっている点を除き、図38に示す計算処理システム300の構成に対応する。
なお、図38と同様に、支援情報提供装置300aと分類器生成装置300bとは、同一のコンピュータシステム上の機能として実装されてもよい。
分類器生成装置300bは、サーバ200´に格納されたデータを「学習データ」(または「発見コホートデータ」)として疾患識別器およびクラスタリング分類器を生成する。なお、サーバ200´に格納された(学習データ以外の)データは、疾患識別器およびクラスタリング分類器に対するバリデーションデータ(または「検証コホートデータ」)として使用することもできる。なお、特に限定されないが、図39の構成においては、サーバ200´の記憶装置210内の健常者・患者データベースには、MRI測定データ3102および被験者の人属性情報・測定パラメータ3104が格納されるだけではなく、このMRI測定データ3102に基づいて算出された機能結合の相関行列のデータ、ハーモナイゼーション処理により補正された相関値のデータが予め算出されて格納されているものとする。疾患識別器生成部3008およびクラスタリング分類器生成部3010は、この補正後の相関値に基づいて、学習処理およびバリデーション処理を実行するものとして説明する。なお、分類器生成装置300b自身が、機能結合の相関行列のデータの算出や、ハーモナイゼーション処理による相関値の補正の処理を実行することとしてもよい。
その他、基本的に、図38の構成と同一部分には、同一符号を付している。
なお、特に限定されないが、MRI装置410で取得されたデータを、記憶装置210の健常者・患者データベースに追加して、疾患識別器およびクラスタリング分類器を再学習する構成としてもよい。
また、支援情報提供装置300aのハーモナイゼーション算出部3020は、図示しないMRI装置100.1~100.Nsに対して、MRI計測装置410についてのハーモナイゼーション処理を実行するものとする。
MRI装置410は、支援情報提供装置300aのクラスタリング分類器の結果を利用する医療機関4に設けられており、特定の被験者について脳活動のデータを計測する。計測された特定の被験者の撮像データ(脳構造画像データおよび脳機能画像データを含む)については、コンピュータ400の匿名化処理部4042において匿名化処理が実行された後に、支援情報提供装置300aに送信される。
特に限定されないが、支援情報提供装置300aについては、いわゆるクラウドコンピュータを利用してもよい。また、これも特に限定されないが、コンピュータ400において、MRI装置410の測定データから特定の被験者の脳の機能結合の相関データを算出し、機能結合の相関データを支援情報提供装置300aに送信する構成としてもよい。この場合は、撮像データについては、機能結合の相関データへ変換されること自体が匿名化に相当する。
治療法情報生成部3200は、治療法情報提供サーバ500内の治療情報DB5100内のデータを治療法情報提供システム5200から受信して、クラスタリング分類器の分類結果に応じて、対応する治療法選択支援データを情報提示部4044に対して返信する。治療法情報生成部3200側で、分類結果と治療法選択支援データとの対応関係の情報を、予め治療法情報提供システム5200から受信して記憶しておく構成としてもよいし、あるいは、治療法情報生成部3200が、分類結果の情報を、都度、治療法情報提供システム5200に問合せ情報として送信して、治療法情報提供システム5200から問合せに対する回答として返信される治療法選択支援データを受信する構成としてもよい。
また、特に限定されないが、たとえば、治療法情報提供サーバ500は、疾患に対する治療薬を開発した製薬メーカや、治療装置を開発した医療装置メーカが管理することとしてもよい。
なお、サーバ200´、支援情報提供装置300a、分類器生成装置300b、コンピュータ400のハードウェアの構成は、基本的に、図5で説明した「データ処理部32」の構成と同様であるので、その説明は繰り返さない。
図39に戻って、相関行列算出部3002、相関値補正処理部3004、疾患識別器生成部3008、クラスタリング分類器生成部3010および判別値算出部3012、ならびに、機能結合の相関行列のデータ3106、測定バイアスデータ3108、補正後相関値データ3110および識別器データ3112については、実施の形態1で説明したのと同様であるので、その説明は、繰り返さない。
ここで、疾患識別器生成部3008、クラスタリング分類器生成部3010及び記憶装置210はクラスタリング装置を構成する。また、疾患識別器生成部3008、及びクラスタリング分類器生成部3010は、クラスタリング装置の演算装置2040bが実行する機能である。
MRI装置410は、クラスタリング分類器による層別化の対象となる被験者の脳活動データを計測し、コンピュータ400は、計測されたMRI測定データを不揮発性記憶装置4100に格納する。
コンピュータ400は、ユーザ(たとえば、医師)の送信の指示に従い、匿名化処理部4042で匿名化された被験者の測定データを脳活動の計測結果として支援情報提供装置300aに送信する。このとき、送信される測定データには、被験者を特定するための仮IDが付与されているものとする。好ましくは、仮IDと、たとえば、患者氏名などの個人情報との対応表は、コンピュータ400からはアクセス不可能な状態で、管理されているものとすることができる。また、少なくとも、支援情報提供装置300aにおいては、対応表の内容については一切アクセスすることができない構成とできる。
これに応じて、支援情報提供装置300aにおいては、必要に応じて、MRI装置410の設置されるサイトに対応したハーモナイゼーション処理を実行することが好ましい。さらに、判別値算出部3012は、入力された被験者の脳活動の計測結果について、クラスタリング分類器による層別化の結果である分類結果を出力する。治療情報出力部3200は、分類結果に応じて、分類結果および対応する治療法選択を支援する情報を、通信インタフェースを介して、コンピュータ400に送信する。
コンピュータ400の情報提示部4044は、支援情報提供装置300aの治療法情報生成部3200から返信されてくる特定の被験者についての治療法の選択を支援する情報および分類結果を、図示しないディスプレイなどの表示装置から医師に提示する。
図40は、治療法情報データベース5100の例を示す図である。
判別値算出部3012は、クラスタリング分類器を使って、被験者を、例えば、クラスタ1から5までに分類できるものとする。治療法情報データベース5100には、各クラスタに応じて、治療方法に関連する所定の治療法情報が格納されている。
治療法情報は、各クラスタに分類された被験者(特に、クラスタリング分類器の生成にそのデータを使用した患者群に属する被験者)の過去の治療履歴に関する情報や、文献等で報告されている効果及び/又は副作用に基づく参照すべき治療情報(参照治療情報)が格納されうる。好ましくは、各クラスタの特定の治療薬に対する応答性を示す情報、特定の物理的治療法に対する応答性を示す情報等が格納され得る。
図40では、治療履歴に関する情報において、推奨される治療方法が、効果の高かった順に第1候補から第3候補まで、クラスタごとに格納されている。また、参照治療情報として、副作用が出る恐れがある等、推奨されない治療方法も、クラスタごとに格納されている。
うつ病の治療薬として、特に制限されないが、アミトリプチリン塩酸塩、アモキサピン、イミプラミン塩酸塩等の三環系抗うつ薬;セチプチリンマレイン酸塩、マプロチリン塩酸塩、ミアンセリン塩酸塩等の四環系抗うつ薬;エスシタロプラムシュウ酸塩、セルトラリン塩酸塩、パロキセチン塩酸塩水和物、フルボキサミンマレイン酸塩等の選択的セロトニン再取り込み阻害薬;デュロキセチン塩酸塩、ベンラファキシン塩酸塩、ミルナシプラン塩酸塩等のセロトニン・ノルアドレナリン再取り込み阻害薬;ミルタザピン等のノルアドレナリン作動性・特異的セロトニン作動薬を挙げることができる。
うつ病の物理的治療法として、経頭蓋磁気刺激法、ニューロフィードバック法、電気けいれん療法や認知行動療法等を挙げることができる。
(治療法選択支援処理)
次に、支援情報提供装置300aの治療法情報生成部3200が行う治療法選択支援処理について説明する。治療法選択支援処理は、治療法選択支援プログラムをコンピュータで治療法情報生成部3200の処理として実行することにより、達成される。
図41は、被験者の治療法選択の支援処理のフローを説明するフローチャートである。
図41に示すステップS502において、支援情報提供装置300aは、コンピュータ400から被験者の脳活動の計測結果を受け付ける。被験者の脳活動の計測結果は、コンピュータ400において匿名化された撮像データ(脳構造画像データおよび脳機能画像データ)である。
図41に示すステップS504において、支援情報提供装置300aは、被験者の脳活動の計測結果に対して、機能結合の相関行列のデータの算出や、ハーモナイゼーション処理による相関値の補正の処理を実行した後、分類器生成装置300bが生成したクラスタリング分類器に入力し、被験者の層別化の結果である所属クラスタ確率を取得する。所属クラスタ確率は、クラスタリング分類器が層別化可能なすべてのクラスタについて、被験者が各クラスタに属する確率として出力される。支援情報提供装置300aは、確率が最も高いクラスタを被験者が属するクラスタである決定することができる。あるいは、確率が高い上位2つのクラスタを被験者が属する可能性があるクラスタであると決定してもよい。このとき、支援情報提供装置300aは、クラスタリング演算装置として機能する。
図41に示すステップS506において、支援情報提供装置300aは、ステップS504において取得した被験者の層別化の結果であるクラスタに基づいて、治療法情報データベース5100から、被験者のクラスタに対応する治療法情報を取得する。このとき、支援情報提供装置300aは、被験者のクラスタに対応する治療法情報について、少なくとも2つの治療法情報、すなわち第1候補、第2候補を取得してもよい。このようにすることで、治療の選択肢を広げることができる。
図41に示すステップS508において、支援情報提供装置300aは、ステップS506において、取得した治療法情報を、通信インタフェースを介して、コンピュータ400に出力する。
以上のような構成とすることで、支援情報提供装置300aは、被験者の脳活動の計測データに基づいて、医師に対して、この被験者に対する治療法の選択を支援する情報を提供することができる。
[スクリーニング支援システム、スクリーニング支援装置]
以下では、図39で説明した支援情報提供装置300aを、創薬における被験者のスクリーニングを実行するための支援装置として使用する態様について説明する。
なお、本明細書では、「治療法」とは、医師が「特定の薬剤を処方し、被験者にこの薬剤を所定の用量・用法で投与すること」や医師が「特定の治療プロセスを選択し、その治療法での治療を実施すること」を意味するものとし、「治療法候補」とは、特定の臨床試験等に基づいて、監督官庁からの承認または認証を受ける前の「治療法の候補」を意味するものとする。対象となる疾患が、精神疾患であれば、「治療プロセス」とは、たとえば、認知行動療法などを意味するものとする。
図42は、一般的な創薬のプロセスを示す図である。
図42に示すように、一般には、創薬の目的に従って、標的を探索し、候補物質についてのスクリーニング、最適化、有効性・安全性・薬物動態の検討を、動物実験や試験管内で細胞培養した細胞に対して創薬候補物質を投与し反応を測定する試験管内試験などで確認する(非臨床試験)とともに、工業化の検討を実施する。
以上の過程を経て、非臨床試験を通過できた候補物質について、いわゆる「治験」が実施される。
「治験」とは、医薬品の製造販売に関して、医薬品医療機器等法上の承認を得るために行われる臨床試験のことである。医薬品の場合は、多くの場合、治験は第I相から第III相までの3段階で行われる。第I相試験は、自由意思に基づき志願した健常成人を対象とし、第II相試験は第I相の結果をうけて、比較的軽度な少数例の患者を対象に、有効性・安全性・薬物動態などの検討を行い、第III相試験は、上市後に実際にその化合物を使用するであろう患者を対象に、有効性の検証や安全性の検討を主な目的として、より大きな規模で行われる。
精神神経系の創薬の場合、問題となるのは、現時点では、治療候補物質について動物モデルでヒトでの有効性を予測することが困難であることや、第I相試験を実施後に、最終的に上市にまで至れる確率が他の疾患の治療候補物質に比べて低いことがある。
そこで、このような問題の解決策の一つとしては、たとえば、第I相試験(または第II相試験)後に、有効性を発揮する可能性のある候補化合物を効率的に第III相へ進めるために、適切な被験者集団の同定を行うことがある。
以下では、支援情報提供装置300aを、このような被験者集団の同定(被験者スクリーニング)の支援装置として使用する場合について説明する。
すなわち、支援情報提供装置300aから、スクリーニング支援データとして、分類されたクラスタに関する情報をコンピュータ400に返信する構成とする。
なお、図40においても示したように、精神神経系の疾患に対する治療としては、治療薬によるものではなく、物理的療法、たとえば、rTMSによるものも知られている。このような物理的療法に使用する医療機器、医療機器プログラムについても、創薬候補物質と同様に、治験を経て、初めて、監督官庁から医療機器として承認される、という手続きが取られる。この場合の治験においても、後述するような被験者のスクリーニングが有効である場合が想定される。
そこで、本実施の形態で説明するようなスクリーニング支援処理は、上述したような「治療法候補」についての臨床試験において適用できるだけではなく、「治療手段候補」についての臨床試験に適用することも可能である。ここで、本明細書では、「治療手段」とは、医師が「特定の治療法を実施するために使用する機器」や「特定の治療法を実施するために使用するプログラム(または、そのプログラムを記録した媒体、あるいはそのプログラムがインストールされた機器)」を意味するものとし、「治療手段候補」とは、特定の臨床試験等に基づいて、監督官庁からの承認または認証を受ける前の「治療機器」や「プログラム(または、そのプログラムが記録された媒体やインストールされた機器)」を意味するものとする。たとえば、対象となる疾患が、精神疾患であれば、「治療手段(または治療機器)」とは、所定の用法で使用される、経頭蓋磁気刺激療法を実施するためのTMS装置、電気けいれん療法に使用されるパルス波治療器、認知行動療法を支援するスマートフォンのアプリケーションプログラム、またはこのようなアプリケーションプログラムが記録された記録媒体や、このようなアプリケーションプログラムがインストールされたスマートフォンなどを意味するものとする。
(スクリーニング支援システムとしての動作)
図44は、スクリーニング支援装置1000´の構成を説明する図である。
スクリーニング支援装置1000´は、治療法選択支援システム1000と同様に、医療機関4において、診断ラベルの推定対象となる被験者の脳活動データを計測するためのMRI装置410からのデータを受信可能なように接続されたコンピュータ400と、支援情報提供装置300aと、クラスタリング分類器の生成処理を実施する分類器生成装置300bと、データサーバ200´を備える。
スクリーニング支援装置1000´としての構成で、図39に示した治療法選択支援システム1000の構成と同一部分には同一符号を付して、その説明は繰り返さないこととし、以下では、スクリーニング支援装置1000´としての主要な動作について説明する。
すなわち、MRI装置410は、クラスタリング分類器による層別化の対象となる被験者の脳活動データを計測し、コンピュータ400は、計測されたMRI測定データを不揮発性記憶装置4100に格納する。
コンピュータ400は、ユーザ(たとえば、医師)の送信の指示に従い、匿名化処理部4042で匿名化された被験者の測定データを脳活動の計測結果として支援情報提供装置300aに送信する。このとき、送信される測定データには、被験者を特定するための仮IDが付与されているものとする。好ましくは、仮IDと、たとえば、患者氏名などの個人情報との対応表は、コンピュータ400からはアクセス不可能な状態で、管理されているものとすることができる。また、少なくとも、支援情報提供装置300aにおいては、対応表の内容については一切アクセスすることができない構成とできる。
これに応じて、支援情報提供装置300aにおいては、必要に応じて、MRI装置410の設置されるサイトに対応したハーモナイゼーション処理を実行することが好ましい。さらに、判別値算出部3012は、入力された被験者の脳活動の計測結果について、クラスタリング分類器による層別化の結果である分類結果を出力する。この分類結果は、被験者を特定する仮IDと関連づけて、たとえば、記憶装置2080´―1に格納される。スクリーニング情報出力部3202は、分類結果に応じて、分類結果に基づくスクリーニングを支援する情報を生成して、通信インタフェースを介して、コンピュータ400に送信する。
コンピュータ400の情報提示部4044は、支援情報提供装置300aのスクリーニング情報出力部3202から返信されてくる特定の被験者について、スクリーニングを支援する情報を、図示しないディスプレイなどの表示装置から医師に提示する。
なお、ここで、スクリーニング情報出力部3202からコンピュータ400に返信される情報を、「スクリーニングを支援する情報」としているのは、たとえば、スクリーニングの結果を用いて、臨床試験や治験を実施する場合に、いわゆる「二重盲検試験」や「ランダム化比較検試験」を実施するなど、医療機関側にも、層別化の結果そのものについては返信せず、これらの試験の特性に応じて、スクリーニングをいかに実施するかを支援する情報としてコンピュータに表示する場合が想定されるからである。
試験の終了後など、試験結果を評価する際には、仮IDに基づいて、試験結果と被験者についての分類結果との照合などが実施されるものとすることができる。
(支援情報提供装置300aにおけるスクリーニング支援処理)
次に、支援情報提供装置300aの演算装置2040aが行うスクリーニング支援処理について説明する。スクリーニング支援処理はスクリーニング支援プログラムをコンピュータで実行することにより、達成される。
図43は、創薬における被験者のスクリーニング支援処理のフローを説明するフローチャートである。
図43に示すステップS602において、支援情報提供装置300aは、コンピュータ400から被験者の脳活動の計測結果を受け付ける。被験者の脳活動の計測結果は、コンピュータ400において匿名化された撮像データ(脳構造画像データおよび脳機能画像データ)である。
図43に示すステップS604において、支援情報提供装置300aは、被験者の脳活動の計測結果に対して、機能結合の相関行列のデータの算出や、ハーモナイゼーション処理による相関値の補正の処理を実行した後、分類器生成装置300bが生成したクラスタリング分類器に入力し、被験者の層別化の結果である所属クラスタ確率を取得する。所属クラスタ確率は、クラスタリング分類器が層別化可能なすべてのクラスタについて、被験者が各クラスタに属する確率として出力される。支援情報提供装置300aは、確率が最も高いクラスタを被験者が属するクラスタである決定することができる。
図43に示すステップS606において、支援情報提供装置300aは、ステップS604において取得した被験者の層別化の結果であるクラスタを示す情報を生成する。
図43に示すステップS608において、支援情報提供装置300aは、ステップS606において生成したクラスタの情報またはスクリーニングを支援する情報を、通信インタフェースを介して、コンピュータ400に出力する。
以上のような構成とすることで、支援情報提供装置300aは、被験者の脳活動の計測データに基づいて、治験のプロセスにおいて、医師に対して、被験者の属するクラスタの情報を提供することができる。
医師は、予め定められた治験プロトコルに従って、対応するクラスタの被験者をスクリーニングして、治験を実施することが可能となる。
[うつ病患者に対する層別化の検証データ]
以上、図39においては、治療法選択支援システム1000について説明し、図42~図44においては、スクリーニング支援装置1000´について説明した。
以下では、以上説明した実施の形態において、判別値算出部3012が、入力された被験者の脳活動の計測結果について、クラスタリング分類器により算出する層別化の結果の臨床的な意義を検証するデータの一例について説明する。
(本検証データに使用したクラスタリング分類器)
本検証データに使用したクラスタリング分類器の生成にあたり疾患識別器生成部3008が生成するMDD識別器の構成は、以下のとおりである。
まず、脳領野のパーセレーション法としては、BSA法(Brainvisa Sulci Atlas)に従って脳領域を分割した。
たとえば、BSA法については、以下の文献に開示がある:
公知文献:Matthieu Perro, Denis Riviere, and Jean-Francois Mangin a, Cortical sulci recognition and spatial normalization. Medical Image Analysis Volume 15, Issue 4, August 2011, Pages 529-550
次に、脳領域間の結合指標としてfMRI信号のピアソン相関係数を算出し、全脳の脳機能結合を基にランダムフォレスト法を用いて、発見コホートデータにより、MDD識別器を学習して生成した。この場合において、トラベリングサブジェクト法による施設間補正は、特に適用していない。
そして、クラスタリング分類器生成部3010におけるクラスタリング分類器の生成処理については、以下のとおりである。
まず、MDD識別器を学習する過程では、入れ子構造の交差検証(Nested Cross Validation)を実施しており、結果として100個のMDD識別器が生成される。
一つ一つの識別器において、ランダムフォレスト法により、各脳機能結合の識別における重要度が決定されており、100個の識別器の総和を計算することで、総合的な重要度を求める。
この総合的重要度を基にして、脳機能結合の順位付けを行い、上位の結合を用いて多重共クラスタリングを実行する。ここで、特に限定されないが、上位何本の結合を用いるかは、所定のパターンの複数の本数についてクラスタリングを実施し、データセット間での安定性が最も高い本数を採用する構成とすることができる。
ここで、「データセット間での安定性」の評価には、図33で説明したARI値を評価指標として用いることができる。
図45は、データセット1についてMDD識別器における重要度上位30本の脳機能結合を用いた多重共クラスタリングによるクラスタリングの結果のビューを示す図である。
図46は、データセット2について多重共クラスタリングによるクラスタリングの結果のビューを示す図である。
図47は、データセット1とデータセット2との間のクラスタリング安定性を示す図である。
図47は、データセット1とデータセット2について、それぞれのビューに対するARIを算出した表であり、図33に示した表に対応するものである。
図47を参照すると、データセット1のビュー1とデータセット2のビュー1、データセット1のビュー2とデータセット2のビュー3は、それぞれ、ARIの値0.67と0.68となっており、互いに類似しているといえる。
図48は、全うつ病患者データ(データセット1+2)に対するクラスタリング結果のビューを示す図である。
図49は、データセット1+2のクラスタリングで生成されたビュー3について、全うつ病患者の数と、臨床データが存在するうつ病患者の数とを、サブタイプごとに示す図である。
ここで、「臨床データ」とは、各患者への「投薬の履歴の情報」、各患者の「うつ病の程度の診断情報」などを指す。
図49を参照すると、臨床データが存在する患者データについては、全体の一部(データセット1とデータセット2の双方の一部)であるのに対して、その一部の患者を取り出したとしても、全体をクラスタリングしたときと、その分布に大きな差がなく、したがって、「臨床データが存在する患者データ」だけについて、クラスタリングの臨床上の意義を検討しても問題はないことがわかる。
以下では、特徴的な臨床的意義、具体的にはサブタイプ間で治療反応性の差が認められたビューとして、図48,49に示したビュー3について、さらに説明する。
なお、各サブタイプが持つ脳機能結合の特徴(結合の平均と標準偏差)についても、臨床データを持つ一部の患者を取り出しても、全体の平均や標準偏差について、大きく変わらないことも確認された。
図50は、クラスタリングに使用された脳機能結合 (ビュー3)を示す図である。
図50(a)は、クラスタリングに使用された脳機能結合の脳内の位置を示し、図50(b)は、機能結合の関心領域を示す。
ビュー3においては、図50(a)(b)に示された3つの結合のみが、クラスタリングに使用されている。
クラスタリング分類器を生成するために使用したMDD識別器における重要度について、上位30本の結合を選択しておき、これらの結合を使って多重共クラスタリングを実施している。
つまり、30本の選択された結合に重みなどは与えずクラスタリングを実施する。そしてその結果、ビューが3つ生成されたことになる。言い換えると、30本の結合が3つのグループに分けられ(ビュー1:22本、ビュー2:5本、ビュー3:3本)、それぞれの結合グループを基に、被検者がクラスタリングされた。
多重共クラスタリングでは、結合をどのように分けて患者(被験者)をどのようにクラスタリングするかはアルゴリズムが自動的に最適解を導くことになる。
なお、図50(b)における関心領域の名称は、以下の内容をそれぞれ示す。
Thalamus_LorR:左or右視床
Precentral_R:右中心前回
Postcentral_L:左中心後回
図51は、図48、図49に示したビュー3の各サブタイプについて、うつ病の重症度及び重症度の改善率との関係を示す図である。
図51(a)は、治療開始後、第0週と第6週のHAMDのスコアを対比して示す図である。
ここで、HAMDとは、ハミルトンうつ病評価尺度のことをいい、うつ病の重症度を評価するための尺度である。うつ病の重症度をあらわす17項目で構成された主要17項目版とこれに追加の4項目を加えた21項目版が主に用いられている。
また、図51(a)において、図49においてサブタイプ1~5とは、図49の「被検者クラスタリングの結果 (ビュー3)」の上位5つのクラスタである。
なお、サブタイプ番号は臨床データが存在する人数の多い順になっている。また、10人以上の被験者データが存在するサブタイプを解析対象としたために、サブタイプ1~5のうち、サブタイプ1,2,4,5がグラフに示されている。
治療開始後0週とは、研究エントリー時及びSSRIによる治療開始を意味し、治療開始後6週とは、研究エントリー及び治療開始から6週間後を意味する。すなわち、研究へのエントリー前からSSRIによる治療開始した者は除かれている。
ここで、SSRIとは、「選択的セロトニン再取り込み阻害薬(Selective Serotonin Reuptake Inhibitor)」のことをいい、たとえば、エスシタロプラムなどが含まれる。
ここでは、研究における治療プロトコルとしては、以下のようなものである。
未治療、もしくは不十分な量・期間の治療しか受けていない大うつ病患者を対象とし、研究にエントリーし、研究エントリー時よりエスシタロプラムなどのSSRIで治療開始、臨床的な判断により増量可能、併用薬の制限はなし、SSRIで良くならない場合は別の薬物療法を実施する。
図51(b)は、第0週と第6週を比較した場合のHAMDの改善率を示す図である。
ここで、改善率とは、以下の式で表される。
(HAMD(0) - HAMD(6W)) / HAMD(0) x 100
図51(a)より、まず初期の段階ではサブタイプ間でHAMDに差がない、つまりうつ重症度は同程度であることが分かる。
図51(a)より、SSRIによる治療開始から6週間後においてはサブタイプ間にHAMDの差が認められることがわかる。図51(b)からは、初期値からの改善率(improvement rate)を見ると分かる通り、サブタイプ1はサブタイプ2や5と比較してSSRIの治療効果が低く、逆にサブタイプ2や5はサブタイプ1に比べてSSRIの治療効果が高いということが示される。
以上の説明より、実施の形態で説明したクラスタリング分類器を用いることで、治療初期の患者に対して、SSRIという薬剤の治療応答性を、クラスタごとに予測することが可能であることがわかる。
今回開示された実施の形態は、本発明を具体的に実施するための構成の例示であって、本発明の技術的範囲を制限するものではない。本発明の技術的範囲は、実施の形態の説明ではなく、特許請求の範囲によって示されるものであり、特許請求の範囲の文言上の範囲および均等の意味の範囲内での変更が含まれることが意図される。
2 被験者、6 ディスプレイ、10 MRI装置、11 磁場印加機構、12 静磁場発生コイル、14 傾斜磁場発生コイル、16 RF照射部、18 寝台、20 受信コイル、21 駆動部、22 静磁場電源、24 傾斜磁場電源、26 信号送信部、28 信号受信部、30寝台駆動部、32 データ処理部、36 記憶部、38 表示部、40 入力部、42 制御部、44 インタフェース部、46 データ収集部、48 画像処理部、50 ネットワークインタフェース、300a 支援情報提供装置、300b 分類器生成装置、 500 治療法情報提供サーバ。

Claims (28)

  1. うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援システムであって、
    複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング装置を備え、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
    前記クラスタリング装置は、
    前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含み、
    前記演算装置は、クラスタリング分類器の生成処理において、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、
    ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成し、
    前記治療法選択支援システムは、さらに、
    前記クラスタリング分類器による層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置と、
    前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果に応じて、対応する治療法情報を出力する支援情報提供装置とを、備える、治療法選択支援システム。
  2. 前記演算装置は、前記識別器モデルを生成する機械学習において、
    前記第1の群と前記第2の群から、アンダーサンプリングおよびサブサンプリングを実行して、複数の学習用サブサンプルを生成し、
    前記学習用サブサンプルのそれぞれについて、機械学習による識別器の生成において使用される特徴量の和集合から、前記和集合に属する特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
    選択された前記クラスタリングのための特徴量に基づいて、前記多重共クラスタリング法により、前記クラスタリング分類器を生成する、請求項1記載の治療法選択支援システム。
  3. 前記支援情報提供装置は、クラスタリング演算装置とインタフェース装置とを備え、
    前記クラスタリング演算装置は、前記クラスタリング分類器により前記第1の被験者が各前記クラスタに属する確率を算出し、前記データベース装置から、前記確率に応じて選択される少なくとも2つの前記治療法情報を読み出し、
    前記インタフェース装置は、前記選択されたクラスタとそれぞれ対応する前記治療法情報とを関連付けて表示するためのデータを出力する、請求項1または2に記載の治療法選択支援システム。
  4. 前記治療法情報は、特定の治療薬に対する応答性を示す情報である、請求項1~3のいずれか1項に記載の治療法選択支援システム。
  5. 前記治療法情報は、特定の物理的治療法に対する応答性を示す情報である、請求項1~3のいずれか1項に記載の治療法選択支援システム。
  6. 前記機械学習による識別器の生成の処理は、前記複数の学習用サブサンプルに対してそれぞれ複数の識別器サブモデルを生成し、前記複数の識別器サブモデルを統合して前記識別器モデルを生成するアンサンブル学習である、請求項2に記載の治療法選択支援システム。
  7. 前記クラスタリング装置は、複数の計測サイトにそれぞれ設けられた複数の脳活動計測装置から、各前記複数の第2の被験者の所定の複数の脳領野ペア間の脳活動の時間相関を表現する情報を受け、
    前記演算装置は、各前記複数の第2の被験者について前記複数の脳機能結合相関値について、前記計測サイトの測定バイアスを除去するように補正することで、補正された調整値を前記特徴量として前記記憶装置に格納するハーモナイゼーション算出手段を含む、請求項1または2に記載の治療法選択支援システム。
  8. うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援装置であって、
    複数の第2の被験者のうち、うつ病の診断ラベルを有する被験者に対する層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置と、
    前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に基づく層別化の結果に応じて、対応する治療法情報を出力する支援情報提供装置とを備え、
    前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
    前記層別化の結果のクラスタは、クラスタリング装置による脳機能結合相関値の計測結果に対するクラスタリング処理により得られたクラスタリング分類器により得られるものであり、
    前記クラスタリング装置は、
    前記第1の群について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含み、
    前記クラスタリング分類器の生成処理において、前記演算装置が、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、
    ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングすることにより、前記クラスタリング分類器を生成する、治療法選択支援装置。
  9. 前記支援情報提供装置は、クラスタリング演算装置とインタフェース装置とを備え、
    前記クラスタリング演算装置は、前記クラスタリング分類器により前記第1の被験者が各前記クラスタに属する確率を算出し、前記データベース装置から、前記確率に応じて選択される少なくとも2つの前記治療法情報を読み出し、
    前記インタフェース装置は、前記選択されたクラスタとそれぞれ対応する前記治療法情報とを関連付けて表示するためのデータを出力する、請求項8に記載の治療法選択支援装置。
  10. 前記治療法情報は、特定の治療薬に対する応答性を示す情報である、請求項8または9に記載の治療法選択支援装置。
  11. 前記治療法情報は、特定の物理的治療法に対する応答性を示す情報である、請求項8~10のいずれか1項に記載の治療法選択支援装置。
  12. うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援方法であって、
    複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング分類器を生成して準備するステップを備え、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
    前記準備するステップは、
    前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算ステップを含み、
    前記演算ステップは、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、
    ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成するステップとを有し、
    前記治療法選択支援方法は、さらに、
    前記第1の被験者の脳活動の計測結果に対する前記クラスタリング分類器による分類結果に応じて、前記クラスタリング分類器による層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベースから、対応する治療法情報を取得し出力する支援情報提供ステップを、備える、治療法選択支援方法。
  13. うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援方法であって、
    前記第1の被験者の脳活動の計測結果に基づく層別化の結果のクラスタに応じて、複数の第2の被験者のうち、うつ病の診断ラベルを有する被験者に対する層別化の結果と対応する所定の治療法情報とを関連付けて格納するためのデータベースから、対応する治療法情報を取得し出力する支援情報提供ステップを備え、
    前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
    前記層別化の結果のクラスタは、脳機能結合相関値の計測結果に対するクラスタリング処理により得られたクラスタリング分類器により得られるものであり、
    前記クラスタリング分類器は、
    前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算ステップにより生成されたものであり、前記演算ステップは、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、
    ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成するステップとを含む、治療法選択支援方法。
  14. うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援プログラムであって、
    前記治療法選択支援プログラムは、コンピュータに実行させたときに、コンピュータに、
    複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング分類器を生成するステップと、
    前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果に応じて、前記クラスタリング分類器による層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置から、対応する治療法情報を取得し出力するステップとを実行させ、
    前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
    前記クラスタリング処理は、
    前記複数の第2の被験者について、クラスタリングの処理を実行するための演算ステップを含み、前記演算ステップは、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を記憶装置に格納するステップと、
    ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成するステップとを含む、治療法選択支援プログラム。
  15. うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援プログラムであって、
    前記治療法選択支援プログラムは、コンピュータに実行させたときに、
    コンピュータに、前記第1の被験者の脳活動の計測結果に基づく層別化の結果のクラスタに応じて、複数の第2の被験者のうち、うつ病の診断ラベルを有する被験者に対する層別化の結果と対応する所定の治療法情報とを関連付けて格納するためのデータベースから、対応する治療法情報を取得し出力する支援情報提供ステップを実行させ、
    前記層別化の結果のクラスタは、脳機能結合相関値の計測結果に対するクラスタリング処理により得られたクラスタリング分類器により得られるものであり、
    前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
    前記クラスタリング分類器は、
    前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算ステップにより生成されたものであり、前記演算ステップは、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、
    ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成するステップとを含む、治療法選択支援プログラム。
  16. うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援システムであって、
    複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング装置を備え、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
    前記クラスタリング装置は、
    前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含み、前記演算装置は、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、
    ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成し、
    前記スクリーニング支援システムは、さらに、
    前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力する支援情報提供装置とを、備える、スクリーニング支援システム。
  17. 前記演算装置は、前記識別器モデルを生成する機械学習において、
    前記第1の群と前記第2の群から、アンダーサンプリングおよびサブサンプリングを実行して、複数の学習用サブサンプルを生成し、
    前記学習用サブサンプルのそれぞれについて、機械学習による識別器の生成において使用される特徴量の和集合から、前記和集合に属する特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
    選択された前記クラスタリングのための特徴量に基づいて、前記多重共クラスタリング法により、前記クラスタリング分類器を生成する、請求項16記載のスクリーニング支援システム。
  18. うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援装置であって、
    クラスタリング分類器を特定する情報を格納するための記憶装置を有する支援情報提供装置であって、前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器に基づく分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力する支援情報提供装置を備え、
    別化の結果のクラスタは、クラスタリング装置による、複数の第2の被験者から取得した脳機能結合相関値の計測結果に対するクラスタリング処理により得られた前記クラスタリング分類器により得られるものであり、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
    前記クラスタリング装置は、
    前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含み、
    前記クラスタリング分類器の生成処理において、前記演算装置が、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、
    ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成する、スクリーニング支援装置。
  19. うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援方法であって、
    記憶装置に格納された情報で特定されるクラスタリング分類器に基づいて、前記脳活動の計測結果により、演算装置が前記第1の被験者の分類を実行するステップと、
    前記分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力するステップと、を備え、
    前記クラスタリング分類器の生成のための処理は、
    うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む複数の第2の被験者について、クラスタリング処理を実行するための演算ステップを含み、前記演算ステップは、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、
    ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成するステップとを含む、スクリーニング支援方法。
  20. うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援プログラムであって、
    前記スクリーニング支援プログラムは、コンピュータに実行させたときに、コンピュータに、
    記憶装置に格納された情報で特定されるクラスタリング分類器に基づいて、前記脳活動の計測結果により、演算装置が前記第1の被験者の分類を実行するステップと、
    前記分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力するステップと、
    を実行させ、
    前記クラスタリング分類器の生成のための処理は、
    うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む複数の第2の被験者について、クラスタリング処理を実行するための演算ステップを含み、前記演算ステップにおいて、
    i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を記憶装置に格納するステップと、
    ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、
    iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、
    iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成するステップを含む、スクリーニング支援プログラム。
  21. 前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である、請求項1に記載の治療法選択支援システム。
  22. 前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である、請求項8に記載の治療法選択支援装置。
  23. 前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である、請求項12または13に記載の治療法選択支援方法。
  24. 前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である、請求項14または15に記載の治療法選択支援プログラム
  25. 前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である、請求項16に記載のスクリーニング支援システム
  26. 前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である、請求項18に記載のスクリーニング支援装置。
  27. 前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である、請求項19に記載のスクリーニング支援方法。
  28. 前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である、請求項20に記載のスクリーニング支援プログラム。
JP2022536448A 2020-07-17 2021-07-15 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム Active JP7737099B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2020123208 2020-07-17
JP2020123208 2020-07-17
PCT/JP2021/026658 WO2022014682A1 (ja) 2020-07-17 2021-07-15 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム

Publications (2)

Publication Number Publication Date
JPWO2022014682A1 JPWO2022014682A1 (ja) 2022-01-20
JP7737099B2 true JP7737099B2 (ja) 2025-09-10

Family

ID=79555705

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022536448A Active JP7737099B2 (ja) 2020-07-17 2021-07-15 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム

Country Status (5)

Country Link
US (1) US20230284983A1 (ja)
EP (1) EP4183333A4 (ja)
JP (1) JP7737099B2 (ja)
CN (1) CN116133587A (ja)
WO (1) WO2022014682A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4308017A4 (en) * 2021-03-18 2025-02-12 The Board of Trustees of the Leland Stanford Junior University Systems and methods to generate a surgical risk score and uses thereof
US20250342960A1 (en) * 2022-05-19 2025-11-06 Shionogi & Co., Ltd. Information processing system and information processing method
US12586360B1 (en) * 2022-12-19 2026-03-24 Whiterabbit.Ai Inc. Methods and systems for monitoring and surveillance of autonomous AI in medical imaging
US20240304330A1 (en) * 2023-03-06 2024-09-12 GE Precision Healthcare LLC System and method for training artificial intelligence models using sub-group training datasets of a majortiy class of samples
CN117224105B (zh) * 2023-10-31 2024-08-23 重庆博恩富克医疗设备有限公司 基于多模式电磁场生物检测的颅脑异常预测方法及装置
CN119557677B (zh) * 2024-11-11 2025-10-24 杭州电子科技大学 一种基于脑电协同聚类的半监督言语想象意图解码方法
CN120260957B (zh) * 2025-03-26 2025-11-18 浙江大学 一种数据处理方法及系统
CN121421537B (zh) * 2025-12-17 2026-03-31 山东第一医科大学(山东省医学科学院) 基于静息态功能磁共振成像的情绪改善评估方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012083136A1 (en) 2010-12-17 2012-06-21 The Trustees Of Columbia University In The City Of New York Apparatus, method and computer-accessible medium for diagnosing and subtyping psychiatric diseases
WO2014178323A1 (ja) 2013-05-01 2014-11-06 株式会社国際電気通信基礎技術研究所 脳活動解析装置、脳活動解析方法およびバイオマーカー装置
US20190090749A1 (en) 2017-09-26 2019-03-28 Washington University Supervised classifier for optimizing target for neuromodulation, implant localization, and ablation
JP2019063478A (ja) 2017-10-03 2019-04-25 株式会社国際電気通信基礎技術研究所 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2245568A4 (en) * 2008-02-20 2012-12-05 Univ Mcmaster EXPERIMENTAL SYSTEM FOR DETERMINING THE RESPONSE OF PATIENTS TO TREATMENT
US9632162B2 (en) 2013-12-06 2017-04-25 Toshiba Medical Systems Corporation Method of, and apparatus for, correcting distortion in medical images
JP6195329B1 (ja) * 2015-11-24 2017-09-13 株式会社国際電気通信基礎技術研究所 脳活動解析装置、脳活動解析方法、プログラムおよびバイオマーカー装置
JP7057995B2 (ja) * 2017-10-03 2022-04-21 株式会社国際電気通信基礎技術研究所 哺乳動物におけるうつ病の病態プロファイルを検出するバイオマーカーとその利用
WO2019182915A1 (en) * 2018-03-19 2019-09-26 The Board Of Trustees Of The Leland Stanford Junior University Treatment of depression
CN108596236A (zh) * 2018-04-18 2018-09-28 东南大学 一种基于全局连接特征和测地距离的丘脑核团分区方法
JP6812022B2 (ja) * 2018-10-11 2021-01-13 株式会社国際電気通信基礎技術研究所 脳機能結合相関値の調整方法、脳機能結合相関値の調整システム、脳活動分類器のハーモナイズ方法、脳活動分類器のハーモナイズシステム、および脳活動バイオマーカシステム
US12383196B2 (en) * 2018-10-15 2025-08-12 The Board Of Trustees Of The Leland Stanford Junior University Treatment of depression using machine learning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012083136A1 (en) 2010-12-17 2012-06-21 The Trustees Of Columbia University In The City Of New York Apparatus, method and computer-accessible medium for diagnosing and subtyping psychiatric diseases
WO2014178323A1 (ja) 2013-05-01 2014-11-06 株式会社国際電気通信基礎技術研究所 脳活動解析装置、脳活動解析方法およびバイオマーカー装置
US20190090749A1 (en) 2017-09-26 2019-03-28 Washington University Supervised classifier for optimizing target for neuromodulation, implant localization, and ablation
JP2019063478A (ja) 2017-10-03 2019-04-25 株式会社国際電気通信基礎技術研究所 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置

Also Published As

Publication number Publication date
EP4183333A4 (en) 2023-12-20
EP4183333A1 (en) 2023-05-24
JPWO2022014682A1 (ja) 2022-01-20
WO2022014682A1 (ja) 2022-01-20
US20230284983A1 (en) 2023-09-14
CN116133587A (zh) 2023-05-16

Similar Documents

Publication Publication Date Title
JP7365496B2 (ja) 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラム、脳活動マーカー分類システムおよび脳機能結合相関値のクラスタリング分類器モデル
JP7737099B2 (ja) 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム
JP7247159B2 (ja) 脳活動分類器のハーモナイズシステム、及び脳活動分類器プログラム
US12390108B2 (en) Brain functional connectivity correlation value adjustment method, brain functional connectivity correlation value adjustment system, brain activity classifier harmonization method, brain activity classifier harmonization system, and brain activity biomarker system
CN111447874B (zh) 判别装置、抑郁症状的判别方法、抑郁症状的水平的判断方法、抑郁症患者的分类方法、抑郁症状的治疗效果的判断方法和脑活动训练装置
JP6497687B2 (ja) 脳活動解析システム及び方法
Schnyer et al. Evaluating the diagnostic utility of applying a machine learning algorithm to diffusion tensor MRI measures in individuals with major depressive disorder
JP6195329B1 (ja) 脳活動解析装置、脳活動解析方法、プログラムおよびバイオマーカー装置
CN105188528B (zh) 脑活动训练装置以及脑活动训练系统
Li et al. Sparse multivariate autoregressive modeling for mild cognitive impairment classification
JP2024003782A (ja) 統合失調スペクトラム障害に関する情報を提供するための診断支援システム、診断支援装置、診断支援方法、及び支援プログラム
Leming Application of deep learning to brain connectivity classification in large mri datasets
Ost Predictive Modeling of Community Acquired MRI After a Diagnosis of Degenerative Cervical Myelopathy: Baseline Disease Severity and Surgical Outcomes
Hansen Irritable Bowel Syndrome and Neuroimaging-based Biomarkers
Chen Machine Learning for the Classification of Dementia in the Presence of Mis-labelled Data
Johnston DOCTOR OF MEDICINE

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240626

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20250311

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20250512

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20250722

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250820

R150 Certificate of patent or registration of utility model

Ref document number: 7737099

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150