JP5856364B2 - 抑うつ病の評価方法、抑うつ病評価装置、抑うつ病評価方法、抑うつ病評価システム、抑うつ病評価プログラム、記録媒体および端末装置 - Google Patents

抑うつ病の評価方法、抑うつ病評価装置、抑うつ病評価方法、抑うつ病評価システム、抑うつ病評価プログラム、記録媒体および端末装置 Download PDF

Info

Publication number
JP5856364B2
JP5856364B2 JP2008527784A JP2008527784A JP5856364B2 JP 5856364 B2 JP5856364 B2 JP 5856364B2 JP 2008527784 A JP2008527784 A JP 2008527784A JP 2008527784 A JP2008527784 A JP 2008527784A JP 5856364 B2 JP5856364 B2 JP 5856364B2
Authority
JP
Japan
Prior art keywords
depression
discriminant
evaluation
value
amino acid
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
JP2008527784A
Other languages
English (en)
Other versions
JPWO2008016101A1 (ja
Inventor
田中 孝幸
孝幸 田中
幸夫 二瓶
幸夫 二瓶
敏彦 安東
敏彦 安東
實 山門
實 山門
慶太 高松
慶太 高松
英孝 高橋
英孝 高橋
典雄 笹森
典雄 笹森
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ajinomoto Co Inc
Original Assignee
Ajinomoto Co Inc
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 Ajinomoto Co Inc filed Critical Ajinomoto Co Inc
Priority to JP2008527784A priority Critical patent/JP5856364B2/ja
Publication of JPWO2008016101A1 publication Critical patent/JPWO2008016101A1/ja
Application granted granted Critical
Publication of JP5856364B2 publication Critical patent/JP5856364B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/68Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
    • G01N33/6803General methods of protein analysis not limited to specific proteins or families of proteins
    • G01N33/6806Determination of free amino acids
    • G01N33/6812Assays for specific amino acids
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • G16B40/20Supervised data analysis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • G16B40/30Unsupervised data analysis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B50/00ICT programming tools or database systems specially adapted for bioinformatics
    • G16B50/30Data warehousing; Computing architectures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/30Psychoses; Psychiatry
    • G01N2800/304Mood disorders, e.g. bipolar, depression
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B50/00ICT programming tools or database systems specially adapted for bioinformatics

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Biophysics (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Biotechnology (AREA)
  • Data Mining & Analysis (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Evolutionary Biology (AREA)
  • Molecular Biology (AREA)
  • Bioethics (AREA)
  • Databases & Information Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Software Systems (AREA)
  • Hematology (AREA)
  • Urology & Nephrology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Epidemiology (AREA)
  • Evolutionary Computation (AREA)
  • Public Health (AREA)
  • Immunology (AREA)
  • Biomedical Technology (AREA)
  • Analytical Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Cell Biology (AREA)
  • Microbiology (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Genetics & Genomics (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Description

本発明は、血液(血漿)中のアミノ酸濃度を利用したストレスの評価方法、ならびにストレス評価装置、ストレス評価方法、ストレス評価システム、ストレス評価プログラムおよび記録媒体に関するものである。
環境の変化や人間関係、過労などに心身の適応が困難な場合、ストレスを引き起こし、生活習慣の偏りや摂食障害など様々な症状が表れることが知られている。このような症状は、肥満、高血圧、高脂血症、糖尿病といった生活習慣病の源流にあたると考えられており、ストレスに、いかに早く気づき、ストレスに対する対応(コーピング)を行っていくのかは、生活習慣病を未然に防ぐためにも、重要な課題である。
また、従来、ストレスが原因で起こる精神疾患:ストレス所見の診断には、問診票が広く用いられており、世界保健機構(WHO)が作成したICD−10(国際疾病分類第10改訂版)や米国精神医学会(APA)が作成したDSM−IV(精神疾患の診断・統計マニュアル)などがある。
また、ストレス所見の客観的計測手法としては脳波や心拍(特許文献1参照)を測定する電気生理的手法やストレスタンパク質(特許文献2参照、特許文献3参照)、各種遺伝子発現(特許文献4、特許文献5、特許文献6、特許文献7、特許文献8参照)を測定する生化学的方法がある。
加えて、アミノ酸を摂取することによってストレスに対する抵抗性が高まること(特許文献9参照)やストレス下において血漿中のアミノ酸濃度が変動していること知られている。
Altamura Cら(非特許文献1参照)は、大うつ病患者の血漿中にて、タウリン、セリン/グリシンの濃度が増加し、グリシンの濃度が減少すること、Mitani H(非特許文献2参照)らは、タウリン、グリシン、グルタミン、グルタミン酸の濃度が増加すること、Mauri MC(非特許文献3、非特許文献4参照)らは、リジン、タウリン、グルタミン酸、チロシン/LNAAs(バリン、ロイシン、イソロイシン、チロシン、フェニルアラニンの総和)の濃度が増加し、トリプトファン/LNAAsの濃度が減少すること、岸本(非特許文献5参照)は、トリプトファン、チロシン、リジン、ロイシン、イソロイシン、フェニルアラニン、メチオニン、アラニン、グルタミン酸、アルギニン、ヒスチジンの濃度が減少することを示している。
加えて、Altamura Cら(非特許文献1参照)は、グリシン、グルタミン酸、タウリンを用いた線形判別式によって大うつ病を判別可能であること、Mauri MC(非特許文献3、非特許文献4参照)らは、チロシン/LNAAs、トリプトファン/LNAAsが、抗うつ剤の薬効指標となり得ることを示している。
特許第3070346号 特許第3429043号 特許第3445038号 特開2003−274949号公報 特開2004−208547号公報 特開2005−143420号公報 特開2005−312435号公報 特開2006−15号公報 国際公開第2005/07408号パンフレット Altamura C, Maes M, Dai J, Meltzer HY(1995) Eur Neuropsychopharmacol.5. Suppl71−5. Mitani H, Shirayama Y, Yamada T, Maeda K, Ashby CR Jr, Kawahara R (2006) Prog Neuropsychopharmacol Biol Psychiatry. 30(6).1155−8. Mauri MC, Ferrara A, Boscati L, Bravin S, Zamberlan F, Alecci M, Invernizzi G(1998) Neuropsychobiology. 37(3).124−9. Mauri MC, Boscati L, Volonteri LS, Scalvini ME, Steinhilber CP, Laini V, Zamberlan F (2001) Neuropsychobiology.44(3).134−8. 岸本英爾(1993) 脳と精神の医学第4巻第2号1993年4月
しかしながら、ストレス所見診断に用いられる問診票の結果には、被験者の主観が入ってしまうことが多く、このことが問診結果の精度を下げる要因となっているという問題があった。
そこで、ストレス所見をより精度の高い診断で抽出するためには、被験者の主観が入ることなく、ストレス所見を客観的に診断できる指標が望まれているのだが、上記ストレス所見の客観的診断手法のうち、脳波や心拍、ストレスタンパク質、各種遺伝子発現を測定する方法は、診断指標の測定や測定結果の解釈が困難であり、測定機器の汎用性が低いことから、普及の目処がたっていないという問題点があった。
本発明は、上記問題点に鑑みてなされたものであって、血液中のアミノ酸濃度のうちストレスの状態と関連するアミノ酸濃度を利用してストレスの状態を精度よく評価することができるストレスの評価方法、ならびにストレス評価装置、ストレス評価方法、ストレス評価システム、ストレス評価プログラムおよび記録媒体を提供することを目的とする。
本発明者らは、上記課題を解決するために鋭意検討し、血液中のアミノ酸濃度によるストレス負荷状態間の判別に有用なアミノ酸変数の同定、ならびにアミノ酸変数を用いる相関式(指標式)とストレス所見を対象とした精神疾患、とりわけ、うつの病態進行に有意な相関があることを見出し、本発明を完成するに至った。具体的には、健常者の血漿中アミノ酸濃度とストレス所見者の血漿中アミノ酸濃度の差異に注目し、血漿中アミノ酸濃度測定によってストレス所見を客観的に高精度で診断する指標を開発した。血漿中アミノ酸であれば、安価・高速・簡便に測定することができ、診断指標として普及する可能性が高く、上記ストレス所見の客観的計測手法の課題を解決することができる。本発明は、以下を包含する。
すなわち、上述した課題を解決し、目的を達成するために、本発明にかかるストレスの評価方法は、評価対象から採取した血液から、アミノ酸の濃度値に関するアミノ酸濃度データを測定する測定ステップと、前記測定ステップで測定した前記評価対象の前記アミノ酸濃度データに含まれるLys,His,ABA(ABAはα−アミノ酪酸を表す。以下同様。),Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの前記濃度値に基づいて、前記評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価する濃度値基準評価ステップとを実行することを特徴とする。
また、本発明にかかるストレスの評価方法は、前記のストレスの評価方法において、前記濃度値基準評価ステップは、前記測定ステップで測定した前記評価対象の前記アミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの前記濃度値に基づいて、前記評価対象につき、前記抑うつ病または非抑うつ病であるか否かを判別、または前記大うつ病または非大うつ病であるか否かを判別する濃度値基準判別ステップをさらに含むことを特徴とする。
また、本発明にかかるストレスの評価方法は、前記のストレスの評価方法において、前記濃度値基準評価ステップは、前記測定ステップで測定した前記評価対象の前記アミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの前記濃度値および前記アミノ酸の濃度を変数とする予め設定した多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出ステップと、前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記ストレスの状態を評価する判別値基準評価ステップとをさらに含み、前記多変量判別式は、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むことを特徴とする。
また、本発明にかかるストレスの評価方法は、前記のストレスの評価方法において、前記判別値基準評価ステップは、前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病または非抑うつ病であるか否かを判別、または前記大うつ病または非大うつ病であるか否かを判別する判別値基準判別ステップをさらに含み、前記多変量判別式は、1つの分数式または複数の前記分数式の和で表され、それを構成する前記分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むこと、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つであることを特徴とする。
また、本発明にかかるストレスの評価方法は、前記のストレスの評価方法において、前記判別値基準判別ステップで前記抑うつ病または前記非抑うつ病であるか否かを判別する場合、前記多変量判別式は、数式1または数式2であること、またはSer,Glu,Cys,Lysを前記変数とする前記ロジスティック回帰式またはSer,Glu,Cys,Lysを前記変数とする前記線形判別式であることを特徴とする。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレスの評価方法は、前記のストレスの評価方法において、前記判別値基準判別ステップで前記大うつ病または前記非大うつ病であるか否かを判別する場合、前記多変量判別式は、数式3、数式4、数式5、数式6または数式7であること、またはSer,Gln,Pro,ABA,Tyr,Trpを前記変数とする前記ロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを前記変数とする前記ロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを前記変数とする前記線形判別式またはAsn,Gln,Tau,Cit,ABAを前記変数とする前記線形判別式であることを特徴とする。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
また、本発明はストレス評価装置に関するものであり、本発明にかかるストレス評価装置は、制御手段と記憶手段とを備え、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価するストレス評価装置であって、前記制御手段は、アミノ酸の濃度値に関する予め取得した前記評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの前記濃度値および前記アミノ酸の濃度を変数とする前記記憶手段で記憶した多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出手段と、前記判別値算出手段で算出した前記判別値に基づいて、前記評価対象につき、前記ストレスの状態を評価する判別値基準評価手段とを備え、前記多変量判別式は、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むことを特徴とする。
また、本発明にかかるストレス評価装置は、前記のストレス評価装置において、前記判別値基準評価手段は、前記判別値算出手段で算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病または非抑うつ病であるか否かを判別、または前記大うつ病または非大うつ病であるか否かを判別する判別値基準判別手段をさらに備え、前記多変量判別式は、1つの分数式または複数の前記分数式の和で表され、それを構成する前記分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むこと、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つであることを特徴とする。
また、本発明にかかるストレス評価装置は、前記のストレス評価装置において、前記判別値基準判別手段で前記抑うつ病または前記非抑うつ病であるか否かを判別する場合、前記多変量判別式は、数式1または数式2であること、またはSer,Glu,Cys,Lysを前記変数とする前記ロジスティック回帰式またはSer,Glu,Cys,Lysを前記変数とする前記線形判別式であることを特徴とする。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価装置は、前記のストレス評価装置において、前記判別値基準判別手段で前記大うつ病または前記非大うつ病であるか否かを判別する場合、前記多変量判別式は、数式3、数式4、数式5、数式6または数式7であること、またはSer,Gln,Pro,ABA,Tyr,Trpを前記変数とする前記ロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを前記変数とする前記ロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを前記変数とする前記線形判別式またはAsn,Gln,Tau,Cit,ABAを前記変数とする前記線形判別式であることを特徴とする。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価装置は、前記のストレス評価装置において、前記制御手段は、前記アミノ酸濃度データと前記ストレスの状態を表す指標に関するストレス状態指標データとを含む前記記憶手段で記憶したストレス状態情報に基づいて、前記記憶手段で記憶する前記多変量判別式を作成する多変量判別式作成手段をさらに備え、前記多変量判別式作成手段は、前記ストレス状態情報から所定の式作成手法に基づいて、前記多変量判別式の候補である候補多変量判別式を作成する候補多変量判別式作成手段と、前記候補多変量判別式作成手段で作成した前記候補多変量判別式を、所定の検証手法に基づいて検証する候補多変量判別式検証手段と、前記候補多変量判別式検証手段での検証結果から所定の変数選択手法に基づいて前記候補多変量判別式の変数を選択することで、前記候補多変量判別式を作成する際に用いる前記ストレス状態情報に含まれる前記アミノ酸濃度データの組み合わせを選択する変数選択手段と、をさらに備え、前記候補多変量判別式作成手段、前記候補多変量判別式検証手段および前記変数選択手段を繰り返し実行して蓄積した前記検証結果に基づいて、複数の前記候補多変量判別式の中から前記多変量判別式として採用する前記候補多変量判別式を選出することで、前記多変量判別式を作成することを特徴とする。
また、本発明はストレス評価方法に関するものであり、本発明にかかるストレス評価方法は、制御手段と記憶手段とを備えた情報処理装置で実行する、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価するストレス評価方法であって、前記制御手段で、アミノ酸の濃度値に関する予め取得した前記評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの前記濃度値および前記アミノ酸の濃度を変数とする前記記憶手段で記憶した多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出ステップと、前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記ストレスの状態を評価する判別値基準評価ステップとを実行し、前記多変量判別式は、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むことを特徴とする。
また、本発明にかかるストレス評価方法は、前記のストレス評価方法において、前記判別値基準評価ステップは、前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病または非抑うつ病であるか否かを判別、または前記大うつ病または非大うつ病であるか否かを判別する判別値基準判別ステップをさらに含み、前記多変量判別式は、1つの分数式または複数の前記分数式の和で表され、それを構成する前記分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むこと、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つであることを特徴とする。
また、本発明にかかるストレス評価方法は、前記のストレス評価方法において、前記判別値基準判別ステップで前記抑うつ病または前記非抑うつ病であるか否かを判別する場合、前記多変量判別式は、数式1または数式2であること、またはSer,Glu,Cys,Lysを前記変数とする前記ロジスティック回帰式またはSer,Glu,Cys,Lysを前記変数とする前記線形判別式であることを特徴とする。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価方法は、前記に記載のストレス評価方法において、前記判別値基準判別ステップで前記大うつ病または前記非大うつ病であるか否かを判別する場合、前記多変量判別式は、数式3、数式4、数式5、数式6または数式7であること、またはSer,Gln,Pro,ABA,Tyr,Trpを前記変数とする前記ロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを前記変数とする前記ロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを前記変数とする前記線形判別式またはAsn,Gln,Tau,Cit,ABAを前記変数とする前記線形判別式であることを特徴とする。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価方法は、前記のストレス評価方法において、前記制御手段で、前記アミノ酸濃度データと前記ストレスの状態を表す指標に関するストレス状態指標データとを含む前記記憶手段で記憶したストレス状態情報に基づいて、前記記憶手段で記憶する前記多変量判別式を作成する多変量判別式作成ステップをさらに実行し、前記多変量判別式作成ステップは、前記ストレス状態情報から所定の式作成手法に基づいて、前記多変量判別式の候補である候補多変量判別式を作成する候補多変量判別式作成ステップと、前記候補多変量判別式作成ステップで作成した前記候補多変量判別式を、所定の検証手法に基づいて検証する候補多変量判別式検証ステップと、前記候補多変量判別式検証ステップでの検証結果から所定の変数選択手法に基づいて前記候補多変量判別式の変数を選択することで、前記候補多変量判別式を作成する際に用いる前記ストレス状態情報に含まれる前記アミノ酸濃度データの組み合わせを選択する変数選択ステップと、をさらに含み、前記候補多変量判別式作成ステップ、前記候補多変量判別式検証ステップおよび前記変数選択ステップを繰り返し実行して蓄積した前記検証結果に基づいて、複数の前記候補多変量判別式の中から前記多変量判別式として採用する前記候補多変量判別式を選出することで、前記多変量判別式を作成することを特徴とする。
また、本発明はストレス評価システムに関するものであり、本発明にかかるストレス評価システムは、制御手段と記憶手段とを備え、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価するストレス評価装置と、アミノ酸の濃度値に関する前記評価対象のアミノ酸濃度データを提供する情報通信端末装置とを、ネットワークを介して通信可能に接続して構成されたストレス評価システムであって、前記情報通信端末装置は、前記評価対象の前記アミノ酸濃度データを前記ストレス評価装置へ送信するアミノ酸濃度データ送信手段と、前記ストレス評価装置から送信された前記ストレスの状態に関する前記評価対象の評価結果を受信する評価結果受信手段とを備え、前記ストレス評価装置の前記制御手段は、前記情報通信端末装置から送信された前記評価対象の前記アミノ酸濃度データを受信するアミノ酸濃度データ受信手段と、前記アミノ酸濃度データ受信手段で受信した前記評価対象の前記アミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの前記濃度値および前記アミノ酸の濃度を変数とする前記記憶手段で記憶した多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出手段と、前記判別値算出手段で算出した前記判別値に基づいて、前記評価対象につき、前記ストレスの状態を評価する判別値基準評価手段と、前記判別値基準評価手段での前記評価対象の前記評価結果を前記情報通信端末装置へ送信する評価結果送信手段と、を備え、前記多変量判別式は、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むことを特徴とする。
また、本発明にかかるストレス評価システムは、前記に記載のストレス評価システムにおいて、前記判別値基準評価手段は、前記判別値算出手段で算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病または非抑うつ病であるか否かを判別、または前記大うつ病または非大うつ病であるか否かを判別する判別値基準判別手段をさらに備え、前記多変量判別式は、1つの分数式または複数の前記分数式の和で表され、それを構成する前記分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むこと、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つであることを特徴とする。
また、本発明にかかるストレス評価システムは、前記に記載のストレス評価システムにおいて、前記判別値基準判別手段で前記抑うつ病または前記非抑うつ病であるか否かを判別する場合、前記多変量判別式は、数式1または数式2であること、またはSer,Glu,Cys,Lysを前記変数とする前記ロジスティック回帰式またはSer,Glu,Cys,Lysを前記変数とする前記線形判別式であることを特徴とする。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価システムは、前記に記載のストレス評価システムにおいて、前記判別値基準判別手段で前記大うつ病または前記非大うつ病であるか否かを判別する場合、前記多変量判別式は、数式3、数式4、数式5、数式6または数式7であること、またはSer,Gln,Pro,ABA,Tyr,Trpを前記変数とする前記ロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを前記変数とする前記ロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを前記変数とする前記線形判別式またはAsn,Gln,Tau,Cit,ABAを前記変数とする前記線形判別式であることを特徴とする。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価システムは、前記に記載のストレス評価システムにおいて、前記ストレス評価装置の前記制御手段は、前記アミノ酸濃度データと前記ストレスの状態を表す指標に関するストレス状態指標データとを含む前記記憶手段で記憶したストレス状態情報に基づいて、前記記憶手段で記憶する前記多変量判別式を作成する多変量判別式作成手段をさらに備え、前記多変量判別式作成手段は、前記ストレス状態情報から所定の式作成手法に基づいて、前記多変量判別式の候補である候補多変量判別式を作成する候補多変量判別式作成手段と、前記候補多変量判別式作成手段で作成した前記候補多変量判別式を、所定の検証手法に基づいて検証する候補多変量判別式検証手段と、前記候補多変量判別式検証手段での検証結果から所定の変数選択手法に基づいて前記候補多変量判別式の変数を選択することで、前記候補多変量判別式を作成する際に用いる前記ストレス状態情報に含まれる前記アミノ酸濃度データの組み合わせを選択する変数選択手段と、をさらに備え、前記候補多変量判別式作成手段、前記候補多変量判別式検証手段および前記変数選択手段を繰り返し実行して蓄積した前記検証結果に基づいて、複数の前記候補多変量判別式の中から前記多変量判別式として採用する前記候補多変量判別式を選出することで、前記多変量判別式を作成することを特徴とする。
また、本発明はストレス評価プログラムに関するものであり、本発明にかかるストレス評価プログラムは、制御手段と記憶手段とを備えた情報処理装置に実行させる、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価するストレス評価プログラムであって、前記制御手段に、アミノ酸の濃度値に関する予め取得した前記評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの前記濃度値および前記アミノ酸の濃度を変数とする前記記憶手段で記憶した多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出ステップと、前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記ストレスの状態を評価する判別値基準評価ステップとを実行させ、前記多変量判別式は、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むことを特徴とする。
また、本発明にかかるストレス評価プログラムは、前記に記載のストレス評価プログラムにおいて、前記判別値基準評価ステップは、前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病または非抑うつ病であるか否かを判別、または前記大うつ病または非大うつ病であるか否かを判別する判別値基準判別ステップをさらに含み、前記多変量判別式は、1つの分数式または複数の前記分数式の和で表され、それを構成する前記分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを前記変数として含むこと、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つであることを特徴とする。
また、本発明にかかるストレス評価プログラムは、前記に記載のストレス評価プログラムにおいて、前記判別値基準判別ステップで前記抑うつ病または前記非抑うつ病であるか否かを判別する場合、前記多変量判別式は、数式1または数式2であること、またはSer,Glu,Cys,Lysを前記変数とする前記ロジスティック回帰式またはSer,Glu,Cys,Lysを前記変数とする前記線形判別式であることを特徴とする。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価プログラムは、前記に記載のストレス評価プログラムにおいて、前記判別値基準判別ステップで前記大うつ病または前記非大うつ病であるか否かを判別する場合、前記多変量判別式は、数式3、数式4、数式5、数式6または数式7であること、またはSer,Gln,Pro,ABA,Tyr,Trpを前記変数とする前記ロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを前記変数とする前記ロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを前記変数とする前記線形判別式またはAsn,Gln,Tau,Cit,ABAを前記変数とする前記線形判別式であることを特徴とする。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価プログラムは、前記に記載のストレス評価プログラムにおいて、前記制御手段に、前記アミノ酸濃度データと前記ストレスの状態を表す指標に関するストレス状態指標データとを含む前記記憶手段で記憶したストレス状態情報に基づいて、前記記憶手段で記憶する前記多変量判別式を作成する多変量判別式作成ステップをさらに実行させ、前記多変量判別式作成ステップは、前記ストレス状態情報から所定の式作成手法に基づいて、前記多変量判別式の候補である候補多変量判別式を作成する候補多変量判別式作成ステップと、前記候補多変量判別式作成ステップで作成した前記候補多変量判別式を、所定の検証手法に基づいて検証する候補多変量判別式検証ステップと、前記候補多変量判別式検証ステップでの検証結果から所定の変数選択手法に基づいて前記候補多変量判別式の変数を選択することで、前記候補多変量判別式を作成する際に用いる前記ストレス状態情報に含まれる前記アミノ酸濃度データの組み合わせを選択する変数選択ステップと、をさらに含み、前記候補多変量判別式作成ステップ、前記候補多変量判別式検証ステップおよび前記変数選択ステップを繰り返し実行して蓄積した前記検証結果に基づいて、複数の前記候補多変量判別式の中から前記多変量判別式として採用する前記候補多変量判別式を選出することで、前記多変量判別式を作成することを特徴とする。
また、本発明は記録媒体に関するものであり、本発明にかかるコンピュータ読み取り可能な記録媒体は、前記のストレス評価プログラムを記録したことを特徴とする。
本発明にかかるストレスの評価方法によれば、評価対象から採取した血液から、アミノ酸の濃度値に関するアミノ酸濃度データを測定し、測定した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価するので、血液中のアミノ酸濃度のうちストレスの状態と関連するアミノ酸濃度を利用してストレスの状態を精度よく評価することができるという効果を奏する。
本発明にかかるストレスの評価方法によれば、測定した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値に基づいて、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別するので、血液中のアミノ酸濃度のうち抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸濃度を利用して抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別を精度よく行うことができるという効果を奏する。
また、本発明にかかるストレスの評価方法によれば、測定した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値およびアミノ酸の濃度を変数とする予め設定した多変量判別式であってLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含むものに基づいて、当該多変量判別式の値である判別値を算出し、算出した判別値に基づいて、評価対象につき、ストレスの状態を評価するので、血液中のアミノ酸濃度のうちストレスの状態と関連するアミノ酸の濃度を変数とする多変量判別式で得られる判別値を利用して、ストレスの状態を精度よく評価することができるという効果を奏する。
また、本発明にかかるストレスの評価方法によれば、算出した判別値に基づいて、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別し、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含む、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つであるので、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸変数を用いる多変量判別式(分数式等)で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
また、本発明にかかるストレスの評価方法によれば、抑うつ病または非抑うつ病であるか否かを判別する場合、多変量判別式は、数式1または数式2である、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式であるので、抑うつ病と非抑うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレスの評価方法によれば、大うつ病または非大うつ病であるか否かを判別する場合、多変量判別式は、数式3、数式4、数式5、数式6または数式7であること、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式であるので、大うつ病と非大うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、大うつ病と非大うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価装置、ストレス評価方法およびストレス評価プログラムによれば、アミノ酸の濃度値に関する予め取得した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値およびアミノ酸の濃度を変数とする記憶手段で記憶した多変量判別式であってLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含むものに基づいて、当該多変量判別式の値である判別値を算出し、算出した判別値に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価するので、血液中のアミノ酸濃度のうちストレスの状態と関連するアミノ酸の濃度を変数とする多変量判別式で得られる判別値を利用して、ストレスの状態を精度よく評価することができるという効果を奏する。
また、本発明にかかるストレス評価装置、ストレス評価方法およびストレス評価プログラムによれば、算出した判別値に基づいて、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別し、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含む、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つであるので、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸変数を用いる多変量判別式(分数式等)で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
また、本発明にかかるストレス評価装置、ストレス評価方法およびストレス評価プログラムによれば、抑うつ病または非抑うつ病であるか否かを判別する場合、多変量判別式は、数式1または数式2である、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式であるので、抑うつ病と非抑うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価装置、ストレス評価方法およびストレス評価プログラムによれば、大うつ病または非大うつ病であるか否かを判別する場合、多変量判別式は、数式3、数式4、数式5、数式6または数式7である、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式であるので、大うつ病と非大うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、大うつ病と非大うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価装置、ストレス評価方法およびストレス評価プログラムによれば、アミノ酸濃度データとストレスの状態を表す指標に関するストレス状態指標データとを含む記憶手段で記憶したストレス状態情報に基づいて、記憶手段で記憶する多変量判別式を作成する。具体的には、(1)ストレス状態情報から所定の式作成手法に基づいて、多変量判別式の候補である候補多変量判別式を作成し、(2)作成した候補多変量判別式を、所定の検証手法に基づいて検証し、(3)(2)での検証結果から所定の変数選択手法に基づいて候補多変量判別式の変数を選択することで、候補多変量判別式を作成する際に用いるストレス状態情報に含まれるアミノ酸濃度データの組み合わせを選択し、(4)(1)、(2)および(3)を繰り返し実行して蓄積した検証結果に基づいて、複数の候補多変量判別式の中から多変量判別式として採用する候補多変量判別式を選出することで、多変量判別式を作成する。これにより、抑うつ病または大うつ病を少なくとも含むストレスの状態の評価に最適な多変量判別式(具体的には、ストレスの状態と有意な相関がある多変量判別式(より具体的には、抑うつ病と非抑うつ病との2群判別に有用な多変量判別式、大うつ病と非大うつ病との2群判別に有用な多変量判別式))を作成することができるという効果を奏する。
また、本発明にかかるストレス評価システムによれば、まず、情報通信端末装置は、評価対象のアミノ酸濃度データをストレス評価装置へ送信する。そして、ストレス評価装置は、情報通信端末装置から送信された評価対象のアミノ酸濃度データを受信し、受信した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値およびアミノ酸の濃度を変数とする記憶手段で記憶した多変量判別式であってLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含むものに基づいて、当該多変量判別式の値である判別値を算出し、算出した判別値に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価し、その評価対象の評価結果を情報通信端末装置へ送信する。そして、情報通信端末装置は、ストレス評価装置から送信されたストレスの状態に関する評価対象の評価結果を受信する。これにより、血液中のアミノ酸濃度のうちストレスの状態と関連するアミノ酸の濃度を変数とする多変量判別式で得られる判別値を利用して、ストレスの状態を精度よく評価することができるという効果を奏する。
また、本発明にかかるストレス評価システムによれば、算出した判別値に基づいて、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別し、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含む、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つであるので、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸変数を用いる多変量判別式(分数式等)で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
また、本発明にかかるストレス評価システムによれば、抑うつ病または非抑うつ病であるか否かを判別する場合、多変量判別式は、数式1または数式2であること、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式であるので、抑うつ病と非抑うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価システムによれば、大うつ病または非大うつ病であるか否かを判別する場合、多変量判別式は、数式3、数式4、数式5、数式6または数式7であること、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式であるので、大うつ病と非大うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、大うつ病と非大うつ病の2群判別をさらに精度よく行うことができるという効果を奏する。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
また、本発明にかかるストレス評価システムによれば、アミノ酸濃度データとストレスの状態を表す指標に関するストレス状態指標データとを含む記憶手段で記憶したストレス状態情報に基づいて、記憶手段で記憶する多変量判別式を作成する。具体的には、(1)ストレス状態情報から所定の式作成手法に基づいて、多変量判別式の候補である候補多変量判別式を作成し、(2)作成した候補多変量判別式を、所定の検証手法に基づいて検証し、(3)(2)での検証結果から所定の変数選択手法に基づいて候補多変量判別式の変数を選択することで、候補多変量判別式を作成する際に用いるストレス状態情報に含まれるアミノ酸濃度データの組み合わせを選択し、(4)(1)、(2)および(3)を繰り返し実行して蓄積した検証結果に基づいて、複数の候補多変量判別式の中から多変量判別式として採用する候補多変量判別式を選出することで、多変量判別式を作成する。これにより、抑うつ病または大うつ病を少なくとも含むストレスの状態の評価に最適な多変量判別式(具体的には、ストレスの状態と有意な相関がある多変量判別式(より具体的には、抑うつ病と非抑うつ病との2群判別に有用な多変量判別式、大うつ病と非大うつ病との2群判別に有用な多変量判別式))を作成することができるという効果を奏する。
また、本発明にかかる記録媒体によれば、当該記録媒体に記録されたストレス評価プログラムをコンピュータに読み取らせて実行することでコンピュータにストレス評価プログラムを実行させるので、ストレス評価プログラムと同様の効果を得ることができるという効果を奏する。
なお、本発明は、ストレスの状態を評価する際(具体的には、抑うつ病または非抑うつ病であるか否かを判別する際、大うつ病または非大うつ病であるか否かを判別する際)、アミノ酸の濃度以外に、その他の代謝物(生体代謝物)の濃度や、メタボロームデータ(有機酸、脂肪酸、核酸など)、トランスクリプトームデータ、プロテオームデータ、タンパク質の発現量、被験者の年齢・性別、生体指標などをさらに用いてもかまわない。また、本発明は、ストレスの状態を評価する際(具体的には、抑うつ病または非抑うつ病であるか否かを判別する際、大うつ病または非大うつ病であるか否かを判別する際)、多変量判別式における変数として、アミノ酸の濃度以外に、その他の代謝物(生体代謝物)の濃度や、メタボロームデータ(有機酸、脂肪酸、核酸など)、トランスクリプトームデータ、プロテオームデータ、タンパク質の発現量、被験者の年齢・性別、生体指標などをさらに用いてもかまわない。
図1は、本発明の基本原理を示す原理構成図である。 図2は、第1実施形態にかかるストレスの評価方法の一例を示すフローチャートである。 図3は、本発明の基本原理を示す原理構成図である。 図4は、本システムの全体構成の一例を示す図である。 図5は、本システムの全体構成の他の一例を示す図である。 図6は、本システムのストレス評価装置100の構成の一例を示すブロック図である。 図7は、利用者情報ファイル106aに格納される情報の一例を示す図である。 図8は、アミノ酸濃度データファイル106bに格納される情報の一例を示す図である。 図9は、ストレス状態情報ファイル106cに格納される情報の一例を示す図である。 図10は、指定ストレス状態情報ファイル106dに格納される情報の一例を示す図である。 図11は、候補多変量判別式ファイル106e1に格納される情報の一例を示す図である。 図12は、検証結果ファイル106e2に格納される情報の一例を示す図である。 図13は、選択ストレス状態情報ファイル106e3に格納される情報の一例を示す図である。 図14は、多変量判別式ファイル106e4に格納される情報の一例を示す図である。 図15は、判別値ファイル106fに格納される情報の一例を示す図である。 図16は、評価結果ファイル106gに格納される情報の一例を示す図である。 図17は、多変量判別式作成部102hの構成を示すブロック図である。 図18は、判別値基準評価部102jの構成を示すブロック図である。 図19は、本システムのクライアント装置200の構成の一例を示すブロック図である。 図20は、本システムのデータベース装置400の構成の一例を示すブロック図である。 図21は、本システムで行うストレス評価サービス処理の一例を示すフローチャートである。 図22は、本システムのストレス評価装置100で行う多変量判別式作成処理の一例を示すフローチャートである。 図23は、健常群と抑うつ病所見群と大うつ病所見群の3群間のアミノ酸変数の分布を示す箱ひげ図である。 図24は、実施例2で示した健常群と抑うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の性能を持つ複数の指標式を示す図である。 図25は、実施例2で示した健常群と抑うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の性能を持つ複数の指標式を示す図である。 図26は、実施例2で示した健常群と抑うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の性能を持つ複数の指標式を示す図である。 図27は、健常群と抑うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図28は、実施例2で示した健常群と大うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の性能を持つ複数の指標式を示す図である。 図29は、実施例2で示した健常群と大うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の性能を持つ複数の指標式を示す図である。 図30は、実施例2で示した健常群と大うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の性能を持つ複数の指標式を示す図である。 図31は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図32は、実施例2で示した健常群と抑うつ病所見群の2群判別性能を最大化する線形判別式と同等の性能を持つ複数の指標式を示す図である。 図33は、実施例2で示した健常群と抑うつ病所見群の2群判別性能を最大化する線形判別式と同等の性能を持つ複数の指標式を示す図である。 図34は、実施例2で示した健常群と抑うつ病所見群の2群判別性能を最大化する線形判別式と同等の性能を持つ複数の指標式を示す図である。 図35は、健常群と抑うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図36は、実施例2で示した健常群と大うつ病所見群の2群判別性能を最大化する線形判別式と同等の性能を持つ複数の指標式を示す図である。 図37は、実施例2で示した健常群と大うつ病所見群の2群判別性能を最大化する線形判別式と同等の性能を持つ複数の指標式を示す図である。 図38は、実施例2で示した健常群と大うつ病所見群の2群判別性能を最大化する線形判別式と同等の性能を持つ複数の指標式を示す図である。 図39は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図40は、指標式1と同等の判別性能を有する複数の指標式を示す図である。 図41は、指標式1と同等の判別性能を有する複数の指標式を示す図である。 図42は、指標式1と同等の判別性能を有する複数の指標式を示す図である。 図43は、健常群と抑うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図44は、指標式3と同等の判別性能を有する複数の指標式を示す図である。 図45は、指標式3と同等の判別性能を有する複数の指標式を示す図である。 図46は、指標式3と同等の判別性能を有する複数の指標式を示す図である。 図47は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図48は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図49は、健常群と抑うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図50は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図51は、実施例7で示した健常群と大うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の判別性能を有する複数の指標式を示す図である。 図52は、実施例7で示した健常群と大うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の判別性能を有する複数の指標式を示す図である。 図53は、実施例7で示した健常群と大うつ病所見群の2群判別性能を最大化するロジスティック回帰式と同等の判別性能を有する複数の指標式を示す図である。 図54は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図55は、実施例7で示した健常群と大うつ病所見群の2群判別性能を最大化する線形判別式と同等の判別性能を有する複数の指標式を示す図である。 図56は、実施例7で示した健常群と大うつ病所見群の2群判別性能を最大化する線形判別式と同等の判別性能を有する複数の指標式を示す図である。 図57は、実施例7で示した健常群と大うつ病所見群の2群判別性能を最大化する線形判別式と同等の判別性能を有する複数の指標式を示す図である。 図58は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図59は、指標式7と同等の判別性能を有する複数の指標式を示す図である。 図60は、指標式7と同等の判別性能を有する複数の指標式を示す図である。 図61は、指標式7と同等の判別性能を有する複数の指標式を示す図である。 図62は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。 図63は、健常群と大うつ病所見群の2群間の診断性能を評価するためのROC曲線を示す図である。
符号の説明
100 ストレス評価装置
102 制御部
102a 要求解釈部
102b 閲覧処理部
102c 認証処理部
102d 電子メール生成部
102e Webページ生成部
102f 受信部
102g ストレス状態情報指定部
102h 多変量判別式作成部
102h1 候補多変量判別式作成部
102h2 候補多変量判別式検証部
102h3 変数選択部
102i 判別値算出部
102j 判別値基準評価部
102j1 判別値基準判別部
102k 結果出力部
102m 送信部
104 通信インターフェース部
106 記憶部
106a 利用者情報ファイル
106b アミノ酸濃度データファイル
106c ストレス状態情報ファイル
106d 指定ストレス状態情報ファイル
106e 多変量判別式関連情報データベース
106e1 候補多変量判別式ファイル
106e2 検証結果ファイル
106e3 選択ストレス状態情報ファイル
106e4 多変量判別式ファイル
106f 判別値ファイル
106g 評価結果ファイル
108 入出力インターフェース部
112 入力装置
114 出力装置
200 クライアント装置(情報通信端末装置)
300 ネットワーク
400 データベース装置
以下に、本発明にかかるストレスの評価方法の実施の形態(第1実施形態)、本発明にかかるストレス評価装置、ストレス評価方法、ストレス評価システム、ストレス評価プログラムおよび記録媒体の実施の形態(第2実施形態)を、図面に基づいて詳細に説明する。なお、本実施の形態により本発明が限定されるものではない。
[第1実施形態]
[1−1.本発明の概要]
ここでは、本発明にかかるストレスの評価方法の概要について図1を参照して説明する。図1は本発明の基本原理を示す原理構成図である。
まず、本発明では、評価対象(例えば動物やヒトなどの個体)から採取した血液から、アミノ酸の濃度値に関するアミノ酸濃度データを測定する(ステップS−11)。ここで、血中アミノ酸濃度の分析は次のように行った。採血した血液サンプルを、ヘパリン処理したチューブに採取し、採取した血液サンプルを遠心することにより血液から血漿を分離した。全ての血漿サンプルは、アミノ酸濃度の測定時まで−70℃で凍結保存した。アミノ酸濃度測定時には、スルホサリチル酸を添加し3%濃度調整により除蛋白処理を行い、測定には、ポストカラムでニンヒドリン反応を用いた高速液体クロマトグラフィー(HPLC)を原理としたアミノ酸分析機と高速液体クロマトグラフ質量分析機(LC/MS)を使用した。なお、アミノ酸濃度の単位は、例えばモル濃度や重量濃度、これらの濃度に任意の定数を加減乗除することで得られるものでもよい。
つぎに、本発明では、ステップS−11で測定した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA(BCAAはVal,Leu,Ileの総和を表す。以下同様。),AAA(AAAはTyr,Phe,Trpの総和を表す。以下同様。),LNAAs(LNAAsはVal,Leu,Ile,Tyr,Phe,Trpの総和を表す。以下同様。)のうち少なくとも1つの濃度値)に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価する(ステップS−12)。
以上、本発明によれば、評価対象から採取した血液からアミノ酸の濃度値に関するアミノ酸濃度データを測定し、測定した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つの濃度値)に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価する。これにより、血液中のアミノ酸濃度のうちストレスの状態と関連するアミノ酸濃度を利用してストレスの状態を精度よく評価することができる。
ここで、ステップS−12を実行する前に、ステップS−11で測定した評価対象のアミノ酸濃度データから欠損値や外れ値などのデータを除去してもよい。これにより、ストレスの状態をさらに精度よく評価することができる。
また、ステップS−12では、ステップS−11で測定した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つの濃度値)に基づいて、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別してもよい。具体的には、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つの濃度値)と予め設定された閾値(カットオフ値)とを比較することで、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別してもよい。これにより、血液中のアミノ酸濃度のうち抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸濃度を利用して抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別を精度よく行うことができる。
また、ステップS−12では、ステップS−11で測定した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つの濃度値)およびアミノ酸の濃度を変数とする予め設定した多変量判別式であってLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含むもの(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つを変数として含むもの)に基づいて、当該多変量判別式の値である判別値を算出し、算出した判別値に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価してもよい。これにより、血液中のアミノ酸濃度のうちストレスの状態と関連するアミノ酸の濃度を変数とする多変量判別式で得られる判別値を利用して、ストレスの状態を精度よく評価することができる。
また、ステップS−12では、算出した判別値に基づいて、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別してもよい。具体的には、判別値と予め設定された閾値(カットオフ値)とを比較することで、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別してもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含んでもよく(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つを変数として含んでもよく)、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つでもよい。これにより、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸変数を用いる多変量判別式(分数式等)で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別をさらに精度よく行うことができる。
なお、ステップS−12で抑うつ病または非抑うつ病であるか否かを判別する場合、多変量判別式は、数式1または数式2でもよく、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にAsn,Tau,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にAsn,Tau,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含んでもよい。これにより、抑うつ病と非抑うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別をさらに精度よく行うことができる。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、ステップS−12で大うつ病または非大うつ病であるか否かを判別する場合、多変量判別式は、数式3、数式4、数式5、数式6または数式7でもよく、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にPhe,Glu,Serのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にTau,Thr,Val,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にTau,Thr,Val,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にPhe,Glu,Serのうち少なくとも1つを変数として含んでもよい。これにより、大うつ病と非大うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、大うつ病と非大うつ病の2群判別をさらに精度よく行うことができる。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
なお、上記した各多変量判別式は、本出願人による国際出願である国際公開WO2004/052191号パンフレットに記載の方法や、本出願人による国際出願である国際出願番号PCT/JP2006/304398に記載の方法(後述する第2実施形態に記載の多変量判別式作成処理)で作成することができる。これら方法で得られた多変量判別式であれば、入力データとしてのアミノ酸濃度データにおけるアミノ酸濃度の単位に因らず、当該多変量判別式をストレスの状態の評価に好適に用いることができる。
また、分数式とは、当該分数式の分子がアミノ酸A,B,C,・・・の和で表わされ及び/又は当該分数式の分母がアミノ酸a,b,c,・・・の和で表わされるものである。また、分数式には、このような構成の分数式α,β,γ,・・・の和(例えばα+βのようなもの)も含まれる。また、分数式には、分割された分数式も含まれる。なお、分子や分母に用いられるアミノ酸にはそれぞれ適当な係数がついてもかまわない。また、分子や分母に用いられるアミノ酸は重複してもかまわない。また、各分数式に適当な係数がついてもかまわない。また、各変数の係数の値や定数項の値は、実数であればかまわない。また、分数式で、分子の変数と分母の変数を入れ替えた組み合わせは、目的変数との相関の正負の符号は概して逆転するが、それらの相関性は保たれるので、判別性では同等と見なせるので、分子の変数と分母の変数を入れ替えた組み合わせも、包含するものである。
また、多変量判別式とは、一般に多変量解析で用いられる式の形式を意味し、例えば重回帰式、多重ロジスティック回帰式、線形判別関数、マハラノビス距離、正準判別関数、サポートベクターマシン、決定木などを包含する。また、異なる形式の多変量判別式の和で示されるような式も含まれる。また、重回帰式、多重ロジスティック回帰式、正準判別関数においては各変数に係数および定数項が付加されるが、この場合の係数および定数項は、好ましくは実数であること、より好ましくはデータから判別を行うために得られた係数および定数項の99%信頼区間の範囲に属する値、さらに好ましくはデータから判別を行うために得られた係数および定数項の95%信頼区間の範囲に属する値であればかまわない。また、各係数の値、及びその信頼区間は、それを実数倍したものでもよく、定数項の値、及びその信頼区間は、それに任意の実定数を加減乗除したものでもよい。
そして、本発明は、ストレスの状態を評価する際(具体的には、抑うつ病または非抑うつ病であるか否かを判別する際、大うつ病または非大うつ病であるか否かを判別する際)、アミノ酸の濃度以外に、その他の代謝物(生体代謝物)の濃度や、メタボロームデータ(有機酸、脂肪酸、核酸など)、トランスクリプトームデータ、プロテオームデータ、タンパク質の発現量、被験者の年齢・性別、生体指標などをさらに用いてもかまわない。また、本発明は、ストレスの状態を評価する際(具体的には、抑うつ病または非抑うつ病であるか否かを判別する際、大うつ病または非大うつ病であるか否かを判別する際)、多変量判別式における変数として、アミノ酸の濃度以外に、その他の代謝物(生体代謝物)の濃度や、メタボロームデータ(有機酸、脂肪酸、核酸など)、トランスクリプトームデータ、プロテオームデータ、タンパク質の発現量、被験者の年齢・性別、生体指標などをさらに用いてもかまわない。
[1−2.第1実施形態にかかるストレスの評価方法]
ここでは、第1実施形態にかかるストレスの評価方法について図2を参照して説明する。図2は、第1実施形態にかかるストレスの評価方法の一例を示すフローチャートである。
まず、動物やヒトなどの個体から採取した血液から、アミノ酸の濃度値に関するアミノ酸濃度データを測定する(ステップSA−11)。なお、アミノ酸の濃度値の測定は、上述した方法で行う。
つぎに、ステップSA−11で測定した個体のアミノ酸濃度データから欠損値や外れ値などのデータを除去する(ステップSA−12)。
つぎに、ステップSA−12で欠損値や外れ値などのデータが除去された個体のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値と予め設定された閾値(カットオフ値)とを比較することで、個体につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別する(ステップSA−13)。
[1−3.第1実施形態のまとめ、およびその他の実施形態]
以上、詳細に説明したように、ストレスの評価方法によれば、(1)個体から採取した血液からアミノ酸濃度データを測定し、(2)測定した個体のアミノ酸濃度データから欠損値や外れ値などのデータを除去し、(3)欠損値や外れ値などのデータが除去された個体のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値と予め設定された閾値(カットオフ値)とを比較することで、個体につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別する。これにより、血液中のアミノ酸濃度のうち抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸濃度を利用して抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別を精度よく行うことができる。
ここで、ストレスの評価方法によれば、ステップSA−12で欠損値や外れ値などのデータが除去された個体のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値およびアミノ酸の濃度を変数とする予め設定した多変量判別式であってLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含むものに基づいて判別値を算出し、算出した判別値と予め設定された閾値(カットオフ値)とを比較することで、個体につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別してもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含んでもよく、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つでもよい。これにより、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸変数を用いる多変量判別式(分数式等)で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別をさらに精度よく行うことができる。
また、ステップSA−13で抑うつ病または非抑うつ病であるか否かを判別する場合、多変量判別式は、数式1または数式2でもよく、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にAsn,Tau,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にAsn,Tau,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含んでもよい。これにより、抑うつ病と非抑うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別をさらに精度よく行うことができる。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、ステップSA−13で大うつ病または非大うつ病であるか否かを判別する場合、多変量判別式は、数式3、数式4、数式5、数式6または数式7でもよく、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にPhe,Glu,Serのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にTau,Thr,Val,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にTau,Thr,Val,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にPhe,Glu,Serのうち少なくとも1つを変数として含んでもよい。これにより、大うつ病と非大うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、大うつ病と非大うつ病の2群判別をさらに精度よく行うことができる。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
なお、上記した各多変量判別式は、本出願人による国際出願である国際公開WO2004/052191号パンフレットに記載の方法や、本出願人による国際出願である国際出願番号PCT/JP2006/304398に記載の方法(後述する第2実施形態に記載の多変量判別式作成処理)で作成することができる。これら方法で得られた多変量判別式であれば、入力データとしてのアミノ酸濃度データにおけるアミノ酸濃度の単位に因らず、当該多変量判別式をストレスの状態の評価に好適に用いることができる。
[第2実施形態]
[2−1.本発明の概要]
ここでは、本発明にかかるストレス評価装置、ストレス評価方法、ストレス評価システム、ストレス評価プログラムおよび記録媒体の概要について、図3を参照して説明する。図3は本発明の基本原理を示す原理構成図である。
まず、本発明は、制御部で、アミノ酸の濃度値に関する予め取得した評価対象(例えば動物やヒトなどの個体)のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つの濃度値)およびアミノ酸の濃度を変数する記憶部で記憶した多変量判別式であってLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含むもの(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つを変数として含むもの)に基づいて、当該多変量判別式の値である判別値を算出する(ステップS−21)。
つぎに、本発明は、制御部で、ステップS−21で算出した判別値に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価する(ステップS−22)。
以上、本発明によれば、アミノ酸の濃度値に関する予め取得した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つの濃度値)およびアミノ酸の濃度を変数する記憶部で記憶した多変量判別式であってLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含むもの(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つを変数として含むもの)に基づいて判別値を算出し、算出した判別値に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価する。これにより、血液中のアミノ酸濃度のうちストレスの状態と関連するアミノ酸の濃度を変数とする多変量判別式で得られる判別値を利用して、ストレスの状態を精度よく評価することができる。
ここで、ステップS−22では、ステップS−21で算出した判別値に基づいて、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別してもよい。具体的には、判別値と予め設定された閾値(カットオフ値)とを比較することで、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別してもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含んでもよく(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つを変数として含んでもよく)、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つでもよい。これにより、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸変数を用いる多変量判別式(分数式等)で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別をさらに精度よく行うことができる。
また、ステップS−22で抑うつ病または非抑うつ病であるか否かを判別する場合、ステップS−21において多変量判別式は、多変量判別式は、数式1または数式2でもよく、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にAsn,Tau,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にAsn,Tau,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含んでもよい。これにより、抑うつ病と非抑うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別をさらに精度よく行うことができる。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、ステップS−22で大うつ病または非大うつ病であるか否かを判別する場合、ステップS−21において多変量判別式は、数式3、数式4、数式5、数式6または数式7でもよく、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にPhe,Glu,Serのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にTau,Thr,Val,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にTau,Thr,Val,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にPhe,Glu,Serのうち少なくとも1つを変数として含んでもよい。これにより、大うつ病と非大うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、大うつ病と非大うつ病の2群判別をさらに精度よく行うことができる。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
なお、上記した各多変量判別式は、本出願人による国際出願である国際公開WO2004/052191号パンフレットに記載の方法や、本出願人による国際出願である国際出願番号PCT/JP2006/304398に記載の方法(後述する多変量判別式作成処理)で作成することができる。これら方法で得られた多変量判別式であれば、入力データとしてのアミノ酸濃度データにおけるアミノ酸濃度の単位に因らず、当該多変量判別式をストレスの状態の評価に好適に用いることができる。
また、分数式とは、当該分数式の分子がアミノ酸A,B,C,・・・の和で表わされ及び/又は当該分数式の分母がアミノ酸a,b,c,・・・の和で表わされるものである。また、分数式には、このような構成の分数式α,β,γ,・・・の和(例えばα+βのようなもの)も含まれる。また、分数式には、分割された分数式も含まれる。なお、分子や分母に用いられるアミノ酸にはそれぞれ適当な係数がついてもかまわない。また、分子や分母に用いられるアミノ酸は重複してもかまわない。また、各分数式に適当な係数がついてもかまわない。また、各変数の係数の値や定数項の値は、実数であればかまわない。また、分数式で、分子の変数と分母の変数を入れ替えた組み合わせは、目的変数との相関の正負の符号は概して逆転するが、それらの相関性は保たれるので、判別性では同等と見なせるので、分子の変数と分母の変数を入れ替えた組み合わせも、包含するものである。
また、多変量判別式とは、一般に多変量解析で用いられる式の形式を意味し、例えば重回帰式、多重ロジスティック回帰式、線形判別関数、マハラノビス距離、正準判別関数、サポートベクターマシン、決定木などを包含する。また、異なる形式の多変量判別式の和で示されるような式も含まれる。また、重回帰式、多重ロジスティック回帰式、正準判別関数においては各変数に係数および定数項が付加されるが、この場合の係数および定数項は、好ましくは実数であること、より好ましくはデータから判別を行うために得られた係数および定数項の99%信頼区間の範囲に属する値、さらに好ましくはデータから判別を行うために得られた係数および定数項の95%信頼区間の範囲に属する値であればかまわない。また、各係数の値、及びその信頼区間は、それを実数倍したものでもよく、定数項の値、及びその信頼区間は、それに任意の実定数を加減乗除したものでもよい。
そして、本発明は、ストレスの状態を評価する際(具体的には、抑うつ病または非抑うつ病であるか否かを判別する際、大うつ病または非大うつ病であるか否かを判別する際)、アミノ酸の濃度以外に、その他の代謝物(生体代謝物)の濃度や、メタボロームデータ(有機酸、脂肪酸、核酸など)、トランスクリプトームデータ、プロテオームデータ、タンパク質の発現量、被験者の年齢・性別、生体指標などをさらに用いてもかまわない。また、本発明は、ストレスの状態を評価する際(具体的には、抑うつ病または非抑うつ病であるか否かを判別する際、大うつ病または非大うつ病であるか否かを判別する際)、多変量判別式における変数として、アミノ酸の濃度以外に、その他の代謝物(生体代謝物)の濃度や、メタボロームデータ(有機酸、脂肪酸、核酸など)、トランスクリプトームデータ、プロテオームデータ、タンパク質の発現量、被験者の年齢・性別、生体指標などをさらに用いてもかまわない。
ここで、多変量判別式作成処理(工程1〜工程4)の概要について詳細に説明する。
まず、本発明は、制御部で、アミノ酸濃度データとストレスの状態を表す指標に関するストレス状態指標データとを含む記憶部で記憶したストレス状態情報から所定の式作成手法に基づいて、多変量判別式の候補である候補多変量判別式(例えば、y=a+a+・・・+a、y:ストレス状態指標データ、x:アミノ酸濃度データ、a:定数、i=1,2,・・・,n)を作成する(工程1)。なお、事前に、ストレス状態情報から欠損値や外れ値などを持つデータを除去してもよい。
なお、工程1において、ストレス状態情報から、複数の異なる式作成手法(主成分分析や判別分析、サポートベクターマシン、重回帰分析、ロジスティック回帰分析、k−means法、クラスター解析、決定木などの多変量解析に関するものを含む。)を併用して複数の候補多変量判別式を作成してもよい。具体的には、多数の健常群およびストレス所見群から得た血液を分析して得たアミノ酸濃度データおよびストレス状態指標データから構成される多変量データであるストレス状態情報に対して、複数の異なるアルゴリズムを利用して複数群の候補多変量判別式を同時並行的に作成してもよい。例えば、異なるアルゴリズムを利用して判別分析およびロジスティック回帰分析を同時に行い、2つの異なる候補多変量判別式を作成してもよい。また、主成分分析を行って作成した候補多変量判別式を利用してストレス状態情報を変換し、変換したストレス状態情報に対して判別分析を行うことで候補多変量判別式を作成してもよい。これにより、最終的に、診断条件に合った適切な多変量判別式を作成することができる。
ここで、主成分分析を用いて作成した候補多変量判別式は、全てのアミノ酸濃度データの分散を最大にするような各アミノ酸変数からなる一次式である。また、判別分析を用いて作成した候補多変量判別式は、各群内の分散の和の全てのアミノ酸濃度データの分散に対する比を最小にするような各アミノ酸変数からなる高次式(指数や対数を含む)である。また、サポートベクターマシンを用いて作成した候補多変量判別式は、群間の境界を最大にするような各アミノ酸変数からなる高次式(カーネル関数を含む)である。また、重回帰分析を用いて作成した候補多変量判別式は、全てのアミノ酸濃度データからの距離の和を最小にするような各アミノ酸変数からなる高次式である。ロジスティック回帰分析を用いて作成した候補多変量判別式は、尤度を最大にするような各アミノ酸変数からなる一次式を指数とする自然対数を項に持つ分数式である。また、k−means法とは、各アミノ酸濃度データのk個近傍を探索し、近傍点の属する群の中で一番多いものをそのデータの所属群と定義し、入力されたアミノ酸濃度データの属する群と定義された群とが最も合致するようなアミノ酸変数を選択する手法である。また、クラスター解析とは、全てのアミノ酸濃度データの中で最も近い距離にある点同士をクラスタリング(群化)する手法である。また、決定木とは、アミノ酸変数に序列をつけて、序列が上位であるアミノ酸変数の取りうるパターンからアミノ酸濃度データの群を予測する手法である。
多変量判別式作成処理の説明に戻り、本発明は、制御部で、工程1で作成した候補多変量判別式を、所定の検証手法に基づいて検証(相互検証)する(工程2)。候補多変量判別式の検証は、工程1で作成した各候補多変量判別式に対して行う。
なお、工程2において、ブートストラップ法やホールドアウト法、リーブワンアウト法などのうち少なくとも1つに基づいて候補多変量判別式の判別率や感度、特異性、情報量基準などのうち少なくとも1つに関して検証してもよい。これにより、ストレス状態情報や診断条件を考慮した予測性または堅牢性の高い候補多変量判別式を作成することができる。
ここで、判別率とは、全入力データの中で、本発明で評価したストレスの状態が正しい割合である。また、感度とは、入力データに記載されたストレス所見になっているものの中で、本発明で評価したストレスの状態が正しい割合である。また、特異性とは、入力データに記載されたストレスの状態が健常になっているものの中で、本発明で評価したストレスの状態が正しい割合である。また、情報量基準とは、工程1で作成した候補多変量判別式のアミノ酸変数の数と、本発明で評価したストレスの状態および入力データに記載されたストレスの状態の差異と、を足し合わせたものである。また、予測性とは、候補多変量判別式の検証を繰り返すことで得られた判別率や感度、特異性を平均したものである。また、堅牢性とは、候補多変量判別式の検証を繰り返すことで得られた判別率や感度、特異性の分散である。
多変量判別式作成処理の説明に戻り、本発明は、制御部で、工程2での検証結果から所定の変数選択手法に基づいて候補多変量判別式の変数を選択することで、候補多変量判別式を作成する際に用いるストレス状態情報に含まれるアミノ酸濃度データの組み合わせを選択する(工程3)。アミノ酸変数の選択は、工程1で作成した各候補多変量判別式に対して行う。これにより、候補多変量判別式のアミノ酸変数を適切に選択することができる。そして、工程3で選択したアミノ酸濃度データを含むストレス状態情報を用いて再び工程1を実行する。
なお、工程3において、工程2での検証結果からステップワイズ法、ベストパス法、近傍探索法、遺伝的アルゴリズムのうち少なくとも1つに基づいて候補多変量判別式のアミノ酸変数を選択してもよい。
ここで、ベストパス法とは、候補多変量判別式に含まれるアミノ酸変数を1つずつ順次減らしていき、候補多変量判別式が与える評価指標を最適化することでアミノ酸変数を選択する方法である。
多変量判別式作成処理の説明に戻り、本発明は、制御部で、上述した工程1、工程2および工程3を繰り返し実行し、これにより蓄積した検証結果に基づいて、複数の候補多変量判別式の中から多変量判別式として採用する候補多変量判別式を選出することで、多変量判別式を作成する(工程4)。なお、候補多変量判別式の選出には、例えば、同じ式作成手法で作成した候補多変量判別式の中から最適なものを選出する場合と、すべての候補多変量判別式の中から最適なものを選出する場合とがある。
以上、説明したように、多変量判別式作成処理では、ストレス状態情報に基づいて、候補多変量判別式の作成、候補多変量判別式の検証および候補多変量判別式の変数の選択に関する処理を一連の流れで体系化(システム化)して実行することにより、ストレスの状態の評価に最適な多変量判別式を作成することができる。換言すると、多変量判別式作成処理では、アミノ酸濃度を多変量の統計解析に用い、最適でロバストな変数の組を選択するために変数選択法とクロスバリデーションとを組み合わせて、診断性能の高い多変量判別式を抽出する。多変量判別式としては、ロジスティック回帰、線形判別、サポートベクターマシン、マハラノビス距離法、重回帰分析、クラスター解析などを用いることができる。
[2−2.システム構成]
ここでは、第2実施形態にかかるストレス評価システム(以下では本システムと記す場合がある。)の構成について、図4から図20を参照して説明する。なお、本システムはあくまでも一例であり、本発明はこれに限定されない。
まず、本システムの全体構成について図4および図5を参照して説明する。図4は本システムの全体構成の一例を示す図である。また、図5は本システムの全体構成の他の一例を示す図である。本システムは、図4に示すように、評価対象につきストレスの状態を評価するストレス評価装置100と、アミノ酸の濃度値に関する評価対象のアミノ酸濃度データを提供するクライアント装置200(本発明の情報通信端末装置に相当)とを、ネットワーク300を介して通信可能に接続して構成されている。
なお、本システムは、図5に示すように、ストレス評価装置100やクライアント装置200の他に、ストレス評価装置100で多変量判別式を作成する際に用いるストレス状態情報やストレスの状態を評価するために用いる多変量判別式などを格納したデータベース装置400を、ネットワーク300を介して通信可能に接続して構成されてもよい。これにより、ネットワーク300を介して、ストレス評価装置100からクライアント装置200やデータベース装置400へ、あるいはクライアント装置200やデータベース装置400からストレス評価装置100へ、ストレスの状態に関する情報などが提供される。ここで、ストレスの状態に関する情報とは、ヒトを含む生物のストレスの状態に関する特定の項目について測定した値に関する情報である。また、ストレスの状態に関する情報は、ストレス評価装置100やクライアント装置200や他の装置(例えば各種の計測装置等)で生成され、主にデータベース装置400に蓄積される。
つぎに、本システムのストレス評価装置100の構成について図6から図18を参照して説明する。図6は、本システムのストレス評価装置100の構成の一例を示すブロック図であり、該構成のうち本発明に関係する部分のみを概念的に示している。
ストレス評価装置100は、当該ストレス評価装置を統括的に制御するCPU等の制御部102と、ルータ等の通信装置および専用線等の有線または無線の通信回線を介して当該ストレス評価装置をネットワーク300に通信可能に接続する通信インターフェース部104と、各種のデータベースやテーブルやファイルなどを格納する記憶部106と、入力装置112や出力装置114に接続する入出力インターフェース部108と、で構成されており、これら各部は任意の通信路を介して通信可能に接続されている。ここで、ストレス評価装置100は、各種の分析装置(例えばアミノ酸アナライザー等)と同一筐体で構成されてもよい。また、ストレス評価装置100の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷等に応じた任意の単位で、機能的または物理的に分散・統合して構成してもよい。例えば、処理の一部をCGI(Common Gateway Interface)を用いて実現してもよい。
記憶部106は、ストレージ手段であり、例えば、RAM・ROM等のメモリ装置や、ハードディスクのような固定ディスク装置、フレキシブルディスク、光ディスク等を用いることができる。記憶部106には、OS(Operating System)と協働してCPUに命令を与え各種処理を行うためのコンピュータプログラムが記録されている。記憶部106は、図示の如く、利用者情報ファイル106aと、アミノ酸濃度データファイル106bと、ストレス状態情報ファイル106cと、指定ストレス状態情報ファイル106dと、多変量判別式関連情報データベース106eと、判別値ファイル106fと、評価結果ファイル106gと、を格納する。
利用者情報ファイル106aは、利用者に関する利用者情報を格納する。図7は、利用者情報ファイル106aに格納される情報の一例を示す図である。利用者情報ファイル106aに格納される情報は、図7に示すように、利用者を一意に識別するための利用者IDと、利用者が正当な者であるか否かの認証を行うための利用者パスワードと、利用者の氏名と、利用者の所属する所属先を一意に識別するための所属先IDと、利用者の所属する所属先の部門を一意に識別するための部門IDと、部門名と、利用者の電子メールアドレスと、を相互に関連付けて構成されている。
図6に戻り、アミノ酸濃度データファイル106bは、アミノ酸の濃度値に関するアミノ酸濃度データを格納する。図8は、アミノ酸濃度データファイル106bに格納される情報の一例を示す図である。アミノ酸濃度データファイル106bに格納される情報は、図8に示すように、評価対象である個体(サンプル)を一意に識別するための個体番号と、アミノ酸濃度データとを相互に関連付けて構成されている。ここで、図8では、アミノ酸濃度データを数値、すなわち連続尺度として扱っているが、アミノ酸濃度データは名義尺度や順序尺度でもよい。なお、名義尺度や順序尺度の場合は、それぞれの状態に対して任意の数値を与えることで解析してもよい。また、アミノ酸濃度データに、他の生体情報(性差、年齢、身長、体重、BMI指数、腹囲、インスリン抵抗性指数、尿酸値、血糖値、中性脂肪、体脂肪率、総コレステロール、HDLコレステロール、LDLコレステロール、収縮期血圧、拡張期血圧、ヘモグロビンA1c、動脈硬化指数、喫煙の有無、飲酒量、心電図の波形を数値化したもの、酵素濃度、遺伝子発現量、脳波、心拍、コルチゾール、海馬体積、脳血流速度、脳由来神経栄養因子、アミラーゼ活性、クロモグラニンA、ドーパミン、アミノ酸以外の代謝物の濃度など)を組み合わせてもよい。
図6に戻り、ストレス状態情報ファイル106cは、多変量判別式を作成する際に用いるストレス状態情報を格納する。図9は、ストレス状態情報ファイル106cに格納される情報の一例を示す図である。ストレス状態情報ファイル106cに格納される情報は、図9に示すように、個体番号と、ストレスの状態を表す指標(指標T、指標T、指標T・・・)に関するストレス状態指標データ(T)と、アミノ酸濃度データと、を相互に関連付けて構成されている。ここで、図9では、ストレス状態指標データおよびアミノ酸濃度データを数値(すなわち連続尺度)として扱っているが、ストレス状態指標データおよびアミノ酸濃度データは名義尺度や順序尺度でもよい。なお、名義尺度や順序尺度の場合は、それぞれの状態に対して任意の数値を与えることで解析してもよい。また、ストレス状態指標データは、ストレスの状態のマーカーとなる既知の単一の状態指標であり、数値データを用いてもよい。
図6に戻り、指定ストレス状態情報ファイル106dは、後述するストレス状態情報指定部102gで指定したストレス状態情報を格納する。図10は、指定ストレス状態情報ファイル106dに格納される情報の一例を示す図である。指定ストレス状態情報ファイル106dに格納される情報は、図10に示すように、個体番号と、指定したストレス状態指標データと、指定したアミノ酸濃度データと、を相互に関連付けて構成されている。
図6に戻り、多変量判別式関連情報データベース106eは、後述する候補多変量判別式作成部102h1で作成した候補多変量判別式を格納する候補多変量判別式ファイル106e1と、後述する候補多変量判別式検証部102h2での検証結果を格納する検証結果ファイル106e2と、後述する変数選択部102h3で選択したアミノ酸濃度データの組み合わせを含むストレス状態情報を格納する選択ストレス状態情報ファイル106e3と、後述する多変量判別式作成部102hで作成した多変量判別式を格納する多変量判別式ファイル106e4と、で構成される。
候補多変量判別式ファイル106e1は、後述する候補多変量判別式作成部102h1で作成した候補多変量判別式を格納する。図11は、候補多変量判別式ファイル106e1に格納される情報の一例を示す図である。候補多変量判別式ファイル106e1に格納される情報は、図11に示すように、ランクと、候補多変量判別式(図11では、F(Gly,Leu,Phe,・・・)やF(Gly,Leu,Phe,・・・)、F(Gly,Leu,Phe,・・・)など)とを相互に関連付けて構成されている。
図6に戻り、検証結果ファイル106e2は、後述する候補多変量判別式検証部102h2での検証結果を格納する。図12は、検証結果ファイル106e2に格納される情報の一例を示す図である。検証結果ファイル106e2に格納される情報は、図12に示すように、ランクと、候補多変量判別式(図12では、F(Gly,Leu,Phe,・・・)やF(Gly,Leu,Phe,・・・)、F(Gly,Leu,Phe,・・・)など)と、各候補多変量判別式の検証結果(例えば各候補多変量判別式の評価値)と、を相互に関連付けて構成されている。
図6に戻り、選択ストレス状態情報ファイル106e3は、後述する変数選択部102h3で選択した変数に対応するアミノ酸濃度データの組み合わせを含むストレス状態情報を格納する。図13は、選択ストレス状態情報ファイル106e3に格納される情報の一例を示す図である。選択ストレス状態情報ファイル106e3に格納される情報は、図13に示すように、個体番号と、後述するストレス状態情報指定部102gで指定したストレス状態指標データと、後述する変数選択部102h3で選択したアミノ酸濃度データと、を相互に関連付けて構成されている。
図6に戻り、多変量判別式ファイル106e4は、後述する多変量判別式作成部102hで作成した多変量判別式を格納する。図14は、多変量判別式ファイル106e4に格納される情報の一例を示す図である。多変量判別式ファイル106e4に格納される情報は、図14に示すように、ランクと、多変量判別式(図14では、F(Phe,・・・)やF(Gly,Leu,Phe)、F(Gly,Leu,Phe,・・・)など)と、各式作成手法に対応する閾値と、各多変量判別式の検証結果(例えば各多変量判別式の評価値)と、を相互に関連付けて構成されている。
図6に戻り、判別値ファイル106fは、後述する判別値算出部102iで算出した判別値を格納する。図15は、判別値ファイル106fに格納される情報の一例を示す図である。判別値ファイル106fに格納される情報は、図15に示すように、評価対象である個体(サンプル)を一意に識別するための個体番号と、ランク(多変量判別式を一意に識別するための番号)と、判別値と、を相互に関連付けて構成されている。
図6に戻り、評価結果ファイル106gは、後述する判別値基準評価部102jでの評価結果(具体的には、後述する判別値基準判別部102j1での判別結果)を格納する。図16は、評価結果ファイル106gに格納される情報の一例を示す図である。評価結果ファイル106gに格納される情報は、評価対象である個体(サンプル)を一意に識別するための個体番号と、予め取得した評価対象のアミノ酸濃度データと、多変量判別式で算出した判別値と、ストレスの状態に関する評価結果(具体的には、抑うつ病または非抑うつ病であるか否かに関する判別結果、大うつ病または非大うつ病であるか否かに関する判別結果)と、を相互に関連付けて構成されている。
図6に戻り、記憶部106には、上述した情報以外にその他情報として、Webサイトをクライアント装置200に提供するための各種のWebデータや、CGIプログラム等が記録されている。Webデータとしては後述する各種のWebページを表示するためのデータ等があり、これらデータは例えばHTMLやXMLで記述されたテキストファイルとして形成されている。また、Webデータを作成するための部品用のファイルや作業用のファイルやその他一時的なファイル等も記憶部106に記憶される。記憶部106には、必要に応じて、クライアント装置200に送信するための音声をWAVE形式やAIFF形式の如き音声ファイルで格納したり、静止画や動画をJPEG形式やMPEG2形式の如き画像ファイルで格納したりすることができる。
通信インターフェース部104は、ストレス評価装置100とネットワーク300(またはルータ等の通信装置)との間における通信を媒介する。すなわち、通信インターフェース部104は、他の端末と通信回線を介してデータを通信する機能を有する。
入出力インターフェース部108は、入力装置112や出力装置114に接続する。ここで、出力装置114には、モニタ(家庭用テレビを含む)の他、スピーカやプリンタを用いることができる(なお、以下では、出力装置114をモニタ114として記載する場合がある。)。入力装置112には、キーボードやマウスやマイクの他、マウスと協働してポインティングデバイス機能を実現するモニタを用いることができる。
制御部102は、OS(Operating System)等の制御プログラム・各種の処理手順等を規定したプログラム・所要データなどを格納するための内部メモリを有し、これらのプログラムに基づいて種々の情報処理を実行する。制御部102は、図示の如く、大別して、要求解釈部102aと閲覧処理部102bと認証処理部102cと電子メール生成部102dとWebページ生成部102eと受信部102fとストレス状態情報指定部102gと多変量判別式作成部102hと判別値算出部102iと判別値基準評価部102jと結果出力部102kと送信部102mとを備えている。制御部102は、データベース装置400から送信されたストレス状態情報やクライアント装置200から送信されたアミノ酸濃度データに対して、欠損値のあるデータの除去・外れ値の多いデータの除去・欠損値のあるデータの多い変数の除去などのデータ処理も行う。
要求解釈部102aは、クライアント装置200やデータベース装置400からの要求内容を解釈し、その解釈結果に応じて制御部102の各部に処理を受け渡す。閲覧処理部102bは、クライアント装置200からの各種画面の閲覧要求を受けて、これら画面のWebデータの生成や送信を行なう。認証処理部102cは、クライアント装置200やデータベース装置400からの認証要求を受けて、認証判断を行う。電子メール生成部102dは、各種の情報を含んだ電子メールを生成する。Webページ生成部102eは、利用者がクライアント装置200で閲覧するWebページを生成する。
受信部102fは、クライアント装置200やデータベース装置400から送信された情報(具体的には、アミノ酸濃度データやストレス状態情報、多変量判別式など)を、ネットワーク300を介して受信する。ストレス状態情報指定部102gは、多変量判別式を作成するにあたり、対象とするストレス状態指標データおよびアミノ酸濃度データを指定する。
多変量判別式作成部102hは、受信部102fで受信したストレス状態情報やストレス状態情報指定部102gで指定したストレス状態情報に基づいて多変量判別式を作成する。具体的には、多変量判別式作成部102hは、ストレス状態情報から、候補多変量判別式作成部102h1、候補多変量判別式検証部102h2および変数選択部102h3を繰り返し実行させることにより蓄積された検証結果に基づいて、複数の候補多変量判別式の中から多変量判別式として採用する候補多変量判別式を選出することで、多変量判別式を作成する。
なお、多変量判別式が予め記憶部106の所定の記憶領域に格納されている場合には、多変量判別式作成部102hは、記憶部106から所望の多変量判別式を選択することで、多変量判別式を作成してもよい。また、多変量判別式作成部102hは、多変量判別式を予め格納した他のコンピュータ装置(例えばデータベース装置400)から所望の多変量判別式を選択しダウンロードすることで、多変量判別式を作成してもよい。
ここで、多変量判別式作成部102hの構成について図17を参照して説明する。図17は、多変量判別式作成部102hの構成を示すブロック図であり、該構成のうち本発明に関係する部分のみを概念的に示している。多変量判別式作成部102hは、候補多変量判別式作成部102h1と、候補多変量判別式検証部102h2と、変数選択部102h3と、をさらに備えている。候補多変量判別式作成部102h1は、ストレス状態情報から所定の式作成手法に基づいて多変量判別式の候補である候補多変量判別式を作成する。なお、候補多変量判別式作成部102h1は、ストレス状態情報から、複数の異なる式作成手法を併用して複数の候補多変量判別式を作成してもよい。候補多変量判別式検証部102h2は、候補多変量判別式作成部102h1で作成した候補多変量判別式を所定の検証手法に基づいて検証する。なお、候補多変量判別式検証部102h2は、ブートストラップ法、ホールドアウト法、リーブワンアウト法のうち少なくとも1つに基づいて候補多変量判別式の判別率、感度、特異性、情報量基準のうち少なくとも1つに関して検証してもよい。変数選択部102h3は、候補多変量判別式検証部102h2での検証結果から所定の変数選択手法に基づいて候補多変量判別式の変数を選択することで、候補多変量判別式を作成する際に用いるストレス状態情報に含まれるアミノ酸濃度データの組み合わせを選択する。なお、変数選択部102h3は、検証結果からステップワイズ法、ベストパス法、近傍探索法、遺伝的アルゴリズムのうち少なくとも1つに基づいて候補多変量判別式の変数を選択してもよい。
図6に戻り、判別値算出部102iは、受信部102fで受信した評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つの濃度値)および多変量判別式作成部102hで作成したLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含む多変量判別式(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つを変数として含む多変量判別式)に基づいて、当該多変量判別式の値である判別値を算出する。
ここで、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含んでもよく(例えば、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsのうち少なくとも1つを変数として含んでもよく)、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つでもよい。
具体的には、後述する判別値基準判別部102j1で抑うつ病または非抑うつ病であるか否かを判別する場合、多変量判別式は、数式1または数式2でもよく、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にAsn,Tau,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にAsn,Tau,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含んでもよい。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、後述する判別値基準判別部102j1で大うつ病または非大うつ病であるか否かを判別する場合、多変量判別式は、数式3、数式4、数式5、数式6または数式7でもよく、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にPhe,Glu,Serのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にTau,Thr,Val,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にTau,Thr,Val,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にPhe,Glu,Serのうち少なくとも1つを変数として含んでもよい。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
図6の説明に戻り、判別値基準評価部102jは、判別値算出部102iで算出した判別値に基づいて、評価対象につき、抑うつ病または大うつ病を少なくとも含むストレスの状態を評価する。判別値基準評価部102jは、判別値基準判別部102j1をさらに備えている。ここで、判別値基準評価部102jの構成について図18を参照して説明する。図18は、判別値基準評価部102jの構成を示すブロック図であり、該構成のうち本発明に関係する部分のみを概念的に示している。判別値基準判別部102j1は、判別値に基づいて評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別する。具体的には、判別値基準判別部102j1は、判別値と予め設定された閾値(カットオフ値)とを比較することで、評価対象につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別する。
図6に戻り、結果出力部102kは、制御部102の各処理部での処理結果(判別値基準評価部102jでの評価結果(具体的には判別値基準判別部102j1での判別結果)を含む)等を出力装置114に出力する。
送信部102mは、評価対象のアミノ酸濃度データの送信元のクライアント装置200に対して評価結果を送信したり、データベース装置400に対して、ストレス評価装置100で作成した多変量判別式や評価結果を送信したりする。
つぎに、本システムのクライアント装置200の構成について図19を参照して説明する。図19は、本システムのクライアント装置200の構成の一例を示すブロック図であり、該構成のうち本発明に関係する部分のみを概念的に示している。
クライアント装置200は、制御部210とROM220とHD230とRAM240と入力装置250と出力装置260と入出力IF270と通信IF280とで構成されており、これら各部は任意の通信路を介して通信可能に接続されている。
制御部210は、Webブラウザ211、電子メーラ212、受信部213、送信部214を備えている。Webブラウザ211は、Webデータを解釈し、解釈したWebデータを後述するモニタ261に表示するブラウズ処理を行う。なお、Webブラウザ211には、ストリーム映像の受信・表示・フィードバック等を行う機能を備えたストリームプレイヤ等の各種のソフトウェアをプラグインしてもよい。電子メーラ212は、所定の通信規約(例えば、SMTP(Simple Mail Transfer Protocol)やPOP3(Post Office Protocol version 3)等)に従って電子メールの送受信を行う。受信部213は、通信IF280を介して、ストレス評価装置100から送信された評価結果などの各種情報を受信する。送信部214は、通信IF280を介して、評価対象のアミノ酸濃度データなどの各種情報をストレス評価装置100へ送信する。
入力装置250はキーボードやマウスやマイク等である。なお、後述するモニタ261もマウスと協働してポインティングデバイス機能を実現する。出力装置260は、通信IF280を介して受信した情報を出力する出力手段であり、モニタ(家庭用テレビを含む)261およびプリンタ262を含む。この他、出力装置260にスピーカ等を設けてもよい。入出力IFは入力装置250や出力装置260に接続する。
通信IF280は、クライアント装置200とネットワーク300(またはルータ等の通信装置)とを通信可能に接続する。換言すると、クライアント装置200は、モデムやTAやルータなどの通信装置および電話回線を介して、または専用線を介してネットワーク300に接続される。これにより、クライアント装置200は、所定の通信規約に従ってストレス評価装置100にアクセスすることができる。
ここで、プリンタ・モニタ・イメージスキャナ等の周辺装置を必要に応じて接続した情報処理装置(例えば、既知のパーソナルコンピュータ・ワークステーション・家庭用ゲーム装置・インターネットTV・PHS端末・携帯端末・移動体通信端末・PDA等の情報処理端末など)に、Webデータのブラウジング機能や電子メール機能を実現させるソフトウェア(プログラム、データ等を含む)を実装することにより、クライアント装置200を実現してもよい。
また、クライアント装置200の制御部210は、制御部210で行う処理の全部または任意の一部を、CPUおよび当該CPUにて解釈して実行するプログラムで実現してもよい。ROM220またはHD230には、OS(Operating System)と協働してCPUに命令を与え、各種処理を行うためのコンピュータプログラムが記録されている。当該コンピュータプログラムは、RAM240にロードされることで実行され、CPUと協働して制御部210を構成する。また、当該コンピュータプログラムは、クライアント装置200と任意のネットワークを介して接続されるアプリケーションプログラムサーバに記録されてもよく、クライアント装置200は、必要に応じてその全部または一部をダウンロードしてもよい。また、制御部210で行う処理の全部または任意の一部を、ワイヤードロジック等によるハードウェアで実現してもよい。
つぎに、本システムのネットワーク300について図4、図5を参照して説明する。ネットワーク300は、ストレス評価装置100とクライアント装置200とデータベース装置400とを相互に通信可能に接続する機能を有し、例えばインターネットやイントラネットやLAN(有線/無線の双方を含む)等である。なお、ネットワーク300は、VANや、パソコン通信網や、公衆電話網(アナログ/デジタルの双方を含む)や、専用回線網(アナログ/デジタルの双方を含む)や、CATV網や、携帯回線交換網または携帯パケット交換網(IMT2000方式、GSM(登録商標)方式またはPDC/PDC−P方式等を含む)や、無線呼出網や、Bluetooth(登録商標)等の局所無線網や、PHS網や、衛星通信網(CS、BSまたはISDB等を含む)等でもよい。
つぎに、本システムのデータベース装置400の構成について図20を参照して説明する。図20は、本システムのデータベース装置400の構成の一例を示すブロック図であり、該構成のうち本発明に関係する部分のみを概念的に示している。
データベース装置400は、ストレス評価装置100または当該データベース装置で多変量判別式を作成する際に用いるストレス状態情報や、ストレス評価装置100で作成した多変量判別式、ストレス評価装置100での評価結果などを格納する機能を有する。図20に示すように、データベース装置400は、当該データベース装置を統括的に制御するCPU等の制御部402と、ルータ等の通信装置および専用線等の有線または無線の通信回路を介して当該データベース装置をネットワーク300に通信可能に接続する通信インターフェース部404と、各種のデータベースやテーブルやファイル(例えばWebページ用ファイル)などを格納する記憶部406と、入力装置412や出力装置414に接続する入出力インターフェース部408と、で構成されており、これら各部は任意の通信路を介して通信可能に接続されている。
記憶部406は、ストレージ手段であり、例えば、RAM・ROM等のメモリ装置や、ハードディスクのような固定ディスク装置や、フレキシブルディスクや、光ディスク等を用いることができる。記憶部406には、各種処理に用いる各種プログラム等を格納する。通信インターフェース部404は、データベース装置400とネットワーク300(またはルータ等の通信装置)との間における通信を媒介する。すなわち、通信インターフェース部404は、他の端末と通信回線を介してデータを通信する機能を有する。入出力インターフェース部408は、入力装置412や出力装置414に接続する。ここで、出力装置414には、モニタ(家庭用テレビを含む)の他、スピーカやプリンタを用いることができる(なお、以下で、出力装置414をモニタ414として記載する場合がある。)。また、入力装置412には、キーボードやマウスやマイクの他、マウスと協働してポインティングデバイス機能を実現するモニタを用いることができる。
制御部402は、OS(Operating System)等の制御プログラム・各種の処理手順等を規定したプログラム・所要データなどを格納するための内部メモリを有し、これらのプログラムに基づいて種々の情報処理を実行する。制御部402は、図示の如く、大別して、要求解釈部402aと閲覧処理部402bと認証処理部402cと電子メール生成部402dとWebページ生成部402eと送信部402fとを備えている。
要求解釈部402aは、ストレス評価装置100からの要求内容を解釈し、その解釈結果に応じて制御部402の各部に処理を受け渡す。閲覧処理部402bは、ストレス評価装置100からの各種画面の閲覧要求を受けて、これら画面のWebデータの生成や送信を行う。認証処理部402cは、ストレス評価装置100からの認証要求を受けて、認証判断を行う。電子メール生成部402dは、各種の情報を含んだ電子メールを生成する。Webページ生成部402eは、利用者がクライアント装置200で閲覧するWebページを生成する。送信部402fは、ストレス状態情報や多変量判別式などの各種情報を、ストレス評価装置100へ送信する。
[2−3.本システムの処理]
ここでは、以上のように構成された本システムで行われるストレス評価サービス処理の一例を、図21を参照して説明する。図21は、ストレス評価サービス処理の一例を示すフローチャートである。
なお、本処理で用いるアミノ酸濃度データは、個体から予め採取した血液を分析して得たアミノ酸の濃度値に関するものである。ここで、血液中のアミノ酸の分析方法について簡単に説明する。まず、採血した血液サンプルを、ヘパリン処理したチューブに採取し、その後、当該チューブに対して遠心分離を行うことで血漿を分離する。なお、分離したすべての血漿サンプルは、アミノ酸濃度の測定時まで−70℃で凍結保存する。そして、アミノ酸濃度の測定時に、血漿サンプルに対してスルホサリチル酸を添加し、3%濃度調整により除蛋白処理を行う。なお、アミノ酸濃度の測定には、ポストカラムでニンヒドリン反応を用いた高速液体クロマトグラフィー(HPLC)を原理としたアミノ酸分析機と高速液体クロマトグラフ質量分析機(LC/MS)を使用した。
まず、Webブラウザ211を表示した画面上で利用者が入力装置250を介してストレス評価装置100が提供するWebサイトのアドレス(URLなど)を指定すると、クライアント装置200はストレス評価装置100へアクセスする。具体的には、利用者がクライアント装置200のWebブラウザ211の画面更新を指示すると、Webブラウザ211は、ストレス評価装置100が提供するWebサイトのアドレスを所定の通信規約でストレス評価装置100へ送信することで、アミノ酸濃度データ送信画面に対応するWebページの送信要求を、当該アドレスに基づくルーティングでストレス評価装置100へ行う。
つぎに、ストレス評価装置100は、要求解釈部102aで、クライアント装置200からの送信を受け、当該送信の内容を解析し、解析結果に応じて制御部102の各部に処理を移す。具体的には、送信の内容がアミノ酸濃度データ送信画面に対応するWebページの送信要求であった場合、ストレス評価装置100は、主として閲覧処理部102bで、記憶部106の所定の記憶領域に格納されている当該Webページを表示するためのWebデータを取得し、取得したWebデータをクライアント装置200へ送信する。より具体的には、利用者からアミノ酸濃度データ送信画面に対応するWebページの送信要求があった場合、ストレス評価装置100は、まず、制御部102で、利用者IDや利用者パスワードの入力を利用者に対して求める。そして、利用者IDやパスワードが入力されると、ストレス評価装置100は、認証処理部102cで、入力された利用者IDやパスワードと利用者情報ファイル106aに格納されている利用者IDや利用者パスワードとの認証判断を行う。そして、ストレス評価装置100は、認証可の場合にのみ、閲覧処理部102bで、アミノ酸濃度データ送信画面に対応するWebページを表示するためのWebデータをクライアント装置200へ送信する。なお、クライアント装置200の特定は、クライアント装置200から送信要求と共に送信されたIPアドレスで行う。
つぎに、クライアント装置200は、ストレス評価装置100から送信されたWebデータ(アミノ酸濃度データ送信画面に対応するWebページを表示するためのもの)を受信部213で受信し、受信したWebデータをWebブラウザ211で解釈し、モニタ261にアミノ酸濃度データ送信画面を表示する。
つぎに、モニタ261に表示されたアミノ酸濃度データ送信画面に対し利用者が入力装置250を介して個体のアミノ酸濃度データなどを入力・選択すると、クライアント装置200は、送信部214で、入力情報や選択事項を特定するための識別子をストレス評価装置100へ送信することで、評価対象の個体のアミノ酸濃度データをストレス評価装置100へ送信する(ステップSA−21)。なお、ステップSA−21におけるアミノ酸濃度データの送信は、FTP等の既存のファイル転送技術等により実現してもよい。
つぎに、ストレス評価装置100は、要求解釈部102aで、クライアント装置200から送信された識別子を解釈することによりクライアント装置200の要求内容を解釈し、Lys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含むストレス評価用(具体的には、抑うつ病と非抑うつ病との2群判別用、大うつ病と非大うつ病との2群判別用)の多変量判別式の送信要求をデータベース装置400へ行う。
つぎに、データベース装置400は、要求解釈部402aで、ストレス評価装置100からの送信要求を解釈し、記憶部406の所定の記憶領域に格納したLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含む多変量判別式(例えばアップデートされた最新のもの)をストレス評価装置100へ送信する(ステップSA−22)。
ここで、ステップSA−22において、ストレス評価装置100へ送信する多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子および/または分母にLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つを変数として含んでもよく、またはロジスティック回帰式、線形判別式、重回帰式、サポートベクターマシンで作成された式、マハラノビス距離法で作成された式、正準判別分析で作成された式、決定木で作成された式のいずれか1つでもよい。
具体的には、後述するステップSA−26で抑うつ病または非抑うつ病であるか否かを判別する場合、多変量判別式は、数式1または数式2でもよく、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にAsn,Tau,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にAsn,Tau,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含んでもよい。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、後述するステップSA−26で大うつ病または非大うつ病であるか否かを判別する場合、多変量判別式は、数式3、数式4、数式5、数式6または数式7でもよく、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にPhe,Glu,Serのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にTau,Thr,Val,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にTau,Thr,Val,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にPhe,Glu,Serのうち少なくとも1つを変数として含んでもよい。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
図21の説明に戻り、ストレス評価装置100は、受信部102fで、クライアント装置200から送信された個体のアミノ酸濃度データおよびデータベース装置400から送信された多変量判別式を受信し、受信したアミノ酸濃度データをアミノ酸濃度データファイル106bの所定の記憶領域に格納すると共に、受信した多変量判別式を多変量判別式ファイル106e4の所定の記憶領域に格納する(ステップSA−23)。
つぎに、ストレス評価装置100は、制御部102で、ステップSA−23で受信した個体のアミノ酸濃度データから欠損値や外れ値などのデータを除去する(ステップSA−24)。
つぎに、ストレス評価装置100は、判別値算出部102iで、ステップSA−24で欠損値や外れ値などのデータが除去された個体のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値およびステップSA−23で受信した多変量判別式に基づいて、判別値を算出する(ステップSA−25)。
つぎに、ストレス評価装置100は、判別値基準判別部102j1で、ステップSA−25で算出した判別値と予め設定された閾値(カットオフ値)とを比較することで、個体につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かの判別し、その判別結果を評価結果ファイル106gの所定の記憶領域に格納する(ステップSA−26)。
つぎに、ストレス評価装置100は、送信部102mで、ステップSA−26で得た判別結果(抑うつ病または非抑うつ病であるか否かに関する判別結果、大うつ病または非大うつ病であるか否かに関する判別結果)を、アミノ酸濃度データの送信元のクライアント装置200とデータベース装置400とへ送信する(ステップSA−27)。具体的には、まず、ストレス評価装置100は、Webページ生成部102eで、判別結果を表示するためのWebページを作成し、作成したWebページに対応するWebデータを記憶部106の所定の記憶領域に格納する。ついで、利用者がクライアント装置200のWebブラウザ211に入力装置250を介して所定のURLを入力し上述した認証を経た後、クライアント装置200は、当該Webページの閲覧要求をストレス評価装置100へ送信する。ついで、ストレス評価装置100は、閲覧処理部102bで、クライアント装置200から送信された閲覧要求を解釈し、判別結果を表示するためのWebページに対応するWebデータを記憶部106の所定の記憶領域から読み出す。そして、ストレス評価装置100は、送信部102mで、読み出したWebデータをクライアント装置200へ送信すると共に、当該Webデータ又は判別結果をデータベース装置400へ送信する。
ここで、ステップSA−27において、ストレス評価装置100は、制御部102で、判別結果を電子メールで利用者のクライアント装置200へ通知してもよい。具体的には、まず、ストレス評価装置100は、電子メール生成部102dで、利用者IDなどを基にして利用者情報ファイル106aに格納されている利用者情報を送信タイミングに従って参照し、利用者の電子メールアドレスを取得する。ついで、ストレス評価装置100は、電子メール生成部102dで、取得した電子メールアドレスを宛て先とし利用者の氏名および判別結果を含む電子メールに関するデータを生成する。ついで、ストレス評価装置100は、送信部102mで、生成した当該データを利用者のクライアント装置200へ送信する。
また、ステップSA−27において、ストレス評価装置100は、FTP等の既存のファイル転送技術等で、判別結果を利用者のクライアント装置200へ送信してもよい。
図21の説明に戻り、データベース装置400は、制御部402で、ストレス評価装置100から送信された判別結果またはWebデータを受信し、受信した判別結果またはWebデータを記憶部406の所定の記憶領域に保存(蓄積)する(ステップSA−28)。
また、クライアント装置200は、受信部213で、ストレス評価装置100から送信されたWebデータを受信し、受信したWebデータをWebブラウザ211で解釈し、個体の判別結果が記されたWebページの画面をモニタ261に表示する(ステップSA−29)。なお、判別結果がストレス評価装置100から電子メールで送信された場合には、クライアント装置200は、電子メーラ212の公知の機能で、ストレス評価装置100から送信された電子メールを任意のタイミングで受信し、受信した電子メールをモニタ261に表示する。
以上により、利用者は、モニタ261に表示されたWebページを閲覧することで、抑うつ病と非抑うつ病との2群判別に関する個体の判別結果、大うつ病と非大うつ病との2群判別に関する個体の判別結果を確認することができる。なお、利用者は、モニタ261に表示されたWebページの表示内容をプリンタ262で印刷してもよい。
また、判別結果がストレス評価装置100から電子メールで送信された場合には、利用者は、モニタ261に表示された電子メールを閲覧することで、抑うつ病と非抑うつ病との2群判別に関する個体の判別結果、大うつ病と非大うつ病との2群判別に関する個体の判別結果を確認することができる。利用者は、モニタ261に表示された電子メールの表示内容をプリンタ262で印刷してもよい。
これにて、ストレス評価サービス処理の説明を終了する。
[2−4.第2実施形態のまとめ、およびその他の実施形態]
以上、詳細に説明したように、ストレス評価システムによれば、クライアント装置200は個体のアミノ酸濃度データをストレス評価装置100へ送信し、データベース装置400はストレス評価装置100からの要求を受けて、抑うつ病と非抑うつ病との2群判別用または大うつ病と非大うつ病との2群判別用の多変量判別式をストレス評価装置100へ送信する。そして、ストレス評価装置100は、(1)クライアント装置200からアミノ酸濃度データを受信すると共にデータベース装置400から多変量判別式を受信し、(2)受信した個体のアミノ酸濃度データから欠損値や外れ値などのデータを除去し、(3)欠損値や外れ値などのデータが除去されたアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも1つの濃度値および受信した多変量判別式に基づいて判別値を算出し、(4)算出した判別値と予め設定した閾値とを比較することで個体につき、抑うつ病または非抑うつ病であるか否かを判別、または大うつ病または非大うつ病であるか否かを判別し、(5)この判別結果をクライアント装置200やデータベース装置400へ送信する。そして、クライアント装置200はストレス評価装置100から送信された判別結果を受信して表示し、データベース装置400はストレス評価装置100から送信された判別結果を受信して格納する。これにより、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別に有用なアミノ酸変数を用いる多変量判別式(分数式等)で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別や大うつ病と非大うつ病の2群判別をさらに精度よく行うことができる。
ここで、ストレス評価システムによれば、ステップSA−26で抑うつ病または非抑うつ病であるか否かを判別する場合、多変量判別式は、数式1または数式2でもよく、またはSer,Glu,Cys,Lysを変数とするロジスティック回帰式またはSer,Glu,Cys,Lysを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にAsn,Tau,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にAsn,Tau,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にCys,Glu,Tyr,Pheのうち少なくとも1つを変数として含んでもよい。これにより、抑うつ病と非抑うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、抑うつ病と非抑うつ病の2群判別をさらに精度よく行うことができる。
×Cys/Asn + b×(Glu+Tyr+Phe)/(Tau+Lys) + c ・・・(数式1)
×Glu/Lys + b×Met/(Tyr+Phe+Trp) + c×Ile/Leu + d×Phe/(Val+Leu+Ile+Tyr+Phe+Trp) + e ・・・(数式2)
(数式1において、a、b、cは任意の実数であり、数式2において、a、b、c、d、eは任意の実数である。)
また、ステップSA−26で大うつ病または非大うつ病であるか否かを判別する場合、多変量判別式は、数式3、数式4、数式5、数式6または数式7でもよく、またはSer,Gln,Pro,ABA,Tyr,Trpを変数とするロジスティック回帰式、Glu,Asn,Gln,Cit,ABAを変数とするロジスティック回帰式、Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpを変数とする線形判別式またはAsn,Gln,Tau,Cit,ABAを変数とする線形判別式でもよい。また、多変量判別式は、1つの分数式または複数の分数式の和で表され、それを構成する分数式の分子にPhe,Glu,Serのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にTau,Thr,Val,Lysのうち少なくとも1つを変数として含む、またはそれを構成する分数式の分子にTau,Thr,Val,Lysのうち少なくとも1つを変数として含み且つそれを構成する分数式の分母にPhe,Glu,Serのうち少なくとも1つを変数として含んでもよい。これにより、大うつ病と非大うつ病の2群判別に特に有用なアミノ酸変数を用いる多変量判別式で得られる判別値を利用して、大うつ病と非大うつ病の2群判別をさらに精度よく行うことができる。
×Phe/(Tau+Thr) + b×(Glu+Ser)/(Val+Lys) + c ・・・(数式3)
×Gly/Gln +b×(Glu+Pro+Ile)/(Orn+Trp+ABA) + c ・・・(数式4)
×Asn + b×Ala/Trp + c×ABA/Tyr + d×Ser/Val + e ・・・(数式5)
×(Glu+Gln)/(Val+Thr+Ser) + b×(Cit+Met)/Asn + c ・・・(数式6)
×Asn/Gln + b×Cit/ABA + c×Glu/Arg +d×Met/(Val+Leu+Ile) + e ・・・(数式7)
(数式3において、a、b、cは任意の実数であり、数式4において、a、b、cは任意の実数であり、数式5において、a、b、c、d、eは任意の実数であり、数式6において、a、b、cは任意の実数であり、数式7において、a、b、c、d、eは任意の実数である。)
なお、上記した各多変量判別式は、本出願人による国際出願である国際公開WO2004/052191号パンフレットに記載の方法や、本出願人による国際出願である国際出願番号PCT/JP2006/304398に記載の方法(後述する多変量判別式作成処理)で作成することができる。これら方法で得られた多変量判別式であれば、入力データとしてのアミノ酸濃度データにおけるアミノ酸濃度の単位に因らず、当該多変量判別式をストレスの状態の評価に好適に用いることができる。
また、本発明にかかるストレス評価装置、ストレス評価方法、ストレス評価システム、ストレス評価プログラムおよび記録媒体は、上述した第2実施形態以外にも、特許請求の範囲の書類に記載した技術的思想の範囲内において種々の異なる実施形態にて実施されてよいものである。例えば、上述した第2実施形態で説明した各処理のうち、自動的に行なわれるものとして説明した処理の全部または一部を手動的に行うこともでき、手動的に行なわれるものとして説明した処理の全部または一部を公知の方法で自動的に行うこともできる。この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種の登録データおよび検索条件等のパラメータを含む情報、画面例、データベース構成については、特記する場合を除いて任意に変更することができる。例えば、ストレス評価装置100に関して、図示の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。また、ストレス評価装置100の各部または各装置が備える処理機能(特に制御部102にて行なわれる各処理機能)については、CPU(Central Processing Unit)および当該CPUにて解釈実行されるプログラムにて、その全部または任意の一部を実現することができ、ワイヤードロジックによるハードウェアとして実現することもできる。
ここで、「プログラム」とは任意の言語や記述方法にて記述されたデータ処理方法であり、ソースコードやバイナリコード等の形式を問わない。なお、「プログラム」は、必ずしも単一的に構成されるものに限られず、複数のモジュールやライブラリとして分散構成されるものや、OS(Operating System)に代表される別個のプログラムと協働してその機能を達成するものを含む。なお、プログラムは、記録媒体に記録されており、必要に応じてストレス評価装置100に機械的に読み取られる。記録媒体に記録されたプログラムを各装置で読み取るための具体的な構成や読み取り手順や読み取り後のインストール手順等については、周知の構成や手順を用いることができる。
また、「記録媒体」とは任意の「可搬用の物理媒体」や任意の「固定用の物理媒体」や「通信媒体」を含むものとする。なお、「可搬用の物理媒体」とはフレキシブルディスクや光磁気ディスクやROMやEPROMやEEPROM(登録商標)やCD−ROMやMOやDVD等である。「固定用の物理媒体」とは各種コンピュータシステムに内蔵されるROMやRAMやHD等である。「通信媒体」とは、LANやWANやインターネット等のネットワークを介してプログラムを送信する場合における通信回線や搬送波のように、短期にプログラムを保持するものである。
最後に、ストレス評価装置100で行う多変量判別式作成処理の一例について図22を参照して詳細に説明する。図22は多変量判別式作成処理の一例を示すフローチャートである。なお、当該多変量判別式作成処理は、ストレス状態情報を管理するデータベース装置400で行ってもよい。
なお、本説明では、ストレス評価装置100は、データベース装置400から事前に取得したストレス状態情報を、ストレス状態情報ファイル106cの所定の記憶領域に格納しているものとする。また、ストレス評価装置100は、ストレス状態情報指定部102gで事前に指定したストレス状態指標データおよびアミノ酸濃度データを含むストレス状態情報を、指定ストレス状態情報ファイル106dの所定の記憶領域に格納しているものとする。
まず、多変量判別式作成部102hは、候補多変量判別式作成部102h1で、指定ストレス状態情報ファイル106dの所定の記憶領域に格納されているストレス状態情報から所定の式作成手法に基づいて候補多変量判別式を作成し、作成した候補多変量判別式を候補多変量判別式ファイル106e1の所定の記憶領域に格納する(ステップSB−21)。具体的には、まず、多変量判別式作成部102hは、候補多変量判別式作成部102h1で、複数の異なる式作成手法(主成分分析や判別分析、サポートベクターマシン、重回帰分析、ロジスティック回帰分析、k−means法、クラスター解析、決定木などの多変量解析に関するものを含む。)の中から所望のものを1つ選択し、選択した式作成手法に基づいて、作成する候補多変量判別式の形(式の形)を決定する。つぎに、多変量判別式作成部102hは、候補多変量判別式作成部102h1で、ストレス状態情報に基づいて、選択した式選択手法に対応する種々(例えば平均や分散など)の計算を実行する。つぎに、多変量判別式作成部102hは、候補多変量判別式作成部102h1で、計算結果および決定した候補多変量判別式のパラメータを決定する。これにより、選択した式作成手法に基づいて候補多変量判別式が作成される。なお、複数の異なる式作成手法を併用して候補多変量判別式を同時並行(並列)的に作成する場合は、選択した式作成手法ごとに上記の処理を並行して実行すればよい。また、複数の異なる式作成手法を併用して候補多変量判別式を直列的に作成する場合は、例えば、主成分分析を行って作成した候補多変量判別式を利用してストレス状態情報を変換し、変換したストレス状態情報に対して判別分析を行うことで候補多変量判別式を作成してもよい。
つぎに、多変量判別式作成部102hは、候補多変量判別式検証部102h2で、ステップSB−21で作成した候補多変量判別式を所定の検証手法に基づいて検証(相互検証)し、検証結果を検証結果ファイル106e2の所定の記憶領域に格納する(ステップSB−22)。具体的には、多変量判別式作成部102hは、候補多変量判別式検証部102h2で、指定ストレス状態情報ファイル106dの所定の記憶領域に格納されているストレス状態情報に基づいて候補多変量判別式を検証する際に用いる検証用データを作成し、作成した検証用データに基づいて候補多変量判別式を検証する。なお、ステップSB−21で複数の異なる式作成手法を併用して候補多変量判別式を複数作成した場合には、多変量判別式作成部102hは、候補多変量判別式検証部102h2で、各式作成手法に対応する候補多変量判別式ごとに所定の検証手法に基づいて検証する。ここで、ステップSB−22において、ブートストラップ法やホールドアウト法、リーブワンアウト法などのうち少なくとも1つに基づいて候補多変量判別式の判別率や感度、特異性、情報量基準などのうち少なくとも1つに関して検証してもよい。これにより、ストレス状態情報や診断条件を考慮した予測性または堅牢性の高い候補指標式を選択することができる。
つぎに、多変量判別式作成部102hは、変数選択部102h3で、ステップSB−22での検証結果から所定の変数選択手法に基づいて、候補多変量判別式の変数を選択することで、候補多変量判別式を作成する際に用いるストレス状態情報に含まれるアミノ酸濃度データの組み合わせを選択し、選択したアミノ酸濃度データの組み合わせを含むストレス状態情報を選択ストレス状態情報ファイル106e3の所定の記憶領域に格納する(ステップSB−23)。なお、ステップSB−21で複数の異なる式作成手法を併用して候補多変量判別式を複数作成し、ステップSB−22で各式作成手法に対応する候補多変量判別式ごとに所定の検証手法に基づいて検証した場合には、ステップSB−23において、多変量判別式作成部102hは、変数選択部102h3で、ステップSB−22での検証結果に対応する候補多変量判別式ごとに所定の変数選択手法に基づいて候補多変量判別式の変数を選択する。ここで、ステップSB−23において、検証結果からステップワイズ法、ベストパス法、近傍探索法、遺伝的アルゴリズムのうち少なくとも1つに基づいて候補多変量判別式の変数を選択してもよい。なお、ベストパス法とは、候補多変量判別式に含まれる変数を1つずつ順次減らしていき、候補多変量判別式が与える評価指標を最適化することで変数を選択する方法である。また、ステップSB−23において、多変量判別式作成部102hは、変数選択部102h3で、指定ストレス状態情報ファイル106dの所定の記憶領域に格納されているストレス状態情報に基づいてアミノ酸濃度データの組み合わせを選択してもよい。
つぎに、多変量判別式作成部102hは、指定ストレス状態情報ファイル106dの所定の記憶領域に格納されているストレス状態情報に含まれるアミノ酸濃度データの全ての組み合わせが終了したか否かを判定し、判定結果が「終了」であった場合(ステップSB−24:Yes)には次のステップ(ステップSB−25)へ進み、判定結果が「終了」でなかった場合(ステップSB−24:No)にはステップSB−21へ戻る。なお、多変量判別式作成部102hは、予め設定した回数が終了したか否かを判定し、判定結果が「終了」であった場合には(ステップSB−24:Yes)次のステップ(ステップSB−25)へ進み、判定結果が「終了」でなかった場合(ステップSB−24:No)にはステップSB−21へ戻ってもよい。また、多変量判別式作成部102hは、ステップSB−23で選択したアミノ酸濃度データの組み合わせが、指定ストレス状態情報ファイル106dの所定の記憶領域に格納されているストレス状態情報に含まれるアミノ酸濃度データの組み合わせまたは前回のステップSB−23で選択したアミノ酸濃度データの組み合わせと同じであるか否かを判定し、判定結果が「同じ」であった場合(ステップSB−24:Yes)には次のステップ(ステップSB−25)へ進み、判定結果が「同じ」でなかった場合(ステップSB−24:No)にはステップSB−21へ戻ってもよい。また、多変量判別式作成部102hは、検証結果が具体的には各候補多変量判別式に関する評価値である場合には、当該評価値と各式作成手法に対応する所定の閾値との比較結果に基づいて、ステップSB−25へ進むかステップSB−21へ戻るかを判定してもよい。
ついで、多変量判別式作成部102hは、検証結果に基づいて、複数の候補多変量判別式の中から多変量判別式として採用する候補多変量判別式を選出することで多変量判別式を決定し、決定した多変量判別式(選出した候補多変量判別式)を多変量判別式ファイル106e4の所定の記憶領域に格納する(ステップSB−25)。ここで、ステップSB−25において、例えば、同じ式作成手法で作成した候補多変量判別式の中から最適なものを選出する場合と、すべての候補多変量判別式の中から最適なものを選出する場合とがある。
これにて、多変量判別式作成処理の説明を終了する。
健常群と問診結果から抑うつ病・大うつ病と判定されたストレス所見群の血液サンプルから、前述のアミノ酸分析法により血中アミノ酸濃度を測定した。図23に、各群毎(1:健常群、2:抑うつ病所見群、3:大うつ病所見群)のアミノ酸変数の分布を(縦軸:nmol/mL、図中のABAはα−ABAを、またCysはCystineを表す)箱ひげ図を示す。
ここで、抑うつ病の診断方法は、職業簡易ストレス調査票で抑うつ感を示す以下のaからfの項目に該当する問診項目の結果(少なくともa〜fの項目を1つ含む)から抑うつ病と判定する。
a 憂鬱だ。
b 何をするのも面倒だ。
c 物事に集中できない。
d 気分が晴れない。
e 仕事が手につかない。
f 悲しいと感じる。
また、大うつ病の診断方法は、以下に示すDSM−IVの問診9項目(AからI)のうち5項目以上(ただし、少なくとも、AまたはBの項目を含む)が同時に2週間以上持続する場合に大うつ病と判定する。
A ほとんど1日中気分が沈んでいる。
B 何に対しても興味がわかず、楽しめない。
C 食欲が低下、または体重の増減が激しい。
D 寝付けない。夜中や早朝に目が覚める。
E 動作や話し方が遅い、又はいらいらしたり落ち着きが無い。
F 疲れを感じたり、気力がわかない。
G 自分に価値が無い、または申し訳ないと感じる。
H 仕事や家事に集中したり、決断することができない。
I この世から消えてしまいたいと思うことがある。
まず、健常群と抑うつ病と判定された抑うつ病所見群の2群間の判別を目的に2群間のノンパラメトリック解析であるMann−Whitney検定を実施した。健常群に比べて抑うつ病所見群では、Lys(有意差確率P<0.05)が有意に減少し、アミノ酸変数Lysが健常群と抑うつ病所見群の2群間の判別能を持つことが判明した。
次に、健常群と大うつ病所見群の2群間の判別を目的に2群間のノンパラメトリック解析であるMann−Whitney検定を実施した。健常群に比べて大うつ病所見群では、Lys,His,ABA,Asn(有意差確率P<0.05)が有意に減少し、アミノ酸変数Lys,His,ABA,Asnが健常群と大うつ病所見群の2群間の判別能を持つことが判明した。
ここでは、多変量判別式の例として、ロジスティック回帰式、線形判別式を挙げ、実施例1で用いたサンプルデータから、健常群とストレス所見群(抑うつ病・大うつ病)の判別式を導出した。本実施の形態にかかる「健常群とストレス所見群(抑うつ病・大うつ病)の多変量判別式」とは、ロジスティック回帰分析、線形判別分析、マハラノビス距離、サポートベクターマシン、正準判別分析、決定木などの判別解析結果から得られる式を指す。
まず、健常群と抑うつ病所見群の2群判別性能を最大化する多変量判別式をロジスティック解析(Wald検定のステップワイズ法による変数選択)により探索し、判別式としてSer,Glu,Cys,Lysから構成されるロジスティック回帰式(アミノ酸変数、Ser,Glu,Cys,Lysの数係数と定数項は順に、−0.03、0.04、0.10、−0.02、2.52)が得られた。なお、このほかに同等の判別性能を有するロジスティック回帰式は複数得られた。それらを図24、図25、図26に示す。また、図24、図25、図26に示す式における各係数の値は、それを実数倍したものでもよい。
この判別式による健常群と抑うつ病所見群の2群判別に関して、ROC曲線(図27)のAUCによる評価を行い、0.796±0.038(95%信頼区間は0.721〜0.872)が得られ診断性能が高く有用な指標であることが判明した。
次に健常群と大うつ病所見群の2群判別性能を最大化する多変量判別式をロジスティック解析(Wald検定のステップワイズ法による変数選択)により探索し、判別式としてSer,Gln,Pro,ABA,Tyr,Trpから構成されるロジスティック回帰式(アミノ酸変数Ser,Gln,Pro,ABA,Tyr,Trpの数係数と定数項は順に、0.2、−0.1、0.4、−2.9、1.1、−1.7、46.8)が得られた。なお、このほかに同等の判別性能を有するロジスティック回帰式は複数得られた。それらを図28、図29、図30に示す。また、図28、図29、図30に示す式における各係数の値は、それを実数倍したものでもよい。
この判別式による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図31)のAUCによる評価を行い、0.997±0.005(95%信頼区間は0.987〜1.000)が得られ診断性能が高く有用な指標であることが判明した。
さらに、健常群と抑うつ病所見群の2群判別性能を最大化する多変量判別式を線形判別解析(ステップワイズ法による変数選択)により探索し、判別式としてSer,Glu,Cys,Lys(アミノ酸変数Ser,Glu,Cys,Lysの数係数と定数項は順に、−0.03、0.03、0.09、−0.02、2.47)から構成される線形判別式が得られた。なお、このほかに同等の判別性能を有する線形判別式は複数得られた。それらを図32、図33、図34に示す。また、図32、図33、図34に示す式における各係数の値は、それを実数倍したものでもよい。
判別式による健常群と抑うつ病所見群の2群判別に関して、ROC曲線(図35)のAUCによる評価を行い、0.798±0.038(95%信頼区間は0.723〜0.873)が得られ診断性能が高く有用な指標であることが判明した。
加えて、健常群と大うつ病所見群の2群判別性能を最大化する多変量判別式を線形判別解析(ステップワイズ法による変数選択)により探索し、判別式としてSer,Gln,Pro,ABA,Cys,Ile,Tyr,Trpから構成される線形判別式(アミノ酸変数Ser,Gln,Pro,ABA,Cys,Ile,Tyr,Trpの数係数と定数項は順に、−0.02、0.02、−0.01、0.11、−0.08、−0.04、−0.04、0.10、−5.77)が得られた。なお、このほかに同等の判別性能を有する線形判別式は複数得られた。それらを図36、図37、図38に示す。また、図36、図37、図38に示す式における各係数の値は、それを実数倍したものでもよい。
この判別式による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図39)のAUCによる評価を行い、0.979±0.014(95%信頼区間は0.952〜1.000)が得られ診断性能が高く有用な指標であることが判明した。
以上のことから、一般的に用いられている多変量解析法から、健常群とストレス所見群を高い精度で判別する式が導出可能であることが示された。また、健常群とストレス所見群の多変量判別式中の変数としては、アミノ酸変数に加えて、アミノ酸以外のメタボロームデータ(有機酸、脂肪酸、核酸など)や、トランスクリプトームデータ、プロテオームデータも用いることにより、判別能は向上すると言える。
実施例1で用いたサンプルデータを用いた。本出願人による国際出願である国際公開WO2004/052191号パンフレットに記載の方法を用いて、健常群と抑うつ病所見群の2群判別性能を最大化する指標を鋭意探索し、同等の性能を持つ複数の指標の中に指標式1が得られた。なお、このほかに指標式1と同等の判別性能を有する多変量判別式は複数得られた。それらを図40、図41、図42に示す。
指標式1:(Cys)/(Asn) + (Glu+Tyr+Phe)/(Tau+Lys)
指標式1による健常群と抑うつ病所見群の2群判別に関して、ROC曲線(図43)のAUCによる評価を行い、0.813±0.037(95%信頼区間は0.741〜0.885)が得られた。また指標式1による健常群と抑うつ病所見群の2群判別のカットオフ値について、抑うつ病所見群の有症率を0.36として最適なカットオフ値を求めると、カットオフ値が1.38となり、感度89%、特異度72%、正診率78%が得られ、診断性能が高く有用な指標であることが判明した。
次に健常群と大うつ病所見群の2群判別性能を最大化する指標を鋭意探索し、同等の性能を持つ複数の指標の中に指標式2、指標式3が得られた。なお、このほかに指標式3と同等の判別性能を有する多変量判別式は複数得られた。それらを図44、図45、図46に示す。
指標式2:(Phe)/(Tau+Thr) + (Glu+Ser)/(Val+Lys)
指標式3:(Gly)/(Gln) + (Glu+Pro+Ile)/(Orn+Trp+ABA)
指標式2による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図47)のAUCによる評価を行い、0.881±0.047(95%信頼区間は0.709〜0.974)が得られた。また指標式2による健常群と大うつ病所見群の2群判別のカットオフ値について、大うつ病所見群の有症率を0.08として最適なカットオフ値を求めると、カットオフ値が0.794となり、感度86%、特異度86%、正診率86%が得られ、診断性能が高く有用な指標であることが判明した。
指標式3による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図48)のAUCによる評価を行い、0.909±0.039(95%信頼区間は0.832〜0.985)が得られた。また指標式3による健常群と大うつ病所見群の2群判別のカットオフ値について、大うつ病所見群の有症率を0.08として最適なカットオフ値を求めると、カットオフ値が2.289となり、感度86%、特異度84%、正診率84%が得られ、診断性能が高く有用な指標であることが判明した。
実施例1で用いたサンプルデータを用いた。本出願人による国際出願である国際公開WO2004/052191号パンフレットに記載の方法を用いて、健常群と抑うつ病所見群の2群判別性能を最大化する指標を鋭意探索し、同等の性能を持つ複数の指標の中に指標式4が得られた。
指標式4:(Glu/Lys) + (−5.169)×(Met/AAA) + 1.6769×(Ile/Leu) + 5.5337×(Phe/LNAAs)
指標式4による健常群と抑うつ病所見群の2群判別に関して、ROC曲線(図49)のAUCによる評価を行い、0.827±0.035(95%信頼区間は0.758〜0.897)が得られた。また指標式4による健常群と抑うつ病所見群の2群判別のカットオフ値について、抑うつ病所見群の有症率を0.36として最適なカットオフ値を求めると、カットオフ値が0.82となり、感度74%、特異度71%、正診率72%が得られ、診断性能が高く有用な指標であることが判明した。
次に健常群と大うつ病所見群の2群判別性能を最大化する指標を鋭意探索し、同等の性能を持つ複数の指標の中に指標式5が得られた。
指標式5:(Asn) + (−2.8556)×(Ala/Trp) + 48.6047×(ABA/Tyr) + (−8.4761)×(Ser/Val)
指標式5による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図50)のAUCによる評価を行い、0.986±0.011(95%信頼区間は0.965〜1.000)が得られた。また指標式5による健常群と大うつ病所見群の2群判別のカットオフ値について、大うつ病所見群の有症率を0.08として最適なカットオフ値を求めると、カットオフ値が27.50となり、感度86%、特異度96%、正診率96%が得られ、診断性能が高く有用な指標であることが判明した。
精神研式ストレス質問紙の試験結果として得られたストレスレベルと相関するアミノ酸変数として、Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAs(有意確率P<0.05)が有意な相関(PEARSON相関)を示し、アミノ酸変数Phe,Leu,Ile,Val,Trp,Tyr,Met,BCAA,AAA,LNAAsが抑うつ病や大うつ病に至るまでのストレスレベルの評価能を持つことが判明した。
ここで、ストレスレベルを示す指標としては、精神研式ストレス質問紙を使用した。この指標は、心身症状、性格特性、タイプA特性、ストレッサーの4つの項目群から成り、各項目群を0点、1点、2点の3段階評価で点数化し集計したものがストレスレベルを示す。
実施例1で用いたサンプルデータを用いた。本出願人による国際出願である国際公開WO2004/052191号パンフレットに記載の方法を用いて、ストレスレベルの状態の診断性能を最大化する指標を鋭意探索し、同等の性能を持つ複数の指標の中に指標式6が得られた。
指標式6:(Phe) + (−21.158)×(Cit/Glu) + 6.9843×(Gly/His)
ここでは、多変量判別式の例として、ロジスティック回帰式、線形判別式を挙げ、健常群と大うつ病所見群の判別式を導出した。本実施の形態にかかる「健常群と大うつ病所見群の多変量判別式」とは、ロジスティック回帰分析、線形判別分析、マハラノビス距離、サポートベクターマシン、正準判別分析、決定木などの判別解析結果から得られる式を指す。
まず、健常群と大うつ病所見群の2群判別性能を最大化する多変量判別式をロジスティック解析(Wald検定のステップワイズ法による変数選択)により探索し、判別式としてGlu,Asn,Gln,Cit,ABAから構成されるロジスティック回帰式(アミノ酸変数Glu,Asn,Gln,Cit,ABAの数係数と定数項は順に、0.03、−0.17、0.01、0.11、−0.12、−0.38)が得られた。なお、このほかに同等の判別性能を有するロジスティック回帰式は複数得られた。それらを図51、図52、図53に示す。また、図51、図52、図53に示す式における各係数の値は、それを実数倍したものでもよい。
この判別式による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図54)のAUCによる評価を行い、0.812±0.055(95%信頼区間は0.703〜0.920)が得られ診断性能が高く有用な指標であることが判明した。
また、健常群と大うつ病所見群の2群判別性能を最大化する多変量判別式を線形判別解析(ステップワイズ法による変数選択)により探索し、判別式としてAsn,Gln,Tau,Cit,ABAから構成される線形判別式(アミノ酸変数Asn,Gln,Tau,Cit,ABAの数係数と定数項は順に、−0.11、0.005、0.04、0.08、−0.06、−0.74)が得られた。なお、このほかに同等の判別性能を有する線形判別式は複数得られた。それらを図55、図56、図57に示す。また、図55、図56、図57に示す式における各係数の値は、それを実数倍したものでもよい。
この判別式による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図58)のAUCによる評価を行い、0.780±0.059(95%信頼区間は0.665〜0.896)が得られ診断性能が高く有用な指標であることが判明した。
以上のことから、一般的に用いられている多変量解析法から、健常群と大うつ病所見群を高い精度で判別する式が導出可能であることが示された。また、健常群とストレス所見群の多変量判別式中の変数としては、アミノ酸変数に加えて、アミノ酸以外のメタボロームデータ(有機酸、脂肪酸、核酸など)や、トランスクリプトームデータ、プロテオームデータも用いることにより、判別能は向上すると言える。
実施例1で用いたサンプルデータを用いた。本出願人による国際出願である国際公開WO2004/052191号パンフレットに記載の方法を用いて、健常群と大うつ病所見群の2群判別性能を最大化する指標を鋭意探索し、同等の性能を持つ複数の指標の中に指標式7が得られた。なお、このほかに指標式7と同等の判別性能を有する多変量判別式は複数得られた。それらを図59、図60、図61に示す。
指標式7:(Glu+Gln)/(Val+Thr+Ser) + (Cit+Met)/(Asn)
指標式7による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図62)のAUCによる評価を行い、0.849±0.071(95%信頼区間は0.668〜0.947)が得られた。また指標式7による健常群と大うつ病所見群の2群判別のカットオフ値について、大うつ病所見群の有症率を0.5として最適なカットオフ値を求めると、カットオフ値が2.45となり、感度70%、特異度84%、正診率77%が得られ、診断性能が高く有用な指標であることが判明した。
次に健常群と大うつ病所見群の2群判別性能を最大化する指標を鋭意探索し、同等の性能を持つ複数の指標の中に指標式8が得られた。
指標式8:(Asn/Gln) + (−0.014363)×(Cit/ABA) +
(−0.045196)×(Glu/Arg) + (−0.80562)×(Met/BCAA)
指標式8による健常群と大うつ病所見群の2群判別に関して、ROC曲線(図63)のAUCによる評価を行い、0.873±0.046(95%信頼区間は0.783〜0.963)が得られた。また指標式8による健常群と大うつ病所見群の2群判別のカットオフ値について、大うつ病所見群の有症率を0.5として最適なカットオフ値を求めると、カットオフ値が−0.020となり、感度77%、特異度87%、正診率82%が得られ、診断性能が高く有用な指標であることが判明した。
以上のように、本発明にかかるストレスの評価方法、ならびにストレス評価装置、ストレス評価方法、ストレス評価システム、ストレス評価プログラムおよび記録媒体は、産業上の多くの分野、特に医薬品や食品、医療などの分野で広く実施することができ、特に、病態予測や疾病リスク予測やプロテオームやメタボローム解析などを行うバイオインフォマティクス分野において極めて有用である。

Claims (23)

  1. 血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む抑うつ病の状態を評価するための多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出ステップ
    を含むことを特徴とする抑うつ病の評価方法。
  2. 前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病の状態を評価する判別値基準評価ステップ
    をさらに含むことを特徴とする請求項1に記載の抑うつ病の評価方法。
  3. 前記判別値基準評価ステップは、前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病であるかを評価すること
    を特徴とする請求項2に記載の抑うつ病の評価方法。
  4. 多変量判別式の値である判別値に基づいて、評価対象につき、抑うつ病の状態を評価する判別値基準評価ステップを含み、
    前記判別値は、血液中のアミノ酸の濃度値に関する前記評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む前記抑うつ病の状態を評価するための前記多変量判別式に基づいて算出されたものであること、
    を特徴とする抑うつ病の評価方法。
  5. 血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値に基づいて、前記評価対象につき、抑うつ病の状態を評価する濃度値基準評価ステップ
    を含むことを特徴とする抑うつ病の評価方法。
  6. 制御手段を備えた抑うつ病評価装置であって、
    前記制御手段は、
    血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む抑うつ病の状態を評価するための多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出手段
    を備えたこと
    を特徴とする抑うつ病評価装置。
  7. 前記制御手段は、
    前記判別値算出手段で算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病の状態を評価する判別値基準評価手段
    をさらに備えたこと
    を特徴とする請求項6に記載の抑うつ病評価装置。
  8. 前記判別値基準評価手段は、前記判別値算出手段で算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病であるかを評価すること
    を特徴とする請求項7に記載の抑うつ病評価装置。
  9. 制御手段を備えた抑うつ病評価装置であって、
    前記制御手段は、
    多変量判別式の値である判別値に基づいて、評価対象につき、抑うつ病の状態を評価する判別値基準評価手段
    を備え、
    前記判別値は、血液中のアミノ酸の濃度値に関する前記評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む前記抑うつ病の状態を評価するための前記多変量判別式に基づいて算出されたものであること、
    を特徴とする抑うつ病評価装置。
  10. 制御手段を備えた抑うつ病評価装置であって、
    前記制御手段は、
    血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値に基づいて、前記評価対象につき、抑うつ病の状態を評価する濃度値基準評価手段
    を備えたこと
    を特徴とする抑うつ病評価装置。
  11. 制御手段を備えた情報処理装置で実行する抑うつ病評価方法であって、
    前記制御手段で実行する、
    血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む抑うつ病の状態を評価するための多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出ステップ
    を含むこと
    を特徴とする抑うつ病評価方法。
  12. 前記制御手段で実行する、
    前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病の状態を評価する判別値基準評価ステップ
    をさらに含むこと
    を特徴とする請求項11に記載の抑うつ病評価方法。
  13. 前記判別値基準評価ステップは、前記判別値算出ステップで算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病であるかを評価すること
    を特徴とする請求項12に記載の抑うつ病評価方法。
  14. 制御手段を備えた情報処理装置で実行する抑うつ病評価方法であって、
    前記制御手段で実行する、
    多変量判別式の値である判別値に基づいて、評価対象につき、抑うつ病の状態を評価する判別値基準評価ステップ
    を含み、
    前記判別値は、血液中のアミノ酸の濃度値に関する前記評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む前記抑うつ病の状態を評価するための前記多変量判別式に基づいて算出されたものであること、
    を特徴とする抑うつ病評価方法。
  15. 制御手段を備えた情報処理装置で実行する抑うつ病評価方法であって、
    前記制御手段で実行する、
    血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値に基づいて、前記評価対象につき、抑うつ病の状態を評価する濃度値基準評価ステップ
    を含むこと
    を特徴とする抑うつ病評価方法。
  16. 制御手段を備えた抑うつ病評価装置と、血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データを提供する端末装置とを、ネットワークを介して通信可能に接続して構成された抑うつ病評価システムであって、
    前記端末装置は、
    前記評価対象の前記アミノ酸濃度データを前記抑うつ病評価装置へ送信するアミノ酸濃度データ送信手段と、
    前記抑うつ病評価装置から送信された、多変量判別式の値である判別値または抑うつ病の状態に関する評価結果を受信する結果受信手段と
    を備え、
    前記抑うつ病評価装置の前記制御手段は、
    前記端末装置から送信された前記評価対象の前記アミノ酸濃度データを受信するアミノ酸濃度データ受信手段と、
    前記アミノ酸濃度データ受信手段で受信した前記評価対象の前記アミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む抑うつ病の状態を評価するための前記多変量判別式に基づいて、前記判別値を算出する判別値算出手段、または、
    前記アミノ酸濃度データ受信手段で受信した前記評価対象の前記アミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む前記抑うつ病の状態を評価するための多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出手段、および、前記判別値算出手段で算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病の状態を評価する判別値基準評価手段と、
    前記判別値算出手段で算出した前記判別値または前記判別値基準評価手段での前記評価対象の前記評価結果を前記端末装置へ送信する結果送信手段と、
    を備えたこと
    を特徴とする抑うつ病評価システム。
  17. 制御手段を備えた情報処理装置に実行させるための抑うつ病評価プログラムであって、
    前記制御手段に実行させるための、
    血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む抑うつ病の状態を評価するための多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出ステップ
    を含むこと
    を特徴とする抑うつ病評価プログラム。
  18. 制御手段を備えた情報処理装置に実行させるための抑うつ病評価プログラムであって、
    前記制御手段に実行させるための、
    多変量判別式の値である判別値に基づいて、評価対象につき、抑うつ病の状態を評価する判別値基準評価ステップ
    を含み、
    前記判別値は、血液中のアミノ酸の濃度値に関する前記評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む前記抑うつ病の状態を評価するための前記多変量判別式に基づいて算出されたものであること、
    を特徴とする抑うつ病評価プログラム。
  19. 制御手段を備えた情報処理装置に実行させるための抑うつ病評価プログラムであって、
    前記制御手段に実行させるための、
    血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値に基づいて、前記評価対象につき、抑うつ病の状態を評価する濃度値基準評価ステップ
    を含むこと
    を特徴とする抑うつ病評価プログラム。
  20. 請求項17から19のいずれか1つに記載の抑うつ病評価プログラムを記録したこと
    を特徴とするコンピュータ読み取り可能な記録媒体。
  21. 多変量判別式の値である判別値または抑うつ病の状態に関する評価結果を取得する結果取得手段を備え、
    前記判別値は、血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む抑うつ病の状態を評価するための前記多変量判別式に基づいて算出されたものである、または、
    前記評価結果は、血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む前記抑うつ病の状態を評価するための多変量判別式に基づいて算出された当該多変量判別式の値である判別値に基づいて、前記評価対象につき、前記抑うつ病の状態を評価した結果であること
    を特徴とする端末装置。
  22. 前記判別値を算出するまたは前記抑うつ病の状態を評価する抑うつ病評価装置とネットワークを介して通信可能に接続されており、
    前記結果取得手段は、前記抑うつ病評価装置から送信された前記判別値または前記評価結果を受信すること、
    を特徴とする請求項21に記載の端末装置。
  23. 血液中のアミノ酸の濃度値に関する評価対象のアミノ酸濃度データを提供する端末装置とネットワークを介して通信可能に接続された、制御手段を備えた抑うつ病評価装置であって、
    前記制御手段は、
    前記端末装置から送信された前記評価対象の前記アミノ酸濃度データを受信するアミノ酸濃度データ受信手段と、
    前記アミノ酸濃度データ受信手段で受信した前記評価対象の前記アミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む抑うつ病の状態を評価するための多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出手段、または、
    前記アミノ酸濃度データ受信手段で受信した前記評価対象の前記アミノ酸濃度データに含まれるLys,His,ABA,Asn,Phe,Leu,Ile,Val,Trp,Tyr,Metのうち少なくとも2つの前記濃度値、および、前記少なくとも2つの濃度値を変数として含む抑うつ病の状態を評価するための多変量判別式に基づいて、当該多変量判別式の値である判別値を算出する判別値算出手段、および、前記判別値算出手段で算出した前記判別値に基づいて、前記評価対象につき、前記抑うつ病の状態を評価する判別値基準評価手段と、
    前記判別値算出手段で算出した前記判別値または前記判別値基準評価手段での前記評価対象の評価結果を前記端末装置へ送信する結果送信手段と、
    を備えたこと
    を特徴とする抑うつ病評価装置。
