CN105163659B - 脑活动分析装置、脑活动分析方法以及生物标记物装置 - Google Patents

脑活动分析装置、脑活动分析方法以及生物标记物装置 Download PDF

Info

Publication number
CN105163659B
CN105163659B CN201480024893.1A CN201480024893A CN105163659B CN 105163659 B CN105163659 B CN 105163659B CN 201480024893 A CN201480024893 A CN 201480024893A CN 105163659 B CN105163659 B CN 105163659B
Authority
CN
China
Prior art keywords
mentioned
cerebration
arbiter
subject
disease
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201480024893.1A
Other languages
English (en)
Other versions
CN105163659A (zh
Inventor
森本淳
川人光男
八幡宪明
桥本龙郎
桥本龙一郎
柴田和久
渡边武郎
佐佐木由香
加藤进昌
笠井清登
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ATR Advanced Telecommunications Research Institute International
Original Assignee
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 ATR Advanced Telecommunications Research Institute International filed Critical ATR Advanced Telecommunications Research Institute International
Publication of CN105163659A publication Critical patent/CN105163659A/zh
Application granted granted Critical
Publication of CN105163659B publication Critical patent/CN105163659B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0263—Measuring blood flow using NMR
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4058—Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
    • A61B5/4064—Evaluating the brain
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
    • A61B5/4088—Diagnosing of monitoring cognitive diseases, e.g. Alzheimer, prion diseases or dementia
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4806—Functional imaging of brain activation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
    • A61B2560/04—Constructional details of apparatus
    • A61B2560/0475—Special features of memory means, e.g. removable memory cards
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis
    • A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/369—Electroencephalography [EEG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4848—Monitoring or testing the effects of treatment, e.g. of medication
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4808—Multimodal MR, e.g. MR combined with positron emission tomography [PET], MR combined with ultrasound or MR combined with computed tomography [CT]
    • 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

Abstract

提供一种用于针对神经/精神疾病实现利用脑功能成像法的生物标记物的脑活动分析方法。根据在健康组、患者组中测定出的静息态功能连接MRI的数据，针对各个受验者导出规定的脑区域间的活动度的相关矩阵(80)。通过正则化典型相关分析(82)对包括受验者的疾病/健康标签在内的受验者的属性以及相关矩阵(80)进行特征提取(84)。基于正则化典型相关分析的结果，通过利用稀疏逻辑回归的判别分析(86)来生成判别器(88)。

Description

脑活动分析装置、脑活动分析方法以及生物标记物装置
技术领域
本发明涉及一种使用脑功能成像法的脑活动分析装置、脑活动分析方法以及生物标记物装置。
背景技术
(生物标记物(Biomarker))
将为了定量地掌握生物体内的生物学变化而将生物体信息数值化、定量化得到的指标称为“生物标记物”。
FDA(美国食品药品局)将生物标记物定位为“作为正常生理、病理生理或者对治疗的药理学反应的指标来客观地测定、评价的项目”。另外，对疾病的状态、变化或治愈的程度赋予特征的生物标记物被用作用于确认新药在临床试验中的有效性的替代标记物(surrogate marker)。血糖值和胆固醇值等是作为生活习惯病的指标而具有代表性的生物标记物。不仅包括尿或血液中包含的来自生物体的物质，还包括心电图、血压、PET图像、骨密度、肺功能等。另外，通过基因组分析(genomic analysis)和蛋白质组分析(proteomeanalysis)的不断进步，发现了与DNA、RNA或生物体蛋白等相关联的各种生物标记物。
关于生物标记物，不仅期待其应用于患病后的治疗效果的测定，还期待其作为用于防疾病于未然的日常指标来应用于疾病的预防，并且期待其在选择避免副作用的有效的治疗方法的个别化医疗中的应用。
但是，在神经/精神疾病的情况下，尽管还对从生物化学或分子遗传学的观点来看能够用作客观指标的分子标记物等进行了研究，但是处于应该称为探讨阶段的状况。
另一方面，也有对如下的疾病判定系统等的报告：使用NIRS(Near-infraRedSpectroscopy：近红外光谱)技术，根据通过生物体光测量而测量出的血红蛋白信号的特征量，来对精神分裂症(schizophrenia)、抑郁症等精神疾病进行分类(专利文献1)。
(实时神经反馈(Real Time Neurofeedback))
例如，以往对于作为神经症的一个类型的强迫性障碍(Obsessive-CompulsiveDisorder：OCD)，作为其治疗方法，已知药物疗法和行为疗法。作为药物疗法，例如使用选择性血清素再摄取抑制剂，作为行为疗法，已知将暴露法与反应阻止法结合而成的暴露反应阻止法等。
另一方面，作为神经/精神疾病的情况下的一种可能的治疗方法，探讨了实时神经反馈法。
以利用核磁共振成像法(MRI：Magnetic Resonance Imaging)将与人脑的活动相关联的血流动态反应视觉化的方法、即功能性磁共振成像法(fMRI：functional MagneticResonance Imaging)为首的脑功能成像法用于检测因感觉刺激或执行认知任务引起的脑活动与静息态或因执行对照任务引起的脑活动之间的不同来确定与所关心的脑功能的构成要素对应的脑活化区域、即明确脑的功能定位。
近年来，有对使用功能性磁共振成像法(fMRI)等脑功能成像法的实时神经反馈技术的报告(非专利文献1)。实时神经反馈技术作为神经疾病、精神疾病的治疗方法开始受到关注。
神经反馈是生物反馈的一种，受验者通过接受与自身的脑活动有关的反馈来学习操作脑活动的方法。
例如有以下报告：当利用fMRI测量前扣带回的活动并将其实时地作为火焰的大小反馈给患者来使患者努力减小火焰时，中枢性的慢性疼痛既实时地得到改善而且也长期地得到改善(参照非专利文献2)。
(静息态fMRI)
另外，在最近的研究中，获知了静息态的脑也在活跃地进行活动。也就是说，脑内存在一种在脑的主动活动时镇静化、在脑休息时活跃地兴奋的神经细胞群。该神经细胞群在解剖学上主要位于左右的大脑相接合的内侧面、即前额叶内侧面、后扣带回、楔前叶、顶叶联合区的后半部、颞中回等。表示该静息态的基线脑活动的区域被命名为默认模式网络(Default Mode Network；DMN)，作为一个网络而同步地活动(参照非专利文献3)。
例如，作为正常人与精神疾病患者的脑活动不同的例子，能够列举出默认模式网络中的脑活动。默认模式网络是指在静息状态下能够观察到比正在执行目标指向性的任务时更活跃的脑活动的部位。做出过以下报告：与正常人相比，在精神分裂症、老年痴呆症等疾病的患者的该默认模式网络中能够观察到异常。例如，有以下报告等：在静息状态下，精神分裂症患者的属于默认模式网络的后扣带回皮质与外侧顶叶皮质、内侧前额叶皮质或小脑皮质之间的活动的相关性下降。
但是，处于以下状况：这种默认模式网络与认知功能的关系如何以及这种脑区域之间的功能性连接的相关性与上述的神经反馈之间的关联等未必明确。
另一方面，通过观察多个脑区域之间因任务等的不同引起的活动的相关关系的变化，来评价这些脑区域间的功能性连接(functional connectivity)。特别是，静息态下的通过fMRI得到的功能性连接的评价也被称为静息态功能连接MRI(rs-fcMRI：resting-state functional connectivity MRI)，也逐渐进行以各种神经/精神疾病为对象的临床研究。但是，以往的rs-fcMRI法用于观察如上所述那样的默认模式网络之类的全局性神经网络的活动，处于不能说对更详细的功能性连接进行了充分研究的状况。
(DecNef法：Decoded NeuroFeedback法)
另一方面，近年来，有对被称为解码神经反馈法(DecNef法)的新类型的神经反馈法的报告(非专利文献4)。
人的感觉和知觉系统根据环绕周围的环境而始终变化。这种变化大半是在发育早期的被确定的阶段、即被称为“临界期”的时期发生。但是，即使是成人，也以能够适应周边环境的重要变化的程度保持着感觉和知觉系统的可塑性。例如，做出过以下报告：成人通过接受使用特定的知觉刺激的训练或者被暴露于特定的知觉刺激下，其训练任务的成绩或对知觉刺激的灵敏度提高，并且，其训练结果会维持几个月到几年(例如参照非专利文献5)。这种变化被称为知觉学习，确认出这种变化在所有感官即视觉、听觉、嗅觉、味觉、触觉上均会发生。
在DecNef法中，不是对受验者直接施加作为学习对象的刺激，而是检测脑的活动，并且对脑活动进行解码，只将与期望的脑活动之间的近似度反馈给受验者，由此能够进行“知觉学习”。
(核磁共振成像法)
简单地说明这种核磁共振成像法则如下。
即，以往以来，在人的临床成像诊断等中使用利用对于生物体中的原子、特别是氢原子的原子核的核磁共振现象的核磁共振成像法来作为对生物体的脑或全身的断面进行成像的方法。
核磁共振成像法在应用于人体的情况下，与同样是人体内断层成像法的“X射线CT”相比，例如存在如下的特征。
(1)能够得到与氢原子的分布及其信号弛豫时间(反映原子结合的强度)对应的浓度的图像。因此，呈现出与组织的特性差异相应的浓淡，从而易于观察组织的不同。
(2)磁场不会被骨骼吸收。因此，易于观察被骨骼包围的部位(头颅内、脊髓等)。
(3)不会像X射线那样对人体有害，因此能够广泛地有效利用。
这种核磁共振成像法利用人体的各细胞中含有量最多的、且具有最大磁性的氢原子核(质子)的磁特性。承担氢原子核的磁性的自旋角动量的磁场内的运动在古典意义上被比喻为陀螺的进动。
下面，为了说明本发明的背景，以该直观的古典模型(classical model)来简单地总结核磁共振的原理。
如上所述的氢原子核的自旋角动量的方向(陀螺的自转轴的方向)在没有磁场的环境中朝向随机的方向，但是若施加静磁场，则会朝向磁力线的方向。
当在该状态下进一步叠加振动磁场时，若该振动磁场的频率是根据静磁场的强度决定的共振频率f0＝γB0/2π(γ：物质所固有的系数)，则通过共振而在原子核侧能量发生移动，磁化矢量的方向发生变化(进动变大)。当在该状态下切断振动磁场时，进动退回倾斜角度来逐渐恢复为在静磁场中的方向。通过从外部利用天线线圈探测该过程，能够得到NMR信号。
在静磁场的强度为B0(T)时，在氢原子的情况下这种共振频率f0为42.6×B0(MHz)。
并且，在核磁共振成像法中，还能够利用与血流量的变化相应地在被检测的信号中出现变化这一情况，来将对于外部刺激等的脑的活动部位视觉化。将这种核磁共振成像法特别称为fMRI(functional MRI)。
在fMRI中，作为装置，使用在通常的MRI装置中进一步配备fMRI测量所需的硬件和软件的装置
在此，血流量的变化对NMR信号强度带来变化是由于血液中的氧合血红蛋白与脱氧血红蛋白的磁特性不同。氧合血红蛋白有抗磁性体的特性，对周围存在的水的氢原子的弛豫时间没有影响。与此相对，脱氧血红蛋白是顺磁性体，改变周围的磁场。因而，当脑受到刺激而局部血流增大、氧合血红蛋白增加时，能够检测其变化量来作为MRI信号。作为对受验者的刺激，例如使用视觉上的刺激或听觉上的刺激、或者规定的任务(task)的执行等(例如专利文献2)。
在此，在脑功能研究中，通过测定与微小静脉或毛细血管中的红细胞中的脱氧血红蛋白的浓度减少的现象(BOLD效应)对应的、氢原子的核磁共振信号(MRI信号)的上升，来进行脑的活动的测定。
特别是，在关于人的运动功能的研究中，一边使受验者进行某些运动，一边通过上述MRI装置测定脑的活动。
另外，在人的情况下，需要无创性的脑活动测量，在该情况下，能够从fMRI数据中提取出更详细的信息的解码技术发展起来(例如非专利文献6)。特别是，fMRI以脑中的像素(volumetric pixel：voxel(体素))为单位对脑活动进行分析，由此能够根据脑活动的空间图案来估计刺激输入和识别状态。因此，这种解码技术对于针对知觉学习的任务的应用即为上述的DecNef法。
专利文献1：日本再表2006-132313号公报
专利文献2：日本特开2011-000184号公报
非专利文献1：Nikolaus Weiskopf，“Real-time fMRI and its application toneurofeedback”，NeuroImage 62(2012)682-692
非专利文献2：deCharms RC,Maeda F,Glover GH et al,“Control overbrainactivation and pain learned by using real-time functional MRI”，ProcNatlAcad Sci USA 102(51),18626-18631,2005
非专利文献3：Raichle ME,Macleod AM,Snyder AZ.et.al.“A default modeofbrain function”，Proc Natl Acad Sci USA 98(2)，676-682，2001
非专利文献4：Kazuhisa Shibata,Takeo Watanabe,Yuka Sasaki,MitsuoKawato,“PerceptualLearning Incepted by Decoded fMRI Neurofeedback WithoutStimulus Presentation”,SCIENCE VOL 3349DECEMBER 2011
非专利文献5：T.Watanabe,J.E.Nanez Sr,S.Koyama,I.Mukai,J.Liederman andY.Sasaki:Greater plasticity in lower-level than higher-level visualmotionprocessing in a passive perceptual learning task.NatureNeuroscience,5,1003-1009,2002.
非专利文献6：Kamitani Y,Tong F.Decoding the visual and subjectivecontents ofthe human brain.Nat Neurosci.2005；8:679-85.
发明内容
发明要解决的问题
如以上那样，指出了利用功能性磁共振成像法等脑功能成像法进行的脑活动分析以及使用该分析的神经反馈的技术被应用于一部分神经/精神疾病的治疗等的可能性，但是现状是尚未达到实用性的应用阶段。
在考虑应用于神经/精神疾病的治疗的情况下，作为如上所述的生物标记物，利用脑功能成像法进行的脑活动分析作为无创性的功能标记物而例如还被期待应用于诊断法的开发、面向用于实现根本治疗的药物研发的靶分子的搜索和识别等。
例如，到目前为止，对于自闭症等精神疾病，使用基因的实用性的生物标记物尚未完成，因此，药物的效果判定等是困难的，因此治疗药物的开发也是困难的。
另一方面，有如下启示：能够基于根据静息态的fMRI数据导出的脑区之间的连接来在某种程度上预测神经疾病的诊断结果。但是，在这些研究中也是关于该预测性能的验证仅使用了在一个设施中测量出的脑活动，作为生物标记物的有用性的验证并不充分。
因而，以往并不清楚如何构成这种利用脑功能成像法的生物标记物更好。
本发明是为了解决如上所述的问题而完成的，其目的在于实现一种能够提供用于客观地判断脑的活动状态是健康状态还是疾病状态的数据的脑活动分析装置和脑活动分析方法。
本发明的其它目的在于提供一种用于实现用于辅助针对神经/精神疾病的诊断的利用脑功能成像法的判别处理的脑活动分析装置和脑活动分析方法。
本发明的其它目的在于提供一种用于实现利用脑功能成像法的生物标记物的脑活动分析装置、脑活动分析方法以及生物标记物装置。
用于解决问题的方案
按照本发明的一个方面，是一种脑活动分析装置，具备：判别器生成单元，其用于根据通过脑活动探测装置按时间序列预先测定表示多个受验者各自的脑内的多个规定区域中的脑活动的信号而得到的信号来生成判别器，判别器生成单元提取多个规定区域中的脑活动的相关关系中的至少在多个受验者的属性上共通的缩略表达(日语：縮約表現)，对于提取出的缩略表达，生成针对受验者的属性中的特定属性的判别器；存储单元，其用于存储用于确定判别器的信息；以及判别单元，其基于根据存储单元中存储的信息而确定的判别器，来对输入数据进行判别处理。
优选的是，特定属性是疾病的判别标签。
优选的是，受验者的属性包括向该受验者给予的药物的标签。
优选的是，脑活动探测装置包括多个脑活动测定装置，判别器生成单元包括提取单元，该提取单元从多个规定区域中的脑活动的相关关系中通过变量选择来提取在多个脑活动测定装置的测定条件和多个受验者的属性上共通的缩略表达。
优选的是，判别器生成单元包括回归单元，该回归单元通过对提取出的缩略表达进一步进行变量选择的回归来生成判别器。
优选的是，提取单元包括相关分析单元，该相关分析单元用于根据由脑活动探测装置探测出的信号来计算多个规定区域间的活动的相关矩阵，在受验者的属性与相关矩阵的非对角元素之间执行正则化典型相关分析来提取缩略表达。
优选的是，回归单元包括回归分析单元，该回归分析单元通过对正则化典型相关分析的结果和受验者的属性进行稀疏逻辑回归来生成判别器。
优选的是，多个脑活动测定装置是分别设置于多个不同的场所的、用于通过脑功能成像法来按时间序列测量脑活动的装置。
优选的是，判别处理是疾病的判别标签的判别，该疾病的判别标签表示与输入数据对应的受验者患有还是不患有神经/精神疾病。
优选的是，受验者的属性包括表示患有还是不患有神经/精神疾病的疾病的判别标签、表示受验者的个人特性的标签、以及对利用脑活动探测装置进行的测定赋予特征的信息，判别处理是疾病的判别标签的判别，该疾病的判别标签表示与输入数据对应的受验者患有还是不患有神经/精神疾病。
优选的是，正则化典型相关分析是基于L1正则化的典型相关分析。
优选的是，脑活动测定装置拍摄静息态功能连接核磁共振影像。
按照本发明的其它方面，是一种脑活动分析方法，用于使具备运算装置和存储装置的计算机进行脑活动的分析，该脑活动分析方法具备以下步骤：根据通过脑活动探测装置按时间序列预先测定表示多个受验者各自的脑内的多个规定区域中的脑活动的信号而得到的信号来生成判别器的步骤，生成判别器的步骤包括以下步骤：提取多个受验者的多个规定区域中的脑活动的相关关系中的至少在多个受验者的属性上共通的缩略表达，对于提取出的缩略表达，生成针对受验者的属性中的特定属性的判别器；将用于确定判别器的信息存储在存储装置中的步骤；以及基于根据存储装置中存储的信息而确定的判别器来对输入数据进行判别处理的步骤。
优选的是，特定属性是疾病的判别标签。
优选的是，受验者的属性包括向该受验者给予的药物的标签。
优选的是，脑活动探测装置包括多个脑活动测定装置，生成判别器的步骤包括以下步骤：从多个规定区域中的脑活动的相关关系中通过变量选择来提取在多个脑活动测定装置的测定条件和多个受验者的属性上共通的缩略表达。
优选的是，生成判别器的步骤包括以下步骤：通过对提取出的缩略表达进一步进行变量选择的回归来生成判别器。
按照本发明的又一方面，是一种生物标记物装置，通过计算机对脑活动的分析来生成作为生物标记物的输出，该生物标记物装置具备：存储装置，其用于存储用于确定判别器的信息，判别器是根据通过脑活动探测装置按时间序列预先测定表示多个受验者各自的脑内的多个规定区域中的脑活动的信号而得到的信号来计算出的判别器，其中，提取多个受验者的多个规定区域中的脑活动的相关关系中的至少在多个受验者的属性上共通的缩略表达，对于提取出的缩略表达，针对受验者的属性中的疾病的判别标签来生成判别器；以及运算装置，其基于根据存储单元中存储的信息而确定的判别器，来对输入数据进行判别处理。
优选的是，受验者的属性包括向该受验者给予的药物的标签。
发明的效果
根据本发明，能够实现能够提供用于客观地判断脑的活动状态是健康状态还是疾病状态的数据的脑活动分析装置和脑活动分析方法。
或者，根据本发明，能够针对神经/精神疾病实现利用脑功能成像法的生物标记物。
附图说明
图1是表示MRI装置10的整体结构的示意图。
图2是数据处理部32的硬件框图。
图3是表示通过本实施方式的rs-fcMRI法拍摄的脑的关注区域(ROI：Region ofInterest)的图。
图4是表示针对关注区域提取表示静息态的功能连接的相关性的相关矩阵的过程的概念图。
图5是说明根据相关矩阵来生成成为生物标记物的判别器的过程的概念图。
图6是用于说明为了生成成为生物标记物的判别器而由数据处理部32进行的处理的流程图。
图7是表示通过稀疏典型相关分析(Sparse Canonical Correlation Analysis：SCCA)得到的特征量(通过SCCA得到的特征空间)以及以此为输入通过稀疏逻辑回归(Sparce Logistic Regression：SLR)来生成判别器的过程的概念的图。
图8是表示通过SCCA得到的特征量(通过SCCA得到的特征空间)以及以此为输入通过SLR来生成判别器的过程的概念的图。
图9是表示通过SCCA得到的特征量(通过SCCA得到的特征空间)以及以此为输入通过SLR来生成判别器的过程的概念的图。
图10是表示生物标记物的生成过程的概念的图。
图11是表示所生成的生物标记物的验证过程的概念图。
图12是表示生物标记物的特性的图。
图13是表示多疾病生物标记物的概念的图。
具体实施方式
下面，按图来说明本发明的实施方式的MRI系统的结构。此外，在下面的实施方式中，标注了相同标记的结构要素和处理步骤是相同或相当的，在不必要的情况下不重复其说明。
[实施方式1]
图1是表示MRI装置10的整体结构的示意图。
如图1所示，MRI装置10具备：磁场施加机构11，其向受验者2的关注区域赋予受控制的磁场并照射RF(射频)波；接收线圈20，其接收来自该受验者2的响应波(NMR信号)并输出模拟信号；驱动部21，其对向该受验者2赋予的磁场进行控制，并且对RF波的发送接收进行控制；以及数据处理部32，其设定该驱动部21的控制时序，并且对各种数据信号进行处理来生成图像。
此外，在此，取载置受验者2的圆筒形状的孔(Bore)的中心轴为Z轴，将与Z轴正交的水平方向定义为X轴，将与Z轴正交的垂直方向定义为Y轴。
由于MRI装置10是这样的结构，因此通过由磁场施加机构11施加的静磁场，构成受验者2的原子核的核自旋顺着磁场方向(Z轴)取向，并且以该原子核所固有的拉莫尔频率(Larmor frequency)进行以该磁场方向为轴的进动。
而且，当照射与该拉莫尔频率相同频率的RF脉冲时，原子发生共振并吸收能量而被激励，从而产生核磁共振现象(NMR现象；Nuclear Magnetic Resonance)。在该共振后，若停止RF脉冲照射，则原子放出能量来恢复为原来的稳定状态。将该过程称为弛豫过程。在弛豫过程中原子输出与拉莫尔频率相同的频率的电磁波(NMR信号)。
利用接收线圈20接收所输出的该NMR信号来作为来自受验者2的响应波，在数据处理部32中，受验者2的关注区域被图像化。
磁场施加机构11具备静磁场发生线圈12、倾斜磁场发生线圈14、RF照射部16以及在孔中载置受验者2的床18。
受验者2在床18上例如仰卧。虽未特别限定，但受验者2例如能够通过棱镜眼镜4看到与Z轴垂直地设置的显示器6上显示的画面。通过该显示器6的图像来对受验者2施加视觉刺激。此外，对受验者2的视觉刺激也可以是通过投影仪将图像投影到受验者2的眼前的结构。
这种视觉刺激在上述的神经反馈中相当于反馈信息的呈现。
驱动部21具备静磁场电源22、倾斜磁场电源24、信号发送部26、信号接收部28以及使床18移动到Z轴方向上的任意位置的床驱动部30。
数据处理部32具备：输入部40，其从操作者(省略图示)受理各种操作和信息输入；显示部38，其对与受验者2的关注区域有关的各种图像和各种信息进行画面显示；显示控制部34，其控制显示部38的显示；存储部36，其存储用于执行各种处理的程序、控制参数、图像数据(构造图像等)以及其它电子数据；控制部42，其对产生使驱动部21驱动的控制时序等的各功能部的动作进行控制；接口部44，其与驱动部21之间执行各种信号的发送接收；数据收集部46，其收集包括源于关注区域的一组NMR信号的数据；图像处理部48，其基于该NMR信号的数据来形成图像；以及网络接口50，其用于执行与网络之间的通信。
另外，数据处理部32除了是专用计算机的情况以外，也包括如下情况：是执行使各功能部动作的功能的通用计算机，基于安装于存储部36的程序来进行所指定的运算、数据处理以及控制时序的产生。下面，设数据处理部32是通用计算机来进行说明。
静磁场发生线圈12向绕Z轴卷绕的螺旋线圈流通从静磁场电源22提供的电流来使螺旋线圈产生感应磁场，使得在孔中产生Z轴方向的静磁场。将受验者2的关注区域设定在形成于该孔的静磁场的均匀性高的区域。在此，更为详细地说，静磁场发生线圈12例如包括四个空心线圈，利用其组合来在内部生成均匀的磁场，对受验者2的体内的规定的原子核、更特定地说是对氢原子核的自旋(spin)赋予取向性。
倾斜磁场发生线圈14包括X线圈、Y线圈以及Z线圈(省略图示)，设置在呈圆筒形状的静磁场发生线圈12的内周面。
这些X线圈、Y线圈以及Z线圈分别顺次地切换X轴方向、Y轴方向以及Z轴方向，来对孔内的均匀磁场叠加倾斜磁场，以对静磁场赋予强度梯度。在Z线圈激励时，使磁场强度向Z方向倾斜来限定共振面，在施加Z方向的磁场之后Y线圈紧接着加以短时间的倾斜地向检测信号施加与Y坐标成比例的相位调制(相位编码)，接着在数据采集时X线圈加以倾斜地向检测信号施加与X坐标成比例的频率调制(频率编码)。
该叠加的倾斜磁场的切换是通过按照控制时序从倾斜磁场电源24向X线圈、Y线圈以及Z线圈分别输出不同的脉冲信号来实现的。由此，能够确定NMR现象所表现的受验者2的位置，能够提供形成受验者2的图像所需的三维坐标上的位置信息。
在此，如上所述，使用三组正交的倾斜磁场，对这些倾斜磁场分别分配切片方向、相位编码方向以及频率编码方向，利用其组合来能够从各种角度进行拍摄。例如，除了能够拍摄与通过X射线CT装置拍摄的方向相同的方向的横向切片(transverse slice)以外，还能够拍摄与其正交的矢状切片(sagittal slice)和冠状切片(coronal slice)、以及垂直于面的方向不与三组正交的倾斜磁场的轴平行的斜切片(oblique slice)等。
RF照射部16用于基于按照控制时序从信号发送部26发送的高频信号，来向受验者2的关注区域照射RF(Radio Frequency：射频)脉冲。
此外，在图1中，RF照射部16内置于磁场施加机构11，但也可以设置于床18，或者与接收线圈20一体化。
接收线圈20用于检测来自受验者2的响应波(NMR信号)，为了高灵敏度地检测该NMR信号，接近受验者2地配置该接收线圈20。
在此，在接收线圈20中，当NMR信号的电磁波切割其线圈线材时基于电磁感应而产生微弱电流。该微弱电流在信号接收部28中被放大，进一步从模拟信号变换为数字信号后被送至数据处理部32。
即，该机制如下。通过RF照射部16向处于在静磁场上加上Z轴倾斜磁场的状态的受验者2施加共振频率的高频电磁场。磁场的强度满足共振条件的部分的规定的原子核、例如氢原子核被选择性地激励而开始共振。处于符合共振条件的部分(例如，受验者2的规定厚度的断层)的规定的原子核被激励，原子核的自旋轴一齐开始进动。当停止激励脉冲时，在接收线圈20中，本次是正在进动的原子核所放射的电磁波感应出信号，在一段时间内检测出该信号。根据该信号，来观察受验者2的体内的包含规定的原子的组织。然后，为了获知信号的发送位置，加上X和Y的倾斜磁场来探测信号。
图像处理部48基于存储部36中构建的数据，一边重复施加激励信号一边测定检测信号，通过第一次傅立叶变换计算来将共振的频率还原为X坐标，通过第二次傅立叶变换来还原Y坐标，从而得到图像，并在显示部38上显示对应的图像。
例如，通过这种MRI系统，实时地拍摄上述的BOLD信号，通过控制部42对按时间序列拍摄的图像进行如后面说明那样的分析处理，由此能够进行静息态功能连接MRI(rs-fcMRI)的拍摄。
图2是数据处理部32的硬件框图。
作为数据处理部32的硬件，如上所述那样不特别限定，能够使用通用计算机。
在图2中，数据处理部32的计算机主体2010除了存储器驱动器2020、盘驱动器2030以外，还包括：CPU 2040；总线2050，其连接在盘驱动器2030和存储器驱动器2020上；ROM2060，其用于存储启动程序等程序；RAM2070，其用于临时存储应用程序的命令并且提供临时存储空间；非易失性存储装置2080，其用于存储应用程序、系统程序以及数据；以及通信接口2090。通信接口2090相当于用于与驱动部21等进行信号的发送接收的接口部44以及用于经由未图示的网络而与其它计算机进行通信的网络接口50。此外，作为非易失性存储装置2080，能够使用硬盘(HDD)或固态硬盘(SSD：Solid State Drive)等。非易失性存储装置2080相当于存储部36。
CPU 2040通过基于程序而执行的运算处理来实现数据处理部32的各功能、例如控制部42、数据收集部46、图像处理部48的各功能。
使数据处理部32执行上述实施方式的功能的程序也可以被存储在CD-ROM 2200或存储器介质2210中并被插入到盘驱动器2030或存储器驱动器2020并进一步被传输至非易失性存储装置2080。程序在执行时被加载至RAM 2070。
数据处理部32还具备作为输入装置的键盘2100和鼠标2110以及作为输出装置的显示器2120。键盘2100和鼠标2110相当于输入部40，显示器2120相当于显示部38。
用于作为如上所述的数据处理部32而发挥功能的程序也可以未必包括使计算机主体2010执行信息处理装置等的功能的操作系统(OS)。程序只要仅包括用于在受控制的方式下调用适当的功能(模块)从而得到期望的结果的命令的部分即可。数据处理部32如何动作这是众所周知的，因此省略详细的说明。
另外，执行上述程序的计算机既可以是单个也可以是多个。即，可以进行集中处理，或者也可以进行分散处理。
图3是表示通过本实施方式的rs-fcMRI法拍摄的脑的关注区域(ROI：Region ofInterest)的图。
在此，设以针对自闭症谱系(Autistic Spectrum Disorder：ASD)的生物标记物为例，采用93个区域作为关注区域。
作为这种关注区域，例如存在如以下那样的区域。
背内侧前额叶皮质(DMPFC：Dorsomedial Prefrontal Cortex)
腹内侧前额叶皮质(VMPFC：Ventromedial Prefrontal Cortex)
前扣带回皮质(ACC：Anterior Cingulate Cortex)
小脑蚓部(Cerebellar Vermis)，
左侧丘脑(Left Thalamus)，
右下顶叶(Right Inferior Parietal Lobe)，
右侧尾状核(Right Caudate Nucleus)，
右中枕叶(Right Middle Occipital Lobe)，
右中扣带回皮质(Right Middle Cingulate Cortex)，
但是，所采用的脑区并不限定于这样的区域。
例如，也可以根据作为对象的神经/精神疾病来变更所选择的区域。
图4是表示针对如图3所示的关注区域提取表示静息态的功能连接的相关性的相关矩阵(correlation matrix)的过程的概念图。
如图4所示，根据实时地测定出的静息态的fMRI的n个(n：自然数)时刻的fMRI数据，计算各关注区域的平均的“活动度”，计算脑区域间(关注区域间)的活动度的相关值。
在此，作为关注区域，如上述那样考虑93个区域，因此，当考虑对称性时，相关矩阵中的独立的非对角元素为：
(93×93-93)/2＝4278(个)。
此外，在图4中，在图上仅示出34×34的相关性。
图5是说明根据如图4中说明那样的相关矩阵来生成成为生物标记物的判别器的过程的概念图。
如图5所示，数据处理部32根据在健康组(在本例中为114人)、患者组(在本例中为74人)中测定出的静息态功能连接MRI的数据，针对各个受验者通过后面说明那样的过程来导出脑区域间(关注区域间)的活动度的相关矩阵80。
接着，由数据处理部32通过正则化典型相关分析82对包括受验者的疾病/健康标签在内的受验者的属性以及相关矩阵进行特征提取84。“正则化”一般是指如下方法：在机器学习和统计学中，对误差函数追加正则化项来对模型的复杂度和自由度加以抑制，从而防止过学习。此外，在正则化典型相关分析的结果、对于解释变量还同时实现了稀疏化的情况下，特别地称为稀疏典型相关分析(SCCA)。下面，作为具体例，设进行SCCA来进行说明。
并且，由数据处理部32基于正则化典型相关分析的结果，通过利用稀疏逻辑回归的判别分析86来生成判别器88。
此外，如后面说明的那样，健康组和患者组的数据不限于由MRI装置10本身来测定，也可以将在其它MRI装置中测定出的数据也合并来进行判别器的生成。另外，更一般地说，数据处理部32无需一定是用于执行MRI装置的控制的计算机，也可以是专用于接收来自多个MRI装置的测定数据来进行判别器的生成处理并通过所生成的判别器进行判别处理的计算机。
图6是用于说明为了生成成为生物标记物的判别器而由数据处理部32进行的处理的流程图。
下面，参照图6来更详细地说明图5中说明的处理。
基于根据静息态的fMRI数据导出的脑区之间的连接以及受验者的疾病的判别标签来制作生物标记物时的最大问题在于，与数据的数量相比，数据的维数压倒性地多。因此，若不进行正则化而使用数据组来进行用于预测疾病的判别标签(在此，将表示受验者是患有疾病还是健康的标签称为“疾病的判别标签”)的判别器的学习，则判别器过拟合(over-fit)，对于未知数据的预测性能显著下降。
在此，一般在机器学习中，将以用更少数量的解释变量来解释观测数据的方式进行的处理称为“变量选择(或特征提取)”。在本实施方式中，将以下处理称为“提取缩略表达”：从“脑区域间(关注区域间)的活动度的多个相关值(多个连接)”中进行变量选择(特征提取)使得在用于预测作为对象的疾病的判别标签的判别器的机器学习中能够以更少的相关值构成判别器，即，选择更重要的相关值来作为解释变量。
而且，在本实施方式中，作为特征提取法，使用正则化法。这样，在典型相关分析中，执行正则化并且执行稀疏化，来留下更重要的变量来作为解释变量，在以这种方式进行处理的意义上，将这种处理称为“稀疏典型相关分析”。更具体地说，例如，作为用于同时实现稀疏化的正则化法，能够使用如以下说明那样的称为“L1正则化”的、对典型相关分析的参数的绝对值的大小加以惩罚(penalty)的方法。
即，参照图6，当开始判别器生成的处理时(S100)，数据处理部32从存储部36读出各受验者的rs-fcMRI数据(S102)，通过SCCA进行特征提取(S104)。
下面，说明L1正则化典型相关分析。在以下的文献中公开了该L1正则化典型相关分析。
文献：Witten DM,Tibshirani R,and T Hastie.A penalized matrixdecomposition,with applications to sparse principal components and canonicalcorrelation analysis.Biostatistics,Vol.10,No.3,pp.515-534,2009.
首先，在一般的典型相关分析(Canonical Correlation Analysis：CCA)中，考虑如以下那样的数据对x1和x2。其中，设各个变量x1和x2被标准化成具有平均值0、标准偏差1。另外，设数据对x1和x2的数据数量分别为n。
[数1]
关于
在CCA中，求出使z1＝w1x1、z2＝w2x2的相关性最大化的参数即，对如下那样的优化问题进行求解：
在的条件下，
与此相对，当导入L1正则化时，对如下那样的优化问题进行求解。
[数2]
对如下那样的优化问题进行求解：
在||w1||2≤1，||w2||2≤1，||w1||1≤c1，||w2||1≤c2的条件下，
在此，c1和c2是表示L1正则化的强度的参数，能够通过周知的各种方法来根据数据而适当地设定。||w1||1和||w2||1的右下的“1”表示||w1||和||w2||是L1范数。
通过附加用于L1正则化的约束条件，参数w1和w2的元素中的重要度低的元素的值变为0，特征量(解释变量)变得稀疏化。
接着，数据处理部32基于SCCA的结果，通过SLR来进行判别分析(S106)。
SLR是将逻辑回归分析扩展到贝叶斯估计的框架内的方法，是将特征矢量的维数压缩与用于判别的权重估计同时进行的方法。SLR在数据的特征矢量的维数非常高、包含大量不需要的特征量的情况下有用。针对不需要的特征量，将线性判别分析中的权重参数设为零(即进行变量选择)，仅取出与判别相关联的极少数的特征量(稀疏性)。
针对通过SLR得到的特征数据，按分出的每个类别求出该特征数据属于该类别的概率p，将该特征数据分配到输出了最大值的类别。p是通过逻辑回归式来输出的。权重的估计是通过ARD(Automatic Relevance determination：自动相关性确定)来进行的，通过使对类别判别的贡献少的特征量的权重接近0来将该特征量从计算中排除。
即，将使用上述的L1正则化CCA提取出的特征量作为输入，使用基于接下来的分层贝叶斯估计的判别器，来进行疾病/健康标签的预测。
在此，例如，当将通过CCA导出的特征量z1＝w1x1作为向SLR的输入、将S＝{0，1}设为疾病(S＝1)/健康(S＝0)的标签时，如下那样定义SLR的输出为S＝1的概率。
[数3]
将此时的参数矢量w的分布设定为如以下那样的正态分布。此外，下面的式子中的α是表现矢量w的正态分布的方差的超参数矢量(hyper parameter vector)。
[数4]
p(w|α)＝N(w|0，diag(α))
并且，如以下那样设定超参数矢量α的分布，由此通过进行分层贝叶斯估计来估计各个参数的分布。
[数5]
在此，符号“Γ”表示伽玛分布，a0、b0是决定超参数矢量α的伽玛分布的参数。矢量α的第i个元素是αi。
此外，在以下的文献中公开了这种SLR。
文献：Okito Yamashita,Masaaki Sato,Taku Yoshioka,Frank Tong,andYukiyasu Kamitani.“Sparse Estimation automatically selects voxels relevantfor the decoding of fMRI activity patterns.”NeuroImage,Vol.42,No.4,pp.1414-1429,2008.
基于这种判别分析的结果，生成对输入的特征量进行疾病/健康标签(疾病的判别标签)的判别的判别器(S108)。用于确定所生成的判别器的信息(与函数形式和参数有关的数据)被保存到存储部36，在之后被输入测试数据时，该信息使用于执行对于测试数据的疾病的判别标签的估计时的判别处理。
即，基于事前的医师的诊断，将受验者分为健康组、患者组。测定这些受验者的脑区域间(关注区域间)的活动度的相关性(连接)。通过对测定结果的机器学习，来生成判别器以判别与不同于这些受验者的新的受验者有关的测试数据是相当于疾病还是相当于健康。该判别器作为精神疾病的生物标记物而发挥功能。此时，由于判别器是通过逻辑回归而生成的，因此作为生物标记物的输出的“疾病的判别标签”也可以包含有疾病的概率(或者健康的概率)。例如，作为“有疾病的概率为○○％”这样的显示而输出。能够将这种概率用作“诊断标记物”。
另外，关于成为判别器的输出的属性，不一定限于疾病的判别，也可以是与其它属性有关的输出。在该情况下也是，既可以输出该属性属于哪一个类别这样的离散性的判别结果，也可以输出该属性属于某一个类别的连续性的值、例如概率。
即，在生物标记物的学习(制作)中，为了制作精神疾病的生物标记物，将静息态功能性连接MRI(rs-fcMRI)数据作为输入，使用上述的L1正则化CA来进行特征提取，将提取出的特征量作为输入，通过SLR来进行疾病/健康标签判别。
此外，在以上的说明中，设将通过CCA导出的特征量z1＝w1x1用作成为向SLR的输入的特征量来进行了说明。
但是，作为成为向SLR的输入的特征量，并不限定于这种情况。
图7～图9是表示通过SCCA得到的特征量(通过SCCA得到的特征空间)以及以此为输入通过SLR来生成判别器的过程的概念的图。
首先，如图7所示，作为第一方法，将健康/疾病的标签用作受验者的属性，在该标签与基于rs-fcMRI的相关矩阵的元素之间执行SCCA。
将由此得到的特征量(中间表达)z1＝w1x1(与二维典型变量对应)作为SLR的输入，来生成判别器(以下，将该过程称为“方式A”)。
在此，“二维典型变量”是指有疾病(1,0)这样的二维矢量、健康(0,1)这样的二维矢量。
在该情况下，通过SCCA，选择相关矩阵的4278个非对角元素中的例如约700个元素(连接)。
即，使用疾病/健康的诊断结果的某个数据，预先利用SCCA、SLR等机器学习法来决定判别器的参数。通过将所决定出的该参数泛化到未知的rs-fcMRI数据来使其作为“生物标记物”而发挥功能。作为判别器的输入，提供rs-fcMRI中的所选择的连接的线性和，判别器的输出对应疾病标签(1,0)、健康标签(0,1)。
另一方面，如图8所示，作为第二方法，也可以如以下那样。作为受验者的属性，使用疾病/健康标签、受验者的年龄、受验者的性别、进行了测定的机关(测定装置的设置场所，称为测定地点。在后述的例子中测定地点的数量为3。)、fMRI摄像装置的施加磁场强度(相当于摄像装置的分辨率等性能的指标之一)、睁眼/闭眼等拍摄时的实验条件。在这些属性与基于rs-fcMRI的相关矩阵的元素之间执行SCCA。此外，作为fMRI摄像装置的性能的指标，也可以不仅使用施加磁场强度，还使用其它指标。
也就是说，作为典型相关分析中的“解释变量”，取基于rs-fcMRI的相关矩阵的独立元素，作为“基准变量”，取如上所述的受验者的属性。
将由此得到的特征量(中间表达)z1＝w1x1(在本例中，与11维典型变量对应)作为SLR的输入，来生成判别器(以下，将该过程称为“方式B”)。
此外，在此，11维典型变量是指疾病/健康(二维)、3个地点(三维：(1,0,0)、(0,1,0)、(0,0,1))、年龄(一维)、性别(二维：(1,0)、(0,1))、性能(一维)、睁眼闭眼(二维：(1,0)、(0,1))共计11维。
并且，在图9中示出第三方法。在第三方法中，作为受验者的属性，使用疾病/健康标签、受验者的年龄、受验者的性别、进行了测定的机关(测定地点)、fMRI摄像装置的性能(例：施加磁场强度)、睁眼/闭眼等拍摄时的实验条件，在该属性与基于rs-fcMRI的相关矩阵的独立元素之间执行SCCA。该过程与方式B相同。
作为SLR的输入，使用作为SCCA的结果的稀疏地选择出的rs-fcMRI的非对角元素，来生成判别器(以下，称为“方式C”)。
在该情况下，设在SCCA中仅选择与疾病/健康标签有关的rs-fcMRI的相关矩阵的非对角元素。由此，过滤掉表现其它因素的rs-fcMRI数据，实现不依赖于获取数据的设施的生物标记物。
通过取入不同的地点、实验条件的影响等来进行SCCA，能够获知哪个连接与哪个因素有关系，排除与疾病/健康标签没有关联的连接，导出仅与疾病/健康标签相关联的连接。其结果，有利于作为生物标记物的“泛化(generalization)”。
并且，通过利用SLR进行的稀疏化，以这种方式C来例如将通过SCCA稀疏化后得到的约700个连接中的21个连接使用于判别器。
此外，在以上所说明的方式B、方式C中，受验者的属性例如也可以包括表示是否给予了规定的药物的标签或者该药物的用量或给予期间等信息。
图10是表示这种生物标记物的生成过程的概念的图。
在地点1～地点3这三个测定机关中，使用各不相同的测定性能的MRI装置，使用健康组114人、患者组74人的rs-fcMRI数据，如上所述，作为相关矩阵存在4000多个连接，在这样的状态下，通过SCCA处理和SLR处理来生成判别器。
(生物标记物的验证)
图11是表示如以上那样生成的生物标记物的验证过程的概念图。
使用基于在日本国内的地点1～地点3得到的数据并通过学习而得到的用于特征提取和判别的参数，对未用在该参数学习中的数据进行疾病/健康标签的预测，评价与实际的疾病的判别标签的一致性。
图12是表示这种生物标记物的特性的图。
上述的一致性评价的结果是，例如在使用了上述的方式C的情况下，对于日本国内的三个不同的地点的数据示出了最高79％以上的判别性能。另外，在使用利用了相同的参数的生物标记物、并将在美国的六个地点处测量出的rs-fcMRI数据用作输入的情况下，示出了66％以上的判别性能。
在此，在图12中，地点4是将美国的六个地点汇总表示的地点。
此外，在图12中，“灵敏度”是对于疾病组的受验者的检查正确地判定为阳性(有疾病)的概率，“特异度”是对于正常组的受验者的检查正确地判定为阴性(无疾病)的概率。“阳性似然比”是真阳性概率与假阳性概率之比，“阴性似然比”是假阴性概率与真阴性概率之比。DOR是灵敏度与特异度之比。
在图12中，对美国地点的泛化在超越了人种和国别的诊断基准这个意义上示出了通过基于有限的受验者数据的学习而得到的生物标记物的性能之高。
(多疾病生物标记物)
此外，在以上的说明中，在SCCA中的基准变量中，疾病/健康标签与一个疾病有关。
但是，若在以rs-fcMRI数据为解释变量时将多个疾病的标签用作疾病标签，则还能够作为多疾病生物标记物来利用。
即，在该情况下，疾病标签包含用于针对多个疾病各自区分疾病和健康的标签。
图13是表示这种多疾病生物标记物的概念的图。
图13所示的只不过是假想的例示，通过利用上述的针对rs-fcMRI数据的(SCCA+SLR)法来研究到目前为止被称为不同病名的疾病以及对各个疾病下的处方药物等的影响与脑活动之间的关系，能够提取各自的相互关系。
例如，作为上述的疾病标签，包含“自闭症”的疾病/健康标签、“精神分裂症”的疾病/健康标签、“抑郁症”的疾病/健康标签、“强迫性障碍”的疾病/健康标签等。
作为药物标签，除了如“精神药物1”“精神药物2”那样药物的种类以外，还包含各药物的用量及给予期间等信息。
作为摄像条件，除了如上所述的fMRI摄像装置的性能(例：施加磁场强度)以外，还包含是在睁眼状态还是在闭眼状态下进行了测定这样的信息。
并且，作为受验者的个别属性，包含性别和年龄等信息。
此外，在以上的说明中，设作为用于通过脑功能成像法按时间序列测量脑活动的脑活动检测装置使用实时fMRI来进行了说明。但是，作为脑活动检测装置，能够使用上述的fMRI、脑磁图仪、近红外光测量装置(NIRS)、脑电图仪或它们的组合。例如，在使用它们的组合的情况下，fMRI和NIRS检测与脑内的血流变化相关联的信号，具有高空间分辨率。另一方面，脑磁图仪和脑电图仪具有以下特征：具有高时间分辨率，用于检测伴随脑活动的电磁场的变化。因而，例如若将fMRI与脑磁图仪相组合，则能够在空间上和时间上均以高分辨率测量脑活动。若将NIRS与脑电图仪组合，则还能够将也同样地在空间上和时间上均以高分辨率测量脑活动的系统构成为小型、能够携带的大小。
通过如以上那样的结构，能够实现对于神经/精神疾病作为利用脑功能成像法的生物标记物而发挥功能的脑活动分析装置和脑活动分析方法。
另外，在以上的说明中，说明了在作为受验者的属性包括“疾病的判别标签”的情况下通过利用机器学习生成判别器来使该判别器作为生物标记物而发挥功能的例子，但是本发明并不一定限定于这种情况。只要能够事前通过客观的方法将得到作为机器学习的对象的测定结果的对象受验者组分为多个类别，测定受验者的脑区域间(关注区域间)的活动度的相关性(连接)，通过对测定结果的机器学习来生成对于类别的判别器，则也可以用于其它判别。
另外，如上所述，这种判别也可以将属于某个属性的可能性显示为概率。
因而，例如，能够客观地评价采取某个“训练”或“行为模式”是否有助于增进受验者的健康。另外，即使处于实际上未达到疾病的状态(“未病”)，也能够客观地评价某个如“食物”“饮料”等那样的摄取物或者某个活动等是否在更接近健康状态的方面有效。
另外，即使在未病的状态下，只要如上所述那样例如输出“健康的概率为○○％”这样的显示，就能够对用户显示关于健康状态的客观数值。此时，所输出的不一定是概率，也可以将“健康度的连续值、例如健康的概率”变换为得分来进行显示。通过进行这种显示，本实施方式的装置除了能够用于诊断的辅助以外，还能够用作用于用户的健康管理的装置。
本次公开的实施方式是用于具体实施本发明的结构的例示，并不限制本发明的技术范围。本发明的技术范围不是由实施方式的说明来表示，而是由权利要求书来表示，意图包括权利要求书的语句上的范围及等同的含义的范围内的变更。
附图标记说明
2：受验者；6：显示器；10：MRI装置；11：磁场施加机构；12：静磁场发生线圈；14：倾斜磁场发生线圈；16：RF照射部；18：床；20：接收线圈；21：驱动部；22：静磁场电源；24：倾斜磁场电源；26：信号发送部；28：信号接收部；30：床驱动部；32：数据处理部；36：存储部；38：显示部；40：输入部；42：控制部；44：接口部；46：数据收集部；48：图像处理部；50：网络接口。

Claims (18)

1.一种脑活动分析装置，具备：
判别器生成单元，其用于根据通过脑活动探测装置按时间序列预先测定表示多个受验者各自的脑内的多个规定区域中的静息态的脑活动的信号而得到的信号来生成判别器，上述判别器生成单元针对上述多个规定区域中的脑活动之间的多个相关性以及包括受验者个人的多种属性和上述脑活动探测装置中的测定条件的属性信息，来提取在上述多个受验者上共通的缩略表达，对于提取出的缩略表达，生成针对上述受验者个人的上述多种属性中的特定属性的判别器；
存储装置，其用于存储用于确定上述判别器的信息；以及
运算装置，其中，上述运算装置构成为：
基于根据上述存储装置中存储的信息而确定的判别器，来对输入数据进行判别处理。
2.根据权利要求1所述的脑活动分析装置，其特征在于，
上述特定属性是疾病的判别标签。
3.根据权利要求2所述的脑活动分析装置，其特征在于，
上述受验者个人的属性包括向该受验者给予的药物的标签。
4.根据权利要求1～3中的任一项所述的脑活动分析装置，其特征在于，
上述脑活动探测装置包括分别设置于多个不同的场所的多个脑活动测定装置，
上述脑活动探测装置中的测定条件包括上述场所的信息，
上述判别器生成单元包括提取单元，该提取单元从上述多个规定区域中的脑活动之间的多个相关性中通过变量选择来提取在上述多个脑活动测定装置的测定条件和上述多个受验者个人的属性上共通的上述缩略表达。
5.根据权利要求1～3中的任一项所述的脑活动分析装置，其特征在于，
上述判别器生成单元包括回归单元，该回归单元通过对提取出的上述缩略表达进一步进行变量选择的回归来生成上述判别器。
6.根据权利要求4所述的脑活动分析装置，其特征在于，
上述判别器生成单元包括回归单元，该回归单元通过对提取出的上述缩略表达进一步进行变量选择的回归来生成上述判别器。
7.根据权利要求4所述的脑活动分析装置，其特征在于，
上述提取单元包括相关分析单元，该相关分析单元用于根据由上述脑活动探测装置探测出的信号来计算上述多个规定区域间的活动的相关矩阵，在上述属性信息与上述相关矩阵的非对角元素之间执行正则化典型相关分析来提取上述缩略表达。
8.根据权利要求1～3中的任一项所述的脑活动分析装置，其特征在于，
上述脑活动探测装置包括多个脑活动测定装置，
上述判别器生成单元包括：
提取单元，其从上述多个规定区域中的脑活动之间的多个相关性中通过变量选择来提取在上述多个脑活动测定装置的测定条件和上述多个受验者个人的属性上共通的上述缩略表达；以及
回归单元，其通过对提取出的上述缩略表达进一步进行变量选择的回归来生成上述判别器，
其中，上述提取单元包括相关分析单元，该相关分析单元用于根据由上述脑活动探测装置探测出的信号来计算上述多个规定区域间的活动的相关矩阵，在上述属性信息与上述相关矩阵的非对角元素之间执行正则化典型相关分析来提取上述缩略表达，
上述回归单元包括回归分析单元，该回归分析单元通过对上述正则化典型相关分析的结果和上述受验者个人的属性进行稀疏逻辑回归来生成判别器。
9.根据权利要求4所述的脑活动分析装置，其特征在于，
上述多个脑活动测定装置是用于通过脑功能成像法来按时间序列测量脑活动的装置。
10.根据权利要求2所述的脑活动分析装置，其特征在于，
上述判别处理是疾病的判别标签的判别，该疾病的判别标签表示与上述输入数据对应的受验者患有还是不患有神经/精神疾病。
11.根据权利要求4所述的脑活动分析装置，其特征在于，
上述属性信息包括表示针对上述受验者个人的患有还是不患有神经/精神疾病的疾病的判别标签、表示上述受验者的个人特性的标签、以及对利用上述脑活动探测装置进行的测定赋予特征的测定条件的信息，
上述判别处理是疾病的判别标签的判别，该疾病的判别标签表示与上述输入数据对应的受验者患有还是不患有神经/精神疾病。
12.根据权利要求7所述的脑活动分析装置，其特征在于，
上述正则化典型相关分析是基于L1正则化的典型相关分析。
13.根据权利要求9所述的脑活动分析装置，其特征在于，
上述脑活动测定装置拍摄静息态功能连接核磁共振影像。
14.一种脑活动分析方法，用于使具备运算装置和存储装置的计算机进行脑活动的分析，该脑活动分析方法具备以下步骤：
根据通过脑活动探测装置按时间序列预先测定表示多个受验者各自的脑内的多个规定区域中的静息态的脑活动的信号而得到的信号来生成判别器的步骤，上述生成判别器的步骤包括以下步骤：针对上述多个规定区域中的脑活动之间的多个相关性以及包括受验者个人的多种属性和上述脑活动探测装置中的测定条件的属性信息，来提取在上述多个受验者上共通的缩略表达，对于提取出的缩略表达，生成针对上述受验者个人的上述多种属性中的特定属性的判别器；
将用于确定上述判别器的信息存储在上述存储装置中的步骤；以及
基于根据上述存储装置中存储的信息而确定的判别器来对输入数据进行判别处理的步骤。
15.根据权利要求14所述的脑活动分析方法，其特征在于，
上述脑活动探测装置包括分别设置于多个不同的场所的多个脑活动测定装置，
上述脑活动探测装置中的测定条件包括上述场所的信息，
上述生成判别器的步骤包括以下步骤：
从上述多个规定区域中的脑活动之间的多个相关性中通过变量选择来提取在上述多个脑活动测定装置的测定条件和上述多个受验者个人的属性上共通的上述缩略表达。
16.根据权利要求14或15所述的脑活动分析方法，其特征在于，
上述生成判别器的步骤包括以下步骤：
通过对提取出的上述缩略表达进一步进行变量选择的回归来生成上述判别器。
17.一种生物标记物装置，通过计算机对脑活动的分析来生成作为生物标记物的输出，该生物标记物装置具备：
存储装置，其用于存储用于确定判别器的信息，
上述判别器是根据通过脑活动探测装置按时间序列预先测定表示多个受验者各自的脑内的多个规定区域中的静息态的脑活动的信号而得到的信号来计算出的判别器，其中，针对上述多个规定区域中的脑活动之间的多个相关性以及包括受验者个人的多种属性和上述脑活动探测装置中的测定条件的属性信息，来提取在上述多个受验者上共通的缩略表达，对于提取出的缩略表达，针对上述受验者个人的上述多种属性中的疾病的判别标签来生成判别器；以及
运算装置，其基于根据上述存储装置中存储的信息而确定的判别器，来对输入数据进行判别处理。
18.根据权利要求17所述的生物标记物装置，其特征在于，
上述受验者个人的属性包括向该受验者给予的药物的标签。
CN201480024893.1A 2013-05-01 2014-04-24 脑活动分析装置、脑活动分析方法以及生物标记物装置 Active CN105163659B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2013-096493 2013-05-01
JP2013096493 2013-05-01
PCT/JP2014/061544 WO2014178323A1 (ja) 2013-05-01 2014-04-24 脳活動解析装置、脳活動解析方法およびバイオマーカー装置

Publications (2)

Publication Number Publication Date
CN105163659A CN105163659A (zh) 2015-12-16
CN105163659B true CN105163659B (zh) 2018-10-23

Family

ID=51843455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480024893.1A Active CN105163659B (zh) 2013-05-01 2014-04-24 脑活动分析装置、脑活动分析方法以及生物标记物装置

Country Status (5)

Country Link
US (2) US10357181B2 (zh)
EP (1) EP2992824B1 (zh)
JP (3) JP5691086B1 (zh)
CN (1) CN105163659B (zh)
WO (1) WO2014178323A1 (zh)

Families Citing this family (42)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9275457B1 (en) 2014-08-28 2016-03-01 International Business Machines Corporation Real-time subject-driven functional connectivity analysis
US20160070879A1 (en) * 2014-09-09 2016-03-10 Lockheed Martin Corporation Method and apparatus for disease detection
US11630512B2 (en) 2015-05-28 2023-04-18 Riken Determination device, determination method, program, and information storage medium
JP6566471B2 (ja) * 2015-07-01 2019-08-28 国立大学法人福井大学 自閉スペクトラム症診断補助のための医療用画像処理方法及び医療用画像処理システム並びにバイオマーカー
JP6590411B2 (ja) * 2015-10-16 2019-10-16 国立大学法人広島大学 感性評価方法
US11382556B2 (en) 2015-11-24 2022-07-12 Advanced Telecommunications Research Institute International Brain activity analyzing apparatus, brain activity analyzing method, program and biomarker apparatus
WO2017109860A1 (ja) * 2015-12-22 2017-06-29 株式会社ニコン 画像処理装置
CN105596004B (zh) * 2015-12-28 2017-08-25 中国人民解放军国防科学技术大学 基于成组典型相关分析的脑功能磁共振成像盲源分离方法
CN108697392B (zh) * 2016-02-29 2022-01-28 大金工业株式会社 判定结果输出装置、判定结果提供装置及判定结果输出系统
KR101796055B1 (ko) * 2016-06-02 2017-11-10 고려대학교 산학협력단 다중 뇌 연결망 구축을 통한 뇌 상태 모니터링 방법 및 장치
WO2018005820A1 (en) * 2016-06-29 2018-01-04 The University Of North Carolina At Chapel Hill Methods, systems, and computer readable media for utlizing brain structural characteristics for predicting a diagnosis of a neurobehavioral disorder
WO2018005814A1 (en) * 2016-06-29 2018-01-04 The University Of North Carolina At Chapel Hill Methods, systems, and computer readable media for utilizing functional connectivity brain imaging for diagnosis of a neurobehavioral disorder
CN106264460B (zh) * 2016-07-29 2019-11-19 北京医拍智能科技有限公司 基于自学习的脑活动多维时间序列信号的解码方法及装置
EP3512415A4 (en) * 2016-09-13 2020-08-19 Ohio State Innovation Foundation NEURAL ARCHITECTURE MODELING SYSTEMS AND METHODS
JP6821171B2 (ja) * 2016-09-20 2021-01-27 学校法人同志社 脳活動状態定量化方法および脳活動状態計測装置
US10034645B1 (en) * 2017-04-13 2018-07-31 The Board Of Trustees Of The Leland Stanford Junior University Systems and methods for detecting complex networks in MRI image data
EP3684463A4 (en) 2017-09-19 2021-06-23 Neuroenhancement Lab, LLC NEURO-ACTIVATION PROCESS AND APPARATUS
JP7057995B2 (ja) * 2017-10-03 2022-04-21 株式会社国際電気通信基礎技術研究所 哺乳動物におけるうつ病の病態プロファイルを検出するバイオマーカーとその利用
WO2019069955A1 (ja) * 2017-10-03 2019-04-11 株式会社国際電気通信基礎技術研究所 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
CN107680677B (zh) * 2017-10-11 2020-09-15 四川大学 基于脑网络分析的神经精神疾病分类方法
US11717686B2 (en) 2017-12-04 2023-08-08 Neuroenhancement Lab, LLC Method and apparatus for neuroenhancement to facilitate learning and performance
JP7069716B2 (ja) * 2017-12-28 2022-05-18 株式会社リコー 生体機能計測解析システム、生体機能計測解析プログラム及び生体機能計測解析方法
WO2019133997A1 (en) 2017-12-31 2019-07-04 Neuroenhancement Lab, LLC System and method for neuroenhancement to enhance emotional response
EP3751473A4 (en) * 2018-02-09 2021-12-08 Axion Research Inc. SYSTEM THAT ESTIMATE THE CONDITION OF A COMPLEX SYSTEM TO INSPECT
EP3763288A4 (en) * 2018-03-09 2021-12-01 Advanced Telecommunications Research Institute International DEVICE, METHOD AND PROGRAM FOR TRAINING BRAIN ACTIVITY
US11364361B2 (en) 2018-04-20 2022-06-21 Neuroenhancement Lab, LLC System and method for inducing sleep by transplanting mental states
CN113382683A (zh) 2018-09-14 2021-09-10 纽罗因恒思蒙特实验有限责任公司 改善睡眠的系统和方法
JP6812022B2 (ja) * 2018-10-11 2021-01-13 株式会社国際電気通信基礎技術研究所 脳機能結合相関値の調整方法、脳機能結合相関値の調整システム、脳活動分類器のハーモナイズ方法、脳活動分類器のハーモナイズシステム、および脳活動バイオマーカシステム
WO2020075737A1 (ja) 2018-10-11 2020-04-16 株式会社国際電気通信基礎技術研究所 脳機能結合相関値の調整方法、脳機能結合相関値の調整システム、脳活動分類器のハーモナイズ方法、脳活動分類器のハーモナイズシステム、及び脳活動バイオマーカシステム
JP7218154B2 (ja) * 2018-11-05 2023-02-06 花王株式会社 評価方法
WO2020100983A1 (ja) * 2018-11-15 2020-05-22 公立大学法人広島市立大学 脳機能計測装置及び脳機能計測方法
JP7147623B2 (ja) * 2019-02-21 2022-10-05 株式会社島津製作所 脳血流量の特徴量の抽出方法
KR102009758B1 (ko) * 2019-03-21 2019-10-21 한국과학기술정보연구원 바이오 시그널 가시화 시스템 및 유효 패턴 추출 방법
US11786694B2 (en) 2019-05-24 2023-10-17 NeuroLight, Inc. Device, method, and app for facilitating sleep
US11443429B2 (en) * 2019-05-30 2022-09-13 Washington University Atlas registration for resting state network mapping in patients with brain tumors
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
JP6999118B2 (ja) * 2020-05-07 2022-01-18 株式会社ニコン 画像処理装置
US11581082B2 (en) * 2020-05-15 2023-02-14 Omniscient Neurotechnology Pty Limited Subsetting brain data
WO2022014682A1 (ja) * 2020-07-17 2022-01-20 株式会社国際電気通信基礎技術研究所 脳機能結合相関値のクラスタリング装置、脳機能結合相関値のクラスタリングシステム、脳機能結合相関値のクラスタリング方法、脳機能結合相関値の分類器プログラムおよび脳活動マーカー分類システム
KR102379132B1 (ko) * 2021-06-30 2022-03-30 액티브레인바이오(주) 디지털 컨텐츠 기반 치료 정보 제공 장치 및 방법
CN113707324A (zh) * 2021-09-17 2021-11-26 四川大学华西医院 一种评估脑功能连接和临床症状之间相关性的分析方法
CN116491960B (zh) * 2023-06-28 2023-09-19 南昌大学第一附属医院 脑瞬态监测设备、电子设备及存储介质

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020103429A1 (en) 2001-01-30 2002-08-01 Decharms R. Christopher Methods for physiological monitoring, training, exercise and regulation
US20040116798A1 (en) * 2002-12-11 2004-06-17 Robert Cancro Method and system for investigation of central nervous system drugs using 3-D brain source localization
JP4025823B2 (ja) * 2004-02-24 2007-12-26 国立精神・神経センター総長 脳疾患の診断支援方法及び装置
US20050215884A1 (en) 2004-02-27 2005-09-29 Greicius Michael D Evaluation of Alzheimer's disease using an independent component analysis of an individual's resting-state functional MRI
JP2006132313A (ja) 2004-10-07 2006-05-25 Yashima:Kk 基礎と架台の接合構造
US7577472B2 (en) 2005-03-18 2009-08-18 The Mcw Research Foundation Inc. MRI method for producing an index indicative of brain disorders
US8386192B2 (en) 2005-06-09 2013-02-26 Hitachi Medical Corporation Disease diagnosis support system
CA2657407A1 (en) * 2006-07-06 2008-01-10 Regents Of The University Of Minnesota Analysis of brain patterns using temporal measures
JP4978860B2 (ja) * 2007-01-24 2012-07-18 株式会社国際電気通信基礎技術研究所 行動予測方法及び行動予測装置
US9480425B2 (en) * 2008-04-17 2016-11-01 Washington University Task-less optical mapping of dynamic brain function using resting state functional connectivity
US20100249573A1 (en) * 2009-03-30 2010-09-30 Marks Donald H Brain function decoding process and system
JP5569714B2 (ja) 2009-06-16 2014-08-13 株式会社国際電気通信基礎技術研究所 磁気共鳴画像化装置および磁気共鳴画像化方法
JP5343790B2 (ja) * 2009-09-28 2013-11-13 株式会社国際電気通信基礎技術研究所 脳活動情報出力装置、脳活動情報出力方法、およびプログラム
WO2011115956A1 (en) 2010-03-15 2011-09-22 Mcw Research Foundation, Inc. Systems and methods for detection and prediction of brain disorders based on neural network interaction
US9612306B2 (en) * 2010-11-26 2017-04-04 Brigham And Women's Hospital, Inc. Method for assessing repetitive head injuries with two-dimensional magnetic resonance spectroscopy
US8451930B2 (en) 2011-02-08 2013-05-28 Mitsubishi Electric Research Laboratories, Inc. Sparse channel estimation using order extension for frequency and time selective MIMO wireless communication networks
JP5959016B2 (ja) 2011-05-31 2016-08-02 国立大学法人 名古屋工業大学 認知機能障害判別装置、認知機能障害判別システム、およびプログラム
US8861815B2 (en) * 2011-08-03 2014-10-14 International Business Machines Corporation Systems and methods for modeling and processing functional magnetic resonance image data using full-brain vector auto-regressive model
JP5320543B2 (ja) 2011-11-08 2013-10-23 株式会社国際電気通信基礎技術研究所 脳機能亢進支援装置および脳機能亢進支援方法
US9480402B2 (en) * 2011-11-11 2016-11-01 Washington University System and method for task-less mapping of brain activity
US9510756B2 (en) * 2012-03-05 2016-12-06 Siemens Healthcare Gmbh Method and system for diagnosis of attention deficit hyperactivity disorder from magnetic resonance images
JP5641531B1 (ja) * 2013-05-01 2014-12-17 株式会社国際電気通信基礎技術研究所 脳活動訓練装置および脳活動訓練システム
US9265441B2 (en) * 2013-07-12 2016-02-23 Siemens Aktiengesellschaft Assessment of traumatic brain injury

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Sparse estimation automatically selects voxels relevant for the decoding of fMRI activity patterns;Okito Yamashita,et al.;《NeuroImage》;20081231;第42卷;第1414-1429页 *

Also Published As

Publication number Publication date
JPWO2014178323A1 (ja) 2017-02-23
JP5691086B1 (ja) 2015-04-01
CN105163659A (zh) 2015-12-16
EP2992824B1 (en) 2021-12-29
JP6497687B2 (ja) 2019-04-10
EP2992824A1 (en) 2016-03-09
US20190298207A1 (en) 2019-10-03
US10357181B2 (en) 2019-07-23
EP2992824A4 (en) 2017-01-11
JP6341513B2 (ja) 2018-06-13
WO2014178323A1 (ja) 2014-11-06
JP2018118144A (ja) 2018-08-02
JP2015062817A (ja) 2015-04-09
US20150272461A1 (en) 2015-10-01

Similar Documents

Publication Publication Date Title
CN105163659B (zh) 脑活动分析装置、脑活动分析方法以及生物标记物装置
CN105188528B (zh) 脑活动训练装置以及脑活动训练系统
JP6722845B2 (ja) 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
JP6903518B2 (ja) 脳活動解析システム、脳活動解析装置、及び脳活動解析プログラム
JP6875054B2 (ja) 脳活動訓練装置、脳活動訓練方法および脳活動訓練プログラム
JP6492356B2 (ja) 脳活動訓練装置および脳活動訓練方法
WO2019069955A1 (ja) 判別装置、うつ症状の判別方法、うつ症状のレベルの判定方法、うつ病患者の層別化方法、うつ症状の治療効果の判定方法及び脳活動訓練装置
JP2020062369A (ja) 脳機能結合相関値の調整方法、脳機能結合相関値の調整システム、脳活動分類器のハーモナイズ方法、脳活動分類器のハーモナイズシステム、および脳活動バイオマーカシステム
JP7057995B2 (ja) 哺乳動物におけるうつ病の病態プロファイルを検出するバイオマーカーとその利用

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant