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

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

Info

Publication number
JP2019155111A
JP2019155111A JP2019061843A JP2019061843A JP2019155111A JP 2019155111 A JP2019155111 A JP 2019155111A JP 2019061843 A JP2019061843 A JP 2019061843A JP 2019061843 A JP2019061843 A JP 2019061843A JP 2019155111 A JP2019155111 A JP 2019155111A
Authority
JP
Japan
Prior art keywords
functional
correlation
classifier
subject
time point
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
JP2019061843A
Other languages
English (en)
Other versions
JP7140348B2 (ja
JP2019155111A5 (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
Publication of JP2019155111A publication Critical patent/JP2019155111A/ja
Publication of JP2019155111A5 publication Critical patent/JP2019155111A5/ja
Application granted granted Critical
Publication of JP7140348B2 publication Critical patent/JP7140348B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/243—Classification techniques relating to the number of classes
    • G06F18/2431—Multiple classes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00—Machine learning
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/70—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mental therapies, e.g. psychological therapy or autogenous training
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H70/00—ICT specially adapted for the handling or processing of medical references
    • G16H70/60—ICT specially adapted for the handling or processing of medical references relating to pathologies
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20081—Training; Learning
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20084—Artificial neural networks [ANN]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30016—Brain

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)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Hospice & Palliative Care (AREA)
  • Psychology (AREA)
  • Social Psychology (AREA)
  • Mathematical Physics (AREA)
  • Psychiatry (AREA)
  • Computing Systems (AREA)
  • Developmental Disabilities (AREA)
  • Child & Adolescent Psychology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)

Abstract

【課題】うつ病患者における治療効果を判定する判別装置を提供する。【解決手段】演算装置は、複数の対象者について、機能的結合から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関とがそれぞれ計測されて得られる情報に基づき、複数の機能的結合の第１の時点と第２の時点の相関の差の各々で張られる相関状態空間において、複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成し、第１の時点の被検者の安静時の脳内における、複数の機能的結合の第１の相関を算出し、第２の時点の被検者と同一の被検者の安静時の脳内における、複数の機能的結合の第２相関を算出し、被検者の前記複数の機能的結合の前記第２の相関と第１の相関との差に基づき、分類器により、被検者に対する治療効果を判別する、うつ病患者における治療効果の判定処理を実行する。【選択図】なし

Description