JP2008527784A 2006-08-04 2007-08-02 抑うつ病の評価方法、抑うつ病評価装置、抑うつ病評価方法、抑うつ病評価システム、抑うつ病評価プログラム、記録媒体および端末装置 Active JP5856364B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008527784A JP5856364B2 (ja) 2006-08-04 2007-08-02 抑うつ病の評価方法、抑うつ病評価装置、抑うつ病評価方法、抑うつ病評価システム、抑うつ病評価プログラム、記録媒体および端末装置

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2006213918 2006-08-04
JP2006213918 2006-08-04
JP2008527784A JP5856364B2 (ja) 2006-08-04 2007-08-02 抑うつ病の評価方法、抑うつ病評価装置、抑うつ病評価方法、抑うつ病評価システム、抑うつ病評価プログラム、記録媒体および端末装置
PCT/JP2007/065139 WO2008016101A1 (fr) 2006-08-04 2007-08-02 Procédé d'évaluation de stress, appareil d'évaluation de stress, système d'évaluation de stress, programme d'évaluation de stress et support d'enregistrement associé

Publications (2)

Publication Number Publication Date
JPWO2008016101A1 JPWO2008016101A1 (ja) 2009-12-24
JP5856364B2 true JP5856364B2 (ja) 2016-02-09

