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
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/10—ICT 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/70—ICT 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features 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/004—Features 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/0042—Features 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
- A61B5/165—Evaluating the state of mind, e.g. depression, anxiety
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/70—ICT 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)。
以下では、人工知能技術による判別やクラスタリングを適用する分野として、医学分野を例にとる。
生体内の生物学的変化を定量的に把握するため、生体情報を数値化・定量化した指標のことを「バイオマーカー」と呼ぶ。
一方で、いわゆる画像診断の分野においては、上述したような「生物学的な分子」というバイオマーカーの概念とは異なり、「画像バイオマーカー」と呼ばれるものも存在する。たとえば、脳神経領域での分子イメージングに、PET(positron emission tomography)を用いて、神経伝達機能や受容体機能解析を行おうとする試みもある。
fMRIでは、装置としては通常のMRI装置に、さらに、fMRI計測に必要なハードおよびソフトを装備したものが使用される。
特に、人の運動機能に関する研究では、被験者に何らかの運動を行わせつつ、上記fMRI計測によって脳の活動を測定することが行われている。
さらに、このような「診断バイオマーカー」について、汎化性能を一層向上させるための技術改良もなされている(特許文献7)。
したがって、たとえば、うつ病に関しては、独立した検証データにおける層別化の精度確認は未実施というのが現状である。
上述したような計測サイト間のサイト間差を補償するための「ハーモナイゼーション法」については、本件の出願人による報告もある(非特許文献8)。
好ましくは、前記治療法情報は、特定の治療薬に対する応答性を示す情報である。
好ましくは、前記治療法情報は、特定の物理的治療法に対する応答性を示す情報である。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
好ましくは、前記治療法情報は、特定の物理的治療法に対する応答性を示す情報である。
好ましくは、前記機械学習による識別器の生成の処理は、前記複数の学習用サブサンプルに対してそれぞれ複数の識別器サブモデルを生成し、前記複数の識別器サブモデルを統合して前記識別器モデルを生成するアンサンブル学習である。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
好ましくは、前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である。
好ましくは、前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である。
好ましくは、前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である。
好ましくは、前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である。
好ましくは、前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である。
以下、本発明の「脳機能結合相関値のクラスタリング装置」、「脳機能結合相関値のクラスタリング方法」等を説明するために、複数の脳活動計測装置からなる計測システムで計測した被験者(精神疾患の患者を含む)の脳機能結合画像データに対する人工知能技術による「クラスタリング」を例として説明することとする。
図1は、複数の計測サイトに設置されたMRI計測システムにより計測されたデータについて、クラスタリング(層別化)処理を説明するための概念図である。
図1を参照して、計測サイトMS.1~MS.Ns(Ns:サイト数)には、それぞれ、MRI装置100.1~100.Nsが設置されているものとする。
そして、各計測サイトにおいては、原則として、MRI装置の仕様上、可能な範囲で、統一した計測プロトコルで各被験者の計測が実施されるものとする。
ここで、特に限定されないが、計測プロトコルとしては、たとえば、以下のような内容が規定されているものとする。
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号明細書
なお、計測プロトコルには、上記の条件の中から適宜、必要となるシーケンス部分を抜き出してもよいし、必要に応じて、他のシーケンスやその条件を追加してもよい。
たとえば、上述した例では、従来の診断基準で、「大うつ病」と診断される患者には、実は、いくつかのサブタイプが含まれることが知られている。
また、MRI装置100.1~100.Nsも、各計測サイトで決して同一計測特性を有するMRI装置が使用されているとは言えない。
以下では、このような「計測関連データ」について説明する。
すなわち、図1における計算処理システム300が、記憶装置210内に格納されたデータに基づいて、後述するような脳活動バイオマーカーを算出する際の基礎となるデータとして、これらのデータが使用される。
なお、特に限定されないが、EPI歪みを画像処理で補正する処理については、各計測サイトで演算処理がされた後に、記憶装置210内にデータが格納される構成とすることができる。
機能的なコネクティビティは、各参加者に対して2つの脳領域間の安静時の機能MRIの血中酸素濃度に依存する(BOLD)信号の時間的な相関性として計算される。
ここで、関心領域としては、上述したようにNr領域を考えているので、相関行列における独立な非対角成分は、対称性を考慮すると、
Nr×(Nr-1)/2(個)
ということになる。
関心領域の設定の方法としては、以下のような方法が想定される。
ここでは、脳活動バイオマーカーのために、たとえば、関心領域として140の領域を採用する。
公知文献2:Perrot et al., Med Image Anal, 15(4), 2011
公知文献3:Tzourio-Mazoyer et al., Neuroimage, 15(1), 2002
このような関心領域としては、たとえば、以下のような領域である。
背内側前頭前皮質(DMPFC)
前頭前野腹内側部(VMPFC)
前帯状皮質(ACC)
小脳虫部、
左視床、
右下頭頂葉、
右尾状核、
右中後頭葉、
右中帯状皮質
ただし、採用する脳の領野は、このような領域に限定されるものではない。
たとえば、対象とする神経・精神疾患に応じて、選択する領域を変更してもよい。
ここで、このような機能的な脳地図の脳領域については、以下のような文献にも開示がされており、特に限定されないが、たとえば、268のノード(脳領域)から成る、という構成とすることができる。
脳領域のパーセレーションについては、ヒューマンコネクトームプロジェクト(HCP)スタイルのマルチモダリティなイメージング(ミエリン・タスク・ファンクショナル)を用いることにより、脳を脳溝に沿ってシート状に変換して作成された脳地図に基づく「表面ベースの方法」で、データを分析することも可能である。
https://edickie.github.io/ciftify/#/
したがって、BOLD信号の時間的変化は、これらの379個の関心領域(ROI)から抽出される。
下記の文献にも開示があるように、事前情報(脳地図)なしに位相のそろったボクセルからネットワークを新規に同定する手法であり、「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.
しかしながら、以下では、特に限定されないが、ピアソンの相関係数を使用するものとして説明する。
前処理されたBOLD信号の可能なノード組の各々の時間経過の間についてのフィッシャーのz変換後のピアソンの相関係数が計算され、要素がそれぞれ2つのノード間のコネクティビティの強度を表わす、379×379の対照的な脳機能結合マトリックスを構築するために使用される。
図3は、「測定パラメータ」および「被験者属性データ」の内容の例を示す概念図である。
「測定パラメータ」は、「計測装置に関する情報」と「計測条件に関する情報」とを含む。
「計測装置に関する情報」には、各計測サイトにおいて被験者の脳活動を計測するためのMRI装置のメーカー名と、型番と、(送)受信コイル数とを含む。
「計測条件に関する情報」も、これらのものに限られるものではない。
なお、「被験者の属性情報」は、必要に応じて、たとえば、計測サイトにおいて、匿名加工処理が行われているものとする。
また、投薬履歴や診断履歴は、日付のランダム化やシフト化(相対化)など、必要に応じて、匿名化のための処理を行う。
図4は、各計測サイトに設置されるMRI装置100.i(1≦i≦Ns)の全体構成を示す模式図である。
図4においては、1番目の計測サイトのMRI装置100.1を、例示的に、詳しく説明している。他のMRI装置100.2~100.Nsについても、基本的な構成は同様である。
なお、ここで、被験者2が載置される円筒形状のボアの中心軸をZ軸にとりZ軸と直交する水平方向にX軸及び鉛直方向にY軸を定義する。
この出力されたNMR信号を被験者2からの応答波として受信コイル20で受信し、データ処理部32において、被験者2の関心領域が画像化される。
磁場印加機構11は、静磁場発生コイル12と、傾斜磁場発生コイル14と、RF照射部16と、被験者2をボア中に載置する寝台18とを備える。
このような視覚刺激は、被験者にニューロフィードバックを行う場合は、フィードバック情報の提示に相当する。
傾斜磁場発生コイル14は、Xコイル、Yコイル及びZコイル(図示省略)から構成され、円筒形状を示す静磁場発生コイル12の内周面に設けられる。
なお、傾斜磁場の均一性を向上させるためにシムコイル(図示せず)が設けられ、「シム調整」が実施される。
(送)受信コイル20については、上述のとおり、SN比の向上のためにマルチアレイコイルが使用される。
図5は、データ処理部32のハードウェアブロック図である。
データ処理部32のハードウェアとしては、上述のとおり、特に限定されないが、汎用コンピュータを使用することが可能である。
(脳機能結合に基づく、疾患/健常の判別器の生成処理、およびクラスタリング処理)
図6は、図2で説明したような相関行列から、診断マーカーとなる判別器を生成する過程およびクラスタリング処理について説明する概念図である。
以下、図6中の各処理について、さらに詳しく説明する。
[前処理から疾患識別器の生成およびクラスタリング処理までの概要]
(安静時機能的結合FCマトリックスの前処理および計算)
計測されたfMRIデータの、たとえば最初の10秒間はT1平衡を考慮に入れるために廃棄される。
このような前処理のパイプライン処理については、たとえば、以下のサイトに開示がある。
http://fmriprep.readthedocs.io/en/latest/workflows.html
脳領域のパーセレーションについては、特に限定されないが、上述した「方法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つのボリュームの間の頭部運動を表わす。
本実施の形態の具体例においては、機能的結合FCは、上述したようなパーセレーション法により領域の分割がされた後に、各参加者につき379個の関心領域(ROI)にわたってBOLD信号の時間的な相関として計算される。
機能的結合の計算では、特に限定されないが、ここでは、ピアソンの相関係数を使用する。
さらに、分析のために、結合マトリックスの下三角行列の71,631(=(379×378)/2)の機能的結合FCの値が使用される。
(脳活動バイオマーカーのためのハーモナイゼーション処理)
適切にそのような異種条件下のデータを管理するために、サイト間のデータをハーモナイズ(調和)させることが必要である。
サイト間差は、本質的に2つのタイプのバイアスを含んでいる。
技術的なバイアス(つまり測定バイアス)および生物学的なバイアス(つまり標本バイアス)である。
(アンサンブル学習による疾患識別器)
(クラスタリングのための特徴量選択とクラスタリング)
特に、限定されないが、たとえば、以下のようにして「重要度」が決定される。
以下では、図6で説明した処理のうち、アンサンブル学習による識別器生成について、より詳しく説明する。
すなわち、2クラス分類のための分類器生成処理、より特定的には、疾患識別器(「健常」または「疾患」についての2クラス分類器)の訓練データとして、学習用データセットを使用して、MDDのためのバイオマーカーを構築する処理を例として説明する。
そこで、以下では、機能的結合FCに基づいて、健常者(HC)およびMDD患者を識別するMDD識別器を構築する手続きについて説明する。
図8は、図7の構成をより詳しく説明する機能ブロック図である。
また、図9は、アンサンブル学習による疾患識別器を生成するための機械学習の手続きを説明するためのフローチャートである。
続いて、ハーモナイゼーション算出部3020が、算出した測定バイアスを使用して、ハーモナイゼーション処理を実施する(S104)。
後述するとおり、ハーモナイゼーション処理は、トラベリングサブジェクトを用いた方法が望ましいが、他の方法であってもよい。
さらに、「サブサンプリング処理」とは、トレーニングデータセットから所定数のサンプルをランダムに抽出する処理を意味する。
ウェイトベクトルwは、以下の評価関数(コスト関数)を最小化するように決定される(LASSO計算)。
ここで、λはハイパーパラメータを表わし、評価に適用される収縮量をコントロールする。
MDD識別器の出力(診断の確率値)が、0.5を超えるときは、MDD患者を示す指標と見なすことができる。
既述したように、機械学習アルゴリズムを使用する信頼できる分類器および回帰モデルの構築については、多数の撮像サイトから集められた大規模なサンプルサイズのデータを用いることが必要である。
図10は、このような学習用のデータセット(データセット1)のデモグラフィック特性を示す図である。
データセット1は、上述したSRPBS中のデータである。
図11は、独立検証のデータセット(データセット2)のデモグラフィック特性を示す図である。
データセット2も、基本的には、上述したSRPBS中のデータである。
すなわち、以下の分析では、2つの安静時機能的MRI(rs-fMRI)データセットを使用する。
データセット1は、「学習用データセット」であり、MDDの識別器およびクラスタリングの分類器を構築するのに使用される。
参加者は、それぞれ10分間の単一の安静時機能的MRI(rs-fMRI)セッションにて計測が実行されている。
安静時機能的MRI(rs-fMRI)のスキャン中に、参加者は、原則として、以下のような教示がされている。
「リラックスしてください。眠らないでください。中央の十字線マーク上を注視して、特定のことに関して考えないでください。」
なお、図10および図11において、括弧中の人数は、BDIスコアのデータを有する参加者の数を示す。
デモグラフィック分布は、すべて学習用データセットの中でMDDとHCの個体群間で統計的に有意な差はない(p>0.05)。
データセット2は、 「独立した検証データセット」であり、MDDの分類器およびクラスタリングの分類器をテストするために使用される。
データセット2に対する撮像を行ったサイトは、データセット1に含まれていない。
また、以下では、機能的結合FCの上のサイト効果を制御するために、後述するような、学習用データセットに対するトラベリングサブジェクトのハーモナイゼーション方法を使用するものとして説明する。
ただし、ハーモナイゼーション方法としては、この方法に限定されるものではなく、たとえば、これもComBat法など他の方法を使用してもよい。
なお、ComBat法については、たとえば、以下の文献に開示がある。
トラベリングサブジェクトのハーモナイゼーション方法を使用することで、純粋なサイト間差(測定バイアス)を除去することが可能となる。
図12は、全撮像サイトについて、学習用データセットに対するMDDの予測性能(出力の確率分布)を示す図である。
識別器モデルは、MDD患者を66%の正確度でHC個体群から分離する。
対応するAUCは、0.77であり、高い識別力を示した。
また、MCCは、約0.33である。
図13は、各撮像サイトについて、学習用データセットに対するMDDの予測性能(識別器の出力の確率分布)を示す図である。
なお、サイト3(SWA)のデータセットには、単に健常者群があるだけであるが、その確率分布は、他のサイトの健常者群のものに相当している。
図14は、独立した検証データセットにおけるMDDの識別器の出力の確率分布を示す図である。
すなわち、独立した検証データセットを使用して、識別器モデルの汎化性能がテストされる。
独立した検証データセットにおいて、約70%の精度で、生成された識別器モデルは、MDD個体群をHC個体群から分離している。
対応するAUCは、0.75となり、高い識別能力となる能力を示した(並べかえ検定p<0.01)。
感度は68%であり、特異度は71%である。これは、0.38という高いMCC値となっている(並べかえ検定 p<0.01)。
図15は、各撮像サイトについて、独立検証データセットに対するMDDの識別器の出力の確率分布を示す図である。
4つの撮像サイトの全データセットだけでなく個々のデータセットに対して高い分類精度が達成されることがわかる。
以下では、図6で説明した処理のうち、クラスタリングのための特徴量選択と、選択された特徴量によるクラスタリング処理について、より詳しく説明する。
図16は、特徴量を選択して、教師なし学習によりクラスタリングを行う処理を説明するためのフローチャートである。
さらに、疾患識別器生成部3008は、サブサンプルされた各サブサンプルに対して、
特徴量選択を伴う学習方法により識別器を生成する(S208)。
なお、ここでは、図9と同様にして、L1正則化(LASSO)で特徴量選択が実行されるものとする。
このようにして生成された(Ns×Ncv)個の識別器の平均を出力とする統合された識別器を、疾患識別器(診断マーカ)として生成する(S210)。
ここまでの処理は、上述したとおり、図9におけるステップS100~S120と同等の処理である。
ここでは、「特徴量として選択された回数」を、このランキングにおける、その特徴量の重要度と呼ぶことにする。
以上の処理により、クラスタリング分類器生成部3010は、うつ病患者群に対するクラスタリング分類器を生成する(S226)。
すなわち、以上の処理により、クラスタリング分類器生成部3010により、各クラスタについて、観測データから、そのような観測データを生成するような確率分布のモデルがそれぞれ特定され、記憶装置2080に各モデルの情報が疾患識別器データ内に格納される。そして、判別値算出部3012は、クラスタリング分類器として、学習データ以外の入力データに対しては、このよう各確率分布のモデルに基づいて、それぞれのクラスタに当該入力データが属する事後確率を算出し、最大の事後確率となるクラスタに、当該入力データが属するとの分類結果を出力することになる(MAP推定法(Maximum A posteriori Probability estimation method))。
なお、以上の説明では、クラスタリング処理は、疾患識別器生成部3008により、健常者群Nh人、うつ病群Nm人の被験者データに対して実施された識別器の生成処理に基づくものとして説明した。しかし、本実施の形態のクラスタリング方法は、このような場合に限定されることなく、「うつ病患者群」以外の疾患群、たとえば、「統合失調症の患者群」「自閉症の患者群」「強迫神経症の患者群」など、他の精神疾患の患者群のクラスタリングに使用することも可能である。
(アンダーサンプリングおよびサブサンプリング処理)
以上説明した処理では、「アンダーサンプリングおよびサブサンプリング処理」を実行しているので、その技術的意義について、簡単に説明する。
まず、「アンダーサンプリング処理」の効果については、識別器における識別の境界が適切に設定されることが挙げられる。
また、サブサンプリングについては、そもそも、これは、以下の理由により、アンダーサンプリングと併せて、複数回実施することが前提である。
このとき、各識別器の生成において、識別に対する特徴量の重要度が決定される。
ここで、「特徴量選択を伴う学習処理」としては、たとえば、上述したLASSO法のような、いわゆる「スパースモデリング」の手法が挙げられる。
図17を参照して、学習処理の対象となる被験者群には、健常者群と患者群とが含まれるものとする。
健常者群の被験者には、ラベルHが対応付けられており、健常者群には、サブタイプh1とサブタイプh2とが含まれるものとする。
患者群の被験者には、ラベルMが対応付けられており、患者群には、サブタイプm1、サブタイプm2、サブタイプm3とが含まれるものとする。
そして「クラスタリング」の目的は、観測量からこれらのサブタイプへのデータ駆動によるクラスタリングを実行することである。
LASSO法では、このような点線の枠内のグループごとに、1つの特徴量が選択されることで、スパース化が実現していることになる。
図19は、アンダーサンプリングおよびサブサンプリングの処理を複数回実施して、識別器を生成される際に、特徴量が選択される様子を示す概念図である。
[多重共クラスタリング処理]
以下では、図16のステップS224における「多重共クラスタリング」について、その概念を説明し、「多重共クラスタリング」との用語の定義を行う。
前提として、「クラスタリング」とは、コンピュータにより実行される教師なし学習によるデータ分類手法をいい、より特定的には、与えられたデータを外的基準なしに自動的に分類する手法をいうものとする。これに対して、「クラス分類」は、一般に、「教師あり学習」による分類手法のことをいう。また、「クラスタ」とは、内的結合と外的分離の性質を持つデータの部分集合である、と定義される。ここで、外的分離とは、違うクラスタにある対象は類似していないという性質をいい、内的結合とは、同じクラスタ内の対象は互いに類似しているという性質をいう。さらに、「類似」の尺度として、集合の要素間の距離が定義される。一般には、いわゆる「距離の公理」を満たすように距離が定義され、距離としては、ユークリッド距離、マハラビノス距離、シティブロック距離、ミンコフスキー距離などが使われることがある。
以下では、まずは、特徴量に依存して複数のクラスタの分け方が存在する場合を説明するために、クラスタリングを実行する複数の対象を含んだ対象群について、各対象は、複数の特徴量で特徴づけられるものとする。
図20は、特徴量に依存して複数のクラスタの分け方が存在する場合を説明するための概念図である。
そして、これらの文字は、異なる背景パターンや、異なる書体(文字スタイル)を有しているものとする。
したがって、同一の文字の集合を考えた場合でも、いずれの特徴量に基づいて、クラスタリングするかによって、異なるクラスタに分けられることになる。
図21は、複数の特徴量により、複数の対象が特徴づけられる場合のクラスタリングの概念を説明するための概念図である。
図22は、多重クラスタリングおよび多重共クラスタリングを説明するための概念図である。
以上のようにして、特徴量の視点ごとに、異なる対象のクラスタリングを実行することを、図22(b)に示すように「多重クラスタリング」と呼ぶ。
図23は、「多重共クラスタリング」において、1つのビューの中に、異なる種類の確率分布の確率モデルが想定されている場合を示す概念図である。
図23(d)においては、白のブロックと、網掛のブロックでは、異なる種類の確率分布の確率モデルに従っている場合を示している。
たとえば、白のブロックは、連続的な確率変数の確率分布であるのに対して、網掛の部分は、離散的な確率変数の確率分布を想定している場合などを示す。
図24は、多重共クラスタリングの学習方法の概要を説明するためのフローチャートである。
そして、クラスタリング分類器生成部3010は、目的関数の大きさを記憶装置2080に格納する(S310)。
(多重共クラスタリングの処理の詳細)
以下では、図24で説明した多重共クラスタリングの学習方法について、さらに詳しく説明する。
なお、多重共クラスタリングの処理の詳細については、下記の文献に開示があるので、以下では、その概要について説明する。
図25は、図24の多重共クラスタリングの学習方法におけるベイズ推定のグラフ表現を示す図である。
(多重共クラスタリングモデル)
特徴量(脳機能結合値)と被験者(ここでは、患者群の被験者)とは、図21(i)に示したようなデータ行列として表現されるものとする。
そして、データ行列Xは、事前に知られているM個の分布からなる分布族で構成されていると仮定する。
分布族に属する確率分布としては、ガウス分布、ポアソン分布、およびカテゴリ分布/多項分布などが含まれうるものとする。
クラスタリング分類器生成部3010は、X(m)に対して、各データサイズが n× d(m)であるように、以下のように分割する。
X= {X(1),…,X(m),…,X(M) }
異なる分布族に対してこれを組み合わせて、Y= {Y(m)}mとする。
Y j.. (m)が、階層棒折りプロセスによって生成される、分布族mの特徴量jのビュー/特徴量クラスタメンバーシップベクトルを示すものとすると、以下の式が成り立つ。
Mul(・|π)は確率パラメータπを持つ1つのサンプルサイズの多項分布である。
Beta(・|a,b)は、事前サンプルサイズ(a,b)を有するベータ分布である。
Yj.. (m)は1×GVベクトル(Yj,1,1 (m), … ,Y j,V,G (m))Tを示す。
ここでは、所定の条件に従って、十分な大きなVのビューの数とGの特徴量クラスタの数を切り捨てる。 この処理については、たとえば、以下の文献を参照。
https://doi.org/10.1214/06-BA104
(オブジェクトクラスタ Zの生成モデル)
ビュー vのオブジェクトiの被験者クラスタメンバーシップベクトルであって、Zi, v.と表記されるベクトルが以下の式によって生成される。
各インスタンスXi,j (m) は、YとZに条件付きで、独立に特定の分布に従うと仮定する。ビューv、特徴量クラスタg、およびオブジェクトクラスタkのクラスタブロックにおける分布族mのパラメータを θv,g,k (m) と表す。
さらに、Θ={θv,g,k (m)} v,g,k,mと表記して、Xの尤度の対数は以下の式に従う:
尤度は、w={wv}v、w´={w´g,v (m)}g,vおよび u={uk,v}k,vに直接関連付けられていない。
未知変数の結合事前分布φ={Y,Z,w,w´,u,Θ}(すなわち、クラスのメンバーシップ変数とモデルパラメータ)は、次のように与えられる。
YとZのMAP(最大事後)推定には変分ベイズEMアルゴリズムを使用する。
このような変分ベイズEMアルゴリズムについては、以下の文献に開示がある。
対数周辺尤度p(X)はジェンセンの不等式を使用して以下のように近似される。
https://doi.org/10.1007/BF02418571
ここでは、異なるパラメータ(平均場近似)に対して因数分解されるq(φ)を選ぶ。
一般に、KL(Πl=1 L ql(φl), p (φ|x))を最小化する分布 q i(φi)は、以下の式によって与えられる。
この性質を現在検討しているモデルに適用すると、以下を示すことができる。
観測モデルとパラメータΘの事前分布は、後述する。
観測モデルでは、ガウス分布、ポアソン分布、およびカテゴリ分布/多項分布を考慮する。クラスタブロックごとに、クラスタブロック内の特徴が独立していることを前提に、これらのファミリーの1変量分布をフィッティングさせる。これらの分布族のパラメータは、共役事前分布を仮定する。
ハイパーパラメータの更新方程式では、変分ベイズEMアルゴリズムを用いて、次のように計算が実行する。
YとZのMAP推定値は、それぞれ argmaxY q* Y(Y)およびargmaxZ q* Z(Z)として評価される。
下限L(q(φ))は以下の式で与えられる。
多重共クラスタモデルは、ビューの数と特徴/オブジェクトクラスタの数がデータ駆動のアプローチで導き出されるため、さまざまなクラスタリング モデルを表現するのに十分な柔軟性を備えている。たとえば、ビューの数が1の場合、モデルは、共クラスタモデルと一致する。特徴量クラスタの数がすべてのビューに対して1つである場合、多重クラスタリングモデルと一致する。さらに、ビューの数が1であり、特徴量クラスタの数が特徴の数と同じである場合、独立した特徴を持つ従来の混合モデルと一致する。さらに、このモデルは、オブジェクトクラスタを区別しない無情報特徴を検出することができる。このような場合、モデルはオブジェクトクラスタの数が 1つのビューを生成する。モデルの利点は、このような基になるデータ構造を自動的に検出することである。
以上のような「多重共クラスタリング法」により、以下のことが可能となる。
2)この手法によって他の手法では見出せなかったクラスターを同定することが可能となる。
3)さらに、それぞれのクラスターの分け方を特徴量によって意味づけることができ、個々のクラスターの解釈を容易にできる。
以下では、上述したようなSRPBSとして公開されている多数の被検者から収集された多施設大規模fMRIデータを2つに分割して、それぞれを上述したような多重共クラスタリング法で使用することで、クラスタリングについての汎化性能を検証する。
図26は、このように2つに分割されたデータセット1とデータセット2とを示す図である。
図27は、各データセットにおいてクラスタリングを実施する概念を説明する概念図である。
図28は、被験者データに対する多重共クラスタリングの例を示す概念図である。
図28(a)に示すように、入力となるデータ行列として、行方向に被験者が並び、列方向に特徴量が並んでいるものとする。
図29は、データセット1およびデータセット2に対して、実際に多重共クラスタリング処理を実施した結果を示す図である。
図29においては、クラスタリングのための特徴量としては、データセット1およびデータセット2とも、それぞれ99個が選択されている。
図30は、データセット1とデータセット2において、それぞれのビューに割り当てられた脳機能結合(FC)の個数を示す表である。
データセット1では、ビュー1に特徴量としてFCが92個、ビュー2にFCが7個割り当てられる。
データセット2では、ビュー1に特徴量としてFCが93個、ビュー2にFCが6個割り当てられる。
(クラスタリング(層別化)の汎化性(データセット間類似性)検証方法)
図31は、このようなクラスタリングの類似度(層別化の汎化性能)の評価手法を説明するための概念図である。
(クラスタリング間の類似度(再現度)を測る評価尺度)
図32は、ARIについて説明するための概念図である。
なお、ARIについては、たとえば、以下の文献に開示がある。
図33は、クラスタリング1とクラスタリング1´の間、および、クラスタリング2とクラスタリング2´の間の類似度の評価結果を示す図である。
図33(a)は、データセット1とデータセット2について、それぞれのビューに対するARIを算出した表である。
図34は、クラスタリング1とクラスタリング1´において、それぞれビュー1の各クラスタに割り当てられる被験者数の分布を示す表である。
以下では、以下の文献に開示され、図6において、ハーモナイゼーション処理と呼んでいる処理の内容について、説明する。
[トラベリングサブジェクト法のハーモナイゼーション]
取得されたトラベリングサブジェクトのデータセットは、データセンター200における記憶装置210に移動被験者データとして格納される。
そして、後述するように、計算処理システム300において、「脳活動バイオマーカーのハーモナイズ方法」のための処理が実行される。
たとえば、GLMMについては、以下の文献に開示がある。
公知文献8:久保 拓弥著、「データ解析のための統計モデリング入門」、岩波書店、2012年 第1刷、2017年 第14刷
i) 各々の機能的なコネクティビティに対する相関値のすべてのサイトにわたる平均からの偏差として、各サイトの測定バイアスを定義する。
iii) 疾患要因は、健常者群の値からの偏差として定義する。
以下では、ベクトルは小文字(例えば、m)によって表示され、ベクトルはすべて列ベクトルであると仮定される。
ベクトルの要素は、mkのように添字によって表示される。
脳領野間のn個の相関値からなる機能的コネクティビティベクトル(列ベクトルとする)の回帰モデルは、以下の式のように表される。
参加者がどんなクラス(健常者、患者、トラベリングサブジェクト)にも属さなければ、ターゲットベクトルは、すべての要素が0に等しいベクトルである。
上付きの添字Tは、マトリックスかベクトルの置換を表示し、xTは行ベクトルを表わす。
各機能的なコネクティビティの相関値に対しては、回帰モデルの計画行列がランク不足であるので、L2正規化による最小二乗回帰を使用して、それぞれのパラメータを評価する。なお、L2正規化による最小二乗回帰法以外にも、たとえば、ベイズ推定法など別の評価方法を用いることも可能である。
上記のような回帰計算の後には、被験者aのb番目のコネクティビティは、以下のように記述できる:
図36においては、1項目と2項目のターゲットベクトルの意義と、測定バイアスベクトルおよび健常者の標本バイアスベクトルについて示している。
第3項目以降についても同様である。
(ハーモナイゼーション処理のフロー)
図37は、ハーモナイゼーションのために、測定バイアスを算出する処理を説明するためのフローチャートである。
判別処理部3000が、被験者に対する疾患あるいは健常ラベルについての疾患識別器を生成する処理における脳機能結合値に対するハーモナイゼーションについて簡単に説明する。
このような疾患識別器は、被験者を診断するための補助情報(支援情報)を提供する。
補正処理が実行された後のデータは、補正後相関値データ3110として、記憶装置2080に格納される。
実施の形態1では、脳活動計測装置(fMRI装置)で、複数の測定場所で計測された脳活動データを計測し、この脳活動データに基づいて、バイオマーカーの生成およびバイオマーカーによる診断ラベルの推定(予測)を行う構成として、分散処理により行う一例の構成を説明した。
図38は、データ収集、推定処理および対象者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。
計算処理システム300は、サーバ200に格納されたデータから識別器を生成する。
MRI装置410は、計算処理システム300上の識別器の結果を利用する別サイトに設けられており、特定の被験者について脳活動のデータを計測する。
コンピュータ400では、図示しない表示装置などを介して、ユーザに対して、判別結果を知らせる。
このような構成とすることで、より多くの被験者について収集したデータに基づいて、識別器による診断ラベルの推定結果を提供することが可能となる。
また、上述のとおり、このような判別は、ある属性に属する可能性を確率として表示するものであってもよい。
[治療法選択支援システム]
(治療法選択支援システムの構成)
本発明のある実施形態は、うつ症状を呈する被験者の脳活動の計測結果に基づいて、この被験者に対する治療法の選択に関する情報を提供するための治療法選択支援システム1000に関する。
図39は、治療法選択支援システム1000の機能構成を示す図である。
なお、図38と同様に、支援情報提供装置300aと分類器生成装置300bとは、同一のコンピュータシステム上の機能として実装されてもよい。
その他、基本的に、図38の構成と同一部分には、同一符号を付している。
図40は、治療法情報データベース5100の例を示す図である。
うつ病の物理的治療法として、経頭蓋磁気刺激法、ニューロフィードバック法、電気けいれん療法や認知行動療法等を挙げることができる。
次に、支援情報提供装置300aの治療法情報生成部3200が行う治療法選択支援処理について説明する。治療法選択支援処理は、治療法選択支援プログラムをコンピュータで治療法情報生成部3200の処理として実行することにより、達成される。
図41は、被験者の治療法選択の支援処理のフローを説明するフローチャートである。
以下では、図39で説明した支援情報提供装置300aを、創薬における被験者のスクリーニングを実行するための支援装置として使用する態様について説明する。
図42は、一般的な創薬のプロセスを示す図である。
以上の過程を経て、非臨床試験を通過できた候補物質について、いわゆる「治験」が実施される。
以下では、支援情報提供装置300aを、このような被験者集団の同定(被験者スクリーニング)の支援装置として使用する場合について説明する。
すなわち、支援情報提供装置300aから、スクリーニング支援データとして、分類されたクラスタに関する情報をコンピュータ400に返信する構成とする。
図44は、スクリーニング支援装置1000´の構成を説明する図である。
スクリーニング支援装置1000´は、治療法選択支援システム1000と同様に、医療機関4において、診断ラベルの推定対象となる被験者の脳活動データを計測するためのMRI装置410からのデータを受信可能なように接続されたコンピュータ400と、支援情報提供装置300aと、クラスタリング分類器の生成処理を実施する分類器生成装置300bと、データサーバ200´を備える。
試験の終了後など、試験結果を評価する際には、仮IDに基づいて、試験結果と被験者についての分類結果との照合などが実施されるものとすることができる。
次に、支援情報提供装置300aの演算装置2040aが行うスクリーニング支援処理について説明する。スクリーニング支援処理はスクリーニング支援プログラムをコンピュータで実行することにより、達成される。
図43は、創薬における被験者のスクリーニング支援処理のフローを説明するフローチャートである。
医師は、予め定められた治験プロトコルに従って、対応するクラスタの被験者をスクリーニングして、治験を実施することが可能となる。
以上、図39においては、治療法選択支援システム1000について説明し、図42~図44においては、スクリーニング支援装置1000´について説明した。
本検証データに使用したクラスタリング分類器の生成にあたり疾患識別器生成部3008が生成するMDD識別器の構成は、以下のとおりである。
まず、脳領野のパーセレーション法としては、BSA法(Brainvisa Sulci Atlas)に従って脳領域を分割した。
たとえば、BSA法については、以下の文献に開示がある:
そして、クラスタリング分類器生成部3010におけるクラスタリング分類器の生成処理については、以下のとおりである。
ここで、「データセット間での安定性」の評価には、図33で説明したARI値を評価指標として用いることができる。
図46は、データセット2について多重共クラスタリングによるクラスタリングの結果のビューを示す図である。
図47は、データセット1とデータセット2との間のクラスタリング安定性を示す図である。
図47は、データセット1とデータセット2について、それぞれのビューに対するARIを算出した表であり、図33に示した表に対応するものである。
図48は、全うつ病患者データ(データセット1+2)に対するクラスタリング結果のビューを示す図である。
ここで、「臨床データ」とは、各患者への「投薬の履歴の情報」、各患者の「うつ病の程度の診断情報」などを指す。
図50は、クラスタリングに使用された脳機能結合 (ビュー3)を示す図である。
図50(a)は、クラスタリングに使用された脳機能結合の脳内の位置を示し、図50(b)は、機能結合の関心領域を示す。
ビュー3においては、図50(a)(b)に示された3つの結合のみが、クラスタリングに使用されている。
多重共クラスタリングでは、結合をどのように分けて患者(被験者)をどのようにクラスタリングするかはアルゴリズムが自動的に最適解を導くことになる。
なお、図50(b)における関心領域の名称は、以下の内容をそれぞれ示す。
Thalamus_LorR:左or右視床
Precentral_R:右中心前回
Postcentral_L:左中心後回
図51は、図48、図49に示したビュー3の各サブタイプについて、うつ病の重症度及び重症度の改善率との関係を示す図である。
図51(a)は、治療開始後、第0週と第6週のHAMDのスコアを対比して示す図である。
また、図51(a)において、図49においてサブタイプ1~5とは、図49の「被検者クラスタリングの結果 (ビュー3)」の上位5つのクラスタである。
ここでは、研究における治療プロトコルとしては、以下のようなものである。
図51(b)は、第0週と第6週を比較した場合のHAMDの改善率を示す図である。
ここで、改善率とは、以下の式で表される。
(HAMD(0) - HAMD(6W)) / HAMD(0) x 100
図51(a)より、まず初期の段階ではサブタイプ間でHAMDに差がない、つまりうつ重症度は同程度であることが分かる。
Claims (28)
- うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援システムであって、
複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング装置を備え、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
前記クラスタリング装置は、
前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含み、
前記演算装置は、クラスタリング分類器の生成処理において、
i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、
ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、
iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成し、
前記治療法選択支援システムは、さらに、
前記クラスタリング分類器による層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置と、
前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果に応じて、対応する治療法情報を出力する支援情報提供装置とを、備える、治療法選択支援システム。 - 前記演算装置は、前記識別器モデルを生成する機械学習において、
前記第1の群と前記第2の群から、アンダーサンプリングおよびサブサンプリングを実行して、複数の学習用サブサンプルを生成し、
前記学習用サブサンプルのそれぞれについて、機械学習による識別器の生成において使用される特徴量の和集合から、前記和集合に属する特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
選択された前記クラスタリングのための特徴量に基づいて、前記多重共クラスタリング法により、前記クラスタリング分類器を生成する、請求項1記載の治療法選択支援システム。 - 前記支援情報提供装置は、クラスタリング演算装置とインタフェース装置とを備え、
前記クラスタリング演算装置は、前記クラスタリング分類器により前記第1の被験者が各前記クラスタに属する確率を算出し、前記データベース装置から、前記確率に応じて選択される少なくとも2つの前記治療法情報を読み出し、
前記インタフェース装置は、前記選択されたクラスタとそれぞれ対応する前記治療法情報とを関連付けて表示するためのデータを出力する、請求項1または2に記載の治療法選択支援システム。 - 前記治療法情報は、特定の治療薬に対する応答性を示す情報である、請求項1~3のいずれか1項に記載の治療法選択支援システム。
- 前記治療法情報は、特定の物理的治療法に対する応答性を示す情報である、請求項1~3のいずれか1項に記載の治療法選択支援システム。
- 前記機械学習による識別器の生成の処理は、前記複数の学習用サブサンプルに対してそれぞれ複数の識別器サブモデルを生成し、前記複数の識別器サブモデルを統合して前記識別器モデルを生成するアンサンブル学習である、請求項2に記載の治療法選択支援システム。
- 前記クラスタリング装置は、複数の計測サイトにそれぞれ設けられた複数の脳活動計測装置から、各前記複数の第2の被験者の所定の複数の脳領野ペア間の脳活動の時間相関を表現する情報を受け、
前記演算装置は、各前記複数の第2の被験者について前記複数の脳機能結合相関値について、前記計測サイトの測定バイアスを除去するように補正することで、補正された調整値を前記特徴量として前記記憶装置に格納するハーモナイゼーション算出手段を含む、請求項1または2に記載の治療法選択支援システム。 - うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援装置であって、
複数の第2の被験者のうち、うつ病の診断ラベルを有する被験者に対する層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置と、
前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に基づく層別化の結果に応じて、対応する治療法情報を出力する支援情報提供装置とを備え、
前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
前記層別化の結果のクラスタは、クラスタリング装置による脳機能結合相関値の計測結果に対するクラスタリング処理により得られたクラスタリング分類器により得られるものであり、
前記クラスタリング装置は、
前記第1の群について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含み、
前記クラスタリング分類器の生成処理において、前記演算装置が、
i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、
ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、
iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングすることにより、前記クラスタリング分類器を生成する、治療法選択支援装置。 - 前記支援情報提供装置は、クラスタリング演算装置とインタフェース装置とを備え、
前記クラスタリング演算装置は、前記クラスタリング分類器により前記第1の被験者が各前記クラスタに属する確率を算出し、前記データベース装置から、前記確率に応じて選択される少なくとも2つの前記治療法情報を読み出し、
前記インタフェース装置は、前記選択されたクラスタとそれぞれ対応する前記治療法情報とを関連付けて表示するためのデータを出力する、請求項8に記載の治療法選択支援装置。 - 前記治療法情報は、特定の治療薬に対する応答性を示す情報である、請求項8または9に記載の治療法選択支援装置。
- 前記治療法情報は、特定の物理的治療法に対する応答性を示す情報である、請求項8~10のいずれか1項に記載の治療法選択支援装置。
- うつ症状を呈する第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援方法であって、
複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング分類器を生成して準備するステップを備え、前記複数の第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の被験者に対する治療法の選択に関する情報を提供するための治療法選択支援プログラムであって、
前記治療法選択支援プログラムは、コンピュータに実行させたときに、コンピュータに、
複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング分類器を生成するステップと、
前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果に応じて、前記クラスタリング分類器による層別化の結果のクラスタと対応する所定の治療法情報とを関連付けて格納するためのデータベース装置から、対応する治療法情報を取得し出力するステップとを実行させ、
前記複数の第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の被験者に対するスクリーニングを支援するためのスクリーニング支援システムであって、
複数の第2の被験者から取得した脳機能結合相関値の計測結果に対して、クラスタリング処理により複数のクラスタに分ける層別化を実行するためのクラスタリング装置を備え、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
前記クラスタリング装置は、
前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含み、前記演算装置は、
i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、
ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、
iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成し、
前記スクリーニング支援システムは、さらに、
前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器による分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力する支援情報提供装置とを、備える、スクリーニング支援システム。 - 前記演算装置は、前記識別器モデルを生成する機械学習において、
前記第1の群と前記第2の群から、アンダーサンプリングおよびサブサンプリングを実行して、複数の学習用サブサンプルを生成し、
前記学習用サブサンプルのそれぞれについて、機械学習による識別器の生成において使用される特徴量の和集合から、前記和集合に属する特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
選択された前記クラスタリングのための特徴量に基づいて、前記多重共クラスタリング法により、前記クラスタリング分類器を生成する、請求項16記載のスクリーニング支援システム。 - うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援装置であって、
クラスタリング分類器を特定する情報を格納するための記憶装置を有する支援情報提供装置であって、前記第1の被験者の脳活動の計測結果を入力として受け、前記計測結果に対する前記クラスタリング分類器に基づく分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力する支援情報提供装置を備え、
層別化の結果のクラスタは、クラスタリング装置による、複数の第2の被験者から取得した脳機能結合相関値の計測結果に対するクラスタリング処理により得られた前記クラスタリング分類器により得られるものであり、前記複数の第2の被験者は、うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含み、
前記クラスタリング装置は、
前記複数の第2の被験者について、前記クラスタリング処理を実行するための演算装置と記憶装置とを含み、
前記クラスタリング分類器の生成処理において、前記演算装置が、
i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を前記記憶装置に格納し、
ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行し、
iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択し、
iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングして前記クラスタリング分類器を生成する、スクリーニング支援装置。 - うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援方法であって、
記憶装置に格納された情報で特定されるクラスタリング分類器に基づいて、前記脳活動の計測結果により、演算装置が前記第1の被験者の分類を実行するステップと、
前記分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力するステップと、を備え、
前記クラスタリング分類器の生成のための処理は、
うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む複数の第2の被験者について、クラスタリング処理を実行するための演算ステップを含み、前記演算ステップは、
i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を取得するステップと、
ii)前記取得された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、
iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、
iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成するステップとを含む、スクリーニング支援方法。 - うつ症状に対する治療手段候補の臨床試験において、第1の被験者の脳活動の計測結果に基づいて、前記第1の被験者に対するスクリーニングを支援するためのスクリーニング支援プログラムであって、
前記スクリーニング支援プログラムは、コンピュータに実行させたときに、コンピュータに、
記憶装置に格納された情報で特定されるクラスタリング分類器に基づいて、前記脳活動の計測結果により、演算装置が前記第1の被験者の分類を実行するステップと、
前記分類結果を前記第1の被験者と関連付けて記録し、前記分類結果に基づく前記第1の被験者のスクリーニングを支援する情報を出力するステップと、
を実行させ、
前記クラスタリング分類器の生成のための処理は、
うつ病の診断ラベルを有する第1の群と前記うつ病の診断ラベルを有さない第2の群とを含む複数の第2の被験者について、クラスタリング処理を実行するための演算ステップを含み、前記演算ステップにおいて、
i)各前記複数の第2の被験者について、所定の複数の脳領野ペア間の脳活動の時間相関をそれぞれ表す複数の脳機能結合相関値に基づく特徴量を記憶装置に格納するステップと、
ii)前記記憶装置に格納された前記特徴量に基づいて、前記診断ラベルの有無を判別するための識別器モデルを生成する機械学習を教師あり学習で実行するステップと、
iii)前記識別器モデルを生成する機械学習において、機械学習による識別器の生成において使用される特徴量の重要度に応じて、クラスタリングのための特徴量を選択するステップと、
iv)選択された前記クラスタリングのための特徴量に基づいて、教師なし学習の多重共クラスタリング法により、前記第1の群をクラスタリングしてクラスタリング分類器を生成するステップを含む、スクリーニング支援プログラム。 - 前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である、請求項1に記載の治療法選択支援システム。
- 前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である、請求項8に記載の治療法選択支援装置。
- 前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である、請求項12または13に記載の治療法選択支援方法。
- 前記所定の治療法情報は、選択的セロトニン再取り込み阻害薬に対する治療応答性に関する情報である、請求項14または15に記載の治療法選択支援プログラム。
- 前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である、請求項16に記載のスクリーニング支援システム。
- 前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である、請求項18に記載のスクリーニング支援装置。
- 前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である、請求項19に記載のスクリーニング支援方法。
- 前記治療手段候補は、選択的セロトニン再取り込み阻害薬を用いた治療法である、請求項20に記載のスクリーニング支援プログラム。
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)
| 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)
| 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)
| 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 |
-
2021
- 2021-07-15 CN CN202180060576.5A patent/CN116133587A/zh active Pending
- 2021-07-15 WO PCT/JP2021/026658 patent/WO2022014682A1/ja not_active Ceased
- 2021-07-15 US US18/016,161 patent/US20230284983A1/en active Pending
- 2021-07-15 EP EP21842607.0A patent/EP4183333A4/en active Pending
- 2021-07-15 JP JP2022536448A patent/JP7737099B2/ja active Active
Patent Citations (4)
| 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 |