本発明は、判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置に関する。
（バイオマーカー）
生体内の生物学的変化を定量的に把握するため、生体情報を数値化・定量化した指標のことを「バイオマーカー」と呼ぶ。
ＦＤＡ(米国食品医薬品局)はバイオマーカーの位置づけを、「正常なプロセスや病的プロセス、あるいは治療に対する薬理学的な反応の指標として客観的に測定・評価される項目」としている。また疾患の状態や変化、治癒の程度を特徴づけるバイオマーカーは、新薬の臨床試験での有効性を確認するためのサロゲートマーカー(代替マーカー)として使われる。血糖値やコレステロール値などは、生活習慣病の指標として代表的なバイオマーカーである。尿や血液中に含まれる生体由来の物質だけでなく、心電図、血圧、PET画像、
骨密度、肺機能なども含まれる。またゲノム解析やプロテオーム解析が進んできたことによって、ＤＮＡやＲＮＡ、生体蛋白等に関連したさまざまなバイオマーカーが見出されている。
バイオマーカーは、疾患にかかった後の治療効果の測定だけでなく、疾患を未然に防ぐための日常的な指標として疾患の予防に、さらに副作用を回避した有効な治療法を選択する個別化医療への応用が期待されている。
ただし、神経・精神疾患の場合、生化学的もしくは分子遺伝学的観点から客観的な指標として利用可能な分子マーカーなども研究されているものの、検討段階というべき状況である。
一方で、ＮＩＲＳ（Near-infraRed Spectroscopy）技術を用いて、生体光計測により計測されたヘモグロビン信号から特徴量に応じて、統合失調症、うつ病などの精神疾患について分類を行う疾患判定システムなども報告がある（特許文献１）。
バイオマーカーは、通常の診療だけでなく、医薬品開発でも利用されている。
一般的な創薬の場合、探索研究の段階では薬力学マーカーが用いられ、非臨床段階になると候補化合物の毒性をみる毒性マーカーも利用されるようになる。
臨床試験の段階になると、評価項目として代替マーカーが利用されたり、診断マーカー、予後マーカー、予測マーカー、モニタリングマーカー、安全性マーカーなど、多様なマーカーが利用される。
「代替マーカー」とは、真のエンドポイントとの科学的な関係が証明されているマーカーのことをいう。
「診断マーカー」とは、疾患の有無や疾患の種類を特定する診断マーカーと，疾患の重症度や進行度を判定する診断マーカーの２種類がある。「予測マーカー」とは、特定の薬剤が効きやすい集団を見分けるバイオマーカーである。「予後マーカー」とは、特定の治
療に関係なく疾患の進行や回復を予想するバイオマーカーである。「モニタリングマーカー」とは、「薬効マーカー」とも呼ばれ、薬剤の効き具合，治療効果を診るマーカーである。このマーカーを長期間モニタリングすることで薬剤耐性の原因解明のヒントにもなりうる、とされる。
たとえば、標的分子が明確である場合などは、たとえば、標的分子などを、特定の薬剤が効きやすい患者の集団を見分けて選択するために、臨床試験で患者層別マーカーとして利用できる場合がある。このような場合、患者層別マーカーを利用する医薬品が承認されると、その患者層別マーカーを測定するための診断薬は「コンパニオン診断薬」と呼ばれる場合がある。あるいは、このような患者層別マーカーにより、患者に対する治療方法を選択するために行われる診断は、「コンパニオン診断」と呼ばれる場合もある。
臨床試験において、後期試験にかかるコストが非常に高いことから、後期試験に適したバイオマーカーの開発が臨床試験の効率を向上させるために必要と考えられる。しかし、このようなバイオマーカーの開発は、一般には容易ではない。
また、たとえば、生物学的指標 （バイオマーカー）を考える際、疾患病態生理におい
て、疾患原因に関連する行動学あるいは生物学的プロセスを反映しているものはtrait markerであり、患者における臨床状態像を反映しているのがstate markerと呼ばれる。trait markerは、遺伝規定性指標とも呼ばれ、state markerは、状態依存性指標とも呼ばれる。
（リアルタイムニューロフィードバック）
たとえば、従来、神経症の一型である強迫性障害 (Obsessive-Compulsive Disorder: OCD)に対しては、その治療法として、薬物療法や行動療法が知られている。薬物療法とし
ては、たとえば、選択的セロトニン再取り込阻害薬が用いられ、行動療法としては、曝露法と反応妨害法を組み合わせた曝露反応妨害法などが知られている。
一方で、神経・精神疾患の場合の治療法の可能性として、リアルタイムニューロフィードバック法が検討されている（特許文献６）。
核磁気共鳴画像法（ＭＲＩ:Magnetic Resonance Imaging）を利用して、ヒトの脳の活
動に関連した血流動態反応を視覚化する方法である機能的磁気共鳴画像法（ｆＭＲＩ：functional Magnetic Resonance Imaging）を始めとする脳機能画像法は、感覚刺激や認知
課題遂行による脳活動と安静時や対照課題遂行による脳活動の違いを検出して、関心のある脳機能の構成要素に対応する脳賦活領域を特定すること、すなわち脳の機能局在を明らかにすることにもちいられてきた。
近年、機能的磁気共鳴画像法（ｆＭＲＩ）などの脳機能画像法を用いたリアルタイムニューロフィードバック技術についての報告がある（非特許文献１）。リアルタイムニューロフィードバック技術は、神経疾患、精神疾患の治療方法として注目され始めている。
ニューロフィードバックとは、バイオフィードバックの一種であり、被検者が自身の脳活動に関するフィードバックを受けることで、脳活動を操作する方法を学ぶというものである。
たとえば、前部帯状回の活動をｆＭＲＩで計測し、それをリアルタイムで炎の大きさとして、患者にフィードバックして、炎を小さくするように努力させると、中枢性の慢性疼痛が実時間でも、また、長期的にも改善したとの報告がある（非特許文献２を参照）。
（安静時ｆＭＲＩ）
また、最近の研究では、安静時の脳も活発に活動していることがわかってきた。つまり、アクティブな活動中には鎮静化しており、休息時には活発に興奮する神経細胞群が脳内に存在する。解剖学的には左右の大脳が合わさった内側面が主な部位で、前頭葉内側面、後部帯状回、楔前部、頭頂連合野の後半部、中側頭回などである。この安静時の脳活動のベースラインを示す領域はデフォルトモードネットワーク（Default Mode Network；DMN
）と名付けられ、ひとつのネットワークとして同期しながら活動する（非特許文献３を参照）。
たとえば、健常者と精神疾患の患者で脳活動が異なるという例として、デフォルトモードネットワークにおける脳活動が挙げられる。デフォルトモードネットワークとは、安静状態において、目標指向的な課題を実行中よりも活発な脳活動がみられる部位のことを指す。統合失調症、アルツハイマー病などの疾患の患者は、健常者と比較したときに、このデフォルトモードネットワークに異常がみられることが報告されている。たとえば、統合失調症の患者は、安静状態において、デフォルトモードネットワークに属する後帯状皮質と、頭頂外側部、前頭前野内側部、小脳皮質の活動の相関が下がっているなどの報告がある。
ただし、このようなデフォルトモードネットワークが、認知機能とどのようにかかわっているのかや、このような脳領域間の機能的結合の相関と上述したニューロフィードバックとの関連などは、必ずしも明らかとはなっていない状況である。
一方で、複数の脳領域間において課題などの違いによる活動の相関関係の変化をみることにより、それら脳領域間の機能的結合（functional connectivity）が評価されている
。とくに、安静時のｆＭＲＩによる機能的結合の評価は安静時機能的結合ＭＲＩ（rs-fcMRI : resting-state functional connectivity MRI）とも呼ばれ、様々な神経・精神疾患を対象とした臨床研究も行われるようになりつつある。ただし、従来のｒｓ−ｆｃＭＲＩ法は、上述したようなデフォルトモードネットワークといったような大域的神経ネットワークの活動をみるものであって、より詳細な機能的結合について、十分な検討がされているとはいえない状況である。
（核磁気共鳴画像法）
このような核磁気共鳴画像法について、簡単に説明すると、以下のとおりである。
すなわち、従来から、生体の脳や全身の断面を画像する方法として、生体中の原子、特に、水素原子の原子核に対する核磁気共鳴現象を利用した核磁気共鳴映像法が、人間の臨床画像診断等に使用されている。
核磁気共鳴映像法は、それを人体に適用する場合、同様の人体内断層画像法である「Ｘ線ＣＴ」に比較して、たとえば、以下のような特徴がある。
（１）水素原子の分布と、その信号緩和時間（原子の結合の強さを反映）に対応した濃度の画像が得られる。このため、組織の性質の差異に応じた濃淡を呈し、組織の違いを観察しやすい。
（２）磁場は、骨による吸収がない。このため、骨に囲まれた部位（頭蓋内、脊髄など）を観察しやすい。
（３）Ｘ線のように人体に害になるということがないので、広範囲に活用できる。
このような核磁気共鳴影像法は、人体の各細胞に最も多く含まれ、かつ最も大きな磁性を有している水素原子核（陽子）の磁気性を利用する。水素原子核の磁性を担うスピン角運動量の磁場内での運動は、古典的には、コマの歳差運動にたとえられる。
以下、本発明の背景の説明のために、この直感的な古典的モデルで、簡単に核磁気共鳴の原理をまとめておく。
上述したような水素原子核のスピン角運動量の方向（コマの自転軸の方向）は、磁場のない環境では、ランダムな方向を向いているものの、静磁場を印加すると、磁力線の方向を向く。
この状態で、さらに振動磁界を重畳すると、この振動磁界の周波数が、静磁界の強さで決まる共鳴周波数ｆ０＝γＢ０／２π（γ：物質に固有の係数）であると、共鳴により原子核側にエネルギーが移動し、磁化ベクトルの方向が変わる（歳差運動が大きくなる）。この状態で、振動磁界を切ると、歳差運動は、傾き角度を戻しながら、静磁界における方向に復帰していく。この過程を外部からアンテナコイルにより検知することで、ＮＭＲ信号を得ることができる。
このような共鳴周波数ｆ０は、静磁界の強度がＢ０（Ｔ）であるとき、水素原子では、４２．６×Ｂ０（ＭＨｚ）となる。
さらに、核磁気共鳴映像法では、血流量の変化に応じて、検出される信号に変化が現れることを用いて、外部刺激等に対する脳の活動部位を視覚化することも可能である。このような核磁気共鳴映像法を、特に、機能的磁気共鳴画像法（ｆＭＲＩ：functional MRI)
と呼ぶ。
ｆＭＲＩでは、装置としては通常のＭＲＩ装置に、さらに、ｆＭＲＩ計測に必要なハードおよびソフトを装備したものが使用される。
ここで、血流量の変化がＮＭＲ信号強度に変化をもたらすのは、血液中の酸素化および脱酸素化ヘモグロビンは磁気的な性質が異なることを利用している。酸素化ヘモグロビンは反磁性体の性質があり、周りに存在する水の水素原子の緩和時間に影響を与えないのに対し、脱酸素化ヘモグロビンは常磁性体であり、周囲の磁場を変化させる。したがって、脳が刺激を受け、局部血流が増大し、酸素化ヘモグロビンが増加すると、その変化分をＭＲＩ信号として検出する事ができる。被検者への刺激は、たとえば、視覚による刺激や聴覚による刺激、あるいは所定の課題（タスク）の実行等が用いられる（たとえば、特許文献２）。
ここで、脳機能研究においては、微小静脈や毛細血管における赤血球中の脱酸素化ヘモグロビンの濃度が減少する現象（ＢＯＬＤ効果）に対応した水素原子の核磁気共鳴信号（ＭＲＩ信号）の上昇を測定することによって脳の活動の測定が行われている。
特に、人の運動機能に関する研究では、被検者に何らかの運動を行わせつつ、上記ＭＲＩ装置によって脳の活動を測定することが行われている。
ところで、ヒトの場合、非侵襲的な脳活動の計測が必要であり、この場合、ｆＭＲＩデータから、より詳細な情報を抽出できるデコーディング技術が発達してきている（たとえば、非特許文献４）。特に、ｆＭＲＩが脳におけるボクセル単位（volumetric pixel : voxel）で脳活動を解析することで、脳活動の空間的パターンから、刺激入力や認識状態を推定することが可能となっている。
さらに、このようなデコーディング技術を発展させた技術として、特許文献３には、神経・精神疾患に対して、脳機能画像法によるバイオマーカーを実現するための脳活動解析方法が開示されている。この方法では、健常群、患者群において測定された安静時機能的結合ＭＲＩのデータから、それぞれの被験者について、所定の脳領域間の活動度の相関行列を導出する。被験者の疾患／健常ラベルを含む被験者の属性と相関行列とについて正則化正準相関解析により特徴抽出が行われる。正則化正準相関解析（スパース正準相関解析）の結果に基づいて、スパースロジスティック回帰(Sparse Logistic Regression: SLR)
による判別分析によりバイオマーカーとして機能する分類器が生成される。このような機械学習の技術により、安静時のｆＭＲＩ データから導き出される脳領野間の結合（機能的結合）にもとづいて自閉症の診断結果を予測可能であることが示された（特許文献４、特許文献５）。しかも、その予測性能の検証は一つの施設において計測された脳活動のみではなく、他の施設で計測された脳活動に対しても、ある程度の汎化が可能であることが示された。
さらに、近年では、多数の研究施設から得られた大規模なマルチサイトサンプル群（ｎ＝１１８８）でｆＭＲＩを用いて、うつ病患者は４つの神経生理学的サブタイプ（バイオタイプ）に分類できること、こうしたバイオタイプは辺縁系および前頭葉-線条体ネット
ワークにおける結合性機能不全を示すそれぞれ別のパターンによって明確に示されることが指摘されている（非特許文献５を参照）。非特許文献５においては、これらのバイオタイプは、臨床的特徴を基盤とするだけでは弁別できないが、異なる臨床症状プロファイルと結び付いており、また経頭蓋磁気刺激治療に対する応答性を予測する（ｎ＝１５４）、とされている。
非特許文献５においては、脳内の２５８の領域からのＢＯＬＤ信号に基づき、脳内の機能的な結合により、上記のサブタイプの分類を行っている。
再表２００６−１３２３１３号公報 特開２０１１−０００１８４号公報 特開２０１５−６２８１７号公報 特許第６１９５３２９号公報 国際公開第２０１７／０９０５９０号 特許第５６４１５３１号公報
Nikolaus Weiskopf ，"Real-time fＭＲＩ 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.
以上のように、機能的磁気共鳴画像法などの脳機能画像法による脳活動の解析の神経・精神疾患の治療への応用を考えた場合、たとえば、上述したようなバイオマーカーとして、脳機能画像法による脳活動の解析は、非侵襲的な機能マーカーとして、診断法の開発、根本治療を実現するための創薬に向けた標的分子の探索・同定などへの応用も期待される。
たとえば、これまで、うつ病に対しては、遺伝子を用いた実用的なバイオマーカーの完成には至っておらず、それゆえに、客観的な薬物の効果判定を行うことが容易ではなく、治療薬の開発も困難であった。
脳活動の計測データに基づいて、機械学習により分類器を生成し、これをバイオマーカーとして実用化させるためには、一つの施設において計測された脳活動に対して機械学習により生成されたバイオマーカーの予測精度の向上が必要である。また、このようにして生成されたバイオマーカーが、他の施設において計測された脳活動に対しても、汎化可能であることが必要である。
すなわち、脳活動の計測データに基づく、機械学習により分類器を構築する際には、主な２つの課題が存在する。
第１の課題は、小標本サイズの問題である。
参加者の人数であるデータ量Ｎが、計測された脳活動計測データの次元Ｍよりはるかに小さいので、分類器のパラメータは、容易に訓練データに過剰適合（「過適応」または「過学習」と呼ばれる）することになる。
この過剰適合のために、構築された分類器は、新しくサンプリングされたテストデータについては、非常に貧弱なパフォーマンスを示すことになる。これは、これらのテストデータは分類器を訓練するのに使用されていないものだからである。
したがって、分類器の望ましい汎化に対して本質的な特徴量だけを識別し利用するために正則化を適切に導入する必要がある。
たとえば、上述した非特許文献５においては、必ずしも、このような正則化について十分な考慮がされているとはいえない。
第２の課題は、構築された分類器が、訓練データが集められたサイトとは異なる撮像サイトで走査されたＭＲＩデータに対しても、よいパフォーマンスを維持する場合のみ、分類器が臨床的に有用で、科学的に信頼できるということである。
これは撮像サイトにわたる、いわゆる汎化能力である。
しかしながら、臨床応用では、異なるサイトで走査されたデータに、特定のサイトで取得したデータを使用して訓練された分類器を汎化することができないことはしばしば観測されている。
この発明は、上記のような問題点を解決するためになされたものであって、その１つの
目的は、脳の活動状態に対して、うつ症状の客観的な疾患ラベルの判別を実現することである。また、この発明の他の目的は、うつ症状に対する治療効果の程度を表す情報の判別を実現することである。この発明の他の目的は、うつ症状の状態を客観的に判別するためのバイオマーカーとしての指標値を出力する分類器を提供することである。この発明の他の目的は、うつ症状のレベルの判別を実現することである。この発明の他の目的は、うつ病患者の層別化を実現することである。
本発明のある実施形態は、被検者におけるうつ症状の有無の判定を補助するための判別装置に関する。本実施形態の判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備る。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。本実施形態の判別装置は、演算装置をさらに備え、前記演算装置は、前記分類器により、被検者のうつ症状に対する分類結果を生成する判別処理を実行する。
本発明のある実施形態は、うつ症状のレベルの判別装置に関する。本実施形態の判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備える。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成する。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。本実施形態の判別装置は、演算装置をさらに備え、前記演算装置は、前記分類器により、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する。前記演算装置は、前記指標値をあらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較する。前記演算装置は、前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する。
本発明のある実施形態は、被検者における治療効果の判定するための判別装置に関する。本実施形態の判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備える。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。本実施形態の判別装置は、演算装置をさらに備え、前記演算装置は、前記分類器により、第１の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第１の値を生成する。前記演算装置は、前記分類器により、
治療開始後であって第１の時点より後の第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素についてうつ症状を評価する第２の値を生成する。前記演算装置は、前記第１の値と前記第２の値とを比較する。本実施形態の判別装置は、前記第２の値が前記第１の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定する、及び／又は前記第２の値が前記第１の値よりも改善していない場合に、前記治療が被検者のうつ症状の改善に有効でないと決定する。
本発明のある実施形態は、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、うつ病患者の層別化するための判別装置に関する。本実施形態の判定装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置を備える。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。本実施形態の判定装置は、演算装置をさらに備え、前記演算装置は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する。前記演算装置は、前記指標値を、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較する。前記演算装置は、前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する。
本発明のある実施形態は、前記判別装置に上記処理を実行させるためのコンピュータプログラムに関する。
本発明のある実施形態は、うつ症状を有する被検者を判定することを補助する判別方法に関する。本実施形態の判別方法は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を含む。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも一つを含む。本実施形態の判別方法は、前記指標値が基準値を超えた場合に、前記被検者がうつ症状を有すると決定する工程を含む。
本発明のある実施形態は、被検者においてうつ症状のレベルを判定することを補助するため判別方法に関する。本実施形態の判別方法は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を備える。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも一つを含む。本実施形態の判別方法は、前記指標値を、あらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較する工程と、前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する工程とを含む。
本発明のある実施の形態は、被検者における治療効果の判定を補助するための判別方法に関する。本実施形態の判別方法は、第１の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第１の値を生成する工程を備える。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質
との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも一つを含む。また、本実施形態の判別方法は、治療開始後であって第１の時点より後の第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素ついてうつ症状を評価する第２の値を生成する工程を含む。本実施形態の判別方法は、前記第１の値と前記第２の値とを比較する工程と、前記第２の値が前記第１の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定する工程と、を含む。
本発明のある実施の形態は、被検者における治療効果の判定を補助するための判別方法に関する。本実施形態の判別方法は、第１の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第１の値を生成する工程を備える。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも一つを含む。本実施形態の判別方法は、治療開始後であって第１の時点より後の第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素ついてうつ症状を評価する第２の値を生成する工程を含む。本実施形態の判別方法は、前記第１の値と前記第２の値とを比較する工程と、前記第２の値が前記第１の値よりも改善していない場合に、前記治療が被検者のうつ症状の改善に有効でないと決定する工程と、を含む。
本発明のある実施の形態は、うつ病患者の層別化を補助するための判別方法に関する。本実施形態の判別方法は、うつ病が予め設定された複数のサブクラスに層別化されている場合に、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程を備える。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも一つを含む。本実施形態の判別方法は、前記指標値を、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較する工程と、前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する工程と、を含む。
本発明のある実施の形態は、うつ病患者の層別化を補助するための判別装置に関する。本実施形態の判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を特定する情報を記憶するための記憶装置と、層別化処理を実行する演算装置とを備える。本実施形態の判別装置において、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ病のサブクラスのラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記選択された機能的結合に対応する判別対象の相関行列の要素に対して前記うつ病のサブクラスのラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。前記演算装置は、前記重み付け和と前記判別対象の相関行列の要素から前記サブクラスを判別するための処理を実行する。。
本発明のある実施の形態は、うつ病患者における治療効果の判定を補助するための判別装置に関する。本実施形態の判別装置は、演算装置を備え、前記演算装置は、複数の対象
者について、下記表１に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関をそれぞれ計測して、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する分類器生成処理を実行する。前記演算装置は、第１の時点の被検者の安静時の脳内における、前記複数の機能的結合の第１の相関を計測し、第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第２相関を計測する。前記演算装置は、前記被検者の前記複数の機能的結合の前記第２の相関と前記第１の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する。
本発明のある実施の形態は、うつ病患者における治療効果の判定を補助するための判別方法に関する。本実施形態の判別方法は、複数の対象者について、下記表１に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関をそれぞれ計測して、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する工程を備える。本実施形態の判別方法は、第１の時点の被検者の安静時の脳内における、前記複数の機能的結合の第１の相関を計測する工程と、第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第２相関を計測する工程を備える。本実施形態の判別方法は、前記被検者の前記複数の機能的結合の前記第２の相関と前記第１の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する工程と、を含む。
本発明のある実施の形態は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、被検者におけるうつ症状の有無の判定またはうつ症状のレベルの判定を補助するために使用する方法に関する。前記使用方法は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別する前記分類器を生成する工程を備える。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程をさらに備える。
本発明のある実施の形態は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、うつ症状の治療の効果の判定を補助するために使用する方法に関する。前記使用方法は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別する前記分類器を生成する工程を備える。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程をさらに備える。
本発明のある実施の形態は、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、安
静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により生成された分類器を、被検者を前記サブクラスに層別化することを補助するために使用する方法に関する。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。前記使用方法は、前記分類器に被検者の前記複数の機能的結合の相関行列の要素から生成された指標値を入力する工程を含む。
本発明のある実施の形態は、 被検者における治療効果の判定を補助するための判別装置に関する。前記判別装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から第１の分類器を生成する第１の演算装置を含む分類器生成装置を備る。前記第１の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。前記第１の演算装置は、複数の対象者について、下記表５に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関をそれぞれ計測して、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第２の分類器を生成する。前記判別装置は、第２の演算装置と記憶装置を含む層別化装置をさらに備る。前記記憶装置は、前記第１および第２の分類器の情報と、うつ症状が予め設定された複数のサブクラスに層別化するための情報とを記憶する。前記第２の演算装置は、前記第１の分類器を用いて、被検者を前記サブクラスに層別化処理を実行する。前記第２の演算装置は、前記層別化処理によって特定のサブクラスに層別化された被検者について、第１の時点の安静時の前記複数の機能的結合の第１の相関を計測するための処理を実行する。前記第２の演算装置は、第２の時点の前記被検者と同一の被検者の安静時の前記複数の機能的結合の第２相関を計測するための処理を実行する。前記第２の演算装置は、前記被検者の前記複数の機能的結合の前記第２の相関と前記第１の相関との差に基づき、前記第２の分類器により、前記被検者に対する治療効果を判別する。
本発明のある実施の形態は、被検者における治療効果の判定を補助するための判別方法に関する。本実施形態の判別方法は、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から第１の分類器生成処理によって生成された第１の分類器を使用して、被検者を前記サブクラスに層別化する層別化工程を含む。本実施形態の判別方法は、前記第１の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機
能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。本実施形態の判別方法は、前記層別化工程によって特定のサブクラスに層別化された被検者について、投薬前の安静時の前記複数の機能的結合の第１の相関を計測する第１の相関計測工程と、投薬開始から所定期間経過後の前記被検者と同一の被検者の安静時の前記複数の機能的結合の第２相関を計測する第２の相関計測工程とを含む。本実施形態の判別方法は、特定のサブクラスに分類された前記被検者に対する薬効を第２の分類器を用いて判別する判別工程であって、前記第２の分類器は、複数の対象者について、下記表１に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と投薬後であって第１の時点よりも後の第２の時点の相関をそれぞれ計測する。そして、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第２の分類器生成処理によって予め生成され前記被検者の前記複数の機能的結合の前記第２の相関と前記第１の相関との差に基づき、前記被検者に対する薬効を判別する、判別工程と、を備える。
本発明のある実施の形態は、演算装置と、記憶装置とを備えた第１の分類器生成装置に関する。前記演算装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器を特定する情報を生成する。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記選択された複数の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも１つを含む。前記記憶装置は、前記演算装置によって生成された前記第１の分類器を特定する情報を記憶する。
本発明のある実施の形態は、第１の分類器の生成方法に関する。本実施形態では、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理により第１の分類器を生成する。前記第１の分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された複数の機能的結合の重み付け和に基づいて、前記うつ症状の疾患ラベルを判別するように生成される。前記機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合、および左下前頭回弁蓋部と、右背内側前頭前皮質および右補足運動野との間の第２の機能的結合から選択される少なくとも一つを含む。
本発明のある実施の形態は、演算装置と、記憶装置とを備えた第２の分類器生成装置に関する。前記演算装置は、複数の対象者について、下記表３に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関をそれぞれ計測する。本実施形態は、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する第２の分類器を生成する。前記記憶装置は、前記演算装置によって生成された前記第２の分類器を特定するための情報を記憶する。
本発明のある実施の形態は、第２の分類器生成方法に関する。本実施形態では、複数の対象者について、下記表４に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関をそれぞれ計測する。本実施形態は、前記複数の機能的結合の第１の時点と第２の時点の
前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する工程を備える。
本発明のある実施の形態は、ニューロフィードバック訓練を実行するための脳活動訓練装置に関する。前記脳活動訓練装置は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を時系列で予め測定した信号に基づいて、前記複数の所定領域間の機能的結合から機械学習により、うつ症状の疾患ラベルを判別するために特徴選択により選択された複数の機能的結合のうち訓練対象の機能的結合を特定する情報を格納する記憶装置を備える。前記訓練対象の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合を含む。前記記憶装置は、ニューロフィードバック訓練における前記訓練対象の機能的結合の目標パターンを記憶する。前記脳活動訓練装置は、前記ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知するための脳活動検知装置をさらに備える。前記被訓練者の脳内の複数の所定領域のそれぞれは、前記複数の参加者における脳内の複数の所定領域にそれぞれ対応する。前記脳活動訓練装置は、提示装置と、演算装置とをさらに備える。前記演算装置は、ｉ）前記脳活動検知装置により検知された信号から、前記訓練対象の機能的結合の時間相関を所定期間について算出し、ｉｉ）算出された前記時間相関に基づいて、前記目標パターンとの近似度に応じて、報酬値を算出し、ｉｉｉ）前記報酬値の大きさを示す情報を前記提示装置により前記被訓練者に対して提示する。
本発明のある実施の形態は、ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知するための脳活動検知装置と、演算装置と、記憶装置と、提示装置と、を備える脳活動訓練装置の制御方法に関する。前記制御方法は、前記記憶装置に、ｉ）安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を時系列で予め測定した信号に基づいて、前記複数の所定領域間の機能的結合から機械学習により、うつ症状の疾患ラベルを判別するために特徴選択により選択された複数の機能的結合のうち訓練対象の機能的結合を特定する情報と、ｉｉ）ニューロフィードバック訓練における前記訓練対象の機能的結合の目標パターンとを格納する工程を備える。前記訓練対象の機能的結合は、左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の第１の機能的結合を含む。前記制御方法は、前記脳活動検知装置により、前記ニューロフィードバック訓練の被訓練者の脳内の複数の所定領域における脳活動を示す信号を時系列で検知する工程をさらに備える。前記被訓練者の脳内の複数の所定領域のそれぞれは、前記複数の参加者における脳内の複数の所定領域にそれぞれ対応する。前記制御方法は、前記演算装置が、前記脳活動検知装置により検知された信号から、前記訓練対象の機能的結合の時間相関を所定期間について算出する工程と、前記演算装置が、算出された前記時間相関に基づいて、前記目標パターンとの近似度に応じて、報酬値を算出する工程と、前記演算装置が、前記報酬値の大きさを示す情報を前記提示装置により前記被訓練者に対して提示する工程と、をさらに備える。
この発明によれば、脳の活動状態に対してうつ症状の状態に対する指標を客観的に判別する判別装置及び判別方法を提供することができる。この発明によれば、うつ症状に対する治療効果の程度を表す情報を客観的に判別することができる。この発明によれば、うつ症状を客観的に判別するためのバイオマーカーとしての指標値を出力する分類器を提供することができる。この発明によれば、うつ症状のレベルを判別することができる。この発明によれば、うつ病患者を客観的に層別化することができる。
この発明によれば、ニューロフィードバック訓練を行うことができる。
ＭＲＩ装置１０の全体構成を示す模式図である。 データ処理部３２のハードウェアブロック図である。 関心領域について、安静時の機能的結合の相関を示す相関行列を抽出する手続きを示す概念図である。 相関行列から、第１の分類器を生成する過程を説明する概念図である。 図４で説明したような第１の分類器の生成処理および生成された第１の分類器による判別処理を行うための機能ブロック図である。 第１の分類器を生成するためにデータ処理部３２が行う処理を説明するためのフローチャートである。 内ループ特徴選択を説明するための概念図である。 内ループ特徴選択処理の概念を説明するための図である。 一例として、内ループ特徴選択を行った際の反復処理のうち、特定のハイパーパラメータλ1およびλ2に対する結果を示す図である。 内ループ特徴選択の処理を、より詳しく説明するためのフローチャートである。 外ループ特徴選択を説明するための概念図である。 外ループ特徴選択の処理を、より詳しく説明するためのフローチャートである。 ステップＳ１０８での判別器の生成処理の概念を示す図である。 被検者についてうつ症状の判別を行うためのフローチャートである。 データ収集、推定処理および対象者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。 被検者についてうつ症状のレベルの判別を行うためのフローチャートである。 被検者について治療効果を判別するためのフローチャートである。 被検者について治療効果を判別するためのフローチャートである。 うつ病のサブクラスを判別するためのフローチャートである。 分類器の生成処理および生成された分類器による層別化処理を行うための機能ブロック図である。 図１９に示す、データ収集、推定処理および被検者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。 第２の分類器を生成するためのフローチャートである。 第２の分類器を使用して薬効を判定するためのフローチャートである。 第２の分類器の生成処理および生成された第２の分類器による治療効果判別処理を行うための機能ブロック図である。 うつ病のサブクラスを決定した被検者に対して第２の分類器を用いて薬効を判定するためのフローチャートである。 選択された１２対の機能的結合を示す。 ５４個の非ゼロのｃi kの大きさの分布を示す図である。 ａは、広島のコホートの関連重み付き和の分布を示す。ｃは、千葉のコホートの関連重み付き和の分布を示す。ｃは、訓練データとテストデータを入れ替えて判別成績を検討した結果を示す。†p < 0.10*p < 0.05, **p < 0.01。 ｄは、メランコリー型ＭＤＤ群の関連重み付き和の平滑化ヒストグラムを示す。ｅは、非メランコリー型ＭＤＤ群の関連重み付き和の平滑化ヒストグラムを示す。†p < 0.10*p < 0.05, **p < 0.01。 ｆは、治療抵抗性ＭＤＤ群の関連重み付き和の平滑化ヒストグラムを示す。ｇは、ＡＳＤ群の関連重み付き和の平滑化ヒストグラムを示す。ｈは、ＳＳＤ群の関連重み付き和の平滑化ヒストグラムを示す。†p < 0.10*p < 0.05, **p < 0.01。 並び替え検定の結果を示す。aは、訓練データのLOOCVの並び替え検定（1000回繰り返し）のヒストグラムを示す。bは、サンプル外のテストデータセットの精度を示し、二項分布を曲線として示す。並び替え無しで訓練され、検証されたメランコリックMDD分類器の精度は、垂直線で示す。並び替え検定の結果は、LOOCV（p=0.049）およびサンプル外検証（p=0.036）について有意であった。 * p <0.05。 ａは、ＢＤＩスコアと第１の分類器を用いた判別結果の相関（全ＭＤＤと健常者全体）を示す。ｂは、ＢＤＩと第１の分類器を用いた判別結果の相関（全ＭＤＤ全体）を示す。ｃは、治療前後の関連重み付き和の平滑化ヒストグラムを示す。ｄは、投薬後のＢＤＩスコアの変化と第１の分類器を用いた判別結果の変化の相関を示す。 ａは、各機能的結合における健常者対照群とＭＤＤ群の相関の差と、投薬前と投薬後の相関の差を示す。ｂは、広島のコホートについての、ＦＣ１とＦＣ２における投薬前と投薬後の相関の差を示す。ｃは、千葉のコホートについての、ＦＣ１とＦＣ２における投薬前と投薬後の相関の差を示す。ｄは、健常者のコホートについての、ＦＣ１とＦＣ２における投薬前と投薬後の相関の差を示す。 ｅは、投薬後寛解群と投薬後非寛解群のＦＣ１とＦＣ２における投薬前と投薬後の相関の差を示す。ｆは、符号付ΔＦＣ１とΔＦＣ２の投薬群での分布を示す。 脳活動訓練装置の構成の概念を示す図である。 指標値と目標値の近さを表すためのモニター表示の一例である。 ニューロフィードバックにおける訓練シーケンスの一例を示す図である。 ａは、ＭＤＤ参加者の訓練期間中のニューロフィードバックスコアを示す。ｂは、ＭＤＤ参加者の訓練前後におけるＨＤＲＳスコアを示す。ｃは、亜臨床的うつ病の参加者の訓練期間中のニューロフィードバックスコアを示す。ｄは、訓練後のＢＤＩスコアとｒｓ−ｆｃＭＲＩ変化量の相関を示す。
以下、本発明の実施の形態について、図に従って説明する。なお、以下の実施の形態において、同じ符号を付した構成要素および処理工程は、同一または相当するものであり、必要でない場合は、その説明は繰り返さない。
１．用語の説明
本発明において、うつ症状には、抑うつ気分；興味や関心の低下；気力の低下；焦燥；制止；思考力、集中力、又は決断力の減退；無価値感、又は罪責感；自殺、自殺念慮、又は自殺企画；病的思考内容；妄想；身体症状（全身倦怠感、頭痛、頭重、又は腰痛等の身体諸所の痛み、動悸、息切れ、食欲減退、又は体重減少等）、睡眠障害から選択される少なくとも一種が含まれる。好ましくは、前記うつ症状には、精神障害の診断と統計マニュアル（ＤＳＭ）−ＩＶの基準の大うつ病性障害（ＭＤＤ）に伴う症状が含まれる。
本発明において、うつ病は、前記うつ症状を伴う限り制限されないが、好ましくはＭＤＤである。ＭＤＤには、メランコリー型ＭＤＤ、非メランコリー型ＭＤＤ及び治療抵抗性ＭＤＤが含まれる。本明細書においては、ＭＤＤを単にうつ病と呼ぶことがある。
抑うつの程度の評価には、たとえば、従来から、自己申告によるＢｅｃｋ Ｄｅｐｒｅｓｓｉｏｎ Ｉｎｖｅｎｔｏｒｙ （ベックうつ病調査表）による評価が、診察時のスクリーニングや面接時の補助資料として利用されている。また、医師用の評価尺度として使用されるハミルトンうつ病評価尺度（ＨＤＲＳ：Hamilton Depression Rating Scale）は、うつ病の指標を提供する多重項目質問票であり、回復を評価する指標となる。略称は、HAM-Dである。
本発明において、被検者は、制限されないが、前記うつ症状を有する者であることが好ましい。年齢及び性別は制限されない。前記被検者は、前記うつ症状を改善するための治
療を受けていない者であっても、前記治療を受けている者であってもよい。
うつ症状を改善するための治療には、投薬療法、ニューロフィードバック療法、修正型電気けいれん療法および反復経頭蓋磁気刺激療法から選択される少なくとも一種が含まれる。投薬される薬剤としては、三環系抗うつ薬（イミプラミン、トリミプラミン、クロミプラミン、アミトリプチリン、ノルトリプチリン、アモキサピン、ロフェプラミン、及びドスレピン等）、四環系抗うつ薬（マプロチリン、アミンセリン、及びセチプチリン等）、トラゾドン、選択的セロトニン再取り込み阻害薬（エスシタロプラム、フルボキサミン、パロキセチン、及びセルトリラリン等）、セロトニン−ノルエピネフリン再取り込み阻害剤（ミルナシプラン、及びデュロキセチン等）、ノルアドレナリン・セロトニン作動性抗うつ剤（ミルタザピン等）、及びベンザミド系薬剤（スルピリド等）から選択される少なくとも一種を挙げることができる。
うつ症状の治療には、「治癒」と言う概念がないため、うつ症状が改善するとは、BDI
等の臨床所見から見て、以前の状態よりもうつ症状がよくなっている状態、または臨床所見から見て「寛解（remitted）」といえる状態になることを意味する。
本発明の実施の形態には、被検者においてうつ症状のレベルを判定することを補助すること、被検者におけるうつ症状の有無の判定を補助すること、被検者における治療効果の判定を補助すること、及び被検者を疾患のサブクラスに層別化することを補助することが含まれる。
本発明において、ドラッグリプロファイリングとは、既にヒトに対して別の疾患または症状に対する薬効が明らかとなっている薬剤、またはヒトに対して毒性も弱いが有意な薬効も認められないことが明らかになっている薬剤であって、うつ症状に対する影響が未知である薬剤について、うつ症状に対する薬効を検出することを意図する。
２．安静時ｆＭＲＩの撮像条件と相関行列の抽出
安静時ｆＭＲＩの撮像装置は、制限されない。撮像条件もｆＭＲＩ画像を取得できる限り制限されない。例えば、ｍａｇｎｅｔｉｃ ｆｉｌｅｄは３．０Ｔ程度、ｆｉｅｌｄ ｏｆ ｖｉｅｗは１９２ｍｍ〜２５６ｍｍ程度、ｍａｔｒｉｘは６４×６４程度、ｎｕｍｂｅｒ ｏｆ ｓｌｉｃｅｓは３０枚〜４０枚程度、ｎｕｍｂｅｒ ｏｆ ｖｏｌｕｍｅｓは１１２〜２４４程度、ｓｌｉｃｅ ｔｈｉｃｋｎｅｓｓは３．０ｍｍ〜４．０ｍｍ程度、ｓｌｉｃｅ ｇａｐは０ｍｍ〜０．８ｍｍ程度、ＴＲは２，０００ｍｓ〜２，７００ｍｓ程度、ＴＥは２５ｍｓ〜３１ｍｓ程度、総撮像時間は５分〜１０分程度、Ｆｌｉｐ ａｎｇｌｅは、７５ｄｅｇ〜９０ｄｅｇ程度、ｓｌｉｃｅ ａｃｑｕｉｓｉｔｉｏｎ ｏｒｄｅｒは、Ａｓｃｅｎｄｉｎｇ （Ｉｎｔｅｒｌｅａｖｅｄ）である。撮像は、暗い照明下で行うことが好ましい。また被検者には、撮像中何も考えず、眠らずにいてもらうことが好ましい。さらに、被検者には、撮像中モニター画面中央の十字マークを見続けてもらうことが好ましい。
撮像されたｆＭＲＩ画像は、Ｙａｈａｔａらの文献（ＮＡＴＵＲＥ ＣＯＭＭＵＮＩＣＡＴＩＯＮＳ ｜ ７：１１２５４ ｜ ＤＯＩ： １０．１０３８／ｎｃｏｍｍｓ１１２５４）に記載の方法で処理することができる。
撮像されたｆＭＲＩ画像データについて、特に限定されないが、たとえば、Ｍａｔｌａｂ Ｒ２０１４ａ（Ｍａｔｈｗｏｒｋｓ ｉｎｃ．， ＵＳＡ）のＳＰＭ８（Ｗｅｌｌｃｏｍｅ Ｔｒｕｓｔ Ｃｅｎｔｅｒ ｆｏｒ Ｎｅｕｒｏｉｍａｇｉｎｇ、Ｕｎｉｖｅｒｓｉｔｙ Ｃｏｌｌｅｇｅ Ｌｏｎｄｏｎ、ＵＫ）を用いてＴ１強調構造画像及び安静状態機能画像における標準的な前処理を施すことができる。例えば、撮像画像に対して、Ｒ
ｅａｌｉｇｎ、Ｓｌｉｃｅ−ｔｉｍｉｎｇ ｃｏｒｒｅｃｔｉｏｎ、Ｃｏｒｅｇｉｓｔｅｒｉｎｇ、Ｎｏｒｍａｌｉｚａｔｉｏｎ及びＳｍｏｏｔｈｉｎｇ（ＦＷＨＭ＝８ｍｍ）を行うことができる。さらに、各被検者の全画像データについて，直前の画像から０．５ｍｍ 以上動いたと判断される画像を解析から除外してもよい。
前処理を行った画像データから、脳内の所定の複数の関心領域（Ｒｅｇｉｏｎ ｏｆ ｉｎｔｅｒｅｓｔ， ＲＯＩ）間について、機能的結合（Ｆｕｎｃｔｉｏｎａｌ ｃｏｎｎｅｃｔｉｖｉｔｙ， ＦＣ）（すなわち、結合強度）を算出する。機能的結合は、安静時脳活動解析において一般的に利用される特徴量（相関行列の要素）であり、異なる関心領域の時系列信号間のピアソンの相関係数により定義される。特に限定されないが、たとえば、所定の時間にわたって計測された値からのピアソンの相関係数の平均値により、結合強度を表すことができる。ただし、所定の時間内のピアソンの相関係数についての他の統計量により結合強度を表してもよい。
なお、後述するように、本実施の形態では、うつ症状の状態を判別するためには、たとえば、下記表１に示す機能的結合識別番号（表１でＩＤ）１〜１２の関心領域間について、その全て、または、一部について機能的結合（結合強度）を算出する。
表１で、Ｌａｔ．の“Ｌ”と“Ｒ”は、左脳と右脳の区別を示す。ＢＳＡは、ブロードマン領野を示し、ＢＡは、ブロードマンの領野の番号を示す。Ｗｅｉｇｈｔは、後述する関連重み付け和（単に「重み付け和」と称することもある）の重みを示す。
相関行列の要素の抽出は、特に限定されないが、たとえば、以下のような手順で実行することができる。
はじめに各関心領域に含まれる全ボクセル（Ｖｏｘｅｌ）の平均信号の時系列を抽出する。次に信号値からのノイズ除去のため、バンドパスフィルタ（０．００８〜０．１Ｈｚ）をかける。その後、９つの説明変数（全脳、白質、脳脊髄液の平均信号及び６つの体動補正パラメータ）による回帰を行う。回帰後の残差系列を機能的結合と関連する時系列信号値と考え、前記時系列信号値を各関心領域の要素とし、関心領域の各対の関心領域の各相関行列の要素について、時系列間でピアソンの相関係数を計算する。前記相関係数が機能的結合の結合強度を表す値であり、各対それぞれに対応する結合強度が求められる。
後述する分類器に入力データとして前記各機能的結合の相関係数を入力することにより、前記相関係数と、分類器の分類処理を実行するための各機能的結合について予め算出された重み（寄与度）を示す係数とに基づいて、機能的結合についてうつ症状の疾患ラベルを判別するための指標値が算出される。あるいは前記相関係数と前記係数とに基づいて、複数の機能的結合についての関連重み付け和がうつ症状の疾患ラベルを判別するための指標値として算出される。ここで、「関連重み付け和」とは、複数の機能的結合にそれぞれ対応する重みを乗算して和をとったものをいう。
このため、前記指標値は、単に被検者の脳活動を測定したデータではなく、各機能的結合の重みを考慮して人為的に算出された値である。前記指標値は、うつ症状のラベルの判別のために、うつ症状のレベルの判別のために、治療効果の判別のために、または、うつ病患者の層別化のために使用される。
前記１２対の機能的結合のうち、最もうつ症状への寄与度が高いのは、機能的結合識別番号１で示される機能的結合である。次に寄与度が高いのは、機能的結合識別番号２で示される機能的結合である。したがって、以下に説明する各実施の形態においては、少なくとも機能的結合識別番号１及び機能的結合識別番号２の一方、あるいは両方を選択して使用してもよい。
前記１２対の機能的結合は、特にメランコリー型ＭＤＤ患者群と健常対照群を判別することに適している。
３．第１の分類器の生成と判別装置
第１の分類器は、安静時において、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号を脳活動検知装置により時系列で予め測定した信号から分類器生成処理によって生成される。ここで、うつ病患者とは、予め医師の診断により、うつ病と診断され、うつ病の「疾患ラベル」と関連づけられた参加者を意味する。前記分類器は、前記複数の所定領域間の機能的結合のうちから機械学習により、うつ症状の疾患ラベルに関連するとして特徴選択により選択された機能的結合の重みに基づいて、前記うつ症状の疾患ラベルを判別するように生成される。
すなわち、後述するように、予め設定された複数の関心領域に対する機能的結合のうち、うつ症状の疾患ラベルに特異的に関連するとしてスパース正準相関分析により抽出された機能的結合から、スパースロジスティック回帰による特徴選択により、うつ病の疾患ラ
ベルの判別に使用される機能的結合が選択される。そして、このようにして選択された機能的結合から、上述した関連重み付け和が算出される。
３−１．判別装置１
例えば、第１の分類器は、図１に示すＭＲＩ装置１０から取得されるｆＭＲＩの画像データから生成される。
図１は、本発明の実施の形態１にかかる、判別装置１であるＭＲＩ装置１０の全体構成を示す模式図である。
図１に示すように、ＭＲＩ装置１０は、ＭＲＩを撮像するＭＲＩ撮像部２５とＭＲＩ撮像部２５の制御シーケンスを設定するとともに各種データ信号を処理して画像を生成するデータ処理部３２とを備える。ＭＲＩ撮像部２５は、被検者２の関心領域に制御された磁場を付与してＲＦ波を照射する磁場印加機構１１と、この被検者(又は被訓練者)２からの応答波（ＮＭＲ信号）を受信してアナログ信号を出力する受信コイル２０と、この被検者２に付与される磁場を制御するとともにＲＦ波の送受信を制御する駆動部２１とを備える。
なお、ここで、被検者２が載置される円筒形状のボアの中心軸をＺ軸にとりＺ軸と直交する水平方向にＸ軸及び鉛直方向にＹ軸を定義する。
ＭＲＩ装置１０は、このような構成であるので、磁場印加機構１１により印加される静磁場により、被検者２を構成する原子核の核スピンは、磁場方向（Ｚ軸）に配向するとともに、この原子核に固有のラーモア周波数でこの磁場方向を軸とする歳差運動を行う。
そして、このラーモア周波数と同じＲＦパルスを照射すると、原子は共鳴しエネルギーを吸収して励起され、核磁気共鳴現象（ＮＭＲ現象；Nuclear Magnetic Resonance）が生じる。この共鳴の後に、ＲＦパルス照射を停止すると、原子はエネルギーを放出して元の定常状態に戻る緩和過程で、ラーモア周波数と同じ周波数の電磁波（ＮＭＲ信号）を出力する。
この出力されたＮＭＲ信号を被検者２からの応答波として受信コイル２０で受信し、データ処理部３２において、被検者２の関心領域が画像化される。
磁場印加機構１１は、静磁場発生コイル１２と、傾斜磁場発生コイル１４と、ＲＦ照射部１６と、被検者２をボア中に載置する寝台１８とを備える。
被検者２は、寝台１８に、たとえば、仰臥する。被検者２は、特に限定されないが、たとえば、プリズムメガネ４により、Ｚ軸に対して垂直に設置されたディスプレイ６に表示される画面を見ることができる。このディスプレイ６の画像により、被検者２に視覚刺激が与えられる。なお、被検者２への視覚刺激は、被検者２の目前にプロジェクタにより画像が投影される構成であってもよい。
このような視覚刺激は、上述したニューロフィードバックでは、フィードバック情報の提示に相当する。
駆動部２１は、静磁場電源２２と、傾斜磁場電源２４と、信号送信部２６と、信号受信部２８と、寝台１８をＺ軸方向の任意位置に移動させる寝台駆動部３０とを備える。
データ処理部３２は、操作者（図示略）から各種操作や情報入力を受け付ける入力部４
０と、被検者２の関心領域に関する各種画像及び各種情報を画面表示する表示部３８と、各種処理を実行させるプログラム・制御パラメータ・画像データ（構造画像等）及びその他の電子データを記憶する記憶部３６と、駆動部２１を駆動させる制御シーケンスを発生させるなどの各機能部の動作を制御する制御部４２と、駆動部２１との間で各種信号の送受信を実行するインタフェース部４４と、関心領域に由来する一群のＮＭＲ信号からなるデータを収集するデータ収集部４６と、このＮＭＲ信号のデータに基づいて画像を形成する画像処理部４８と、ネットワークとの間で通信を実行するためのネットワークインタフェース部５０を備える。
また、データ処理部３２は、専用コンピュータである場合の他、各機能部を動作させる機能を実行する汎用コンピュータであって、記憶部３６にインストールされたプログラムに基づいて、指定された演算やデータ処理や制御シーケンスの発生をさせるものである場合も含まれる。以下では、データ処理部３２は、汎用コンピュータであるものとして説明する。
静磁場発生コイル１２は、Ｚ軸周りに巻回される螺旋コイルに静磁場電源２２から供給される電流を流して誘導磁場を発生させ、ボアにＺ軸方向の静磁場を発生させるものである。このボアに形成される静磁場の均一性の高い領域に被検者２の関心領域を設定することになる。ここで、静磁場コイル１２は、より詳しくは、たとえば、４個の空芯コイルから構成され、その組み合わせで内部に均一な磁界を作り、被検者２の体内の所定の原子核、より特定的には水素原子核のスピンに配向性を与える。
傾斜磁場発生コイル１４は、Ｘコイル、Ｙコイル及びＺコイル（図示省略）から構成され、円筒形状を示す静磁場発生コイル１２の内周面に設けられる。
これらＸコイル、Ｙコイル及びＺコイルは、それぞれＸ軸方向、Ｙ軸方向及びＺ軸方向を順番に切り替えながら、ボア内の均一磁場に対し傾斜磁場を重畳させ、静磁場に強度勾配を付与する。Ｚコイルは励起時に、磁界強度をＺ方向に傾斜させて共鳴面を限定し、Ｙコイルは、Ｚ方向の磁界印加の直後に短時間の傾斜を加えて検出信号にＹ座標に比例した位相変調を加え（位相エンコーディング）、Ｘコイルは、続いてデータ採取時に傾斜を加えて、検出信号にＸ座標に比例した周波数変調を与える（周波数エンコーディング）。
この重畳される傾斜磁場の切り替えは、制御シーケンスに従って、Ｘコイル、Ｙコイル及びＺコイルにそれぞれ異なるパルス信号が送信部２４から出力されることにより実現される。これにより、ＮＭＲ現象が発現する被検者２の位置を特定することができ、被検者２の画像を形成するのに必要な三次元座標上の位置情報が与えられる。
ここで、上述のように、３組の直交する傾斜磁場を用いて、それぞれにスライス方向、位相エンコード方向、および周波数エンコード方向を割り当ててその組み合わせにより様々な角度から撮影を行える。たとえば、Ｘ線ＣＴ装置で撮像されるものと同じ方向のトランスバーススライスに加えて、それと直交するサジタルスライスやコロナルスライス、更には面と垂直な方向が３組の直交する傾斜磁場の軸と平行でないオブリークスライス等について撮像することができる。
ＲＦ照射部１６は、制御シーケンスに従って信号送信部３３から送信される高周波信号に基づいて、被検者２の関心領域にＲＦ（Radio Frequency）パルスを照射するものであ
る。
なお、ＲＦ照射部１６は、図１において、磁場印加機構１１に内蔵されているが、寝台１８に設けられたり、あるいは、受信コイル２０と一体化されていてもよい。
受信コイル２０は、被検者２からの応答波（ＮＭＲ信号）を検出するものであって、このＮＭＲ信号を高感度で検出するために、被検者２に近接して配置されている。
ここで、受信コイル２０には、ＮＭＲ信号の電磁波がそのコイル素線を切ると電磁誘導に基づき微弱電流が生じる。この微弱電流は、信号受信部２８において増幅され、さらにアナログ信号からデジタル信号に変換されデータ処理部３２に送られる。
すなわち、静磁界にＺ軸傾斜磁界を加えた状態にある被検者２に、共鳴周波数の高周波電磁界を、ＲＦ照射部１６を通じて印加すると、磁界の強さが共鳴条件になっている部分の所定の原子核、たとえば、水素原子核が、選択的に励起されて共鳴し始める。共鳴条件に合致した部分（たとえば、被検者２の所定の厚さの断層）にある所定の原子核が励起され、スピンがいっせいに回転する。励起パルスを止めると、受信コイル２０には、今度は、回転しているスピンが放射する電磁波が信号を誘起し、しばらくの間、この信号が検出される。この信号によって、被検者２の体内の、所定の原子を含んだ組織を観察する。そして、信号の発信位置を知るために、ＸとＹの傾斜磁界を加えて信号を検知する、という構成になっている。
画像処理部４８は、記憶部３６に構築されているデータに基づき、励起信号を繰り返し与えつつ検出信号を測定し、１回目のフーリエ変換計算により、共鳴の周波数をＸ座標に還元し、２回目のフーリエ変換でＹ座標を復元して画像を得て、表示部３８に対応する画像を表示する。
たとえば、このようなＭＲＩシステムにより、上述したＢＯＬＤ信号をリアルタイムで撮像し、制御部４２により、時系列に撮像される画像について、後に説明するような解析処理を行うことで、安静時機能的結合ＭＲＩ（ｒｓ−ｆｃＭＲＩ）の撮像を行うことが可能となる。
図２は、データ処理部３２のハードウェアブロック図である。
データ処理部３２のハードウェアとしては、上述のとおり、特に限定されないが、汎用コンピュータを使用することが可能である。
図２において、データ処理部３２のコンピュータ本体２０１０は、メモリドライブ２０２０、ディスクドライブ２０３０に加えて、演算装置（ＣＰＵ）２０４０と、ディスクドライブ２０３０及びメモリドライブ２０２０に接続されたバス２０５０と、ブートアッププログラム等のプログラムを記憶するためのＲＯＭ２０６０とに接続され、アプリケーションプログラムの命令を一時的に記憶するとともに一時記憶空間を提供するためのＲＡＭ２０７０と、アプリケーションプログラム、システムプログラム、およびデータを記憶するための不揮発性記憶装置２０８０と、通信インタフェース２０９０とを含む。通信インタフェース２０９０は、駆動部２１等と信号の授受を行うためのインタフェース部４４および図示しないネットワークを介して他のコンピュータと通信するためのネットワークインタフェース５０に相当する。なお、不揮発性記憶装置２０８０としては、ハードディスク（ＨＤＤ）やソリッドステートドライブ（ＳＳＤ：Solid State Drive）などを使用す
ることが可能である。不揮発性記憶装置２０８０が、記憶部３６に相当する。データ処理部３２のコンピュータ本体２０１０は、第１の分類器生成装置、及び／又は第２の分類器生成装置としても機能する。
ＣＰＵ２０４０が、プログラムに基づいて実行する演算処理により、データ処理部３２の各機能、たとえば、制御部４２、データ収集部４６、画像処理部４８の各機能が実現さ
れる。
データ処理部３２に、上述した実施の形態の機能を実行させるプログラムは、ＣＤ−ＲＯＭ２２００、またはメモリ媒体２２１０に記憶されて、ディスクドライブ２０３０またはメモリドライブ２０２０に挿入され、さらに不揮発性記憶装置２０８０に転送されても良い。プログラムは実行の際にＲＡＭ２０７０にロードされる。
データ処理部３２は、さらに、入力装置としてのキーボード２１００およびマウス２１１０と、出力装置としてのディスプレイ２１２０とを備える。キーボード２１００およびマウス２１１０が入力部４０に相当し、ディスプレイ２１２０が表示部３８に相当する。
上述したようなデータ処理部３２として機能するためのプログラムは、コンピュータ本体２０１０に、情報処理装置等の機能を実行させるオペレーティングシステム（ＯＳ）は、必ずしも含まなくても良い。プログラムは、制御された態様で適切な機能（モジュール）を呼び出し、所望の結果が得られるようにする命令の部分のみを含んでいれば良い。データ処理部３２がどのように動作するかは周知であり、詳細な説明は省略する。
また、上記プログラムを実行するコンピュータは、単数であってもよく、複数であってもよい。すなわち、集中処理を行ってもよく、あるいは分散処理を行ってもよい。
図３は、うつ症状の疾患ラベルに関連する安静時の機能的結合の相関を示す相関行列を抽出する手続きを示す概念図である。
図３に示すように、リアルタイムで測定した安静時のｆＭＲＩのｎ個（ｎ：自然数）の時刻分のｆＭＲＩデータから、各関心領域の平均的な「活動度」を算出し、脳領域間（関心領域間）の活動度の相関値を算出する。
ここで、関心領域としては、小脳領域を除く１３７領域を考えているので、相関行列における独立な非対角成分は、対称性を考慮すると、
（１３７×１３７−１３７）／２＝９３１６（個）
ということになる。なお、図３では、図の上では、相関は３４×３４のみを示している。
このような相関行列の要素の算出については、特に限定されないが、たとえば、以下のように実行することが可能である。
安静時脳活動データから、各参加者について、異なる関心領域間の機能的結合（Functional connectivity， ＦＣ）を計算する。ＦＣは安静時脳活動解析において一般的に利用される特徴量であり、異なる関心領域の時系列信号間のピアソンの相関係数により定義される。
まず、各関心領域に含まれる全ボクセルの平均信号の時系列を抽出する。
次に信号値からのノイズ除去のため、バンドパスフィルタをかけ、その後、９つの説明変数（全脳，白質，脳脊，髄液の平均信号及び６つの体動補正パラメータ）による回帰を行う。
回帰後の残差系列を機能的結合と関連する時系列信号値と考え、異なるＲＯＩの時系列間でピアソンの相関係数を計算する。
関心領域については、Brain Sulci Atlas （BAL）に含まれる１３７の関心領域を用い
る。これら１３７の関心領域間の機能的結合ＦＣを特徴量として用いる。
ここで、Brain Sulci Atlas （BAL）については、以下に開示がある。
文献：Perrot et al., Med Image Anal, 15(4), 2011
文献：Tzourio-Mazoyer et al., Neuroimage, 15(1), 2002
図４は、図３で説明したような相関行列から、バイオマーカーとなる第１の分類器を生成する過程を説明する概念図である。
以下、説明するように、第１の分類器は、たとえば、「うつ病」に対する診断マーカーとしての指標として機能し、疾患の重症度や進行度を医師が診断する際の支援のために使用することができる。したがって、第１の分類器の演算と結果の出力を実行する装置は、診断支援装置と呼ぶこともできる。
図４に示すように、ＭＲＩ撮像部２５で撮像された、健常群、及びうつ病患者群を含む参加者において測定されたｆＭＲＩのデータから、それぞれの参加者について後に説明するような手続により、データ処理部３２が、脳領域間（関心領域間）の活動度の相関行列の要素を導出する。
続いて、データ処理部３２により、参加者の疾患／健常ラベルを含む参加者の属性と相関行列とについて正則化正準相関解析により特徴抽出が行われる。「正則化」とは、一般に、機械学習や統計学において、誤差関数に、ハイパーパラメータにより重み付けされた正則化項を追加して、モデルの複雑度・自由度に抑制を加え、過学習を防ぐ方法である。なお、正則化正準相関解析の結果、説明変数についてスパース化も併せて実現される場合は、特に、スパース正準相関解析（Sparse Canonical Correlation Analysis : SCCA）と呼ぶ。以下では、具体例としては、スパース正準相関解析を行うものとして説明する。
そして、このようなスパース正準相関解析において、後述するように、ハイパーパラメータの値の調節により「疾患ラベル」とのみ結合する正準変数が存在するようにしたときに、これに対応する正準変数に接続する機能的結合ＦＣを抽出する。所定範囲でハイパーパラメータを変化させるとき、このような条件を満たす正準変数が存在する範囲で抽出される機能的結合ＦＣの和集合を「第１和集合」と呼ぶ。
さらに、データ処理部３２により、スパース正則化正準相関解析の結果得られた「第１和集合」を説明変数とし、たとえば、１個抜き交差検証（Leave-One-Out Cross Validation：ＬＯＯＣＶ）を行いながら、各交差検証のステップでスパースロジスティック回帰による判別分析を行い、すべての交差検証にわたって、説明変数として抽出された機能的結合ＦＣの和集合を「第２和集合」と呼ぶ。
最後に、全ての参加者のデータに対して、「第２和集合」を説明変数として、目的変数である「疾患ラベル」に対して、スパースロジスティック回帰による判別分析を行って、第１の分類器が生成される。
図５は、図４で説明したような第１の分類器の生成処理および生成された第１の分類器による判別処理を行うための機能ブロック図である。
まず、不揮発性記憶装置２０８０には、健常者とうつ病患者を含む複数の参加者の各々の脳内の複数の所定領域における脳活動を示す信号をＭＲＩ装置により時系列で予め測定した信号の情報である参加者ｒｓ−ｆｃＭＲＩ測定データ３１０２と、ＭＲＩ測定データを計測した各参加者に関連づけられた、複数の人属性情報３１０４とが格納されている。
ここで、「人属性情報」とは、参加者個人を特定するための「人特性情報」と、各参加者に対する測定条件を特定するための「測定条件情報」とを含む。
「人特性情報」とは、参加者に対する疾患ラベル、年齢、性別、投薬履歴などの情報を意味する。
「測定条件情報」とは、参加者を測定した測定サイト情報（測定した施設および／または測定した装置を特定する情報を含む）や、測定中に、開眼／閉眼のいずれで計測したか、測定磁場強度などの測定条件を特定する条件である。
演算装置２０４０は、参加者ｒｓ−ｆｃＭＲＩ測定データ３１０２および対応する人属性情報３１０４とに基づいて、疾患ラベルに対する分類器を生成する処理を行う。
相関行列算出部３００２は、参加者ｒｓ−ｆｃＭＲＩ測定データ３１０２から複数の所定領域間の脳活動の機能的結合の相関行列を参加者ごとに算出する。算出された機能的結合の相関行列のデータは、参加者ごとに、不揮発性記憶装置２０８０に、機能的結合の相関行列のデータ３１０６として格納される。
第１特徴選択部３００４は、複数の参加者から抽出されたＫ個（Ｋ：２以上の自然数）の異なる部分集合に対して、順次１つの部分集合を選択し、選択された部分集合を除く（Ｋ−１）個の部分集合に対して、複数の属性情報と相関行例の要素とに対するスパース正準相関分析により、複数の人属性情報のうちの特定の属性情報、たとえば、疾患ラベルのみに対応する正準変数と接続する相関行例の要素を抽出する。さらに、第１特徴選択部３００４は、順次選択される部分集合にわたって、抽出された相関行例の要素の和集合である第１和集合を取得し、不揮発性記憶装置２０８０に、第１の機能的結合和集合データ３１０８として格納する。なお、「第１の機能的結合和集合データ」は、機能的結合の相関行列のデータ３１０６のうち、第１の和集合に対応する要素を特定するためのインデックスであってもよい。
第２特徴選択部３００６は、複数の参加者のうち上記Ｋ個の部分集合を除く残りの参加者をテスト集合とし、テスト集合を異なるＮ個のグループに分けるとき、複数の参加者からＮ個のグループうちの選択された１つのグループを除く参加者の集合を対象として、第１和集合に基づいて、上記特定の属性情報（たとえば、疾患ラベル）を推定するためのテスト分類器をスパースロジスティック回帰により算出し、スパース化に伴いテスト分類器の説明変数となる相関行例の要素を抽出する。第２特徴選択部３００６は、さらに、Ｎ個のグループうちから１つのグループを順次選択しつつ、特徴抽出を繰り返して、テスト分類器の説明変数として抽出された相関行例の要素の和集合である第２和集合を取得して、不揮発性記憶装置２０８０に、第２の機能的結合和集合データ３１１０として格納する。なお、「第２の機能的結合和集合データ」も、機能的結合の相関行列のデータ３１０６のうち、第２の和集合に対応する要素を特定するためのインデックスであってもよい。
なお、Ｎ個グループの各々の要素の個数は１個であってもよい。
分類器生成部３００８は、第２和集合を説明変数として、特定の属性情報（たとえば、疾患ラベル）を推定するための第１の分類器をスパースロジスティック回帰により算出する。分類器生成部３００８は、生成された第１の分類器を特定するための情報を、不揮発性記憶装置２０８０に、分類器データ３１１２として格納する。
判別処理部３０１０は、分類器データ３１１２により特定される分類器に基づいて、入
力データに対して判別処理を行う。
なお、上記の説明では、分類器生成部３００８は、第２和集合を説明変数として、分類器を生成するものとして説明したが、たとえば、分類器生成部３００８は、直接、第１和集合を説明変数として、第１の分類器を生成してもよい。ただし、後に説明するように、次元の低減や、汎化性能の観点からは、第２和集合を説明変数とすることが、望ましい。
なお、第１の分類器の生成にあたっては、正則化法（たとえば、Ｌ１正則化やＬ２正則化）を用いたロジスティック回帰である正則化ロジスティック回帰を用いることができ、より具体的には、たとえば、上述した「スパースロジスティック回帰」を使用することもできる。また、第１の分類器の生成には、サポートベクトルマシンや、ＬＤＡ(Linear discriminant analysis)などを使用してもよい。以下では、スパースロジスティック回帰を例にとって説明することにする。
後に説明するように、第２特徴選択部３００６による特徴選択処理と並行して、Ｎ個のグループうちから除かれる１つのグループを順次選択しつつ、第２特徴量選択部３００６により算出されたテスト分類器により、除かれたグループをテストサンプルとして判別結果を算出することにより交差検証を行ってもよい。
このように説明変数の次元の縮小を、入れ子構造となった特徴選択手続きにより行うことにより、第２特徴選択部３００６の行う処理において、ほとんどの参加者のデータを使用する一方、時間的に効率的な次元の数の縮小を実行できる。
さらに、第２特徴選択部３００６の処理におけるテスト集合を、次元の数を縮小するために使用されるデータセットとは、独立となるようにしておくことにより、極端に楽観的な結果が得られることを回避することができる。
図６は、バイオマーカーとなる分類器を生成するためにデータ処理部３２が行う処理を説明するためのフローチャートである。
以下、図４で説明した処理を、図６を参照しながら、より詳しく説明する。
安静時のｆＭＲＩデータから導き出される脳領野間の結合と参加者の疾患の判別ラベル（疾患ラベル）にもとづいてバイオマーカーを作成する場合の最大の問題点は、データ数にくらべてデータの次元が圧倒的に多いことである。そのため、正則化を行うことなしに、疾患ラベル（ここでは、参加者が疾患であるか健常であるかを示すラベルを「疾患ラベル」と呼ぶ）の予測のための分類器の学習をデータセットを用いておこなうと、オーバーフィット（過剰適合）することで未知のデータに対する予測性能が著しく低下する。
ここで、一般に、機械学習において、観測データをより少ない数の説明変数で説明するように行われる処理のことを「特徴選択（または特徴抽出）」と呼ぶ。本実施の形態では、「脳領域間（関心領域間）の活動度の複数個の相関値（複数個の結合）」のうち、対象となる疾患ラベルの予測のための第１の分類器の機械学習において、より少ない相関値により第１の分類器を構成できるように特徴選択（特徴抽出）をすること、すなわち、説明変数として、より重要な相関係数を選択することを「縮約表現を抽出する」と呼ぶこともある。
そして、本実施の形態では、特徴抽出法として、正則化法を用いる。このように、正準相関分析において、正則化とともにスパース化を実行し、説明変数としてより重要なものを残すように処理するという意味で、このような処理を「スパース正準相関分析」と呼ぶ
。より具体的には、たとえば、スパース化を併せて実現するための正則化法としては、以下に説明するような「Ｌ１正則化」という正準相関分析のパラメータの絶対値の大きさにペナルティーを課す手法を用いることが可能である。
すなわち、図６を参照して、データ処理部３２は、入力部４０からの処理開始の入力等により、分類器生成の処理が開始されると（Ｓ１００）、記憶部３６から各参加者についてのＭＲＩ測定データを読み出し（Ｓ１０２）、スパース正準相関分析（SCCA）により第１の特徴抽出の処理を行う（Ｓ１０４）。
以下では、ステップＳ１０４の処理を「内ループ特徴選択」と呼ぶ。
ここで、ステップＳ１０４での内ループ特徴選択の処理を説明するために、「スパース正準相関分析」と、「内ループ特徴選択」について、以下、それぞれ説明する。
（スパース正準相関分析）
以下では、スパース正準相関分析として、Ｌ１正則化正準相関分析について説明する。このＬ１正則化正準相関分析については、以下の文献に開示がある。
文献：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）では、データ対
ｘ1およびｘ2を考える際に、それぞれの変数ｘ1およびｘ2は、平均０、標準偏差１を持つように標準化される。
一般的に言えば、正準相関分析(CCA)の使用によって、ペアになった観測量の間の潜在
している関係を識別することができる。
具体的には、正準相関分析(CCA)では、対となる射影された変数（正準変数）に最大の
相関があるように、投影ベクトルを見つけ出す。
これに対して、正準相関分析にＬ１正則化を導入した場合、すなわち、Ｌ１ノルム正規化によるスパースＣＣＡ、すなわち、Ｌ１−ＳＣＣＡを使用する場合は、以下のような最適化問題を解くことになる。
以下のような変数の組を考える。
ここでも、行列Ｘ1およびＸ2を構成する列については、それぞれ平均が０で分散が１となるように正規化されているものとする。
そうすると、Ｌ１−ＳＣＣＡは、以下の式（１）のように定式化することができる。
ここで、ハイパーパラメータλ1およびλ2が、重みパラメータw1およびw2の（w1およびw2により対応する変数がスパース化されることから「スパース投影ベクトル」と呼ぶ）スパース化の度合いをそれぞれ示す。
具体的に、本実施の形態の例としては、人属性情報および機能的な結合ＦＣの潜在している関係を識別するために、上記の変数に相当する２つのデータ行列を構築するものとする。
以下では、参加者には、うつ病の患者群（疾患ラベルを「MDD（Major depressive disorder）」で表す）と健常者対照群（疾患ラベルを「HC（Healthy control）」で表す）第
１のデータ行列Ｘ1の行は、１人の参加者の人属性情報（人特性情報および測定条件情報
）を表しており、たとえば、この特性情報や測定条件情報は、以下のものを含む：
ｉ）疾患ラベル（ＭＤＤあるいはＨＣ），
ｉｉ）サイト情報（参加者の脳活動がどこで計測されたか、サイトＡ，Ｂ，Ｃ）
ｉｉｉ）年齢
ｉｖ）性別
ｖ）撮像条件（開眼あるいは閉眼）
ｖｉ）薬剤投与の状態１（抗精神病薬）
ｖｉｉ）薬剤投与の状態２（抗うつ薬）
ｖｉｉｉ）薬剤投与の状態３（精神安定剤）
具体的には、人属性情報データ行列Ｘ1の列数は、１０であり、つまりｐ1＝１０である。
最初の列は、１（＝ＭＤＤ）あるいは０（＝ＨＣ）を含んでいる。
次の３つの列は、脳活動の計測された施設（場所）を示しており、この場合は、３つの異なるサイトで計測が行われたことに相当して、［１００］（サイトＡ）、［０１０］（サイトＢ）あるいは［００１］（サイトＣ）のいずれかを含んでいる。
第５列は、参加者の年齢の値を含んでおり、第６列は、参加者の性別の情報であって、１（男性）あるいは０（女性）のいずれかの値を含んでおり、第７列は、測定条件のうち、測定中に参加者が、目を開いていたか、あるいは、目を閉じていたかを示し、１（開眼）あるいは０（閉眼）の値を含んでいる。また、最後の３つの列は、３つの薬剤の投薬履歴に関するステータス情報を含んでおり、それぞれの列は、１(薬物治療あり)あるい０（薬物治療履歴無し)の値を含んでいる。
なお、人特性情報および測定条件情報は、上記のようなものに限定されず、たとえば、参加者の他の薬剤の投薬履歴の情報などの他の特性や、他の測定条件、たとえば、ＭＲＩ装置の印加磁場の大きさなどの条件の情報を含んでいてもよい。
第２のデータ行列Ｘ2は、参加者の機能的結合の相関（ＦＣ）を表わす相関行列の非対
角の下三角の部分の要素を行ベクトル形式で表現したものであるとする。
そして、Ｌ１−ＳＣＣＳが、行列のペアＸ1およびＸ2に適用され、それから、スパース投影ベクトルw1およびw2が導き出される。
そして、上述のとおり、ハイパーパラメータλ1およびλ2を所定の値に設定したときに、スパース化により第１のデータ行列Ｘ1が、スパース投影ベクトルw1により、特定の人
属性情報、この場合は、「疾患ラベル」のみとなる正準変数に投影される条件が満たされるようにする。この場合に、対応する第２のデータ行列Ｘ2において、疾患ラベルにのみ
関連する機能的結合の相関行列要素を特定するためのインデックス（要素）を特定するために、スパース投影ベクトルw2を使用する。
すなわち、内ループ特徴選択では、Ｌ１−ＳＣＣＡのハイパーパラメータλ1およびλ2は、たとえば、０．１の大きさのステップで、０．１と０．９の間で独立して変更される。
ただし、ハイパーパラメータλ1およびλ2を可変とする範囲および変化させるステップの大きさは、このような例に限定されるものではない。
Ｌ１−ＳＣＣＡの処理に対して、「疾患」ラベルにのみ接続した正準変数が存在するようなハイパーパラメータλ1およびλ2の範囲を見つけ出す。
導出されたスパース投影ベクトルw2のゼロでない要素によって規定される部分空間へ元の相関行列要素を投影することを以下のようにして表現する。
ここで、投影ベクトルw2のｋ番目のゼロでない要素のインデックスを示すために変数ｉkを定義する。ここで、１≦ｋ≦ｍであり、ｍはゼロでないエレメントの数を示す。
そして、以下のような部分空間への投影マトリックスＥを考える。
最後に、元の相関行列要素ベクトルｘ2を以下のように投影することにより部分空間ｚ
∈Ｒmへのベクトルが導かれる。
結果として、疾患ラベル（ＭＤＤ／ＨＣ）に関連する特定の数だけの特徴量（相関行列の要素）を選択することができる。
診断のラベルとのみ関係のある正準変数に相当する相関行列要素を選ぶことによって、分類に本質的な相関行列要素が選択される。
（内ループ特徴抽出）
図７は、内ループ特徴選択を説明するための概念図である。
図６のステップＳ１０４の内ループ特徴選択では、図７に説明するようにして、参加者の集合をＫ個（Ｋ：２以上の自然数）に分割して、（Ｋ−１）個の部分集合のうち、１つの集合を除いた残りの部分集合について、上述したＬ１−ＳＣＣＡを実行することを、（Ｋ−１）回繰り返す。
すなわち、図７（ａ）に示すように、参加者の集合のうちの（Ｋ−１）／Ｋ（たとえば、Ｋ＝９であれば８／９）に相当する部分集合（これを、以下「内ループ訓練データ」と呼ぶ）が、このような内ループ特徴選択に使用される。Ｋ個の部分集合のうちの残りの１つは、後に説明する外ループ特徴選択において、訓練対象となるテストデータを含む「テストプール」として使用されるため、内ループ特徴選択では使用されない。
図７（ｂ）に示すように、（Ｋ−１）個の部分集合のうち、１つの集合（図７（ｂ）で斜線で示す）を除いた残りの部分集合について、ハイパーパラメータλ1およびλ2を所定の範囲で、所定ステップずつ変化させ、Ｌ１−ＳＣＣＡを実行する。特定のハイパーパラメータλ1およびλ2の範囲で、「診断」ラベルとのみ関連がある正準変数に関係している機能的結合要素ＦＣを特徴量として抽出する。
このような抽出処理を、（Ｋ−１）個の部分集合のうち、除くことになる１つの集合を順次変えた部分集合について、反復する。
各反復において、抽出された機能的結合の相関行列の要素ＦＣの和集合を「内ループで選択された機能的結合要素ＦＣの和集合」（第１の和集合）とする。
このような処理を行うことで、攪乱変数ＮＶに相当する、異なる撮像サイトでの人属性情報や異なる撮像サイトでの撮像条件の相違によって引き起こされた不適当な影響は縮小されることになる。
この手続きは、後に説明するように、たとえば、日本国内の複数の撮像サイトからのＭＲＩ計測データに基づいて生成した第１の分類器を、外国、たとえば、米国にまで汎化するロバストに構築するのに役立つことになる。
再び、図６に戻って、続いて、データ処理部３２は、内ループ特徴選択の結果に基づいて、スパースロジスティック回帰を用いて、第２の特徴抽出の処理を実行する（Ｓ１０６）。
以下では、ステップＳ１０６の処理を「外ループ特徴選択」と呼ぶ。
ここで、ステップＳ１０６での外ループ特徴選択の処理を説明するために、「スパースロジスティック回帰」と、「外ループ特徴選択」について、以下、それぞれ説明する。
（スパースロジスティック回帰）
スパースロジスティック回帰は、ロジスティック回帰分析をベイズ推定の枠組みに拡張した手法であり、特徴ベクトルの次元圧縮を判別のための重み推定と同時に行う手法である。データの特徴ベクトルの次元数が非常に高く、不要な特徴量が多く含まれている場合に有用である。不要な特徴量に対しては線形判別分析における重みパラメータをゼロにし（すなわち、特徴選択を行い）、判別に関連するごく少数の特徴量だけを取り出す（スパース性）。
スパースロジスティック回帰では得られた特徴データに対して、分類するクラスへの所属する確率ｐをクラス毎に求め、最大値を出力したクラスに割り当てる。ｐはロジスティック回帰式によって出力される。重みの推定はＡＲＤ（Automatic Relevance determination）によって行われ、クラス判別への貢献が少ない特徴量は、重みが０に近づくことで
計算から除外される。
すなわち、上述のＬ１正則化ＣＣＡを用いて抽出された特徴量の第１の和集合を入力として、次の階層ベイズ推定にもとづく第１の分類器を用いて、疾患ラベルの予測を行う。
ここで、上述した式（２）で抽出された特徴入力ｚ（選択されたFC）から診断のラベルが疾患（ここでは、自閉症との診断）となる確率を予測するために、分類器としてロジスティック回帰を使用する。
ここで、ｙは、診断クラス・ラベルを表わし、つまり、ｙ＝１はＭＤＤクラスを、ｙ＝０はＨＣクラスをそれぞれ示す。
また、以下のｚハット（文字の頭部に＾が付けられたものを、「ハット」と呼ぶ）は、拡張入力を備えた特徴ベクトルである。
ここで、特徴ベクトルｚは、１人の参加者の安静時ＭＲＩサンプルの結合相関行列から式（２）に従って、抽出される。
拡張入力“１”を使用することは、第１の分類器に対する一定の(バイアス)入力を導入するための、標準的なアプローチである。
以下のようなθは、ロジスティック関数のパラメータベクトルである。
この際のパラメータθの分布を以下のような正規分布に設定する。
さらに、パラメータｗの分布のハイパーパラメータαの分布を以下のように設定することにより、階層ベイズ推定を行うことによりそれぞれのパラメータの分布を推定する。
ここで、ａ0, ｂ0はハイパーパラメータのガンマ分布を決定するパラメータである。αはベクトルθの正規分布の分散を表現するパラメータベクトルで、そのベクトルの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.
（外ループ特徴選択）
図１１は、外ループ特徴選択を説明するための概念図である。
図６のステップＳ１０６の外ループ特徴選択では、図１１に説明するようにして、内ループ特徴選択の際に、参加者の集合をＫ個（Ｋ：２以上の自然数）に分割したうち、内ループ特徴選択に使用した（Ｋ−１）個の部分集合を除く残りの１つの部分集合を、訓練対象となるテストデータを含む「テストプール」として使用する。
すなわち、図１１に示すように、たとえば、このテストプールをＬ個のグループに分割し、このＬ個のグループから順次１つのグループ（図１１の点線四角内で、斜線で示す）を選択し、参加者の集合この選択されたグループを除く参加者について、第１の和集合である機能的結合の相関行列の要素ＦＣを説明変数として、疾患ラベルを予測する第１の分類器をスパースロジスティック回帰により生成する。このときに、スパースロジスティック回帰により「機能的結合の相関行列の要素ＦＣ」が、さらに選択される。これを「外ループで選択された機能的結合の相関行列の要素ＦＣ」と呼ぶ。
このようにして、Ｌ個のグループから順次１つのグループを選択して、このような処理をＬ回反復することで、抽出された機能的結合の相関行列の要素ＦＣの和集合を「外ループで選択された機能的結合要素ＦＣの和集合」（第２の和集合）とする。
このように、外ループ特徴抽出において使用されるテスト集合は、内ループ特徴抽出において次元数を縮小するために使用されたデータセットから常に独立している。
なお、ここでは、一例として、特に、Ｌ個のグループの各々は、１人の参加者を含むものとする。この場合は、上記反復処理において、除かれた１人の参加者について、生成された第１の分類器で予測処理を行い、その除かれた１人の参加者についての疾患ラベルと予測結果との誤差を累積処理することは、いわゆる１個抜き交差検証（Leave-One-Out Cross Validation：ＬＯＯＣＶ）を行うことに相当する。したがって、図１１の処理において、各反復処理で、上述の誤差を算出し、これを累積してＬ回について平均することで、生成された第１の分類器の交差検証、言い換えると、汎化能力の評価を行うことも可能である。
（第１の分類器の生成）
再び図６に戻って、続いて、データ処理部３２は、外ループ特徴選択の結果（第２の和集合）に基づいて、スパースロジスティック回帰により、分類器の生成処理を実行する（Ｓ１０８）。
図１３は、ステップＳ１０８での分類器の生成処理の概念を示す図である。
図１３に示すように、ステップＳ１０８では、第１の和集合から外ループ特徴選択により抽出された第２の和集合を説明変数として、全ての参加者についてのスパースロジスティック回帰により、分類器が生成される。
生成された分類器を特定するための分類器データ（関数形およびパラメータに関するデータ）３１１２は、不揮発性記憶装置２０８０に格納され、後に、以上のような訓練に使用されたのとは異なるＭＲＩ測定データ（テストデータ）が入力された際に、テストデータに対する疾患ラベルの推定を実行する際の判別処理に使用される。
すなわち、事前の医師の診断に基づいて、健常群、うつ病患者群に分けられた参加者の脳領域間（関心領域間）の活動度の相関（結合）を測定し、測定結果に対する機械学習により、これらとは異なる新たな参加者についてのテストデータがうつ症状または健常のいずれに相当するかを判別するように生成された分類器が、うつ症状のバイオマーカーとして機能する。
このとき、バイオマーカーの出力である「疾患ラベル」としては、分類器がロジスティック回帰により生成されることから、疾患である確率（または健常である確率）を含んでもよい。たとえば、「疾患である確率は、○○％」といような表示として出力する。このような確率を、「疾患マーカー」として使用することが可能である。
また、分類器の出力となる属性については、必ずしも、疾患の判別に限られず、他の属性についての出力であってもよい。この場合も、その属性について、いずれのクラスに属するかという離散的な判別結果を出力してもよいし、その属性について、あるクラスに属する連続的な値、たとえば確率を出力してもよい。
すなわち、分類器の学習（作成）においては、うつ症状のバイオマーカー作成のため，ｒｓ−ｆｃＭＲＩデータを入力として，上述の内ループ特徴選択および外ループ特徴選択により特徴抽出を行い、抽出された特徴量を説明変数として生成された分類器によりうつ症状・健常ラベル判別を行う。
（内ループ特徴選択での処理）
図８は、内ループ特徴選択処理の概念を説明するための図である。
図９は、一例として、内ループ特徴選択を行った際の反復処理のうち、特定のハイパーパラメータλ1およびλ2に対する結果を示す図である。
図８を参照して、入れ子構造になっている特徴選択において、Ｌ１−ＳＣＣＡの反復のうちの１つを表わし、ここでは、正準変数が「疾患ラベル」にのみ接続されている。
図８では、正準変数が「疾患ラベル」にのみ接続されていることを、疾患ラベルがただ１つの正準変数w1 Tx1と、実線で接続されることで表している。
また、正準変数ｗ１ Ｔｘ１と正準変数ｗ２ Ｔｘ２とを接続する点線の上の記号ｃｊは、それぞれの正準変数間の相関係数を表す。
正準変数ｗ１ Ｔｘ１が「疾患ラベル」にのみ接続されているときに、対応する正準変数ｗ２ Ｔｘ２と接続する機能的結合の相関行列の非対角の下三角行列の要素が、特徴量として選択される。
図９では、例として、最初の外ループの中で、各人属性情報につき少なくとも１つの正準相関を生成する最も小さなハイパーパラメータの組合せを示す。
正準変数は、丸によって表わされる。
左の列の丸は人特性情報と測定条件に由来する、正準変数ｗ１ Ｔｘ１を示す。一方で、右の列の丸は、機能的結合（ＦＣ）に由来する正準変数ｗ２ Ｔｘ２を示す。
上述のとおり、正準変数を接続する点線の上の数字は、正準変数ｗ１ Ｔｘ１とｗ２ Ｔｘ２の間の相関係数を表わす。
人属性情報のラベルおよび正準変数ｗ２ Ｔｘ２の間の接続は実線または点線で表わされる。
（内ループ特徴選択での処理）
図１０は、内ループ特徴選択の処理を、より詳しく説明するためのフローチャートである。
内ループ特徴選択にあたっては、参加者の集合は、図７（ａ）に示したように、Ｋ個の部分集合に分けられ、このうちの特定の１個の部分集合を除く（Ｋ−１）個の部分集合が、使用される。
図１０を参照して、内ループ特徴選択の処理が開始されると、演算処理装置（ＣＰＵ）２０４０は、ハイパーパラメータλ1およびλ2を初期値である（λ1，λ2）＝（０．１，０．１）に設定する（Ｓ２００）。
続いて、変数ｉの値が１に設定され（Ｓ２０２）、ＣＰＵ２０４０は、ｉの値が、内ループ特徴選択の反復回数（Ｋ−１）を超えない場合（Ｓ２０４でＮ）は、内ループ訓練データのうち、i番目のデータブロックを除き、不揮発性記憶装置２０８０中に格納された
、参加者の人属性情報３１０４および機能的結合の相関行列データ３１０６に基づいて、スパース正準相関分析を実行する（Ｓ２０６）。
現在の（λ1，λ2）の値に組について、疾患ラベルとのみ結合する正準変数が存在する場合（Ｓ２０８でＹ）、ＣＰＵ２０４０は、疾患ラベルとのみ結合する正準変数に対応する機能的結合の相関行列の要素（ＦＣ）を抽出し、不揮発性記憶装置２０８０に、記憶する（Ｓ２１０）。ステップＳ２１０の処理の後、または、現在の（λ1，λ2）の値に組について、疾患ラベルとのみ結合する正準変数が存在しない場合（Ｓ２０８でＮ）は、ｉの値を１だけインクリメントして、処理は、ステップＳ２０４に復帰する。
したがって、ＣＰＵ２０４０は、（Ｋ−１）回にわたって、ステップＳ２０６〜Ｓ２１２の処理を繰り返す。
ステップＳ２０４において、ｉの値が、内ループ特徴選択の反復回数（Ｋ−１）を超え
た場合（Ｓ２０４でＹ）は、ＣＰＵ２０４０は、（λ1，λ2）について、所定のルールにしがたい、λ1およびλ2のいずれか一方を所定のステップ量だけ変化させる（Ｓ２１４）。ステップＳ２１６において、（λ1，λ2）の値が可変範囲内であるときは（Ｓ２１６でＹ）、ＣＰＵ２０４０は、処理をステップＳ２０２に復帰させる。
（λ1，λ2）の可能な全ての値の組について処理が終了すると（Ｓ２１６でＮ）、ＣＰＵ２０４０は、これまでの処理で抽出された機能的結合の相関行列の要素（ＦＣ）の和集合を求め、不揮発性記憶装置２０８０に、第１の機能的結合和集合データ３１０８として格納することで処理が終了する。
（外ループ特徴選択での処理）
図１２は、外ループ特徴選択の処理を、より詳しく説明するためのフローチャートである。
外ループ特徴選択にあたっては、参加者の集合は、図１１で示したように、Ｋ個の部分集合に分けられ、このうちの内ループ特徴選択では使用されなかった１個の部分集合がテストプールとして使用される。このテストプールに含まれる参加者の人数をＮＴ人とする。
また、以下では、一例として、１個抜き交差検証（ＬＯＯＣＶ）処理を行うものとして説明する。
図１２を参照して、外ループ特徴選択の処理が開始されると、演算処理装置（ＣＰＵ）２０４０は、変数ｉの値を１に設定し（Ｓ３００）、ＣＰＵ２０４０は、ｉの値が、外ループ特徴選択の反復回数ＮＴを超えない場合（Ｓ３０２でＮ）は、テストプールのうち、i番目のデータを除き、不揮発性記憶装置２０８０中に格納された、参加者の人属性情報
３１０４および機能的結合の相関行列データ３１０６に基づいて、スパースロジスティック回帰（ＳＬＲ）を実行して、テスト分類器を生成する（Ｓ３０４）。
続いて、ＣＰＵ２０４０は、生成されたテスト分類器により、除かれたi番目のデータ
に対するＳＬＲによる予測値を算出する（Ｓ３０６）。
さらに、ＣＰＵ２０４０は、テスト分類器を生成する際のスパース化により抽出された機能的結合ＦＣを、選択された特徴量として不揮発性記憶装置２０８０に格納し、変数ｉの値を１だけインクリメントする（Ｓ３１０）。
Ｓ３０４〜Ｓ３０８の処理がＮＴ回繰り返された後、変数ｉの値がＮＴを超えると（Ｓ３０２でＹ）、ＣＰＵ２０４０は、ＬＯＯＣＶとしての２乗誤差見積の算出する（Ｓ３１２）。
続いて、ＣＰＵ２０４０は、ＮＴ回の繰り返し処理により、抽出された機能的結合ＦＣの和集合（第２の和集合）を求め、第２の機能的結合和集合データ３１１０として、不揮発性記憶装置２０８０に格納して処理が終了する。
（分類器の生成処理）
図１３は、最終的に第１の分類器を生成する手続きを説明するための概念図である。
図１３に示すように、分類器生成処理では、内ループ特徴選択処理および外ループ特徴選択処理により、最終的に抽出された第２の機能的結合和集合データ３１１０を説明変数とし、参加者全員について、スパースロジスティック回帰によって、第１の分類器を生成
する。生成された第１の分類器を特定するための情報（ロジスティック関数のパラメータベクトルθ等）は、不揮発性記憶装置２０８０に記憶される。
上述したように、ここまで説明した手続きでは、「うつ病」を疾患ラベルとした場合は、表１に記載の１２個の機能的結合が、最終的に抽出された第２の機能的結合和集合データ３１１０の説明変数として使用され、第１の分類器が生成される。
ただし、ロジスティック回帰式で表現される分類器の生成には表１に記載の全ての機能的結合を使用してもよいが、少なくとも機能的結合識別番号１及び機能的結合識別番号２の一方、あるいは両方を選択して使用してもよい。また、寄与度が高い複数の機能的結合を選択してもよい。このような１２個の機能的結合からの分類器に使用する前記機能的結合の選択は、データ処理部３２または後述するコンピュータ３００のＣＰＵが行ってもよいが、ヒトが行ってもよい。
（判別処理）
図１４は、生成した第１の分類器を用いて被検者のｒｓ−ｆｃＭＲＩのデータから疾患ラベルを判別するためにデータ処理部３２が行う処理を説明するためのフローチャート図である。
ＣＰＵ２０４０は、ステップＳ４０１において、ＭＲＩ撮像部２５から安静時の被検者ｒｓ−ｆｃＭＲＩ測定データ３１１３をインタフェース部４４を介して取得する。
ＣＰＵ２０４０は、ステップＳ４０２〜Ｓ４０３おいて上記２．に記載の前処理を行い、表１に記載の機能的結合の全て、あるいは一部について相関行列の要素を抽出する。
ステップＳ４０４おいて、ＣＰＵ２０４０は、抽出した相関行列の要素に対してピアソン相関係数を求める。ステップＳ４０５〜Ｓ４０６おいて、ＣＰＵ２０４０は、前記相関係数を分類器に入力し、少なくとも１つの機能的結合についてうつ症状の疾患ラベルを判別するための指標値を生成する。あるいは、望ましくはＣＰＵ２０４０は、前記相関係数を分類器に入力し、複数の機能的結合についての関連重み付け和をうつ症状の疾患ラベルを判別するための指標値として生成する。
次にＣＰＵ２０４０は、ステップＳ４０７において指標値を基準値と比較する。例えば数６で表されるような数５によるロジスティック回帰の式では、うつ症状の疾患ラベルと健常者のラベルの閾値は０となるように設定されているため、基準値は０となる。たとえば、「関連重み付け和」をバイオマーカーの指標値として使用する場合は、第１の分類器の感度や特異度を調整するために、基準値を０とは異なる値に設定してもよい。たとえば、基準値を０よりも小さくすれば、「うつ病」の疾患レベルと判別される感度は向上するが、特異度は低下する。上述した分類器の判別結果を、後続する他の情報に基づく医師による診断のための参考情報（支援情報）として使用することを意図する場合は、感度の向上を優先させるという設定もあり得る。
ステップＳ４０８にお
いてＣＰＵ２０４０は、前記指標値が０よりも高い場合（「ＹＥＳ：Ｙ」）には、前記被検者がうつ症状のラベルであると判別することができる（Ｓ４０９）。また、前記指標値が０よりも低い場合（「ＮＯ：Ｎ」）には、前記被検者がうつ症状のラベルに相当すると判別することができる（Ｓ４１０）。
３−２．判別装置２
実施の形態１では、脳活動計測装置（ｆＭＲＩ装置）で、１つの測定場所で計測された
脳活動データを計測し、この脳活動データに基づいて、分類器の生成および分類器による疾患ラベルの推定（予測）を、同一のコンピュータによる処理または分散処理により行う構成を説明した。
ただし、ｉ）分類器を機械学習により訓練するための脳活動データの計測（データ収集）、ｉｉ）分類器の機械学習による生成処理および特定の被験者についての分類器による疾患ラベルの推定（予測）の処理（推定処理）、ｉｉｉ）上記特定の参加者についての脳活動データの計測（対象者の脳活動計測）を、それぞれ、異なる施設で分散して実行する構成とすることも可能である。
すなわち、健常群およびうつ病患者群のデータは、ＭＲＩ装置１０自身で測定される場合に限られず、他のＭＲＩ装置において測定されたデータをも総合して、分類器の生成が行われてもよい。また、より一般には、データ処理部３２は、ＭＲＩ装置の制御を実行するためのコンピュータである必要は必ずしもなく、１つまたは複数のＭＲＩ装置からの測定データを受け取って、分類器の生成処理および生成された分類器により、判別処理を行うことに特化したコンピュータであってもよい。
図１５は、本発明の実施の形態１の別実施形態である、データ収集、推定処理および対象者の脳活動計測を、分散して処理する場合の一例を示す機能ブロック図である。
図１５を参照して、サイト１００．１〜１００．Ｎは、うつ病患者群、健常者群を含む参加者のデータを脳活動計測装置により計測する施設であり、管理サーバ２００は、サイト１００．１〜１００．Ｎからの計測データを管理する。
判別装置２に相当するコンピュータ３００は、管理サーバ２００に格納されたデータから分類器を生成する。
ＭＲＩ装置４１０は、コンピュータ３００上の分類器の結果を利用する別サイトに設けられており、特定の被検者について脳活動のデータを計測する。
コンピュータ４００は、ＭＲＩ装置４１０が設けられる別サイトに設置され、ＭＲＩ装置４１０の測定データから特定被検者の脳の機能的結合の相関データを算出し、機能的結合の相関データをコンピュータ３００に送信して、返信されてくる分類器による判別結果を利用する。
サーバ２００は、サイト１００．１〜１００．Ｎから送信されてくるうつ病患者群および健常者群の参加者ＦＭＲＩ測定データ３１０２と、参加者ＦＭＲＩ測定データ３１０２に関連付けられた参加者の人属性情報３１０４とを格納し、コンピュータ３００からのアクセスに従って、これらのデータをコンピュータ３００に送信する。
コンピュータ３００は、通信インタフェース２０９０を介して、サーバ２００からの参加者ＦＭＲＩ測定データ３１０２および被験者の人属性情報３１０４を受信する。
なお、サーバ２００、コンピュータ３００、コンピュータ４００のハードウェアの構成は、基本的に、図２で説明した「データ処理部３２」の構成と同様であるので、その説明は繰り返さない。
図１５に戻って、相関行列算出部３００２、第１特徴選択部３００４、第２特徴選択部３００６、分類器生成部３００８および判別処理部３０１０、ならびに、機能的結合の相関行列のデータ３１０６、第１の機能的結合和集合データ３１０８、第２の機能的結合和
集合データ３１１０および分類器データ３１１２については、実施の形態１で説明したのと同様であるので、その説明は、繰り返さない。
ＭＲＩ装置４１０は、診断ラベルの推定対象となる被検者の脳活動データを計測し、コンピュータ４００の処理装置４０４０は、計測されたＭＲＩ測定データ４１０２を不揮発性記憶装置４１００に格納する。
さらに、コンピュータ４００の処理装置４０４０は、ＭＲＩ測定データ４１０２とＭＲＩ装置４１０の計測対象の被験者の人属性情報とに基づいて、相関行列算出部３００２と同様にして、機能的結合の相関行列のデータ４１０６を算出し、不揮発性記憶装置４１００に格納する。
コンピュータ４００のユーザから診断の対象となる疾患が指定され、当該ユーザの送信の指示に従い、コンピュータ４００が、機能的結合の相関行列のデータ４１０６をコンピュータ３００に送信する。これに応じて、判別処理部３０１０は、指定された診断ラベルについての判別結果を算出し、コンピュータ３００は、通信インタフェース２０９０を介して、コンピュータ４００に送信する。
コンピュータ４００では、図示しない表示装置などを介して、ユーザに対して、判別結果を知らせる。
このような構成とすることで、より多くの被験者について収集したデータに基づいて、より多くのユーザに対して、分類器による診断ラベルの推定結果を提供することが可能となる。
また、サーバ２００とコンピュータ３００とを別個の管理者が管理する形態とすることも可能で、その場合、サーバ２００にアクセスできるコンピュータを制限することで、サーバ２００に格納される被験者の情報のセキュリティを向上させることも可能となる。
さらに、コンピュータ３００の運営主体からみると、「分類器による判別のサービスを受ける側（コンピュータ４００）」に対して、分類器についての情報は、一切提供しなくても、「判別結果を提供するサービス」を行うことが可能となる。
３−４．分類器及び判別装置の用途
前記判別装置１及び２は、うつ症状に対する疾患ラベルを判別する他、うつ症状のレベルを判別するために、うつ症状の治療の効果の程度を判別した情報を生成するために、又はうつ症状が予め設定された複数のサブクラスに層別化されている場合に、又は被検者を前記サブクラスに層別化する際の支援情報を生成するために使用することができる。それぞれの使用方法については、後述する。
第１の分類器は、うつ症状を判別するための診断バイオマーカーとして使用することができる。また、第１の分類器は、うつ症状のレベルを判別するための診断バイオマーカーとして使用することができる。第１の分類器は、うつ症状の治療の効果を判別するための薬効マーカーとして使用することができる。第１の分類器は、被検者をうつ病のサブクラスに層別化するためのバイオマーカーとして使用することができる。
３−５．分類器を生成するためのコンピュータプログラムとうつ症状を判別するためのコンピュータプログラム
実施の形態１の別形態には、ステップＳ１００〜Ｓ１０８（具体的には、Ｓ２００〜Ｓ２１８とＳ３００〜Ｓ３１４を含む）を含む処理を実行させて、分類器生成処理の機能を
コンピュータに実行させるためのプログラムと、ステップＳ４０１〜Ｓ４１０を含む処理を実行させて、判別処理の機能をコンピュータに実行させるためのプログラムとを含む。また実施の形態１の別形態には、ステップＳ１００〜Ｓ１０８（具体的には、Ｓ２００〜Ｓ２１８とＳ３００〜Ｓ３１４を含む）と、ステップＳ４０１〜Ｓ４１０を含む処理を実行させて、判別装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
３−６．うつ症状の判別方法
実施の形態２は、うつ症状の判別方法に関する。前記判別方法は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程と、前記指標値が基準値を超えた場合に、前記被検者がうつ症状であると決定する工程を備える。具体的な工程は、上記Ｓ４０１〜Ｓ４１０に準ずるが、前記工程は、指標値を生成する工程を除いた工程の全部または一部をヒトが行ってもよい。
４．うつ症状のレベルの判別装置及び判別方法
４−１．判別装置
判別装置１及び判別装置２は、うつ症状のレベルを判別するために使用することができる。すなわち、実施の形態１の別形態であるうつ症状のレベルの判別装置は、演算装置と記憶装置とを有する。この演算装置は、記憶装置に格納されたプログラムに基づいて、前記第１の分類器により、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成し、前記指標値をあらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較し、前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する、という処理を実行する。この場合、判別装置１及び判別装置２の構成は、上記３．に記載した構成と同様であるため、説明は省略する。判別処理部３０１０は、分類器データ３１１２により特定される分類器に基づいて、入力データに対して判別処理を行い、前記被検者のうつ症状のレベルを判別する。
この場合、「うつ症状のレベル」の判別には、たとえば、指標値である「関連重み付け和」の値が取り得る値を、予め、複数のレベルの基準範囲に分類しておくことで、「うつ症状のレベル」に対応した判別結果を出力することが可能となる。
あるいは、この場合、分類器の生成に用いる「疾患ラベル」を予めうつ病のレベルに応じた複数の「レベル別疾患ラベル」として、分類器の機械学習を実行してもよい。
４−２．判別方法
実施の形態３は、うつ症状のレベルの判別方法に関する。前記判別方法は、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程と、あらかじめ前記機能的結合ごとにうつ症状のレベルに応じて設定された指標値の基準範囲と比較する工程と、前記被検者が、前記指標値が含まれる基準範囲に対応するうつ症状のレベルであると決定する工程とを備える。
図１６は、生成した分類器を用いて被検者のｒｓ−ｆｃＭＲＩのデータからうつ症状のレベルを判別するためにデータ処理部３２が行う処理を説明するためのフローチャート図である。
ＣＰＵ２０４０は、入力部４０からの処理開始の入力等により、図１５に記載のステッ
プＳ４０１〜Ｓ４０６を行うことにより、指標値を生成する。次にＣＰＵ２０４０は、図１６に示すステップＳ５０１において前記指標値を予めうつ症状のレベルに応じて設定された指標値の基準範囲と比較する。前記１２対の機能的結合は、たとえば、ＢＤＩのうつ症状の重症度と相関するため、うつ症状のレベルを決定するための指標値の基準範囲は、ＢＤＩのうつ症状の重症度との相関から予め決定しておくことができる。ステップＳ５０２においてＣＰＵ２０４０は、図１５のステップＳ４０６で生成された指標値が、どのレベルのうつ症状に対応しているかを決定する（Ｓ５０２）。そして、ＣＰＵ２０４０は、ステップＳ５０３において、前記被検者がステップＳ５０２で決定されたうつ症状のレベルであると判別する。前記ＣＰＵ２０４０が行う各ステップは、指標値を生成する工程を除いた工程の全部または一部をヒトが行ってもよい。
また、実施の形態３の別態様として、上記ステップＳ５０１〜Ｓ５０３を含む処理を実行させて、判別装置としての機能をコンピュータに実行させるプログラムを含む。実施の形態３の別態様として、上記ステップＳ４０１〜Ｓ４０６およびステップＳ５０１〜Ｓ５０３を含む処理を実行させて、判別装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
５．第１の分類器を使用する治療効果の判別装置及び判別方法
５−１．判別装置
判別装置１及び判別装置２は、治療効果を判別するための情報の生成のために使用することができる。すなわち、実施の形態１の別形態である被検者における治療効果の判定するための判別装置は、演算装置と記憶装置とを有する。演算装置は、記憶装置に格納されたプログラムに基づいて、前記第１の分類器により、第１の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第１の値を生成し、前記分類器により、治療開始後であって第１の時点より後の第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素についてうつ症状を評価する第２の値を生成し、前記第１の値と前記第２の値とを比較し、ｉ）前記第２の値が前記第１の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると決定する処理、及び／又は前記第２の値が前記第１の値よりも改善していない場合に、前記治療が被検者のうつ症状の改善に有効でないと決定する処理を実行する。この場合、判別装置１及び判別装置２の構成は、上記３．に記載した構成と同様であるため、説明は省略する。判別処理部３０１０は、分類器データ３１１２により特定される分類器に基づいて、入力データに対して判別処理を行い、前記被検者における治療効果を判別する情報を生成する。
本実施の態様の別の態様は、前記判別装置を、ドラッグリプロファイリングのために使用することを含む。
５−２．判別方法
実施の形態４は、前記第１の分類器に、少なくとも２つの時点において撮像された被検者の安静時ｆＭＲＩのデータから取得される表１に記載の機能的結合についての相関行列の要素を入力し、治療効果を判別する情報を生成する。具体的には、第１の時点の被検者において安静時に計測された、機能的結合の相関行列の要素について、うつ症状を評価するための第１の値を生成する工程と、第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記機能的結合と同一の機能的結合の相関行列の要素ついてうつ症状を評価する第２の値を生成する工程と、第１の値と第２の値とを比較する工程と、第２の値が第１の値よりも改善している場合に、前記治療が被検者のうつ症状の改善に有効であると
決定し、有効である可能性を示す情報を生成する工程とを備える。第１の時点は、治療前であっても、治療開始から所定の時間が経過していてもよいが、好ましくは、治療開始前である。また、一度治療を終了した後、治療が再度開始される場合も治療終了から治療再開前までの間は、「治療前」としてもよい。第２の時点は、治療開始後であって第１の時点よりも後である限り制限されない。
図１７は、第１の分類器を用いて被検者のｆＭＲＩのデータから治療効果を判別するための処理を説明するためのフローチャート図である。
ＣＰＵ２０４０は、入力部４０からの処理開始の入力等により、図１７に記載のステップＳ６０１において、第１の時点においてＭＲＩ撮像部２５によって被検者から撮像された安静時の被検者ｆＭＲＩ測定データ３１１３をインタフェース部４４を介して取得する。
ステップＳ６０２〜Ｓ６０３おいて上記２．に記載の前処理を行い、表１に記載の１２対の機能的結合から選択した全て、あるいは一部の機能的結合について第１の時点での相関行列の要素を抽出する。
ステップＳ６０４において、ＣＰＵ２０４０は、ステップＳ６０３で抽出した第１の時点での相関行列の要素に対してピアソン相関係数を算出する。
ステップＳ６０５〜Ｓ６０６おいて、ＣＰＵ２０４０は、ステップＳ６０４において算出された相関係数を第１の分類器に入力し、前記選択された機能的結合について第１の値を生成する。前記第１の値は、例えば、上記２．で述べた指標値そして算出される。
ＣＰＵ２０４０は、第２の時点において、ステップＳ６０７に示すＭＲＩ撮像部２５によって被検者から撮像された安静時の被検者ｆＭＲＩ測定データ３１１３をインタフェース部４４を介して取得する。
ステップＳ６０８〜Ｓ６０９おいて上記２．に記載の前処理を行い、表１に記載の１２対の機能的結合から選択した全て、あるいは一部の機能的結合について第２の時点での相関行列の要素を抽出する。
ステップＳ６１０において、ＣＰＵ２０４０は、ステップＳ６０９で抽出した第１の時点での相関行列の要素に対してピアソン相関係数を算出する。
ステップＳ６１１〜Ｓ６１２おいて、ＣＰＵ２０４０は、ステップＳ６１０において算出された相関係数を第１の分類器に入力し、前記選択された機能的結合について第１の値を生成する。前記第２の値は、例えば、上記２．で述べた指標値として生成される。
ＣＰＵ２０４０は、ステップＳ６０６において生成された第１の値とステップＳ６１２おいて生成された第２の値を比較する（Ｓ６１３）。
ＣＰＵ２０４０は、ステップＳ６１４において、第１の値よりも第２の値が改善していると判別した場合には（「ＹＥＳ：Ｙ」）、前記被検者のうつ病症状の改善に対する前記治療の状態を示す情報を生成する（Ｓ６１５）。この場合、次にＣＰＵ２０４０は、ステップＳ６１６に進み、前記治療の継続の可能性を示す情報を提示してもよい。医師は、このような治療継続の可能性を示す情報を元に、治療継続をするか否かの判断を行う。特に限定されないが、たとえば、出力画面上には、治療の有効性があると予め設定された領域についての情報が区別可能に表示されており、「治療の状態を示す情報」とは、分類器の
出力する指標値であり、「治療を継続の可能性を示す情報の提示」とは、上記領域を表示した状態で、併せて、この指標値を表示することであってもよい。
また、ＣＰＵ２０４０は、ステップＳ６１４において、第１の値よりも第２の値が改善していないと判別した場合には（「ＮＯ：Ｎ」）、前記被検者のうつ病症状の改善に対する前記治療の状態を示す情報を生成する（Ｓ６１７）。この場合、次にＣＰＵ２０４０は、ステップＳ６１８に進み、前記治療の終了の可能性を示す情報を提示してもよい。
医師は、このような治療終了の可能性を示す情報を元に、治療を終了をするか否かの判断を行う。特に限定されないが、たとえば、出力画面上には、治療の有効性がないと予め設定された領域についての情報が区別可能に表示されており、「治療の状態を示す情報」とは、分類器の出力する指標値であり、「治療の終了の可能性を示す情報の提示」とは、上記領域を表示した状態で、併せて、この指標値を表示することであってもよい。
さらに、ＣＰＵ２０４０は、ステップＳ６１９に進み、他の治療方法に変更可能性を示す情報を提示してもよい。この場合、さらに具体的な治療方法の可能性を提示してもよい。この場合、医師は、表示された治療方法の可能性のうちから、現在の治療方法から変更すべき治療方法を判断することができる。
例えば数６で表される式では、うつ症状の疾患ラベルと健常者のラベルの閾値は０となるように設定されているため、基準値は０となる。
たとえば、「関連重み付け和」をバイオマーカーの指標値として使用する場合は、ステップＳ６１４において、うつ症状が重くなるにつれ、前記指標値は、正方向に大きくなる。したがって、ＣＰＵ２０４０は、前記第１の値よりも前記第２の値が小さくなっている場合に、第１の値よりも第２の値が改善していると判別することができる。また、ＣＰＵ２０４０は、前記第１の値よりも前記第２の値が大きくなっている場合、または第１の値と第２の値に顕著な差が認められない場合に、第１の値よりも第２の値が改善していないと判別することができる。
前記第１の値及び第２の値は、指標値であってもよいが、ステップＳ６０４及びステップＳ６１０で算出される相関係数を、それぞれ第１の値及び第２の値としてもよい。
また、実施の態様３では、ステップＳ６０１は、ステップＳ６０７よりも前に行われている必要があるが、ステップＳ６０２〜Ｓ６０６は、ステップＳ６０７よりも前に行われている必要はない。ステップＳ６０２〜Ｓ６０６は、ステップＳ６０１の後、少なくともステップＳ６１３よりも前のどこかで行われていればよい。
ステップＳ６０１〜Ｓ６１９は、その全部あるいは一部をヒトが行ってもよい。
実施の形態３は、投薬治療の効果判定に使用することも可能である。第１の分類器を使って投薬効果を判定する場合には、少なくとも機能的結合識別番号１及び２の機能的結合を指標とすることが好ましい。
実施の態様４の別の態様は、前記判定方法を、ドラッグリプロファイリングのために使用することを含む。
また、実施の形態４の別態様として、上記ステップＳ６０１〜Ｓ６１９を含む処理を実行させて、判別装置としての機能をコンピュータに実行させるプログラムを含む。このプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の
記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
６．うつ病患者の層別化装置及び層別化方法
６−１．層別化装置
実施の形態１の別形態は、うつ病患者の層別化装置に関する。うつ病が予め設定された複数のサブクラスに層別化されている場合に、判別装置１及び判別装置２は、それぞれうつ病患者を医師が層別化するための支援をする層別化装置１及び層別化装置２（以下、総称するときは、「層別化装置」という）として使用することができる。層別化装置は、演算装置と記憶装置とを含む。すなわち、うつ症状が予め設定された複数のサブクラスに層別化されている場合に、うつ病患者の層別化するための判別装置の演算装置は、記憶装置に格納されたプログラムに基づいて、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成し、、前記指標値を、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較し、、前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する、という処理を実行する。この場合、層別化装置１及び層別化装置２の構成は、上記３．に記載した判別装置１及び判別装置２の構成と同様であるため、説明は省略する。図１９と図２０には、層別化装置１及び層別化装置２の機能ブロック図を示す。図１９と図２０では、判別処理部３０１０が、層別化処理部３０２０となる点を除き、図５及び図１５と同じである。層別化処理部３０２０は、分類器データ３１１２により特定される分類器に基づいて、入力データに対して層別化処理を行い、うつ病患者の層別化のための判別を行い結果を出力する。
医師は、層別化装置１または層別化装置２にから出力される、うつ病患者の層別化のための判別結果を支援情報として、患者の層別化を行う。この意味で、このような装置を、層別化支援装置と呼ぶこともできる。
６−２．層別化方法
実施の形態５は、うつ病患者得層別化するための層別化する方法に関する。前記判別方法は、うつ病が予め設定された複数のサブクラスに層別化されている場合に、安静時の被検者において計測された機能的結合の相関行列の要素について、うつ症状を評価するための指標値を生成する工程と、あらかじめ前記機能的結合ごとに前記サブクラスに応じて設定された指標値の基準範囲と比較する工程と、前記被検者が、前記指標値が含まれる基準範囲に対応する前記サブクラスであると決定する工程を備える。
図１８は、生成した分類器を用いて被検者のｆＭＲＩのデータからうつ症状のサブクラスを判別するための情報を生成するデータ処理部３２が行う処理を説明するためのフローチャート図である。
ＣＰＵ２０４０は、入力部４０からの処理開始の入力等により、図１４に記載のステップＳ４０１〜Ｓ４０６を行うことにより、指標値を生成する。次にＣＰＵ２０４０は、図１８に示すステップＳ７０１において前記指標値を、予めうつ病のサブクラスに応じて設定された指標値の基準範囲と比較する。うつ病のサブクラスは、臨床所見からメランコリー型ＭＤＤ、非メランコリー型ＭＤＤ及び治療抵抗性ＭＤＤのサブクラスに分類される。このため、たとえば、メランコリー型ＭＤＤを「疾患ラベル」として生成した分類器を用いることで、脳活動を計測した対象者からの計測データが、メランコリー型ＭＤＤの疾患ラベルに対応するか否かを判別できる。同様にして、非メランコリー型ＭＤＤまたは治療抵抗性ＭＤＤを「疾患ラベル」として生成した分類器を用いることで、脳活動を計測した対象者からの計測データが、非メランコリー型ＭＤＤまたは治療抵抗性ＭＤＤの疾患ラベ
ルに対応するか否かを判別できる。
ステップＳ７０２においてＣＰＵ２０４０は、図１５のステップＳ４０６で生成された指標値が、どのサブクラスに対応しているかを決定する。そして、ＣＰＵ２０４０は、ステップＳ７０３において、前記被検者がステップＳ７０２で決定されたうつ病のサブクラスに属すると判別する。医師は、このようにして出力されるサブクラスの情報を、うつ病患者の層別化のための支援情報として、患者の層別化を行う。この意味で、このような方法を、層別化支援方法と呼ぶこともできる。前記ＣＰＵ２０４０が行う各ステップは、指標値を生成する工程を除いた全部または一部をヒトが行ってもよい。
また、実施の形態５の別態様として、上記ステップＳ７０１〜Ｓ７０３を含む処理を実行させて、層別化装置としての機能をコンピュータに実行させるプログラムを含む。実施の形態５の別態様として、上記ステップＳ４０１〜Ｓ４０６およびステップＳ７０１〜Ｓ７０３を含む処理を実行させて、層別化装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
７．第２の分類器の生成と第２の分類器を使用する治療効果の判別
７−１．治療効果判別装置
実施の形態６は、複数の対象者について、治療のために投与した治療効果がみられる群とみられない群とを判別するための第２の分類器を生成する分類器生成処理を実行し、医師が治療効果の判断をする支援のための指標となる情報を生成するする、治療効果判別装置、あるいは、治療支援装置に関する。
治療効果判別装置の構成は、基本的に図１と同じであるが、図１におけるデータ処理部３２を実施の形態６ではデータ処理部６２と読み替える。また、治療効果判別装置において、データ処理部６２の構成は、図１及び図２と基本的に同様であるが、実施の形態６では、図１において記憶部３６を記憶部６６と、制御部４２を制御部６２と、入力部４０を入力部７０と、インタフェース部４４をインタフェース部７４と、データ収集部４６をデータ収集部７６と、画像処理部４８を画像処理部７８と、表示部３８を表示部６８と、表示制御部３４を表示制御部６４と読み替える。また、図２において、データ処理部３２のコンピュータ本体２０１０をコンピュータ本体６０１０と、メモリドライブ２０２０をメモリドライブ６０２０、ディスクドライブ２０３０をディスクドライブ６０３０と、演算装置（ＣＰＵ）２０４０を演算装置（ＣＰＵ）６０４０と、ディスクドライブ２０３０をディスクドライブ６０３０と、メモリドライブ２０２０をメモリドライブ６０２０と、バス２０５０をバス６０５０と、ＲＯＭ２０６０をＲＯＭ６０６０と、ＲＡＭ２０７０をＲＡＭ６０７０と、不揮発性記憶装置２０８０を不揮発性記憶装置６０８０と、通信インタフェース２０９０を通信インタフェース６０９０と読み替える。データ処理部６２のハードウェアとしては、上述のとおり、特に限定されないが、汎用コンピュータを使用することが可能である。
（分類器の生成）
図２３は、第２の分類器の生成処理および生成された第２の分類器による判別処理を行うための機能ブロック図である。
不揮発性記憶装置６０８０には、健常者、治療を受けたうつ病患者を含む複数の対象者の脳内の複数の所定領域における脳活動を示す信号をＭＲＩ装置により時系列で予め測定した信号の情報であって、第１の時点で測定した対象者のｒｓ−ｆｃＭＲＩ測定データ６
１０２と、第２の時点で測定した対象者のｒｓ−ｆｃＭＲＩ測定データ６１０４と、前記ｒｓ−ｆｃＭＲＩ測定データを計測した各対象者に関連づけられた、複数の治療データ６１０５とが格納されている。
ここで、治療データ６１０５には、対象者の疾患ラベル、治療履歴（治療方法、投薬した薬剤名、用量、投与期間等）の情報が含まれる。ＣＰＵ６０４０は、第１の時点の対象者のｒｓ−ｆｃＭＲＩ測定データ６１０２と、第２の時点の対象者のｒｓ−ｆｃＭＲＩ測定データ６１０４と、対象者の治療データ６１０５に基づいて、複数の対象者において治療効果がみられた群（「寛解群(remitted)」）と治療効果がみられない群（「非寛解群(non-remitted)」）とを識別するための分類器を生成する処理を行う。したがって、疾患ラベルには、この場合、「寛解ラベル」と「非寛解ラベル」とが含まれる。
第１の相関計測部６００２は、第１の時点の対象者のｒｓ−ｆｃＭＲＩ測定データ６１０２から、複数の所定領域間の脳活動の機能的結合の相関行列を対象者ごとに算出する。算出された機能的結合の相関行列のデータは、対象者ごとに、不揮発性記憶装置６０８０に、機能的結合の相関行列の第１の相関計測データ６１０６として格納される。
第２の相関計測部６００４は、第２の時点の対象者のｒｓ−ｆｃＭＲＩ測定データ６１０４から、複数の所定領域間の脳活動の機能的結合の相関行列を対象者ごとに算出する。算出された機能的結合の相関行列のデータは、対象者ごとに、不揮発性記憶装置６０８０に、機能的結合の相関行列の第２の相関計測データ６１０８として格納される。
第２の分類器生成部６００８は、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差に基づいて、回帰式を求める。前記回帰式は、第２の分類器を特定するための第２の分類器データ６０２０として、不揮発性記憶装置６０８０に格納する。治療効果判別処理部６１００は、第２の分類器データ６０２０により特定される第２の分類器に基づいて、入力データに対して治療効果の判別処理を行う。
図２１は、第２の分類器を生成するためにデータ処理部６２が行う処理を説明するためのフローチャートを示す。
以下、図２３で説明した処理を図２１を参照しながら、より詳しく説明する。
図２１において、データ処理部６２は、入力部７０からの処理開始の入力等により、分類器生成の処理が開始されると（スタート）、第１の時点の各対象者のｒｓ−ｆｃＭＲＩ測定データ６１０２を読み出し、各表１に示される複数の機能的結合それぞれについて、第１の時点の相関を計測する（Ｓ８０２）。この相関の計測は、各機能的結合のピアソン相関係数を算出することによって行われる。
次にデータ処理部６２は、記憶部６６から第２の時点の各対象者のｒｓ−ｆｃＭＲＩ測定データ６１０４を読み出し、各対象者について表１に示される複数の機能的結合それぞれについて、第２の時点の相関を計測する（Ｓ８０３）。この相関の計測は、各機能的結合のピアソン相関係数を算出することによって行われる。
次に、データ処理部６２は、各対象者の表１に示される複数の機能的結合のうち、たとえば、第１の機能的結合および第２の機能的結合のそれぞれについて、第１の時点の相関と第２の時点の相関との差を算出する（Ｓ８０４）。前記差は、例えば、下式（４）に基づいて算出される。
式（４）において、ＦＣ１は、表１における第１番目の結合強度を示し、ＦＣ２は、表１における第２番目の結合強度を示す。すなわち、関連重み付け和の計算の結果、重みの絶対値が最も大きな結合がＦＣ１であり、２番目に大きな結合がＦＣ２である。上付きの“post”は、第２の時点を示し、“pre”は、第１の時点を示す。典型的には、第１の時
点は、治療前（投薬前）の時点であり、第２の時点は、治療開始後（投薬開始後）、所定の時点である。また、sign(w1)は、ＦＣ１の重みの符号を示し、sign(w2)は、ＦＣ２の重みの符号を示す。たとえば、後述するように図３３ｆにおいて作成したこの分類器は、正確性：０．７５， ＡＵＣ：０．７９， 特異度０．８８、感度０．４３を示す。
さらに、データ処理部６２は、各対象者の各機能的結合における前記２つの相関の差と、治療データ６１０５に含まれる説明変数（「寛解群」または「非寛解群」のラベル）との関係を説明するための回帰式を求める。回帰式は、特に制限されないが、線形回帰分析によって求めることができる。この回帰式に基づいて、相関状態空間における治療効果あり群と、治療効果なし群を識別する分類器を生成することができる（Ｓ８０４）。
ここで、「相関状態空間」とは、このようにして複数の機能的結合の第１および第２の時点の間の強度差（必要により、さらに強度差に重みの符号をかけたもの）をそれぞれ空間の軸として張られる空間のことをいう。上記の例では、２次元の空間となる。
上述のとおり、例えば、各対象者の表１に示す機能的結合識別番号１における第１の時点と第２の時点の前記相関の差（Δsign(W)ＦＣ１とする）と、表１に示す機能的結合識
別番号２における第１の時点と第２の時点の前記相関の差（Δsign(W)ＦＣ２とする）の
相関状態をプロットすると図３３ｆのようになる。Δsign(W)ＦＣ１とΔsign(W)ＦＣ２の相関と、各対象者における治療効果の有無との関係をみると負の傾きをもつ線形回帰式を求めることができる。そして、この線形回帰式に基づく回帰直線を基準として、回帰直線よりも上に分布する対象者を治療効果が認められない群、回帰直線よりも下に分布する対象者を治療効果が認められる群として識別することができる。
医師は、このようにして出力される治療効果に関する情報を、支援情報として、うつ病患者に対する治療効果の判断を行う。
（説明変数としてのΔsign(W)ＦＣ１とΔsign(W)ＦＣ２）
上記の例では、各結合強度の重みの符号sign(W)は、図として表示する際の便宜のため
であり、本質的には、治療（投薬）前と治療（投薬開始後）後におけるＦＣ１とＦＣ２の強度変化が重要である。
図３３ｆのように、寛解と非寛解を区別せずに見れば、第１の機能的結合の強度変化Δsign(W)ＦＣ１と第２の機能的結合の強度変化Δsign(W)ＦＣ２との間には、治療開始前と治療開始後所定期間経過後の強度変化の間には相関が見られず、各々の変化に対する補償のような現象は見られない。言い換えれば、これらの結合強度の変化は、相互に独立であることになり、治療の効果を説明するための変数の選択としては適切であることを示している。
したがって、上記のように、治療開始の前後における強度変化が相互に独立と見なせる
機能的結合を用いて、上記のような第１および第２の機能的結合に限られず、他の機能的結合を用いることも可能である。また、治療効果を説明するための変数としては、より多くの機能的結合を、表１に示した機能的結合から選択してもよい。
（判別処理）
図２２は、実施の形態７生成した第２の分類器を用いて被検者のｒｓ−ｆｃＭＲＩのデータから治療効果を判別するデータ処理部６２が行う処理を説明するためのフローチャート図である。
ＣＰＵ６０４０は、入力部７０からの処理開始の入力等により、ステップＳ９０１において、ＭＲＩ撮像部２５から第１の時点における安静時の被検者ｒｓ−ｆｃＭＲＩ測定データ６１１３をインタフェース部６４を介して取得する。
ＣＰＵ６０４０は、ステップＳ９０２〜Ｓ９０３おいて上記２．に記載の前処理を行い、表１に記載の機能的結合の全て、あるいは一部について相関行列の要素を抽出し、第１の時点の各機能的結合の相関を計測する。
ＣＰＵ６０４０は、ステップＳ９０４において、ＭＲＩ撮像部２５から第２の時点における安静時の被検者ｒｓ−ｆｃＭＲＩ測定データ６１１３をインタフェース部６４を介して取得する。
ＣＰＵ６０４０は、ステップＳ９０５〜Ｓ９０６おいて上記２．に記載の前処理を行い、表１に記載の機能的結合の全て、あるいは一部について相関行列の要素を抽出し、第２の時点の各機能的結合の相関を計測する。
次に、ＣＰＵ６０４０は、ステップＳ９０７おいて、前記複数の機能的結合それぞれについて、ステップＳ９０３で計測された第１の時点の各機能的結合の相関とステップＳ９０６で計測された第２の時点の各機能的結合の相関の差を算出する。
ＣＰＵ６０４０は、ステップＳ９０７で算出された差を第２の分類器に入力し（Ｓ９０８）、第２の分類器によって前記被検者において治療効果が認められているか、認められていないかを判別する（Ｓ９０８）。
各相関の計測方法、各相関の差の算出方法、及び治療効果の判別方法は、第２の分類器の生成方法で記載した内容と同様である。
ここで、第１の時点は、治療開始前であっても、治療開始から所定の時間が経過していてもよいが、好ましくは、治療開始前である。また、一度治療を終了した後、治療が再度開始される場合も治療終了から治療再開前までの間は、「治療開始前」としてもよい。第２の時点は、治療開始後であって第１の時点よりも後である限り制限されない。
また、実施の形態７は、ステップＳ９０１は、ステップＳ６０６よりも前に行われている必要があるが、ステップＳ９０２〜Ｓ９０３は、ステップＳ９０４よりも前に行われている必要はない。ステップＳ９０２〜Ｓ９０３は、ステップＳ９０１の後、少なくともステップＳ９０７よりも前のどこかで行われていればよい。
また、実施の形態６の別形態は、前記第２の分類器、治療効果判別装置、あるいは、治療支援装置を、ドラッグリプロファイリングのために使用することを含む。
７−２．治療効果判別方法
実施の形態７は、第２の分類器を用いた治療効果の判別方法に関する。前記判別方法は、第１の時点の被検者の安静時の脳内における、前記複数の機能的結合の第１の相関を計測する工程と、第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第２相関を計測する工程と、前記被検者の前記複数の機能的結合の前記第２の相関と前記第１の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する工程とを備える。具体的な工程は、上記Ｓ９０１〜Ｓ９０９に準ずるが、前記工程は、全部または一部をヒトが行ってもよい。
ここでも、このようにして出力される治療効果に関する判別情報を、支援情報として、うつ病患者に対する治療効果の判断を行う。
実施の形態７の別の形態は、前記治療効果判別方法を、ドラッグリプロファイリングのために使用することを含む。
また、実施の形態７の別形態として、上記ステップＳ８０１〜Ｓ８０４を含む処理を実行させて、第２の分類器としての機能をコンピュータに実行させるプログラムを含む。実施の形態７の別態様として、上記ステップＳ８０１〜Ｓ８０４およびステップＳ９０１〜Ｓ９０９を含む処理を実行させて、治療効果判別装置としての機能をコンピュータに実行させるプログラムを含む。実施の形態７の別態様として、上記ステップＳ９０１〜Ｓ９０９を含む処理を実行させて、治療効果判別装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
８．層別化されたうつ病患者における第２の分類器を使用した治療効果の判定
８−１．治療効果判別装置
実施の形態８は、実施の態様４に記載のうつ病患者の層別化装置において、層別化されたうつ病患者からのｒｓ−ｆｃＭＲＩデータを使って、実施の形態６の治療効果判別装置において、治療効果を判別する装置である。
すなわち、実施の態様８は、上記３．で述べた第１の分類器生成処理と、上記６．で述べた層別化処理と上記７．で述べた第２の分類器生成処理と、第１の相関計測処理と、第２の相関計測処理と、前記第２の分類器により、前記被検者に対する治療効果を判別する判別処理とを実行する治療効果判別装置である。
したがって、上記３．の記載、上記６−１．の記載と、上記７−１．記載は、ここに援用される。
実施の形態８の別の形態は、前記治療効果判別装置を、ドラッグリプロファイリングのために使用することを含む。
８−２．治療効果判別方法
実施の形態９は、層別化されたうつ病患者における第２の分類器を使用した治療効果の判別方法に関する。
図２４に、層別化されたうつ病患者における第２の分類器を使用した治療効果を判別するためのフローチャートを示す。
図２４において、ステップＳ７０１〜Ｓ７０３は、図１８と同じである。したがって、
上記６−２．の説明は、ここに援用される。
ＣＰＵ６０４０は、図２４に示すステップＳ７０３においてサブクラスが決定された被検者について、ステップＳ９２１〜Ｓ９２９を行う。ステップＳ９２１〜Ｓ９２９の各ステップは、上記７−２．及び図２２に記載のステップＳ９０１〜Ｓ９０９の各ステップに対応する。したがって、上記７−２．の記載は、ここに援用される。
なお、特に限定されないが、以上の説明において、治療効果判別装置の支援により、医師が、現在継続中の治療方法（特定薬の投薬など）に効果が見られないと判断した場合は、それ以後に、医師は、別の治療方法（別の特定薬の投与、ニューロフィードバック法、修正型電気けいれん療法、反復経頭蓋磁気刺激治療法など）の治療方法を併用する、あるいは、別の治療方法へと治療方法を変更する、という判断をすることができる。
実施の形態９の別の形態は、前記治療効果判別方法を、ドラッグリプロファイリングのために使用することを含む。
また、実施の形態８の別形態として、上記ステップＳ９２１〜Ｓ９２９を含む処理を実行させて、第２の分類器としての機能をコンピュータに実行させるプログラムを含む。実施の形態８の別形態として、上記ステップＳ７０１〜Ｓ７０３およびステップＳ９２１〜Ｓ９２９を含む処理を実行させて、治療効果判別装置としての機能をコンピュータに実行させるプログラムを含む。実施の形態７の別形態として、上記ステップＳ４０１〜Ｓ４０６、ステップＳ７０１〜Ｓ７０３およびステップＳ９０１〜Ｓ９０９を含む処理を実行させて、治療効果判別装置としての機能をコンピュータに実行させるプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
９．ニューロフィードバック
９−１．脳活動訓練装置
実施の形態１０は、被訓練者の脳活動に対してフィードバックを与え、健常者の相関（結合）状態に近づける訓練（ニューロフィードバック訓練、あるいは単に訓練ともいう）ための脳活動訓練装置に関する。
なお、より一般的には、このような脳活動訓練装置は、健常群と患者群の脳活動という関係において、被訓練者の脳活動の機能的結合状態を健常者の脳活動の機能的結合状態に近づける訓練に留まらず、目標とする脳活動の機能的結合状態に、被訓練者の現在の脳活動の機能的結合状態を近づけるという意味での訓練に使用することも可能である。健常者の脳活動の機能的結合状態及び目標とする脳活動の結合状態を、後述するように「目標パターン」と呼ぶ。
ここで、被訓練者は、上記３．記載のうつ症状の判別装置１、うつ症状の判別装置２又はうつ症状の判別方法により、「うつ症状」のラベルを付された被検体（好ましくはＭＤＤ、より好ましくはメランコリー型ＭＤＤの被検体）、上記４．に記載のうつ症状のレベルの判別装置及び判別方法により、うつ症状のレベルを判定された被検体（好ましくはＭＤＤ、より好ましくはメランコリー型ＭＤＤの被検体）、上記６．に記載のうつ病患者の層別化装置及び層別化方法により、層別化された被検体（好ましくはＭＤＤに層別化された被検体、より好ましくはメランコリー型ＭＤＤに層別化された被検体）であることが好ましい。従って、上記記３．記載のうつ症状の判別装置１、うつ症状の判別装置２又はうつ症状の判別方法に関する説明、上記４．に記載のうつ症状のレベルの判別装置及び判別方法に関する説明及び上記６．に記載のうつ病患者の層別化装置及び層別化方法に関する
説明は、ここに援用される。
図３４は、このような脳活動訓練装置の構成の概念を示す図である。
脳活動訓練装置のハードウェア構成としては、たとえば、上述した図１に示されるＭＲＩ装置１０の構成を用いることができる。脳活動訓練装置のハードウェア構成は、基本的に図１と同じであるが、図１におけるデータ処理部３２を実施の態様１０ではデータ処理部１０２と読み替える。また、脳活動訓練装置において、データ処理部１０２の構成は、図１及び図２と基本的に同様であるが、実施の形態１０では、図１において記憶部３６を記憶部１０６と、制御部４２を制御部１１２と、入力部４０を入力部１１０と、インタフェース部４４をインタフェース部１１４と、データ収集部４６をデータ収集部１１６と、画像処理部４８を画像処理部１１８と、表示部３８を表示部１０８と、表示制御部３４を表示制御部１０４と読み替える。また、図２において、データ処理部３２のコンピュータ本体２０１０をコンピュータ本体１００１０と、メモリドライブ２０２０をメモリドライブ１００２０、ディスクドライブ２０３０をディスクドライブ１００３０と、演算装置（ＣＰＵ）２０４０を演算装置（ＣＰＵ）１００４０と、ディスクドライブ２０３０をディスクドライブ６０３０と、メモリドライブ２０２０をメモリドライブ６０２０と、バス２０５０をバス１００５０と、ＲＯＭ２０６０をＲＯＭ１００６０と、ＲＡＭ２０７０をＲＡＭ１００７０と、不揮発性記憶装置２０８０を不揮発性記憶装置１００８０と、通信インタフェース２０９０を通信インタフェース１００９０と読み替える。データ処理部６２のハードウェアとしては、上述のとおり、特に限定されないが、汎用コンピュータを使用することが可能である。
以下に説明する脳活動訓練装置の駆動方法では、脳機能画像法により脳活動を示す信号を時系列に計測するための脳活動検出装置としては、リアルタイムｆＭＲＩを用いるものとして説明する。
図３４を参照して、まず、ＭＲＩ装置１０により被訓練者の脳活動が所定時間の間について時系列に測定される。ｆＭＲＩとして撮像するためにＥＰＩ（Echo-planar imaging
）撮像が実行される。
続いて、データ処理部１０２により、撮像された画像がリアルタイムで再構成処理される。
ここで、実施の形態１で説明したとおり、表１に示す機能的結合識別番号（表１でＩＤ）１〜１２の関心領域間の機能的結合が、「うつ症状」のラベルの判別のために特徴選択により選択された。
ここで、被訓練者においては、上記１２個の機能的結合のうち、少なくとも１つの特定の機能結合を「訓練対象の機能的結合」として、さらに選択・抽出する。ここで、特に限定されないが、このような「訓練対象の機能的結合」としては、たとえば、表１で示した機能結合のうち、うつ病ラベルの判別に最も寄与度の大きい「左背外側前頭前皮質と、左楔前部および左後部帯状皮質との間の機能的結合」（表１のＩＤ＝１の結合）を選択するものとする。ここで、「寄与度」とは、「関連重み付け和における重みの絶対値」を意味するものとする。
図３２において、説明したように、この結合は、上述したような「抗うつ薬」の投与によって、健常者側とは反対方向に変化する結合のうち、最大の寄与度を有する結合である。
なお、「訓練対象の機能的結合」としては、このＩＤ＝１の結合のみを選択してもよいし、ＩＤ＝１の結合に加えて、他の結合を対象として選択してもよい。このとき、ＩＤ＝１の結合に加えて選択する機能的結合としては、２番目に寄与度の大きなＩＤ＝２の結合というような選択をしてもよい。あるいは、ＩＤ＝１の結合に加えて、「抗うつ薬」の投与によって、健常者側とは反対方向に変化する他の結合、たとえば、ＩＤ＝３の結合を選択してもよい。ＩＤ＝１の結合に加えて選択する他の機能結合の個数は、１つでもそれ以上でもよい。
さらに、データ処理部１０２は、「訓練対象の機能的結合」に対応するの関心領域について、機能的結合の所定時間内の時間的な相関を算出する手続きを行う。
すなわち、データ処理部１０２は、リアルタイムで測定した安静時のｆＭＲＩの連続するｎ個（ｎ：自然数、ｎ≧１）の計測時刻分のｆＭＲＩデータ（信号）から、まず、各関心領域の「活動度」を算出する。この活動度は、たとえば、連続するｎ個の計測時刻の計測値についての平均値とすることができる。ここで、連続するｎ個の計測時刻からなる期間を「サンプルステップ」と呼ぶことにする。
そして、このような各サンプルステップの活動度について、データ処理部１０２は、所定の期間の計測ウィンドウについて、以下の数式により、活動度の時間軸上の相関（以下、「時間相関」という）を算出する。なお、計測ウィンドウは、時間軸上で連続する複数のサンプルステップの活動度（たとえば、ｍ個（ｍ：自然数、ｍ≧２））を含むような期間である。
ここで、ｘは、「訓練対象の機能的結合」のうちの一方の関心領域の「活動度」を意味し、ｙは、「訓練対象の機能的結合」のうちの他方の関心領域の「活動度」を意味する。また、ｘ（バー）（以下、文字ｘの頭部に−を付して「平均」を表す記号を、「ｘ（バー）」と呼ぶ。他の文字についても同様である。）とは、上記のような計測ウィンドウ内での活動度ｘの平均を表す。「ｙ（バー）」についても同様である。上記の式において、Σは、たとえば、計測ウィンドウ内のｍ個のサンプルステップの活動度について和を取ることを意味する。また、計測ウィンドウは、特に限定されないが、たとえば、１０秒〜２０秒程度に設定することができる。
ここで、機能結合の時間相関は、符号として、正または負の値を取りうる。
そして、機能結合の時間相関の符号は、たとえば、「健常者」であれば、負であるのに対して、「うつ症状」有りと判別された被験者では、正である、というような状況がありうる。もちろん、その逆の場合もありうるし、あるいは、両者で同一の符号であることもありうる。機能結合の時間相関の符号およびその絶対値の大きさを併せて、「時間相関のパターン」と呼ぶことにする。
そこで、たとえば、「健常者」において、時間相関の符号が負であるような機能的結合を「訓練対象の機能的結合」として選択した場合、被訓練者の訓練対象の機能的結合の時間相関の健常者の時間相関に対する「近似度」を以下の関数Ｆ１のようにして算出できる。
ｉ）被訓練者の訓練対象の機能的結合の時間相関の符号が正であれば、一定値Ｖ１をとり、前記符号が負であれば、値がＶ２（＞Ｖ１≧０）となり、時間相関の絶対値の増加にとともに、値Ｖ２が増加するような関数Ｆ１。
逆に、「健常者」において、時間相関の符号が正であるような機能的結合を「訓練対象の機能的結合」として選択した場合、被訓練者の訓練対象の機能的結合の時間相関の健常者の時間相関に対する「近似度」を以下の関数Ｆ２のようにして算出できる。
ｉｉ）被訓練者の訓練対象の機能的結合の時間相関の符号が負であれば、一定値Ｖ１をとり、前記符号が正であれば、値がＶ２（＞Ｖ１≧０）となり、時間相関の絶対値の増加にとともに、値Ｖ２が増加するような関数Ｆ２。
なお、「近似度」としては、上記のような例に限られず、符号を含めて、健常者の機能結合の時間相関に近似するパターンであるほど、「近似度」が高くなる関数であれば、他の関数でもよい。
また、特に限定されないが、「訓練対象の機能的結合」が２つ以上選択されている場合は、各々の近似度を平均した値を「近似度」としてもよいし、あるいは、近似度を各結合の寄与度で重み付け平均した値を「近似度」としてもよい。
さらに、データ処理部１０２は、上記のようにして算出した近似度に応じて、報酬値を算出する。ここで、報酬値（スコアＳＣ）は、目標の脳活動の結合強度に対する近似度に応じて、より近似度が高いほど、高い得点となる。
データ処理部１０２は、ディスプレイ６に算出されたスコアＳＣを表示することで、被訓練者２にフィードバックする。なお、フィードバックされる情報としては、スコアの数値そのものでもよいし、あるいは、例えば図３５に示すような、前記数値の大きさに応じて、大きさの変化する図形でもよい。その他、被訓練者が、スコアの大きさを認識可能な画像であれば、他の構成を用いてもよい。
以後は、所定の期間について、ＥＰＩ撮像から、スコアのフィードバックの処理をリアルタイムで繰り返す。
図３６は、ニューロフィードバック訓練シーケンスの一例を示す図である。
図３６に示すように、まず、第１日において、被訓練者は、事前の安静時の脳活動測定される。
その後、第１日において、ニューロフィードバックによる訓練を行う。
以後、第２日、第３日についても、ニューロフィードバック訓練を行う。
第４日（最終日）には、ニューロフィードバックによる訓練を行った後に、事後の安静時の脳活動の計測を実施する。訓練にかける日数は、これより少なくとも、あるいは、多くともよい。
さらに、被訓練者に対して、所定期間の経過後、たとえば、２か月後において、安静時
の脳活動の計測を実施する。
ニューロフィードバック訓練の治療効果は、上記５．に記載の第１の分類器を使用する治療効果の判別装置及び判別方法、又は上記７．に記載の第２の分類器の生成と第２の分類器を使用する治療効果の判別方法によって、判定することができる。
従って、上記５．に記載の第１の分類器を使用する治療効果の判別装置及び判別方法、又は上記７．に記載の第２の分類器の生成と第２の分類器を使用する治療効果の判別方法の説明は、ここに援用される。
９−２．脳活動訓練装置を制御するコンピュータプログラム
実施の形態１０の別形態として、上記処理を実行させて、脳活動訓練装置を駆動するコンピュータプログラムを含む。これらのプログラムは、ハードディスク、フラッシュメモリ等の半導体メモリ素子、光ディスク等の記憶媒体に記憶されていてもよい。前記記憶媒体へのプログラムの記憶形式は、前記提示装置が前記プログラムを読み取り可能である限り制限されない。前記記憶媒体への記憶は、不揮発性であることが好ましい。
（実施例）
I．データ収集及び評価方法
１．データ収集対象
広島大学病院及び広島市内のクリニックで募集した105名の患者について、精神障害の
診断と統計マニュアル-IVの基準の大うつ病性障害（MDD）の患者を選別するために、M.I.N.I.（精神疾患簡易構造化面接法）を使用してスクリーニングした。現在及び過去に躁病症状、精神病症状、アルコール依存症又は乱用、薬物（substance）依存症又は乱用、反
社会性人格障害を除外基準とした。そして、最終的には、自己申告でうつ症状を伴う93名の患者をMDD分類器の訓練データセットとして選択した。患者は、投薬を開始する前に、
又は投薬開始から０〜２週間以内に後述する条件で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未満である）数を等しくした（表２a）。知能指標値（IQ）を評価するためのJapanese Adult Reading Test（JART）の得点については、訓練データセットに3つの欠失
データ（メランコリー型MDD群に1つ、全MDD群及びメタンコリー型MDD群の両方に対する健常対照群うち２つ）と、テストデータセットの治療抵抗型MDD群に2つの欠失データがあった。BDIスコアでは、訓練データセットに、全MDD群のみに対する健常対照群に2つの欠失
データと抗うつ剤治療後のデータに1つの欠失データがあった。
表２bには、本実施例で使用するテストデータの内訳を示す。
２．独立外部バリデーションコホートに対する汎化
日本放射線医学研究所において独立した検証コホート（千葉のコホート）を作成した。参加者について、M.I.N.I.に基づいて精神疾患の生涯履歴を評価した。
MDD患者には併存症の精神障害はなく、健常対照群には身体的、神経学的、精神医学的
障害はなく、現在または過去の薬物乱用の病歴はなかった。すべての参加者に対して実験前に書面によるインフォームドコンセントを実施した。今回の実験は、1964年ヘルシンキ宣言及びその後の改正に規定された倫理基準に従って、放射性医薬品安全委員会および放
射線医学研究所の機関レビュー委員会によって承認された。
３．うつ病及びうつ症状の評価
うつ病及びうつ病の症状を評価するためにBeck Depression Inventory（BDI）を使用した。さらに、治療効果の評価には、Hamilton Rating Scale for Depression（HAMD）も使用した。
４．非メランコリー型MDD及び治療抵抗型MDD
非メランコリー型MDD群は、BDIスコア17以上のすべてのMDDを含む。治療抵抗型MDDは、MDDと診断された中で、抗うつ薬を2種類以上投与しても抑うつ症状が改善しない者である。
５．他の精神疾患における分類器の評価
自閉症（ASD）患者および統合失調症（SSD）患者のfMRIデータは、Yahataらの文献（NATURE COMMUNICATIONS | 7:11254 | DOI: 10.1038/ncomms11254）に記載のデータを使用した。うつ病の合併症による影響を最小限に抑えるために、ASD群は、抗うつ薬が有効でな
かった者に限定した。fMRIデータは、MDD患者健常者と同様に、目を開いた安静状態で取
得した。
６．fMRIデータの取得
fMRIデータの取得の際、被験者に、暗い照明のあるスキャンルームでは、特に何も考えず、眠らずに、モニター画面中央の十字マークを見続けてもらうように依頼した。各施設におけるfMRIデータ取得条件の詳細は表３に示す。
８．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´n〜ez 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全幅半値ガウスカーネルで平滑化した。これらの前処理ス
テップの後、時系列データのフレーム間の相対的な変化に基づいて、スクラビング手順を実行して余分な頭部運動を伴う任意のボリューム（すなわち、機能的画像）を除外した。（頭部運動の要約については表４参照）。
各参加者について、Brainvisa Sulci Atlas（BSA; http：// brainvisa）に解剖学的に定義された大脳皮質全体を覆う137の関心領域（ROI）のそれぞれについて、fMRIデータを経時的に抽出した。本実施例において、施設１において、参加者の小脳の構造及び機能画像が含まれていなかったため、ROIには小脳は含まれていない。バンドパスフィルタ（0.008〜0.1Hz）をかけた後、以下の9つのパラメータ（再配列からの6つの頭部運動パラメー
タ、白質の時間的ゆらぎ、脳脊髄液の時間的ゆらぎ、脳全体のゆらぎ）を直線的に回帰した。
各参加者の9,316機能的結合（FC）それぞれの行列を得るため、137のROI間のペアワイ
ズピアソン相関を算出した。
II.実施例１．メランコリー型MDD分類のための12対の機能的結合の選択
メランコリー型MDDを分類するための12対の機能的結合を選択するため、表１aに示す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対の機能的結合であった（図２６）。さらに、LOOCVを行い、最終的に12対の機
能的結合を選択した（表５及び図２６）。
選択された１２個の機能的結合ＦＣは、はたして、外ループの手続きの全体にわたって、大きな重みで、頻繁に選択されたかものかどうかについて検討した。これは、最後に選択された１２個の機能的結合ＦＣの安定性およびロバスト性に関して重要である。
このような検討のために、以下の式
でｋ番目のＦＣ(ｋ＝１，２，…，９３１６)のための累積的な絶対的な重みを定義する。
ここで、ＮがＬＯＯＣＶ回数(つまり被験者の数)であり、ｗi kは、i番目のＬＯＯＣＶ
のｋ番目の機能的結合ＦＣに関連した重みである。
累積的な重みｃi kがより大きな値であることは、ＬＯＯＣＶの全体にわたって、MDDとHCへの分類に対してｋ番目の機能的結合ＦＣが、より重要な貢献をすることを意味する。
表２で、Lat.の“Ｌ”と“Ｒ”は、左脳と右脳の区別を示す。BSAは、ブロードマン領
野を示し、BAは、ブロードマンの領野の番号を示す。rControlは、健常者対照群に対する相関係数を示す。rMDDは、メランコリー型MDD群に対する相関係数を示す。Weightは、関
連重み付け和の重みを示す。
図２５に選択された12対の機能的結合とその重みを示す。
特に重みが大きかった２つの機能的結合は、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の分布を図２７aに示す。黒棒がメランコリー
型MDD群であり、白棒が健常対照群である。正確性は70 %であった（感度64 %, 特異度77 % AUC 0.77; 並べ替え検定においてp=.049：図３０a）。この結果は、選択された12対の
機能的結合から構成される分類器（以下、「本発明の第１の分類器」とする）は、メランコリー型MDDを健常者と識別できることを示している。
また、千葉において集められたコホート（慶応大学病院にて収集した。11名のメランコリー型MDDと40名の健常者を含む）についても同様に上記方法で取得した分類器を使って
検討した。その結果を図２７bに示す。千葉から集めたコホートの正確性は65％（感度64 %,特異度65 %, 及びAUC 0.62; 並べ替え検定においてp=0.036（図３０b））であった。本発明の分類器は、訓練データからは独立したコホートについてもメランコリー型MDDと健
常者を分類可能であった。
以上の結果から本発明の分類器は、メランコリー型MDDを分類可能であり、かつ汎化さ
れていると考えられた。
III.実施例２．本発明の分類器の非メランコリー型MDD群及び全MDD群への適用 図２７c
には、全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のサブタイプを含んでいると考えられた。
図２８d、e及び図２９f、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. 実施例３．WLSスコアを用いたうつ病の重症度の評価と治療効果の評価
WLSスコアがうつ病の重症度と相関するか否かを調べるため、BDスコアとWLSスコアの相関を検証した。図３１aに示すように、全MDD群と健常対照群を合わせた群についてBDIス
コアとWLSスコアの相関をとったところ、n=186 r = 0.655, 並び替え検定p = 0.001とな
り、両者は相関していることが示された。また、全MDD群のみをみても（図３１b）、n=93でr=0.188、並び替え検定p =0.046となり、比較的相関していた。
次に、選択的セロトニン再取り込み阻害薬（SSRI）であるエスシタロプラムの効果をWLSスコアで評価できるか検討した。エスシタロプラムを6〜8週間投与した、表１aに示す寛解した24名のメランコリー型MDD患者各々のWLSスコアの平滑化ヒストグラムを図３１cに
示す。抗うつ剤を投与するとWLSスコアの分布は健常対照群よりにシフトした。AUC 0.72
、Benjamini-Hochberg修正Kolmogorov-Smirnov検定p=0.008であった。
ただし、この結果は、集団についての結果であり、WLSと個人の改善との関係を示すも
のではない。そこで、ΔWLS（具体的には WLSpost - WLSpre）とΔBDI（具体的には BDIpost - BDIpre)の関係、ΔWLSとΔHAMD (具体的にはHAMDpost - HAMDpre)を検討した。有
意差は並び替え検定で求めた。図３１dに示すように、ΔBDIとΔWLSは有意に相関してい
た（r=0.373、p=0.040）が、ΔHAMDとΔWLSは相関していなかった（r=0.154、p=0.237）
。BDIとHAMDが乖離することは既に報告されている。
実施例５：12対の機能的結合それぞれの治療前後での変化
次に12対の機能的結合それぞれに対するＳＳＲＩ（エスシタロプラム）の影響について検討するため、個々の機能的結合についてΔWLSに対する寄与を検討した。
具体的には、投与前後解析において、式:
（式中、Nは治療を受け患者数を示す。wiは、分類器の重みを示す。）
を用いて各機能的結合の寄与スコアを算出した。
さらに、参考点を得るために、同様の寄与スコアをMDDと健常対照群の間でも式：
（式中、MはMDDの数を示し、Ｎは健常対照の数を示す。）
によって寄与スコアを算出した。
すなわち、siは分類器の重みで重み付けされた健常対照群とMDD群の平均結合強度の差
を示す。分類器において、MDDは陽性（正のスコア）を示し、健常者は陰性（負のスコア
）示すため、siは常に陰性（負のスコア）となる。また、式２において
は減数 (減算において負の符号)となる。ciとsiが大きな負の値をとることは治療後−治
療前健常者−MDDにおいてFCiの寄与が大きいことを示す。言い換えると、ciが負になることは 治療後のFCiが健常者の対応するFCにより近づいたことを示す。これに対して、ciが正になることは治療後のFCiがMDDとなったことを示す。健常対照群とMDD群で人数が異な
るため、それぞれのFCiについてのciとsiの有意差は、ウェルチのt-検定で評価した。ウ
ェルチのt-検定は、下式：
にしたがって、算出した。
各検定のp値は、Benjamini-Hochberg法によって複数の比較に批准させた。図３２aに示すように、FC1 (i.e. DLPFC-Precuneus/PCC)は、ci及びsiの差が最も大きく顕著であった。これは、抗うつ剤の投与後に FC1が健常対照群とは反対のMDD群の方に推移したことを
示す。
メランコリー型MDD分類のための12対の機能的結合（表Cに記載のこの機能的結合識別番号を「FC♯」と標記することがある）の中で8つの機能的結合については、健常対照群の
方向へ向かって推移し、いわゆる「正常化」した。しかし、FC# 1、3、9、12の４対の機
能的結合については、健常対照群とは逆の方向に向かって推移した（図３２a）。この中
で、FC# 1（Left DLPFC - Left Precuneus / PCC FC）は、12対の中で最も寄与度が高い
機能的結合である。メランコリー型MDD群と健常対照群の分類においてFC#1の治療前後の
変化の有意差は、ペアワイズ比較でp=0.009であった（図３２a）。
最も寄与度が高い、２つの主要な機能的結合（FC#1, FC#2）に着目した結果、治療後に有意な差が観察されたのに対して(対応のあるｔ検定, 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）。ここで、図３２
、３３の全ての図において、正方向への推移はメランコリー型MDDの方向へ向かっている
ことを示し、負方向への推移は健常者方向へ向かっていることを示す。これは、メランコリー型MDD分類器の重みの符号に、FCまたはFCの差を乗じているためである。上述のよう
に、エスシタロプラム治療後にFC# 1で正方向へ推移したため、完全に独立したコホート
でも同じ結果となるか検証した。千葉において収集された11名のメランコリー型MDDの患
者（セロトニン-ノルエピネフリン再取り込み阻害剤（SNRI）デュロキセチンを6〜8週間
投与）について検討した（図３２c）。治療前には有意差がなかった（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）、有意な関連は認められなかった（図３２d）。
さらに、広島コホート内の、臨床医によって治療効果を判定された患者について、非寛解患者（n=7）と寛解患者（n=7）に分け、抗うつ治療の効果を調べた。その結果、この結果、寛解・非寛解の群間 ｘ FC#1・FC#2の機能的結合間 ｘ 治療前後の時間に有意な交互作用が見られ (Fig. ３３e, 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次元プロットに表した（図３３f）。
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. 実施例４．ニューロフィードバックV-1.FCニューロフィードバック訓練に使用する標
的機能的結合の選択
以下の工程に従って、ニューロフィードバックの標的として、前記表５に示す第１の機能的結合を選択した。
はじめに、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を患う３名と、亜臨床的うつ病の７名の計10名が参加した。ニューロフ
ィードバック訓練に先立って、行われた２回の異なる時点における亜臨床的うつ病の参加者のBDI-IIの平均値は>10であった。
ニューロフィードバック訓練は、参加者にｆＭＲＩ装置に仰向けになってもらい、プリズム眼鏡を介して装置内のモニターを見せ、図３５に示す機能的結合の結合指標が負の方向に動くようにイメージさせることにより行った。
ニューロフィードバックスコアは、公知の方法に従って算出した。
図３７aにMDD参加者３名の各訓練日のニューロフィードバックスコアの集計を示す。スコアは、訓練の4日間にわたって増加傾向を示した。結果の有意性を、2つの説明変数（訓練日それぞれと、訓練項目）と１つの応答変数（ニューロフィードバックスコア）を含む重回帰モデルを用いて確認したところ、訓練日において有意に良い効果を示した（95% 信頼区間（CI）計数 1.9-9.1）。事後比較として行ったt-検定では、MDD参加者３名全てに
おいて、訓練の初日と比較して、最終日で有意に高いニューロフィードバックスコアが得られた。(t = 4.01, P < .001). これらの結果は、MDDの参加者全てが、訓練を通して標
的FCの負の相関を誘導することを学んだこと示している。さらに、図３７bに示す用に、MDD参加者３名全てにおいて、訓練後に、うつ症状の重症度を示すHamilton Rating Scale for Depression（HAMD）が減少した。
MDD参加者と同様に、亜臨床的うつ病の参加者７名においても、訓練期間中、ニューロ
フィードバックスコアが増加する傾向にあった（図３７c）。この結果について有意性を
解析したところ、一元配置分散分析（One-way ANOVA）で訓練日の主効果が有意であった(p = 0.0046)。事後比較ペアt検定では、訓練の最終日にニューロフィードバックスコアが訓練の初日より有意に高かった。訓練後のBDIスコアも減少傾向を示した（p = .07）。また、図３７dに示すように、亜臨床的うつ病の参加者７名中５名については、正常方向に
標的rs-fcMRIを変化させ、訓練前後間でのrs-fcMRI変化量は、BDIスコアと有意に相関し
た（r = 0.87, p = 0.011）。
上述のように、半数以上の参加者が、訓練を通じてBDIスコアと相関する表５に示す機
能的結合を減少させたことから、表５に示す機能的結合は、うつ症状を改善するためのニューロフィードバック訓練の標的機能的結合として使用できることが示された。さら、本発明の表５に示す機能的結合は、うつ症状を改善するためのニューロフィードバック訓練法を開発するための治療方法バイオマーカーとして使用できることが示された。
V-2.他の抗うつ治療との比較
ニューロフィードバック訓練法によるMDDの治療効果を、他の治療方法（抗うつ薬投与
、反復経頭蓋磁気刺激療法：rTMS、修正型電気けいれん療法：ECT intervention）と比較した。ニューロフィードバック訓練法は、他のうつ治療よりも治療効果が改善していた。
２ 被験者、６ ディスプレイ、１０ ＭＲＩ装置、１１ 磁場印加機構、１２ 静磁場発生コイル、１４ 傾斜磁場発生コイル、１６ ＲＦ照射部、１８ 寝台、２０ 受信コイル、２１ 駆動部、２２ 静磁場電源、２４ 傾斜磁場電源、２６ 信号送信部、２８ 信号受信部、３０寝台駆動部、３２ データ処理部、３６ 記憶部、３８ 表示部、４０ 入力部、４２ 制御部、４４ インタフェース部、４６ データ収集部、４８ 画像処理部、５０ ネットワークインタフェース。