Family

ID=38997275

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008527784A Active JP5856364B2 (ja) 2006-08-04 2007-08-02 抑うつ病の評価方法、抑うつ病評価装置、抑うつ病評価方法、抑うつ病評価システム、抑うつ病評価プログラム、記録媒体および端末装置

Country Status (4)

Country Link
US (1) US8673647B2 (ja)
EP (1) EP2053405A4 (ja)
JP (1) JP5856364B2 (ja)
WO (1) WO2008016101A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2794161C2 (ru) * 2017-10-10 2023-04-12 Сисейдо Компани, Лтд. Программа и система обработки информации

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2008015929A1 (ja) * 2006-08-04 2009-12-24 味の素株式会社 メタボリック・シンドロームの評価方法、メタボリック・シンドローム評価装置、メタボリック・シンドローム評価方法、メタボリック・シンドローム評価システム、メタボリック・シンドローム評価プログラムおよび記録媒体、ならびにメタボリック・シンドロームの予防・改善物質の探索方法
US20110245092A1 (en) * 2008-03-04 2011-10-06 John Bilello Diagnosing and monitoring depression disorders based on multiple serum biomarker panels
CA2740736A1 (en) * 2008-10-15 2010-04-22 Ridge Diagnostics, Inc. Human biomarker hypermapping for depressive disorders
CN102301234B (zh) * 2008-11-18 2015-06-17 里奇诊断学股份有限公司 针对重度抑郁疾病的代谢综合症状及hpa轴生物标志物
KR102204498B1 (ko) * 2013-04-09 2021-01-19 아지노모토 가부시키가이샤 생활습관병 지표의 평가 방법, 생활습관병 지표 평가 장치, 생활습관병 지표 평가 방법, 생활습관병 지표 평가 프로그램 제품, 생활습관병 지표 평가 시스템, 및 정보통신 단말장치
JP6598425B2 (ja) * 2014-02-27 2019-10-30 花王株式会社 疲労の評価方法
US9801553B2 (en) 2014-09-26 2017-10-31 Design Interactive, Inc. System, method, and computer program product for the real-time mobile evaluation of physiological stress
EP3779453A4 (en) * 2018-04-10 2022-01-05 Ajinomoto Co., Inc. METHOD, DEVICE AND PROGRAM FOR EVALUATING THE HEALTH AND / OR GROWTH STATUS OF LIVESTOCK
JP7132568B2 (ja) * 2018-05-17 2022-09-07 Cyberdyne株式会社 生体情報計測装置及び生体情報計測方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004052191A1 (ja) * 2002-12-09 2004-06-24 Ajinomoto Co., Inc. 生体状態情報処理装置、生体状態情報処理方法、生体状態情報管理システム、プログラム、および、記録媒体

