JP2019063478A - 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置 - Google Patents

判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置 Download PDF

Info

Publication number
JP2019063478A
JP2019063478A JP2018005110A JP2018005110A JP2019063478A JP 2019063478 A JP2019063478 A JP 2019063478A JP 2018005110 A JP2018005110 A JP 2018005110A JP 2018005110 A JP2018005110 A JP 2018005110A JP 2019063478 A JP2019063478 A JP 2019063478A
Authority
JP
Japan
Prior art keywords
functional
classifier
subject
correlation
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.)
Granted
Application number
JP2018005110A
Other languages
English (en)
Other versions
JP2019063478A5 (ja
JP6722845B2 (ja
Inventor
ジュゼッペ リシ
Lisi Giuseppe
ジュゼッペ リシ
森本 淳
Jun Morimoto
淳 森本
光男 川人
Mitsuo Kawato
光男 川人
貴志 山田
Takashi Yamada
貴志 山田
奈穂 市川
Nao Ichikawa
奈穂 市川
泰昌 岡本
Yasumasa Okamoto
泰昌 岡本
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
Original Assignee
Hiroshima University NUC
ATR Advanced Telecommunications Research Institute International
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 filed Critical Hiroshima University NUC
Priority to PCT/JP2018/036952 priority Critical patent/WO2019069955A1/ja
Priority to US16/753,291 priority patent/US12020427B2/en
Priority to EP18864859.6A priority patent/EP3692909A4/en
Priority to CN201880078144.5A priority patent/CN111447874A/zh
Publication of JP2019063478A publication Critical patent/JP2019063478A/ja
Publication of JP2019063478A5 publication Critical patent/JP2019063478A5/ja
Application granted granted Critical
Publication of JP6722845B2 publication Critical patent/JP6722845B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/243Classification techniques relating to the number of classes
    • G06F18/2431Multiple classes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • 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
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • 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/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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
    • 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
    • G16H70/00ICT specially adapted for the handling or processing of medical references
    • G16H70/60ICT specially adapted for the handling or processing of medical references relating to pathologies
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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; CALCULATING OR 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; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Public Health (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Quality & Reliability (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Mathematical Physics (AREA)
  • Social Psychology (AREA)
  • Psychology (AREA)
  • Computing Systems (AREA)
  • Psychiatry (AREA)
  • Hospice & Palliative Care (AREA)
  • Developmental Disabilities (AREA)
  • Child & Adolescent Psychology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)

Abstract

【課題】脳の活動状態に対して、うつ症状の客観的な疾患ラベルを判別できる判別装置を提供する。【解決手段】被検者におけるうつ症状の有無の判定を補助するための判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備える。分類器は、複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、うつ症状の疾患ラベルを判別するように生成される。判別装置は、演算装置をさらに備え、演算装置は、分類器により、被検者のうつ症状に対する分類結果を生成する判別処理を実行する。【選択図】図15

Description

本発明は、判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置に関する。
(バイオマーカー)
生体内の生物学的変化を定量的に把握するため、生体情報を数値化・定量化した指標のことを「バイオマーカー」と呼ぶ。
FDA(米国食品医薬品局)はバイオマーカーの位置づけを、「正常なプロセスや病的プロセス、あるいは治療に対する薬理学的な反応の指標として客観的に測定・評価される項目」としている。また疾患の状態や変化、治癒の程度を特徴づけるバイオマーカーは、新薬の臨床試験での有効性を確認するためのサロゲートマーカー(代替マーカー)として使われる。血糖値やコレステロール値などは、生活習慣病の指標として代表的なバイオマーカーである。尿や血液中に含まれる生体由来の物質だけでなく、心電図、血圧、PET画像、骨密度、肺機能なども含まれる。またゲノム解析やプロテオーム解析が進んできたことによって、DNAやRNA、生体蛋白等に関連したさまざまなバイオマーカーが見出されている。
バイオマーカーは、疾患にかかった後の治療効果の測定だけでなく、疾患を未然に防ぐための日常的な指標として疾患の予防に、さらに副作用を回避した有効な治療法を選択する個別化医療への応用が期待されている。
ただし、神経・精神疾患の場合、生化学的もしくは分子遺伝学的観点から客観的な指標として利用可能な分子マーカーなども研究されているものの、検討段階というべき状況である。
一方で、NIRS(Near-infraRed Spectroscopy)技術を用いて、生体光計測により計測されたヘモグロビン信号から特徴量に応じて、統合失調症、うつ病などの精神疾患について分類を行う疾患判定システムなども報告がある(特許文献1)。
バイオマーカーは、通常の診療だけでなく、医薬品開発でも利用されている。
一般的な創薬の場合、探索研究の段階では薬力学マーカーが用いられ、非臨床段階になると候補化合物の毒性をみる毒性マーカーも利用されるようになる。
臨床試験の段階になると、評価項目として代替マーカーが利用されたり、診断マーカー、予後マーカー、予測マーカー、モニタリングマーカー、安全性マーカーなど、多様なマーカーが利用される。
「代替マーカー」とは、真のエンドポイントとの科学的な関係が証明されているマーカーのことをいう。
「診断マーカー」とは、疾患の有無や疾患の種類を特定する診断マーカーと,疾患の重症度や進行度を判定する診断マーカーの2種類がある。「予測マーカー」とは、特定の薬剤が効きやすい集団を見分けるバイオマーカーである。「予後マーカー」とは、特定の治療に関係なく疾患の進行や回復を予想するバイオマーカーである。「モニタリングマーカー」とは、「薬効マーカー」とも呼ばれ、薬剤の効き具合,治療効果を診るマーカーである。このマーカーを長期間モニタリングすることで薬剤耐性の原因解明のヒントにもなりうる、とされる。
たとえば、標的分子が明確である場合などは、たとえば、標的分子などを、特定の薬剤が効きやすい患者の集団を見分けて選択するために、臨床試験で患者層別マーカーとして利用できる場合がある。このような場合、患者層別マーカーを利用する医薬品が承認されると、その患者層別マーカーを測定するための診断薬は「コンパニオン診断薬」と呼ばれる場合がある。あるいは、このような患者層別マーカーにより、患者に対する治療方法を選択するために行われる診断は、「コンパニオン診断」と呼ばれる場合もある。
臨床試験において、後期試験にかかるコストが非常に高いことから、後期試験に適したバイオマーカーの開発が臨床試験の効率を向上させるために必要と考えられる。しかし、このようなバイオマーカーの開発は、一般には容易ではない。
また、たとえば、生物学的指標 (バイオマーカー)を考える際、疾患病態生理において、疾患原因に関連する行動学あるいは生物学的プロセスを反映しているものはtrait markerであり、患者における臨床状態像を反映しているのがstate markerと呼ばれる。trait markerは、遺伝規定性指標とも呼ばれ、state markerは、状態依存性指標とも呼ばれる。
(リアルタイムニューロフィードバック)
たとえば、従来、神経症の一型である強迫性障害 (Obsessive-Compulsive Disorder: OCD)に対しては、その治療法として、薬物療法や行動療法が知られている。薬物療法としては、たとえば、選択的セロトニン再取り込阻害薬が用いられ、行動療法としては、曝露法と反応妨害法を組み合わせた曝露反応妨害法などが知られている。
一方で、神経・精神疾患の場合の治療法の可能性として、リアルタイムニューロフィードバック法が検討されている(特許文献6)。
核磁気共鳴画像法(MRI:Magnetic Resonance Imaging)を利用して、ヒトの脳の活動に関連した血流動態反応を視覚化する方法である機能的磁気共鳴画像法(fMRI:functional Magnetic Resonance Imaging)を始めとする脳機能画像法は、感覚刺激や認知課題遂行による脳活動と安静時や対照課題遂行による脳活動の違いを検出して、関心のある脳機能の構成要素に対応する脳賦活領域を特定すること、すなわち脳の機能局在を明らかにすることにもちいられてきた。
近年、機能的磁気共鳴画像法(fMRI)などの脳機能画像法を用いたリアルタイムニューロフィードバック技術についての報告がある(非特許文献1)。リアルタイムニューロフィードバック技術は、神経疾患、精神疾患の治療方法として注目され始めている。
ニューロフィードバックとは、バイオフィードバックの一種であり、被検者が自身の脳活動に関するフィードバックを受けることで、脳活動を操作する方法を学ぶというものである。
たとえば、前部帯状回の活動をfMRIで計測し、それをリアルタイムで炎の大きさとして、患者にフィードバックして、炎を小さくするように努力させると、中枢性の慢性疼痛が実時間でも、また、長期的にも改善したとの報告がある(非特許文献2を参照)。
(安静時fMRI)
また、最近の研究では、安静時の脳も活発に活動していることがわかってきた。つまり、アクティブな活動中には鎮静化しており、休息時には活発に興奮する神経細胞群が脳内に存在する。解剖学的には左右の大脳が合わさった内側面が主な部位で、前頭葉内側面、後部帯状回、楔前部、頭頂連合野の後半部、中側頭回などである。この安静時の脳活動のベースラインを示す領域はデフォルトモードネットワーク(Default Mode Network;DMN)と名付けられ、ひとつのネットワークとして同期しながら活動する(非特許文献3を参照)。
たとえば、健常者と精神疾患の患者で脳活動が異なるという例として、デフォルトモードネットワークにおける脳活動が挙げられる。デフォルトモードネットワークとは、安静状態において、目標指向的な課題を実行中よりも活発な脳活動がみられる部位のことを指す。統合失調症、アルツハイマー病などの疾患の患者は、健常者と比較したときに、このデフォルトモードネットワークに異常がみられることが報告されている。たとえば、統合失調症の患者は、安静状態において、デフォルトモードネットワークに属する後帯状皮質と、頭頂外側部、前頭前野内側部、小脳皮質の活動の相関が下がっているなどの報告がある。
ただし、このようなデフォルトモードネットワークが、認知機能とどのようにかかわっているのかや、このような脳領域間の機能的結合の相関と上述したニューロフィードバックとの関連などは、必ずしも明らかとはなっていない状況である。
一方で、複数の脳領域間において課題などの違いによる活動の相関関係の変化をみることにより、それら脳領域間の機能的結合(functional connectivity)が評価されている。とくに、安静時のfMRIによる機能的結合の評価は安静時機能的結合MRI(rs-fcMRI : resting-state functional connectivity MRI)とも呼ばれ、様々な神経・精神疾患を対象とした臨床研究も行われるようになりつつある。ただし、従来のrs−fcMRI法は、上述したようなデフォルトモードネットワークといったような大域的神経ネットワークの活動をみるものであって、より詳細な機能的結合について、十分な検討がされているとはいえない状況である。
(核磁気共鳴画像法)
このような核磁気共鳴画像法について、簡単に説明すると、以下のとおりである。
すなわち、従来から、生体の脳や全身の断面を画像する方法として、生体中の原子、特に、水素原子の原子核に対する核磁気共鳴現象を利用した核磁気共鳴映像法が、人間の臨床画像診断等に使用されている。
核磁気共鳴映像法は、それを人体に適用する場合、同様の人体内断層画像法である「X線CT」に比較して、たとえば、以下のような特徴がある。
(1)水素原子の分布と、その信号緩和時間(原子の結合の強さを反映)に対応した濃度の画像が得られる。このため、組織の性質の差異に応じた濃淡を呈し、組織の違いを観察しやすい。
(2)磁場は、骨による吸収がない。このため、骨に囲まれた部位(頭蓋内、脊髄など)を観察しやすい。
(3)X線のように人体に害になるということがないので、広範囲に活用できる。
このような核磁気共鳴影像法は、人体の各細胞に最も多く含まれ、かつ最も大きな磁性を有している水素原子核(陽子)の磁気性を利用する。水素原子核の磁性を担うスピン角運動量の磁場内での運動は、古典的には、コマの歳差運動にたとえられる。
以下、本発明の背景の説明のために、この直感的な古典的モデルで、簡単に核磁気共鳴の原理をまとめておく。
上述したような水素原子核のスピン角運動量の方向(コマの自転軸の方向)は、磁場のない環境では、ランダムな方向を向いているものの、静磁場を印加すると、磁力線の方向を向く。
この状態で、さらに振動磁界を重畳すると、この振動磁界の周波数が、静磁界の強さで決まる共鳴周波数f0=γB0/2π(γ:物質に固有の係数)であると、共鳴により原子核側にエネルギーが移動し、磁化ベクトルの方向が変わる(歳差運動が大きくなる)。この状態で、振動磁界を切ると、歳差運動は、傾き角度を戻しながら、静磁界における方向に復帰していく。この過程を外部からアンテナコイルにより検知することで、NMR信号を得ることができる。
このような共鳴周波数f0は、静磁界の強度がB0(T)であるとき、水素原子では、42.6×B0(MHz)となる。
さらに、核磁気共鳴映像法では、血流量の変化に応じて、検出される信号に変化が現れることを用いて、外部刺激等に対する脳の活動部位を視覚化することも可能である。このような核磁気共鳴映像法を、特に、機能的磁気共鳴画像法(fMRI:functional MRI)と呼ぶ。
fMRIでは、装置としては通常のMRI装置に、さらに、fMRI計測に必要なハードおよびソフトを装備したものが使用される。
ここで、血流量の変化がNMR信号強度に変化をもたらすのは、血液中の酸素化および脱酸素化ヘモグロビンは磁気的な性質が異なることを利用している。酸素化ヘモグロビンは反磁性体の性質があり、周りに存在する水の水素原子の緩和時間に影響を与えないのに対し、脱酸素化ヘモグロビンは常磁性体であり、周囲の磁場を変化させる。したがって、脳が刺激を受け、局部血流が増大し、酸素化ヘモグロビンが増加すると、その変化分をMRI信号として検出する事ができる。被検者への刺激は、たとえば、視覚による刺激や聴覚による刺激、あるいは所定の課題(タスク)の実行等が用いられる(たとえば、特許文献2)。
ここで、脳機能研究においては、微小静脈や毛細血管における赤血球中の脱酸素化ヘモグロビンの濃度が減少する現象(BOLD効果)に対応した水素原子の核磁気共鳴信号(MRI信号)の上昇を測定することによって脳の活動の測定が行われている。
特に、人の運動機能に関する研究では、被検者に何らかの運動を行わせつつ、上記MRI装置によって脳の活動を測定することが行われている。
ところで、ヒトの場合、非侵襲的な脳活動の計測が必要であり、この場合、fMRIデータから、より詳細な情報を抽出できるデコーディング技術が発達してきている(たとえば、非特許文献4)。特に、fMRIが脳におけるボクセル単位(volumetric pixel : voxel)で脳活動を解析することで、脳活動の空間的パターンから、刺激入力や認識状態を推定することが可能となっている。
さらに、このようなデコーディング技術を発展させた技術として、特許文献3には、神経・精神疾患に対して、脳機能画像法によるバイオマーカーを実現するための脳活動解析方法が開示されている。この方法では、健常群、患者群において測定された安静時機能的結合MRIのデータから、それぞれの被験者について、所定の脳領域間の活動度の相関行列を導出する。被験者の疾患/健常ラベルを含む被験者の属性と相関行列とについて正則化正準相関解析により特徴抽出が行われる。正則化正準相関解析(スパース正準相関解析)の結果に基づいて、スパースロジスティック回帰(Sparse Logistic Regression: SLR)による判別分析によりバイオマーカーとして機能する分類器が生成される。このような機械学習の技術により、安静時のfMRI データから導き出される脳領野間の結合(機能的結合)にもとづいて自閉症の診断結果を予測可能であることが示された(特許文献4、特許文献5)。しかも、その予測性能の検証は一つの施設において計測された脳活動のみではなく、他の施設で計測された脳活動に対しても、ある程度の汎化が可能であることが示された。
さらに、近年では、多数の研究施設から得られた大規模なマルチサイトサンプル群(n=1188)でfMRIを用いて、うつ病患者は4つの神経生理学的サブタイプ(バイオタイプ)に分類できること、こうしたバイオタイプは辺縁系および前頭葉-線条体ネットワークにおける結合性機能不全を示すそれぞれ別のパターンによって明確に示されることが指摘されている(非特許文献5を参照)。非特許文献5においては、これらのバイオタイプは、臨床的特徴を基盤とするだけでは弁別できないが、異なる臨床症状プロファイルと結び付いており、また経頭蓋磁気刺激治療に対する応答性を予測する(n=154)、とされている。
非特許文献5においては、脳内の258の領域からのBOLD信号に基づき、脳内の機能的な結合により、上記のサブタイプの分類を行っている。
再表2006−132313号公報 特開2011−000184号公報 特開2015−62817号公報 特許第6195329号公報 国際公開第2017/090590号 特許第5641531号公報
Nikolaus Weiskopf ,"Real-time fMRI and its application to neurofeedback",NeuroImage 62 (2012) 682-692 deCharms RC, Maeda F, Glover GH et al, "Control over brain activation and pain learned by using real-time functional MRI",Proc Natl Acad Sci USA 102(51), 18626-18631, 2005 Raichle ME, Macleod AM, Snyder AZ. et. al. "A default mode of brain function",Proc Natl Acad Sci USA 98(2),676-682,2001 Kamitani Y, Tong F. Decoding the visual and subjective contents of the human brain. Nat Neurosci. 2005; 8: 679-85. Drysdale AT, Grosenick L, Downar J, Dunlop K, Mansouri F, Meng Y, Fetcho RN, Zebley B, Oathes DJ, Etkin A, Schatzberg AF, Sudheimer K, Keller J, Mayberg HS, Gunning FM, Alexopoulos GS, Fox MD, Pascual-Leone A, Voss HU, Casey BJ, Dubin MJ, Liston C. Resting-state connectivity biomarkers define neurophysiological subtypes of depression. Nat Med. 2017 23(1):28-38.
以上のように、機能的磁気共鳴画像法などの脳機能画像法による脳活動の解析の神経・精神疾患の治療への応用を考えた場合、たとえば、上述したようなバイオマーカーとして、脳機能画像法による脳活動の解析は、非侵襲的な機能マーカーとして、診断法の開発、根本治療を実現するための創薬に向けた標的分子の探索・同定などへの応用も期待される。
たとえば、これまで、うつ病に対しては、遺伝子を用いた実用的なバイオマーカーの完成には至っておらず、それゆえに、客観的な薬物の効果判定を行うことが容易ではなく、治療薬の開発も困難であった。
脳活動の計測データに基づいて、機械学習により分類器を生成し、これをバイオマーカーとして実用化させるためには、一つの施設において計測された脳活動に対して機械学習により生成されたバイオマーカーの予測精度の向上が必要である。また、このようにして生成されたバイオマーカーが、他の施設において計測された脳活動に対しても、汎化可能であることが必要である。
すなわち、脳活動の計測データに基づく、機械学習により分類器を構築する際には、主な2つの課題が存在する。
第1の課題は、小標本サイズの問題である。
参加者の人数であるデータ量Nが、計測された脳活動計測データの次元Mよりはるかに小さいので、分類器のパラメータは、容易に訓練データに過剰適合(「過適応」または「過学習」と呼ばれる)することになる。
この過剰適合のために、構築された分類器は、新しくサンプリングされたテストデータについては、非常に貧弱なパフォーマンスを示すことになる。これは、これらのテストデータは分類器を訓練するのに使用されていないものだからである。
したがって、分類器の望ましい汎化に対して本質的な特徴量だけを識別し利用するために正則化を適切に導入する必要がある。
たとえば、上述した非特許文献5においては、必ずしも、このような正則化について十分な考慮がされているとはいえない。
第2の課題は、構築された分類器が、訓練データが集められたサイトとは異なる撮像サイトで走査されたMRIデータに対しても、よいパフォーマンスを維持する場合のみ、分類器が臨床的に有用で、科学的に信頼できるということである。
これは撮像サイトにわたる、いわゆる汎化能力である。
しかしながら、臨床応用では、異なるサイトで走査されたデータに、特定のサイトで取得したデータを使用して訓練された分類器を汎化することができないことはしばしば観測されている。
この発明は、上記のような問題点を解決するためになされたものであって、その1つの目的は、脳の活動状態に対して、うつ症状の客観的な疾患ラベルの判別を実現することである。また、この発明の他の目的は、うつ症状に対する治療効果の程度を表す情報の判別を実現することである。この発明の他の目的は、うつ症状の状態を客観的に判別するためのバイオマーカーとしての指標値を出力する分類器を提供することである。この発明の他の目的は、うつ症状のレベルの判別を実現することである。この発明の他の目的は、うつ病患者の層別化を実現することである。
本発明のある実施形態は、被検者におけるうつ症状の有無の判定を補助するための判別装置に関する。本実施形態の判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備る。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。本実施形態の判別装置は、演算装置をさらに備え、前記演算装置は、前記分類器により、被検者のうつ症状に対する分類結果を生成する判別処理を実行する。
本発明のある実施形態は、うつ症状のレベルの判別装置に関する。本実施形態の判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備える。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成する。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。本実施形態の判別装置は、演算装置をさらに備え、前記演算装置は、前記分類器により、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する。前記演算装置は、前記指標値をあらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較する。前記演算装置は、前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する。
本発明のある実施形態は、被検者における治療効果の判定するための判別装置に関する。本実施形態の判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備える。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。本実施形態の判別装置は、演算装置をさらに備え、前記演算装置は、前記分類器により、第1の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第1の値を生成する。前記演算装置は、前記分類器により、治療開始後であって第1の時点より後の第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素についてうつ症状を評価する第2の値を生成する。前記演算装置は、前記第1の値と前記第2の値とを比較する。本実施形態の判別装置は、前記第2の値が前記第1の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定する、及び/又は前記第2の値が前記第1の値よりも改善していない場合に、前記治療が被検者のうつ症状の改善に有効でないと決定する。
本発明のある実施形態は、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、うつ病患者の層別化するための判別装置に関する。本実施形態の判定装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備える。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。本実施形態の判定装置は、演算装置をさらに備え、前記演算装置は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する。前記演算装置は、前記指標値を、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較する。前記演算装置は、前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する。
本発明のある実施形態は、前記判別装置に上記処理を実行させるためのコンピュータプログラムに関する。
本発明のある実施形態は、うつ症状を有する被検者を判定することを補助する判別方法に関する。本実施形態の判別方法は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を含む。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含む。本実施形態の判別方法は、前記指標値が基準値を超えた場合に、前記被検者がうつ症状を有すると決定する工程を含む。
本発明のある実施形態は、被検者においてうつ症状のレベルを判定することを補助するため判別方法に関する。本実施形態の判別方法は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を備える。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含む。本実施形態の判別方法は、前記指標値を、あらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較する工程と、前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する工程とを含む。
本発明のある実施の形態は、被検者における治療効果の判定を補助するための判別方法に関する。本実施形態の判別方法は、第1の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第1の値を生成する工程を備える。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含む。また、本実施形態の判別方法は、治療開始後であって第1の時点より後の第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素ついてうつ症状を評価する第2の値を生成する工程を含む。本実施形態の判別方法は、前記第1の値と前記第2の値とを比較する工程と、前記第2の値が前記第1の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定する工程と、を含む。
本発明のある実施の形態は、被検者における治療効果の判定を補助するための判別方法に関する。本実施形態の判別方法は、第1の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第1の値を生成する工程を備える。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含む。本実施形態の判別方法は、治療開始後であって第1の時点より後の第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素ついてうつ症状を評価する第2の値を生成する工程を含む。本実施形態の判別方法は、前記第1の値と前記第2の値とを比較する工程と、前記第2の値が前記第1の値よりも改善していない場合に、前記治療が被検者のうつ症状の改善に有効でないと決定する工程と、を含む。
本発明のある実施の形態は、うつ病患者の層別化を補助するための判別方法に関する。本実施形態の判別方法は、うつ病が予め設定された複数のサブクラスに層別化されている場合に、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を備える。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含む。本実施形態の判別方法は、前記指標値を、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較する工程と、前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する工程と、を含む。
本発明のある実施の形態は、うつ病患者の層別化を補助するための判別装置に関する。本実施形態の判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置と、層別化処理を実行する演算装置とを備える。本実施形態の判別装置において、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ病のサブクラスのラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記選択された機能的結合に対応する判別対象の相関行列の要素に対して前記うつ病のサブクラスのラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。前記演算装置は、前記重み付け和と前記判別対象の相関行列の要素から前記サブクラスを判別するための処理を実行する。。
本発明のある実施の形態は、うつ病患者における治療効果の判定を補助するための判別装置に関する。本実施形態の判別装置は、演算装置を備え、前記演算装置は、複数の対象者について、下記表1に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する分類器生成処理を実行する。前記演算装置は、第1の時点の被検者の安静時の脳内における、前記複数の機能的結合の第1の相関を計測し、第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第2相関を計測する。前記演算装置は、前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する。
本発明のある実施の形態は、うつ病患者における治療効果の判定を補助するための判別方法に関する。本実施形態の判別方法は、複数の対象者について、下記表1に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する工程を備える。本実施形態の判別方法は、第1の時点の被検者の安静時の脳内における、前記複数の機能的結合の第1の相関を計測する工程と、第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第2相関を計測する工程を備える。本実施形態の判別方法は、前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する工程と、を含む。
本発明のある実施の形態は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、被検者におけるうつ症状の有無の判定またはうつ症状のレベルの判定を補助するために使用する方法に関する。前記使用方法は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別する前記分類器を生成する工程を備える。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程をさらに備える。
本発明のある実施の形態は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、うつ症状の治療の効果の判定を補助するために使用する方法に関する。前記使用方法は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別する前記分類器を生成する工程を備える。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程をさらに備える。
本発明のある実施の形態は、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、安
静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、被検者を前記サブクラスに層別化することを補助するために使用する方法に関する。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程を含む。
本発明のある実施の形態は、 被検者における治療効果の判定を補助するための判別装置に関する。前記判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から第1の分類器を生成する第1の演算装置を含む分類器生成装置を備る。前記第1の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。前記第1の演算装置は、複数の対象者について、下記表5に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第2の分類器を生成する。前記判別装置は、第2の演算装置と記憶装置を含む層別化装置をさらに備る。前記記憶装置は、前記第1および第2の分類器の情報と、うつ症状が予め設定された複数のサブクラスに層別化するための情報とを記憶する。前記第2の演算装置は、前記第1の分類器を用いて、被検者を前記サブクラスに層別化処理を実行する。前記第2の演算装置は、前記層別化処理によって特定のサブクラスに層別化された被検者について、第1の時点の安静時の前記複数の機能的結合の第1の相関を計測するための処理を実行する。前記第2の演算装置は、第2の時点の前記被検者と同一の被検者の安静時の前記複数の機能的結合の第2相関を計測するための処理を実行する。前記第2の演算装置は、前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記第2の分類器により、前記被検者に対する治療効果を判別する。
本発明のある実施の形態は、被検者における治療効果の判定を補助するための判別方法に関する。本実施形態の判別方法は、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から第1の分類器生成処理によって生成された第1の分類器を使用して、被検者を前記サブクラスに層別化する層別化工程を含む。本実施形態の判別方法は、前記第1の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。本実施形態の判別方法は、前記層別化工程によって特定のサブクラスに層別化された被検者について、投薬前の安静時の前記複数の機能的結合の第1の相関を計測する第1の相関計測工程と、投薬開始から所定期間経過後の前記被検者と同一の被検者の安静時の前記複数の機能的結合の第2相関を計測する第2の相関計測工程とを含む。本実施形態の判別方法は、特定のサブクラスに分類された前記被検者に対する薬効を第2の分類器を用いて判別する判別工程であって、前記第2の分類器は、複数の対象者について、下記表1に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と投薬後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測する。そして、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第2の分類器生成処理によって予め生成され前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記被検者に対する薬効を判別する、判別工程と、を備える。
本発明のある実施の形態は、演算装置と、記憶装置とを備えた第1の分類器生成装置に関する。前記演算装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器を特定する情報を生成する。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む。前記記憶装置は、前記演算装置によって生成された前記第1の分類器を特定する情報を記憶する。
本発明のある実施の形態は、第1の分類器の生成方法に関する。本実施形態では、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により第1の分類器を生成する。前記第1の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含む。
本発明のある実施の形態は、演算装置と、記憶装置とを備えた第2の分類器生成装置に関する。前記演算装置は、複数の対象者について、下記表3に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測する。本実施形態は、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第2の分類器を生成する。前記記憶装置は、前記演算装置によって生成された前記第2の分類器を特定するための情報を記憶する。
本発明のある実施の形態は、第2の分類器生成方法に関する。本実施形態では、複数の対象者について、下記表4に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測する。本実施形態は、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する工程を備える。
本発明のある実施の形態は、ニューロフィードバック訓練を実行するための脳活動訓練装置に関する。前記脳活動訓練装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を時系列で予め測定した信号に基づいて、前記複数の所定領域間の機能的結合から機械学習により、うつ症状の疾患ラベルを判別するために特徴選択により選択された複数の機能的結合のうち訓練対象の機能的結合を特定する情報を格納する記憶装置を備える。前記訓練対象の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合を含む。前記記憶装置は、ニューロフィードバック訓練における前記訓練対象の機能的結合の目標パターンを記憶する。前記脳活動訓練装置は、前記ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知するための脳活動検知装置をさらに備える。前記被訓練者の脳内の複数の所定領域のそれぞれは、前記複数の参加者における脳内の複数の所定領域にそれぞれ対応する。前記脳活動訓練装置は、提示装置と、演算装置とをさらに備える。前記演算装置は、i)前記脳活動検知装置により検知された信号から、前記訓練対象の機能的結合の時間相関を所定期間について算出し、ii)算出された前記時間相関に基づいて、前記目標パターンとの近似度に応じて、報酬値を算出し、iii)前記報酬値の大きさを示す情報を前記提示装置により前記被訓練者に対して提示する。
本発明のある実施の形態は、ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知するための脳活動検知装置と、演算装置と、記憶装置と、提示装置と、を備える脳活動訓練装置の制御方法に関する。前記制御方法は、前記記憶装置に、i)安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を時系列で予め測定した信号に基づいて、前記複数の所定領域間の機能的結合から機械学習により、うつ症状の疾患ラベルを判別するために特徴選択により選択された複数の機能的結合のうち訓練対象の機能的結合を特定する情報と、ii)ニューロフィードバック訓練における前記訓練対象の機能的結合の目標パターンとを格納する工程を備える。前記訓練対象の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合を含む。前記制御方法は、前記脳活動検知装置により、前記ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知する工程をさらに備える。前記被訓練者の脳内の複数の所定領域のそれぞれは、前記複数の参加者における脳内の複数の所定領域にそれぞれ対応する。前記制御方法は、前記演算装置が、前記脳活動検知装置により検知された信号から、前記訓練対象の機能的結合の時間相関を所定期間について算出する工程と、前記演算装置が、算出された前記時間相関に基づいて、前記目標パターンとの近似度に応じて、報酬値を算出する工程と、前記演算装置が、前記報酬値の大きさを示す情報を前記提示装置により前記被訓練者に対して提示する工程と、をさらに備える。
この発明によれば、脳の活動状態に対してうつ症状の状態に対する指標を客観的に判別する判別装置及び判別方法を提供することができる。この発明によれば、うつ症状に対する治療効果の程度を表す情報を客観的に判別することができる。この発明によれば、うつ症状を客観的に判別するためのバイオマーカーとしての指標値を出力する分類器を提供することができる。この発明によれば、うつ症状のレベルを判別することができる。この発明によれば、うつ病患者を客観的に層別化することができる。
この発明によれば、ニューロフィードバック訓練を行うことができる。
MRI装置10の全体構成を示す模式図である。 データ処理部32のハードウェアブロック図である。 関心領域について、安静時の機能的結合の相関を示す相関行列を抽出する手続きを示す概念図である。 相関行列から、第1の分類器を生成する過程を説明する概念図である。 図4で説明したような第1の分類器の生成処理および生成された第1の分類器による判別処理を行うための機能ブロック図である。 第1の分類器を生成するためにデータ処理部32が行う処理を説明するためのフローチャートである。 内ループ特徴選択を説明するための概念図である。 内ループ特徴選択処理の概念を説明するための図である。 一例として、内ループ特徴選択を行った際の反復処理のうち、特定のハイパーパラメータλ1およびλ2に対する結果を示す図である。 内ループ特徴選択の処理を、より詳しく説明するためのフローチャートである。 外ループ特徴選択を説明するための概念図である。 外ループ特徴選択の処理を、より詳しく説明するためのフローチャートである。 ステップS108での判別器の生成処理の概念を示す図である。 被検者についてうつ症状の判別を行うためのフローチャートである。 データ収集、推定処理および対象者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。 被検者についてうつ症状のレベルの判別を行うためのフローチャートである。 被検者について治療効果を判別するためのフローチャートである。 被検者について治療効果を判別するためのフローチャートである。 うつ病のサブクラスを判別するためのフローチャートである。 分類器の生成処理および生成された分類器による層別化処理を行うための機能ブロック図である。 図19に示す、データ収集、推定処理および被検者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。 第2の分類器を生成するためのフローチャートである。 第2の分類器を使用して薬効を判定するためのフローチャートである。 第2の分類器の生成処理および生成された第2の分類器による治療効果判別処理を行うための機能ブロック図である。 うつ病のサブクラスを決定した被検者に対して第2の分類器を用いて薬効を判定するためのフローチャートである。 選択された12対の機能的結合を示す。 54個の非ゼロのci kの大きさの分布を示す図である。 aは、広島のコホートの関連重み付き和の分布を示す。cは、千葉のコホートの関連重み付き和の分布を示す。cは、訓練データとテストデータを入れ替えて判別成績を検討した結果を示す。†p < 0.10*p < 0.05, **p < 0.01。 dは、メランコリー型MDD群の関連重み付き和の平滑化ヒストグラムを示す。eは、非メランコリー型MDD群の関連重み付き和の平滑化ヒストグラムを示す。†p < 0.10*p < 0.05, **p < 0.01。 fは、治療抵抗性MDD群の関連重み付き和の平滑化ヒストグラムを示す。gは、ASD群の関連重み付き和の平滑化ヒストグラムを示す。hは、SSD群の関連重み付き和の平滑化ヒストグラムを示す。†p < 0.10*p < 0.05, **p < 0.01。 並び替え検定の結果を示す。aは、訓練データのLOOCVの並び替え検定(1000回繰り返し)のヒストグラムを示す。bは、サンプル外のテストデータセットの精度を示し、二項分布を曲線として示す。並び替え無しで訓練され、検証されたメランコリックMDD分類器の精度は、垂直線で示す。並び替え検定の結果は、LOOCV(p=0.049)およびサンプル外検証(p=0.036)について有意であった。 * p <0.05。 aは、BDIスコアと第1の分類器を用いた判別結果の相関(全MDDと健常者全体)を示す。bは、BDIと第1の分類器を用いた判別結果の相関(全MDD全体)を示す。cは、治療前後の関連重み付き和の平滑化ヒストグラムを示す。dは、投薬後のBDIスコアの変化と第1の分類器を用いた判別結果の変化の相関を示す。 aは、各機能的結合における健常者対照群とMDD群の相関の差と、投薬前と投薬後の相関の差を示す。bは、広島のコホートについての、FC1とFC2における投薬前と投薬後の相関の差を示す。cは、千葉のコホートについての、FC1とFC2における投薬前と投薬後の相関の差を示す。dは、健常者のコホートについての、FC1とFC2における投薬前と投薬後の相関の差を示す。 eは、投薬後寛解群と投薬後非寛解群のFC1とFC2における投薬前と投薬後の相関の差を示す。fは、符号付ΔFC1とΔFC2の投薬群での分布を示す。 脳活動訓練装置の構成の概念を示す図である。 指標値と目標値の近さを表すためのモニター表示の一例である。 ニューロフィードバックにおける訓練シーケンスの一例を示す図である。 aは、MDD参加者の訓練期間中のニューロフィードバックスコアを示す。bは、MDD参加者の訓練前後におけるHDRSスコアを示す。cは、亜臨床的うつ病の参加者の訓練期間中のニューロフィードバックスコアを示す。dは、訓練後のBDIスコアとrs−fcMRI変化量の相関を示す。
以下、本発明の実施の形態について、図に従って説明する。なお、以下の実施の形態において、同じ符号を付した構成要素および処理工程は、同一または相当するものであり、必要でない場合は、その説明は繰り返さない。
1.用語の説明
本発明において、うつ症状には、抑うつ気分;興味や関心の低下;気力の低下;焦燥;制止;思考力、集中力、又は決断力の減退;無価値感、又は罪責感;自殺、自殺念慮、又は自殺企画;病的思考内容;妄想;身体症状(全身倦怠感、頭痛、頭重、又は腰痛等の身体諸所の痛み、動悸、息切れ、食欲減退、又は体重減少等)、睡眠障害から選択される少なくとも一種が含まれる。好ましくは、前記うつ症状には、精神障害の診断と統計マニュアル(DSM)−IVの基準の大うつ病性障害(MDD)に伴う症状が含まれる。
本発明において、うつ病は、前記うつ症状を伴う限り制限されないが、好ましくはMDDである。MDDには、メランコリー型MDD、非メランコリー型MDD及び治療抵抗性MDDが含まれる。本明細書においては、MDDを単にうつ病と呼ぶことがある。
抑うつの程度の評価には、たとえば、従来から、自己申告によるBeck Depression Inventory (ベックうつ病調査表)による評価が、診察時のスクリーニングや面接時の補助資料として利用されている。また、医師用の評価尺度として使用されるハミルトンうつ病評価尺度(HDRS:Hamilton Depression Rating Scale)は、うつ病の指標を提供する多重項目質問票であり、回復を評価する指標となる。略称は、HAM-Dである。
本発明において、被検者は、制限されないが、前記うつ症状を有する者であることが好ましい。年齢及び性別は制限されない。前記被検者は、前記うつ症状を改善するための治療を受けていない者であっても、前記治療を受けている者であってもよい。
うつ症状を改善するための治療には、投薬療法、ニューロフィードバック療法、修正型電気けいれん療法および反復経頭蓋磁気刺激療法から選択される少なくとも一種が含まれる。投薬される薬剤としては、三環系抗うつ薬(イミプラミン、トリミプラミン、クロミプラミン、アミトリプチリン、ノルトリプチリン、アモキサピン、ロフェプラミン、及びドスレピン等)、四環系抗うつ薬(マプロチリン、アミンセリン、及びセチプチリン等)、トラゾドン、選択的セロトニン再取り込み阻害薬(エスシタロプラム、フルボキサミン、パロキセチン、及びセルトリラリン等)、セロトニン−ノルエピネフリン再取り込み阻害剤(ミルナシプラン、及びデュロキセチン等)、ノルアドレナリン・セロトニン作動性抗うつ剤(ミルタザピン等)、及びベンザミド系薬剤(スルピリド等)から選択される少なくとも一種を挙げることができる。
うつ症状の治療には、「治癒」と言う概念がないため、うつ症状が改善するとは、BDI等の臨床所見から見て、以前の状態よりもうつ症状がよくなっている状態、または臨床所見から見て「寛解(remitted)」といえる状態になることを意味する。
本発明の実施の形態には、被検者においてうつ症状のレベルを判定することを補助すること、被検者におけるうつ症状の有無の判定を補助すること、被検者における治療効果の判定を補助すること、及び被検者を疾患のサブクラスに層別化することを補助することが含まれる。
本発明において、ドラッグリプロファイリングとは、既にヒトに対して別の疾患または症状に対する薬効が明らかとなっている薬剤、またはヒトに対して毒性も弱いが有意な薬効も認められないことが明らかになっている薬剤であって、うつ症状に対する影響が未知である薬剤について、うつ症状に対する薬効を検出することを意図する。
2.安静時fMRIの撮像条件と相関行列の抽出
安静時fMRIの撮像装置は、制限されない。撮像条件もfMRI画像を取得できる限り制限されない。例えば、magnetic filedは3.0T程度、field of viewは192mm〜256mm程度、matrixは64×64程度、number of slicesは30枚〜40枚程度、number of volumesは112〜244程度、slice thicknessは3.0mm〜4.0mm程度、slice gapは0mm〜0.8mm程度、TRは2,000ms〜2,700ms程度、TEは25ms〜31ms程度、総撮像時間は5分〜10分程度、Flip angleは、75deg〜90deg程度、slice acquisition orderは、Ascending (Interleaved)である。撮像は、暗い照明下で行うことが好ましい。また被検者には、撮像中何も考えず、眠らずにいてもらうことが好ましい。さらに、被検者には、撮像中モニター画面中央の十字マークを見続けてもらうことが好ましい。
撮像されたfMRI画像は、Yahataらの文献(NATURE COMMUNICATIONS | 7:11254 | DOI: 10.1038/ncomms11254)に記載の方法で処理することができる。
撮像されたfMRI画像データについて、特に限定されないが、たとえば、Matlab R2014a(Mathworks inc., USA)のSPM8(Wellcome Trust Center for Neuroimaging、University College London、UK)を用いてT1強調構造画像及び安静状態機能画像における標準的な前処理を施すことができる。例えば、撮像画像に対して、Realign、Slice−timing correction、Coregistering、Normalization及びSmoothing(FWHM=8mm)を行うことができる。さらに、各被検者の全画像データについて,直前の画像から0.5mm 以上動いたと判断される画像を解析から除外してもよい。
前処理を行った画像データから、脳内の所定の複数の関心領域(Region of interest, ROI)間について、機能的結合(Functional connectivity, FC)(すなわち、結合強度)を算出する。機能的結合は、安静時脳活動解析において一般的に利用される特徴量(相関行列の要素)であり、異なる関心領域の時系列信号間のピアソンの相関係数により定義される。特に限定されないが、たとえば、所定の時間にわたって計測された値からのピアソンの相関係数の平均値により、結合強度を表すことができる。ただし、所定の時間内のピアソンの相関係数についての他の統計量により結合強度を表してもよい。
なお、後述するように、本実施の形態では、うつ症状の状態を判別するためには、たとえば、下記表1に示す機能的結合識別番号(表1でID)1〜12の関心領域間について、その全て、または、一部について機能的結合(結合強度)を算出する。
表1で、Lat.の“L”と“R”は、左脳と右脳の区別を示す。BSAは、ブロードマン領野を示し、BAは、ブロードマンの領野の番号を示す。Weightは、後述する関連重み付け和(単に「重み付け和」と称することもある)の重みを示す。
相関行列の要素の抽出は、特に限定されないが、たとえば、以下のような手順で実行することができる。
はじめに各関心領域に含まれる全ボクセル(Voxel)の平均信号の時系列を抽出する。次に信号値からのノイズ除去のため、バンドパスフィルタ(0.008〜0.1Hz)をかける。その後、9つの説明変数(全脳、白質、脳脊髄液の平均信号及び6つの体動補正パラメータ)による回帰を行う。回帰後の残差系列を機能的結合と関連する時系列信号値と考え、前記時系列信号値を各関心領域の要素とし、関心領域の各対の関心領域の各相関行列の要素について、時系列間でピアソンの相関係数を計算する。前記相関係数が機能的結合の結合強度を表す値であり、各対それぞれに対応する結合強度が求められる。
後述する分類器に入力データとして前記各機能的結合の相関係数を入力することにより、前記相関係数と、分類器の分類処理を実行するための各機能的結合について予め算出された重み(寄与度)を示す係数とに基づいて、機能的結合についてうつ症状の疾患ラベルを判別するための指標値が算出される。あるいは前記相関係数と前記係数とに基づいて、複数の機能的結合についての関連重み付け和がうつ症状の疾患ラベルを判別するための指標値として算出される。ここで、「関連重み付け和」とは、複数の機能的結合にそれぞれ対応する重みを乗算して和をとったものをいう。
このため、前記指標値は、単に被検者の脳活動を測定したデータではなく、各機能的結合の重みを考慮して人為的に算出された値である。前記指標値は、うつ症状のラベルの判別のために、うつ症状のレベルの判別のために、治療効果の判別のために、または、うつ病患者の層別化のために使用される。
前記12対の機能的結合のうち、最もうつ症状への寄与度が高いのは、機能的結合識別番号1で示される機能的結合である。次に寄与度が高いのは、機能的結合識別番号2で示される機能的結合である。したがって、以下に説明する各実施の形態においては、少なくとも機能的結合識別番号1及び機能的結合識別番号2の一方、あるいは両方を選択して使用してもよい。
前記12対の機能的結合は、特にメランコリー型MDD患者群と健常対照群を判別することに適している。
3.第1の分類器の生成と判別装置
第1の分類器は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理によって生成される。ここで、うつ病患者とは、予め医師の診断により、うつ病と診断され、うつ病の「疾患ラベル」と関連づけられた参加者を意味する。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された機能的結合の重みに基づいて、前記うつ症状の疾患ラベルを判別するように生成される。
すなわち、後述するように、予め設定された複数の関心領域に対する機能的結合のうち、うつ症状の疾患ラベルに特異的に関連するとしてスパース正準相関分析により抽出された機能的結合から、スパースロジスティック回帰による特徴選択により、うつ病の疾患ラベルの判別に使用される機能的結合が選択される。そして、このようにして選択された機能的結合から、上述した関連重み付け和が算出される。
3−1.判別装置1
例えば、第1の分類器は、図1に示すMRI装置10から取得されるfMRIの画像データから生成される。
図1は、本発明の実施の形態1にかかる、判別装置1であるMRI装置10の全体構成を示す模式図である。
図1に示すように、MRI装置10は、MRIを撮像するMRI撮像部25とMRI撮像部25の制御シーケンスを設定するとともに各種データ信号を処理して画像を生成するデータ処理部32とを備える。MRI撮像部25は、被検者2の関心領域に制御された磁場を付与してRF波を照射する磁場印加機構11と、この被検者(又は被訓練者)2からの応答波(NMR信号)を受信してアナログ信号を出力する受信コイル20と、この被検者2に付与される磁場を制御するとともにRF波の送受信を制御する駆動部21とを備える。
なお、ここで、被検者2が載置される円筒形状のボアの中心軸をZ軸にとりZ軸と直交する水平方向にX軸及び鉛直方向にY軸を定義する。
MRI装置10は、このような構成であるので、磁場印加機構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は、制御シーケンスに従って信号送信部33から送信される高周波信号に基づいて、被検者2の関心領域にRF(Radio Frequency)パルスを照射するものである。
なお、RF照射部16は、図1において、磁場印加機構11に内蔵されているが、寝台18に設けられたり、あるいは、受信コイル20と一体化されていてもよい。
受信コイル20は、被検者2からの応答波(NMR信号)を検出するものであって、このNMR信号を高感度で検出するために、被検者2に近接して配置されている。
ここで、受信コイル20には、NMR信号の電磁波がそのコイル素線を切ると電磁誘導に基づき微弱電流が生じる。この微弱電流は、信号受信部28において増幅され、さらにアナログ信号からデジタル信号に変換されデータ処理部32に送られる。
すなわち、静磁界にZ軸傾斜磁界を加えた状態にある被検者2に、共鳴周波数の高周波電磁界を、RF照射部16を通じて印加すると、磁界の強さが共鳴条件になっている部分の所定の原子核、たとえば、水素原子核が、選択的に励起されて共鳴し始める。共鳴条件に合致した部分(たとえば、被検者2の所定の厚さの断層)にある所定の原子核が励起され、スピンがいっせいに回転する。励起パルスを止めると、受信コイル20には、今度は、回転しているスピンが放射する電磁波が信号を誘起し、しばらくの間、この信号が検出される。この信号によって、被検者2の体内の、所定の原子を含んだ組織を観察する。そして、信号の発信位置を知るために、XとYの傾斜磁界を加えて信号を検知する、という構成になっている。
画像処理部48は、記憶部36に構築されているデータに基づき、励起信号を繰り返し与えつつ検出信号を測定し、1回目のフーリエ変換計算により、共鳴の周波数をX座標に還元し、2回目のフーリエ変換でY座標を復元して画像を得て、表示部38に対応する画像を表示する。
たとえば、このようなMRIシステムにより、上述したBOLD信号をリアルタイムで撮像し、制御部42により、時系列に撮像される画像について、後に説明するような解析処理を行うことで、安静時機能的結合MRI(rs−fcMRI)の撮像を行うことが可能となる。
図2は、データ処理部32のハードウェアブロック図である。
データ処理部32のハードウェアとしては、上述のとおり、特に限定されないが、汎用コンピュータを使用することが可能である。
図2において、データ処理部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に相当する。データ処理部32のコンピュータ本体2010は、第1の分類器生成装置、及び/又は第2の分類器生成装置としても機能する。
CPU2040が、プログラムに基づいて実行する演算処理により、データ処理部32の各機能、たとえば、制御部42、データ収集部46、画像処理部48の各機能が実現される。
データ処理部32に、上述した実施の形態の機能を実行させるプログラムは、CD−ROM2200、またはメモリ媒体2210に記憶されて、ディスクドライブ2030またはメモリドライブ2020に挿入され、さらに不揮発性記憶装置2080に転送されても良い。プログラムは実行の際にRAM2070にロードされる。
データ処理部32は、さらに、入力装置としてのキーボード2100およびマウス2110と、出力装置としてのディスプレイ2120とを備える。キーボード2100およびマウス2110が入力部40に相当し、ディスプレイ2120が表示部38に相当する。
上述したようなデータ処理部32として機能するためのプログラムは、コンピュータ本体2010に、情報処理装置等の機能を実行させるオペレーティングシステム(OS)は、必ずしも含まなくても良い。プログラムは、制御された態様で適切な機能(モジュール)を呼び出し、所望の結果が得られるようにする命令の部分のみを含んでいれば良い。データ処理部32がどのように動作するかは周知であり、詳細な説明は省略する。
また、上記プログラムを実行するコンピュータは、単数であってもよく、複数であってもよい。すなわち、集中処理を行ってもよく、あるいは分散処理を行ってもよい。
図3は、うつ症状の疾患ラベルに関連する安静時の機能的結合の相関を示す相関行列を抽出する手続きを示す概念図である。
図3に示すように、リアルタイムで測定した安静時のfMRIのn個(n:自然数)の時刻分のfMRIデータから、各関心領域の平均的な「活動度」を算出し、脳領域間(関心領域間)の活動度の相関値を算出する。
ここで、関心領域としては、小脳領域を除く137領域を考えているので、相関行列における独立な非対角成分は、対称性を考慮すると、
(137×137−137)/2=9316(個)
ということになる。なお、図3では、図の上では、相関は34×34のみを示している。
このような相関行列の要素の算出については、特に限定されないが、たとえば、以下のように実行することが可能である。
安静時脳活動データから、各参加者について、異なる関心領域間の機能的結合(Functional connectivity, FC)を計算する。FCは安静時脳活動解析において一般的に利用される特徴量であり、異なる関心領域の時系列信号間のピアソンの相関係数により定義される。
まず、各関心領域に含まれる全ボクセルの平均信号の時系列を抽出する。
次に信号値からのノイズ除去のため、バンドパスフィルタをかけ、その後、9つの説明変数(全脳,白質,脳脊,髄液の平均信号及び6つの体動補正パラメータ)による回帰を行う。
回帰後の残差系列を機能的結合と関連する時系列信号値と考え、異なるROIの時系列間でピアソンの相関係数を計算する。
関心領域については、Brain Sulci Atlas (BAL)に含まれる137の関心領域を用いる。これら137の関心領域間の機能的結合FCを特徴量として用いる。
ここで、Brain Sulci Atlas (BAL)については、以下に開示がある。
文献:Perrot et al., Med Image Anal, 15(4), 2011
文献:Tzourio-Mazoyer et al., Neuroimage, 15(1), 2002
図4は、図3で説明したような相関行列から、バイオマーカーとなる第1の分類器を生成する過程を説明する概念図である。
以下、説明するように、第1の分類器は、たとえば、「うつ病」に対する診断マーカーとしての指標として機能し、疾患の重症度や進行度を医師が診断する際の支援のために使用することができる。したがって、第1の分類器の演算と結果の出力を実行する装置は、診断支援装置と呼ぶこともできる。
図4に示すように、MRI撮像部25で撮像された、健常群、及びうつ病患者群を含む参加者において測定されたfMRIのデータから、それぞれの参加者について後に説明するような手続により、データ処理部32が、脳領域間(関心領域間)の活動度の相関行列の要素を導出する。
続いて、データ処理部32により、参加者の疾患/健常ラベルを含む参加者の属性と相関行列とについて正則化正準相関解析により特徴抽出が行われる。「正則化」とは、一般に、機械学習や統計学において、誤差関数に、ハイパーパラメータにより重み付けされた正則化項を追加して、モデルの複雑度・自由度に抑制を加え、過学習を防ぐ方法である。なお、正則化正準相関解析の結果、説明変数についてスパース化も併せて実現される場合は、特に、スパース正準相関解析(Sparse Canonical Correlation Analysis : SCCA)と呼ぶ。以下では、具体例としては、スパース正準相関解析を行うものとして説明する。
そして、このようなスパース正準相関解析において、後述するように、ハイパーパラメータの値の調節により「疾患ラベル」とのみ結合する正準変数が存在するようにしたときに、これに対応する正準変数に接続する機能的結合FCを抽出する。所定範囲でハイパーパラメータを変化させるとき、このような条件を満たす正準変数が存在する範囲で抽出される機能的結合FCの和集合を「第1和集合」と呼ぶ。
さらに、データ処理部32により、スパース正則化正準相関解析の結果得られた「第1和集合」を説明変数とし、たとえば、1個抜き交差検証(Leave-One-Out Cross Validation:LOOCV)を行いながら、各交差検証のステップでスパースロジスティック回帰による判別分析を行い、すべての交差検証にわたって、説明変数として抽出された機能的結合FCの和集合を「第2和集合」と呼ぶ。
最後に、全ての参加者のデータに対して、「第2和集合」を説明変数として、目的変数である「疾患ラベル」に対して、スパースロジスティック回帰による判別分析を行って、第1の分類器が生成される。
図5は、図4で説明したような第1の分類器の生成処理および生成された第1の分類器による判別処理を行うための機能ブロック図である。
まず、不揮発性記憶装置2080には、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号をMRI装置により時系列で予め測定した信号の情報である参加者rs−fcMRI測定データ3102と、MRI測定データを計測した各参加者に関連づけられた、複数の人属性情報3104とが格納されている。
ここで、「人属性情報」とは、参加者個人を特定するための「人特性情報」と、各参加者に対する測定条件を特定するための「測定条件情報」とを含む。
「人特性情報」とは、参加者に対する疾患ラベル、年齢、性別、投薬履歴などの情報を意味する。
「測定条件情報」とは、参加者を測定した測定サイト情報(測定した施設および/または測定した装置を特定する情報を含む)や、測定中に、開眼/閉眼のいずれで計測したか、測定磁場強度などの測定条件を特定する条件である。
演算装置2040は、参加者rs−fcMRI測定データ3102および対応する人属性情報3104とに基づいて、疾患ラベルに対する分類器を生成する処理を行う。
相関行列算出部3002は、参加者rs−fcMRI測定データ3102から複数の所定領域間の脳活動の機能的結合の相関行列を参加者ごとに算出する。算出された機能的結合の相関行列のデータは、参加者ごとに、不揮発性記憶装置2080に、機能的結合の相関行列のデータ3106として格納される。
第1特徴選択部3004は、複数の参加者から抽出されたK個(K:2以上の自然数)の異なる部分集合に対して、順次1つの部分集合を選択し、選択された部分集合を除く(K−1)個の部分集合に対して、複数の属性情報と相関行例の要素とに対するスパース正準相関分析により、複数の人属性情報のうちの特定の属性情報、たとえば、疾患ラベルのみに対応する正準変数と接続する相関行例の要素を抽出する。さらに、第1特徴選択部3004は、順次選択される部分集合にわたって、抽出された相関行例の要素の和集合である第1和集合を取得し、不揮発性記憶装置2080に、第1の機能的結合和集合データ3108として格納する。なお、「第1の機能的結合和集合データ」は、機能的結合の相関行列のデータ3106のうち、第1の和集合に対応する要素を特定するためのインデックスであってもよい。
第2特徴選択部3006は、複数の参加者のうち上記K個の部分集合を除く残りの参加者をテスト集合とし、テスト集合を異なるN個のグループに分けるとき、複数の参加者からN個のグループうちの選択された1つのグループを除く参加者の集合を対象として、第1和集合に基づいて、上記特定の属性情報(たとえば、疾患ラベル)を推定するためのテスト分類器をスパースロジスティック回帰により算出し、スパース化に伴いテスト分類器の説明変数となる相関行例の要素を抽出する。第2特徴選択部3006は、さらに、N個のグループうちから1つのグループを順次選択しつつ、特徴抽出を繰り返して、テスト分類器の説明変数として抽出された相関行例の要素の和集合である第2和集合を取得して、不揮発性記憶装置2080に、第2の機能的結合和集合データ3110として格納する。なお、「第2の機能的結合和集合データ」も、機能的結合の相関行列のデータ3106のうち、第2の和集合に対応する要素を特定するためのインデックスであってもよい。
なお、N個グループの各々の要素の個数は1個であってもよい。
分類器生成部3008は、第2和集合を説明変数として、特定の属性情報(たとえば、疾患ラベル)を推定するための第1の分類器をスパースロジスティック回帰により算出する。分類器生成部3008は、生成された第1の分類器を特定するための情報を、不揮発性記憶装置2080に、分類器データ3112として格納する。
判別処理部3010は、分類器データ3112により特定される分類器に基づいて、入力データに対して判別処理を行う。
なお、上記の説明では、分類器生成部3008は、第2和集合を説明変数として、分類器を生成するものとして説明したが、たとえば、分類器生成部3008は、直接、第1和集合を説明変数として、第1の分類器を生成してもよい。ただし、後に説明するように、次元の低減や、汎化性能の観点からは、第2和集合を説明変数とすることが、望ましい。
なお、第1の分類器の生成にあたっては、正則化法(たとえば、L1正則化やL2正則化)を用いたロジスティック回帰である正則化ロジスティック回帰を用いることができ、より具体的には、たとえば、上述した「スパースロジスティック回帰」を使用することもできる。また、第1の分類器の生成には、サポートベクトルマシンや、LDA(Linear discriminant analysis)などを使用してもよい。以下では、スパースロジスティック回帰を例にとって説明することにする。
後に説明するように、第2特徴選択部3006による特徴選択処理と並行して、N個のグループうちから除かれる1つのグループを順次選択しつつ、第2特徴量選択部3006により算出されたテスト分類器により、除かれたグループをテストサンプルとして判別結果を算出することにより交差検証を行ってもよい。
このように説明変数の次元の縮小を、入れ子構造となった特徴選択手続きにより行うことにより、第2特徴選択部3006の行う処理において、ほとんどの参加者のデータを使用する一方、時間的に効率的な次元の数の縮小を実行できる。
さらに、第2特徴選択部3006の処理におけるテスト集合を、次元の数を縮小するために使用されるデータセットとは、独立となるようにしておくことにより、極端に楽観的な結果が得られることを回避することができる。
図6は、バイオマーカーとなる分類器を生成するためにデータ処理部32が行う処理を説明するためのフローチャートである。
以下、図4で説明した処理を、図6を参照しながら、より詳しく説明する。
安静時のfMRIデータから導き出される脳領野間の結合と参加者の疾患の判別ラベル(疾患ラベル)にもとづいてバイオマーカーを作成する場合の最大の問題点は、データ数にくらべてデータの次元が圧倒的に多いことである。そのため、正則化を行うことなしに、疾患ラベル(ここでは、参加者が疾患であるか健常であるかを示すラベルを「疾患ラベル」と呼ぶ)の予測のための分類器の学習をデータセットを用いておこなうと、オーバーフィット(過剰適合)することで未知のデータに対する予測性能が著しく低下する。
ここで、一般に、機械学習において、観測データをより少ない数の説明変数で説明するように行われる処理のことを「特徴選択(または特徴抽出)」と呼ぶ。本実施の形態では、「脳領域間(関心領域間)の活動度の複数個の相関値(複数個の結合)」のうち、対象となる疾患ラベルの予測のための第1の分類器の機械学習において、より少ない相関値により第1の分類器を構成できるように特徴選択(特徴抽出)をすること、すなわち、説明変数として、より重要な相関係数を選択することを「縮約表現を抽出する」と呼ぶこともある。
そして、本実施の形態では、特徴抽出法として、正則化法を用いる。このように、正準相関分析において、正則化とともにスパース化を実行し、説明変数としてより重要なものを残すように処理するという意味で、このような処理を「スパース正準相関分析」と呼ぶ。より具体的には、たとえば、スパース化を併せて実現するための正則化法としては、以下に説明するような「L1正則化」という正準相関分析のパラメータの絶対値の大きさにペナルティーを課す手法を用いることが可能である。
すなわち、図6を参照して、データ処理部32は、入力部40からの処理開始の入力等により、分類器生成の処理が開始されると(S100)、記憶部36から各参加者についてのMRI測定データを読み出し(S102)、スパース正準相関分析(SCCA)により第1の特徴抽出の処理を行う(S104)。
以下では、ステップS104の処理を「内ループ特徴選択」と呼ぶ。
ここで、ステップS104での内ループ特徴選択の処理を説明するために、「スパース正準相関分析」と、「内ループ特徴選択」について、以下、それぞれ説明する。
(スパース正準相関分析)
以下では、スパース正準相関分析として、L1正則化正準相関分析について説明する。このL1正則化正準相関分析については、以下の文献に開示がある。
文献:Witten DM, Tibshirani R, and T Hastie. A penalized matrix decomposition, with applications to sparse principal components and canonical correlation analysis. Biostatistics, Vol. 10, No. 3, pp. 515-534, 2009.
まず、一般の正準相関分析(Canonical Correlation Analysis: CCA)では、データ対x1およびx2を考える際に、それぞれの変数x1およびx2は、平均0、標準偏差1を持つように標準化される。
一般的に言えば、正準相関分析(CCA)の使用によって、ペアになった観測量の間の潜在している関係を識別することができる。
具体的には、正準相関分析(CCA)では、対となる射影された変数(正準変数)に最大の相関があるように、投影ベクトルを見つけ出す。
これに対して、正準相関分析にL1正則化を導入した場合、すなわち、L1ノルム正規化によるスパースCCA、すなわち、L1−SCCAを使用する場合は、以下のような最適化問題を解くことになる。
以下のような変数の組を考える。
ここでも、行列X1およびX2を構成する列については、それぞれ平均が0で分散が1となるように正規化されているものとする。
そうすると、L1−SCCAは、以下の式(1)のように定式化することができる。
ここで、ハイパーパラメータλ1およびλ2が、重みパラメータw1およびw2の(w1およびw2により対応する変数がスパース化されることから「スパース投影ベクトル」と呼ぶ)スパース化の度合いをそれぞれ示す。
具体的に、本実施の形態の例としては、人属性情報および機能的な結合FCの潜在している関係を識別するために、上記の変数に相当する2つのデータ行列を構築するものとする。
以下では、参加者には、うつ病の患者群(疾患ラベルを「MDD(Major depressive disorder)」で表す)と健常者対照群(疾患ラベルを「HC(Healthy control)」で表す)第1のデータ行列X1の行は、1人の参加者の人属性情報(人特性情報および測定条件情報)を表しており、たとえば、この特性情報や測定条件情報は、以下のものを含む:
i)疾患ラベル(MDDあるいはHC),
ii)サイト情報(参加者の脳活動がどこで計測されたか、サイトA,B,C) iii)年齢
iv)性別
v)撮像条件(開眼あるいは閉眼)
vi)薬剤投与の状態1(抗精神病薬)
vii)薬剤投与の状態2(抗うつ薬)
viii)薬剤投与の状態3(精神安定剤)
具体的には、人属性情報データ行列X1の列数は、10であり、つまりp1=10である。
最初の列は、1(=MDD)あるいは0(=HC)を含んでいる。
次の3つの列は、脳活動の計測された施設(場所)を示しており、この場合は、3つの異なるサイトで計測が行われたことに相当して、[100](サイトA)、[010](サイトB)あるいは[001](サイトC)のいずれかを含んでいる。
第5列は、参加者の年齢の値を含んでおり、第6列は、参加者の性別の情報であって、1(男性)あるいは0(女性)のいずれかの値を含んでおり、第7列は、測定条件のうち、測定中に参加者が、目を開いていたか、あるいは、目を閉じていたかを示し、1(開眼)あるいは0(閉眼)の値を含んでいる。また、最後の3つの列は、3つの薬剤の投薬履歴に関するステータス情報を含んでおり、それぞれの列は、1(薬物治療あり)あるい0(薬物治療履歴無し)の値を含んでいる。
なお、人特性情報および測定条件情報は、上記のようなものに限定されず、たとえば、参加者の他の薬剤の投薬履歴の情報などの他の特性や、他の測定条件、たとえば、MRI装置の印加磁場の大きさなどの条件の情報を含んでいてもよい。
第2のデータ行列X2は、参加者の機能的結合の相関(FC)を表わす相関行列の非対角の下三角の部分の要素を行ベクトル形式で表現したものであるとする。
そして、L1−SCCSが、行列のペアX1およびX2に適用され、それから、スパース投影ベクトルw1およびw2が導き出される。
そして、上述のとおり、ハイパーパラメータλ1およびλ2を所定の値に設定したときに、スパース化により第1のデータ行列X1が、スパース投影ベクトルw1により、特定の人属性情報、この場合は、「疾患ラベル」のみとなる正準変数に投影される条件が満たされるようにする。この場合に、対応する第2のデータ行列X2において、疾患ラベルにのみ関連する機能的結合の相関行列要素を特定するためのインデックス(要素)を特定するために、スパース投影ベクトルw2を使用する。
すなわち、内ループ特徴選択では、L1−SCCAのハイパーパラメータλ1およびλ2は、たとえば、0.1の大きさのステップで、0.1と0.9の間で独立して変更される。
ただし、ハイパーパラメータλ1およびλ2を可変とする範囲および変化させるステップの大きさは、このような例に限定されるものではない。
L1−SCCAの処理に対して、「疾患」ラベルにのみ接続した正準変数が存在するようなハイパーパラメータλ1およびλ2の範囲を見つけ出す。
導出されたスパース投影ベクトルw2のゼロでない要素によって規定される部分空間へ元の相関行列要素を投影することを以下のようにして表現する。
ここで、投影ベクトルw2のk番目のゼロでない要素のインデックスを示すために変数ikを定義する。ここで、1≦k≦mであり、mはゼロでないエレメントの数を示す。
そして、以下のような部分空間への投影マトリックスEを考える。
最後に、元の相関行列要素ベクトルx2を以下のように投影することにより部分空間z∈Rmへのベクトルが導かれる。
結果として、疾患ラベル(MDD/HC)に関連する特定の数だけの特徴量(相関行列の要素)を選択することができる。
診断のラベルとのみ関係のある正準変数に相当する相関行列要素を選ぶことによって、分類に本質的な相関行列要素が選択される。
(内ループ特徴抽出)
図7は、内ループ特徴選択を説明するための概念図である。
図6のステップS104の内ループ特徴選択では、図7に説明するようにして、参加者の集合をK個(K:2以上の自然数)に分割して、(K−1)個の部分集合のうち、1つの集合を除いた残りの部分集合について、上述したL1−SCCAを実行することを、(K−1)回繰り返す。
すなわち、図7(a)に示すように、参加者の集合のうちの(K−1)/K(たとえば、K=9であれば8/9)に相当する部分集合(これを、以下「内ループ訓練データ」と呼ぶ)が、このような内ループ特徴選択に使用される。K個の部分集合のうちの残りの1つは、後に説明する外ループ特徴選択において、訓練対象となるテストデータを含む「テストプール」として使用されるため、内ループ特徴選択では使用されない。
図7(b)に示すように、(K−1)個の部分集合のうち、1つの集合(図7(b)で斜線で示す)を除いた残りの部分集合について、ハイパーパラメータλ1およびλ2を所定の範囲で、所定ステップずつ変化させ、L1−SCCAを実行する。特定のハイパーパラメータλ1およびλ2の範囲で、「診断」ラベルとのみ関連がある正準変数に関係している機能的結合要素FCを特徴量として抽出する。
このような抽出処理を、(K−1)個の部分集合のうち、除くことになる1つの集合を順次変えた部分集合について、反復する。
各反復において、抽出された機能的結合の相関行列の要素FCの和集合を「内ループで選択された機能的結合要素FCの和集合」(第1の和集合)とする。
このような処理を行うことで、攪乱変数NVに相当する、異なる撮像サイトでの人属性情報や異なる撮像サイトでの撮像条件の相違によって引き起こされた不適当な影響は縮小されることになる。
この手続きは、後に説明するように、たとえば、日本国内の複数の撮像サイトからのMRI計測データに基づいて生成した第1の分類器を、外国、たとえば、米国にまで汎化するロバストに構築するのに役立つことになる。
再び、図6に戻って、続いて、データ処理部32は、内ループ特徴選択の結果に基づいて、スパースロジスティック回帰を用いて、第2の特徴抽出の処理を実行する(S106)。
以下では、ステップS106の処理を「外ループ特徴選択」と呼ぶ。
ここで、ステップS106での外ループ特徴選択の処理を説明するために、「スパースロジスティック回帰」と、「外ループ特徴選択」について、以下、それぞれ説明する。
(スパースロジスティック回帰)
スパースロジスティック回帰は、ロジスティック回帰分析をベイズ推定の枠組みに拡張した手法であり、特徴ベクトルの次元圧縮を判別のための重み推定と同時に行う手法である。データの特徴ベクトルの次元数が非常に高く、不要な特徴量が多く含まれている場合に有用である。不要な特徴量に対しては線形判別分析における重みパラメータをゼロにし(すなわち、特徴選択を行い)、判別に関連するごく少数の特徴量だけを取り出す(スパース性)。
スパースロジスティック回帰では得られた特徴データに対して、分類するクラスへの所属する確率pをクラス毎に求め、最大値を出力したクラスに割り当てる。pはロジスティック回帰式によって出力される。重みの推定はARD(Automatic Relevance determination)によって行われ、クラス判別への貢献が少ない特徴量は、重みが0に近づくことで計算から除外される。
すなわち、上述のL1正則化CCAを用いて抽出された特徴量の第1の和集合を入力として、次の階層ベイズ推定にもとづく第1の分類器を用いて、疾患ラベルの予測を行う。
ここで、上述した式(2)で抽出された特徴入力z(選択されたFC)から診断のラベルが疾患(ここでは、自閉症との診断)となる確率を予測するために、分類器としてロジスティック回帰を使用する。
ここで、yは、診断クラス・ラベルを表わし、つまり、y=1はMDDクラスを、y=0はHCクラスをそれぞれ示す。
また、以下のzハット(文字の頭部に^が付けられたものを、「ハット」と呼ぶ)は、拡張入力を備えた特徴ベクトルである。
ここで、特徴ベクトルzは、1人の参加者の安静時MRIサンプルの結合相関行列から式(2)に従って、抽出される。
拡張入力“1”を使用することは、第1の分類器に対する一定の(バイアス)入力を導入するための、標準的なアプローチである。
以下のようなθは、ロジスティック関数のパラメータベクトルである。
この際のパラメータθの分布を以下のような正規分布に設定する。
さらに、パラメータwの分布のハイパーパラメータαの分布を以下のように設定することにより、階層ベイズ推定を行うことによりそれぞれのパラメータの分布を推定する。
ここで、a0, b0はハイパーパラメータのガンマ分布を決定するパラメータである。αはベクトルθの正規分布の分散を表現するパラメータベクトルで、そのベクトルのi番目の要素がαiである。
なお、このようなスパースロジスティック回帰については、以下の文献に開示がある。
文献:Okito Yamashita, Masa aki Sato, Taku Yoshioka, Frank Tong, and Yukiyasu Kamitani. “Sparse Estimation automatically selects voxels relevant for the decoding of fMRI activity patterns.” NeuroImage, Vol. 42, No. 4, pp. 1414-1429, 2008.
(外ループ特徴選択)
図11は、外ループ特徴選択を説明するための概念図である。
図6のステップS106の外ループ特徴選択では、図11に説明するようにして、内ループ特徴選択の際に、参加者の集合をK個(K:2以上の自然数)に分割したうち、内ループ特徴選択に使用した(K−1)個の部分集合を除く残りの1つの部分集合を、訓練対象となるテストデータを含む「テストプール」として使用する。
すなわち、図11に示すように、たとえば、このテストプールをL個のグループに分割し、このL個のグループから順次1つのグループ(図11の点線四角内で、斜線で示す)を選択し、参加者の集合この選択されたグループを除く参加者について、第1の和集合である機能的結合の相関行列の要素FCを説明変数として、疾患ラベルを予測する第1の分類器をスパースロジスティック回帰により生成する。このときに、スパースロジスティック回帰により「機能的結合の相関行列の要素FC」が、さらに選択される。これを「外ループで選択された機能的結合の相関行列の要素FC」と呼ぶ。
このようにして、L個のグループから順次1つのグループを選択して、このような処理をL回反復することで、抽出された機能的結合の相関行列の要素FCの和集合を「外ループで選択された機能的結合要素FCの和集合」(第2の和集合)とする。
このように、外ループ特徴抽出において使用されるテスト集合は、内ループ特徴抽出において次元数を縮小するために使用されたデータセットから常に独立している。
なお、ここでは、一例として、特に、L個のグループの各々は、1人の参加者を含むものとする。この場合は、上記反復処理において、除かれた1人の参加者について、生成された第1の分類器で予測処理を行い、その除かれた1人の参加者についての疾患ラベルと予測結果との誤差を累積処理することは、いわゆる1個抜き交差検証(Leave-One-Out Cross Validation:LOOCV)を行うことに相当する。したがって、図11の処理において、各反復処理で、上述の誤差を算出し、これを累積してL回について平均することで、生成された第1の分類器の交差検証、言い換えると、汎化能力の評価を行うことも可能である。
(第1の分類器の生成)
再び図6に戻って、続いて、データ処理部32は、外ループ特徴選択の結果(第2の和集合)に基づいて、スパースロジスティック回帰により、分類器の生成処理を実行する(S108)。
図13は、ステップS108での分類器の生成処理の概念を示す図である。
図13に示すように、ステップS108では、第1の和集合から外ループ特徴選択により抽出された第2の和集合を説明変数として、全ての参加者についてのスパースロジスティック回帰により、分類器が生成される。
生成された分類器を特定するための分類器データ(関数形およびパラメータに関するデータ)3112は、不揮発性記憶装置2080に格納され、後に、以上のような訓練に使用されたのとは異なるMRI測定データ(テストデータ)が入力された際に、テストデータに対する疾患ラベルの推定を実行する際の判別処理に使用される。
すなわち、事前の医師の診断に基づいて、健常群、うつ病患者群に分けられた参加者の脳領域間(関心領域間)の活動度の相関(結合)を測定し、測定結果に対する機械学習により、これらとは異なる新たな参加者についてのテストデータがうつ症状または健常のいずれに相当するかを判別するように生成された分類器が、うつ症状のバイオマーカーとして機能する。
このとき、バイオマーカーの出力である「疾患ラベル」としては、分類器がロジスティック回帰により生成されることから、疾患である確率(または健常である確率)を含んでもよい。たとえば、「疾患である確率は、○○%」といような表示として出力する。このような確率を、「疾患マーカー」として使用することが可能である。
また、分類器の出力となる属性については、必ずしも、疾患の判別に限られず、他の属性についての出力であってもよい。この場合も、その属性について、いずれのクラスに属するかという離散的な判別結果を出力してもよいし、その属性について、あるクラスに属する連続的な値、たとえば確率を出力してもよい。
すなわち、分類器の学習(作成)においては、うつ症状のバイオマーカー作成のため,rs−fcMRIデータを入力として,上述の内ループ特徴選択および外ループ特徴選択により特徴抽出を行い、抽出された特徴量を説明変数として生成された分類器によりうつ症状・健常ラベル判別を行う。
(内ループ特徴選択での処理)
図8は、内ループ特徴選択処理の概念を説明するための図である。
図9は、一例として、内ループ特徴選択を行った際の反復処理のうち、特定のハイパーパラメータλ1およびλ2に対する結果を示す図である。
図8を参照して、入れ子構造になっている特徴選択において、L1−SCCAの反復のうちの1つを表わし、ここでは、正準変数が「疾患ラベル」にのみ接続されている。
図8では、正準変数が「疾患ラベル」にのみ接続されていることを、疾患ラベルがただ1つの正準変数w1 Tx1と、実線で接続されることで表している。
また、正準変数w と正準変数w とを接続する点線の上の記号cは、それぞれの正準変数間の相関係数を表す。
正準変数w が「疾患ラベル」にのみ接続されているときに、対応する正準変数w と接続する機能的結合の相関行列の非対角の下三角行列の要素が、特徴量として選択される。
図9では、例として、最初の外ループの中で、各人属性情報につき少なくとも1つの正準相関を生成する最も小さなハイパーパラメータの組合せを示す。
正準変数は、丸によって表わされる。
左の列の丸は人特性情報と測定条件に由来する、正準変数w を示す。一方で、右の列の丸は、機能的結合(FC)に由来する正準変数w を示す。
上述のとおり、正準変数を接続する点線の上の数字は、正準変数w とw の間の相関係数を表わす。
人属性情報のラベルおよび正準変数w の間の接続は実線または点線で表わされる。
(内ループ特徴選択での処理)
図10は、内ループ特徴選択の処理を、より詳しく説明するためのフローチャートである。
内ループ特徴選択にあたっては、参加者の集合は、図7(a)に示したように、K個の部分集合に分けられ、このうちの特定の1個の部分集合を除く(K−1)個の部分集合が、使用される。
図10を参照して、内ループ特徴選択の処理が開始されると、演算処理装置(CPU)2040は、ハイパーパラメータλ1およびλ2を初期値である(λ1,λ2)=(0.1,0.1)に設定する(S200)。
続いて、変数iの値が1に設定され(S202)、CPU2040は、iの値が、内ループ特徴選択の反復回数(K−1)を超えない場合(S204でN)は、内ループ訓練データのうち、i番目のデータブロックを除き、不揮発性記憶装置2080中に格納された、参加者の人属性情報3104および機能的結合の相関行列データ3106に基づいて、スパース正準相関分析を実行する(S206)。
現在の(λ1,λ2)の値に組について、疾患ラベルとのみ結合する正準変数が存在する場合(S208でY)、CPU2040は、疾患ラベルとのみ結合する正準変数に対応する機能的結合の相関行列の要素(FC)を抽出し、不揮発性記憶装置2080に、記憶する(S210)。ステップS210の処理の後、または、現在の(λ1,λ2)の値に組について、疾患ラベルとのみ結合する正準変数が存在しない場合(S208でN)は、iの値を1だけインクリメントして、処理は、ステップS204に復帰する。
したがって、CPU2040は、(K−1)回にわたって、ステップS206〜S212の処理を繰り返す。
ステップS204において、iの値が、内ループ特徴選択の反復回数(K−1)を超えた場合(S204でY)は、CPU2040は、(λ1,λ2)について、所定のルールにしがたい、λ1およびλ2のいずれか一方を所定のステップ量だけ変化させる(S214)。ステップS216において、(λ1,λ2)の値が可変範囲内であるときは(S216でY)、CPU2040は、処理をステップS202に復帰させる。
(λ1,λ2)の可能な全ての値の組について処理が終了すると(S216でN)、CPU2040は、これまでの処理で抽出された機能的結合の相関行列の要素(FC)の和集合を求め、不揮発性記憶装置2080に、第1の機能的結合和集合データ3108として格納することで処理が終了する。
(外ループ特徴選択での処理)
図12は、外ループ特徴選択の処理を、より詳しく説明するためのフローチャートである。
外ループ特徴選択にあたっては、参加者の集合は、図11で示したように、K個の部分集合に分けられ、このうちの内ループ特徴選択では使用されなかった1個の部分集合がテストプールとして使用される。このテストプールに含まれる参加者の人数をNT人とする。
また、以下では、一例として、1個抜き交差検証(LOOCV)処理を行うものとして説明する。
図12を参照して、外ループ特徴選択の処理が開始されると、演算処理装置(CPU)2040は、変数iの値を1に設定し(S300)、CPU2040は、iの値が、外ループ特徴選択の反復回数NTを超えない場合(S302でN)は、テストプールのうち、i番目のデータを除き、不揮発性記憶装置2080中に格納された、参加者の人属性情報3104および機能的結合の相関行列データ3106に基づいて、スパースロジスティック回帰(SLR)を実行して、テスト分類器を生成する(S304)。
続いて、CPU2040は、生成されたテスト分類器により、除かれたi番目のデータに対するSLRによる予測値を算出する(S306)。
さらに、CPU2040は、テスト分類器を生成する際のスパース化により抽出された機能的結合FCを、選択された特徴量として不揮発性記憶装置2080に格納し、変数iの値を1だけインクリメントする(S310)。
S304〜S308の処理がNT回繰り返された後、変数iの値がNTを超えると(S302でY)、CPU2040は、LOOCVとしての2乗誤差見積の算出する(S312)。
続いて、CPU2040は、NT回の繰り返し処理により、抽出された機能的結合FCの和集合(第2の和集合)を求め、第2の機能的結合和集合データ3110として、不揮発性記憶装置2080に格納して処理が終了する。
(分類器の生成処理)
図13は、最終的に第1の分類器を生成する手続きを説明するための概念図である。
図13に示すように、分類器生成処理では、内ループ特徴選択処理および外ループ特徴選択処理により、最終的に抽出された第2の機能的結合和集合データ3110を説明変数とし、参加者全員について、スパースロジスティック回帰によって、第1の分類器を生成する。生成された第1の分類器を特定するための情報(ロジスティック関数のパラメータベクトルθ等)は、不揮発性記憶装置2080に記憶される。
上述したように、ここまで説明した手続きでは、「うつ病」を疾患ラベルとした場合は、表1に記載の12個の機能的結合が、最終的に抽出された第2の機能的結合和集合データ3110の説明変数として使用され、第1の分類器が生成される。
ただし、ロジスティック回帰式で表現される分類器の生成には表1に記載の全ての機能的結合を使用してもよいが、少なくとも機能的結合識別番号1及び機能的結合識別番号2の一方、あるいは両方を選択して使用してもよい。また、寄与度が高い複数の機能的結合を選択してもよい。このような12個の機能的結合からの分類器に使用する前記機能的結合の選択は、データ処理部32または後述するコンピュータ300のCPUが行ってもよいが、ヒトが行ってもよい。
(判別処理)
図14は、生成した第1の分類器を用いて被検者のrs−fcMRIのデータから疾患ラベルを判別するためにデータ処理部32が行う処理を説明するためのフローチャート図である。
CPU2040は、ステップS401において、MRI撮像部25から安静時の被検者rs−fcMRI測定データ3113をインタフェース部44を介して取得する。
CPU2040は、ステップS402〜S403おいて上記2.に記載の前処理を行い、表1に記載の機能的結合の全て、あるいは一部について相関行列の要素を抽出する。
ステップS404おいて、CPU2040は、抽出した相関行列の要素に対してピアソン相関係数を求める。ステップS405〜S406おいて、CPU2040は、前記相関係数を分類器に入力し、少なくとも1つの機能的結合についてうつ症状の疾患ラベルを判別するための指標値を生成する。あるいは、望ましくはCPU2040は、前記相関係数を分類器に入力し、複数の機能的結合についての関連重み付け和をうつ症状の疾患ラベルを判別するための指標値として生成する。
次にCPU2040は、ステップS407において指標値を基準値と比較する。例えば数6で表されるような数5によるロジスティック回帰の式では、うつ症状の疾患ラベルと健常者のラベルの閾値は0となるように設定されているため、基準値は0となる。たとえば、「関連重み付け和」をバイオマーカーの指標値として使用する場合は、第1の分類器の感度や特異度を調整するために、基準値を0とは異なる値に設定してもよい。たとえば、基準値を0よりも小さくすれば、「うつ病」の疾患レベルと判別される感度は向上するが、特異度は低下する。上述した分類器の判別結果を、後続する他の情報に基づく医師による診断のための参考情報(支援情報)として使用することを意図する場合は、感度の向上を優先させるという設定もあり得る。
ステップS408にお
いてCPU2040は、前記指標値が0よりも高い場合(「YES:Y」)には、前記被検者がうつ症状のラベルであると判別することができる(S409)。また、前記指標値が0よりも低い場合(「NO:N」)には、前記被検者がうつ症状のラベルに相当すると判別することができる(S410)。
3−2.判別装置2
実施の形態1では、脳活動計測装置(fMRI装置)で、1つの測定場所で計測された脳活動データを計測し、この脳活動データに基づいて、分類器の生成および分類器による疾患ラベルの推定(予測)を、同一のコンピュータによる処理または分散処理により行う構成を説明した。
ただし、i)分類器を機械学習により訓練するための脳活動データの計測(データ収集)、ii)分類器の機械学習による生成処理および特定の被験者についての分類器による疾患ラベルの推定(予測)の処理(推定処理)、iii)上記特定の参加者についての脳活動データの計測(対象者の脳活動計測)を、それぞれ、異なる施設で分散して実行する構成とすることも可能である。
すなわち、健常群およびうつ病患者群のデータは、MRI装置10自身で測定される場合に限られず、他のMRI装置において測定されたデータをも総合して、分類器の生成が行われてもよい。また、より一般には、データ処理部32は、MRI装置の制御を実行するためのコンピュータである必要は必ずしもなく、1つまたは複数のMRI装置からの測定データを受け取って、分類器の生成処理および生成された分類器により、判別処理を行うことに特化したコンピュータであってもよい。
図15は、本発明の実施の形態1の別実施形態である、データ収集、推定処理および対象者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。
図15を参照して、サイト100.1〜100.Nは、うつ病患者群、健常者群を含む参加者のデータを脳活動計測装置により計測する施設であり、管理サーバ200は、サイト100.1〜100.Nからの計測データを管理する。
判別装置2に相当するコンピュータ300は、管理サーバ200に格納されたデータから分類器を生成する。
MRI装置410は、コンピュータ300上の分類器の結果を利用する別サイトに設けられており、特定の被検者について脳活動のデータを計測する。
コンピュータ400は、MRI装置410が設けられる別サイトに設置され、MRI装置410の測定データから特定被検者の脳の機能的結合の相関データを算出し、機能的結合の相関データをコンピュータ300に送信して、返信されてくる分類器による判別結果を利用する。
サーバ200は、サイト100.1〜100.Nから送信されてくるうつ病患者群および健常者群の参加者FMRI測定データ3102と、参加者FMRI測定データ3102に関連付けられた参加者の人属性情報3104とを格納し、コンピュータ300からのアクセスに従って、これらのデータをコンピュータ300に送信する。
コンピュータ300は、通信インタフェース2090を介して、サーバ200からの参加者FMRI測定データ3102および被験者の人属性情報3104を受信する。
なお、サーバ200、コンピュータ300、コンピュータ400のハードウェアの構成は、基本的に、図2で説明した「データ処理部32」の構成と同様であるので、その説明は繰り返さない。
図15に戻って、相関行列算出部3002、第1特徴選択部3004、第2特徴選択部3006、分類器生成部3008および判別処理部3010、ならびに、機能的結合の相関行列のデータ3106、第1の機能的結合和集合データ3108、第2の機能的結合和集合データ3110および分類器データ3112については、実施の形態1で説明したのと同様であるので、その説明は、繰り返さない。
MRI装置410は、診断ラベルの推定対象となる被検者の脳活動データを計測し、コンピュータ400の処理装置4040は、計測されたMRI測定データ4102を不揮発性記憶装置4100に格納する。
さらに、コンピュータ400の処理装置4040は、MRI測定データ4102とMRI装置410の計測対象の被験者の人属性情報とに基づいて、相関行列算出部3002と同様にして、機能的結合の相関行列のデータ4106を算出し、不揮発性記憶装置4100に格納する。
コンピュータ400のユーザから診断の対象となる疾患が指定され、当該ユーザの送信の指示に従い、コンピュータ400が、機能的結合の相関行列のデータ4106をコンピュータ300に送信する。これに応じて、判別処理部3010は、指定された診断ラベルについての判別結果を算出し、コンピュータ300は、通信インタフェース2090を介して、コンピュータ400に送信する。
コンピュータ400では、図示しない表示装置などを介して、ユーザに対して、判別結果を知らせる。
このような構成とすることで、より多くの被験者について収集したデータに基づいて、より多くのユーザに対して、分類器による診断ラベルの推定結果を提供することが可能となる。
また、サーバ200とコンピュータ300とを別個の管理者が管理する形態とすることも可能で、その場合、サーバ200にアクセスできるコンピュータを制限することで、サーバ200に格納される被験者の情報のセキュリティを向上させることも可能となる。
さらに、コンピュータ300の運営主体からみると、「分類器による判別のサービスを受ける側(コンピュータ400)」に対して、分類器についての情報は、一切提供しなくても、「判別結果を提供するサービス」を行うことが可能となる。
3−4.分類器及び判別装置の用途
前記判別装置1及び2は、うつ症状に対する疾患ラベルを判別する他、うつ症状のレベルを判別するために、うつ症状の治療の効果の程度を判別した情報を生成するために、又はうつ症状が予め設定された複数のサブクラスに層別化されている場合に、又は被検者を前記サブクラスに層別化する際の支援情報を生成するために使用することができる。それぞれの使用方法については、後述する。
第1の分類器は、うつ症状を判別するための診断バイオマーカーとして使用することができる。また、第1の分類器は、うつ症状のレベルを判別するための診断バイオマーカーとして使用することができる。第1の分類器は、うつ症状の治療の効果を判別するための薬効マーカーとして使用することができる。第1の分類器は、被検者をうつ病のサブクラスに層別化するためのバイオマーカーとして使用することができる。
3−5.分類器を生成するためのコンピュータプログラムとうつ症状を判別するためのコンピュータプログラム
実施の形態1の別形態には、ステップS100〜S108(具体的には、S200〜S218とS300〜S314を含む)を含む処理を実行させて、分類器生成処理の機能をコンピュータに実行させるためのプログラムと、ステップS401〜S410を含む処理を実行させて、判別処理の機能をコンピュータに実行させるためのプログラムとを含む。また実施の形態1の別形態には、ステップS100〜S108(具体的には、S200〜S218とS300〜S314を含む)と、ステップS401〜S410を含む処理を実行させて、判別装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
3−6.うつ症状の判別方法
実施の形態2は、うつ症状の判別方法に関する。前記判別方法は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程と、前記指標値が基準値を超えた場合に、前記被検者がうつ症状であると決定する工程を備える。具体的な工程は、上記S401〜S410に準ずるが、前記工程は、指標値を生成する工程を除いた工程の全部または一部をヒトが行ってもよい。
4.うつ症状のレベルの判別装置及び判別方法
4−1.判別装置
判別装置1及び判別装置2は、うつ症状のレベルを判別するために使用することができる。すなわち、実施の形態1の別形態であるうつ症状のレベルの判別装置は、演算装置と記憶装置とを有する。この演算装置は、記憶装置に格納されたプログラムに基づいて、前記第1の分類器により、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成し、前記指標値をあらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較し、前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する、という処理を実行する。この場合、判別装置1及び判別装置2の構成は、上記3.に記載した構成と同様であるため、説明は省略する。判別処理部3010は、分類器データ3112により特定される分類器に基づいて、入力データに対して判別処理を行い、前記被検者のうつ症状のレベルを判別する。
この場合、「うつ症状のレベル」の判別には、たとえば、指標値である「関連重み付け和」の値が取り得る値を、予め、複数のレベルの基準範囲に分類しておくことで、「うつ症状のレベル」に対応した判別結果を出力することが可能となる。
あるいは、この場合、分類器の生成に用いる「疾患ラベル」を予めうつ病のレベルに応じた複数の「レベル別疾患ラベル」として、分類器の機械学習を実行してもよい。
4−2.判別方法
実施の形態3は、うつ症状のレベルの判別方法に関する。前記判別方法は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程と、あらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較する工程と、前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する工程とを備える。
図16は、生成した分類器を用いて被検者のrs−fcMRIのデータからうつ症状のレベルを判別するためにデータ処理部32が行う処理を説明するためのフローチャート図である。
CPU2040は、入力部40からの処理開始の入力等により、図15に記載のステップS401〜S406を行うことにより、指標値を生成する。次にCPU2040は、図16に示すステップS501において前記指標値を予めうつ症状のレベルに応じて設定された指標値の基準範囲と比較する。前記12対の機能的結合は、たとえば、BDIのうつ症状の重症度と相関するため、うつ症状のレベルを決定するための指標値の基準範囲は、BDIのうつ症状の重症度との相関から予め決定しておくことができる。ステップS502においてCPU2040は、図15のステップS406で生成された指標値が、どのレベルのうつ症状に対応しているかを決定する(S502)。そして、CPU2040は、ステップS503において、前記被検者がステップS502で決定されたうつ症状のレベルであると判別する。前記CPU2040が行う各ステップは、指標値を生成する工程を除いた工程の全部または一部をヒトが行ってもよい。
また、実施の形態3の別態様として、上記ステップS501〜S503を含む処理を実行させて、判別装置としての機能をコンピュータに実行させるプログラムを含む。実施の形態3の別態様として、上記ステップS401〜S406およびステップS501〜S503を含む処理を実行させて、判別装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
5.第1の分類器を使用する治療効果の判別装置及び判別方法
5−1.判別装置
判別装置1及び判別装置2は、治療効果を判別するための情報の生成のために使用することができる。すなわち、実施の形態1の別形態である被検者における治療効果の判定するための判別装置は、演算装置と記憶装置とを有する。演算装置は、記憶装置に格納されたプログラムに基づいて、前記第1の分類器により、第1の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第1の値を生成し、前記分類器により、治療開始後であって第1の時点より後の第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素についてうつ症状を評価する第2の値を生成し、前記第1の値と前記第2の値とを比較し、i)前記第2の値が前記第1の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定する処理、及び/又は前記第2の値が前記第1の値よりも改善していない場合に、前記治療が被検者のうつ症状の改善に有効でないと決定する処理を実行する。この場合、判別装置1及び判別装置2の構成は、上記3.に記載した構成と同様であるため、説明は省略する。判別処理部3010は、分類器データ3112により特定される分類器に基づいて、入力データに対して判別処理を行い、前記被検者における治療効果を判別する情報を生成する。
本実施の態様の別の態様は、前記判別装置を、ドラッグリプロファイリングのために使用することを含む。
5−2.判別方法
実施の形態4は、前記第1の分類器に、少なくとも2つの時点において撮像された被検者の安静時fMRIのデータから取得される表1に記載の機能的結合についての相関行列の要素を入力し、治療効果を判別する情報を生成する。具体的には、第1の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第1の値を生成する工程と、第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素ついてうつ症状を評価する第2の値を生成する工程と、第1の値と第2の値とを比較する工程と、第2の値が第1の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定し、有効である可能性を示す情報を生成する工程とを備える。第1の時点は、治療前であっても、治療開始から所定の時間が経過していてもよいが、好ましくは、治療開始前である。また、一度治療を終了した後、治療が再度開始される場合も治療終了から治療再開前までの間は、「治療前」としてもよい。第2の時点は、治療開始後であって第1の時点よりも後である限り制限されない。
図17は、第1の分類器を用いて被検者のfMRIのデータから治療効果を判別するための処理を説明するためのフローチャート図である。
CPU2040は、入力部40からの処理開始の入力等により、図17に記載のステップS601において、第1の時点においてMRI撮像部25によって被検者から撮像された安静時の被検者fMRI測定データ3113をインタフェース部44を介して取得する。
ステップS602〜S603おいて上記2.に記載の前処理を行い、表1に記載の12対の機能的結合から選択した全て、あるいは一部の機能的結合について第1の時点での相関行列の要素を抽出する。
ステップS604において、CPU2040は、ステップS603で抽出した第1の時点での相関行列の要素に対してピアソン相関係数を算出する。
ステップS605〜S606おいて、CPU2040は、ステップS604において算出された相関係数を第1の分類器に入力し、前記選択された機能的結合について第1の値を生成する。前記第1の値は、例えば、上記2.で述べた指標値そして算出される。
CPU2040は、第2の時点において、ステップS607に示すMRI撮像部25によって被検者から撮像された安静時の被検者fMRI測定データ3113をインタフェース部44を介して取得する。
ステップS608〜S609おいて上記2.に記載の前処理を行い、表1に記載の12対の機能的結合から選択した全て、あるいは一部の機能的結合について第2の時点での相関行列の要素を抽出する。
ステップS610において、CPU2040は、ステップS609で抽出した第1の時点での相関行列の要素に対してピアソン相関係数を算出する。
ステップS611〜S612おいて、CPU2040は、ステップS610において算出された相関係数を第1の分類器に入力し、前記選択された機能的結合について第1の値を生成する。前記第2の値は、例えば、上記2.で述べた指標値として生成される。
CPU2040は、ステップS606において生成された第1の値とステップS612おいて生成された第2の値を比較する(S613)。
CPU2040は、ステップS614において、第1の値よりも第2の値が改善していると判別した場合には(「YES:Y」)、前記被検者のうつ病症状の改善に対する前記治療の状態を示す情報を生成する(S615)。この場合、次にCPU2040は、ステップS616に進み、前記治療の継続の可能性を示す情報を提示してもよい。医師は、このような治療継続の可能性を示す情報を元に、治療継続をするか否かの判断を行う。特に限定されないが、たとえば、出力画面上には、治療の有効性があると予め設定された領域についての情報が区別可能に表示されており、「治療の状態を示す情報」とは、分類器の出力する指標値であり、「治療を継続の可能性を示す情報の提示」とは、上記領域を表示した状態で、併せて、この指標値を表示することであってもよい。
また、CPU2040は、ステップS614において、第1の値よりも第2の値が改善していないと判別した場合には(「NO:N」)、前記被検者のうつ病症状の改善に対する前記治療の状態を示す情報を生成する(S617)。この場合、次にCPU2040は、ステップS618に進み、前記治療の終了の可能性を示す情報を提示してもよい。
医師は、このような治療終了の可能性を示す情報を元に、治療を終了をするか否かの判断を行う。特に限定されないが、たとえば、出力画面上には、治療の有効性がないと予め設定された領域についての情報が区別可能に表示されており、「治療の状態を示す情報」とは、分類器の出力する指標値であり、「治療の終了の可能性を示す情報の提示」とは、上記領域を表示した状態で、併せて、この指標値を表示することであってもよい。
さらに、CPU2040は、ステップS619に進み、他の治療方法に変更可能性を示す情報を提示してもよい。この場合、さらに具体的な治療方法の可能性を提示してもよい。この場合、医師は、表示された治療方法の可能性のうちから、現在の治療方法から変更すべき治療方法を判断することができる。
例えば数6で表される式では、うつ症状の疾患ラベルと健常者のラベルの閾値は0となるように設定されているため、基準値は0となる。
たとえば、「関連重み付け和」をバイオマーカーの指標値として使用する場合は、ステップS614において、うつ症状が重くなるにつれ、前記指標値は、正方向に大きくなる。したがって、CPU2040は、前記第1の値よりも前記第2の値が小さくなっている場合に、第1の値よりも第2の値が改善していると判別することができる。また、CPU2040は、前記第1の値よりも前記第2の値が大きくなっている場合、または第1の値と第2の値に顕著な差が認められない場合に、第1の値よりも第2の値が改善していないと判別することができる。
前記第1の値及び第2の値は、指標値であってもよいが、ステップS604及びステップS610で算出される相関係数を、それぞれ第1の値及び第2の値としてもよい。
また、実施の態様3では、ステップS601は、ステップS607よりも前に行われている必要があるが、ステップS602〜S606は、ステップS607よりも前に行われている必要はない。ステップS602〜S606は、ステップS601の後、少なくともステップS613よりも前のどこかで行われていればよい。
ステップS601〜S619は、その全部あるいは一部をヒトが行ってもよい。
実施の形態3は、投薬治療の効果判定に使用することも可能である。第1の分類器を使って投薬効果を判定する場合には、少なくとも機能的結合識別番号1及び2の機能的結合を指標とすることが好ましい。
実施の態様4の別の態様は、前記判定方法を、ドラッグリプロファイリングのために使用することを含む。
また、実施の形態4の別態様として、上記ステップS601〜S619を含む処理を実行させて、判別装置としての機能をコンピュータに実行させるプログラムを含む。このプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
6.うつ病患者の層別化装置及び層別化方法
6−1.層別化装置
実施の形態1の別形態は、うつ病患者の層別化装置に関する。うつ病が予め設定された複数のサブクラスに層別化されている場合に、判別装置1及び判別装置2は、それぞれうつ病患者を医師が層別化するための支援をする層別化装置1及び層別化装置2(以下、総称するときは、「層別化装置」という)として使用することができる。層別化装置は、演算装置と記憶装置とを含む。すなわち、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、うつ病患者の層別化するための判別装置の演算装置は、記憶装置に格納されたプログラムに基づいて、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成し、、前記指標値を、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較し、、前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する、という処理を実行する。この場合、層別化装置1及び層別化装置2の構成は、上記3.に記載した判別装置1及び判別装置2の構成と同様であるため、説明は省略する。図19と図20には、層別化装置1及び層別化装置2の機能ブロック図を示す。図19と図20では、判別処理部3010が、層別化処理部3020となる点を除き、図5及び図15と同じである。層別化処理部3020は、分類器データ3112により特定される分類器に基づいて、入力データに対して層別化処理を行い、うつ病患者の層別化のための判別を行い結果を出力する。
医師は、層別化装置1または層別化装置2にから出力される、うつ病患者の層別化のための判別結果を支援情報として、患者の層別化を行う。この意味で、このような装置を、層別化支援装置と呼ぶこともできる。
6−2.層別化方法
実施の形態5は、うつ病患者得層別化するための層別化する方法に関する。前記判別方法は、うつ病が予め設定された複数のサブクラスに層別化されている場合に、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程と、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較する工程と、前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する工程を備える。
図18は、生成した分類器を用いて被検者のfMRIのデータからうつ症状のサブクラスを判別するための情報を生成するデータ処理部32が行う処理を説明するためのフローチャート図である。
CPU2040は、入力部40からの処理開始の入力等により、図14に記載のステップS401〜S406を行うことにより、指標値を生成する。次にCPU2040は、図18に示すステップS701において前記指標値を、予めうつ病のサブクラスに応じて設定された指標値の基準範囲と比較する。うつ病のサブクラスは、臨床所見からメランコリー型MDD、非メランコリー型MDD及び治療抵抗性MDDのサブクラスに分類される。このため、たとえば、メランコリー型MDDを「疾患ラベル」として生成した分類器を用いることで、脳活動を計測した対象者からの計測データが、メランコリー型MDDの疾患ラベルに対応するか否かを判別できる。同様にして、非メランコリー型MDDまたは治療抵抗性MDDを「疾患ラベル」として生成した分類器を用いることで、脳活動を計測した対象者からの計測データが、非メランコリー型MDDまたは治療抵抗性MDDの疾患ラベルに対応するか否かを判別できる。
ステップS702においてCPU2040は、図15のステップS406で生成された指標値が、どのサブクラスに対応しているかを決定する。そして、CPU2040は、ステップS703において、前記被検者がステップS702で決定されたうつ病のサブクラスに属すると判別する。医師は、このようにして出力されるサブクラスの情報を、うつ病患者の層別化のための支援情報として、患者の層別化を行う。この意味で、このような方法を、層別化支援方法と呼ぶこともできる。前記CPU2040が行う各ステップは、指標値を生成する工程を除いた全部または一部をヒトが行ってもよい。
また、実施の形態5の別態様として、上記ステップS701〜S703を含む処理を実行させて、層別化装置としての機能をコンピュータに実行させるプログラムを含む。実施の形態5の別態様として、上記ステップS401〜S406およびステップS701〜S703を含む処理を実行させて、層別化装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
7.第2の分類器の生成と第2の分類器を使用する治療効果の判別
7−1.治療効果判別装置
実施の形態6は、複数の対象者について、治療のために投与した治療効果がみられる群とみられない群とを判別するための第2の分類器を生成する分類器生成処理を実行し、医師が治療効果の判断をする支援のための指標となる情報を生成するする、治療効果判別装置、あるいは、治療支援装置に関する。
治療効果判別装置の構成は、基本的に図1と同じであるが、図1におけるデータ処理部32を実施の形態6ではデータ処理部62と読み替える。また、治療効果判別装置において、データ処理部62の構成は、図1及び図2と基本的に同様であるが、実施の形態6では、図1において記憶部36を記憶部66と、制御部42を制御部62と、入力部40を入力部70と、インタフェース部44をインタフェース部74と、データ収集部46をデータ収集部76と、画像処理部48を画像処理部78と、表示部38を表示部68と、表示制御部34を表示制御部64と読み替える。また、図2において、データ処理部32のコンピュータ本体2010をコンピュータ本体6010と、メモリドライブ2020をメモリドライブ6020、ディスクドライブ2030をディスクドライブ6030と、演算装置(CPU)2040を演算装置(CPU)6040と、ディスクドライブ2030をディスクドライブ6030と、メモリドライブ2020をメモリドライブ6020と、バス2050をバス6050と、ROM2060をROM6060と、RAM2070をRAM6070と、不揮発性記憶装置2080を不揮発性記憶装置6080と、通信インタフェース2090を通信インタフェース6090と読み替える。データ処理部62のハードウェアとしては、上述のとおり、特に限定されないが、汎用コンピュータを使用することが可能である。
(分類器の生成)
図23は、第2の分類器の生成処理および生成された第2の分類器による判別処理を行うための機能ブロック図である。
不揮発性記憶装置6080には、健常者、治療を受けたうつ病患者を含む複数の対象者の脳内の複数の所定領域における脳活動を示す信号をMRI装置により時系列で予め測定した信号の情報であって、第1の時点で測定した対象者のrs−fcMRI測定データ6102と、第2の時点で測定した対象者のrs−fcMRI測定データ6104と、前記rs−fcMRI測定データを計測した各対象者に関連づけられた、複数の治療データ6105とが格納されている。
ここで、治療データ6105には、対象者の疾患ラベル、治療履歴(治療方法、投薬した薬剤名、用量、投与期間等)の情報が含まれる。CPU6040は、第1の時点の対象者のrs−fcMRI測定データ6102と、第2の時点の対象者のrs−fcMRI測定データ6104と、対象者の治療データ6105に基づいて、複数の対象者において治療効果がみられた群(「寛解群(remitted)」)と治療効果がみられない群(「非寛解群(non-remitted)」)とを識別するための分類器を生成する処理を行う。したがって、疾患ラベルには、この場合、「寛解ラベル」と「非寛解ラベル」とが含まれる。
第1の相関計測部6002は、第1の時点の対象者のrs−fcMRI測定データ6102から、複数の所定領域間の脳活動の機能的結合の相関行列を対象者ごとに算出する。算出された機能的結合の相関行列のデータは、対象者ごとに、不揮発性記憶装置6080に、機能的結合の相関行列の第1の相関計測データ6106として格納される。
第2の相関計測部6004は、第2の時点の対象者のrs−fcMRI測定データ6104から、複数の所定領域間の脳活動の機能的結合の相関行列を対象者ごとに算出する。算出された機能的結合の相関行列のデータは、対象者ごとに、不揮発性記憶装置6080に、機能的結合の相関行列の第2の相関計測データ6108として格納される。
第2の分類器生成部6008は、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差に基づいて、回帰式を求める。前記回帰式は、第2の分類器を特定するための第2の分類器データ6020として、不揮発性記憶装置6080に格納する。治療効果判別処理部6100は、第2の分類器データ6020により特定される第2の分類器に基づいて、入力データに対して治療効果の判別処理を行う。
図21は、第2の分類器を生成するためにデータ処理部62が行う処理を説明するためのフローチャートを示す。
以下、図23で説明した処理を図21を参照しながら、より詳しく説明する。
図21において、データ処理部62は、入力部70からの処理開始の入力等により、分類器生成の処理が開始されると(スタート)、第1の時点の各対象者のrs−fcMRI測定データ6102を読み出し、各表1に示される複数の機能的結合それぞれについて、第1の時点の相関を計測する(S802)。この相関の計測は、各機能的結合のピアソン相関係数を算出することによって行われる。
次にデータ処理部62は、記憶部66から第2の時点の各対象者のrs−fcMRI測定データ6104を読み出し、各対象者について表1に示される複数の機能的結合それぞれについて、第2の時点の相関を計測する(S803)。この相関の計測は、各機能的結合のピアソン相関係数を算出することによって行われる。
次に、データ処理部62は、各対象者の表1に示される複数の機能的結合のうち、たとえば、第1の機能的結合および第2の機能的結合のそれぞれについて、第1の時点の相関と第2の時点の相関との差を算出する(S804)。前記差は、例えば、下式(4)に基づいて算出される。
式(4)において、FC1は、表1における第1番目の結合強度を示し、FC2は、表1における第2番目の結合強度を示す。すなわち、関連重み付け和の計算の結果、重みの絶対値が最も大きな結合がFC1であり、2番目に大きな結合がFC2である。上付きの“post”は、第2の時点を示し、“pre”は、第1の時点を示す。典型的には、第1の時点は、治療前(投薬前)の時点であり、第2の時点は、治療開始後(投薬開始後)、所定の時点である。また、sign(w1)は、FC1の重みの符号を示し、sign(w2)は、FC2の重みの符号を示す。たとえば、後述するように図33fにおいて作成したこの分類器は、正確性:0.75, AUC:0.79, 特異度0.88、感度0.43を示す。
さらに、データ処理部62は、各対象者の各機能的結合における前記2つの相関の差と、治療データ6105に含まれる説明変数(「寛解群」または「非寛解群」のラベル)との関係を説明するための回帰式を求める。回帰式は、特に制限されないが、線形回帰分析によって求めることができる。この回帰式に基づいて、相関状態空間における治療効果あり群と、治療効果なし群を識別する分類器を生成することができる(S804)。
ここで、「相関状態空間」とは、このようにして複数の機能的結合の第1および第2の時点の間の強度差(必要により、さらに強度差に重みの符号をかけたもの)をそれぞれ空間の軸として張られる空間のことをいう。上記の例では、2次元の空間となる。
上述のとおり、例えば、各対象者の表1に示す機能的結合識別番号1における第1の時点と第2の時点の前記相関の差(Δsign(W)FC1とする)と、表1に示す機能的結合識別番号2における第1の時点と第2の時点の前記相関の差(Δsign(W)FC2とする)の相関状態をプロットすると図33fのようになる。Δsign(W)FC1とΔsign(W)FC2の相関と、各対象者における治療効果の有無との関係をみると負の傾きをもつ線形回帰式を求めることができる。そして、この線形回帰式に基づく回帰直線を基準として、回帰直線よりも上に分布する対象者を治療効果が認められない群、回帰直線よりも下に分布する対象者を治療効果が認められる群として識別することができる。
医師は、このようにして出力される治療効果に関する情報を、支援情報として、うつ病患者に対する治療効果の判断を行う。
(説明変数としてのΔsign(W)FC1とΔsign(W)FC2)
上記の例では、各結合強度の重みの符号sign(W)は、図として表示する際の便宜のためであり、本質的には、治療(投薬)前と治療(投薬開始後)後におけるFC1とFC2の強度変化が重要である。
図33fのように、寛解と非寛解を区別せずに見れば、第1の機能的結合の強度変化Δsign(W)FC1と第2の機能的結合の強度変化Δsign(W)FC2との間には、治療開始前と治療開始後所定期間経過後の強度変化の間には相関が見られず、各々の変化に対する補償のような現象は見られない。言い換えれば、これらの結合強度の変化は、相互に独立であることになり、治療の効果を説明するための変数の選択としては適切であることを示している。
したがって、上記のように、治療開始の前後における強度変化が相互に独立と見なせる機能的結合を用いて、上記のような第1および第2の機能的結合に限られず、他の機能的結合を用いることも可能である。また、治療効果を説明するための変数としては、より多くの機能的結合を、表1に示した機能的結合から選択してもよい。
(判別処理)
図22は、実施の形態7生成した第2の分類器を用いて被検者のrs−fcMRIのデータから治療効果を判別するデータ処理部62が行う処理を説明するためのフローチャート図である。
CPU6040は、入力部70からの処理開始の入力等により、ステップS901において、MRI撮像部25から第1の時点における安静時の被検者rs−fcMRI測定データ6113をインタフェース部64を介して取得する。
CPU6040は、ステップS902〜S903おいて上記2.に記載の前処理を行い、表1に記載の機能的結合の全て、あるいは一部について相関行列の要素を抽出し、第1の時点の各機能的結合の相関を計測する。
CPU6040は、ステップS904において、MRI撮像部25から第2の時点における安静時の被検者rs−fcMRI測定データ6113をインタフェース部64を介して取得する。
CPU6040は、ステップS905〜S906おいて上記2.に記載の前処理を行い、表1に記載の機能的結合の全て、あるいは一部について相関行列の要素を抽出し、第2の時点の各機能的結合の相関を計測する。
次に、CPU6040は、ステップS907おいて、前記複数の機能的結合それぞれについて、ステップS903で計測された第1の時点の各機能的結合の相関とステップS906で計測された第2の時点の各機能的結合の相関の差を算出する。
CPU6040は、ステップS907で算出された差を第2の分類器に入力し(S908)、第2の分類器によって前記被検者において治療効果が認められているか、認められていないかを判別する(S908)。
各相関の計測方法、各相関の差の算出方法、及び治療効果の判別方法は、第2の分類器の生成方法で記載した内容と同様である。
ここで、第1の時点は、治療開始前であっても、治療開始から所定の時間が経過していてもよいが、好ましくは、治療開始前である。また、一度治療を終了した後、治療が再度開始される場合も治療終了から治療再開前までの間は、「治療開始前」としてもよい。第2の時点は、治療開始後であって第1の時点よりも後である限り制限されない。
また、実施の形態7は、ステップS901は、ステップS606よりも前に行われている必要があるが、ステップS902〜S903は、ステップS904よりも前に行われている必要はない。ステップS902〜S903は、ステップS901の後、少なくともステップS907よりも前のどこかで行われていればよい。
また、実施の形態6の別形態は、前記第2の分類器、治療効果判別装置、あるいは、治療支援装置を、ドラッグリプロファイリングのために使用することを含む。
7−2.治療効果判別方法
実施の形態7は、第2の分類器を用いた治療効果の判別方法に関する。前記判別方法は、第1の時点の被検者の安静時の脳内における、前記複数の機能的結合の第1の相関を計測する工程と、第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第2相関を計測する工程と、前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する工程とを備える。具体的な工程は、上記S901〜S909に準ずるが、前記工程は、全部または一部をヒトが行ってもよい。
ここでも、このようにして出力される治療効果に関する判別情報を、支援情報として、うつ病患者に対する治療効果の判断を行う。
実施の形態7の別の形態は、前記治療効果判別方法を、ドラッグリプロファイリングのために使用することを含む。
また、実施の形態7の別形態として、上記ステップS801〜S804を含む処理を実行させて、第2の分類器としての機能をコンピュータに実行させるプログラムを含む。実施の形態7の別態様として、上記ステップS801〜S804およびステップS901〜S909を含む処理を実行させて、治療効果判別装置としての機能をコンピュータに実行させるプログラムを含む。実施の形態7の別態様として、上記ステップS901〜S909を含む処理を実行させて、治療効果判別装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
8.層別化されたうつ病患者における第2の分類器を使用した治療効果の判定
8−1.治療効果判別装置
実施の形態8は、実施の態様4に記載のうつ病患者の層別化装置において、層別化されたうつ病患者からのrs−fcMRIデータを使って、実施の形態6の治療効果判別装置において、治療効果を判別する装置である。
すなわち、実施の態様8は、上記3.で述べた第1の分類器生成処理と、上記6.で述べた層別化処理と上記7.で述べた第2の分類器生成処理と、第1の相関計測処理と、第2の相関計測処理と、前記第2の分類器により、前記被検者に対する治療効果を判別する判別処理とを実行する治療効果判別装置である。
したがって、上記3.の記載、上記6−1.の記載と、上記7−1.記載は、ここに援用される。
実施の形態8の別の形態は、前記治療効果判別装置を、ドラッグリプロファイリングのために使用することを含む。
8−2.治療効果判別方法
実施の形態9は、層別化されたうつ病患者における第2の分類器を使用した治療効果の判別方法に関する。
図24に、層別化されたうつ病患者における第2の分類器を使用した治療効果を判別するためのフローチャートを示す。
図24において、ステップS701〜S703は、図18と同じである。したがって、上記6−2.の説明は、ここに援用される。
CPU6040は、図24に示すステップS703においてサブクラスが決定された被検者について、ステップS921〜S929を行う。ステップS921〜S929の各ステップは、上記7−2.及び図22に記載のステップS901〜S909の各ステップに対応する。したがって、上記7−2.の記載は、ここに援用される。
なお、特に限定されないが、以上の説明において、治療効果判別装置の支援により、医師が、現在継続中の治療方法(特定薬の投薬など)に効果が見られないと判断した場合は、それ以後に、医師は、別の治療方法(別の特定薬の投与、ニューロフィードバック法、修正型電気けいれん療法、反復経頭蓋磁気刺激治療法など)の治療方法を併用する、あるいは、別の治療方法へと治療方法を変更する、という判断をすることができる。
実施の形態9の別の形態は、前記治療効果判別方法を、ドラッグリプロファイリングのために使用することを含む。
また、実施の形態8の別形態として、上記ステップS921〜S929を含む処理を実行させて、第2の分類器としての機能をコンピュータに実行させるプログラムを含む。実施の形態8の別形態として、上記ステップS701〜S703およびステップS921〜S929を含む処理を実行させて、治療効果判別装置としての機能をコンピュータに実行させるプログラムを含む。実施の形態7の別形態として、上記ステップS401〜S406、ステップS701〜S703およびステップS901〜S909を含む処理を実行させて、治療効果判別装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
9.ニューロフィードバック
9−1.脳活動訓練装置
実施の形態10は、被訓練者の脳活動に対してフィードバックを与え、健常者の相関(結合)状態に近づける訓練(ニューロフィードバック訓練、あるいは単に訓練ともいう)ための脳活動訓練装置に関する。
なお、より一般的には、このような脳活動訓練装置は、健常群と患者群の脳活動という関係において、被訓練者の脳活動の機能的結合状態を健常者の脳活動の機能的結合状態に近づける訓練に留まらず、目標とする脳活動の機能的結合状態に、被訓練者の現在の脳活動の機能的結合状態を近づけるという意味での訓練に使用することも可能である。健常者の脳活動の機能的結合状態及び目標とする脳活動の結合状態を、後述するように「目標パターン」と呼ぶ。
ここで、被訓練者は、上記3.記載のうつ症状の判別装置1、うつ症状の判別装置2又はうつ症状の判別方法により、「うつ症状」のラベルを付された被検体(好ましくはMDD、より好ましくはメランコリー型MDDの被検体)、上記4.に記載のうつ症状のレベルの判別装置及び判別方法により、うつ症状のレベルを判定された被検体(好ましくはMDD、より好ましくはメランコリー型MDDの被検体)、上記6.に記載のうつ病患者の層別化装置及び層別化方法により、層別化された被検体(好ましくはMDDに層別化された被検体、より好ましくはメランコリー型MDDに層別化された被検体)であることが好ましい。従って、上記記3.記載のうつ症状の判別装置1、うつ症状の判別装置2又はうつ症状の判別方法に関する説明、上記4.に記載のうつ症状のレベルの判別装置及び判別方法に関する説明及び上記6.に記載のうつ病患者の層別化装置及び層別化方法に関する説明は、ここに援用される。
図34は、このような脳活動訓練装置の構成の概念を示す図である。
脳活動訓練装置のハードウェア構成としては、たとえば、上述した図1に示されるMRI装置10の構成を用いることができる。脳活動訓練装置のハードウェア構成は、基本的に図1と同じであるが、図1におけるデータ処理部32を実施の態様10ではデータ処理部102と読み替える。また、脳活動訓練装置において、データ処理部102の構成は、図1及び図2と基本的に同様であるが、実施の形態10では、図1において記憶部36を記憶部106と、制御部42を制御部112と、入力部40を入力部110と、インタフェース部44をインタフェース部114と、データ収集部46をデータ収集部116と、画像処理部48を画像処理部118と、表示部38を表示部108と、表示制御部34を表示制御部104と読み替える。また、図2において、データ処理部32のコンピュータ本体2010をコンピュータ本体10010と、メモリドライブ2020をメモリドライブ10020、ディスクドライブ2030をディスクドライブ10030と、演算装置(CPU)2040を演算装置(CPU)10040と、ディスクドライブ2030をディスクドライブ6030と、メモリドライブ2020をメモリドライブ6020と、バス2050をバス10050と、ROM2060をROM10060と、RAM2070をRAM10070と、不揮発性記憶装置2080を不揮発性記憶装置10080と、通信インタフェース2090を通信インタフェース10090と読み替える。データ処理部62のハードウェアとしては、上述のとおり、特に限定されないが、汎用コンピュータを使用することが可能である。
以下に説明する脳活動訓練装置の駆動方法では、脳機能画像法により脳活動を示す信号を時系列に計測するための脳活動検出装置としては、リアルタイムfMRIを用いるものとして説明する。
図34を参照して、まず、MRI装置10により被訓練者の脳活動が所定時間の間について時系列に測定される。fMRIとして撮像するためにEPI(Echo-planar imaging)撮像が実行される。
続いて、データ処理部102により、撮像された画像がリアルタイムで再構成処理される。
ここで、実施の形態1で説明したとおり、表1に示す機能的結合識別番号(表1でID)1〜12の関心領域間の機能的結合が、「うつ症状」のラベルの判別のために特徴選択により選択された。
ここで、被訓練者においては、上記12個の機能的結合のうち、少なくとも1つの特定の機能結合を「訓練対象の機能的結合」として、さらに選択・抽出する。ここで、特に限定されないが、このような「訓練対象の機能的結合」としては、たとえば、表1で示した機能結合のうち、うつ病ラベルの判別に最も寄与度の大きい「左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の機能的結合」(表1のID=1の結合)を選択するものとする。ここで、「寄与度」とは、「関連重み付け和における重みの絶対値」を意味するものとする。
図32において、説明したように、この結合は、上述したような「抗うつ薬」の投与によって、健常者側とは反対方向に変化する結合のうち、最大の寄与度を有する結合である。
なお、「訓練対象の機能的結合」としては、このID=1の結合のみを選択してもよいし、ID=1の結合に加えて、他の結合を対象として選択してもよい。このとき、ID=1の結合に加えて選択する機能的結合としては、2番目に寄与度の大きなID=2の結合というような選択をしてもよい。あるいは、ID=1の結合に加えて、「抗うつ薬」の投与によって、健常者側とは反対方向に変化する他の結合、たとえば、ID=3の結合を選択してもよい。ID=1の結合に加えて選択する他の機能結合の個数は、1つでもそれ以上でもよい。
さらに、データ処理部102は、「訓練対象の機能的結合」に対応するの関心領域について、機能的結合の所定時間内の時間的な相関を算出する手続きを行う。
すなわち、データ処理部102は、リアルタイムで測定した安静時のfMRIの連続するn個(n:自然数、n≧1)の計測時刻分のfMRIデータ(信号)から、まず、各関心領域の「活動度」を算出する。この活動度は、たとえば、連続するn個の計測時刻の計測値についての平均値とすることができる。ここで、連続するn個の計測時刻からなる期間を「サンプルステップ」と呼ぶことにする。
そして、このような各サンプルステップの活動度について、データ処理部102は、所定の期間の計測ウィンドウについて、以下の数式により、活動度の時間軸上の相関(以下、「時間相関」という)を算出する。なお、計測ウィンドウは、時間軸上で連続する複数のサンプルステップの活動度(たとえば、m個(m:自然数、m≧2))を含むような期間である。
ここで、xは、「訓練対象の機能的結合」のうちの一方の関心領域の「活動度」を意味し、yは、「訓練対象の機能的結合」のうちの他方の関心領域の「活動度」を意味する。また、x(バー)(以下、文字xの頭部に−を付して「平均」を表す記号を、「x(バー)」と呼ぶ。他の文字についても同様である。)とは、上記のような計測ウィンドウ内での活動度xの平均を表す。「y(バー)」についても同様である。上記の式において、Σは、たとえば、計測ウィンドウ内のm個のサンプルステップの活動度について和を取ることを意味する。また、計測ウィンドウは、特に限定されないが、たとえば、10秒〜20秒程度に設定することができる。
ここで、機能結合の時間相関は、符号として、正または負の値を取りうる。
そして、機能結合の時間相関の符号は、たとえば、「健常者」であれば、負であるのに対して、「うつ症状」有りと判別された被験者では、正である、というような状況がありうる。もちろん、その逆の場合もありうるし、あるいは、両者で同一の符号であることもありうる。機能結合の時間相関の符号およびその絶対値の大きさを併せて、「時間相関のパターン」と呼ぶことにする。
そこで、たとえば、「健常者」において、時間相関の符号が負であるような機能的結合を「訓練対象の機能的結合」として選択した場合、被訓練者の訓練対象の機能的結合の時間相関の健常者の時間相関に対する「近似度」を以下の関数F1のようにして算出できる。
i)被訓練者の訓練対象の機能的結合の時間相関の符号が正であれば、一定値V1をとり、前記符号が負であれば、値がV2(>V1≧0)となり、時間相関の絶対値の増加にとともに、値V2が増加するような関数F1。
逆に、「健常者」において、時間相関の符号が正であるような機能的結合を「訓練対象の機能的結合」として選択した場合、被訓練者の訓練対象の機能的結合の時間相関の健常者の時間相関に対する「近似度」を以下の関数F2のようにして算出できる。
ii)被訓練者の訓練対象の機能的結合の時間相関の符号が負であれば、一定値V1をとり、前記符号が正であれば、値がV2(>V1≧0)となり、時間相関の絶対値の増加にとともに、値V2が増加するような関数F2。
なお、「近似度」としては、上記のような例に限られず、符号を含めて、健常者の機能結合の時間相関に近似するパターンであるほど、「近似度」が高くなる関数であれば、他の関数でもよい。
また、特に限定されないが、「訓練対象の機能的結合」が2つ以上選択されている場合は、各々の近似度を平均した値を「近似度」としてもよいし、あるいは、近似度を各結合の寄与度で重み付け平均した値を「近似度」としてもよい。
さらに、データ処理部102は、上記のようにして算出した近似度に応じて、報酬値を算出する。ここで、報酬値(スコアSC)は、目標の脳活動の結合強度に対する近似度に応じて、より近似度が高いほど、高い得点となる。
データ処理部102は、ディスプレイ6に算出されたスコアSCを表示することで、被訓練者2にフィードバックする。なお、フィードバックされる情報としては、スコアの数値そのものでもよいし、あるいは、例えば図35に示すような、前記数値の大きさに応じて、大きさの変化する図形でもよい。その他、被訓練者が、スコアの大きさを認識可能な画像であれば、他の構成を用いてもよい。
以後は、所定の期間について、EPI撮像から、スコアのフィードバックの処理をリアルタイムで繰り返す。
図36は、ニューロフィードバック訓練シーケンスの一例を示す図である。
図36に示すように、まず、第1日において、被訓練者は、事前の安静時の脳活動測定される。
その後、第1日において、ニューロフィードバックによる訓練を行う。
以後、第2日、第3日についても、ニューロフィードバック訓練を行う。
第4日(最終日)には、ニューロフィードバックによる訓練を行った後に、事後の安静時の脳活動の計測を実施する。訓練にかける日数は、これより少なくとも、あるいは、多くともよい。
さらに、被訓練者に対して、所定期間の経過後、たとえば、2か月後において、安静時の脳活動の計測を実施する。
ニューロフィードバック訓練の治療効果は、上記5.に記載の第1の分類器を使用する治療効果の判別装置及び判別方法、又は上記7.に記載の第2の分類器の生成と第2の分類器を使用する治療効果の判別方法によって、判定することができる。
従って、上記5.に記載の第1の分類器を使用する治療効果の判別装置及び判別方法、又は上記7.に記載の第2の分類器の生成と第2の分類器を使用する治療効果の判別方法の説明は、ここに援用される。
9−2.脳活動訓練装置を制御するコンピュータプログラム
実施の形態10の別形態として、上記処理を実行させて、脳活動訓練装置を駆動するコンピュータプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
(実施例)
I.データ収集及び評価方法
1.データ収集対象
広島大学病院及び広島市内のクリニックで募集した105名の患者について、精神障害の診断と統計マニュアル-IVの基準の大うつ病性障害(MDD)の患者を選別するために、M.I.N.I.(精神疾患簡易構造化面接法)を使用してスクリーニングした。現在及び過去に躁病症状、精神病症状、アルコール依存症又は乱用、薬物(substance)依存症又は乱用、反社会性人格障害を除外基準とした。そして、最終的には、自己申告でうつ症状を伴う93名の患者をMDD分類器の訓練データセットとして選択した。患者は、投薬を開始する前に、又は投薬開始から0〜2週間以内に後述する条件でfMRIデータを取得した。145人の健常者は、地域から募集し、M.I.N.I.面接を行い、いかなる精神疾患の既往もない者を選択した。本実施例における実験は、広島大学の倫理委員会の承認を得て行った。また、実験を開始するに先立ち、全ての参加者に対して書面でインフォームドコンセントを行った。
また、完全に独立した4つ施設からもうつ病コホートを収集した。
前記MDDには、メランコリー型MDD、非メランコリー型MDD、治療抵抗型MDDの3つのサブタイプが含まれる。以下、実施例において、全てのMDD患者を含む群を全MDD群、メランコリー型MDD患者を含む群をメランコリー型MDD群、非メランコリー型うつ病患者を含む群を非メランコリー型MDD群、治療抵抗型MDD患者を含む群を治療抵抗型MDD群、健常者の群を健常対照群とする。
メランコリー型MDD分類器を生成するための訓練データセットは、Beck Depression Inventory(BDI)に基づいて、年齢及び性別が健常対照群と一致した、中程度のうつ病症状を伴ううつ病の(M.I.N.I.に基づく)サブタイプとなるよう限定した。各グループに偏りがないようにするために、患者数及び健常者(年齢及び性別が患者と一致し、かつBDI-IIスコアが10未満である)数を等しくした(表2a)。知能指標値(IQ)を評価するためのJapanese Adult Reading Test(JART)の得点については、訓練データセットに3つの欠失データ(メランコリー型MDD群に1つ、全MDD群及びメタンコリー型MDD群の両方に対する健常対照群うち2つ)と、テストデータセットの治療抵抗型MDD群に2つの欠失データがあった。BDIスコアでは、訓練データセットに、全MDD群のみに対する健常対照群に2つの欠失データと抗うつ剤治療後のデータに1つの欠失データがあった。
表2bには、本実施例で使用するテストデータの内訳を示す。
2.独立外部バリデーションコホートに対する汎化
日本放射線医学研究所において独立した検証コホート(千葉のコホート)を作成した。参加者について、M.I.N.I.に基づいて精神疾患の生涯履歴を評価した。
MDD患者には併存症の精神障害はなく、健常対照群には身体的、神経学的、精神医学的障害はなく、現在または過去の薬物乱用の病歴はなかった。すべての参加者に対して実験前に書面によるインフォームドコンセントを実施した。今回の実験は、1964年ヘルシンキ宣言及びその後の改正に規定された倫理基準に従って、放射性医薬品安全委員会および放射線医学研究所の機関レビュー委員会によって承認された。
3.うつ病及びうつ症状の評価
うつ病及びうつ病の症状を評価するためにBeck Depression Inventory(BDI)を使用した。さらに、治療効果の評価には、Hamilton Rating Scale for Depression(HAMD)も使用した。
4.非メランコリー型MDD及び治療抵抗型MDD
非メランコリー型MDD群は、BDIスコア17以上のすべてのMDDを含む。治療抵抗型MDDは、MDDと診断された中で、抗うつ薬を2種類以上投与しても抑うつ症状が改善しない者である。
5.他の精神疾患における分類器の評価
自閉症(ASD)患者および統合失調症(SSD)患者のfMRIデータは、Yahataらの文献(NATURE COMMUNICATIONS | 7:11254 | DOI: 10.1038/ncomms11254)に記載のデータを使用した。うつ病の合併症による影響を最小限に抑えるために、ASD群は、抗うつ薬が有効でなかった者に限定した。fMRIデータは、MDD患者健常者と同様に、目を開いた安静状態で取得した。
6.fMRIデータの取得
fMRIデータの取得の際、被験者に、暗い照明のあるスキャンルームでは、特に何も考えず、眠らずに、モニター画面中央の十字マークを見続けてもらうように依頼した。各施設におけるfMRIデータ取得条件の詳細は表3に示す。
8.fMRIイメージングデータの前処理及び領域間相関
すべてのfMRIデータは、Yahataらの以下の文献に記載されている同一の手法を用いて前処理した。
文献:Noriaki Yahata, Jun Morimoto, Ryuichiro Hashimoto, Giuseppe Lisi, Kazuhisa Shibata,Yuki Kawakubo, Hitoshi Kuwabara, Miho Kuroda, Takashi Yamada, Fukuda Megumi, Hiroshi Imamizu, Jose´ E. Na´nez Sr, Hidehiko Takahashi, Yasumasa Okamoto, Kiyoto Kasai, Nobumasa Kato, Yuka Sasaki, Takeo Watanabe & Mitsuo Kawato, “A small number of abnormal brain connections predicts adult autism spectrum disorder”, NATURE COMMUNICATIONS, DOI: 10.1038/ncomms11254 Matlab R2014a(Mathworks inc., USA)のSPM8(Wellcome Trust Center for Neuroimaging、University College London、UK)を用いてT1強調構造画像及び安静状態機能画像を前処理した。機能イメージは、スライスタイミング補正および平均イメージへのアライメントで前処理された。次に、平均機能画像により同調した構造画像のセグメント化によって得られた正規化パラメータを用いて、fMRIデータを正規化し、2×2×2mm3ボクセルで再サンプリングした。最後に、機能的画像を等方性6mm全幅半値ガウスカーネルで平滑化した。これらの前処理ステップの後、時系列データのフレーム間の相対的な変化に基づいて、スクラビング手順を実行して余分な頭部運動を伴う任意のボリューム(すなわち、機能的画像)を除外した。(頭部運動の要約については表4参照)。
各参加者について、Brainvisa Sulci Atlas(BSA; http:// brainvisa)に解剖学的に定義された大脳皮質全体を覆う137の関心領域(ROI)のそれぞれについて、fMRIデータを経時的に抽出した。本実施例において、施設1において、参加者の小脳の構造及び機能画像が含まれていなかったため、ROIには小脳は含まれていない。バンドパスフィルタ(0.008〜0.1Hz)をかけた後、以下の9つのパラメータ(再配列からの6つの頭部運動パラメータ、白質の時間的ゆらぎ、脳脊髄液の時間的ゆらぎ、脳全体のゆらぎ)を直線的に回帰した。
各参加者の9,316機能的結合(FC)それぞれの行列を得るため、137のROI間のペアワイズピアソン相関を算出した。
II.実施例1.メランコリー型MDD分類のための12対の機能的結合の選択
メランコリー型MDDを分類するための12対の機能的結合を選択するため、表1aに示す66名のメランコリー型MDD患者と66名の健常者のrs-fMRIのデータを使用した。メランコリー型MDD群を分類するための機能的結合の選択は、上述したYahataらの文献で報告した自閉症(ASD)を分類する分類器の作成方法に順じて、実施の形態において説明した手順で行った。
このシステムは、L1正規化されたスパース正準相関分析(L1-SCCA)とスパースロジスティック回帰(SLR)を使用している。SLRは、MDD分類の目的には有用ではないが、客観的に各機能的結合を枝切りしながら、ロジスティック回帰モデルを訓練する能力がある。SLRで訓練する前に、L1-SCCAにより、入力量をある程度減らし、同時に致命的な過適応(過学習)を引き起こす可能性のある攪乱変数(NV:nuisance variable))の影響を減らした。本実施例において、施設、性別、年齢は確率変数に含められているため、これらのうち不要な要素をL1-SCCAにより排除した。この方法では、後述するような内ループと外ループによる入れ子状の(nested)特徴選択とLOOCV(Leave-One-Out Cross Validation)の逐次工程を使用して、情報リークや楽観的すぎる結果を回避した。SLRにより出力されたのは54対の機能的結合であった(図26)。さらに、LOOCVを行い、最終的に12対の機能的結合を選択した(表5及び図26)。
選択された12個の機能的結合FCは、はたして、外ループの手続きの全体にわたって、大きな重みで、頻繁に選択されたかものかどうかについて検討した。これは、最後に選択された12個の機能的結合FCの安定性およびロバスト性に関して重要である。
このような検討のために、以下の式
でk番目のFC(k=1,2,…,9316)のための累積的な絶対的な重みを定義する。
ここで、NがLOOCV回数(つまり被験者の数)であり、wi kは、i番目のLOOCVのk番目の機能的結合FCに関連した重みである。
累積的な重みci kがより大きな値であることは、LOOCVの全体にわたって、MDDとHCへの分類に対してk番目の機能的結合FCが、より重要な貢献をすることを意味する。
表2で、Lat.の“L”と“R”は、左脳と右脳の区別を示す。BSAは、ブロードマン領野を示し、BAは、ブロードマンの領野の番号を示す。rControlは、健常者対照群に対する相関係数を示す。rMDDは、メランコリー型MDD群に対する相関係数を示す。Weightは、関連重み付け和の重みを示す。
図25に選択された12対の機能的結合とその重みを示す。
特に重みが大きかった2つの機能的結合は、left dorsolateral prefrontal cortex (DLPFC, BA46) - left posterior cingulate cortex (PCC) / Precuneusと、left inferior frontal gyrus (IFG opecular, BA44) - right DLPFC (BA9) / frontal eye field (FEF, BA8) /supplementary motor area (SMA, BA6) であった。これらの機能的結合は、MDDの反復経頭蓋磁気刺激(rTMS)治療において標的とされるLeft DLPFC/IFGと重複していた。
このようにして、メランコリー型MDD患者と健常者を分類するための回帰式、すなわち分類器を作成した。さらに、各被験者がメランコリー型MDDか健常であるかを判断するための指標として関連重み付け和(the associated weighted linear sum:WLS)を算出した。
メランコリー型MDD群と健常対照群のWLSの分布を図27aに示す。黒棒がメランコリー型MDD群であり、白棒が健常対照群である。正確性は70 %であった(感度64 %, 特異度77 % AUC 0.77; 並べ替え検定においてp=.049:図30a)。この結果は、選択された12対の機能的結合から構成される分類器(以下、「本発明の第1の分類器」とする)は、メランコリー型MDDを健常者と識別できることを示している。
また、千葉において集められたコホート(慶応大学病院にて収集した。11名のメランコリー型MDDと40名の健常者を含む)についても同様に上記方法で取得した分類器を使って検討した。その結果を図27bに示す。千葉から集めたコホートの正確性は65%(感度64 %,特異度65 %, 及びAUC 0.62; 並べ替え検定においてp=0.036(図30b))であった。本発明の分類器は、訓練データからは独立したコホートについてもメランコリー型MDDと健常者を分類可能であった。
以上の結果から本発明の分類器は、メランコリー型MDDを分類可能であり、かつ汎化されていると考えられた。
III.実施例2.本発明の分類器の非メランコリー型MDD群及び全MDD群への適用 図27cには、全MDD群、メランコリー型MDD群、及び非メランコリー型MDD群の各群を訓練データ(縦方向)として本発明の分類器を作成し、それぞれの群をテストデータ(横方向)として正確性を検討したデータを示す。LOOCVを行った結果、たとえば、メランコリー型MDD群で生成した分類器は、非メランコリー型MDDにおける正確性は、54%(感度42 %、特異度67 %、AUC 0.65)であった。全MDD群については、LOOCVにおける正確性は66 % (感度58 %、特異度74 %、AUC 0.74)であった。
この結果は、全MDD群を訓練データとする分類器は、メランコリー型MDD群を訓練データとした分類器やメランコリー型MDD群を訓練データとした分類器と比較して分類精度が悪くなることを示している。そして、全MDD群が、より多様なMDDのサブタイプを含んでいると考えられた。
図28d、e及び図29f、g、hには、メランコリー型MDD群(d)、非メランコリー型MDD群(e)、治療抵抗型MDD群(f)、ASD群(g)、SSD群(h)の各被験者について求められたWLSの平滑化ヒストグラムとAUC値を示す。有意差は、Benjamini-Hochberg修正Kolmogorov-Smirnov検定により算出した。各群のWLSデータは、健常対照群の中央値と標準偏差が一致するように標準化した。なお、この標準化は定量分析には適用されない。メランコリー型MDD群−健常対照群(d)についてはAUC 0.77、p=1.5×10-5であった。非メランコリー型MDD群−健常対照群(e)では、AUC 0.65、p=.051であった。治療抵抗型MDD群−健常対照群(f)ではAUC 0.46、p=0.54であった。自閉症群−健常者対象群(g) では、AUC 0.51、p=0.74であった。統合失調症群−健常者群(h)ではAUC 0.43、p=0.038であった。
IV. 実施例3.WLSスコアを用いたうつ病の重症度の評価と治療効果の評価
WLSスコアがうつ病の重症度と相関するか否かを調べるため、BDスコアとWLSスコアの相関を検証した。図31aに示すように、全MDD群と健常対照群を合わせた群についてBDIスコアとWLSスコアの相関をとったところ、n=186 r = 0.655, 並び替え検定p = 0.001となり、両者は相関していることが示された。また、全MDD群のみをみても(図31b)、n=93でr=0.188、並び替え検定p =0.046となり、比較的相関していた。
次に、選択的セロトニン再取り込み阻害薬(SSRI)であるエスシタロプラムの効果をWLSスコアで評価できるか検討した。エスシタロプラムを6〜8週間投与した、表1aに示す寛解した24名のメランコリー型MDD患者各々のWLSスコアの平滑化ヒストグラムを図31cに示す。抗うつ剤を投与するとWLSスコアの分布は健常対照群よりにシフトした。AUC 0.72、Benjamini-Hochberg修正Kolmogorov-Smirnov検定p=0.008であった。
ただし、この結果は、集団についての結果であり、WLSと個人の改善との関係を示すものではない。そこで、ΔWLS(具体的には WLSpost - WLSpre)とΔBDI(具体的には BDIpost - BDIpre)の関係、ΔWLSとΔHAMD (具体的にはHAMDpost - HAMDpre)を検討した。有意差は並び替え検定で求めた。図31dに示すように、ΔBDIとΔWLSは有意に相関していた(r=0.373、p=0.040)が、ΔHAMDとΔWLSは相関していなかった(r=0.154、p=0.237)。BDIとHAMDが乖離することは既に報告されている。
実施例5:12対の機能的結合それぞれの治療前後での変化
次に12対の機能的結合それぞれに対するSSRI(エスシタロプラム)の影響について検討するため、個々の機能的結合についてΔWLSに対する寄与を検討した。
具体的には、投与前後解析において、式:
(式中、Nは治療を受け患者数を示す。wiは、分類器の重みを示す。)
を用いて各機能的結合の寄与スコアを算出した。
さらに、参考点を得るために、同様の寄与スコアをMDDと健常対照群の間でも式:
(式中、MはMDDの数を示し、Nは健常対照の数を示す。)
によって寄与スコアを算出した。
すなわち、siは分類器の重みで重み付けされた健常対照群とMDD群の平均結合強度の差を示す。分類器において、MDDは陽性(正のスコア)を示し、健常者は陰性(負のスコア)示すため、siは常に陰性(負のスコア)となる。また、式2において
は減数 (減算において負の符号)となる。ciとsiが大きな負の値をとることは治療後−治療前健常者−MDDにおいてFCiの寄与が大きいことを示す。言い換えると、ciが負になることは 治療後のFCiが健常者の対応するFCにより近づいたことを示す。これに対して、ciが正になることは治療後のFCiがMDDとなったことを示す。健常対照群とMDD群で人数が異なるため、それぞれのFCiについてのciとsiの有意差は、ウェルチのt-検定で評価した。ウェルチのt-検定は、下式:
にしたがって、算出した。
各検定のp値は、Benjamini-Hochberg法によって複数の比較に批准させた。図32aに示すように、FC1 (i.e. DLPFC-Precuneus/PCC)は、ci及びsiの差が最も大きく顕著であった。これは、抗うつ剤の投与後に FC1が健常対照群とは反対のMDD群の方に推移したことを示す。
メランコリー型MDD分類のための12対の機能的結合(表Cに記載のこの機能的結合識別番号を「FC♯」と標記することがある)の中で8つの機能的結合については、健常対照群の方向へ向かって推移し、いわゆる「正常化」した。しかし、FC# 1、3、9、12の4対の機能的結合については、健常対照群とは逆の方向に向かって推移した(図32a)。この中で、FC# 1(Left DLPFC - Left Precuneus / PCC FC)は、12対の中で最も寄与度が高い機能的結合である。メランコリー型MDD群と健常対照群の分類においてFC#1の治療前後の変化の有意差は、ペアワイズ比較でp=0.009であった(図32a)。
最も寄与度が高い、2つの主要な機能的結合(FC#1, FC#2)に着目した結果、治療後に有意な差が観察されたのに対して(対応のあるt検定, p = .002)、治療前での差異は観察されなかった(p= 0.96)。二要因分散分析による検定では、FC#1・FC#2の機能的結合と、治療前後の時間の交互作用が有意であり(FC #number x time, F(1,23)=6.70, p =0.016)、各主効果に関しては有意な差は見られなかった。しかし、治療後のFC# 1とFC#2との間でのペアワイズt-検定では、推移に有意差が認められた(p=0.002)。ここで、図32、33の全ての図において、正方向への推移はメランコリー型MDDの方向へ向かっていることを示し、負方向への推移は健常者方向へ向かっていることを示す。これは、メランコリー型MDD分類器の重みの符号に、FCまたはFCの差を乗じているためである。上述のように、エスシタロプラム治療後にFC# 1で正方向へ推移したため、完全に独立したコホートでも同じ結果となるか検証した。千葉において収集された11名のメランコリー型MDDの患者(セロトニン-ノルエピネフリン再取り込み阻害剤(SNRI)デュロキセチンを6〜8週間投与)について検討した(図32c)。治療前には有意差がなかった(n.s.、p = 0.77)のに対して、治療後の対t-検定では、FC# 1とFC#2のトレンドレベルの差が認められた(p=0.10)。次に、健常対象者(19名)にSSRIを投与した際のFC# 1とFC#2の変化を観察した。パロキセチン投与の単回投与後にFC#1とFC#2との間に有意差が認められた(p=0.033 )。治療前には差はなく(p = 0.18)、有意な関連は認められなかった(図32d)。
さらに、広島コホート内の、臨床医によって治療効果を判定された患者について、非寛解患者(n=7)と寛解患者(n=7)に分け、抗うつ治療の効果を調べた。その結果、この結果、寛解・非寛解の群間 x FC#1・FC#2の機能的結合間 x 治療前後の時間に有意な交互作用が見られ (Fig. 33e, F(1,22) = 8.86, p=.007)、FC#1とFC#2における治療反応のパターンに差があることが示された。具体的には、FC#2の値は、寛解群においてのみ、抗うつ剤の治療後に減少し(p=0.001)、これらの値には有意差が認められた(p=.019)。さらに、非寛解群では、治療前には差がないにも関わらず、抗うつ剤治療後にFC#1とFC#2に有意差が認められた(p=0.033)。
寛解群と非寛解群でFC#1とFC#2がどのように変化するかを検討するため、24名の患者についてFC#1の変化とFC#2の変化を2次元プロットに表した(図33f)。
FC#1及びFC#2は、下式によりもとめた。
FC1postは、FC1の投与後を示し、FC1preは、FC1の投与前を示し、 FC2postは、FC2の投与後を示し、FC2preは、FC2の投与前を示し、signは各機能的結合の重みの符号を示す。この式の精度は、Accuracy:0.75, AUC:0.79, Specificity:0.88, Sensitivity:0.43であった。
その結果、FC#1の変化とFC#2の変化は、寛解に寄与しているように見える。一方で、FC#1の変化とFC#2の変化との間には、有意差や相関は認められなかった(p=0.57)。
V. 実施例4.ニューロフィードバックV-1.FCニューロフィードバック訓練に使用する標的機能的結合の選択
以下の工程に従って、ニューロフィードバックの標的として、前記表5に示す第1の機能的結合を選択した。
はじめに、rs-fcMRIのバイオマーカーとして、本発明の分類器とBDIスコアに基づくうつ症状の重症度を予測するためのバイオマーカー(Yamashita A, Hayasaka S, Lisi G, Ichikawa N, Takamura M, Okada G, Morimoto J, Yahata N, Okamoto Y, Kawato M, Imamizu H(2015) Common functional connectivity between depression and depressed mood. The 37th annual meeting of theJapanese society of biological psychiatry; 24-26 September; Tokyo, Japan.)を構築した。
標的機能的結合には、上記のバイオマーカーを含むものとした。
ニューロフィードバック訓練の間、参加者は、標的FCの相関を減少させることを目的とした。
訓練には、MDDを患う3名と、亜臨床的うつ病の7名の計10名が参加した。ニューロフィードバック訓練に先立って、行われた2回の異なる時点における亜臨床的うつ病の参加者のBDI-IIの平均値は>10であった。
ニューロフィードバック訓練は、参加者にfMRI装置に仰向けになってもらい、プリズム眼鏡を介して装置内のモニターを見せ、図35に示す機能的結合の結合指標が負の方向に動くようにイメージさせることにより行った。
ニューロフィードバックスコアは、公知の方法に従って算出した。
図37aにMDD参加者3名の各訓練日のニューロフィードバックスコアの集計を示す。スコアは、訓練の4日間にわたって増加傾向を示した。結果の有意性を、2つの説明変数(訓練日それぞれと、訓練項目)と1つの応答変数(ニューロフィードバックスコア)を含む重回帰モデルを用いて確認したところ、訓練日において有意に良い効果を示した(95% 信頼区間(CI)計数 1.9-9.1)。事後比較として行ったt-検定では、MDD参加者3名全てにおいて、訓練の初日と比較して、最終日で有意に高いニューロフィードバックスコアが得られた。(t = 4.01, P < .001). これらの結果は、MDDの参加者全てが、訓練を通して標的FCの負の相関を誘導することを学んだこと示している。さらに、図37bに示す用に、MDD参加者3名全てにおいて、訓練後に、うつ症状の重症度を示すHamilton Rating Scale for Depression(HAMD)が減少した。
MDD参加者と同様に、亜臨床的うつ病の参加者7名においても、訓練期間中、ニューロフィードバックスコアが増加する傾向にあった(図37c)。この結果について有意性を解析したところ、一元配置分散分析(One-way ANOVA)で訓練日の主効果が有意であった(p = 0.0046)。事後比較ペアt検定では、訓練の最終日にニューロフィードバックスコアが訓練の初日より有意に高かった。訓練後のBDIスコアも減少傾向を示した(p = .07)。また、図37dに示すように、亜臨床的うつ病の参加者7名中5名については、正常方向に標的rs-fcMRIを変化させ、訓練前後間でのrs-fcMRI変化量は、BDIスコアと有意に相関した(r = 0.87, p = 0.011)。
上述のように、半数以上の参加者が、訓練を通じてBDIスコアと相関する表5に示す機能的結合を減少させたことから、表5に示す機能的結合は、うつ症状を改善するためのニューロフィードバック訓練の標的機能的結合として使用できることが示された。さら、本発明の表5に示す機能的結合は、うつ症状を改善するためのニューロフィードバック訓練法を開発するための治療方法バイオマーカーとして使用できることが示された。
V-2.他の抗うつ治療との比較
ニューロフィードバック訓練法によるMDDの治療効果を、他の治療方法(抗うつ薬投与、反復経頭蓋磁気刺激療法:rTMS、修正型電気けいれん療法:ECT intervention)と比較した。ニューロフィードバック訓練法は、他のうつ治療よりも治療効果が改善していた。
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 ネットワークインタフェース。

Claims (51)

  1. 被検者におけるうつ症状の有無の判定を補助するための判別装置であって、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備え、
    前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    演算装置をさらに備え、前記演算装置は、
    前記分類器により、前記被検者のうつ症状に対する分類結果を生成する判別処理を実行する、判別装置。
  2. 前記分類器生成処理において、前記機能的結合の重み付け和は、前記うつ症状の疾患ラベルに特異的に関連するとしてスパース正準相関分析により抽出された機能的結合から、スパースロジスティック回帰による特徴選択により算出される、請求項1に記載の判別装置。
  3. 前記選択された複数の機能的結合は、第1の機能的結合と第2の機能的結合の両方を含み、
    前記分類器は、各機能的結合の相関行列の要素の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される、請求項1または2に記載の判別装置。
  4. 前記第1および第2の機能的結合が、下記表1に示す機能的結合識別番号1および2にそれぞれ対応するとき、
    前記選択された複数の機能的結合は、下記表1に示す機能的結合識別番号3〜12から選択される少なくとも一つの機能的結合をさらに含む、請求項3に記載の判別装置:
  5. 前記うつ症状は、メランコリー型うつ病に伴う症状である、請求項1〜4のいずれか一項に記載の判別装置。
  6. うつ症状のレベルの判別装置であって、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備え、
    前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    演算装置をさらに備え、前記演算装置は、
    前記分類器により、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成し、
    前記指標値をあらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較し、
    前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する、ように構成される、判別装置。
  7. 被検者における治療効果の判定するための判別装置であって、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備え、
    前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    演算装置をさらに備え、前記演算装置は、
    前記分類器により、第1の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第1の値を生成し、
    前記分類器により、治療開始後であって第1の時点より後の第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素についてうつ症状を評価する第2の値を生成し、
    前記第1の値と前記第2の値とを比較し、
    i)前記第2の値が前記第1の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定する、及び/又は、ii)前記第2の値が前記第1の値よりも改善していない場合に、前記治療が被検者のうつ症状の改善に有効でないと決定する、ように構成される、判別装置。
  8. うつ病患者を層別化するための判別装置であって、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備え、
    前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    演算装置をさらに備え、前記演算装置は、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、
    安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成し、
    前記指標値を、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較し、
    前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する、ように構成される、判別装置。
  9. 請求項1〜8のいずれか一項に記載の判別装置の機能をコンピュータに実行させるためのプログラム。
  10. 安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を備え、前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含み、
    前記指標値が基準値を超えた場合に、前記被検者がうつ症状を有すると決定する工程を備える、うつ症状を有する被検者を判定することを補助する判別方法。
  11. 前記指標値は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により算出された重み付け和に基づいて生成され、前記分類器生成処理においては、前記所定領域間の機能的結合の重み付け和は、機械学習によりうつ症状の疾患ラベルに関連する特徴を選択する際に算出される、請求項10に記載の判別方法。
  12. 前記分類器生成処理においては、前記機能的結合の重み付け和は、前記うつ症状の疾患ラベルに特異的に関連するとしてスパース正準相関分析により抽出された機能的結合から、スパースロジスティック回帰による特徴選択により算出される、請求項11に記載の判別方法。
  13. 前記機能的結合は、第1の機能的結合と第2の機能的結合の両方を含み、
    前記指標値は、各機能的結合の相関係数の重み付け和である、請求項10〜12のいずれか一項に記載の判別方法。
  14. 前記第1及び第2の機能的結合が、下記表2に示す機能的結合識別番号1および2にそれぞれ対応するとき、前記機能的結合は、下記表2に示す機能的結合識別番号3〜12から選択される少なくとも一つの機能的結合をさらに含む、請求項13に記載の判別方法:
  15. 前記うつ症状は、メランコリー型うつ病に伴う症状である、請求項10〜14のいずれか一項に記載の判別方法。
  16. 前記指標値は、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器であって、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成された分類器により生成される、請求項10〜15のいずれか一項に記載の判別方法。
  17. 安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を備え、機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含み、
    前記指標値を、あらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較する工程と、
    前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する工程と、を備える、被検者においてうつ症状のレベルを判定することを補助するため判別方法。
  18. 第1の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第1の値を生成する工程を備え、前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含み、
    治療開始後であって第1の時点より後の第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素についてうつ症状を評価する第2の値を生成する工程と、
    前記第1の値と前記第2の値とを比較する工程と、
    前記第2の値が前記第1の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定する工程と、
    を含む、被検者における治療効果の判定を補助するための判別方法。
  19. 前記決定工程の後に、前記治療を継続することを提示する工程をさらに含む、請求項18に記載の判別方法。
  20. 前記治療が、投薬療法、ニューロフィードバック療法、修正型電気けいれん療法および反復経頭蓋磁気刺激療法から選択される少なくとも一種である、請求項18または19に記載の判別方法。
  21. 前記治療方法が投薬療法であって、前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合を含む、請求項18及び19に記載の判別方法。
  22. 第1の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第1の値を生成する工程を備え、前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含み、
    治療開始後であって第1の時点より後の第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素ついてうつ症状を評価する第2の値を生成する工程と、
    前記第1の値と前記第2の値とを比較する工程と、
    前記第2の値が前記第1の値よりも改善していない場合に、前記治療が被検者のうつ症状の改善に有効でないと決定する工程と、
    を含む、被検者における治療効果の判定を補助するための判別方法。
  23. 前記決定工程の後に、前記治療の終了を提示する工程を含む、請求項22に記載の判別方法。
  24. 前記決定工程の後に、前記治療を他の治療へ変更することを提示する工程を含む、請求項22に記載の判別方法。
  25. 前記治療が、投薬療法、ニューロフィードバック療法、修正型電気けいれん療法および反復経頭蓋磁気刺激療法から選択される少なくとも一種である、請求項22〜24のいずれか一項に記載の判別方法。
  26. 前記治療方法が投薬療法であって、前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合を含む、請求項22〜24のいずれか一項に記載の判別方法。
  27. うつ病が予め設定された複数のサブクラスに層別化されている場合に、
    安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を備え、前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含み、
    前記指標値を、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較する工程と、
    前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する工程と、を含む、うつ病患者の層別化を補助するための判別方法。
  28. うつ病患者の層別化を補助するための判別装置であって、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置と、層別化処理を実行する演算装置とを備え、
    うつ症状が予め設定された複数のサブクラスに層別化されている場合に、前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ病のサブクラスのラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記選択された機能的結合に対応する判別対象の相関行列の要素に対して前記うつ病のサブクラスのラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    前記演算装置は、前記重み付け和と前記判別対象の相関行列の要素から前記サブクラスを判別するための処理を実行する、判別装置。
  29. うつ病患者における治療効果の判定を補助するための判別装置であって、
    演算装置を備え、前記演算装置は、
    複数の対象者について、下記表3に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する分類器生成処理を実行し、
    第1の時点の被検者の安静時の脳内における、前記複数の機能的結合の第1の相関を計測し、
    第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第2相関を計測し、
    前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する、
    ように構成された、判別装置:
  30. 前記複数の機能的結合は、
    機能的結合識別番号1の機能的結合、及び機能的結合識別番号2の機能的結合である、請求項29の判別装置。
  31. 複数の対象者について、下記表4に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する工程と、
    第1の時点の被検者の安静時の脳内における、前記複数の機能的結合の第1の相関を計測する工程と、
    第2の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第2相関を計測する工程と、
    前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する工程と、
    を含む、うつ病患者における治療効果の判定を補助するための判別方法:
  32. 安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、被検者におけるうつ症状の有無の判定またはうつ症状のレベルの判定を補助するために使用する方法であって、
    前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程を含む、使用方法。
  33. 安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、うつ症状の治療の効果の判定を補助するために使用する方法であって、
    前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように前記分類器を生成する工程を備え、、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程をさらに備える、使用方法。
  34. うつ症状が予め設定された複数のサブクラスに層別化されている場合に、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、被検者を前記サブクラスに層別化することを補助するために使用する方法であって、
    前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように前記分類器を生成する工程を備え、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程をさらに備える、使用方法。
  35. 被検者における治療効果の判定を補助するための判別装置であって、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から第1の分類器を生成する第1の演算装置を含む分類器生成装置を備え、
    前記第1の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    前記第1の演算装置は、複数の対象者について、下記表5に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第2の分類器を生成し、
    第2の演算装置と記憶装置を含む層別化装置をさらに備え、
    前記記憶装置は、前記第1および第2の分類器の情報と、うつ症状が予め設定された複数のサブクラスに層別化するための情報とを記憶し、
    前記第2の演算装置は、
    前記第1の分類器を用いて、被検者を前記サブクラスに層別化処理を実行し、
    前記層別化処理によって特定のサブクラスに層別化された被検者について、第1の時点の安静時の前記複数の機能的結合の第1の相関を計測するための処理を実行し、
    第2の時点の前記被検者と同一の被検者の安静時の前記複数の機能的結合の第2相関を計測するための処理を実行し、
    前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記第2の分類器により、前記被検者に対する治療効果を判別する、ように構成される、判別装置:
  36. うつ症状が予め設定された複数のサブクラスに層別化されている場合に、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から第1の分類器生成処理によって生成された第1の分類器を使用して、被検者を前記サブクラスに層別化する層別化工程であって、
    前記第1の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含む、層別化工程と、
    前記層別化工程によって特定のサブクラスに層別化された被検者について、治療開始前の安静時の前記複数の機能的結合の第1の相関を計測する第1の相関計測工程と、
    治療開始から所定期間経過後の前記被検者と同一の被検者の安静時の前記複数の機能的結合の第2相関を計測する第2の相関計測工程と、
    特定のサブクラスに分類された前記被検者に対する治療効果を第2の分類器を用いて判別する判別工程であって、
    前記第2の分類器は、複数の対象者について、下記表6に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第2の分類器生成処理によって予め生成され、
    前記被検者の前記複数の機能的結合の前記第2の相関と前記第1の相関との差に基づき、前記被検者に対する治療効果を判別する、判別工程と、
    を備える、被検者における治療効果の判定を補助するための判別方法:
  37. さらに、ドラッグリプロファイリングのために使用される、請求項18〜21、または請求項22〜26に記載の判別方法。
  38. さらに、ドラッグリプロファイリングのために使用される、請求項29または30に記載の判別装置。
  39. ドラッグリプロファイリングのために使用される、請求項31に記載の判別方法。
  40. ドラッグリプロファイリングのために使用される、請求項35に記載の判別装置。
  41. ドラッグリプロファイリングのために使用される、請求項36に記載の判別方法。
  42. 演算装置と、記憶装置とを備えた第1の分類器生成装置であって、
    前記演算装置は、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器を特定する情報を生成し、
    前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記選択された複数の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および
    左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも1つを含み、
    前記記憶装置は、前記演算装置によって生成された前記第1の分類器を特定する情報を記憶する、
    第1の分類器生成装置。
  43. 安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により第1の分類器を生成する方法であって、前記第1の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成され、
    前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第2の機能的結合から選択される少なくとも一つを含む、
    第1の分類器の生成方法。
  44. 演算装置と、記憶装置とを備えた第2の分類器生成装置であって、
    前記演算装置は、
    複数の対象者について、下記表3に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第2の分類器を生成し、
    前記記憶部は、前記演算装置によって生成された前記第2の分類器を特定するための情報を記憶する、第2の分類器生成装置:
  45. 複数の対象者について、下記表4に示す機能的結合識別番号1〜12から選択される複数の機能的結合の第1の時点の相関と治療開始後であって第1の時点よりも後の第2の時点の相関をそれぞれ計測して、前記複数の機能的結合の第1の時点と第2の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する工程を備える、第2の分類器生成方法:
  46. ニューロフィードバック訓練を実行するための脳活動訓練装置であって、
    安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を時系列で予め測定した信号に基づいて、前記複数の所定領域間の機能的結合から機械学習により、うつ症状の疾患ラベルを判別するために特徴選択により選択された複数の機能的結合のうち訓練対象の機能的結合を特定する情報を格納する記憶装置を備え、
    前記訓練対象の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合を含み、
    前記記憶装置は、ニューロフィードバック訓練における前記訓練対象の機能的結合の目標パターンを記憶し、
    前記ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知するための脳活動検知装置をさらに備え、前記被訓練者の脳内の複数の所定領域のそれぞれは、前記複数の参加者における脳内の複数の所定領域にそれぞれ対応し、
    提示装置と、
    演算装置とをさらに備え、前記演算装置は、
    i)前記脳活動検知装置により検知された信号から、前記訓練対象の機能的結合の時間相関を所定期間について算出し、
    ii)算出された前記時間相関に基づいて、前記目標パターンとの近似度に応じて、報酬値を算出し、
    iii)前記報酬値の大きさを示す情報を前記提示装置により前記被訓練者に対して提示する、ように構成される、脳活動訓練装置。
  47. 前記分類器生成処理において、前記機能的結合の重み付け和は、前記うつ症状の疾患ラベルに特異的に関連するとしてスパース正準相関分析により抽出された機能的結合から、スパースロジスティック回帰による特徴選択により算出される、請求項46に記載の脳活動訓練装置。
  48. 前記うつ症状は、メランコリー型うつ病に伴う症状である、請求項46または47のいずれか一項に記載の脳活動訓練装置。
  49. 前記被訓練者が、
    請求項1〜5に記載の判別装置、又は請求項10〜16に記載の判別方法によりうつ症状があると判別された被検者であるか、
    請求項6に記載の判別装置、又は請求項17に記載の判別方法によりうつ症状のレベルを決定された被検者であるか、
    請求項8に記載の判別装置、又は請求項27の判別方法により層別化されたうつ疾患患者である、
    請求項46〜48のいずれかに記載の脳活動訓練装置。
  50. 請求項46〜49のいずれか一項に記載の脳活動訓練装置を制御するコンピュータプログラム。
  51. ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知するための脳活動検知装置と、演算装置と、記憶装置と、提示装置と、を備える脳活動訓練装置の制御方法であって、
    前記記憶装置に、i)安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を時系列で予め測定した信号に基づいて、前記複数の所定領域間の機能的結合から機械学習により、うつ症状の疾患ラベルを判別するために特徴選択により選択された複数の機能的結合のうち訓練対象の機能的結合を特定する情報と、ii)ニューロフィードバック訓練における前記訓練対象の機能的結合の目標パターンとを格納する工程を備え、
    前記訓練対象の機能的結合は、
    左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第1の機能的結合を含み、 前記脳活動検知装置により、前記ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知する工程をさらに備え、前記被訓練者の脳内の複数の所定領域のそれぞれは、前記複数の参加者における脳内の複数の所定領域にそれぞれ対応し、
    前記演算装置が、前記脳活動検知装置により検知された信号から、前記訓練対象の機能的結合の時間相関を所定期間について算出する工程と、
    前記演算装置が、算出された前記時間相関に基づいて、前記目標パターンとの近似度に応じて、報酬値を算出する工程と、
    前記演算装置が、前記報酬値の大きさを示す情報を前記提示装置により前記被訓練者に対して提示する工程と、をさらに備える、脳活動訓練装置の制御方法。