Claims (11)

  1. 判別装置であって、
    演算装置を備え、前記演算装置は、
    複数の対象者について、下記表１に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関とがそれぞれ計測されて得られる情報に基づき、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する分類器生成処理を実行し、
    第１の時点の被検者の安静時の脳内における、前記複数の機能的結合の第１の相関を算出し、
    第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第２相関を算出し、
    前記被検者の前記複数の機能的結合の前記第２の相関と前記第１の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する、うつ病患者における治療効果の判定処理を実行する、ように構成された、判別装置：
    。
  2. 前記複数の機能的結合は、
    機能的結合識別番号１の機能的結合、及び機能的結合識別番号２の機能的結合である、請求項１の判別装置。
  3. 前記うつ病は、メランコリー型うつ病である、請求項１又は２に記載の脳活動訓練装置。
  4. さらに、ドラッグリプロファイリングのために使用される、請求項１から３のいずれか一項に記載の判別装置。
  5. 請求項１から４のいずれか一項に記載の判別装置の機能をコンピュータに実行させるためのプログラム。
  6. 複数の対象者について、下記表２に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関とがそれぞれ計測されて得られる情報に基づき、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する工程と、
    第１の時点の被検者の安静時の脳内における、前記複数の機能的結合の第１の相関を計測する工程と、
    第２の時点の前記被検者と同一の被検者の安静時の脳内における、前記複数の機能的結合の第２相関を計測する工程と、
    前記被検者の前記複数の機能的結合の前記第２の相関と前記第１の相関との差に基づき、前記分類器により、前記被検者に対する治療効果を判別する工程と、
    を含む、うつ病患者における治療効果の判定を補助するための判別方法：
    。
  7. ドラッグリプロファイリングのために使用される、請求項６に記載の判別方法。
  8. 分類器生成装置であって、
    演算装置を備え、前記演算装置は、
    複数の対象者について、下記表３に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関をそれぞれ計測して、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を生成し、
    記憶装置をさらに備え、前記記憶装置は、前記演算装置によって生成された前記分類器を特定するための情報を記憶する、分類器生成装置：
    。
  9. 前記対象者は、うつ病に伴う症状を有する、請求項８に記載の脳活動訓練装置。
  10. 前記うつ症状は、メランコリー型うつ病に伴う症状である、請求項９に記載の脳活動訓練装置。
  11. 複数の対象者について、下記表４に示す機能的結合識別番号１〜１２から選択される複数の機能的結合の第１の時点の相関と治療開始後であって第１の時点よりも後の第２の時点の相関をそれぞれ計測して、前記複数の機能的結合の第１の時点と第２の時点の前記相関の差の各々で張られる相関状態空間において、前記複数の対象者において治療効果がみられた群と治療効果がみられない群とを識別する分類器を予め生成する工程を備える、分類器生成方法：
    。