Family Cites Families (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59215673A (ja) 1983-05-20 1984-12-05 Japan Storage Battery Co Ltd 液循環式電池用バイポーラ極板
DE69214628T2 (de) 1991-01-18 1997-05-22 Esa Inc Verfahren zur Diagnose von Störungen der Körperfunktionen unter Verwendung von biochemischen Pryfilen
JP3215924B2 (ja) 1992-01-24 2001-10-09 久信 菅野 唾液の粘性測定とストレス表示装置
JP3070346B2 (ja) 1993-08-09 2000-07-31 セイコーエプソン株式会社 ストレス評価装置および生理的年齢評価装置
JP3429043B2 (ja) 1993-12-24 2003-07-22 兼房 加藤 新規ストレスタンパク質p20、その精製方法、合成方法、抗体、測定方法及び免疫学的試薬
JP3445038B2 (ja) 1995-09-26 2003-09-08 サントリー株式会社 精神的ストレス性障害関連抗原の測定法
JP3667511B2 (ja) 1997-10-31 2005-07-06 ポーラ化成工業株式会社 ストレスの鑑別法
US7329489B2 (en) 2000-04-14 2008-02-12 Matabolon, Inc. Methods for drug discovery, disease treatment, and diagnosis using metabolomics
US20020009740A1 (en) 2000-04-14 2002-01-24 Rima Kaddurah-Daouk Methods for drug discovery, disease treatment, and diagnosis using metabolomics
JP3754611B2 (ja) 2000-10-03 2006-03-15 旭テクノグラス株式会社 ヒト老化マーカー及びストレスマーカーの検定方法
CA2399169A1 (en) 2001-09-07 2003-03-07 Queen's University At Kingston Diagnostic methods for determining susceptibility to convulsive conditions
JP3764663B2 (ja) 2001-09-14 2006-04-12 独立行政法人電子航法研究所 心身診断システム
JP2003274959A (ja) 2002-03-22 2003-09-30 Adgene Co Ltd ポリメラーゼ連鎖反応を応用した非特異核酸増幅法
JP2003274949A (ja) 2002-03-22 2003-09-30 Inst Of Physical & Chemical Res うつ病薬の標的遺伝子及びその利用
JP3978721B2 (ja) 2002-07-31 2007-09-19 株式会社 矢内原研究所 心的外傷後ストレス障害の判定法
JP2004208547A (ja) 2002-12-27 2004-07-29 Hitachi Ltd うつ病の評価方法
JP2005004398A (ja) 2003-06-11 2005-01-06 Hirokazu Tachikawa インターネットを利用した精神疾患の診断・予防システム
US7064104B2 (en) * 2003-06-13 2006-06-20 The Procter & Gamble Company Methods of managing the symptoms of premenstrual syndrome
ES2318301T3 (es) 2003-07-11 2009-05-01 KOENIG & BAUER AKTIENGESELLSCHAFT Rotativo de bobinas.
JP4442197B2 (ja) 2003-11-17 2010-03-31 ソニー株式会社 ストレス評価方法
US20090297494A1 (en) 2004-01-15 2009-12-03 Michel Cuenod Diagnostic and treatment of a mental disorder
JP3923507B2 (ja) 2004-02-17 2007-06-06 株式会社総合医科学研究所 疲労度評価装置、疲労度評価方法及びその利用
JP2005312435A (ja) 2004-03-29 2005-11-10 Kazuhito Rokutan うつ病の評価方法
JP2006000015A (ja) 2004-06-15 2006-01-05 Kazuhito Rokutan ストレスの評価方法
JPWO2006098192A1 (ja) 2005-03-16 2008-08-21 味の素株式会社 生体状態評価装置、生体状態評価方法、生体状態評価システム、生体状態評価プログラム、評価関数作成装置、評価関数作成方法、評価関数作成プログラムおよび記録媒体
DE602006014379D1 (de) 2005-04-06 2010-07-01 Klinik Und Poliklinik Fuer Psy Neurodegenerative marker für depression.
EP1887362A4 (en) * 2005-05-30 2010-01-27 Ajinomoto Kk HEPATHOPATHY EVALUATION APPARATUS, HEPATHOPATHY EVALUATION METHOD, HEPATHOPATHY EVALUATION SYSTEM, HEPATHOPATHY EVALUATION PROGRAM, AND RECORDING MEDIUM
WO2007100782A2 (en) 2006-02-28 2007-09-07 Metabolon, Inc. Biomarkers for amyotrophic lateral sclerosis and methods using the same
US8394354B2 (en) 2006-05-09 2013-03-12 Metabolon, Inc. Biomarkers for depression and methods using the same
JP2008203250A (ja) * 2007-01-24 2008-09-04 Ajinomoto Co Inc Ibdの評価方法、ならびにアミノ酸情報処理装置、アミノ酸情報処理方法、アミノ酸情報処理システム、アミノ酸情報処理プログラムおよび記録媒体
EP2172775A4 (en) * 2007-06-25 2010-12-01 Ajinomoto Kk METHOD OF EVALUATING VISCERAL GREASE ACCUMULATION
JPWO2009054351A1 (ja) * 2007-10-25 2011-03-03 味の素株式会社 生体状態評価装置、生体状態評価方法、生体状態評価システム、生体状態評価プログラムおよび記録媒体
JP5976987B2 (ja) * 2008-02-06 2016-08-24 味の素株式会社 胃癌の評価方法、ならびに胃癌評価装置、胃癌評価方法、胃癌評価システム、胃癌評価プログラム、記録媒体および情報通信端末装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004052191A1 (ja) * 2002-12-09 2004-06-24 Ajinomoto Co., Inc. 生体状態情報処理装置、生体状態情報処理方法、生体状態情報管理システム、プログラム、および、記録媒体

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN6012067145; 河北英詮: 'うつ病の血漿遊離アミノ酸' 東京医科大学雑誌 Vol.52, No.2, 1994, P.135-144 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2794161C2 (ru) * 2017-10-10 2023-04-12 Сисейдо Компани, Лтд. Программа и система обработки информации

Also Published As

Publication number Publication date
EP2053405A1 (en) 2009-04-29
JPWO2008016101A1 (ja) 2009-12-24
WO2008016101A1 (fr) 2008-02-07
US20090234586A1 (en) 2009-09-17
EP2053405A4 (en) 2009-11-11
US8673647B2 (en) 2014-03-18

Similar Documents

Publication Publication Date Title
JP5856364B2 (ja) 抑うつ病の評価方法、抑うつ病評価装置、抑うつ病評価方法、抑うつ病評価システム、抑うつ病評価プログラム、記録媒体および端末装置
JP5746811B2 (ja) 大腸癌の評価方法、ならびに大腸癌評価装置、大腸癌評価方法、大腸癌評価システム、大腸癌評価プログラムおよび記録媒体
JP5976987B2 (ja) 胃癌の評価方法、ならびに胃癌評価装置、胃癌評価方法、胃癌評価システム、胃癌評価プログラム、記録媒体および情報通信端末装置
JP5470848B2 (ja) 肺癌の評価方法、肺癌評価装置、肺癌評価方法、肺癌評価システム、肺癌評価プログラム、記録媒体、および、情報通信端末装置
JP5817527B2 (ja) 肥満の評価方法、肥満評価装置、肥満評価方法、肥満評価システム、肥満評価プログラム、記録媒体、および情報通信端末装置
JP5754137B2 (ja) 前立腺癌の評価方法、前立腺癌評価装置、前立腺癌評価方法、前立腺癌評価プログラム、記録媒体、前立腺癌評価システム、および情報通信端末装置
JP2016029398A (ja) 内臓脂肪面積の評価方法、内臓脂肪面積評価装置、内臓脂肪面積評価方法、内臓脂肪面積評価プログラム、記録媒体、内臓脂肪面積評価システム、及び端末装置
JPWO2008015929A1 (ja) メタボリック・シンドロームの評価方法、メタボリック・シンドローム評価装置、メタボリック・シンドローム評価方法、メタボリック・シンドローム評価システム、メタボリック・シンドローム評価プログラムおよび記録媒体、ならびにメタボリック・シンドロームの予防・改善物質の探索方法
WO2014168125A1 (ja) 生活習慣病指標の評価方法、生活習慣病指標評価装置、生活習慣病指標評価方法、生活習慣病指標評価プログラム、生活習慣病指標評価システム、および情報通信端末装置
WO2013011919A1 (ja) Nashの評価方法、nash評価装置、nash評価方法、nash評価プログラム、nash評価システム、情報通信端末装置、およびnashの予防・改善物質の探索方法
JPWO2009054351A1 (ja) 生体状態評価装置、生体状態評価方法、生体状態評価システム、生体状態評価プログラムおよび記録媒体
JP2014002174A (ja) 乳癌の評価方法、ならびに乳癌評価装置、乳癌評価方法、乳癌評価システム、乳癌評価プログラムおよび記録媒体
WO2016056631A1 (ja) 評価方法、評価装置、評価プログラム、評価システム、及び端末装置
JP6614215B2 (ja) 取得方法、算出方法、脳卒中評価装置、算出装置、脳卒中評価プログラム、算出プログラムおよび脳卒中評価システム
JP6002117B2 (ja) 取得方法、評価装置、評価プログラム、記録媒体および評価システム
JP6260275B2 (ja) 脂肪肝の評価方法、脂肪肝評価装置、脂肪肝評価方法、脂肪肝評価プログラム、脂肪肝評価システム、および端末装置
JPWO2009054350A1 (ja) 耐糖能異常の評価方法、耐糖能異常評価装置、耐糖能異常評価方法、耐糖能異常評価システム、耐糖能異常評価プログラムおよび記録媒体、ならびに耐糖能異常の予防・改善物質の探索方法
WO2013115283A1 (ja) 早期腎症の評価方法、早期腎症評価装置、早期腎症評価方法、早期腎症評価プログラム、早期腎症評価システムおよび情報通信端末装置
JP6270257B2 (ja) 取得方法、心血管イベント評価装置、心血管イベント評価プログラム、心血管イベント評価システムおよび端末装置
JP2011080957A (ja) 線維筋痛症とうつ病の評価方法、線維筋痛症・うつ病評価装置および線維筋痛症・うつ病評価方法
WO2013005790A1 (ja) 生体酸化の評価方法、生体酸化評価装置、生体酸化評価方法、生体酸化評価プログラム、生体酸化評価システム、情報通信端末装置、および生体酸化の予防・改善物質の探索方法
JP6886241B2 (ja) 骨格筋面積の評価方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100713

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20121225

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130225

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20131203

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140228

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20140307

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20140509

A912 Re-examination (zenchi) completed and case transferred to appeal board

Free format text: JAPANESE INTERMEDIATE CODE: A912

Effective date: 20140704

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150624

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150918

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20151211

R150 Certificate of patent or registration of utility model

Ref document number: 5856364

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250