JP2018005110A 2017-10-03 2018-01-16 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置 Active JP6722845B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
PCT/JP2018/036952 WO2019069955A1 (ja) 2017-10-03 2018-10-02 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
US16/753,291 US12020427B2 (en) 2017-10-03 2018-10-02 Differentiation device, differentiation method for depression symptoms, determination method for level of depression symptoms, stratification method for depression patients, determination method for effects of treatment of depression symptoms, and brain activity training device
EP18864859.6A EP3692909A4 (en) 2017-10-03 2018-10-02 DIFFERENTIATION DEVICE, DEPRESSION SYMPTOMS DIFFERENTIATION PROCESS, DEPRESSION SYMPTOM LEVEL DETERMINATION METHOD, DEPRESSION PATIENT STRATIFICATION METHOD, DEPRESSION SYMPTOM TREATMENT EFFECTS DETERMINATION AND DEPRESSION SYMPTOMS TREATMENT PROCESS, DETERMINATION DEVICE '' BRAIN ACTIVITY
CN201880078144.5A CN111447874A (zh) 2017-10-03 2018-10-02 判别装置、抑郁症状的判别方法、抑郁症状的水平的判断方法、抑郁症患者的分类方法、抑郁症状的治疗效果的判断方法和脑活动训练装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017193897 2017-10-03
JP2017193897 2017-10-03

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2019061843A Division JP7140348B2 (ja) 2017-10-03 2019-03-27 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置