JP2019061843A 2017-10-03 2019-03-27 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置 Active JP7140348B2 (ja)

Applications Claiming Priority (2)

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

Related Parent Applications (1)

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

Publications (3)

Publication Number Publication Date
JP2019155111A true JP2019155111A (ja) 2019-09-19
JP2019155111A5 JP2019155111A5 (ja) 2021-02-25
JP7140348B2 JP7140348B2 (ja) 2022-09-21

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 Before (1)

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

Country Status (4)

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

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111588351A (zh) * 2020-06-09 2020-08-28 天津大学 一种基于功能近红外脑网络的中医针刺得气智能检测系统

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200410890A1 (en) * 2018-03-09 2020-12-31 Advanced Telecommunications Research Institute International Brain activity training apparatus, brain activity training method and brain activity training program
JP7147623B2 (ja) * 2019-02-21 2022-10-05 株式会社島津製作所 脳血流量の特徴量の抽出方法
BR112021018770A2 (pt) * 2019-03-22 2022-02-15 Cognoa Inc Métodos e dispositivos de terapia digital personalizada
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
US20230107263A1 (en) * 2020-04-06 2023-04-06 Advanced Telecommunications Research Institute International Brain functional connectivity correlation value clustering device, brain functional connectivity correlation value clustering system, brain functional connectivity correlation value clustering method, brain functional connectivity correlation value classifier program, brain activity marker classification system and clustering classifier model for brain functional connectivity correlation values
EP4183333A4 (en) * 2020-07-17 2023-12-20 Advanced Telecommunications Research Institute International BRAIN FUNCTIONAL CONNECTIVITY CORRELATION VALUES GROUPING DEVICE, BRAIN FUNCTIONAL CONNECTIVITY CORRELATION VALUES GROUPING SYSTEM, BRAIN FUNCTIONAL CONNECTIVITY CORRELATION VALUES GROUPING METHOD, VAL CLASSIFICATION PROGRAM BRAIN FUNCTIONAL CONNECTIVITY CORRELATION ORS, AND SYSTEM OF BRAIN ACTIVITY MARKER CLASSIFICATION
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
CN113947580B (zh) * 2021-10-18 2022-07-19 上海市精神卫生中心（上海市心理咨询培训中心） 针对精神疾病实现神经影像辅助诊断处理的系统
CN114860922B (zh) * 2022-03-25 2023-04-18 南京脑科医院 心理测评量表的分类模型获取方法、筛查方法及系统
CN115147417B (zh) 2022-09-02 2022-11-15 之江实验室 一种基于过滤法特征选择的功能连接矩阵处理系统及装置
CN116469541B (zh) * 2023-03-24 2024-01-30 深圳市第二人民医院(深圳市转化医学研究院) 一种抑郁症标志物及其在抑郁症预后中的应用、评估装置
CN116504363B (zh) * 2023-05-10 2024-03-08 中国科学院心理研究所 基于组水平差异统计图的抑郁症tms个体化靶点定位方法及系统
CN116543042B (zh) * 2023-05-10 2024-01-23 中国科学院心理研究所 基于组水平平均统计图的抑郁症tms个体化靶点定位方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008531158A (ja) * 2005-03-04 2008-08-14 メンティス キュラ イーエイチエフ． 神経学的状態を評価する方法および系
JP2011520095A (ja) * 2008-03-12 2011-07-14 リッジ ダイアグノスティックス，インコーポレイテッド うつ障害をモニタリングするための炎症バイオマーカー
JP2012520499A (ja) * 2009-03-12 2012-09-06 コミッサリア ア レネルジー アトミーク エ オ ゼネルジ ザルタナテイヴ 情報予測デバイスを構築する方法、情報予測デバイスの使用方法、ならびに対応する記憶媒体および記憶装置
WO2014178322A1 (ja) * 2013-05-01 2014-11-06 株式会社国際電気通信基礎技術研究所 脳活動訓練装置および脳活動訓練システム
US20150119689A1 (en) * 2012-05-16 2015-04-30 Beth Israel Deaconess Medical Center, Inc. Identifying individual target sites for transcranial magnetic stimulation applications
JP6195329B1 (ja) * 2015-11-24 2017-09-13 株式会社国際電気通信基礎技術研究所 脳活動解析装置、脳活動解析方法、プログラムおよびバイオマーカー装置

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006102370A2 (en) * 2005-03-22 2006-09-28 University Of Massachusetts Functional brain mri mapping as a marker in cns diseases and disorders
US8597614B2 (en) * 2005-08-03 2013-12-03 Centre For Addiction And Mental Health Diagnosis and treatment of mood disorders
EP2158574B1 (en) * 2007-06-21 2018-12-12 Koninklijke Philips N.V. Model-based differential diagnosis of dementia and interactive setting of level of significance
WO2009111595A2 (en) * 2008-03-04 2009-09-11 Ridge Diagnostics, Inc. Diagnosing and monitoring depression disorders based on multiple biomarker panels
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
US20130102918A1 (en) * 2011-08-16 2013-04-25 Amit Etkin System and method for diagnosing and treating psychiatric disorders
WO2014190414A1 (en) * 2013-05-28 2014-12-04 Laszlo Osvath Systems and methods for diagnosis of depression and other medical conditions
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
WO2015134762A1 (en) * 2014-03-07 2015-09-11 The Research Foundation For The State University Of New York Method of diagnosing depression by pet imaging
WO2015164477A1 (en) * 2014-04-25 2015-10-29 The General Hospital Corporation Hybrid system for treating mental and emotional disorders with responsive brain stimulation
DE102015109853A1 (de) * 2015-06-19 2016-12-22 Nclogics Ag Assistenz- und Entscheidungssystem und Verfahren zur Auswertung von Elektroenzephalogrammen
WO2018023056A1 (en) * 2016-07-29 2018-02-01 Beth Israel Deaconess Medical Center, Inc. Methods and apparatus for network localization of neurological symptoms from focal brain lesions
CN106407733A (zh) * 2016-12-12 2017-02-15 兰州大学 基于虚拟现实场景脑电信号的抑郁症风险筛查系统和方法
CN106725532B (zh) * 2016-12-13 2018-04-24 兰州大学 基于语音特征与机器学习的抑郁症自动评估系统和方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008531158A (ja) * 2005-03-04 2008-08-14 メンティス キュラ イーエイチエフ． 神経学的状態を評価する方法および系
JP2011520095A (ja) * 2008-03-12 2011-07-14 リッジ ダイアグノスティックス，インコーポレイテッド うつ障害をモニタリングするための炎症バイオマーカー
JP2012520499A (ja) * 2009-03-12 2012-09-06 コミッサリア ア レネルジー アトミーク エ オ ゼネルジ ザルタナテイヴ 情報予測デバイスを構築する方法、情報予測デバイスの使用方法、ならびに対応する記憶媒体および記憶装置
US20150119689A1 (en) * 2012-05-16 2015-04-30 Beth Israel Deaconess Medical Center, Inc. Identifying individual target sites for transcranial magnetic stimulation applications
WO2014178322A1 (ja) * 2013-05-01 2014-11-06 株式会社国際電気通信基礎技術研究所 脳活動訓練装置および脳活動訓練システム
JP6195329B1 (ja) * 2015-11-24 2017-09-13 株式会社国際電気通信基礎技術研究所 脳活動解析装置、脳活動解析方法、プログラムおよびバイオマーカー装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111588351A (zh) * 2020-06-09 2020-08-28 天津大学 一种基于功能近红外脑网络的中医针刺得气智能检测系统

Also Published As

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

Similar Documents

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

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210115

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210115

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20211221

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220208

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220408

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220826

R150 Certificate of patent or registration of utility model

Ref document number: 7140348

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150