(57)【要約】      （修正有）【課題】健常なヒトの日々の生活や労働で生じる疲労を客観的に評価できる新たな評価方法及び評価システムを提供する。【解決手段】被験者の睡眠中の心拍間隔から、周波数解析によって高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝を算出する工程と、該ＨＦと、［ＨＦ／（ＬＦ＋ＨＦ）］について、被験者の睡眠の第１〜第３の睡眠周期毎のＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］の平均値を算出する工程とを含む、疲労の評価方法。【選択図】図１
2015109888.tif 000002
【特許請求の範囲】
【請求項１】
被験者の睡眠中の心拍間隔から、周波数解析によって高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝を算出する工程と、
該ＨＦと、［ＨＦ／（ＬＦ＋ＨＦ）］について、被験者の睡眠の第１〜第３の睡眠周期毎のＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］の平均値を算出する工程と、
次の（ａ）〜（ｈ）から選択される１又は２以上の工程と、
（ａ）第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第２睡眠周期）｝を算出する工程、
（ｂ）第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の差｛△ＨＦ（第２睡眠周期−第１睡眠周期）｝を算出する工程、
（ｃ）第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する工程、
（ｄ）第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）｝を算出する工程、
（ｅ）第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の差｛△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）｝を算出する工程、
（ｆ）第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の比｛［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）｝を算出する工程、
（ｇ）第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第１睡眠周期）｝を算出する工程、
（ｈ）第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する工程、
を含む、疲労の評価方法。
【請求項２】
工程（ａ）〜（ｈ）は、前記△ＨＦ（第３睡眠周期−第２睡眠周期）、△ＨＦ（第２睡眠周期−第１睡眠周期）、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）、［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）、△ＨＦ（第３睡眠周期−第１睡眠周期）、又はＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）を、予め設定した各々の所定の閾値と比較し、その結果に基づいて疲労度が判定される工程を含む、請求項１記載の疲労の評価方法。
【請求項３】
次の１〜５４いずれかに該当する場合に疲労度が高いと判定される工程を含む、請求項１又は２に記載の疲労の評価方法。
【表１】
2015109888.tif 000014
【請求項４】
被験者の睡眠中の心拍間隔から、周波数解析によって高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝を算出する手段と、
該ＨＦと、［ＨＦ／（ＬＦ＋ＨＦ）］について、被験者の睡眠の第１〜第３の睡眠周期毎のＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］の平均値を算出する手段と、
次の（ｉ）〜(vii)から選択される１又は２以上の手段と、
（ｉ）第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第２睡眠周期）｝を算出する手段、
(ii)第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の差｛△ＨＦ（第２睡眠周期−第１睡眠周期）｝を算出する手段、
(iii)第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する手段、
(iv)第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）｝を算出する手段、
（ｖ）第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の差｛△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）｝を算出する手段、
(vi)第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の比｛［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）｝を算出する手段、
(vii)第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第１睡眠周期）｝を算出する手段、
(vii)第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する手段、
を備える、疲労の評価システム。
【請求項５】
更に、前記△ＨＦ（第３睡眠周期−第２睡眠周期）、△ＨＦ（第２睡眠周期−第１睡眠周期）、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）、［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）、△ＨＦ（第３睡眠周期−第１睡眠周期）、又はＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）に関する所定の閾値を記憶する手段を備え、前記（ｉ）〜(vii)の手段は、前記△ＨＦ（第３睡眠周期−第２睡眠周期）、△ＨＦ（第２睡眠周期−第１睡眠周期）、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）、［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）、△ＨＦ（第３睡眠周期−第１睡眠周期）、又はＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）を、前記記憶手段に記憶した各々の所定の閾値と比較し、その結果に基づいて疲労度を判定する手段を備える、請求項４記載の疲労の評価システム。
【請求項６】
次の１〜５４いずれかに該当する場合に疲労度が高いと判定する手段を備える、請求項４又は５に記載の疲労の評価システム。
【表１】
2015109888.tif 000015
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、疲労の評価方法及び疲労の評価システムに関する。
【背景技術】
【０００２】
疲労は、現代の社会に生きる人が日常向き合っている現象である。我が国で疲労感を自覚しているヒトの割合は、就労人口の約６０％以上とも云われている。
このような状況の中、疲労の指標化が望まれ、疲労を非侵襲的な手法で客観的に評価する方法の探索研究が精力的になされてきた。
【０００３】
例えば、非特許文献１では、覚醒時平均活動量や睡眠時間等の睡眠関連指標から三つの異なる線形／非線形分析を行って慢性疲労病態を判別する方法や、指尖加速度脈波から得られた心拍時系列データによる自律神経機能解析で低周波成分（ＬＦ）と高周波成分（ＨＦ）を求め慢性疲労症候群患者の疲労を評価する方法が報告されている。
また、特許文献１では、被験者の交感神経の働き値（ＬＦ値）と、被験者の副交感神経の働き値（ＨＦ値）からＬＦ／ＨＦ値を計算してＬＦ／ＨＦ値を得、得られたＬＦ／ＨＦ値と疲労度判定基準値とを対比して、被験者の客観的疲労度を判定する、被験者の疲労度の判定処理を行う方法が報告されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−２０１１１３号公報
【非特許文献】
【０００５】
【非特許文献１】厚生労働科学研究費補助金（障害者対策総合研究事業）（神経・筋疾患分野）、「自律神経機能異常を伴い慢性的な疲労を訴える患者に対する客観的な疲労診断法の確立と慢性疲労診断指針の作成  平成２１〜２３年度総合研究報告書」ｐ．２０−２８、［ｏｎｌｉｎｅ］、平成２４年（２０１２年）３月、［平成２５年８月２６日検索］、インターネット＜ＵＲＬ：ｈｔｔｐ：／／ｗｗｗ．ｆｕｋｓｉ−ｋａｇｋ−ｕ．ａｃ．ｊｐ／ｇｕｉｄｅ／ｅｆｆｏｒｔｓ／ｒｅｓｅａｒｃｈ／ｋｕｒａｔｓｕｎｅ／ｈ２３／ｐｄｆ／ｈ２３ｈｏｕｋｏｋｕ．ｐｄｆ＞
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来の技術は、半年以上続く慢性的な疲労を抱える者、すなわち慢性疲労症候群の病的な疲労の抽出を目的とするもので、健常な労働者の日々の疲労を評価しようとしてもできないか、評価できても感度が低く偽陰性の確率が高いといった問題がある。
すなわち、今までのところ、健常なヒトの日々の生活や労働で生じる疲労を客観的に評価できる方法は確立されていなかった。
そこで、健常なヒトの疲労を客観的に評価する評価方法及び評価システムを提供するのが本発明の課題である。
【課題を解決するための手段】
【０００７】
本発明者らは、睡眠中の自律神経活動動態を知るために心拍の変化を解析していたところ、疲労の有無によって、睡眠周期毎の高周波数成分（ＨＦ）の変化と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比の変化に差があることを見出し、当該変化の違いを指標とすれば、健常者の日常生活で生じる疲労であっても客観的に評価できることに想到し本発明を完成した。
【０００８】
すなわち、本発明は、被験者の睡眠中の心拍間隔から、周波数解析によって高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝を算出する工程と、
該ＨＦと、［ＨＦ／（ＬＦ＋ＨＦ）］について、被験者の睡眠の第１〜第３の睡眠周期毎のＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］の平均値を算出する工程と、
次の（ａ）〜（ｈ）から選択される１又は２以上の工程と、
（ａ）第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第２睡眠周期）｝を算出する工程、
（ｂ）第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の差｛△ＨＦ（第２睡眠周期−第１睡眠周期）｝を算出する工程、
（ｃ）第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する工程、
（ｄ）第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）｝を算出する工程、
（ｅ）第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の差｛△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）｝を算出する工程、
（ｆ）第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の比｛［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）｝を算出する工程、
（ｇ）第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第１睡眠周期）｝を算出する工程、
（ｈ）第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する工程、
を含む、疲労の評価方法を提供するものである。
また、本発明は、被験者の睡眠中の心拍間隔から、周波数解析によって高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝を算出する手段と、
該ＨＦと、［ＨＦ／（ＬＦ＋ＨＦ）］について、被験者の睡眠の第１〜第３の睡眠周期毎のＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］の平均値を算出する手段と、
次の（ｉ）〜(vii)から選択される１又は２以上の手段と、
（ｉ）第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第２睡眠周期）｝を算出する手段、
(ii)第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の差｛△ＨＦ（第２睡眠周期−第１睡眠周期）｝を算出する手段、
(iii)第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する手段、
(iv)第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）｝を算出する手段、
（ｖ）第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の差｛△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）｝を算出する手段、
(vi)第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の比｛［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）｝を算出する手段、
(vii)第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第１睡眠周期）｝を算出する手段、
(vii)第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する手段、
を備える、疲労の評価システムを提供するものである。
【発明の効果】
【０００９】
本発明によれば、健常者の日常の生活や労働で生じる疲労を客観的に評価することができる。これにより、疲労度を的確に捉え、適切な人材を選択することや、状態に応じたアドバイスを与えることが可能である。
【図面の簡単な説明】
【００１０】
【図１】本発明の疲労の評価方法の概略フローチャートである。
【図２】本発明の疲労の評価システムの概略構成図である。
【図３】疲労の有無による１睡眠周期毎の高周波数成分（ＨＦ）の平均値の変化と［ＨＦ／（ＬＦ＋ＨＦ）］％の平均値の変化を示す図である。
【発明を実施するための形態】
【００１１】
以下、本発明の実施の形態を図面と共に説明する。
図１は、本発明の疲労の評価方法を示す概略フローチャートである。また、図２は、当該疲労の評価方法を実施する疲労の評価システムの一例の概略構成図である。
この疲労の評価システムは、被験者の睡眠中の心拍間隔から、周波数解析によって高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝を算出する手段１の一つとして、心拍間隔を測定する測定計６と、データを入力・周波数解析を行う解析装置８を備えている。睡眠周期を、心拍間隔を測定する測定計６を用いずに測定する場合は、睡眠周期の測定計７を備えていてもよい。
また、被験者の睡眠の第１〜第３の睡眠周期毎の自律神経活動解析を行って、ＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］の平均値を算出する手段２と、
次の（ｉ）〜(vii)から選択される１又は２以上の手段により、疲労を評価する手段３を備えており、これに算出・評価結果等を出力・記録する装置５が接続されている。手段３は、後述する閾値を記憶する手段４を備えていてもよい。
（ｉ）第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第２睡眠周期）｝を算出する手段、
(ii)第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の差｛△ＨＦ（第２睡眠周期−第１睡眠周期）｝を算出する手段、
(iii)第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する手段、
(iv)第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）｝を算出する手段、
（ｖ）第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の差｛△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）｝を算出する手段、
(vi)第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の比｛［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）｝を算出する手段、
(vii)第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第１睡眠周期）｝を算出する手段、
(vii)第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）｝を算出する手段。
【００１２】
なお、本発明において、疲労とは、一般的に、肉体的及び精神的活動、又は疾病によって生じた独特の不快感と休養の願望を伴う身体の活動能力の減退状態を云い、疲労度とは該疲労の度合いをいう。
疲労は、前記のとおり、病的な疲労とそれ以外の生理的な疲労に分類されるが、本発明における疲労は病的な疲労ではなく、生理的な、肉体的及び精神的な疲労を含む複合疲労であることが好ましい。
すなわち、本発明における被験者は、健常者、更に健常労働者、更に健常なディスクワーカーであることが好ましい。
【００１３】
手段１では、被験者の睡眠中の心拍間隔から、周波数解析によって高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝を算出する。低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比はパーセンテージで表すことが可能である。
心拍間隔は、心臓の拍動１拍ごとの間隔で、被験者の心拍を表す心拍信号から振幅ピークを検出することで取得できる。心拍信号としては、例えば、心電信号、脈波信号等が挙げられ、心電計、脈波計等によって測定される。
【００１４】
心拍間隔の周波数解析としては、従来公知の高速フーリエ変換法、最大エントロピー法、ウェーブレット法等が利用できる。例えば、所定時間における心拍間隔の変動をフーリエ変換して得られるパワースペクトルから、特定周波数よりも高い高周波帯域の積分値、すなわち高周波数成分（ＨＦ）と、特定周波数よりも低い低周波帯域の積分値、すなわち低周波数成分（ＬＦ）を算出する。
ここで、高周波帯域は０．１５〜０．４０Ｈｚ、低周波帯域は０．０４〜０．１５Ｈｚまでである。
【００１５】
次に手段２では、この高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝について、被験者の睡眠の第１〜第３の睡眠周期毎のＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］の平均値を算出する。
正常な睡眠では、覚醒期から睡眠状態に入ると、睡眠周期（ウルトラディアンリズム）と呼ばれる典型的には９０分、一般的には６０分〜１２０分の周期で睡眠状態の遷移が一晩に４周期〜５周期繰り返される（睡眠学、２００９年、ｐ．１７１、日本睡眠学会編、朝倉書店東京）。睡眠周期では、レム睡眠とノンレム睡眠とが交互に出現する。本発明では、被験者の一晩の睡眠周期（第１睡眠周期〜第４乃至第５睡眠周期）のうち、第１、第２及び第３の睡眠周期毎のＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］の平均値を求める。
【００１６】
睡眠周期は、従来公知の睡眠ポリグラフ検査（ＰＳＧ、ｐｏｌｙｓｏｍｎｏｇｒａｐｈｙ）による睡眠ポリグラフから算出できる。また、心拍間隔（Ｒ−Ｒ間隔）の周期がレム睡眠とノンレム睡眠の週間と相関すること（Ｎｅｕｒｏｓｃｉ  Ｌｅｔｔｅｒ，１９９７，２２９，ｐ．１７３−１７６）を利用して、心拍間隔の１分間のＲ−Ｒ間隔の平均、又は１分間の心拍数（回／分）の自己相関係数から算出することもできる。
また、より簡易的に、睡眠周期は典型的には９０分周期で繰り返されることから、睡眠時間を９０分間間隔に分割して求めてもよい。この場合、周波数解析に先だってデータを９０分間間隔に分割して、その後手段１により、高周波数成分（ＨＦ）と、低周波数成分（ＬＦ）と、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比｛［ＨＦ／（ＬＦ＋ＨＦ）］｝を算出することもできる。  
【００１７】
次に手段３では、第１〜第３睡眠周期間における、高周波数成分（ＨＦ）の平均値の差、ＨＦの平均値の比、低周波数成分（ＬＦ）及び高周波数成分（ＨＦ）の総和に対する高周波数成分（ＨＦ）の比［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の差、又は［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の比を求め、疲労を評価する。
具体的には、
第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第２睡眠周期）｝、
第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の差｛△ＨＦ（第２睡眠周期−第１睡眠周期）｝、
第２睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）｝、
第１睡眠周期のＨＦの平均値と第２睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）｝、
第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の差｛△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）｝、
第２睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］の平均値の比｛［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）｝、
第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の差｛△ＨＦ（第３睡眠周期−第１睡眠周期）｝、
第１睡眠周期のＨＦの平均値と第３睡眠周期のＨＦの平均値の比｛ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）｝
の１又は２以上を求める。
【００１８】
図３に示すように、被験者の疲労度が高いと、第１〜第３睡眠周期毎、ことに第２睡眠周期間と第３睡眠周期の間でＨＦの平均値が大きく変化することが見出された。したがって、この変化は、被験者の疲労を反映する指標となり得る。
そのため、前記△ＨＦ（第３睡眠周期−第２睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）及び△ＨＦ（第３睡眠周期−第１睡眠周期）については、その値が大きいと被験者の疲労度は高いと判定され、他方、その値が小さいと被験者の疲労度は低いと判定される。
また、△ＨＦ（第２睡眠周期−第１睡眠周期）、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）、［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）については、その値が小さいと被験者の疲労度は高いと判定され、他方、その値が大きいと被験者の疲労度は低いと判定される。
【００１９】
本発明では、前記△ＨＦ（第３睡眠周期−第２睡眠周期）、△ＨＦ（第２睡眠周期−第１睡眠周期）、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）、［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）、△ＨＦ（第３睡眠周期−第１睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）のうち、好ましくは４つ以上、より好ましくは５つ以上を求め、疲労の評価に用いるのが好ましい。更に好ましくは、△ＨＦ（第３睡眠周期−第２睡眠周期）、△ＨＦ（第２睡眠周期−第１睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）を組み合わせた評価が好ましい。
【００２０】
更に、評価システムでは、前記△ＨＦ（第３睡眠周期−第２睡眠周期）、△ＨＦ（第２睡眠周期−第１睡眠周期）、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）、［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）、△ＨＦ（第３睡眠周期−第１睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）に関する所定の閾値を記憶する手段４を備えるのが好ましい。
すなわち、手段３は、算出した前記△ＨＦ（第３睡眠周期−第２睡眠周期）、△ＨＦ（第２睡眠周期−第１睡眠周期）、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期−第２睡眠周期）、［ＨＦ／（ＬＦ＋ＨＦ）］（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］（第３睡眠周期）、△ＨＦ（第３睡眠周期−第１睡眠周期）、又はＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）を、手段４において記憶した各々の所定の閾値と比較し、その結果に基づいて被験者の疲労度を判定する手段を備えることが好ましい。
閾値は、この値を超えると或いは満たないと疲労度が高い、疲労が蓄積していると判断される基準となる値である。
【００２１】
閾値は、例えば、後述するＶＡＳ評価で「疲労あり」の被験者、「疲労なし」の被験者毎に、ＨＦの平均値、［ＨＦ／（ＬＦ＋ＨＦ）］の平均値から、ＲＯＣ曲線（Ｒｅｃｅｉｖｅｒ  Ｏｐｅｒａｔｉｎｇ  Ｃｈａｒａｃｔｅｒｉｓｔｉｃ  ｃｕｒｖｅ、受信者動作特性曲線）を作成し予め設定することができる。
この際、判定の信頼性を図る尺度である感度・特異度・有効度をそれぞれ任意に設定し、その割合に応じた閾値を設定することができる。
感度は、疲労感がある被験者が陽性となる割合（真陽性率）を意味し、特異度は疲労感のない被験者が陰性となる割合（真陰性率）を意味する。特異度が高くなるように閾値を設定すると偽陰性が多くなり、感度が低くなる。一方、感度が高くなるように設定すると偽陽性が多くなり、特異度は低くなる。そこで、閾値は、目的に応じて決定される。
【００２２】
閾値は、例えば、疲労度が低いことが明らかな被験者の場合は、△ＨＦ（第３睡眠周期−第２睡眠周期）の閾値は０．２以上、△ＨＦ（第２睡眠周期−第１睡眠周期）の閾値は−０．２以下、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）の閾値は１以下、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）の閾値は１以上、△［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）−［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）の閾値は−１．１以上、［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）の閾値は１以下、△ＨＦ（第３睡眠周期−第１睡眠周期）の閾値は−０．１以上、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）の閾値は1以下で設定することが好ましい。
より好ましくは、△ＨＦ（第３睡眠周期−第２睡眠周期）の閾値は２．３以上、△ＨＦ（第２睡眠周期−第１睡眠周期）の閾値は−０．２２以下、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）の閾値は０．９２以下、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）の閾値は１．０１以上、△［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）−［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）の閾値は−１．１１以上、［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）の閾値は１．０２以下、△ＨＦ（第３睡眠周期−第１睡眠周期）の閾値は２．２以上、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）の閾値は０．９１以下で設定することが好ましい。
さらに好ましくは、△ＨＦ（第３睡眠周期−第２睡眠周期）の閾値は１２以上、△ＨＦ（第２睡眠周期−第１睡眠周期）の閾値は−６以下、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）の閾値は０．７以下、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）の閾値は１．２以上、△［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）−［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）の閾値は６以上、［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）の閾値は０．８以下、△ＨＦ（第３睡眠周期−第１睡眠周期）の閾値は−０．１２以上、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）の閾値は１以下で設定することがより好ましい。
【００２３】
また、疲労度が高いことが明らかな被験者の場合は、△ＨＦ（第３睡眠周期−第２睡眠周期）の閾値は０．２以上、△ＨＦ（第２睡眠周期−第１睡眠周期）の閾値は−０．２以下、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）の閾値は１以下、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）の閾値は１以上、△［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）−［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）の閾値は−１．１以上、［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）の閾値は１以下、△ＨＦ（第３睡眠周期−第１睡眠周期）の閾値は−０．１以上、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）の閾値は1以下で設定することが好ましい。
より好ましくは、△ＨＦ（第３睡眠周期−第２睡眠周期）の閾値は−１．３２以上、△ＨＦ（第２睡眠周期−第１睡眠周期）の閾値は０．２９以下、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）の閾値は１．０５以下、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）の閾値は０．９９以上、△［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）−［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）の閾値は−２．２６以上、［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）の閾値は１．０５以下、△ＨＦ（第３睡眠周期−第１睡眠周期）の閾値は−０．２１以上、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）の閾値は１以下で設定することが好ましい。
さらに好ましくは、△ＨＦ（第３睡眠周期−第２睡眠周期）の閾値は−３．６以上、△ＨＦ（第２睡眠周期−第１睡眠周期）の閾値は７以下、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）の閾値は１．２以下、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）の閾値は０．８以上、△［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）−［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）の閾値は−１０以上、［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）の閾値は１．２以下、△ＨＦ（第３睡眠周期−第１睡眠周期）の閾値は−６以上、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）の閾値は１．２５以下で設定することがより好ましい。
【００２４】
本発明では、手段３で、各項目について次の表１に示す閾値を記憶する手段４を備え、次の１〜５４いずれかに該当する場合に疲労度が高いと判定する手段を備えることが好ましい。
【００２５】
【表１】
2015109888.tif 000003
【００２６】
後述する実施例５に示すように、前記△ＨＦ（第３睡眠周期−第２睡眠周期）、△ＨＦ（第２睡眠周期−第１睡眠周期）、ＨＦ（第２睡眠周期）／ＨＦ（第３睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第２睡眠周期）、△［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期−第２睡眠周期）、［ＨＦ／（ＬＦ＋ＨＦ）］％（第２睡眠周期）／［ＨＦ／（ＬＦ＋ＨＦ）］％（第３睡眠周期）、△ＨＦ（第３睡眠周期−第１睡眠周期）、ＨＦ（第１睡眠周期）／ＨＦ（第３睡眠周期）の閾値を、例えば表１の４６〜５３のとおり設定すると、感度が９１〜９６％となり、偽陰性が極めて少なくなるため、疲労度が高い被験者を確実に抽出したい場合に有効な評価方法となる。
このように、本発明の評価方法を実施することで、健常者の日常の生活や労働で生じる疲労を客観的に評価することが可能である。
【実施例】
【００２７】
発明評価用データの収集
〔対象者〕
週５日（月曜日〜金曜日勤務、週休２日（土日））の日勤労働者（標準労働時間８：３０〜１７：００、１時間の昼休み含む）のうち４５名の健常男性労働者（２５〜５３歳）を対象として、連続７日間の就寝時の主観的な疲労と、睡眠中の第１〜第３睡眠周期までの心拍変動を測定した。のべ２１７日分のデータを解析した。
【００２８】
〔疲労評価〕
１．ＶＡＳ評価
就寝時の主観的な疲労の評価は、日本疲労学会の抗疲労臨床評価ガイドライン第５版（ｈｔｔｐ：／／ｗｗｗ．ｈｉｒｏｕｇａｋｋａｉ．ｃｏｍ／ｇｕｉｄｅｌｉｎｅ．ｐｄｆ）の疲労感ＶＡＳ（Ｖｉｓｕａｌ  Ａｎａｌｏｇｕｅ  Ｓｃａｌｅ）検査方法（ｈｔｔｐ：／／ｗｗｗ．ｈｉｒｏｕｇａｋｋａｉ．ｃｏｍ／ＶＡＳ．ｐｄｆ）を利用した。
疲労感ＶＡＳ検査は、１００ｍｍの直線の左端（０ｍｍ）を「疲れを全く感じない」、右端（１００ｍｍ）を「何もできないほど疲れきっている」感覚として、今、感じている疲労感を、直線の左右両端に示した感覚を参考に１００ｍｍの直線上に示してもらう方法である。
【００２９】
疲労の有無の判定は、疲労感ＶＡＳ（Ｖｉｓｕａｌ  Ａｎａｌｏｇｕｅ  Ｓｃａｌｅ）検査を３分位に分割し、下位３分の１（平均±標準偏差２７．８±１０ｍｍ、最小値〜最大値２〜４４ｍｍ）を「疲労なし」、上位３分の１（平均±標準偏差７２．７±６．９ｍｍ、最小値〜最大値６４〜９４ｍｍ）を「疲労あり」とした。なお、中位３分の１は平均±標準偏差５３．５±５．３ｍｍ、最小値〜最大値４５〜６３ｍｍである。
【００３０】
２．睡眠中の心拍変動の測定
睡眠時の心拍変動は、心臓の拍動とともに発生する指尖容積脈波から脈波間隔（Ｒ−Ｒ間隔）を測定し（東芝体動計ＮＥＭ−Ｔ１、（株）パラマテック）、１分間の脈波間隔を高速フーリエ変化で周波数解析を行った（東芝汎用睡眠解析ソフト「ねむリズム」ＮＥＭ−ＳＳ１、東芝ホームアプライアンス（株））。周波数成分の０．０４〜０．１５Ｈｚの成分を低周波数成分（ＬＦ）とし、０．１５から０．４Ｈｚの成分を高周波数成分（ＨＦ）とし、領域内のピーク値を各成分の１分間の値とした。また、ＬＦとＨＦの総和に対するＨＦのパーセンテージ［ＨＦ／（ＬＦ＋ＨＦ）］％を算出した。
【００３１】
睡眠周期は、１分間の心拍数、又は１分間の脈波間隔（Ｒ−Ｒ間隔）の平均（ｍｓｅｃ）の自己相関係数から算出した。具体的には、次式（１）から自己相関係数を計算し、自己相関係数が低下、増加し最大値となった時間を各睡眠周期の終了時とした。
１分毎のＨＦ、［ＨＦ／（ＬＦ＋ＨＦ）］％を各睡眠周期で平均し、睡眠周期毎の平均値を算出した。
【００３２】
【数１】
2015109888.tif 000004
【００３３】
次いで、第１睡眠周期と第２睡眠周期、第１睡眠周期と第３睡眠周期、又は第２睡眠周期と第３睡眠周期のＨＦの差と比、第２睡眠周期と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］％の差と比を算出した。
【００３４】
３．閾値の設定
得られたデータから、ＶＡＳ評価での「疲労なし」と「疲労あり」それぞれの対象者についてＲＯＣ曲線（Ｒｅｃｅｉｖｅｒ  Ｏｐｅｒａｔｉｎｇ  Ｃｈａｒａｃｔｅｒｉｓｔｉｃ  ｃｕｒｖｅ、受信者動作特性曲線）を作成し、表３に示す閾値を設定した。なお、感度と特異度の優れた独立変数のＲＯＣ曲線は、左上隅に近づいていくという事実から、本発明においてはこの左上隅との距離が最小となる点を閾値とした。
該閾値に基づく、疲労評価の感度、特異度、有効度を求め、表３に示した。なお、表３〜表１０では、第１睡眠周期はＣｙｃｌｅ１、第２睡眠周期はＣｙｃｌｅ２、第３睡眠周期はＣｙｃｌｅ３と表記した。
感度・特異度・有効度の算出方法は表２のとおりである。
【００３５】
【表２】
2015109888.tif 000005
【００３６】
【表３】
2015109888.tif 000006
【００３７】
実施例１
発明評価用データの収集で得たデータを利用し、各項目について表４に示す閾値を設定し対象者の疲労を評価した。その結果、感度が５０％以上で特異度が最大であった。
【００３８】
【表４】
2015109888.tif 000007
【００３９】
実施例２
発明評価用データの収集で得たデータを利用し、各項目について表５に示す閾値を設定し対象者の疲労を評価した。その結果、特異度が５０％以上で感度が最大であった。
【００４０】
【表５】
2015109888.tif 000008
【００４１】
実施例３
発明評価用データの収集で得たデータを利用し、各項目について表６に示す閾値を設定し対象者の疲労を評価した。その結果、有効度が高かった。
【００４２】
【表６】
2015109888.tif 000009
【００４３】
実施例４
発明評価用データの収集で得たデータを利用し、各項目について表７に示す閾値を設定し対象者の疲労を評価した。その結果、特異度が９０％以上と高く偽陽性が少ないため、疲労度が低い対象者の抽出に有効と考えられた。
【００４４】
【表７】
2015109888.tif 000010
【００４５】
実施例５
発明評価用データの収集で得たデータを利用し、各項目について表８に示す閾値を設定し対象者の疲労を評価した。その結果、感度が９０％以上と高く偽陰性が少ないため、疲労度が高い対象者の抽出に有効と考えられた。
【００４６】
【表８】
2015109888.tif 000011
【００４７】
実施例６
対象者の睡眠時間を９０分間で分割し、睡眠周期を９０分間隔で固定して第１〜第３の睡眠周期毎のＨＦの平均値、及び［ＨＦ／（ＬＦ＋ＨＦ）］％の平均値を求めた。
各項目について表９に示す閾値を設定し、疲労評価の感度、特異度、有効度を求めた。その結果、感度、特異度、有効度共に５０％前後と評価可能であることが確認された。
【００４８】
【表９】
2015109888.tif 000012
【００４９】
実施例７
本発明の頑健性と汎用性の確認のために、発明評価用データの収集とは異なる対象者の疲労判定を行った。
〔対象者〕
週５日（月曜日〜金曜日勤務、週休２日（土日））の日勤労働者（標準労働時間８：３０〜１７：００、１時間の昼休み含む）のうち１６名の健常男性労働者（２６〜５４歳）を対象として、連続１０日間の就寝時の主観的な疲労と、睡眠中の第１〜第３睡眠周期までの心拍変動を測定した。のべ７６日分のデータを解析した。
【００５０】
〔疲労評価〕
１．ＶＡＳ評価
対象者の主観的な疲労の評価は、実施例１と同様にＶＡＳ評価により行った。
疲労感ＶＡＳ（Ｖｉｓｕａｌ  Ａｎａｌｏｇｕｅ  Ｓｃａｌｅ）検査を３分位に分割し、下位３分の１（平均±標準偏差３１．０±１２．８ｍｍ、最小値〜最大値２〜４７ｍｍ）を「疲労なし」、上位３分の１（平均±標準偏差７６．３±８．５ｍｍ、最小値〜最大値６４〜９２ｍｍ）を「疲労あり」とした。なお、中位３分の１は平均±標準偏差５４．９±４．８ｍｍ、最小値〜最大値４８〜６３ｍｍである。
【００５１】
２．睡眠中の心拍変動の測定
発明評価用データの収集と同様に、睡眠時の心拍変動は、心臓の拍動とともに発生する指尖容積脈波から脈波間隔（Ｒ−Ｒ間隔）を測定し（東芝体動計ＮＥＭ−Ｔ１  （株）パラマテック）、１分間の脈波間隔を高速フーリエ変化で周波数解析を行った（東芝汎用睡眠解析ソフト「ねむリズム」ＮＥＭ−ＳＳ１、東芝ホームアプライアンス（株））。
また、実施例１と同様に睡眠周期を算出し、第１睡眠周期と第２睡眠周期、第１睡眠周期と第３睡眠周期、又は第２睡眠周期と第３睡眠周期のＨＦの差と比、第２睡眠周期と第３睡眠周期の［ＨＦ／（ＬＦ＋ＨＦ）］％の差と比を算出した。
【００５２】
３．閾値
閾値は発明評価用データの収集で得たデータを利用し各項目について表１０に示す閾値を設定した。該閾値に基づく、疲労評価の感度、特異度、有効度を求め、表１０に示した。
【００５３】
【表１０】
2015109888.tif 000013
【００５４】
その結果、良好な感度、特異度、有効度が示され、本発明の頑健性、汎用性が確認された。
【符号の説明】
【００５５】
１  算出手段１
２  算出手段２
３  算出手段３
４  閾値記憶手段
５  出力・記録装置
６  心拍間隔測定計
７  睡眠周期測定計
８  解析装置
【図１】
2015109888.tif 000016
【図２】
2015109888.tif 000017
【図３】
2015109888.tif 000018