Publications (3)

Publication Number Publication Date
JP2019063478A true JP2019063478A (ja) 2019-04-25
JP2019063478A5 JP2019063478A5 (ja) 2019-10-10
JP6722845B2 JP6722845B2 (ja) 2020-07-15

Family

ID=66337446

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2018005110A Active JP6722845B2 (ja) 2017-10-03 2018-01-16 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
JP2019061843A Active JP7140348B2 (ja) 2017-10-03 2019-03-27 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2019061843A Active JP7140348B2 (ja) 2017-10-03 2019-03-27 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置

Country Status (4)

Country Link
US (1) US12020427B2 (ja)
EP (1) EP3692909A4 (ja)
JP (2) JP6722845B2 (ja)
CN (1) CN111447874A (ja)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019172245A1 (ja) * 2018-03-09 2019-09-12 株式会社国際電気通信基礎技術研究所 脳活動訓練装置、脳活動訓練方法および脳活動訓練プログラム
WO2021205996A1 (ja) * 2020-04-06 2021-10-14 株式会社国際電気通信基礎技術研究所 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラム、脳活動マーカー分類システムおよび脳機能結合相関値のクラスタリング分類器モデル
WO2022014682A1 (ja) * 2020-07-17 2022-01-20 株式会社国際電気通信基礎技術研究所 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム
CN116469541A (zh) * 2023-03-24 2023-07-21 深圳市第二人民医院(深圳市转化医学研究院) 一种抑郁症标志物及其在抑郁症预后中的应用、评估装置
JP2023143663A (ja) * 2022-03-24 2023-10-06 之江実験室 磁気共鳴及びシャムグラフニューラルネットワークに基づく脳アトラスの個別化方法及びシステム
US12020427B2 (en) 2017-10-03 2024-06-25 Advanced Telecommunications Research Institute International Differentiation device, differentiation method for depression symptoms, determination method for level of depression symptoms, stratification method for depression patients, determination method for effects of treatment of depression symptoms, and brain activity training device
KR102710036B1 (ko) * 2024-01-31 2024-09-25 이화여자대학교 산학협력단 산후우울증 발생 위험도를 예측하는 전자 장치 및 이의 동작 방법
US12127826B2 (en) 2022-03-24 2024-10-29 Zhejiang Lab Brain atlas individualization method and system based on magnetic resonance and twin graph neural network

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20240135558A1 (en) * 2014-04-25 2024-04-25 Thornhill Scientific Inc. Imaging abnormalities in vascular response
CA3123192A1 (en) * 2018-12-14 2020-06-18 Keio University Device and method for inferring depressive state and program for same
JP7147623B2 (ja) * 2019-02-21 2022-10-05 株式会社島津製作所 脳血流量の特徴量の抽出方法
CA3134521A1 (en) * 2019-03-22 2020-10-01 Cognoa, Inc. Personalized digital therapy methods and devices
US11412933B2 (en) * 2019-07-19 2022-08-16 Brainsight Technology Private Limited System for analysing an activity of brain using magnetic resonance imaging (MRI) data
CN111588351B (zh) * 2020-06-09 2023-04-11 天津大学 一种基于功能近红外脑网络的中医针刺得气智能检测系统
WO2022056652A1 (zh) * 2020-09-15 2022-03-24 洪硕宏 评估跨颅磁刺激对忧郁症患者是否有效的辅助判断装置
CN112370659B (zh) * 2020-11-10 2023-03-14 四川大学华西医院 基于机器学习的头部刺激训练装置的实现方法
US11948303B2 (en) * 2021-02-26 2024-04-02 Capsovision Inc. Method and apparatus for objective assessment of gastrointestinal conditions based on images captured in the GI tract
CN113052800B (zh) * 2021-03-09 2022-02-22 山东大学 一种阿尔茨海默症影像分析方法与装置
CN112991335B (zh) * 2021-04-23 2024-04-26 合肥中聚源智能科技有限公司 一种想像能力评估方法、系统及可读存储介质
CN113331131B (zh) * 2021-05-28 2022-10-14 南方医科大学 一种可卡因成瘾戒断后焦虑抑郁样动物模型的构建及其应用
US11379981B1 (en) * 2021-06-07 2022-07-05 The Florida International University Board Of Trustees Systems and methods for diagnosing autism spectrum disorder using fMRI data
CN113505721A (zh) * 2021-07-22 2021-10-15 北京脑陆科技有限公司 基于两级分类器的抑郁识别方法、装置、终端及介质
CN113947580B (zh) * 2021-10-18 2022-07-19 上海市精神卫生中心(上海市心理咨询培训中心) 针对精神疾病实现神经影像辅助诊断处理的系统
CN114860922B (zh) * 2022-03-25 2023-04-18 南京脑科医院 心理测评量表的分类模型获取方法、筛查方法及系统
CN115147417B (zh) 2022-09-02 2022-11-15 之江实验室 一种基于过滤法特征选择的功能连接矩阵处理系统及装置
CN116543042B (zh) * 2023-05-10 2024-01-23 中国科学院心理研究所 基于组水平平均统计图的抑郁症tms个体化靶点定位方法及系统
CN116504363B (zh) * 2023-05-10 2024-03-08 中国科学院心理研究所 基于组水平差异统计图的抑郁症tms个体化靶点定位方法及系统

Family Cites Families (55)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020103429A1 (en) 2001-01-30 2002-08-01 Decharms R. Christopher Methods for physiological monitoring, training, exercise and regulation
US20040116798A1 (en) 2002-12-11 2004-06-17 Robert Cancro Method and system for investigation of central nervous system drugs using 3-D brain source localization
JP4555230B2 (ja) 2003-09-11 2010-09-29 株式会社日立メディコ 生体光計測装置
US20050215884A1 (en) 2004-02-27 2005-09-29 Greicius Michael D Evaluation of Alzheimer's disease using an independent component analysis of an individual's resting-state functional MRI
JP2008531158A (ja) * 2005-03-04 2008-08-14 メンティス キュラ イーエイチエフ. 神経学的状態を評価する方法および系
US7577472B2 (en) 2005-03-18 2009-08-18 The Mcw Research Foundation Inc. MRI method for producing an index indicative of brain disorders
WO2006102370A2 (en) * 2005-03-22 2006-09-28 University Of Massachusetts Functional brain mri mapping as a marker in cns diseases and disorders
WO2006132313A1 (ja) 2005-06-09 2006-12-14 Hitachi Medical Corporation 疾患判定支援システム
US8597614B2 (en) 2005-08-03 2013-12-03 Centre For Addiction And Mental Health Diagnosis and treatment of mood disorders
US9101276B2 (en) 2006-07-06 2015-08-11 Regents Of The University Of Minnesota Analysis of brain patterns using temporal measures
JP4978860B2 (ja) 2007-01-24 2012-07-18 株式会社国際電気通信基礎技術研究所 行動予測方法及び行動予測装置
CN101681515A (zh) * 2007-06-21 2010-03-24 皇家飞利浦电子股份有限公司 基于模型的痴呆的辨别诊断和显著性水平的交互设置
US9119549B2 (en) 2007-11-12 2015-09-01 Siemens Aktiengesellschaft Method for developing test for neuropsychiatric disease
JP5663314B2 (ja) * 2008-03-04 2015-02-04 リッジ ダイアグノスティックス,インコーポレイテッド 多数のバイオマーカーパネルに基づくうつ障害の診断およびモニタリング
EP2272044A4 (en) * 2008-03-12 2011-07-06 Ridge Diagnostics Inc INFLAMMATION BIOMARKERS FOR MONITORING DEPRESSION DISORDERS
FR2943156B1 (fr) * 2009-03-12 2014-04-11 Commissariat Energie Atomique Procede d'elaboration d'un dispositif de prediction, utilisation, support de stockage d'information et appareil correspondants.
US20100249573A1 (en) 2009-03-30 2010-09-30 Marks Donald H Brain function decoding process and system
JP5569714B2 (ja) 2009-06-16 2014-08-13 株式会社国際電気通信基礎技術研究所 磁気共鳴画像化装置および磁気共鳴画像化方法
WO2011038124A2 (en) 2009-09-24 2011-03-31 President And Fellows Of Harvard College Systems and methods for treating depression, including the use of neurobehavioral tasks
WO2011115956A1 (en) 2010-03-15 2011-09-22 Mcw Research Foundation, Inc. Systems and methods for detection and prediction of brain disorders based on neural network interaction
EP2642919A4 (en) 2010-11-26 2014-10-08 Brigham & Womens Hospital METHOD FOR EVALUATING REPETITIVE HEAD INJURY LESIONS USING TWO-DIMENSIONAL MAGNETIC RESONANCE SPECTROSCOPY
US10238329B2 (en) 2010-12-17 2019-03-26 The Trustees Of Columbia University In The City Of New York Apparatus, method and computer-accessible medium for diagnosing and subtyping psychiatric diseases
US8451930B2 (en) 2011-02-08 2013-05-28 Mitsubishi Electric Research Laboratories, Inc. Sparse channel estimation using order extension for frequency and time selective MIMO wireless communication networks
JP5959016B2 (ja) 2011-05-31 2016-08-02 国立大学法人 名古屋工業大学 認知機能障害判別装置、認知機能障害判別システム、およびプログラム
US8861815B2 (en) 2011-08-03 2014-10-14 International Business Machines Corporation Systems and methods for modeling and processing functional magnetic resonance image data using full-brain vector auto-regressive model
US20130102918A1 (en) 2011-08-16 2013-04-25 Amit Etkin System and method for diagnosing and treating psychiatric disorders
JP5320543B2 (ja) 2011-11-08 2013-10-23 株式会社国際電気通信基礎技術研究所 脳機能亢進支援装置および脳機能亢進支援方法
US9480402B2 (en) 2011-11-11 2016-11-01 Washington University System and method for task-less mapping of brain activity
US9072496B2 (en) 2012-02-02 2015-07-07 International Business Machines Corporation Method and system for modeling and processing fMRI image data using a bag-of-words approach
US9510756B2 (en) 2012-03-05 2016-12-06 Siemens Healthcare Gmbh Method and system for diagnosis of attention deficit hyperactivity disorder from magnetic resonance images
US10137307B2 (en) * 2012-05-16 2018-11-27 Beth Israel Deaconess Medical Center, Inc. Identifying individual target sites for transcranial magnetic stimulation applications
WO2014059234A1 (en) 2012-10-11 2014-04-17 The Research Foundation Of The City University Of New York Predicting response to stimulus
JP5200311B1 (ja) 2012-10-19 2013-06-05 株式会社神島組 破砕方法
US20150294074A1 (en) * 2013-05-01 2015-10-15 Advanced Telecommunications Research Institute International Brain activity training apparatus and brain activity training system
US10357181B2 (en) 2013-05-01 2019-07-23 Advanced Telecommunications Research Institute International Brain activity analyzing apparatus, brain activity analyzing method and biomarker apparatus
CA2913945A1 (en) * 2013-05-28 2014-12-04 Laszlo Osvath Systems and methods for diagnosis of depression and other medical conditions
JP5538597B2 (ja) 2013-06-19 2014-07-02 株式会社日立製作所 異常検知方法及び異常検知システム
US9265441B2 (en) 2013-07-12 2016-02-23 Siemens Aktiengesellschaft Assessment of traumatic brain injury
US9632162B2 (en) 2013-12-06 2017-04-25 Toshiba Medical Systems Corporation Method of, and apparatus for, correcting distortion in medical images
US9694155B2 (en) * 2013-12-17 2017-07-04 Juliana Stoianova Panova Adjuvant method for the interface of psychosomatic approaches and technology for improving medical outcomes
US20170071522A1 (en) 2014-03-07 2017-03-16 The Research Foundation For The State University Of New York Method of diagnosing depression by pet imaging
WO2016069058A1 (en) * 2014-04-25 2016-05-06 The General Hospital Corporation Method for cross-diagnostic identification and treatment of neurologic features underpinning mental and emotional disorders
DE102015109853A1 (de) * 2015-06-19 2016-12-22 Nclogics Ag Assistenz- und Entscheidungssystem und Verfahren zur Auswertung von Elektroenzephalogrammen
EP3381361B1 (en) 2015-11-24 2023-12-27 Advanced Telecommunications Research Institute International Brain activity analysis device, brain activity analysis method, program, and biomarker device
DK3417074T3 (da) 2016-03-23 2019-09-16 Senzagen Ab Analytiske fremgangsmåder og arrays til bestemmelse af den hudsensibiliserende virkning af midler in vitro
US11666219B2 (en) * 2016-07-29 2023-06-06 Beth Israel Deaconess Medical Center, Inc. Methods and apparatus for network localization of neurological symptoms from focal brain lesions
JP2018089142A (ja) 2016-12-05 2018-06-14 学校法人同志社 脳機能イメージングデータからヒトの脳活動状態を推定する方法
CN106407733A (zh) 2016-12-12 2017-02-15 兰州大学 基于虚拟现实场景脑电信号的抑郁症风险筛查系统和方法
CN106725532B (zh) * 2016-12-13 2018-04-24 兰州大学 基于语音特征与机器学习的抑郁症自动评估系统和方法
WO2018147193A1 (ja) 2017-02-08 2018-08-16 日本電信電話株式会社 モデル学習装置、推定装置、それらの方法、およびプログラム
US11704790B2 (en) 2017-09-26 2023-07-18 Washington University Supervised classifier for optimizing target for neuromodulation, implant localization, and ablation
WO2019069955A1 (ja) 2017-10-03 2019-04-11 株式会社国際電気通信基礎技術研究所 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
JP6722845B2 (ja) 2017-10-03 2020-07-15 株式会社国際電気通信基礎技術研究所 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
JP2019198376A (ja) 2018-05-14 2019-11-21 キヤノンメディカルシステムズ株式会社 医用画像処理装置、医用画像処理方法、及び医用画像処理システム
JP7148311B2 (ja) 2018-08-07 2022-10-05 ファナック株式会社 製品検査装置

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US12020427B2 (en) 2017-10-03 2024-06-25 Advanced Telecommunications Research Institute International Differentiation device, differentiation method for depression symptoms, determination method for level of depression symptoms, stratification method for depression patients, determination method for effects of treatment of depression symptoms, and brain activity training device
WO2019172245A1 (ja) * 2018-03-09 2019-09-12 株式会社国際電気通信基礎技術研究所 脳活動訓練装置、脳活動訓練方法および脳活動訓練プログラム
JPWO2019172245A1 (ja) * 2018-03-09 2021-04-01 株式会社国際電気通信基礎技術研究所 脳活動訓練装置、脳活動訓練方法および脳活動訓練プログラム
JPWO2021205996A1 (ja) * 2020-04-06 2021-10-14
JP7365496B2 (ja) 2020-04-06 2023-10-19 株式会社国際電気通信基礎技術研究所 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラム、脳活動マーカー分類システムおよび脳機能結合相関値のクラスタリング分類器モデル
WO2021205996A1 (ja) * 2020-04-06 2021-10-14 株式会社国際電気通信基礎技術研究所 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラム、脳活動マーカー分類システムおよび脳機能結合相関値のクラスタリング分類器モデル
WO2022014682A1 (ja) * 2020-07-17 2022-01-20 株式会社国際電気通信基礎技術研究所 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム
JP2023143663A (ja) * 2022-03-24 2023-10-06 之江実験室 磁気共鳴及びシャムグラフニューラルネットワークに基づく脳アトラスの個別化方法及びシステム
JP7524279B2 (ja) 2022-03-24 2024-07-29 之江実験室 磁気共鳴及びシャムグラフニューラルネットワークに基づく脳アトラスの個別化方法及びシステム
US12127826B2 (en) 2022-03-24 2024-10-29 Zhejiang Lab Brain atlas individualization method and system based on magnetic resonance and twin graph neural network
CN116469541A (zh) * 2023-03-24 2023-07-21 深圳市第二人民医院(深圳市转化医学研究院) 一种抑郁症标志物及其在抑郁症预后中的应用、评估装置
CN116469541B (zh) * 2023-03-24 2024-01-30 深圳市第二人民医院(深圳市转化医学研究院) 一种抑郁症标志物及其在抑郁症预后中的应用、评估装置
KR102710036B1 (ko) * 2024-01-31 2024-09-25 이화여자대학교 산학협력단 산후우울증 발생 위험도를 예측하는 전자 장치 및 이의 동작 방법

Also Published As

Publication number Publication date
EP3692909A1 (en) 2020-08-12
US20210034912A1 (en) 2021-02-04
EP3692909A4 (en) 2021-10-27
US12020427B2 (en) 2024-06-25
JP2019155111A (ja) 2019-09-19
CN111447874A (zh) 2020-07-24
JP7140348B2 (ja) 2022-09-21
JP6722845B2 (ja) 2020-07-15

Similar Documents

Publication Publication Date Title
JP7140348B2 (ja) 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
JP6497687B2 (ja) 脳活動解析システム及び方法
JP6903518B2 (ja) 脳活動解析システム、脳活動解析装置、及び脳活動解析プログラム
WO2020075737A1 (ja) 脳機能結合相関値の調整方法、脳機能結合相関値の調整システム、脳活動分類器のハーモナイズ方法、脳活動分類器のハーモナイズシステム、及び脳活動バイオマーカシステム
WO2019069955A1 (ja) 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
JP7247159B2 (ja) 脳活動分類器のハーモナイズシステム、及び脳活動分類器プログラム
JP5641531B1 (ja) 脳活動訓練装置および脳活動訓練システム
JP6875054B2 (ja) 脳活動訓練装置、脳活動訓練方法および脳活動訓練プログラム
JP2024003782A (ja) 統合失調スペクトラム障害に関する情報を提供するための診断支援システム、診断支援装置、診断支援方法、及び支援プログラム
JP7057995B2 (ja) 哺乳動物におけるうつ病の病態プロファイルを検出するバイオマーカーとその利用
WO2022014682A1 (ja) 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム
US20240354947A1 (en) Differentiation device, differentiation method for depression symptoms, determination method for level of depression symptoms, stratification method for depression patients, determination method for effects of treatment of depression symptoms, and brain activity training device
Abdulkadir Brain MRI analysis and machine learning for diagnosis of neurodegeneration

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190327

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190327

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20190327

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190412

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20191122

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20191203

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20200131

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200318

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: 20200422

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200521

R150 Certificate of patent or registration of utility model

Ref document number: 6722845

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250