WO2019230235A1 - ストレス評価装置、ストレス評価方法及びプログラム - Google Patents

ストレス評価装置、ストレス評価方法及びプログラム Download PDF

Info

Publication number
WO2019230235A1
WO2019230235A1 PCT/JP2019/016437 JP2019016437W WO2019230235A1 WO 2019230235 A1 WO2019230235 A1 WO 2019230235A1 JP 2019016437 W JP2019016437 W JP 2019016437W WO 2019230235 A1 WO2019230235 A1 WO 2019230235A1
Authority
WO
WIPO (PCT)
Prior art keywords
change
amount
heart rate
stress
threshold value
Prior art date
Application number
PCT/JP2019/016437
Other languages
English (en)
French (fr)
Inventor
武央 頭川
幸弘 森田
Original Assignee
パナソニックＩｐマネジメント株式会社
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
Priority claimed from JP2019067219A external-priority patent/JP7262006B2/ja
Application filed by パナソニックＩｐマネジメント株式会社 filed Critical パナソニックＩｐマネジメント株式会社
Priority to CN201980017146.8A priority Critical patent/CN111818850A/zh
Publication of WO2019230235A1 publication Critical patent/WO2019230235A1/ja
Priority to US17/029,128 priority patent/US20210000355A1/en

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle

Definitions

  • the present disclosure relates to a stress evaluation apparatus, a stress evaluation method, and a program for determining a factor of stress of a measurement subject.
  • Patent Document 1 calculates an activity intensity of a measurement subject based on a detection value of an acceleration sensor, and a biometric index such as heart rate, heartbeat waveform, blood pressure, blood oxygen saturation, body temperature, or sweating degree, and the like.
  • a biometric index such as heart rate, heartbeat waveform, blood pressure, blood oxygen saturation, body temperature, or sweating degree, and the like.
  • Patent Document 2 provides a stress relieving method and the like to the subject by analyzing and judging the subject's stress state together with surrounding conditions based on the subject's biological index and behavior information.
  • a life support device and a life support method are disclosed.
  • the present disclosure provides a stress evaluation apparatus, a stress evaluation method, and a program capable of determining a factor of stress of a measurement subject.
  • a stress evaluation device includes a first sensor unit that measures a heart rate and heart rate fluctuation of a measurement person, (i) a change amount of the heart rate, and (ii) a change amount of the heart rate fluctuation.
  • a calculation unit that calculates, and (ii) a determination unit that determines a factor of the stress of the person to be measured based on the change amount of the heart rate and (ii) the change amount of the heartbeat fluctuation, and outputs information based on the determination result
  • the amount of change in the heart rate is the amount of change to the heart rate measured by the first sensor unit from a reference that is the heart rate at rest of the subject
  • the determination unit includes (I) the heart rate A comparison of the magnitude relationship between the amount of change in the value and the first threshold, and , By comparing the magnitude relation between (II) the
  • the stress evaluation method includes an acquisition step of acquiring the measured heart rate and heart rate fluctuation of the measured person, (i) a change amount of the heart rate, and (ii) heart rate fluctuation.
  • the change amount of the heart rate is a change amount to the heart rate measured by the first sensor unit from a reference which is a heart rate at rest of the subject, and a change amount of the heart rate fluctuation Is a change amount of the heart rate fluctuation measured by the first sensor unit from a reference which is a heart rate fluctuation of the measurement subject at rest, and in the determination step, (I) a change amount of the heart rate And the first threshold It compares the magnitude relation, and, by comparing the magnitude relation between (II) the variation of the heartbeat fluctuation and the second threshold value, determining the cause of the stress.
  • the stress evaluation device According to the stress evaluation device, the stress evaluation method, and the program of the present disclosure, it is possible to evaluate the stress factor of the measurement subject.
  • FIG. 1 is a diagram in which the amount of change in the biometric index for each stress factor of each of the 20 subjects is plotted.
  • FIG. 2 is a diagram showing an average value of the amount of change in the biometric index for each stress factor plotted in FIG.
  • FIG. 3 is a schematic configuration diagram showing an example of the configuration of the stress evaluation apparatus according to the first embodiment.
  • FIG. 4 is a block diagram showing a specific example of a stress evaluation apparatus based on the configuration of FIG.
  • FIG. 5 is a flowchart for explaining the stress evaluation method according to the first embodiment.
  • FIG. 6 is a diagram illustrating an example of heart rate information obtained by the stress evaluation device according to the first embodiment.
  • FIG. 7 is a diagram for explaining a method of calculating the fluctuation amount of the heartbeat interval (RRI).
  • RRI heartbeat interval
  • FIG. 8 is a diagram for explaining an example of use of the stress evaluation apparatus according to the first embodiment.
  • FIG. 9A is a diagram in which the amount of change in the biometric index for each stress factor of each of the 20 subjects is plotted.
  • FIG. 9B is a diagram when FIG. 9A is viewed from the positive side of the axis indicating the amount of change in RRI.
  • FIG. 9C is a diagram when FIG. 9A is viewed from the negative side of the axis indicating the amount of change in CvRR.
  • FIG. 9D is a diagram when FIG. 9A is viewed from the negative side of the axis indicating the amount of change in SC.
  • FIG. 9A is a diagram in which the amount of change in the biometric index for each stress factor of each of the 20 subjects is plotted.
  • FIG. 9B is a diagram when FIG. 9A is viewed from the positive side of the axis indicating the amount of change in RRI.
  • FIG. 9C is a diagram when FIG
  • FIG. 10A is a diagram showing an average value of the amount of change in the biometric index for each stress factor plotted in FIG. 9A.
  • FIG. 10B is a diagram when FIG. 10A is viewed from the positive side of the axis indicating the amount of change in RRI.
  • FIG. 10C is a diagram when FIG. 10A is viewed from the negative side of the axis indicating the amount of change in CvRR.
  • FIG. 10D is a diagram when FIG. 10A is viewed from the negative side of the axis indicating the amount of change in SC.
  • FIG. 11 is a schematic configuration diagram illustrating an example of the configuration of the stress evaluation apparatus according to the embodiment.
  • FIG. 12 is a block diagram showing a specific example of a stress evaluation apparatus based on the configuration of FIG.
  • FIG. 13 is a flowchart for explaining the stress evaluation method according to the second embodiment.
  • FIG. 14 is a diagram for explaining an example of use of the stress evaluation apparatus according to the second embodiment.
  • the stress determination system described in Patent Document 1 calculates the activity intensity of the measurement subject based on information obtained from the acceleration sensor, and calculates the heart rate, heartbeat waveform, blood pressure, blood oxygen saturation, body temperature,
  • the stress state of the person to be measured is determined based on a biological index such as the degree of sweating and the activity intensity.
  • the stress state in the daily life of the measurement subject is determined by measuring the biometric index only when the activity intensity is a certain value or less.
  • the life support system described in Patent Document 2 acquires not only biological information such as electrocardiograms and pulse waves but also behavior information of the measurement subject, and analyzes and determines the situation around the measurement subject.
  • a stress relieving method and the like are provided to the person being measured.
  • the present inventors conducted the following monitor test in order to find the relationship between stress factors and multiple types of biological indices obtained from biological information such as heartbeat information.
  • test subjects were 20- to 30-year-old adults and men and women who were university students, and 20 subjects were selected who did not show abnormal values in the results of questionnaires regarding their health and mental status.
  • the biological signal at the time of the subject's rest is a biological signal measured for 5 minutes in the same posture as the task executing posture before performing each of the tasks [1] to [4].
  • a biometric index was calculated from this biosignal and used as a reference value for calculating the amount of change in the biometric index.
  • the amount of change of the biometric index is a biometric index calculated from the biometric signal of the subject measured during the task execution based on the biometric index of the subject at rest.
  • the measured biological signals are an electrocardiogram (ECG), a breathing interval, a fingertip temperature (Skin Temperature: SKT), and a skin conductance (Skin Conductance: SC) of the fingertip. These biosignals were measured simultaneously. A plurality of types of biological indices were obtained from each biological signal. Hereinafter, the result of examination using ECG will be described.
  • ECG electrocardiogram
  • SKT fingertip temperature
  • SC skin conductance
  • a heartbeat interval (RR interval: RRI), which is the interval between the R wave peaks of two consecutive heartbeats, was calculated (see (a) of FIG. 7).
  • RRI is one index of heart rate.
  • CvRR coefficient of Variation of RR interval: CvRR
  • CvRR is one of indices of heartbeat fluctuation.
  • CvRR was calculated from RRI by standardizing RRI standard deviation SD in an arbitrary time zone with an average value of RRI in an arbitrary time zone, as shown in the following formula (1).
  • each continuous RRI is converted into a biaxial relationship between time and RRI, and further converted into RRI equidistant time series data (see FIG. 7B), and then fast Fourier transform (Fast Fourier Transform).
  • Frequency analysis was performed using Transform (FFT) (see FIG. 7C).
  • FFT Fast Fourier transform
  • HF High Frequency
  • LF Low Frequency
  • HF and LF are indexes of heartbeat fluctuation.
  • HF is an integral value of the power spectrum in the high frequency region of 0.14 Hz to 0.4 Hz, and is considered to reflect the amount of activity of the parasympathetic nerve.
  • LF is an integral value of the power spectrum in the low frequency range of 0.04 Hz to 0.14 Hz, and is considered to reflect the amount of activity of the sympathetic nerve and the parasympathetic nerve. Note that data for performing frequency analysis using FFT was 60-second heart rate fluctuation data, and frequency conversion was performed at intervals of 5 seconds.
  • the biometric index when the measurement subject is at rest and the biometric index measured while the measurement subject is performing the task are average values of the biometric index for 240 seconds after 60 seconds from the start of measurement.
  • the amount of change of the biometric index is the amount of change from the reference, which is the average value of the biometric index when the subject is at rest, to the average value of the biometric index measured while the measured subject is performing the task. It is.
  • the amount of change is expressed as a ratio or a difference. When the change amount of the biometric index is expressed as a ratio, the change amount of the biometric index is calculated using the following equation (2).
  • the determination accuracy was 46.3%. Therefore, when the change amount of LF and the change amount of HF are used, the determination accuracy is greatly reduced as compared with the combination including the change amount of RRI. From the above examination, it has been found that when the amount of change in RRI and the amount of change in CvRR are used, the cause of stress can be determined with relatively high accuracy.
  • FIG. 1 is a diagram in which the amount of change in the biometric index for each stress factor of each of the 20 subjects is plotted. Since stress 1 and 2 related to thought fatigue both showed similar results, they were illustrated as stress related to thought fatigue. From FIG. 1, it was found that the change amount of the biometric index varies depending on the type of task to be executed. In order to make the change tendency clearer, an average value of the amount of change of the biometric index of 20 subjects was obtained.
  • FIG. 2 is a diagram showing an average value of the amount of change in the biometric index for each stress factor of 20 subjects. From FIG. 2, it was found that the change amount of the biometric index has the following characteristic change tendency depending on the factor of stress.
  • the amount of change in RRI greatly shifts to the negative side (that is, the heart rate increases), and the amount of change in CvRR tends to shift to the positive side.
  • the RRI change amount tends to shift to the plus side (that is, the heart rate becomes smaller), and the CvRR change amount tends to shift slightly to the minus side.
  • the stress factor is thought fatigue, the amount of change in RRI shifts slightly to the minus side (that is, the heart rate does not change much), and the amount of change in CvRR tends to greatly shift to the minus side. I found out.
  • the present inventors have a tendency that the amount of change of each biometric index changes due to stress factors, and in particular, both the amount of change of the biometric index related to heart rate and heart rate fluctuations.
  • the stress factor can be determined more accurately than when using either one as the determination index.
  • the invention of the apparatus for determining the stress factor and the intensity of the stress of the person to be measured by comparing the amount of change of a plurality of types of biological indexes obtained from the person to be measured and the threshold value. I thought.
  • the present disclosure provides a stress evaluation apparatus, a stress evaluation method, and a program that can determine a factor of stress of a measurement subject.
  • a stress evaluation device includes a first sensor unit that measures a heart rate and heart rate fluctuation of a measurement person, (i) a change amount of the heart rate, and (ii) a change amount of the heart rate fluctuation.
  • a calculation unit that calculates, and (ii) a determination unit that determines a factor of the stress of the person to be measured based on the change amount of the heart rate and (ii) the change amount of the heartbeat fluctuation, and outputs information based on the determination result
  • the amount of change in the heart rate is the amount of change to the heart rate measured by the first sensor unit from a reference that is the heart rate at rest of the subject
  • the determination unit includes (I) the heart rate A comparison of the magnitude relationship between the amount of change in the value and the first threshold, and , By comparing the magnitude relation between (II) the
  • the factor of stress can be determined by comparing the magnitude relationship between the amount of change of each biological index and the threshold value of each biological index.
  • the change amount of the heart rate is the change amount of the heart rate measured at a first time
  • the change amount of the heart rate fluctuation is the second change amount.
  • the amount of change in the heart rate fluctuation measured over time, and the first threshold is an arbitrary value different from the first time and the second time based on the heart rate of the measurement subject at rest
  • the heart rate measured at time may be the heart rate fluctuation measured at the arbitrary time on the basis of the heart rate fluctuation at rest of the measurement subject.
  • the arbitrary time refers to, for example, the time when the measurement subject is in a state before feeling stress.
  • the first threshold value and the second threshold value can be set accurately.
  • each biometric index measured at a predetermined time such as during sleep or immediately before going to bed may be set as the threshold value of each biometric index.
  • the factor of stress can be determined more accurately.
  • the heartbeat fluctuation may be obtained by performing frequency analysis on the heartbeat interval of the measurement subject.
  • the stress evaluation apparatus can obtain information on the breathing interval and blood pressure from the frequency component of heartbeat fluctuation. Therefore, the stress evaluation apparatus can use a biometric index including detailed information of the measurement subject as an index (determination index) for determining the stress, and thus more accurately determine the factor of the measurement subject's stress. be able to.
  • the determination unit includes a change amount of the heart rate that is greater than the first threshold value, and a change amount of the heart rate fluctuation is greater than the second threshold value. May be determined to be a factor related to the person.
  • the stress factor is an interpersonal factor by comparing the magnitude relationship between the amount of change of each biometric index and the threshold value of each biometric index.
  • the determination unit includes a change amount of the heart rate that is greater than the first threshold value, and a change amount of the heart rate fluctuation is greater than the second threshold value. May be determined to be pain.
  • the stress factor is pain by comparing the magnitude relationship between the amount of change of each biological index and the threshold value of each biological index.
  • the determination unit includes a change amount of the heart rate that is smaller than the first threshold value, and a change amount of the heart rate fluctuation is less than the second threshold value. May be determined to be fatigue due to thought.
  • the determination unit further includes a difference between the change amount of the heart rate and the first threshold value, and a change amount of the heart rate fluctuation and the second amount.
  • the intensity of the stress may be determined according to a difference from a threshold value, and a determination result may be output as the information based on the determination result.
  • the person being measured can know the strength of his / her stress. Thereby, it becomes easy to have consciousness about control of stress, and it becomes easy to grasp
  • the person to be measured can recognize that the strength of the stress that can be withstood among a plurality of types of stress is different. As a result, the person to be measured can determine whether or not stress control is required immediately in accordance with the state of stress. Therefore, the measurement subject can efficiently control the stress, and thus can continuously control the stress.
  • the stress evaluation device further includes a presentation unit that presents the information based on the determination result output by the determination unit, wherein the information includes the stress factor, the stress It may include at least one selected from the group consisting of strength and stress reduction measures.
  • the presenting unit may present by voice.
  • the presenting unit may present an image.
  • the person being measured can visually know the state of stress and the control method of his / her own stress, so that he / she can be clearly aware of the control of his / her stress. Therefore, the measurement subject can continue to control his own stress.
  • the stress evaluation method includes an acquisition step of acquiring the measured heart rate and heart rate fluctuation of the measured person, (i) a change amount of the heart rate, and (ii) heart rate fluctuation.
  • the change amount of the heart rate is a change amount to the heart rate measured by the first sensor unit from a reference which is a heart rate at rest of the subject, and a change amount of the heart rate fluctuation Is a change amount of the heart rate fluctuation measured by the first sensor unit from a reference which is a heart rate fluctuation of the measurement subject at rest, and in the determination step, (I) a change amount of the heart rate And the first threshold It compares the magnitude relation, and, by comparing the magnitude relation between (II) the variation of the heartbeat fluctuation and the second threshold value, determining the cause of the stress.
  • the factor of stress can be determined by comparing the magnitude relationship between the amount of change of each biological index and the threshold value of each biological index.
  • FIG. 3 is a schematic configuration diagram of the stress evaluation apparatus 100 according to the present embodiment.
  • the stress evaluation device 100 includes a first sensor unit 11 a, a calculation unit 12, a determination unit 13, a presentation unit 14, and a storage unit 15.
  • the first sensor unit 11a includes a wearable first biological sensor 111a (see FIG. 4) that measures a biological signal of the measurement subject.
  • the first sensor unit 11a calculates a plurality of types of biological indices from the biological signal measured by the first biological sensor 111a, and outputs the calculated biological indices to the computing unit 12.
  • the calculation unit 12 calculates an average value (hereinafter also referred to as a reference value) of each biological index when the measurement subject is at rest and a threshold value of each biological index, and stores the calculated values in the storage unit 15. In addition, the calculation unit 12 calculates the average value of each measured biometric index and the amount of change of each biometric index, and outputs it to the determination unit 13.
  • the determination unit 13 determines a factor of stress of the measurement subject based on the amount of change of each biological index. More specifically, the determination unit 13 determines the factor of stress by comparing the magnitude relationship between the change amount of each biological index and the threshold value of each biological index. Further, the determination unit 13 determines the strength of the stress according to the difference between the change amount of each biometric index and the threshold value of each biometric index.
  • the determination unit 13 outputs information based on these determination results to the presentation unit 14.
  • the determination unit 13 causes the storage unit 15 to store information based on the determination result.
  • the presentation unit 14 presents information based on the determination result.
  • the stress evaluation apparatus 100 may include an input unit 16 (see FIG. 4) for inputting an instruction of the measurement subject (user).
  • the determination unit 13 causes the presentation unit 14 to present determination result information based on the measurement subject's instruction input to the input unit 16.
  • FIG. 4 is a block diagram showing a specific example of a stress evaluation apparatus based on the configuration of FIG.
  • the stress evaluation apparatus 100 includes a first sensor unit 11a including a first biological sensor 111a and a first signal processing unit 112a, a calculation unit 12, a determination unit 13, a presentation unit 14, A storage unit 15 and an input unit 16 are provided.
  • the first biological sensor 111a measures a biological signal of the measurement subject.
  • a biological signal is a signal of biological information.
  • the biological information is physiological information that is affected by stress such as heart rate, pulse, respiratory rate, blood oxygen saturation, blood pressure, or body temperature. From the ease of measurement, the biological information is, for example, heart rate information.
  • Heartbeat information is information obtained from a heartbeat.
  • the biological information may be pulse wave information.
  • the first biological sensor 111a is a sensor that acquires heartbeat information or pulse wave information.
  • the heart rate sensor is, for example, a sensor that includes a pair of detection electrodes that contact the surface of the body of the measurement subject.
  • the heartbeat information obtained by the heartbeat sensor is an electrical signal obtained by the heartbeat, for example, an electrocardiogram.
  • the heart rate sensor may be a conductive adhesive gel electrode or a dry electrode composed of conductive fibers.
  • the wearing part of the heart rate sensor is the chest, and the shape of the heart rate sensor is, for example, a wear type in which wear and electrodes are integrated.
  • the pulse wave sensor When the first biological sensor 111a is a sensor that acquires pulse wave information (hereinafter referred to as a pulse wave sensor), the pulse wave sensor reflects, for example, a change in blood volume in a blood vessel using a phototransistor and a photodiode. It is a sensor which measures by.
  • the pulse wave sensor is attached to the user's wrist and measures pulse wave information with the attached shape.
  • the attachment site of the pulse wave sensor may be an ankle, a finger, an upper arm, or the like.
  • the shape of the pulse wave sensor is not limited to a band type (for example, a wrist watch type), and may be a pasting type attached to a neck or the like, an eyeglass type, or the like.
  • the pulse wave sensor may be an image sensor that calculates pulse by measuring pulse wave information from a change in chromaticity of skin such as a face or a hand.
  • the biological signal measured by the first biological sensor 111a is output to the first signal processing unit 112a.
  • the first signal processing unit 112a calculates a plurality of types of biological indices from one biological signal measured by the first biological sensor 111a.
  • two types of biometric indicators biometric index 1 and biometric index 2 are calculated.
  • the plurality of types of biological indices are RRI, CvRR, HF, LF, and the like.
  • RRI is a heart rate index
  • CvRR, HF, and LF are heart rate fluctuation indexes.
  • the first signal processing unit 112a may calculate a biological index of fluctuations in respiratory rate and blood pressure from the frequency component of heartbeat fluctuation.
  • RRI and CvRR are combinations having relatively high determination accuracy among these plural types of biological indices.
  • the first signal processing unit 112 a outputs the calculated biometric index 1 and biometric index 2 to the calculation unit 12.
  • the calculation unit 12 acquires the biometric index 1 and the biometric index 2 output from the first signal processing unit 112a, and calculates the change amount of the biometric index 1 and the change amount of the biometric index 2 from the acquired biometric index 1 and biometric index 2. To do.
  • the amount of change of the biometric index is a measured biometric index based on a biometric index (hereinafter sometimes referred to as a reference value) measured when the measurement subject is at rest, and is represented by a difference or a ratio.
  • the reference value of each biometric index is stored in the storage unit 15.
  • the calculation unit 12 reads the reference values of the biometric index 1 and the biometric index 2 stored in the storage unit 15 and calculates the amount of change of the biometric index 1 and the biometric index 2 with respect to the reference value.
  • the calculation unit 12 outputs the calculated change amount of each biological index to the determination unit 13. Note that the reference value may vary depending on the season or the physiological period of the measurement subject, and may be updated every predetermined period.
  • the calculation unit 12 calculates a threshold value for each biological index.
  • the biological index 1 is, for example, a heart rate
  • the amount of change in heart rate is the amount of change in heart rate measured at the first time.
  • the first threshold value is a threshold value of the biometric index 1, for example, an RRI threshold value that is an index of heart rate.
  • the first threshold is a heart rate measured at an arbitrary time on the basis of the heart rate at rest of the measurement subject.
  • the biological index 2 is, for example, heart rate fluctuation
  • the amount of change in heart rate fluctuation is the amount of change in heart rate fluctuation measured at the second time.
  • the second threshold value is a threshold value of the biological index 2, for example, a CvRR threshold value that is an index of heartbeat fluctuation.
  • the second threshold is a heart rate fluctuation measured at an arbitrary time on the basis of the heart rate of the subject at rest. That is, these threshold values are changes in the biological index that is the difference or ratio between the measured value of the biological index measured at an arbitrary time different from the first time and the second time and the reference value.
  • the arbitrary time indicates, for example, the time when the measurement subject is in a state before feeling stress.
  • the first signal processing unit 112a may calculate a plurality of types of heart rates and heart rate fluctuations in a time division manner from one biological signal measured by the first biological sensor 111a.
  • the calculation unit 12 calculates the amount of change in heart rate measured at the first time, and calculates the amount of change in heart rate fluctuation measured at a second time different from the first time.
  • the calculation unit 12 reads the threshold value of each biometric index stored in the storage unit 15 and compares the magnitude relationship between the change amount of each biometric index and the threshold value. And the calculating part 12 determines the period when at least 1 of the variation
  • the stress generation period is a period during which the measurement subject feels stress.
  • the calculation unit 12 calculates a representative value of the change amount of each biometric index from the change amount of each biometric index during the stress occurrence period.
  • the average value of the change amount of each biometric index during the stress occurrence period may be used, and the value having the largest difference from the reference value (maximum value) May be used.
  • the determination unit 13 acquires the representative value of the change amount of the biometric index 1 and the biometric index 2 output from the calculation unit 12 and reads the first threshold value and the second threshold value stored in the storage unit 15.
  • the determination unit 13 compares the magnitude relationship between the representative value of the change amount of the biometric index 1 and the first threshold during the stress occurrence period, and compares the magnitude of the representative value of the change amount of the biometric index 2 with the second threshold value. By comparing the relationship, a factor of the stress of the measurement subject is determined. That is, the determination unit 13 determines a factor of stress for each stress occurrence period. Since the representative value of the change amount of the biometric index is an example of the change amount of the biometric index, hereinafter, the representative value of the change amount of the biometric index is also simply referred to as a change amount of the biometric index.
  • the determination unit 13 has a change amount of the biometric index 1 (here, heart rate) larger than the first threshold and a change amount of the biometric index 2 (here, heart rate fluctuation) is the second. If the threshold value is larger than the threshold value, it is determined that the stress factor is a factor related to the person. The determination unit 13 determines that the stress factor is pain when the change amount of the biometric index 1 is larger than the first threshold value and the change amount of the biometric index 2 is smaller than the second threshold value. To do. In addition, when the change amount of the biometric index 1 is smaller than the first threshold value and the change amount of the biometric index 2 is larger than the second threshold value, the determination unit 13 causes fatigue due to thought. Is determined.
  • the biometric index 1 here, heart rate
  • a change amount of the biometric index 2 here, heart rate fluctuation
  • the determination unit 13 determines the strength of the stress according to the difference between the change amount of the biometric index 1 and the first threshold value and the difference between the change amount of the biometric index 2 and the second threshold value.
  • the result is output as information based on the determination result.
  • the information based on the determination result includes, for example, at least one of a stress factor, a stress intensity, and a stress reduction measure.
  • the stress reduction measure is, for example, a stress elimination method or a stress avoidance method.
  • the stress reduction measure is included in a presentation information table described later.
  • the determination unit 13 reads out an appropriate stress reduction measure from the presentation information table stored in the storage unit 15 and outputs it to the presentation unit 14.
  • the determination unit 13 stores information based on the determination result in the storage unit 15. At this time, the determination unit 13 may associate the information on the time when the measurement subject feels stress with the information based on the determination result and store the information in the storage unit 15.
  • the presentation unit 14 presents information based on the determination result output by the determination unit 13.
  • the presentation unit 14 may present information based on the determination result by voice or an image.
  • the presentation unit 14 is, for example, a speaker.
  • the presentation part 14 presents the said information with an image, the presentation part 14 is a display, for example.
  • the storage unit 15 stores a reference value for each biometric index, a threshold value for each biometric index, a presentation information table, and the like.
  • the presentation information table is a table of presentation information such as stress reduction measures presented according to the stress factor and the intensity of the stress. As described above, the reference value and threshold value of each biometric index may be updated in a predetermined period. Similarly, the presentation information table may be updated in a predetermined period.
  • the storage unit 15 stores information based on the determination results such as the stress factor, the stress intensity, and the stress reduction measure output by the determination unit 13. At this time, the storage unit 15 may store information based on the determination result and the stress occurrence period in association with each other. Thereby, the measurement subject can call information based on the determination result at a desired timing. At this time, the determination unit 13 causes the presentation unit 14 to present information based on the determination result based on the measurement subject's operation input by the input unit 16.
  • the input unit 16 outputs an operation signal indicating an operation by the measurement subject to the determination unit 13.
  • the input unit 16 is, for example, a keyboard, a mouse, a touch panel, or a microphone.
  • the operation signal is a signal for setting an information extraction method based on the determination result or a presentation method in the presentation unit 14.
  • various types of determination results are presented to the presentation unit 14. For example, a change in stress during a predetermined period, a factor of stress that the measurement subject is likely to be influenced by, a stress reduction measure suitable for the measurement subject, and the like.
  • the measurement subject can grasp not only the short-term stress tendency but also the medium-term and long-term stress tendency. In this way, the person to be measured can know an effective stress reduction measure suitable for himself / herself, and therefore can control mid- to long-term stress.
  • FIG. 5 is a flowchart for explaining the stress evaluation method according to the embodiment.
  • the stress evaluation method includes an acquisition step S10 for acquiring the measured heart rate and heart rate fluctuation of the measured person, (i) a change amount of the heart rate, and (ii) a change amount of the heart rate fluctuation.
  • Calculation step S20, and determination step S30 for determining a factor of stress of the measurement subject based on the amount of change in heart rate and the amount of change in heart rate fluctuation, and outputting information based on the determination result.
  • the amount of change in heart rate is the amount of change to the heart rate measured by the first sensor unit 11a from the reference, which is the heart rate at rest of the subject, and the amount of change in heart rate fluctuation is the amount of change in the subject.
  • the present embodiment further includes a presentation step S40 for presenting information based on the determination result of determination step S30.
  • the calculation unit 12 acquires a plurality of types of biological indices (here, heart rate and heart rate fluctuation) of the measurement subject measured by the first sensor unit 11a.
  • heart rate information here, an electrocardiogram
  • biological indexes such as a heart rate index and a heart rate fluctuation index are calculated by the first signal processing unit 112a.
  • the biological information is not limited to heartbeat information, and may be physiological information that is affected by stress such as pulse wave information.
  • the heart rate information is simple in a state where the burden on the subject is less than other biological information such as pulse, respiratory rate, blood pressure, and blood oxygen saturation, and It can be measured in real time. Therefore, by using the heartbeat information of the person to be measured as the biological information, it is possible to appropriately evaluate the stress state of the person to be measured.
  • the biological indices obtained from the heartbeat information are RRI that is an index of heart rate, CvRR, LF, HF, and LF / HF that are indices of heartbeat fluctuation.
  • RRI an index of heart rate
  • CvRR cardiac rate
  • LF LF
  • HF HF
  • LF / HF indices of heartbeat fluctuation.
  • the combination of these biomarkers can determine the cause of stress with relatively high determination accuracy, so that highly reliable evaluation can be obtained.
  • FIG. 6 is a diagram illustrating an example of heart rate information obtained by the first sensor unit 11a of the stress evaluation apparatus 100 according to the present embodiment.
  • the heart rate information is, for example, an electrocardiogram and has an electrocardiogram waveform shown in FIG.
  • the ECG waveform reflects the P wave that reflects the electrical excitation of the atrium, the Q wave, the R wave, and the S wave that reflect the electrical excitation of the ventricle, and the process of repolarization of the excited ventricular cardiomyocytes. It consists of T waves.
  • the wave height (potential difference) of the R wave is the largest and is most robust against noise such as myoelectric potential. Therefore, the R wave peak interval between two consecutive heartbeats in these electrocardiographic waveforms, that is, the heartbeat interval (RRI) is calculated.
  • the heart rate is calculated by multiplying the reciprocal of RRI by 60.
  • the CvRR is calculated by normalizing the standard deviation SD of the RRI in an arbitrary time zone from the RRI using the above formula (2) with the average value of the heartbeat interval.
  • the first signal processing unit 112a detects an electrical signal (R wave) generated when the left ventricle contracts suddenly and pumps blood from the heart from the heartbeat information obtained by the first biological sensor 111a, and calculates the RRI. calculate.
  • R wave an electrical signal generated when the left ventricle contracts suddenly and pumps blood from the heart from the heartbeat information obtained by the first biological sensor 111a
  • calculates the RRI. calculate for detecting the R wave, for example, a known method such as the Pan & Topkins method is used.
  • FIG. 7 is a diagram for explaining a technique for calculating the fluctuation amount of the heartbeat interval (RRI).
  • the first signal processing unit 112a calculates the RRI fluctuation amount from the obtained R-wave detection data as follows.
  • the first signal processing unit 112a calculates the RRI that is the interval between the peaks of the R waves of two consecutive heartbeats.
  • the first signal processing unit 112a converts each calculated RRI into a biaxial relationship between time and RRI. Since the converted data is discrete data with unequal intervals, the calculation unit 12 converts the converted RRI time-series data into equal-interval time-series data shown in FIG. 7B.
  • the computing unit 12 obtains a frequency component of heartbeat fluctuation shown in FIG. 7C by performing frequency analysis on the equally-spaced time series data using, for example, fast Fourier transform (FFT).
  • FFT fast Fourier transform
  • the frequency component of heart rate variability can be divided into, for example, a high frequency component HF and a low frequency component LF.
  • HF is considered to reflect the amount of parasympathetic nerve activity.
  • LF is also considered to reflect the amount of sympathetic and parasympathetic activity. Therefore, LF / HF, which is the ratio of LF to HF, is considered to indicate the amount of sympathetic nerve activity.
  • a plurality of types of biological indices are calculated from the heartbeat information.
  • the calculation unit 12 acquires two types of biological indexes (here, heart rate and heart rate fluctuation) from these biological indexes.
  • the calculation unit 12 calculates the amount of change in the two types of biological indices acquired in the acquisition step S10.
  • the amount of change of each biological index is the ratio or difference between the reference value of each biological index and the value of each acquired biological index, with the value of each biological index when the subject is at rest as a reference value. Obtained by calculation.
  • the calculation unit 12 reads out and uses the reference value of each biological index stored in the storage unit 15.
  • the amount of change of each biological index is calculated by, for example, subtracting the reference value of each biological index from the value of each biological index acquired in the acquisition step S10 when the amount of change is represented by a difference.
  • the change amount of the heart rate is calculated by subtracting the reference value of the heart rate from the value of the heart rate of the measurement subject acquired in the acquisition step S10.
  • the change amount is expressed as a ratio, it is calculated by dividing the value of each biological index acquired in the acquisition step S10 by the reference value of each biological index.
  • the amount of change in the heart rate is calculated by dividing the value of the heart rate of the measured person acquired in the acquisition step S10 by the reference value of the heart rate.
  • the calculation unit 12 calculates the change amount of each biological index.
  • the determination unit 13 determines a stress factor based on the change amount of each biological index calculated in the calculation step S20.
  • the determination unit 13 determines the factor of the stress of the measurement subject by comparing the magnitude relationship between the amount of change of each biological index and the threshold value of each biological index. Specifically, in determination step S ⁇ b> 30, when the amount of change in heart rate is greater than the first threshold and the amount of change in heart rate fluctuation is greater than the second threshold, the determination unit 13 determines that the stress factor is It is determined that this is a factor related to the person.
  • the determination unit 13 determines that the stress factor is pain when the change amount of the biometric index 1 is larger than the first threshold value and the change amount of the biometric index 2 is smaller than the second threshold value. To do.
  • the determination unit 13 causes fatigue due to thought. Is determined.
  • the determination unit 13 determines the strength of the stress according to the difference between the change amount of the biometric index 1 and the first threshold value and the difference between the change amount of the biometric index 2 and the second threshold value. The result is output as information based on the determination result.
  • the first threshold value is a heart rate threshold value, and the heart rate measured at an arbitrary time different from the first time and the second time on the basis of the heart rate at rest of the measurement subject. It is.
  • the second threshold value is a threshold value of heart rate fluctuation, and is a heart rate fluctuation measured at an arbitrary time different from the first time and the second time with reference to the heart rate fluctuation at rest of the measurement subject. .
  • These threshold values are calculated by the calculation unit 12 and stored in the storage unit 15.
  • the determination unit 13 reads out and uses the threshold value of each biological index stored in the storage unit 15.
  • the arbitrary time indicates, for example, the time when the measurement subject is in a state before feeling stress.
  • a threshold value when the amount of change of each biological index is a positive value and a threshold value when the amount of change of each biological index is a negative value are set.
  • the reference value is the zero point of change.
  • the magnitude relationship between the change amount of each biological index and the threshold value is compared as follows. When the change amount of the biometric index is a positive value, the magnitude relationship between the change amount of the biometric index and the positive threshold value is compared. When the change amount of the biometric index is a negative value, the magnitude relationship between the absolute value of the change amount of the biometric index and the absolute value of the negative threshold is compared.
  • the threshold value of each biometric index may be a fixed value, may be updated in a predetermined period, or may be updated each time based on daily measurement.
  • the threshold value may be calculated by relatively simple machine learning such as linear discrimination or decision tree. Therefore, since the determination reference value and threshold value suitable for the person to be measured can be set, the cause of stress can be determined with higher accuracy.
  • the factor of stress of the measurement subject is determined by comparing the magnitude relationship between the change amount of each biological index and the threshold value of each biological index.
  • the presentation unit 14 presents information based on the determination result determined by the determination unit 13.
  • the presentation unit 14 may present information based on the determination result by voice, or may present it by an image.
  • the information based on the determination result includes at least one of a stress factor, a stress intensity, and a stress reduction measure.
  • the presentation unit 14 displays various types of determination results based on the settings input by the measurement subject using the input unit 16.
  • FIG. 8 is a diagram for explaining an example of use of the stress evaluation apparatus 100 according to the present embodiment.
  • the stress evaluation apparatus 100 includes a first biological sensor 111a that is a part of the first sensor unit 11a and an evaluation terminal 20 that includes components other than the first biological sensor 111a.
  • the measurement subject wears the first living body sensor 111a so as to contact the skin of the chest, and measures an electrocardiogram (ECG).
  • ECG electrocardiogram
  • the first biological sensor 111a may be a conductive adhesive gel electrode or a dry electrode composed of conductive fibers or the like.
  • the first biological sensor 111a transmits the measured heartbeat electrical signal to the evaluation terminal 20 by communication.
  • the communication method may be wireless communication such as Bluetooth (registered trademark) or wired communication.
  • the evaluation terminal 20 includes a first signal processing unit 112a, a calculation unit 12, a determination unit 13, a presentation unit 14, a storage unit 15, and an input unit 16 of the first sensor unit 11a.
  • the 1st signal processing part 112a receives the electrical signal of the heartbeat transmitted by communication from the 1st living body sensor 111a.
  • the first signal processing unit 112 a calculates RRI that is a heart rate index and CvRR that is a heart rate fluctuation index from the received heartbeat electrical signal, and outputs these biological indices to the calculation unit 12.
  • the calculation unit 12 acquires the RRI and CvRR output from the first signal processing unit 112a, and reads the RRI reference value and the CvRR reference value stored in the storage unit 15.
  • the calculation unit 12 calculates the amount of change of the biometric index, which is a biometric index, based on the read reference value.
  • the amount of change in the biometric index is represented by a difference or a ratio. In the present embodiment, the amount of change is expressed as a ratio.
  • the calculation unit 12 calculates the threshold value of each biological index and outputs it to the storage unit 15.
  • a threshold value when the change amount of each biometric index becomes a positive value and a threshold value when the change amount of each biometric index becomes a negative value are set.
  • the reference value is zero change.
  • the positive threshold value is a value larger than the reference value, and the first threshold value 1a in the change amount graph 120 (hereinafter, positive value).
  • the threshold value 1a) and the second threshold value 2a hereinafter, positive threshold value 2a).
  • the negative threshold value is smaller than the reference value, and the first threshold value 1b (hereinafter, negative threshold value 1b) and the first threshold value in the change amount graph 120 are the same. 2 threshold value 2b (hereinafter, negative threshold value 2b).
  • the calculation unit 12 calculates a reference value for each biometric index and outputs it to the storage unit 15.
  • the reference value of each biometric index the amount of change of each biometric index is zero.
  • the reference value is a solid line between the positive threshold value 1a and the negative threshold value 1b.
  • the positive threshold value and the negative threshold value may be set at equal intervals with the reference value (change amount zero) interposed therebetween, or may not be set at equal intervals with the reference value interposed therebetween. These threshold values may be appropriately set according to the magnitude of the change amount of each biological index.
  • the determination unit 13 acquires the amount of change of each biometric index output by the calculation unit 12 and reads the threshold value of each biometric index stored in the storage unit 15.
  • the determination unit 13 compares the magnitude of the change amount of each biological index and the threshold value of each biological index, and determines a factor of stress. For example, when the change amount of each biological index is a positive value, the determination unit 13 compares the magnitude relationship between the change amount of each biological index and the positive threshold value. When the change amount of each biological index is a negative value, the determination unit 13 compares the magnitude relationship between the absolute value of the change amount of each biological index and the absolute value of the negative threshold value.
  • the change amount graph 120 and the determination table 130 will be described in more detail.
  • the determination unit 13 determines that the stress factor felt by the measured person in the period A1 is a factor related to the person.
  • the change amount of RRI is larger than the positive threshold value 1a
  • the absolute value of change amount of CvRR is smaller than the absolute value of the negative threshold value 2b. Therefore, the determination unit 13 determines that the stress factor felt by the measurement subject during the period B1 is pain.
  • the determination unit 13 determines that the stress factor felt by the measured person in the period C1 is fatigue due to thinking (thinking fatigue).
  • the change amount of each biological index based on the reference value is indicated by the direction and number of arrows.
  • a horizontal arrow indicates that the change amount of the biometric index is not accompanied by a change exceeding the threshold value.
  • the determination unit 13 determines the difference between the absolute value of the RRI change amount and the absolute value of the first threshold value, and the difference between the absolute value of the CvRR change amount and the absolute value of the second threshold value, Determine the strength of the stress.
  • the determination unit 13 outputs information based on these determination results to the presentation unit 14.
  • the presentation unit 14 is, for example, a smartphone display.
  • the determination unit 13 can call the information based on the determination result at a desired timing.
  • the determination unit 13 causes the presentation unit 14 to present information based on the determination result based on the operation of the measurement subject input by the input unit 16 such as a touch panel.
  • the determination unit 13 presents the presentation information 140 to the presentation unit 14 based on the instruction of the measurement subject.
  • the presentation information 140 includes a time when the measurement subject feels stress, a factor of stress, and a measure for reducing stress.
  • the stress reduction measure is, for example, a message for proposing a stress relieving method or a stress avoiding method according to a stress factor.
  • the message is, for example, if the stress factor is thought fatigue, let's take a break or stretch, and if it's a factor related to the person, meditate a little or take a deep breath. And so on.
  • the measurement subject can grasp his / her stress state and appropriate stress reduction measures more accurately than in the past. As a result, the measurement subject can appropriately and efficiently control their own stress, and therefore can continuously control the stress.
  • the present inventors conducted the following monitor test in order to find the relationship between the stress factor and the biological index obtained from the biological information such as heart rate information and sweating information.
  • test subjects were 20- to 30-year-old adults and men and women who were university students, and 20 subjects were selected who did not show abnormal values in the results of questionnaires regarding their health and mental status.
  • the biological signal at the time of the subject's rest is a biological signal measured for 5 minutes in the same posture as the task executing posture before performing each of the tasks [1] to [4].
  • a biometric index was calculated from this biosignal and used as a reference value for calculating the amount of change in the biometric index.
  • the amount of change of the biometric index is a biometric index calculated from the biometric signal of the subject measured during the task execution based on the biometric index of the subject at rest.
  • the measured biological signals are an electrocardiogram (ECG), a breathing interval, a fingertip temperature (Skin Temperature: SKT), and a skin conductance (Skin Conductance: SC) of the fingertip. These biosignals were measured simultaneously. A plurality of types of biological indices were obtained from each biological signal.
  • ECG electrocardiogram
  • SKT fingertip temperature
  • SC skin conductance
  • the calculation method of the biometric index varies depending on each biometric index.
  • the biometric index is SKT
  • the SKT is obtained by averaging the fingertip temperature over an arbitrary interval. Since CvRR, HF, and LF are also as described above, description thereof is omitted here.
  • the determination accuracy was 81.3%, and in the determination by the decision tree, the determination accuracy was 66.3%. Therefore, it was found that even if the number of changes in the biometric index used for determining the stress factor is reduced to three, relatively high determination accuracy is maintained.
  • the determination accuracy was 62.5%. Therefore, it was found that when the change amount of RRI that is an index of heart rate is excluded from the change amount of the biometric index used for the determination of the stress factor, the determination accuracy is remarkably lowered.
  • FIG. 9A is a diagram in which the amount of change in the biometric index for each factor of stress for each of the 20 subjects is plotted.
  • FIG. 9B is a diagram when FIG. 9A is viewed from the positive side of the axis indicating the amount of change in RRI.
  • FIG. 9C is a diagram when FIG. 9A is viewed from the negative side of the axis indicating the amount of change in CvRR.
  • FIG. 9D is a diagram when FIG. 9A is viewed from the negative side of the axis indicating the amount of change in SC.
  • FIG. 10A is a diagram illustrating an average value of the amount of change in the biometric index for each stress factor of the 20 subjects plotted in FIG. 9A.
  • FIG. 10B is a diagram when FIG. 10A is viewed from the positive side of the axis indicating the amount of change in RRI.
  • FIG. 10C is a diagram when FIG. 10A is viewed from the negative side of the axis indicating the amount of change in CvRR.
  • FIG. 10D is a diagram when FIG. 10A is viewed from the negative side of the axis indicating the amount of change in SC. From FIG. 10A to FIG. 10D, it was found that the change amount of the biometric index has the following characteristic change tendency depending on the factor of stress.
  • the amount of change in RRI greatly shifts to the minus side (that is, the heart rate increases), the amount of change in CvRR shifts to the plus side, and the amount of change in SC increases to the plus side.
  • the stress factor is pain, the RRI change amount shifts to the plus side (that is, the heart rate decreases), the CvRR change amount slightly shifts to the minus side, and the SC change amount increases to the plus side. There is a tendency to shift greatly.
  • the present inventors have a tendency of a predetermined change in the amount of change of each biometric index due to stress factors, and in particular, (i) heart rate, (ii) heart rate fluctuation, and (Iii) The knowledge that the cause of stress can be discriminated with relatively high accuracy by using the amount of change of the biometric index relating to skin conductance or skin temperature as a determination index. And based on this examination result, it resulted in invention of the apparatus which determines the factor of a to-be-measured person's stress by comparing the variation
  • the present disclosure provides a stress evaluation apparatus, a stress evaluation method, and a program that can determine a factor of stress of a measurement subject.
  • the stress evaluation device further includes a second sensor unit that measures at least one of the skin conductance or the skin temperature of the subject, and the calculation unit further includes (iii) a skin conductance.
  • the amount of change or the amount of change in skin temperature is calculated, and the amount of change in skin conductance is calculated based on the skin conductance measured by the second sensor unit from a reference that is the skin conductance at rest of the subject.
  • the amount of change in the skin temperature is the amount of change to the skin temperature measured by the second sensor unit from a reference that is the skin temperature at rest of the subject, and the determination
  • the unit includes (III) the magnitude of the change in the skin conductance or the change in the skin temperature and the third threshold.
  • the factor of stress can be determined by comparing the magnitude relationship between the amount of change of each biological index and the threshold value of each biological index.
  • the change amount of the heart rate is the change amount of the heart rate measured at a first time
  • the change amount of the heart rate fluctuation is the second change amount.
  • the amount of change in the heart rate fluctuation measured in time is the amount of change in skin conductance or the skin temperature measured in a third time
  • the first threshold value is the heart rate measured at an arbitrary time different from the first, second, and third times, based on the resting heart rate of the subject.
  • the second threshold is the heartbeat fluctuation measured at the arbitrary time with reference to the heartbeat fluctuation at rest of the subject
  • the third threshold is the rest of the subject.
  • the aforementioned skin conductance as a standard It said skin conductance was measured time, or the referenced skin temperature at rest of the subject, may be the skin temperature measured on the arbitrary time.
  • the arbitrary time refers to, for example, the time when the measurement subject is in a state before feeling stress.
  • a 1st threshold value, a 2nd threshold value, and a 3rd threshold value can be set correctly.
  • each biometric index measured at a predetermined time such as during sleep or immediately before going to bed may be set as the threshold value of each biometric index.
  • the threshold value which considered the menstrual change of a woman, a secular variation, etc. can be set, without measuring a subject every time, the factor of stress can be determined more correctly.
  • the heartbeat fluctuation may be obtained by performing frequency analysis on the heartbeat interval of the measurement subject.
  • the stress evaluation apparatus can obtain information on the breathing interval and blood pressure from the frequency component of heartbeat fluctuation. Thereby, since the stress evaluation apparatus can use the biometric index including the detailed information of the measurement subject as the determination index, the stress factor of the measurement target can be determined more accurately.
  • the stress evaluation apparatus can obtain information on the breathing interval and blood pressure from the frequency component of heartbeat fluctuation. Therefore, the stress evaluation apparatus can use a biometric index including the detailed state of the person being measured as an index (determination index) for determining the stress, and thus more accurately determine the cause of the stress of the person being measured. be able to.
  • the determination unit includes a change amount of the heart rate that is greater than a first threshold value, and a change amount of the heart rate fluctuation is greater than a second threshold value.
  • change_quantity of the said skin temperature is larger than a 3rd threshold value, you may determine with the factor of the said stress being a factor regarding an interpersonal.
  • the stress factor is an interpersonal factor by comparing the magnitude relationship between the amount of change of each biometric index and the threshold value of each biometric index.
  • the determination unit has a change amount of the heart rate larger than a first threshold value and a change amount of the heart rate fluctuation smaller than a second threshold value.
  • the amount of change in skin conductance or the amount of change in skin temperature is greater than a third threshold value, it may be determined that the stress factor is pain.
  • the stress factor is pain by comparing the magnitude relationship between the amount of change of each biological index and the threshold value of each biological index.
  • the determination unit includes a change amount of the heart rate that is smaller than a first threshold value and a change amount of the heart rate fluctuation is larger than a second threshold value.
  • the determination unit further includes a difference between the change amount of the heart rate and the first threshold value, a change amount of the heart rate fluctuation, and the second threshold value. And the difference between the skin conductance change amount or the skin temperature change amount and the third threshold value, the strength of the stress is determined, and the determination result is used as the information based on the determination result. It may be output.
  • the person being measured can know the strength of his / her stress. Thereby, it becomes easy to have consciousness about control of stress, and it becomes easy to grasp
  • the person to be measured can recognize that the strength of the stress that can be withstood among a plurality of types of stress is different. As a result, the person to be measured can determine whether or not stress control is required immediately in accordance with the state of stress. Therefore, the measurement subject can efficiently control the stress, and thus can continuously control the stress.
  • the stress evaluation device further includes a presentation unit that presents the information based on the determination result output by the determination unit, wherein the information includes the stress factor, the stress It may include at least one selected from the group consisting of strength and stress reduction measures.
  • the presenting unit may present by voice.
  • the presenting unit may present an image.
  • the person being measured can visually know the state of stress and the control method of his / her own stress, so that he / she can be clearly aware of the control of his / her stress. Therefore, the measurement subject can continue to control his own stress.
  • the obtaining step further obtains at least one of the skin conductance or skin temperature of the subject
  • the calculating step further comprises (iii) skin conductance.
  • the amount of change in skin temperature or the amount of change in skin temperature is calculated, and the amount of change in skin conductance is the skin conductance measured by the second sensor unit from a reference that is the skin conductance at rest of the subject.
  • the amount of change in the skin temperature is the skin temperature measured by the second sensor unit from a reference that is the skin temperature of the subject at rest, and the determination step includes: (I) and (II) and (III) the amount of change in the skin conductance or the amount of change in the skin temperature and the third threshold value By comparing the engagement, the determining factors of stress of the subject, and outputs the information based on the determination result.
  • the factor of stress can be determined by comparing the magnitude relationship between the amount of change of each biological index and the threshold value of each biological index.
  • FIG. 11 is a schematic configuration diagram of a stress evaluation apparatus 100a according to the present embodiment.
  • the stress evaluation apparatus 100a includes a first sensor unit 11a, a second sensor unit 11b, a calculation unit 12a, a determination unit 13a, a presentation unit 14a, and a storage unit 15a.
  • the first sensor unit 11a and the second sensor unit 11b respectively include wearable first and second biological sensors 111a and 111b (see FIG. 12) that measure the biological signal of the measurement subject. Including.
  • the first sensor unit 11a calculates a plurality of types of biological indices from the biological signal measured by the first biological sensor 111a, and outputs the calculated biological indices to the computing unit 12a.
  • the second sensor unit 11b calculates at least one type of biological index from the biological signal measured by the second biological sensor 111b, and outputs the calculated biological index to the computing unit 12a.
  • the computing unit 12a calculates an average value (hereinafter also referred to as a reference value) of each biological index when the measurement subject is at rest and a threshold value of each biological index, and stores the calculated values in the storage unit 15a.
  • the calculating part 12a calculates the average value of each measured biomarker and the variation
  • the determination unit 13a determines a factor of stress of the measurement subject based on the amount of change of each biological index. More specifically, the determination unit 13a determines the cause of stress by comparing the magnitude relationship between the amount of change of each biometric index and the threshold value of each biometric index. Further, the determination unit 13a determines the strength of the stress according to the difference between the change amount of each biometric index and the threshold value of each biometric index. And the determination part 13a outputs the information based on these determination results to the presentation part 14a.
  • the determination unit 13a stores information based on the determination result in the storage unit 15a.
  • the presentation unit 14a presents information based on the determination result.
  • the stress evaluation apparatus 100a may include an input unit 16a (see FIG. 12) for inputting an instruction of the measurement subject (user).
  • the determination unit 13a causes the presentation unit 14a to present determination result information based on the measurement subject's instruction input to the input unit 16a.
  • FIG. 12 is a block diagram showing a specific example of a stress evaluation apparatus based on the configuration of FIG.
  • the stress evaluation apparatus 100a includes a first sensor unit 11a including a first biological sensor 111a and a first signal processing unit 112a, a second biological sensor 111b, and a second signal processing unit 112b.
  • storage part 15a, and the input part 16a are provided.
  • the first biological sensor 111a and the second biological sensor 111b measure a biological signal of the measurement subject.
  • a biological signal is a signal of biological information.
  • the biological information is physiological information that is affected by stress such as heart rate, pulse, respiratory rate, blood oxygen saturation, blood pressure, or body temperature. From the ease of measurement, the biological information is, for example, heart rate information.
  • Heartbeat information is information obtained from a heartbeat.
  • the biological information may be pulse wave information.
  • the first biological sensor 111a and the second biological sensor 111b use sensors corresponding to each biological information.
  • the biological sensor is a sensor (heart rate sensor) that acquires heart rate information
  • the heart rate sensor is, for example, a sensor including a pair of detection electrodes that are in contact with the surface of the body of the measurement subject.
  • the heartbeat information obtained by the heartbeat sensor is an electrical signal obtained by the heartbeat, for example, an electrocardiogram.
  • the heart rate sensor may be a conductive adhesive gel electrode or a dry electrode composed of conductive fibers.
  • the wearing part of the heart rate sensor is the chest, and the shape of the heart rate sensor is, for example, a wear type in which wear and electrodes are integrated.
  • the pulse wave sensor measures, for example, a change in blood volume in the blood vessel using reflected light or transmitted light using a phototransistor and a photodiode. It is a sensor.
  • the pulse wave sensor is attached to the user's wrist and measures pulse wave information with the attached shape.
  • the attachment site of the pulse wave sensor may be an ankle, a finger, an upper arm, or the like.
  • the shape of the pulse wave sensor is not limited to a band type (for example, a wrist watch type), and may be a pasting type attached to a neck or the like, an eyeglass type, or the like.
  • the pulse wave sensor may be an image sensor that calculates pulse by measuring pulse wave information from a change in chromaticity of skin such as a face or a hand.
  • the biological sensor is, for example, a belt-type sensor having a pressure sensor wound around the chest or abdomen, or a temperature sensor attached under the nose.
  • the biological sensor measures, for example, a change in saturated oxygen concentration contained in blood in the blood vessel by reflected light or transmitted light using a phototransistor and two types of photodiodes. It is a sensor.
  • the biological sensor is, for example, a sensor that wraps a belt with a pressure sensor around the upper arm and a fingertip or a rib.
  • the biological sensor is a thermocouple sensor that is attached to a site where capillary contraction is likely to occur due to stress, such as the palm or the head of the nose.
  • the biological sensor is a sensor including a pair of detection electrodes that come into contact with a site where sweating is likely to occur due to stress such as a palm or a face.
  • Each biological signal measured by the first biological sensor 111a and the second biological sensor 111b is output to the first signal processing unit 112a and the second signal processing unit 112b.
  • the first signal processing unit 112a calculates a plurality of types of biological indices from one biological signal measured by the first biological sensor 111a.
  • the first sensor 111a is a heart rate sensor.
  • the plurality of types of biomarkers are RRI, CvRR, HF, LF, and the like.
  • RRI is a heart rate index
  • CvRR, HF, and LF are heart rate fluctuation indexes.
  • the first signal processing unit 112a may calculate a biological index of fluctuations in respiratory rate and blood pressure from the frequency component of heartbeat fluctuation.
  • RRI and CvRR are combinations having relatively high determination accuracy among these plural types of biological indices.
  • the first signal processing unit 112a outputs the calculated biometric index 1 and biometric index 2 to the calculation unit 12a.
  • the second signal processing unit 112b calculates at least one type of biological index from one piece of biological information measured by the second biological sensor 111b.
  • the biological index 3 is calculated.
  • the second biological sensor 111b is a sensor including a pair of detection electrodes.
  • the second biological sensor 111b is, for example, a thermocouple sensor.
  • the second biological sensor 111b is attached by being wound around the finger of the measurement subject.
  • the second signal processing unit 112b calculates skin conductance.
  • the biological index 3 is skin conductance or skin temperature.
  • the second signal processing unit 112b outputs the calculated biological index 3 to the calculation unit 12a.
  • the calculation unit 12a acquires the biometric index 1 and the biometric index 2 output from the first signal processing unit 112a, and calculates the change amount of the biometric index 1 and the change amount of the biometric index 2 from the acquired biometric index 1 and biometric index 2. To do. In addition, the calculation unit 12a acquires the biometric index 3 output from the second signal processing unit 112b, and calculates a change amount of the biometric index 3 from the acquired biometric index 3.
  • the amount of change of the biometric index is a measured biometric index based on a biometric index (hereinafter sometimes referred to as a reference value) measured when the measurement subject is at rest, and is represented by a difference or a ratio.
  • the reference value of each biometric index is stored in the storage unit 15a.
  • the calculation unit 12a reads the reference value of each biometric index stored in the storage unit 15a, and calculates the change amount of each biometric index with respect to the reference value.
  • the calculation unit 12a outputs the calculated change amount of each biological index to the determination unit 13a.
  • the reference value may vary depending on the season or the physiological period of the measurement subject, and may be updated every predetermined period.
  • the calculation unit 12a calculates a threshold value for each biological index.
  • the biological index 1 is, for example, a heart rate
  • the amount of change in heart rate is the amount of change in heart rate measured at the first time.
  • the first threshold value is a threshold value of the biometric index 1, for example, an RRI threshold value that is an index of heart rate.
  • the first threshold is a heart rate measured at an arbitrary time on the basis of the heart rate at rest of the measurement subject.
  • the biological index 2 is, for example, heart rate fluctuation
  • the amount of change in heart rate fluctuation is the amount of change in heart rate fluctuation measured at the second time.
  • the second threshold value is a threshold value of the biological index 2, for example, a CvRR threshold value that is an index of heartbeat fluctuation.
  • the second threshold is a heart rate fluctuation measured at an arbitrary time on the basis of the heart rate of the subject at rest.
  • the biological index 3 is, for example, skin conductance or skin temperature
  • the amount of change in skin conductance or skin temperature is the amount of change in skin conductance or skin temperature measured at the third time.
  • the third threshold is a threshold of the biometric index 3, and is, for example, a threshold of skin conductance or a threshold of skin temperature.
  • the third threshold value is the skin conductance measured at any time based on the skin conductance of the subject at rest, or at any time based on the skin temperature of the subject at rest. The measured skin temperature.
  • These threshold values are changes in the biometric index that is the difference or ratio between the measured value of the biometric index measured at an arbitrary time different from the first, second, and third times and the reference value.
  • the arbitrary time indicates, for example, the time when the measurement subject is in a state before feeling stress.
  • the first time, the second time, and the third time are the same time.
  • the times may be different times.
  • the first signal processing unit 112a may calculate a plurality of types of heart rates and heart rate fluctuations in a time division manner from one biological signal measured by the first biological sensor 111a.
  • the calculation unit 12 calculates the amount of change in heart rate measured at the first time, and calculates the amount of change in heart rate fluctuation measured at a second time different from the first time.
  • the second signal processing unit 112b may measure sweat or skin temperature with the second biological sensor 112b at the third time. At this time, the calculation unit 12 calculates the amount of change in skin conductance or the amount of change in skin temperature measured at the third time.
  • the third time may be the same time as either the first time or the second time.
  • the calculation unit 12a reads the threshold value of each biological index stored in the storage unit 15a, and compares the magnitude relationship between the change value of each biological index and the threshold value of each biological index. And the calculating part 12a determines the period when at least 1 of the variation
  • the stress generation period is a period during which the measurement subject feels stress.
  • the calculation unit 12a calculates a representative value of the change amount of each biological index from the value of the change amount of each biometric index during the stress occurrence period.
  • the average value of the change amount of each biometric index during the stress occurrence period may be used, and the value having the largest difference from the reference value (maximum value) May be used.
  • the determination unit 13 acquires a representative value of the change amount of each biological index output by the calculation unit 12a, and reads the first threshold value, the second threshold value, and the third threshold value stored in the storage unit 15a.
  • the determination unit 13a compares the magnitude relationship between the representative value of the change amount of the biometric index 1 and the first threshold value, and compares the magnitude relationship between the representative value of the change amount of the biometric index 2 and the second threshold value.
  • the factor of the stress of the person to be measured is determined by comparing the magnitude relationship between the representative value of the change amount of the biological index 3 and the third threshold value. That is, the determination unit 13a determines a stress factor for each stress occurrence period. Since the representative value of the change amount of the biometric index is an example of the change amount of the biometric index, hereinafter, the representative value of the change amount of the biometric index is also simply referred to as a change amount of the biometric index.
  • the determination unit 13a has a change amount of the biometric index 1 (here, heart rate) larger than the first threshold value, and a change amount of the biometric index 2 (here, heart rate fluctuation) is the second. And the change amount of the biometric index 3 (here, skin conductance or skin temperature) is greater than the third threshold value, it is determined that the stress factor is a factor related to the person. In addition, the determination unit 13a determines that the change amount of the biometric index 1 is larger than the first threshold value, the change amount of the biometric index 2 is smaller than the second threshold value, and the change amount of the biometric index 3 is the third change amount. If the threshold value is larger than the threshold value, it is determined that the stress factor is pain.
  • the determination unit 13a determines that the change amount of the biometric index 1 is smaller than the first threshold value, the change amount of the biometric index 2 is larger than the second threshold value, and the change amount of the biometric index 3 is the third change amount. If it is smaller than the threshold value, it is determined that the stress factor is fatigue due to thought.
  • the determination unit 13a includes the difference between the change amount of the biometric index 1 and the first threshold value, the difference between the change amount of the biometric index 2 and the second threshold value, and the change amount of the biometric index 3 and the third threshold value.
  • the strength of the stress is determined according to the difference between and the determination result is output as information based on the determination result.
  • the information based on the determination result includes, for example, at least one of a stress factor, a stress intensity, and a stress reduction measure.
  • the stress reduction measure is, for example, a stress elimination method or a stress avoidance method.
  • the stress reduction measure is included in a presentation information table described later.
  • the determination unit 13a reads out an appropriate stress reduction measure from the presentation information table stored in the storage unit 15a and outputs it to the presentation unit 14a.
  • the determination unit 13a stores information based on the determination result in the storage unit 15a. At this time, the determination unit 13a may associate the information on the time when the measurement subject feels stress with the information based on the determination result and store the information in the storage unit 15a.
  • the presentation unit 14a presents information based on the determination result output by the determination unit 13a.
  • the presentation unit 14a may present information based on the determination result by voice or an image.
  • the presentation unit 14a is, for example, a speaker.
  • the presentation part 14a presents the said information with an image, the presentation part 14a is a display, for example.
  • the storage unit 15a stores a reference value of each biological index, a threshold value of each biological index, a presentation information table, and the like.
  • the presentation information table is a table of presentation information such as stress reduction measures presented according to the stress factor and the intensity of the stress. As described above, the reference value and threshold value of each biometric index may be updated in a predetermined period. Similarly, the presentation information table may be updated in a predetermined period.
  • the storage unit 15a stores information based on the determination results such as the stress factor, the stress intensity, and the stress reduction measure output by the determination unit 13a. At this time, the storage unit 15a may store information based on the determination result and the stress occurrence period in association with each other. Thereby, the measurement subject can call information based on the determination result at a desired timing. At this time, the determination unit 13a causes the presentation unit 14 to present information based on the determination result based on the measurement subject's operation input by the input unit 16a.
  • the input unit 16a outputs an operation signal indicating an operation by the measurement subject to the determination unit 13a.
  • the input unit 16a is, for example, a keyboard, a mouse, a touch panel, or a microphone.
  • the operation signal is a signal for setting an information extraction method based on the determination result or a presentation method in the presentation unit 14a.
  • the presentation unit 14a presents various types of determination results based on the settings input to the input unit 16a. For example, a change in stress during a predetermined period, a factor of stress that the measurement subject is likely to be influenced by, a stress reduction measure suitable for the measurement subject, and the like. As a result, the measurement subject can grasp not only the short-term stress tendency but also the medium-term and long-term stress tendency. In this way, the person to be measured can know an effective stress reduction measure suitable for himself / herself, and therefore can control mid- to long-term stress.
  • FIG. 13 is a flowchart illustrating the stress evaluation method according to the embodiment.
  • the stress evaluation method includes (i) heart rate, (ii) heart rate fluctuation, and (iii) skin conductance or skin temperature of the measured person to be measured.
  • Calculation step S200 for calculating (i) heart rate change amount, (ii) heart rate fluctuation change amount, and (iii) skin conductance change amount or skin temperature change amount; and (i) heart rate change amount.
  • the amount of change in heart rate is the amount of change to the heart rate measured by the first sensor unit 11a from the reference, which is the heart rate at rest of the subject, and the amount of change in heart rate fluctuation is the amount of change in the subject. This is the amount of change to heartbeat fluctuation measured by the first sensor unit 11a from the reference, which is heartbeat fluctuation at rest.
  • the amount of change in skin conductance is the amount of change in skin conductance measured by the second sensor unit 11b from the reference, which is the skin conductance at rest of the subject, and the amount of change in skin temperature is It is the skin temperature measured by the 2nd sensor part 11b from the reference
  • the determination step S300 (I) the magnitude relationship between the heart rate change amount and the first threshold value is compared, and (II) the magnitude relationship between the heart rate fluctuation change amount and the second threshold value is compared, and (III) The factor of the stress is determined by comparing the magnitude relationship between the change amount of the skin conductance or the change amount of the skin temperature and the third threshold value.
  • the present embodiment further includes a presentation step S400 for presenting information based on the determination result of determination step S300.
  • the calculation unit 12a acquires a plurality of biological indices of the measurement subject measured by the first sensor unit 11a and the second sensor unit 11b.
  • heart rate information here, an electrocardiogram
  • a heart rate index and a heart rate fluctuation index are calculated by the first signal processing unit 112a.
  • the biological information on the temperature or sweating is measured by the second biological sensor 111b, and the skin temperature (SKT) or skin conductance (SC) is calculated by the second signal processing unit 112b.
  • the biological information may be physiological information that is affected by stress such as heart rate, pulse rate, respiratory rate, blood oxygen saturation, blood pressure, body temperature, and sweating.
  • the heart rate information is simple in a state where the burden on the subject is less than other biological information such as pulse, respiratory rate, blood pressure, and blood oxygen saturation, and It can be measured in real time. Therefore, by using the heartbeat information of the person to be measured as the biological information, it is possible to appropriately evaluate the stress state of the person to be measured.
  • the biological index obtained from the heart rate information includes RRI that is an index of heart rate, CvRR that is an index of heart rate fluctuation, LF, HF, LF / HF, and the like.
  • RRI an index of heart rate
  • CvRR an index of heart rate fluctuation
  • LF HF
  • LF / HF a plurality of types of biological indices
  • the combination of these biomarkers can determine the cause of stress with relatively high determination accuracy, so that highly reliable evaluation can be obtained.
  • the heart rate information is, for example, an electrocardiogram and has an electrocardiogram waveform shown in FIG.
  • the ECG waveform reflects the P wave that reflects the electrical excitation of the atrium, the Q wave, the R wave, and the S wave that reflect the electrical excitation of the ventricle, and the process of repolarization of the excited ventricular cardiomyocytes. It consists of T waves.
  • the wave height (potential difference) of the R wave is the largest and is most robust against noise such as myoelectric potential. For this reason, the interval between the R wave peaks of two consecutive heartbeats in these electrocardiographic waveforms, that is, the heartbeat interval (RRI) is calculated.
  • the heart rate is calculated by multiplying the reciprocal of RRI by 60.
  • CvRR is obtained by normalizing the standard deviation SD of RRI in an arbitrary time zone from the RRI by the average value of the heartbeat interval using the above formula (2). Calculated.
  • the first signal processing unit 112a detects an electrical signal (R wave) generated when the left ventricle contracts suddenly and pumps blood from the heart from the heartbeat information obtained by the first biological sensor 111a, and calculates the RRI. calculate.
  • R wave an electrical signal generated when the left ventricle contracts suddenly and pumps blood from the heart from the heartbeat information obtained by the first biological sensor 111a
  • calculates the RRI. calculate for detecting the R wave, for example, a known method such as the Pan & Topkins method is used.
  • the first signal processing unit 112a calculates the RRI fluctuation amount from the obtained R-wave detection data as follows.
  • the first signal processing unit 112a calculates the RRI that is the interval between the R wave peaks of two consecutive heartbeats.
  • the first signal processing unit 112a converts each calculated RRI into a biaxial relationship between time and RRI. Since the converted data is discrete data with unequal intervals, the arithmetic unit 12a converts the converted RRI time-series data into the equally-spaced time-series data shown in FIG. 7B.
  • the computing unit 12a obtains a frequency component of heartbeat variability shown in FIG. 7C by performing frequency analysis on the equidistant time series data using, for example, fast Fourier transform (FFT).
  • FFT fast Fourier transform
  • the frequency component of heart rate variability can be divided into, for example, a high frequency component HF and a low frequency component LF.
  • HF is considered to reflect the amount of parasympathetic nerve activity.
  • LF is also considered to reflect the amount of sympathetic and parasympathetic activity. Therefore, LF / HF, which is the ratio of LF to HF, is considered to indicate the amount of sympathetic nerve activity.
  • a plurality of types of biological indices are calculated from the heartbeat information.
  • the calculation unit 12a outputs the two types of biological indexes (here, heart rate and heart rate fluctuation) output from the first sensor unit 11a and the second sensor unit 11b.
  • One type of biological index here, skin conductance
  • the calculation unit 12a calculates the amount of change of each biological index acquired in acquisition step S100.
  • the amount of change of each biological index is, for example, the ratio between the reference value of each biological index and the value of each acquired biological index, with the value of each biological index resting on the measurement subject as a reference value. Obtained by calculating the difference.
  • the calculation unit 12a reads out and uses the reference value of each biological index stored in the storage unit 15a.
  • the change amount of each biometric index is calculated by subtracting each biometric index reference value from the value of each biometric index acquired in the acquisition step S100 when the change amount is represented by a difference.
  • the change amount of the heart rate is calculated by subtracting the reference value of the heart rate from the value of the heart rate of the measurement subject acquired in the acquisition step S100.
  • the amount of change is expressed as a ratio, it is calculated by dividing the value of each biological index acquired in acquisition step S100 by the reference value of each biological index.
  • the amount of change in heart rate is calculated by dividing the value of the heart rate of the measurement subject acquired in acquisition step S100 by the reference value of the heart rate.
  • the change amount of each biological index is calculated by the calculation unit 12a.
  • the determination unit 13a determines a stress factor based on the amount of change of each biological index calculated in calculation step S200.
  • the determination unit 13a determines the factor of stress of the measurement subject by comparing the magnitude relationship between the change amount of each biometric index and the threshold value of each biometric index. Specifically, in the determination step S300, the determination unit 13a determines that the change amount of the heart rate is greater than the first threshold value, the change amount of the heart rate fluctuation is greater than the second threshold value, and the skin conductance. When the amount of change or the amount of change in skin temperature is greater than the third threshold, it is determined that the stress factor is a factor related to the person.
  • the determination unit 13a determines that the change amount of the biometric index 1 is larger than the first threshold value, the change amount of the biometric index 2 is smaller than the second threshold value, and the change amount of the skin conductance or the skin temperature. When the amount of change is larger than the third threshold, it is determined that the stress factor is pain. In addition, the determination unit 13a determines that the change amount of the biometric index 1 is smaller than the first threshold value, the change amount of the biometric index 2 is larger than the second threshold value, and the change amount of the skin conductance or the skin temperature. When the amount of change is smaller than the third threshold, it is determined that the stress factor is fatigue due to thought.
  • the determination unit 13a includes a difference between the change amount of the biological index 1 and the first threshold value, a difference between the change amount of the heartbeat fluctuation and the second threshold value, and a change amount of the skin conductance or a change amount of the skin temperature.
  • the strength of the stress is determined according to the difference from the third threshold value, and the determination result is output as information based on the determination result.
  • the first threshold value is a heart rate threshold value, which is a heart rate measured at any time with respect to the measurement subject based on the heart rate of the measurement subject at rest.
  • the second threshold value is a heart rate fluctuation threshold value, which is a heart rate fluctuation value measured at an arbitrary time based on the heart rate fluctuation value of the measurement subject at rest.
  • the third threshold value is a skin conductance or skin temperature threshold value, and is a skin conductance or skin temperature measured at an arbitrary time based on the skin conductance or skin temperature of the measurement subject at rest.
  • a threshold value when the amount of change of each biological index is a positive value and a threshold value when the amount of change of each biological index is a negative value are set.
  • the reference value is the zero point of change.
  • the magnitude relationship between the change amount of each biological index and the threshold value is compared as follows. When the change amount of the biometric index is a positive value, the magnitude relationship between the change amount of the biometric index and the positive threshold value is compared. When the change amount of the biometric index is a negative value, the magnitude relationship between the absolute value of the change amount of the biometric index and the absolute value of the negative threshold is compared.
  • the threshold value of each biometric index may be a fixed value, may be updated in a predetermined period, or may be updated each time based on daily measurement.
  • the threshold value may be calculated by relatively simple machine learning such as linear discrimination or decision tree. Therefore, since the determination reference value and threshold value suitable for the person to be measured can be set, the cause of stress can be determined with higher accuracy.
  • the factor of stress of the measurement subject is determined by comparing the magnitude relationship between the amount of change of each biological index and the threshold value of each biological index.
  • the presentation unit 14a presents information based on the determination result determined by the determination unit 13a.
  • the presentation unit 14a may present information based on the determination result by voice or may be presented by an image.
  • the information based on the determination result includes at least one of a stress factor, a stress intensity, and a stress reduction measure.
  • the presentation unit 14a displays various types of determination results based on the settings input by the measurement subject using the input unit 16a.
  • FIG. 14 is a diagram for explaining an example of use of the stress evaluation apparatus 100a according to the present embodiment.
  • the stress evaluation apparatus 100a includes a first biosensor 111a that is a part of the first sensor unit 11a, a second biosensor 111b that is a part of the second sensor unit 11b, and a first biosensor.
  • 111a and the evaluation terminal 20 including configurations other than the second biological sensor 111b.
  • the measurement subject wears the first living body sensor 111a so as to contact the skin of the chest, and measures an electrocardiogram (ECG).
  • ECG electrocardiogram
  • the first biological sensor 111a may be a conductive adhesive gel electrode or a dry electrode composed of conductive fibers or the like.
  • the first biological sensor 111a transmits the measured heartbeat electrical signal to the evaluation terminal 20 by communication.
  • the second biological sensor 111b is a wristwatch type sensor and includes a sensor electrode that is used by being attached to a palm.
  • the second biological sensor 111b measures the skin potential of the palm measured by the sensor electrode and transmits it to the evaluation terminal 20 by communication.
  • the second biological sensor 111b may include a thermocouple sensor that is used by being attached to a fingertip. Thereby, the 2nd biological sensor 111b can measure the temperature of a fingertip with a thermocouple type sensor.
  • the communication method between the first biological sensor 111a and the second biological sensor 111b and the evaluation terminal 20 may be wireless communication such as Bluetooth (registered trademark) or wired communication.
  • the evaluation terminal 20 includes a first signal processing unit 112a of the first sensor unit 11a, a second signal processing unit 112b of the second sensor unit 11b, a calculation unit 12a, a determination unit 13a, a presentation unit 14a, a storage unit 15a, and an input unit 16a. Is provided.
  • the 1st signal processing part 112a and the 2nd signal processing part 112b receive the living body signal transmitted by communication from the 1st living body sensor 111a and the 2nd living body sensor 111b, respectively.
  • the first signal processing unit 112a calculates RRI that is a heart rate index and CvRR that is a heart rate fluctuation index from the received heartbeat electrical signal, and outputs these biological indices to the calculation unit 12a.
  • the second signal processing unit 112b calculates skin conductance (SC) that is an index of sweating from the received skin potential signal, and outputs SC to the calculation unit 12a.
  • SC skin conductance
  • the second signal processing unit 112b receives the skin temperature signal from the second biological sensor 111b and calculates the skin temperature (SKT) that is an index of body temperature. Then, the SKT is output to the calculation unit 12a.
  • the calculation unit 12a acquires the RRI and CvRR output from the first signal processing unit 112a, and reads the RRI reference value and the CvRR reference value stored in the storage unit 15. In addition, the calculation unit 12a acquires the SC output from the second signal processing unit 112b, and reads the reference value of the SC stored in the storage unit 15a. The calculation unit 12a calculates the amount of change of each of these biometric indices with reference to the read reference value. The amount of change in the biometric index is represented by a difference or a ratio. In the present embodiment, the amount of change is expressed as a ratio.
  • the calculation unit 12a calculates the threshold value of each biological index and outputs it to the storage unit 15a.
  • a threshold value when the change amount of each biometric index becomes a positive value and a threshold value when the change amount of each biometric index becomes a negative value are set.
  • the reference value is zero change.
  • the positive threshold value is a value larger than the reference value, and the first threshold value 1a (hereinafter, positive value) in the change amount graph 120a.
  • the negative threshold value is smaller than the reference value, and the first threshold value 1b (hereinafter, negative threshold value 1b) in the change amount graph 120, the first threshold value.
  • a second threshold 2b (hereinafter, negative threshold 2b) and a third threshold 3b (hereinafter, negative threshold 3b).
  • the calculating part 12a calculates the reference value of each biological index, and outputs it to the storage part 15a.
  • the reference value of each biometric index the amount of change of each biometric index is zero.
  • the reference value is indicated by a solid line between the positive threshold value 1a and the negative threshold value 1b.
  • the positive threshold value and the negative threshold value may be set at equal intervals with the reference value (change amount zero) interposed therebetween, or may not be set at equal intervals with the reference value interposed therebetween. These threshold values may be appropriately set according to the magnitude of the change amount of each biological index.
  • the determination unit 13a acquires the amount of change of each biological index output from the calculation unit 12a and reads the threshold value of each biological index stored in the storage unit 15a.
  • the determination unit 13a compares the magnitude of the change amount of each biometric index with the threshold value of each biometric index, and determines a factor of stress. For example, when the change amount of each biological index is a positive value, the determination unit 13a compares the magnitude relationship between the change amount of each biological index and a positive threshold value. When the change amount of each biological index is a negative value, the determination unit 13a compares the magnitude relationship between the absolute value of the change amount of each biological index and the absolute value of the negative threshold value.
  • the change amount graph 120a and the determination table 130a will be described in more detail.
  • the determination unit 13a determines that the stress factor felt by the measurement subject during the period A2 is a factor related to the person.
  • the change amount of RRI is larger than the positive threshold value 1a
  • the absolute value of change amount of CvRR is smaller than the absolute value of the negative threshold value 2b
  • the change amount of skin conductance is , Greater than the positive threshold 3a.
  • the determination unit 13a determines that the stress factor felt by the measurement subject during the period B2 is pain.
  • the absolute value of the change amount of RRI is smaller than the absolute value of the negative threshold value 1b
  • the absolute value of the change amount of CvRR is larger than the absolute value of the negative threshold value 2b
  • the absolute value of the change amount of the skin conductance is smaller than the absolute value of the negative threshold 3b. Therefore, the determination unit 13a determines that the factor of the stress felt by the measurement person in the period C2 is fatigue due to thought (thinking fatigue).
  • the change amount of each biological index based on the reference value is indicated by the direction and number of arrows.
  • a horizontal arrow indicates that the change amount of the biometric index is not accompanied by a change exceeding the threshold value.
  • the determination unit 13a determines the difference between the absolute value of the RRI change amount and the absolute value of the first threshold, the difference between the absolute value of the CvRR change amount and the absolute value of the second threshold, and the change of the SC.
  • the intensity of stress is determined according to the difference between the absolute value of the quantity and the absolute value of the third threshold value.
  • the determination unit 13a outputs information based on these determination results to the presentation unit 14a.
  • the presentation unit 14a is, for example, a smartphone display.
  • the determination unit 13a stores information based on the determination result in the storage unit 15a. Thereby, the measurement subject can call information based on the determination result at a desired timing.
  • the determination unit 13a causes the presentation unit 14a to present information based on the determination result based on the operation of the measurement subject input by the input unit 16a such as a touch panel. For example, when the measurement subject inputs an instruction to extract necessary information at the input unit 16a of the evaluation terminal 20, the determination unit 13a presents the presentation information 140a to the presentation unit 14a based on the instruction of the measurement subject.
  • the presentation information 140a includes a time when the measurement subject feels stress, a factor of stress, and a measure for reducing stress.
  • the stress reduction measure is, for example, a message for proposing a stress relieving method or a stress avoiding method according to a stress factor.
  • the message is, for example, if the stress factor is thought fatigue, let's take a break or stretch, and if it's a factor related to the person, meditate a little or take a deep breath. And so on.
  • the measurement subject can grasp his / her stress state and appropriate stress reduction measures more accurately than in the past. As a result, the measurement subject can appropriately and efficiently control their own stress, and therefore can continuously control the stress.
  • heartbeat information is used as biological information
  • a heart rate index and a heartbeat fluctuation index are used as a plurality of types of biological indices obtained from the heartbeat information.
  • entropy E that is the degree of autonomic nerve activity
  • tone T that is the autonomic nerve balance
  • the stress evaluation apparatus 100 showed the example comprised from the biosensor 111 and the evaluation terminal 20 in Embodiment 1, it comprises the structure other than the 1st sensor part 11a and the 1st sensor part 11a, for example. You may comprise from an evaluation terminal.
  • the stress evaluation apparatus 100a showed the example comprised from the biosensor 111a and the evaluation terminal 20, the 1st sensor part 11a, the 2nd sensor part 11b, and the 1st sensor part were shown, for example. You may comprise from the evaluation terminal provided with structures other than 11a and the 2nd sensor part 11b.
  • the stress evaluation apparatus may be an integrated apparatus in which all components are incorporated into one device.
  • the biological sensor is a heartbeat sensor
  • the biological sensor may be a pulse wave sensor.
  • the stress evaluation apparatus may be a wristwatch type apparatus including a display.
  • the evaluation terminal 20 showed the example which is a smart phone or a tablet terminal
  • the smart phone or the tablet terminal is provided with the presentation part 14 and the input part 16, and the 1st signal processing part 112a and the calculating part 12 are provided.
  • the determination unit 13 and the storage unit 15 may be provided on a server connected via a communication network such as the Internet.
  • the evaluation terminal 20 showed the example which is a smart phone or a tablet terminal
  • the smart phone or the tablet terminal is provided with the presentation part 14a and the input part 16a, the 1st signal processing part 112a, the 2nd signal You may provide the process part 112b calculating part 12a, the determination part 13a, and the memory
  • the reference value and threshold value of each biometric index are shown as an example of a form stored in the storage unit provided in the evaluation terminal, the reference value and threshold value are stored in a server on the Internet and sent to the evaluation terminal as needed. It may be in the form that comes.
  • skin conductance is cited as one of the indices for determining the factor of stress, but it is not particularly limited as long as it is an index that can measure mental sweating.
  • it may be an index obtained by measuring a skin potential or current value such as skin resistance, or may be an index obtained by measuring a moisture content such as humidity on the skin surface.
  • the skin conductance or skin temperature is measured with the palm. However, it may be measured with a part of the face where mental sweating is likely to occur or may be measured with the sole of the foot. .
  • a simulated job interview in the monitor test is given as a specific example of factors related to interpersonalness, which is one of the factors of stress, but is not limited thereto.
  • the interpersonal factor may be a factor that causes the person to be measured to feel anxiety or tension in matters related to people, such as human relations at work and private, talking in front of people, or negotiating with people.
  • the pain may be pain that causes fear or patience, such as physical pain such as bruise, headache, toothache, laceration, or pain associated with physical stimulation such as rubbing, stab, cut, and strike.
  • fatigue due to thought which is one of the factors of stress
  • the task of mental arithmetic and voice-based soap is cited as a task that requires thought, but the present invention is not limited thereto.
  • fatigue due to thinking may be a factor that causes fatigue by continuing thinking work such as work on a personal computer or intellectual activity such as experiments that require concentration as work that requires thought. .
  • the present disclosure is useful as a stress evaluation apparatus that can easily and accurately determine a stress factor of a measurement subject from a plurality of types of biological indices.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Psychiatry (AREA)
  • Physiology (AREA)
  • Child & Adolescent Psychology (AREA)
  • Developmental Disabilities (AREA)
  • Educational Technology (AREA)
  • Hospice & Palliative Care (AREA)
  • Signal Processing (AREA)
  • Psychology (AREA)
  • Social Psychology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

ストレス評価装置１００は、被測定者の心拍数及び心拍揺らぎを測定する第１センサ部１１ａと、（ｉ）心拍数の変化量及び（ｉｉ）心拍揺らぎの変化量を算出する演算部１２と、（ｉ）心拍数の変化量及び（ｉｉ）心拍揺らぎの変化量に基づいて被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定部１３と、を備える。心拍数の変化量は、被測定者の安静時の心拍数である基準からの第１センサ部１１ａによって測定された心拍数への変化量である。心拍揺らぎの変化量は、被測定者の安静時の心拍揺らぎである基準からの第１センサ部１１ａによって測定された心拍揺らぎへの変化量である。判定部１３は、（Ｉ）心拍数の変化量と第１の閾値との大小関係の比較及び（ＩＩ）心拍揺らぎの変化量と第２の閾値との大小関係の比較を行うことによりストレスの要因を判定する。

Description

ストレス評価装置、ストレス評価方法及びプログラム
　本開示は、被測定者のストレスの要因を判定するストレス評価装置、ストレス評価方法及びプログラムに関する。
　近年のウエアラブルデバイスの発達により、日常生活での生体指標の測定が可能な生体指標測定装置が普及している。例えば、ストレスの評価用デバイスでは、デバイスに搭載した加速度センサによって被測定者の動きを検知し、安静時のストレス測定を行うことが試みられている。
　例えば、特許文献１は、加速度センサの検出値を基に被測定者の活動強度等を算出し、心拍数、鼓動波形、血圧、血中酸素飽和度、体温、又は発汗度などの生体指標及び活動強度に基づいて、被測定者のストレス状態を判定できるシステムを開示している。
　また、特許文献２は、被測定者の生体指標及び行動情報に基づいて被測定者のストレス状態を周囲の状況と併せて分析して判断することにより、被測定者にストレス解消方法などを提供する生活支援装置及び生活支援方法を開示している。
特開２００９－１４８３７２号公報 特開２００１－３４４３５２号公報
　本開示は、被測定者のストレスの要因を判定することができるストレス評価装置、ストレス評価方法及びプログラムを提供する。
　本開示の一態様に係るストレス評価装置は、被測定者の心拍数及び心拍揺らぎを測定する第１センサ部と、（ｉ）心拍数の変化量、及び、（ｉｉ）心拍揺らぎの変化量を算出する演算部と、（ｉ）前記心拍数の変化量及び（ｉｉ）前記心拍揺らぎの変化量に基づいて前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定部と、を備え、前記心拍数の変化量は、前記被測定者の安静時の心拍数である基準からの前記第１センサ部によって測定された前記心拍数への変化量であり、前記心拍揺らぎの変化量は、前記被測定者の安静時の心拍揺らぎである基準からの前記第１センサ部によって測定された前記心拍揺らぎへの変化量であり、前記判定部は、（Ｉ）前記心拍数の変化量と第１の閾値との大小関係の比較、及び、（ＩＩ）前記心拍揺らぎの変化量と第２の閾値との大小関係の比較を行うことにより、前記ストレスの要因を判定する。
　また、本開示の一態様に係るストレス評価方法は、測定された被測定者の心拍数及び心拍揺らぎを取得する取得ステップと、（ｉ）心拍数の変化量、及び、（ｉｉ）心拍揺らぎの変化量を算出する算出ステップと、前記心拍数の変化量及び前記心拍揺らぎの変化量に基づいて前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定ステップと、を含み、前記心拍数の変化量は、前記被測定者の安静時の心拍数である基準からの前記第１センサ部によって測定された前記心拍数への変化量であり、前記心拍揺らぎの変化量は、前記被測定者の安静時の心拍揺らぎである基準からの前記第１センサ部によって測定された前記心拍揺らぎへの変化量であり、前記判定ステップでは、（Ｉ）前記心拍数の変化量と第１の閾値との大小関係を比較し、かつ、（ＩＩ）前記心拍揺らぎの変化量と第２の閾値との大小関係を比較することにより、前記ストレスの要因を判定する。
　なお、これらの包括的又は具体的な態様は、システム、装置、集積回路、コンピュータプログラム又はコンピュータで読み取り可能なＣＤ－ＲＯＭなどの記録媒体で実現されてもよく、システム、装置、集積回路、コンピュータプログラム及び記録媒体の任意な組み合わせで実現されてもよい。
　本開示のストレス評価装置、ストレス評価方法及びプログラムによれば、被測定者のストレスの要因を評価することができる。
図１は、被験者２０名それぞれのストレス要因毎の生体指標の変化量をプロットした図である。 図２は、図１でプロットされたストレスの要因毎の生体指標の変化量の平均値を示す図である。 図３は、実施の形態１に係るストレス評価装置の構成の一例を示す概略構成図である。 図４は、図３の構成に基づくストレス評価装置の具体例を示す構成図である。 図５は、実施の形態１に係るストレス評価方法を説明するフローチャートである。 図６は、実施の形態１に係るストレス評価装置で得られる心拍情報の一例を示す図である。 図７は、心拍間隔（ＲＲＩ）の変動量を算出する手法を説明する図である。 図８は、実施の形態１に係るストレス評価装置の使用例を説明する図である。 図９Ａは、被験者２０名それぞれのストレス要因毎の生体指標の変化量をプロットした図である。 図９Ｂは、図９ＡをＲＲＩの変化量を示す軸のプラス側から見た図である。 図９Ｃは、図９ＡをＣｖＲＲの変化量を示す軸のマイナス側から見た図である。 図９Ｄは、図９ＡをＳＣの変化量を示す軸のマイナス側から見た図である。 図１０Ａは、図９Ａでプロットされたストレスの要因毎の生体指標の変化量の平均値を示す図である。 図１０Ｂは、図１０ＡをＲＲＩの変化量を示す軸のプラス側から見た図である。 図１０Ｃは、図１０ＡをＣｖＲＲの変化量を示す軸のマイナス側から見た図である。 図１０Ｄは、図１０ＡをＳＣの変化量を示す軸のマイナス側から見た図である。 図１１は、実施の形態に係るストレス評価装置の構成の一例を示す概略構成図である。 図１２は、図１１の構成に基づくストレス評価装置の具体例を示す構成図である。 図１３は、実施の形態２に係るストレス評価方法を説明するフローチャートである。 図１４は、実施の形態２に係るストレス評価装置の使用例を説明する図である。
　（本開示の基礎となった第１の知見）
　現代社会における鬱等のストレス障害は、日々の生活で蓄積されたストレスが原因で重症化することが多い。このような問題を回避するために、日常生活の中でストレスの蓄積を低減させることが重要となる。つまり、人々が自身のストレス状態をコントロールできることが望ましい。そのために、日常生活におけるストレスの状態をセンシングして、ストレスの強度及びストレスの要因に応じて適切なストレス解消方法及びストレス回避方法などのストレス低減策をユーザに提供することが望ましい。
　例えば、特許文献１に記載のストレス判定システムは、加速度センサから得られた情報を基に被測定者の活動強度等を算出し、心拍数、鼓動波形、血圧、血中酸素飽和度、体温、発汗度などの生体指標及び活動強度に基づいて、被測定者のストレス状態を判定する。当該システムでは、活動強度が一定値以下の場合にのみ生体指標を測定することにより、被測定者の日常生活におけるストレス状態を判定している。
　しかしながら、特許文献１に記載のストレス判定システムでは、ストレスの有無の判定は可能であるが、ストレスの要因についての情報が得られない。人がストレスを受ける要因、つまり、ストレスの要因は様々である。また、ストレスの要因に応じて最適なストレス解消方法及びストレス回避方法は異なる。特許文献１に記載のストレス判定システムでは、ストレスの要因についての情報が得られないため、ユーザに適切なストレス解消方法及びストレス回避方法を提供することができず、ユーザのストレスの制御を行うには不十分である。
　また、特許文献２に記載の生活支援システムは、心電及び脈波等の生体情報だけではなく、被測定者の行動情報を取得して、被測定者の周囲の状況を分析して判断することにより、被測定者にストレス解消方法などを提供する。
　しかしながら、特許文献２に記載の生活支援システムでは、被測定者の周囲の状況が同じであっても、被測定者によってストレスの要因が異なる場合があるため、被測定者が実際に感じているストレスの要因を判定することは難しい。そのため、特許文献２に記載の生活支援システムでは、被測定者に不適切なストレス解消方法及びストレス対処行動を提示する危険性がある。
　本発明者らは、上記課題に鑑みて鋭意検討をした。検討内容を以下に記す。
　本発明者らは、ストレスの要因と、心拍情報などの生体情報から得られる複数種類の生体指標との関連性を見出すために、以下のモニター試験を実施した。
　［モニター試験］
　２０名の被験者に対してストレスの要因が異なる４つのタスクを与え、タスクを実行している被験者の生体信号を測定した。
　被験者は、２０代から３０代の社会人又は大学生の男女であり、健康状態及び精神状態に関するアンケ―トの結果が異常値を示さない２０名が選出された。
　タスクは、［１］対人に関するストレス、［２］痛みに関するストレス、［３］思考による疲労（以下、思考疲労）に関するストレス１、［４］思考疲労に関するストレス２の４種類である。各タスクは各被験者に対して個別に実施された。タスクの詳細は、以下の通りである。
　［１］対人に関するストレス
　被験者と初対面の男性１名及び女性１名の合計２名のタスク説明者が、被験者にタスクの説明をした後、被験者にタスクを実行させ、タスク実行中の被験者の生体信号を測定した。具体的には、タスク説明者は、５分間後に模擬の就職面談を行うこと、及び、面談開始までの５分間で話す内容を決定することを被験者に伝えた。生体信号の測定は、会話による動き及びノイズを考慮し、被験者が話す内容を考える５分間に実施した。
　［２］痛みに関するストレス
　被験者が十分に痛みを感じる程度に調整した電気刺激を、被験者の前腕部に１０分間与えた。電気刺激は、約１分間に、ランダムに１０回程度実施した。これを１０分間繰り返した。生体信号の測定は、電気刺激を開始してから最初の５分間に実施した。
　［３］思考疲労に関するストレス１
　ディスプレイに表示された２桁又は３桁の掛け算の問題を制限時間内に被験者に解答させた。被験者は、掛け算の問題を暗算し、ディスプレイに表示された３つの選択肢から解答を選択した。問題の難易度及び１問当たりの制限時間は、被験者の暗算能力を事前に測定することにより、決定された。被験者は、このタスクを１５分間実行した。生体信号の測定は、被験者がタスクを開始してから最初の５分間に実施した。
　［４］思考疲労に関するストレス２
　スピーカーから指示されるじゃんけんの問題に対してディスプレイに表示された３つの選択肢から正しいものを制限時間内に被験者に選択させた。１問当たりの制限時間は、被験者の解答能力を事前に測定することにより、決定された。被験者は、このタスクを１５分間実行した。生体信号の測定は、被験者がタスクを開始してから最初の５分間に実施した。
　上記のモニター試験は、日内変動を考慮し、被験者毎に別日の同時刻に実施した。
　被験者の安静時の生体信号は、上記［１］～［４］の各タスクを実施する前に、タスクを実行する姿勢と同じ姿勢で、５分間測定した生体信号である。この生体信号から生体指標を算出し、生体指標の変化量を算出するための基準値とした。生体指標の変化量は、被験者の安静時の生体指標を基準とするタスク実行中に測定された被験者の生体信号から算出された生体指標である。
　測定された生体信号は、心電図（Ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｍ：ＥＣＧ）、呼吸間隔、指先温度（Ｓｋｉｎ　Ｔｅｍｐｅｒａｔｕｒｅ：ＳＫＴ）、及び、指先の皮膚コンダクタンス（Ｓｋｉｎ　Ｃｏｎｄｕｃｔａｎｃｅ：ＳＣ）である。これらの生体信号は同時に測定された。そして、各生体信号から複数種類の生体指標を得た。以下、ＥＣＧを用いて検討した結果について説明する。
　測定されたＥＣＧから、連続する２つの心拍のＲ波のピークの間隔である心拍間隔（Ｒ－Ｒ　ｉｎｔｅｒｖａｌｓ：ＲＲＩ）が算出された（図７の（ａ）参照）。ＲＲＩは、心拍数の指標の１つである。さらに、算出されたＲＲＩから、心拍変動の変動係数（Ｃｏｅｆｆｉｃｉｅｎｔ　ｏｆ　Ｖａｒｉａｔｉｏｎ　ｏｆ　Ｒ－Ｒ　ｉｎｔｅｒｖａｌｓ：ＣｖＲＲ）が算出された。ＣｖＲＲは、心拍揺らぎの指標の１つである。ＣｖＲＲは、下記式（１）に示すように、ＲＲＩから、任意時間帯におけるＲＲＩの標準偏差ＳＤを任意時間帯におけるＲＲＩの平均値で規格化することにより算出された。
　ＣｖＲＲ＝任意時間帯における心拍間隔のＳＤ／任意時間帯における心拍間隔の平均　・・・式（１）
　また、連続する各ＲＲＩを、時間とＲＲＩとの２軸の関係に変換し、さらに、ＲＲＩの等間隔時系列データ（図７の（ｂ）参照）に変換した後に、高速フーリエ変換（Ｆａｓｔ　Ｆｏｕｒｉｅｒ　Ｔｒａｎｓｆｏｒｍ：ＦＦＴ）を用いて周波数解析した（図７の（ｃ）参照）。これにより、心拍変動の周波数成分を示す生体指標であるＨＦ（Ｈｉｇｈ　Ｆｒｅｑｕｅｎｃｙ）とＬＦ（Ｌｏｗ　Ｆｒｅｑｕｅｎｃｙ）とが算出された。ＨＦ及びＬＦは、心拍揺らぎの指標である。ＨＦは、０．１４Ｈｚ～０．４Ｈｚの高周波数領域のパワースペクトルの積分値であり、副交感神経の活動量を反映していると考えられている。また、ＬＦは、０．０４Ｈｚ～０．１４Ｈｚの低周波数領域のパワースペクトルの積分値であり、交感神経及び副交感神経の活動量を反映すると考えられている。なお、ＦＦＴを用いた周波数解析を行うデータは、６０秒間の心拍変動のデータであり、周波数変換は、５秒間隔で行われた。
　被測定者の安静時の生体指標、及び、被測定者がタスクを実行している間に測定された生体指標は、それぞれ測定開始６０秒後から２４０秒間の生体指標の平均値である。また、生体指標の変化量は、被測定者の安静時の生体指標の平均値である基準からの被測定者がタスクを実行している間に測定された生体指標の平均値への変化量である。なお、変化量は、比又は差で表される。生体指標の変化量が比で表される場合、生体指標の変化量は、下記式（２）を用いて算出される。
　生体指標の変化量＝（タスク実行中の生体指標の平均値－安静時の生体指標の平均値）／安静時の生体指標の平均値　・・・式（２）
　続いて、ストレスの要因を判定する性能の高い生体指標の変化量の組み合わせを検討した。具体的には、算出されたＲＲＩ、ＣｖＲＲ、ＬＦ、及び、ＨＦのそれぞれの変化量を用いて線形判別分析を行った。
　ＲＲＩ及びＣｖＲＲの変化量を用いて線形判別分析を行った結果、判定精度は７５．０％であった。したがって、ＲＲＩの変化量及びＣｖＲＲの変化量を用いると、比較的高い精度でストレスの要因を判定できることが分かった。
　また、ＲＲＩ、ＬＦ及びＨＦの変化量を用いて線形判別分析を行った結果、判定精度は６７．５％であった。したがって、ＲＲＩの変化量、ＬＦの変化量及びＨＦの変化量を用いると、比較的良好な精度でストレスの要因を判定できることが分かった。
　一方、ＬＦ及びＨＦの変化量を用いて線形判別分析を行った結果、判定精度は４６．３％であった。したがって、ＬＦの変化量及びＨＦの変化量を用いると、ＲＲＩの変化量を含む組み合わせに比べて、判定精度が大きく低下した。以上の検討により、ＲＲＩの変化量及びＣｖＲＲの変化量を用いると、比較的高い精度でストレスの要因を判定できることが分かった。
　そこで、生体指標の変化量としてＲＲＩの変化量及びＣｖＲＲの変化量を用いてストレスの要因を判定した。図１は、被験者２０名それぞれのストレスの要因毎の生体指標の変化量をプロットした図である。思考疲労に関するストレス１及び２は、ともに同様の結果を示したため、思考疲労に関するストレスとして図示した。図１から、生体指標の変化量は、実行されるタスクの種類によって変化の傾向が異なることが分かった。変化の傾向をより明確にするために、被験者２０名の生体指標の変化量の平均値を求めた。図２は、被験者２０名のストレス要因毎の生体指標の変化量の平均値を示す図である。図２から、ストレスの要因によって生体指標の変化量は、以下の特徴的な変化の傾向を有することが分かった。
　ストレスの要因が対人に関する要因である場合、ＲＲＩの変化量はマイナス側に大きく移行し（すなわち、心拍数が大きくなり）、ＣｖＲＲの変化量はプラス側に移行する傾向がある。また、ストレスの要因が痛みである場合、ＲＲＩの変化量はプラス側に移行し（すなわち、心拍数が小さくなり）、ＣｖＲＲの変化量はマイナス側にわずかに移行する傾向がある。また、ストレスの要因が思考疲労である場合、ＲＲＩの変化量はマイナス側にごくわずかに移行し（すなわち、心拍数はあまり変化せず）、ＣｖＲＲの変化量はマイナス側に大きく移行する傾向があることが分かった。
　以上の結果により、ＲＲＩの変化量とＣｖＲＲの変化量とを用いてストレスの要因を判定すると、比較的高い判定精度が得られることが分かった。また、ＲＲＩの変化量及びＣｖＲＲの変化量は、ストレスの要因によって変化の傾向があることが分かった。これらの変化量の変化の傾向に基づいて、被験者のストレスの要因を容易に、かつ、精度良く判定できることが分かった。
　以上の検討結果から、本発明者らは、ストレスの要因によって各生体指標の変化量は所定の変化の傾向を有しており、特に、心拍数及び心拍揺らぎに関する生体指標の変化量の両方を判定の指標に用いることにより、いずれか一方を判定の指標に用いる場合よりも、より正確にストレスの要因を判定できるとの知見を得た。そして、この検討結果を基に、被測定者から得られる複数種類の生体指標の変化量と閾値とを比較することにより、被測定者のストレスの要因及びストレスの強度を判定する装置の発明に想到した。
　そこで、本開示は、被測定者のストレスの要因を判定できるストレス評価装置、ストレス評価方法及びプログラムを提供する。
　本開示の一態様の概要は、以下の通りである。
　本開示の一態様に係るストレス評価装置は、被測定者の心拍数及び心拍揺らぎを測定する第１センサ部と、（ｉ）心拍数の変化量、及び、（ｉｉ）心拍揺らぎの変化量を算出する演算部と、（ｉ）前記心拍数の変化量及び（ｉｉ）前記心拍揺らぎの変化量に基づいて前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定部と、を備え、前記心拍数の変化量は、前記被測定者の安静時の心拍数である基準からの前記第１センサ部によって測定された前記心拍数への変化量であり、前記心拍揺らぎの変化量は、前記被測定者の安静時の心拍揺らぎである基準からの前記第１センサ部によって測定された前記心拍揺らぎへの変化量であり、前記判定部は、（Ｉ）前記心拍数の変化量と第１の閾値との大小関係の比較、及び、（ＩＩ）前記心拍揺らぎの変化量と第２の閾値との大小関係の比較を行うことにより、前記ストレスの要因を判定する。
　上記構成によれば、被測定者の安静時の各生体指標を基準として各生体指標の変化量を算出するため、各生体指標の推移をより正確に把握することができる。したがって、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因を判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記心拍数の変化量は、第１の時間に測定された前記心拍数の変化量であり、前記心拍揺らぎの変化量は、第２の時間に測定された前記心拍揺らぎの変化量であり、前記第１の閾値は、前記被測定者の安静時の心拍数を基準とする、前記第１及び前記第２の時間とは異なる任意の時間に測定された前記心拍数であり、前記第２の閾値は、前記被測定者の安静時の心拍揺らぎを基準とする、前記任意の時間に測定された前記心拍揺らぎであってもよい。
　ここで、任意の時間とは、例えば、被測定者がストレスを感じる手前の状態にある時を指す。これにより、第１の閾値及び第２の閾値を正確に設定することができる。例えば、各生体指標の変化量と閾値との大小関係を比較する場合、被測定者の睡眠中又は就寝直前等の所定の時刻に測定した各生体指標を各生体指標の閾値に設定してもよい。これにより、被測定者が任意の時間を都度設定することなく女性の月経変動、又は、経年変動等を考慮した閾値を設定できるため、より正確にストレスの要因を判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記心拍揺らぎは、前記被測定者の心拍間隔を周波数分析して求められてもよい。
　これにより、ストレス評価装置は、心拍揺らぎの周波数成分から呼吸間隔及び血圧の情報を得ることができる。したがって、ストレス評価装置は、被測定者の詳細な情報を含む生体指標をストレスの判定のための指標（判定指標）に用いることができるため、被測定者のストレスの要因をより正確に判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記判定部は、前記心拍数の変化量が前記第１の閾値よりも大きく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも大きい場合、前記ストレスの要因は、対人に関する要因であると判定してもよい。
　上記構成によれば、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因が対人に関する要因であると判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記判定部は、前記心拍数の変化量が前記第１の閾値よりも大きく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも小さい場合、前記ストレスの要因は、痛みであると判定してもよい。
　上記構成によれば、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因が痛みであると判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記判定部は、前記心拍数の変化量が前記第１の閾値よりも小さく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも大きい場合、思考による疲労であると判定してもよい。
　上記構成によれば、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因が思考による疲労であると判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、さらに、前記判定部は、前記心拍数の変化量と前記第１の閾値との差、及び、前記心拍揺らぎの変化量と前記第２の閾値との差に応じて、前記ストレスの強度を判定し、判定結果を前記判定結果に基づく前記情報として出力してもよい。
　これにより、被測定者は、自身のストレスの強度を知ることができる。これにより、ストレスの制御について意識を持ちやすくなり、自身のストレスに対する傾向を把握しやすくなる。例えば、被測定者は、複数種類のストレスの要因の中でも耐えうるストレスの強度が異なることを認識することができる。これにより、被測定者は、ストレスの状況に応じてストレスの制御がすぐに必要かどうかを判断することができるようになる。そのため、被測定者は、ストレスの制御を効率良く行うことができるため、ストレスの制御を継続して行うことができる。
　例えば、本開示の一態様に係るストレス評価装置は、さらに、前記判定部によって出力された前記判定結果に基づく前記情報を提示する提示部を備え、前記情報は、前記ストレスの要因、前記ストレスの強度及び前記ストレスの低減策からなる群から選択される少なくとも１つを含んでもよい。
　これにより、被測定者は、ストレスを受けた直後に、自身のストレスの状況及びストレスの制御方法を知ることができるため、ストレスの蓄積をより低減することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記提示部は、音声で提示してもよい。
　これにより、被測定者は、日常生活を送りながら簡便に自身のストレスの状況及び制御方法を知ることができるため、自身のストレスの制御に対する意識を維持しやすくなる。そのため、被測定者は、自身のストレスの制御を継続して行うことができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記提示部は、画像で提示してもよい。
　これにより、被測定者は、視覚的に自身のストレスの状況及び制御方法を知ることができるため、自身のストレスの制御について明確に意識することができる。そのため、被測定者は、自身のストレスの制御を継続して行うことができる。
　また、本開示の一態様に係るストレス評価方法は、測定された被測定者の心拍数及び心拍揺らぎを取得する取得ステップと、（ｉ）心拍数の変化量、及び、（ｉｉ）心拍揺らぎの変化量を算出する算出ステップと、前記心拍数の変化量及び前記心拍揺らぎの変化量に基づいて前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定ステップと、を含み、前記心拍数の変化量は、前記被測定者の安静時の心拍数である基準からの前記第１センサ部によって測定された前記心拍数への変化量であり、前記心拍揺らぎの変化量は、前記被測定者の安静時の心拍揺らぎである基準からの前記第１センサ部によって測定された前記心拍揺らぎへの変化量であり、前記判定ステップでは、（Ｉ）前記心拍数の変化量と第１の閾値との大小関係を比較し、かつ、（ＩＩ）前記心拍揺らぎの変化量と第２の閾値との大小関係を比較することにより、前記ストレスの要因を判定する。
　上記方法によれば、被測定者の安静時の各生体指標を基準として各生体指標の変化量を算出するため、各生体指標の推移をより正確に把握することができる。そのため、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因を判定することができる。
　なお、これらの包括的又は具体的な態様は、システム、方法、集積回路、コンピュータプログラム又はコンピュータで読み取り可能なＣＤ－ＲＯＭなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラム及び記録媒体の任意な組み合わせで実現されてもよい。
　以下、本開示の実施の形態について、図面を参照しながら具体的に説明する。
　なお、以下で説明する実施の形態は、いずれも包括的又は具体的な例を示すものである。以下の実施の形態で示される数値、形状、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、本開示を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、最上位概念を示す独立請求項に記載されていない構成要素については、任意の構成要素として説明される。また、各図は、必ずしも厳密に図示したものではない。各図において、実質的に同一の構成については同一の符号を付し、重複する説明は省略又は簡略化することがある。
　（実施の形態１）
　以下、本実施の形態に係るストレス評価装置、ストレス評価方法及びプログラムについて具体例を挙げて説明する。
　［ストレス評価装置の概要］
　図３は、本実施の形態に係るストレス評価装置１００の概略構成図である。図３に示すように、ストレス評価装置１００は、第１センサ部１１ａと、演算部１２と、判定部１３と、提示部１４と、記憶部１５と、を備える。ストレス評価装置１００では、例えば、第１センサ部１１ａは、被測定者の生体信号を測定するウエアラブルの第１生体センサ１１１ａ（図４参照）を含む。第１センサ部１１ａは、第１生体センサ１１１ａで測定された生体信号から複数種類の生体指標を算出し、測定された生体指標として演算部１２に出力する。演算部１２は、被測定者の安静時の各生体指標の平均値（以下、基準値ともいう）及び各生体指標の閾値を算出し、記憶部１５に格納させる。また、演算部１２は、測定された各生体指標の平均値及び各生体指標の変化量を算出し、判定部１３に出力する。判定部１３は、各生体指標の変化量に基づいて被測定者のストレスの要因を判定する。より具体的には、判定部１３は、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因を判定する。また、判定部１３は、各生体指標の変化量と各生体指標の閾値との差に応じて、ストレスの強度を判定する。そして、判定部１３は、これらの判定結果に基づく情報を提示部１４に出力する。このとき、判定部１３は、判定結果に基づく情報を記憶部１５に格納させる。提示部１４は、判定結果に基づく情報を提示する。さらに、ストレス評価装置１００は、被測定者（ユーザ）の指示を入力する入力部１６（図４参照）を備えてもよい。判定部１３は、入力部１６に入力された被測定者の指示に基づいて判定結果の情報を提示部１４に提示させる。
　［ストレス評価装置の構成］
　本実施の形態に係るストレス評価装置１００の構成についてより具体的に説明する。図４は、図３の構成に基づくストレス評価装置の具体例を示す構成図である。
　図４に示すように、ストレス評価装置１００は、第１生体センサ１１１ａと第１信号処理部１１２ａとを含む第１センサ部１１ａと、演算部１２と、判定部１３と、提示部１４と、記憶部１５と、入力部１６と、を備える。
　第１生体センサ１１１ａは、被測定者の生体信号を測定する。生体信号は、生体情報の信号である。生体情報は、例えば、心拍、脈拍、呼吸数、血中酸素飽和度、血圧、又は、体温などのストレスにより影響を受ける生理学的な情報である。測定の容易性から、生体情報は、例えば、心拍情報である。心拍情報とは、心拍から得られる情報である。また、生体情報は、脈波情報であってもよい。
　第１生体センサ１１１ａは、心拍情報又は脈波情報を取得するセンサである。第１生体センサ１１１ａが心拍情報を取得するセンサ（以下、心拍センサ）である場合、心拍センサは、例えば、被測定者の体の表面に接触する一対の検出電極を備えるセンサである。心拍センサにより得られる心拍情報は、心臓の拍動により得られる電気信号であり、例えば、心電図である。心拍センサは、導電性粘着ゲル電極であってもよく、導電性繊維などで構成されるドライ電極であってもよい。心拍センサの装着部位は、胸部であり、心拍センサの形状は、例えば、ウエアと電極とが一体となったウエア型である。
　第１生体センサ１１１ａが脈波情報を取得するセンサ（以下、脈波センサ）である場合、脈波センサは、例えば、フォトトランジスタ及びフォトダイオードにより血管中の血液量の変化を反射光又は透過光により測定するセンサである。脈波センサは、ユーザの手首に装着され、当該装着された形状で脈波情報を測定する。脈波センサの装着部位は、足首、指、上腕などでもよい。脈波センサの形状は、バンド型（例えば、腕時計型）に限定されず、頸部等に貼り付ける貼付型、メガネ型などであってもよい。また、脈波センサは、顔又は手などの皮膚の色度の変化から脈波情報を測定して脈拍を算出する画像センサであってもよい。
　第１生体センサ１１１ａで測定された生体信号は、第１信号処理部１１２ａに出力される。
　第１信号処理部１１２ａは、第１生体センサ１１１ａで測定された１つの生体信号から複数種類の生体指標を算出する。本実施の形態では、生体指標１及び生体指標２の２種類の生体指標が算出される。上述したように、生体信号が心電図の場合、複数種類の生体指標は、ＲＲＩ、ＣｖＲＲ、ＨＦ及びＬＦなどである。ＲＲＩは、心拍数の指標であり、ＣｖＲＲ、ＨＦ及びＬＦは、心拍揺らぎの指標である。さらに、第１信号処理部１１２ａは、心拍揺らぎの周波数成分から呼吸数及び血圧の変動の生体指標を算出してもよい。また、これらの複数種類の生体指標のうち判定精度が比較的高い組み合わせは、ＲＲＩ及びＣｖＲＲである。したがって、本実施の形態では、生体指標１及び生体指標２は、それぞれＲＲＩ及びＣｖＲＲである例について説明する。なお、ＲＲＩ及びＣｖＲＲの算出方法については、モニター試験にて上述した通りである。第１信号処理部１１２ａは、算出された生体指標１及び生体指標２を演算部１２に出力する。
　演算部１２は、第１信号処理部１１２ａが出力した生体指標１及び生体指標２を取得し、取得した生体指標１及び生体指標２から生体指標１の変化量及び生体指標２の変化量を算出する。生体指標の変化量は、被測定者の安静時に測定された生体指標（以下、基準値と称する場合がある。）を基準とする測定された生体指標であり、差又は比で表される。各生体指標の基準値は、記憶部１５に格納されている。演算部１２は、記憶部１５に格納された生体指標１及び生体指標２の基準値を読み出し、当該基準値に対する生体指標１及び生体指標２の変化量を算出する。演算部１２は、算出した各生体指標の変化量を判定部１３に出力する。なお、基準値は、季節又は被測定者の生理周期などにより変動する場合があるため、所定の期間毎に更新されてもよい。
　また、演算部１２は、各生体指標の閾値を算出する。生体指標１が、例えば、心拍数である場合、心拍数の変化量は、第１の時間に測定された心拍数の変化量である。第１の閾値は、生体指標１の閾値であり、例えば、心拍数の指標であるＲＲＩの閾値である。第１の閾値は、被測定者の安静時の心拍数を基準とする、任意の時間に測定された心拍数である。また、生体指標２が、例えば、心拍揺らぎである場合、心拍揺らぎの変化量は、第２の時間に測定された心拍揺らぎの変化量である。第２の閾値は、生体指標２の閾値であり、例えば、心拍揺らぎの指標であるＣｖＲＲの閾値である。第２の閾値は、被測定者の安静時の心拍数を基準とする、任意の時間に測定された心拍揺らぎである。つまり、これらの閾値は、第１の時間及び第２の時間とは異なる任意の時間に測定された生体指標の測定値と基準値との差又は比である生体指標の変化量である。ここで、任意の時間とは、例えば、被測定者がストレスを感じる手前の状態にある時を指す。
　以下、本実施の形態では、第１の時間及び第２の時間は、同じ時間である場合について説明するが、第１の時間及び第２の時間は異なる時間であってもよい。例えば、第１信号処理部１１２ａは、第１生体センサ１１１ａで測定された１つの生体信号から時分割で複数種類の心拍数及び心拍揺らぎを算出してもよい。このとき、演算部１２は、第１の時間に測定された心拍数の変化量を算出し、第１の時間とは異なる第２の時間に測定された心拍揺らぎの変化量を算出する。
　演算部１２は、記憶部１５に格納された各生体指標の閾値を読み出し、各生体指標の変化量と閾値との大小関係を比較する。そして、演算部１２は、各生体指標の変化量の少なくとも１つが閾値を一定時間超えている期間をストレス発生期間と判定する。ストレス発生期間とは、被測定者がストレスを感じた期間である。演算部１２は、ストレス発生期間中の各生体指標の変化量から各生体指標の変化量の代表値を算出する。例えば、ストレス発生期間における各生体指標の変化量の代表値は、ストレス発生期間中の各生体指標の変化量の平均値を用いてもよく、基準値からの差分が最も大きい値（最大値）を用いてもよい。
　判定部１３は、演算部１２が出力した生体指標１及び生体指標２の変化量の代表値を取得し、記憶部１５に格納された第１の閾値及び第２の閾値を読み出す。判定部１３は、ストレス発生期間における生体指標１の変化量の代表値と第１の閾値との大小関係を比較し、かつ、生体指標２の変化量の代表値と第２の閾値との大小関係を比較することにより、被測定者のストレスの要因を判定する。つまり、判定部１３は、ストレス発生期間毎にストレスの要因を判定する。生体指標の変化量の代表値は、生体指標の変化量の一例であると言えるため、以下、生体指標の変化量の代表値を単に生体指標の変化量とも呼ぶ。
　具体的には、判定部１３は、生体指標１（ここでは、心拍数）の変化量が第１の閾値よりも大きく、かつ、生体指標２（ここでは、心拍揺らぎ）の変化量が第２の閾値よりも大きい場合、ストレスの要因は、対人に関する要因であると判定する。また、判定部１３は、生体指標１の変化量が第１の閾値よりも大きく、かつ、生体指標２の変化量が第２の閾値よりも小さい場合、ストレスの要因は、痛みであると判定する。また、判定部１３は、生体指標１の変化量が第１の閾値よりも小さく、かつ、生体指標２の変化量が第２の閾値よりも大きい場合、ストレスの要因は、思考による疲労であると判定する。
　さらに、判定部１３は、生体指標１の変化量と第１の閾値との差、及び、生体指標２の変化量と第２の閾値との差に応じて、ストレスの強度を判定し、判定結果を当該判定結果に基づく情報として出力する。判定結果に基づく情報は、例えば、ストレスの要因、ストレスの強度及びストレスの低減策の少なくとも１つを含む。ストレスの低減策は、例えば、ストレスの解消方法又はストレスの回避方法などである。ストレスの低減策は、後述する提示情報テーブルに含まれる。判定部１３は、記憶部１５に格納された提示情報テーブルから適切なストレス低減策を読み出し、提示部１４に出力する。
　また、判定部１３は、判定結果に基づく情報を記憶部１５に格納する。このとき、判定部１３は、被測定者がストレスを感じた時間の情報と上記判定結果に基づく情報とを紐づけして記憶部１５に格納してもよい。
　提示部１４は、判定部１３によって出力された上記判定結果に基づく情報を提示する。提示部１４は、上記判定結果に基づく情報を音声で提示してもよく、画像で提示してもよい。提示部１４が上記情報を音声で提示する場合は、提示部１４は、例えば、スピーカーである。また、提示部１４が上記情報を画像で提示する場合は、提示部１４は、例えば、ディスプレイである。
　記憶部１５は、各生体指標の基準値、各生体指標の閾値、及び、提示情報テーブルなどを格納する。提示情報テーブルは、ストレスの要因及び当該ストレスの強度に応じて提示されるストレス低減策などの提示情報のテーブルである。上述したように、各生体指標の基準値及び閾値は、所定の期間で更新されてもよい。なお、提示情報テーブルも同様に、所定の期間で更新されてもよい。
　また、記憶部１５は、判定部１３が出力したストレスの要因、ストレスの強度及びストレス低減策などの判定結果に基づく情報を格納する。このとき、記憶部１５は、判定結果に基づく情報とストレス発生期間とを紐付けて格納してもよい。これにより、被測定者は、所望のタイミングで判定結果に基づく情報を呼び出すことができる。このとき、判定部１３は、入力部１６により入力された被測定者の操作に基づいて、判定結果に基づく情報を提示部１４に提示させる。
　入力部１６は、被測定者による操作を示す操作信号を判定部１３に出力する。入力部１６は、例えば、キーボード、マウス、タッチパネル、又は、マイクなどである。操作信号とは、判定結果に基づく情報の抽出方法又は提示部１４における提示方法などの設定を行う信号である。提示部１４には、入力部１６に入力された設定に基づき、様々な形式の判定結果が提示される。例えば、所定の期間におけるストレスの変化、被測定者が影響を受けやすいストレスの要因、及び、被測定者に適したストレス低減策などである。これにより、被測定者は、短期的なストレスの傾向を把握できるだけでなく、中期的及び長期的なストレスの傾向を把握することができる。このように、被測定者は、自己に適した効果的なストレス低減策を知ることができるため、中長期的なストレスを制御することができる。
　［ストレス評価方法］
　次に、本実施の形態に係るストレス評価方法について図５を用いて具体的に説明する。図５は、実施の形態に係るストレス評価方法を説明するフローチャートである。
　本実施の形態に係るストレス評価方法は、測定された被測定者の心拍数及び心拍揺らぎを取得する取得ステップＳ１０と、（ｉ）心拍数の変化量、及び、（ｉｉ）心拍揺らぎの変化量を算出する算出ステップＳ２０と、心拍数の変化量及び心拍揺らぎの変化量に基づいて被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定ステップＳ３０と、を含む。心拍数の変化量は、被測定者の安静時の心拍数である基準からの第１センサ部１１ａによって測定された心拍数への変化量であり、心拍揺らぎの変化量は、被測定者の安静時の心拍揺らぎである基準からの第１センサ部１１ａによって測定された心拍揺らぎへの変化量である。判定ステップＳ３０では、（Ｉ）前記心拍数の変化量と第１の閾値との大小関係を比較し、かつ、（ＩＩ）前記心拍揺らぎの変化量と第２の閾値との大小関係を比較することにより、前記ストレスの要因を判定する。本実施の形態では、さらに、判定ステップＳ３０の判定結果に基づく情報を提示する提示ステップＳ４０を含む。
　以下、各ステップについてより具体的に説明する。
　まず取得ステップＳ１０では、演算部１２は、第１センサ部１１ａで測定された被測定者の複数種類の生体指標（ここでは、心拍数及び心拍揺らぎ）を取得する。第１センサ部１１ａでは、第１生体センサ１１１ａで心拍情報（ここでは、心電図）が測定され、第１信号処理部１１２ａで、心拍数の指標及び心拍揺らぎの指標などの生体指標が算出される。なお、上述した通り、生体情報は、心拍情報に限られず、脈波情報などのストレスで影響を受ける生理学的な情報であってもよい。特に、心拍情報は、ウエアラブルな生体センサを用いた場合、脈拍、呼吸数、血圧、及び血中酸素飽和度などの他の生体情報よりも被測定者の負担が少ない状態で簡便に、かつ、リアルタイムに測定することができる。そのため、生体情報として被測定者の心拍情報を用いることにより、被測定者のストレスの状態を適切に評価することができる。
　心拍情報から得られる生体指標は、心拍数の指標であるＲＲＩ、心拍揺らぎの指標であるＣｖＲＲ、ＬＦ、ＨＦ、及び、ＬＦ／ＨＦなどである。このように１つの生体情報から複数種類の生体指標が得られる。また、上述したように、これらの生体指標の組み合わせにより、比較的高い判定精度でストレスの要因を判定することができるため、信頼性の高い評価が得られる。
　図６は、本実施の形態に係るストレス評価装置１００の第１センサ部１１ａで得られる心拍情報の一例を示す図である。心拍情報は、例えば、心電図であり、図６に示す心電波形となる。心電波形は、心房の電気的興奮を反映するＰ波と、心室の電気的興奮を反映するＱ波、Ｒ波、及びＳ波と、興奮した心室の心筋細胞が再分極する過程を反映するＴ波とから構成されている。これらの心電波形のうち、Ｒ波の波高（電位差）が最も大きく、筋電位などのノイズに対して最も頑健である。そのため、これらの心電波形における連続する２つの心拍のＲ波のピークの間隔、つまり、心拍間隔（ＲＲＩ）を算出する。心拍数は、ＲＲＩの逆数に６０を乗じて算出される。
　さらに、モニター試験で上述した通り、ＣｖＲＲは、上記式（２）を用いて、ＲＲＩから、任意時間帯におけるＲＲＩの標準偏差ＳＤを心拍間隔の平均値で規格化することにより算出される。
　第１信号処理部１１２ａは、第１生体センサ１１１ａで得られた心拍情報から、左心室が急激に収縮して心臓から血液を送り出す際に発生する電気信号（Ｒ波）を検出し、ＲＲＩを算出する。なお、Ｒ波の検出には、例えば、Ｐａｎ＆Ｔｏｍｐｋｉｎｓ法などの公知の手法が用いられる。
　次に、演算部１２において検出されたＲ波から心拍間隔（ＲＲＩ）の変動量を算出する方法について説明する。
　図７は、心拍間隔（ＲＲＩ）の変動量を算出する手法を説明する図である。第１信号処理部１１２ａは、得られたＲ波の検出データから、以下のように、ＲＲＩの変動量を算出する。
　図７の（ａ）に示すように、第１信号処理部１１２ａは、連続する２つの心拍のＲ波のピークの間隔であるＲＲＩを算出する。第１信号処理部１１２ａは、算出された各ＲＲＩを時間とＲＲＩとの２軸の関係に変換する。変換されたデータは不等間隔の離散的なデータであるため、演算部１２は、変換されたＲＲＩの時系列データを、図７の（ｂ）に示す等間隔時系列データに変換する。次いで、演算部１２は、この等間隔時系列データに対して、例えば、高速フーリエ変換（ＦＦＴ）を用いて周波数解析することにより、図７の（ｃ）に示す心拍変動の周波数成分を求める。
　心拍変動の周波数成分は、例えば、高周波成分ＨＦと低周波成分ＬＦとに分けることができる。モニター試験で上述した通り、ＨＦは、副交感神経活動量を反映していると考えられる。また、ＬＦは、交感神経及び副交感神経の活動量を反映すると考えられている。そのため、ＬＦとＨＦとの比であるＬＦ／ＨＦは、交感神経活動量を示すと考えられる。
　このように、第１センサ部１１ａでは、心拍情報から複数種類の生体指標が算出される。
　取得ステップＳ１０では、演算部１２にて、これらの生体指標から２種類の生体指標（ここでは、心拍数及び心拍揺らぎ）を取得する。
　次いで、算出ステップＳ２０では、演算部１２にて、取得ステップＳ１０で取得された２種類の生体指標の変化量を算出する。各生体指標の変化量は、上述した通り、被測定者の安静時の各生体指標の値を基準値として、各生体指標の基準値と取得された各生体指標の値との比又は差を算出して得られる。演算部１２は、記憶部１５に格納された各生体指標の基準値を読み出して使用する。
　なお、各生体指標の変化量は、例えば、変化量が差で表される場合は、取得ステップＳ１０で取得された各生体指標の値から各生体指標の基準値を差し引くことにより算出される。例えば、心拍数の変化量は、取得ステップＳ１０で取得された被測定者の心拍数の値から心拍数の基準値を差し引くことにより算出される。また、変化量が比で表される場合は、取得ステップＳ１０で取得された各生体指標の値を各生体指標の基準値で割ることにより算出される。例えば、心拍数の変化量は、取得ステップＳ１０で取得された被測定者の心拍数の値を心拍数の基準値で割ることにより算出される。
　以上のように、算出ステップＳ２０では、演算部１２にて、各生体指標の変化量を算出する。
　次いで、判定ステップＳ３０では、判定部１３にて、算出ステップＳ２０で算出された各生体指標の変化量に基づいてストレスの要因を判定する。判定部１３は、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、被測定者のストレスの要因を判定する。具体的には、判定ステップＳ３０では、判定部１３は、心拍数の変化量が第１の閾値よりも大きく、かつ、心拍揺らぎの変化量が第２の閾値よりも大きい場合、ストレスの要因は、対人に関する要因であると判定する。また、判定部１３は、生体指標１の変化量が第１の閾値よりも大きく、かつ、生体指標２の変化量が第２の閾値よりも小さい場合、ストレスの要因は、痛みであると判定する。また、判定部１３は、生体指標１の変化量が第１の閾値よりも小さく、かつ、生体指標２の変化量が第２の閾値よりも大きい場合、ストレスの要因は、思考による疲労であると判定する。
　さらに、判定部１３は、生体指標１の変化量と第１の閾値との差、及び、生体指標２の変化量と第２の閾値との差に応じて、ストレスの強度を判定し、判定結果を当該判定結果に基づく情報として出力する。
　なお、第１の閾値は、心拍数の閾値であり、被測定者の安静時の心拍数を基準とする、第１の時間及び第２の時間とは異なる任意の時間に測定された心拍数である。第２の閾値は、心拍揺らぎの閾値であり、被測定者の安静時の心拍揺らぎを基準とする、第１の時間及び第２の時間とは異なる任意の時間に測定された心拍揺らぎである。これらの閾値は、演算部１２にて算出され、記憶部１５に格納される。判定部１３は、記憶部１５に格納された各生体指標の閾値を読み出して使用する。上述したように、任意の時間とは、例えば、被測定者がストレスを感じる手前の状態にある時を指す。
　各生体指標の閾値は、各生体指標の変化量が正の値である場合の閾値と、各生体指標の変化量が負の値である場合の閾値とが設定される。基準値は変化量のゼロ点である。各生体指標の変化量と閾値との大小関係は、以下のように比較される。生体指標の変化量が正の値である場合、生体指標の変化量と正の閾値との大小関係を比較する。また、生体指標の変化量が負の値である場合、生体指標の変化量の絶対値と負の閾値の絶対値との大小関係を比較する。なお、各生体指標の閾値は、固定値であってもよく、所定の期間で更新されてもよく、日々の測定に基づいて都度更新されてもよい。
　なお、閾値は、線形判別又は決定木等の比較的単純な機械学習によって算出されてもよい。これにより、被測定者に適した判定基準値及び閾値を設定できるため、ストレスの要因をより精度良く判定することができる。
　以上のように、判定ステップＳ３０では、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、被測定者のストレスの要因を判定する。
　次いで、提示ステップＳ４０では、提示部１４にて、判定部１３で判定された判定結果に基づく情報を提示する。提示部１４は、判定結果に基づく情報を音声で提示してもよく、画像で提示してもよい。判定結果に基づく情報は、ストレスの要因、ストレスの強度及びストレスの低減策の少なくとも１つを含む。提示部１４は、被測定者が入力部１６で入力した設定に基づき、様々な形式の判定結果を表示する。
　［ストレス評価装置の使用例］
　次に、本実施の形態に係るストレス評価装置１００の使用例について具体的に説明する。図８は、本実施の形態に係るストレス評価装置１００の使用例を説明する図である。
　図８に示すように、ストレス評価装置１００は、第１センサ部１１ａの一部である第１生体センサ１１１ａと、第１生体センサ１１１ａ以外の構成を含む評価端末２０とから構成される。被測定者は、第１生体センサ１１１ａを胸部の肌に接触するように装着し、心電図（ＥＣＧ）を測定する。第１生体センサ１１１ａは、導電性粘着ゲル電極であってもよく、導電性繊維などで構成されるドライ電極であってもよい。第１生体センサ１１１ａは、測定した心拍の電気信号を、通信により評価端末２０に送信する。通信方法は、Ｂｌｕｅｔｏｏｈ（登録商標）などの無線通信であってもよく、有線通信であってもよい。
　評価端末２０は、第１センサ部１１ａの第１信号処理部１１２ａ、演算部１２、判定部１３、提示部１４、記憶部１５及び入力部１６を備える。第１信号処理部１１２ａは、第１生体センサ１１１ａから通信により送信された心拍の電気信号を受信する。第１信号処理部１１２ａは、受信した心拍の電気信号から心拍数の指標であるＲＲＩ及び心拍揺らぎの指標であるＣｖＲＲを算出し、これらの生体指標を演算部１２に出力する。
　演算部１２は、第１信号処理部１１２ａが出力したＲＲＩ及びＣｖＲＲを取得し、記憶部１５に格納されたＲＲＩの基準値及びＣｖＲＲの基準値を読み出す。演算部１２は、読み出した基準値を基準とする、これらの生体指標である生体指標の変化量をそれぞれ算出する。生体指標の変化量は、差又は比で表される。本実施の形態では、当該変化量は、比で表される。
　また、上述した通り、演算部１２は、各生体指標の閾値を算出し、記憶部１５に出力する。各生体指標の閾値は、各生体指標の変化量が正の値になる場合の閾値と、各生体指標の変化量が負の値になる場合の閾値とが設定される。基準値は変化量ゼロである。具体的には、各生体指標の変化量が正の値になる場合、正の閾値は、基準値よりも大きい値であり、変化量のグラフ１２０中の第１の閾値１ａ（以下、正の閾値１ａ）及び第２の閾値２ａ（以下、正の閾値２ａ）である。各生体指標の変化量が負の値になる場合、負の閾値は、基準値よりも小さい値であり、変化量のグラフ１２０中の第１の閾値１ｂ（以下、負の閾値１ｂ）及び第２の閾値２ｂ（以下、負の閾値２ｂ）である。また、演算部１２は、各生体指標の基準値を算出し、記憶部１５に出力する。各生体指標の基準値は、各生体指標の変化量がゼロである。例えば、変化量のグラフ１２０では、基準値は、正の閾値１ａ及び負の閾値１ｂの間の実線である。なお、正の閾値及び負の閾値は、基準値（変化量ゼロ）を挟んで等間隔に設定されてもよく、基準値を挟んで等間隔に設定されなくてもよい。これらの閾値は、各生体指標の変化量の大きさに応じて、適宜設定されてもよい。
　判定部１３は、演算部１２が出力した各生体指標の変化量を取得し、記憶部１５に格納された各生体指標の閾値を読み出す。判定部１３は、各生体指標の変化量と各生体指標の閾値との大小関係を比較し、ストレスの要因を判定する。例えば、各生体指標の変化量が正の値である場合、判定部１３は、各生体指標の変化量と正の閾値との大小関係を比較する。また、各生体指標の変化量が負の値である場合、判定部１３は、各生体指標の変化量の絶対値と負の閾値の絶対値との大小関係を比較する。以下、変化量のグラフ１２０及び判定表１３０を用いて、より具体的に説明する。
　変化量のグラフ１２０に示すように、期間Ａ１では、ＲＲＩの変化量の絶対値は、負の閾値１ｂの絶対値よりも大きく、かつ、ＣｖＲＲの変化量は、正の閾値２ａよりも大きい。よって、判定部１３は、被測定者が期間Ａ１で感じたストレスの要因は、対人に関する要因であると判定する。また、期間Ｂ１では、ＲＲＩの変化量は、正の閾値１ａよりも大きく、かつ、ＣｖＲＲの変化量の絶対値は、負の閾値２ｂの絶対値よりも小さい。よって、判定部１３は、被測定者が期間Ｂ１で感じたストレスの要因は、痛みであると判定する。また、期間Ｃ１では、ＲＲＩの変化量の絶対値は、負の閾値１ｂの絶対値よりも小さく、かつ、ＣｖＲＲの変化量の絶対値は、負の閾値２ｂの絶対値よりも大きい。よって、判定部１３は、被測定者が期間Ｃ１で感じたストレスの要因は、思考による疲労（思考疲労）であると判定する。
　判定表１３０では、矢印の向き及び本数で、基準値（変化量ゼロ）に基づく各生体指標の変化量の推移を示している。横向きの矢印は、生体指標の変化量が閾値を超える変化を伴わないことを示している。
　さらに、判定部１３は、ＲＲＩの変化量の絶対値と第１の閾値の絶対値との差、及び、ＣｖＲＲの変化量の絶対値と第２の閾値の絶対値との差に応じて、ストレスの強度を判定する。
　判定部１３は、これらの判定結果に基づく情報を提示部１４に出力する。提示部１４は、例えば、スマートフォンのディスプレイである。また、判定部１３は、被測定者は、所望のタイミングで判定結果に基づく情報を呼び出すことができる。このとき、判定部１３は、タッチパネルなどの入力部１６により入力された被測定者の操作に基づいて、判定結果に基づく情報を提示部１４に提示させる。例えば、被測定者が評価端末２０の入力部１６で必要な情報を抽出する指示を入力すると、判定部１３は、被測定者の指示に基づいて提示部１４に提示情報１４０を提示する。提示情報１４０は、被測定者がストレスを感じた時間、ストレスの要因、及びストレスの低減策を含んでいる。ストレスの低減策は、例えば、ストレスの要因に応じたストレス解消方法又はストレス回避方法を提案するメッセージである。当該メッセージは、例えば、ストレスの要因が思考疲労である場合、少し休憩しましょう、又は、ストレッチをしましょう、などであり、対人に関わる要因である場合、少し瞑想しましょう、又は、深呼吸をしましょう、などである。
　以上のように、本実施の形態によれば、被測定者が日常生活を送りながら簡便に、かつ、正確に、ストレスの要因を判定することができる。そのため、被測定者は、従来よりも正確に自身のストレス状態及び適切なストレス低減策を把握することができる。これにより、被測定者は、適切に、かつ、効率良く、自身のストレスの制御を行うことができるため、ストレスの制御を継続して行うことができる。
　（本開示の基礎となった第２の知見）
　本発明者らは、本開示の基礎となった第１の知見に記載の上記課題に鑑みて鋭意検討をした。検討内容を以下に記す。
　本発明者らは、ストレスの要因と、心拍情報及び発汗情報などの生体情報から得られる生体指標との関連性を見出すために、以下のモニター試験を実施した。
　［モニター試験］
　２０名の被験者に対してストレスの要因が異なる４つのタスクを与え、タスクを実行している被験者の生体信号を測定した。
　被験者は、２０代から３０代の社会人又は大学生の男女であり、健康状態及び精神状態に関するアンケ―トの結果が異常値を示さない２０名が選出された。
　タスクは、［１］対人に関するストレス、［２］痛みに関するストレス、［３］思考による疲労（以下、思考疲労）に関するストレス１、［４］思考疲労に関するストレス２の４種類である。各タスクは各被験者に対して個別に実施された。タスクの詳細は、第１の知見に記載のモニター試験と同様であるため、ここでの記載を省略する。
　上記のモニター試験は、日内変動を考慮し、被験者毎に別日の同時刻に実施した。
　被験者の安静時の生体信号は、上記［１］～［４］の各タスクを実施する前に、タスクを実行する姿勢と同じ姿勢で、５分間測定した生体信号である。この生体信号から生体指標を算出し、生体指標の変化量を算出するための基準値とした。生体指標の変化量は、被験者の安静時の生体指標を基準とするタスク実行中に測定された被験者の生体信号から算出された生体指標である。
　測定された生体信号は、心電図（Ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａｍ：ＥＣＧ）、呼吸間隔、指先温度（Ｓｋｉｎ　Ｔｅｍｐｅｒａｔｕｒｅ：ＳＫＴ）、及び、指先の皮膚コンダクタンス（Ｓｋｉｎ　Ｃｏｎｄｕｃｔａｎｃｅ：ＳＣ）である。これらの生体信号は同時に測定された。そして、各生体信号から複数種類の生体指標を得た。
　生体指標の算出方法は、各生体指標により様々である。例えば、生体指標がＳＫＴである場合、ＳＫＴは、指先の温度を任意の区間で平均化して得られる。なお、ＣｖＲＲ、ＨＦ、ＬＦについても、上述した通りであるため、ここでの記載を省略する。
　続いて、ストレスの要因を判定する性能の高い生体指標の変化量の組み合わせを検討した。具体的には、算出されたＲＲＩ、ＣｖＲＲ、ＬＦ、ＨＦ、ＳＣ、及び、ＳＫＴのそれぞれの変化量を用いて線形判別分析を行った。これら全ての生体指標の変化量を用いて線形判別分析を行った結果、判定精度は約８１．３％であった。また、より単純な決定木による判別では、判定精度は７７．５％であった。
　また、ＲＲＩ、ＣｖＲＲ、及びＳＣの変化量を用いて線形判別分析を行った結果、判定精度は８１．３％であり、決定木による判別では、判定精度は６６．３％であった。したがって、ストレス要因の判定に用いる生体指標の変化量の数を３つに減らしても、比較的高い判定精度を保つことが分かった。
　一方、例えば心拍数の生体指標であるＲＲＩを除き、ＣｖＲＲ及びＳＣの変化量を用いて線形判別分析を行った結果、判定精度は６２．５％であった。したがって、ストレス要因の判定に用いる生体指標の変化量から心拍数の指標であるＲＲＩの変化量を除くと、判定精度が著しく低下することが分かった。
　そこで、生体指標の変化量としてＲＲＩの変化量、ＣｖＲＲの変化量及びＳＣの変化量を用いてストレスの要因を判定した。図９Ａは、被験者２０名それぞれのストレスの要因毎の生体指標の変化量をプロットした図である。図９Ｂは、図９ＡをＲＲＩの変化量を示す軸のプラス側から見た図である。図９Ｃは、図９ＡをＣｖＲＲの変化量を示す軸のマイナス側から見た図である。図９Ｄは、図９ＡをＳＣの変化量を示す軸のマイナス側から見た図である。
　図９Ａ～図９Ｄから、生体指標の変化量は、実行されるタスクの種類によって変化の傾向が異なることが分かった。変化の傾向をより明確にするために、被験者２０名の生体指標の変化量の平均値を求めた。図１０Ａは、図９Ａにプロットした被験者２０名のストレス要因毎の生体指標の変化量の平均値を示す図である。図１０Ｂは、図１０ＡをＲＲＩの変化量を示す軸のプラス側から見た図である。図１０Ｃは、図１０ＡをＣｖＲＲの変化量を示す軸のマイナス側から見た図である。図１０Ｄは、図１０ＡをＳＣの変化量を示す軸のマイナス側から見た図である。図１０Ａ～図１０Ｄから、ストレスの要因によって生体指標の変化量は、以下の特徴的な変化の傾向を有することが分かった。
　ストレスの要因が対人に関する要因である場合、ＲＲＩの変化量はマイナス側に大きく移行し（すなわち心拍数が大きくなり）、ＣｖＲＲの変化量はプラス側に移行し、ＳＣの変化量はプラス側に移行する傾向がある。また、ストレスの要因が痛みである場合、ＲＲＩの変化量はプラス側に移行し（すなわち心拍数が小さくなり）、ＣｖＲＲの変化量はマイナス側にわずかに移行し、ＳＣの変化量はプラス側に大きく移行する傾向がある。また、ストレスの要因が思考疲労である場合、ＲＲＩの変化量はマイナス側にごくわずかに移行し（すなわち心拍数はあまり変化せず）、ＣｖＲＲの変化量はマイナス側に大きく移行し、ＳＣの変化量はプラス側に移行する傾向があることが分かった。
　以上の結果により、ＲＲＩの変化量、ＣｖＲＲの変化量、及び、ＳＣの変化量を用いてストレスの要因を判定すると、比較的高い判定精度が得られることが分かった。また、これらの変化量は、ストレスの要因によって変化の傾向があることが分かった。これらの変化量の変化の傾向に基づいて、被験者のストレスの要因を容易に、かつ、精度良く判定できることが分かった。
　以上の検討結果から、本発明者らは、ストレスの要因によって各生体指標の変化量は所定の変化の傾向を有しており、特に、（ｉ）心拍数、（ｉｉ）心拍揺らぎ、及び、（ｉｉｉ）皮膚コンダクタンス又は皮膚温度に関する生体指標の変化量を判定の指標に用いることにより、比較的高い精度でストレスの要因を判別できるとの知見を得た。そして、この検討結果を基に、被測定者から得られる複数種類の生体指標の変化量と閾値とを比較することにより、被測定者のストレスの要因を判定する装置の発明に結実した。
　そこで、本開示は、被測定者のストレスの要因を判定できるストレス評価装置、ストレス評価方法及びプログラムを提供する。
　本開示の一態様の概要は、以下の通りである。
　本開示の一態様に係るストレス評価装置は、さらに、前記被測定者の皮膚コンダクタンス又は皮膚温度の少なくとも一方を測定する第２センサ部を備え、前記演算部は、さらに、（ｉｉｉ）皮膚コンダクタンスの変化量、又は、皮膚温度の変化量を算出し、前記皮膚コンダクタンスの変化量は、前記被測定者の安静時の皮膚コンダクタンスである基準からの前記第２センサ部によって測定された前記皮膚コンダクタンスへの変化量であり、前記皮膚温度の変化量は、前記被測定者の安静時の皮膚温度である基準からの前記第２センサ部によって測定された前記皮膚温度への変化量であり、前記判定部は、前記（Ｉ）及び前記（ＩＩ）に加えて、（ＩＩＩ）前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量と第３の閾値との大小関係の比較も行うことにより、前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する。
　上記構成によれば、被測定者の安静時の各生体指標を基準として各生体指標の変化量を算出するため、各生体指標の推移をより正確に把握することができる。したがって、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因を判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記心拍数の変化量は、第１の時間に測定された前記心拍数の変化量であり、前記心拍揺らぎの変化量は、第２の時間に測定された前記心拍揺らぎの変化量であり、前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量は、第３の時間に測定された皮膚コンダクタンス又は前記皮膚温度の変化量であり、前記第１の閾値は、前記被測定者の安静時の心拍数を基準とする、前記第１、前記第２及び前記第３の時間とは異なる任意の時間に測定された前記心拍数であり、前記第２の閾値は、前記被測定者の安静時の心拍揺らぎを基準とする、前記任意の時間に測定された前記心拍揺らぎであり、前記第３の閾値は、前記被測定者の安静時の皮膚コンダクタンスを基準とする、前記任意の時間に測定された前記皮膚コンダクタンス、又は、前記被測定者の安静時の皮膚温度を基準とする、前記任意の時間に測定された前記皮膚温度であってもよい。
　ここで、任意の時間とは、例えば、被測定者がストレスを感じる手前の状態にある時を指す。これにより、第１の閾値、第２の閾値及び第３の閾値を正確に設定することができる。
　例えば、各生体指標の変化量と閾値との大小関係を比較する場合、被測定者の睡眠中又は就寝直前等の所定の時刻に測定した各生体指標を各生体指標の閾値に設定してもよい。これにより、被測定者が任意の時間を都度設定することなく、女性の月経変動、又は、経年変動等を考慮した閾値を設定できるため、より正確にストレスの要因を判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記心拍揺らぎは、前記被測定者の心拍間隔を周波数分析して求められてもよい。
　これにより、ストレス評価装置は、心拍揺らぎの周波数成分から呼吸間隔及び血圧の情報を得ることができる。これにより、ストレス評価装置は、被測定者の詳細な情報を含む生体指標を判定指標に用いることができるため、被測定者のストレスの要因をより正確に判定することができる。
　これにより、ストレス評価装置は、心拍揺らぎの周波数成分から呼吸間隔及び血圧の情報を得ることができる。したがって、ストレス評価装置は、被測定者の詳細な状態を含む生体指標をストレスの判定のための指標（判定指標）に用いることができるため、被測定者のストレスの要因をより正確に判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記判定部は、前記心拍数の変化量が第１の閾値よりも大きく、かつ、前記心拍揺らぎの変化量が第２の閾値よりも大きく、かつ、前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量が第３の閾値よりも大きい場合、前記ストレスの要因は、対人に関する要因であると判定してもよい。
　上記構成によれば、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因が対人に関する要因であると判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記判定部は、前記心拍数の変化量が第１の閾値よりも大きく、かつ、前記心拍揺らぎの変化量が第２の閾値よりも小さく、かつ、前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量が第３の閾値よりも大きい場合、ストレスの要因は、痛みであると判定してもよい。
　上記構成によれば、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因が痛みであると判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記判定部は、前記心拍数の変化量が第１の閾値よりも小さく、かつ、前記心拍揺らぎの変化量が第２の閾値よりも大きく、かつ、前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量が第３の閾値よりも小さい場合、ストレスの要因は、思考による疲労であると判定してもよい。
　上記構成によれば、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因が思考による疲労であると判定することができる。
　例えば、本開示の一態様に係るストレス評価装置では、さらに、前記判定部は、前記心拍数の変化量と前記第１の閾値との差、前記心拍揺らぎの変化量と前記第２の閾値との差、及び、前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量と前記第３の閾値との差に応じて、前記ストレスの強度を判定し、判定結果を前記判定結果に基づく前記情報として出力してもよい。
　これにより、被測定者は、自身のストレスの強度を知ることができる。これにより、ストレスの制御について意識を持ちやすくなり、自身のストレスに対する傾向を把握しやすくなる。例えば、被測定者は、複数種類のストレスの要因の中でも耐えうるストレスの強度が異なることを認識することができる。これにより、被測定者は、ストレスの状況に応じてストレスの制御がすぐに必要かどうかを判断することができるようになる。そのため、被測定者は、ストレスの制御を効率良く行うことができるため、ストレスの制御を継続して行うことができる。
　例えば、本開示の一態様に係るストレス評価装置は、さらに、前記判定部によって出力された前記判定結果に基づく前記情報を提示する提示部を備え、前記情報は、前記ストレスの要因、前記ストレスの強度及び前記ストレスの低減策からなる群から選択される少なくとも１つを含んでもよい。
　これにより、被測定者は、ストレスを受けた直後に、自身のストレスの状況及びストレスの制御方法を知ることができるため、ストレスの蓄積をより低減することができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記提示部は、音声で提示してもよい。
　これにより、被測定者は、日常生活を送りながら簡便に自身のストレスの状況及び制御方法を知ることができるため、自身のストレスの制御に対する意識を維持しやすくなる。そのため、被測定者は、自身のストレスの制御を継続して行うことができる。
　例えば、本開示の一態様に係るストレス評価装置では、前記提示部は、画像で提示してもよい。
　これにより、被測定者は、視覚的に自身のストレスの状況及び制御方法を知ることができるため、自身のストレスの制御について明確に意識することができる。そのため、被測定者は、自身のストレスの制御を継続して行うことができる。
　また、本開示の一態様に係るストレス評価方法は、前記取得ステップは、さらに、前記被測定者の皮膚コンダクタンス又は皮膚温度の少なくとも一方を取得し、前記算出ステップは、さらに、（ｉｉｉ）皮膚コンダクタンスの変化量、又は、皮膚温度の変化量を算出し、前記皮膚コンダクタンスの変化量は、前記被測定者の安静時の皮膚コンダクタンスである基準からの前記第２センサ部によって測定された前記皮膚コンダクタンスへの変化量であり、前記皮膚温度の変化量は、前記被測定者の安静時の皮膚温度である基準からの前記第２センサ部によって測定された前記皮膚温度であり、前記判定ステップは、前記（Ｉ）かつ前記（ＩＩ）かつ（ＩＩＩ）前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量と第３の閾値との大小関係を比較することにより、前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する。
　上記方法によれば、被測定者の安静時の各生体指標を基準として各生体指標の変化量を算出するため、各生体指標の推移をより正確に把握することができる。そのため、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因を判定することができる。
　なお、これらの包括的又は具体的な態様は、システム、方法、集積回路、コンピュータプログラム又はコンピュータで読み取り可能なＣＤ－ＲＯＭなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラム及び記録媒体の任意な組み合わせで実現されてもよい。
　以下、本開示の実施の形態２について、図面を参照しながら具体的に説明する。
　（実施の形態２）
　以下、本実施の形態に係るストレス評価装置、ストレス評価方法及びプログラムについて具体例を挙げて説明する。
　［ストレス評価装置の概要］
　図１１は、本実施の形態に係るストレス評価装置１００ａの概略構成図である。図１１に示すように、ストレス評価装置１００ａは、第１センサ部１１ａと、第２センサ部１１ｂと、演算部１２ａと、判定部１３ａと、提示部１４ａと、記憶部１５ａと、を備える。ストレス評価装置１００ａでは、例えば、第１センサ部１１ａ及び第２センサ部１１ｂはそれぞれ、被測定者の生体信号を測定するウエアラブルの第１生体センサ１１１ａ及び第２生体センサ１１１ｂ（図１２参照）を含む。第１センサ部１１ａは、第１生体センサ１１１ａで測定された生体信号から複数種類の生体指標を算出し、測定された生体指標として演算部１２ａに出力する。第２センサ部１１ｂは、第２生体センサ１１１ｂで測定された生体信号から少なくとも１種類の生体指標を算出し、測定された生体指標として演算部１２ａに出力する。演算部１２ａは、被測定者の安静時の各生体指標の平均値（以下、基準値ともいう）及び各生体指標の閾値を算出し、記憶部１５ａに格納させる。また、演算部１２ａは、測定された各生体指標の平均値及び各生体指標の変化量を算出し、判定部１３ａに出力する。判定部１３ａは、各生体指標の変化量に基づいて被測定者のストレスの要因を判定する。より具体的には、判定部１３ａは、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、ストレスの要因を判定する。また、判定部１３ａは、各生体指標の変化量と各生体指標の閾値との差に応じて、ストレスの強度を判定する。そして、判定部１３ａは、これらの判定結果に基づく情報を提示部１４ａに出力する。このとき、判定部１３ａは、判定結果に基づく情報を記憶部１５ａに格納させる。提示部１４ａは、判定結果に基づく情報を提示する。さらに、ストレス評価装置１００ａは、被測定者（ユーザ）の指示を入力する入力部１６ａ（図１２参照）を備えてもよい。判定部１３ａは、入力部１６ａに入力された被測定者の指示に基づいて判定結果の情報を提示部１４ａに提示させる。
　［ストレス評価装置の構成］
　本実施の形態に係るストレス評価装置１００ａの構成についてより具体的に説明する。図１２は、図１１の構成に基づくストレス評価装置の具体例を示す構成図である。
　図１２に示すように、ストレス評価装置１００ａは、第１生体センサ１１１ａと第１信号処理部１１２ａとを含む第１センサ部１１ａと、第２生体センサ１１１ｂと第２信号処理部１１２ｂとを含む第２センサ部１１ｂと、演算部１２ａと、判定部１３ａと、提示部１４ａと、記憶部１５ａと、入力部１６ａと、を備える。
　第１生体センサ１１１ａ及び第２生体センサ１１１ｂは、被測定者の生体信号を測定する。生体信号は、生体情報の信号である。生体情報は、例えば、心拍、脈拍、呼吸数、血中酸素飽和度、血圧、又は、体温などのストレスにより影響を受ける生理学的な情報である。測定の容易性から、生体情報は、例えば、心拍情報である。心拍情報とは、心拍から得られる情報である。また、生体情報は、脈波情報であってもよい。
　第１生体センサ１１１ａ及び第２生体センサ１１１ｂ（以下、単に「生体センサ」と称する）は、各生体情報に応じたセンサを用いる。例えば、生体センサが心拍情報を取得するセンサ（心拍センサ）である場合、心拍センサは、例えば、被測定者の体の表面に接触する一対の検出電極を備えるセンサである。心拍センサにより得られる心拍情報は、心臓の拍動により得られる電気信号であり、例えば、心電図である。心拍センサは、導電性粘着ゲル電極であってもよく、導電性繊維などで構成されるドライ電極であってもよい。心拍センサの装着部位は、胸部であり、心拍センサの形状は、例えば、ウエアと電極とが一体となったウエア型である。
　生体センサが脈波情報を取得するセンサ（以下、脈波センサ）である場合、脈波センサは、例えば、フォトトランジスタ及びフォトダイオードにより血管中の血液量の変化を反射光又は透過光により測定するセンサである。脈波センサは、ユーザの手首に装着され、当該装着された形状で脈波情報を測定する。脈波センサの装着部位は、足首、指、上腕などでもよい。脈波センサの形状は、バンド型（例えば、腕時計型）に限定されず、頸部等に貼り付ける貼付型、メガネ型などであってもよい。また、脈波センサは、顔又は手などの皮膚の色度の変化から脈波情報を測定して脈拍を算出する画像センサであってもよい。
　また、生体情報が呼吸数である場合、生体センサは、例えば、胸部又は腹部に巻き付ける圧力センサを備えるベルト型のセンサ、又は、鼻の下に取り付ける温度センサである。
　また、生体情報が血中酸素飽和度である場合、生体センサは、例えば、フォトトランジスタ及び２種類のフォトダイオードにより血管中の血液に含まれる飽和酸素濃度の変化を反射光又は透過光により測定するセンサである。
　また、生体情報が血圧である場合、生体センサは、例えば、圧力センサが付いたベルトを上腕部と、指先又は橈骨と、に巻き付けるセンサである。
　また、生体情報が体温である場合、生体センサは、例えば、掌又は鼻の頭などストレスにより毛細血管の収縮が起こりやすい部位に貼り付ける熱電対のセンサである。
　また、生体情報が発汗である場合、生体センサは、例えば、掌又は顔などストレスにより発汗が起こりやすい部位に接触する一対の検出電極を備えるセンサである。
　第１生体センサ１１１ａ及び第２生体センサ１１１ｂで測定された各生体信号は、第１信号処理部１１２ａ及び第２信号処理部１１２ｂに出力される。
　第１信号処理部１１２ａは、第１生体センサ１１１ａで測定された１つの生体信号から複数種類の生体指標を算出する。本実施の形態では、第１センサ１１１ａは、心拍センサである。上述したように、心拍の生体信号が心電図の場合、複数種類の生体指標は、ＲＲＩ、ＣｖＲＲ、ＨＦ及びＬＦなどである。ＲＲＩは、心拍数の指標であり、ＣｖＲＲ、ＨＦ及びＬＦは、心拍揺らぎの指標である。さらに、第１信号処理部１１２ａは、心拍揺らぎの周波数成分から呼吸数及び血圧の変動の生体指標を算出してもよい。また、これらの複数種類の生体指標のうち判定精度が比較的高い組み合わせは、ＲＲＩ及びＣｖＲＲである。したがって、本実施の形態では、生体指標１及び生体指標２は、それぞれＲＲＩ及びＣｖＲＲである例について説明する。なお、ＲＲＩ及びＣｖＲＲの算出方法については、モニター試験にて上述した通りである。第１信号処理部１１２ａは、算出された生体指標１及び生体指標２を演算部１２ａに出力する。
　また、第２信号処理部１１２ｂは、第２生体センサ１１１ｂで測定された１つの生体情報から少なくとも１種類の生体指標を算出する。本実施の形態では、生体指標３が算出される。上述したように、生体情報が発汗である場合、第２生体センサ１１１ｂは一対の検出電極を備えるセンサである。また、生体情報が体温である場合、第２生体センサ１１１ｂは、例えば、熱電対のセンサである。第２生体センサ１１１ｂは、例えば、被測定者の指に巻き付けて装着される。生体情報が発汗である場合、第２信号処理部１１２ｂは、皮膚コンダクタンスを算出する。また、第２生体センサ１１１ｂから出力された生体情報が体温である場合、第２信号処理部１１２ｂは、皮膚温度を算出する。したがって、本実施の形態では、生体指標３は、皮膚コンダクタンス又は皮膚温度である。第２信号処理部１１２ｂは、算出された生体指標３を演算部１２ａに出力する。
　演算部１２ａは、第１信号処理部１１２ａが出力した生体指標１及び生体指標２を取得し、取得した生体指標１及び生体指標２から生体指標１の変化量及び生体指標２の変化量を算出する。また、演算部１２ａは、第２信号処理部１１２ｂが出力した生体指標３を取得し、取得した生体指標３から生体指標３の変化量を算出する。生体指標の変化量は、被測定者の安静時に測定された生体指標（以下、基準値と称する場合がある。）を基準とする測定された生体指標であり、差又は比で表される。各生体指標の基準値は、記憶部１５ａに格納されている。演算部１２ａは、記憶部１５ａに格納された各生体指標の基準値を読み出し、当該基準値に対する各生体指標の変化量を算出する。演算部１２ａは、算出した各生体指標の変化量を判定部１３ａに出力する。なお、基準値は、季節又は被測定者の生理周期などにより変動する場合があるため、所定の期間毎に更新されてもよい。
　また、演算部１２ａは、各生体指標の閾値を算出する。生体指標１が、例えば、心拍数である場合、心拍数の変化量は、第１の時間に測定された心拍数の変化量である。第１の閾値は、生体指標１の閾値であり、例えば、心拍数の指標であるＲＲＩの閾値である。第１の閾値は、被測定者の安静時の心拍数を基準とする、任意の時間に測定された心拍数である。また、生体指標２が、例えば、心拍揺らぎである場合、心拍揺らぎの変化量は、第２の時間に測定された心拍揺らぎの変化量である。第２の閾値は、生体指標２の閾値であり、例えば、心拍揺らぎの指標であるＣｖＲＲの閾値である。第２の閾値は、被測定者の安静時の心拍数を基準とする、任意の時間に測定された心拍揺らぎである。また、生体指標３が、例えば、皮膚コンダクタンス又は皮膚温度である場合、皮膚コンダクタンス又は皮膚温度の変化量は、第３の時間に測定された皮膚コンダクタンス又は皮膚温度の変化量である。第３の閾値は、生体指標３の閾値であり、例えば、皮膚コンダクタンスの閾値又は皮膚温度の閾値である。第３の閾値は、被測定者の安静時の皮膚コンダクタンスを基準とする、任意の時間に測定された皮膚コンダクタンス、又は、被測定者の安静時の皮膚温度を基準とする、任意の時間に測定された皮膚温度である。これらの閾値は、第１、第２及び第３の時間とは異なる任意の時間に測定された生体指標の測定値と基準値との差又は比である生体指標の変化量である。ここで、任意の時間とは、例えば、被測定者がストレスを感じる手前の状態にある時を指す。
　以下、本実施の形態では、第１の時間、第２の時間、及び、第３の時間は、同じ時間である場合について説明するが、第１の時間、第２の時間、及び、第３の時間は、それぞれ異なる時間であってもよい。例えば、第１信号処理部１１２ａは、第１生体センサ１１１ａで測定された１つの生体信号から時分割で複数種類の心拍数及び心拍揺らぎを算出してもよい。このとき、演算部１２は、第１の時間に測定された心拍数の変化量を算出し、第１の時間とは異なる第２の時間に測定された心拍揺らぎの変化量を算出する。また、第２信号処理部１１２ｂは、第３の時間に、第２生体センサ１１２ｂで発汗又は皮膚温度を測定してもよい。このとき、演算部１２は、第３の時間に測定された皮膚コンダクタンスの変化量又は皮膚温度の変化量を算出する。なお、第３の時間は、第１の時間及び第２の時間のいずれかと同じ時間であってもよい。
　演算部１２ａは、記憶部１５ａに格納された各生体指標の閾値を読み出し、各生体指標の変化量の値と各生体指標の閾値との大小関係を比較する。そして、演算部１２ａは、各生体指標の変化量の少なくとも１つが閾値を一定時間超えている期間をストレス発生期間と判定する。ストレス発生期間とは、被測定者がストレスを感じた期間である。演算部１２ａは、ストレス発生期間中の各生体指標の変化量の値から各生体指標の変化量の代表値を算出する。例えば、ストレス発生期間における各生体指標の変化量の代表値は、ストレス発生期間中の各生体指標の変化量の平均値を用いてもよく、基準値からの差分が最も大きい値（最大値）を用いてもよい。
　判定部１３は、演算部１２ａが出力した各生体指標の変化量の代表値を取得し、記憶部１５ａに格納された第１の閾値、第２の閾値及び第３の閾値を読み出す。判定部１３ａは、生体指標１の変化量の代表値と第１の閾値との大小関係を比較し、かつ、生体指標２の変化量の代表値と第２の閾値との大小関係を比較し、かつ、生体指標３の変化量の代表値と第３の閾値との大小関係を比較することにより、被測定者のストレスの要因を判定する。つまり、判定部１３ａは、ストレス発生期間毎にストレスの要因を判定する。生体指標の変化量の代表値は、生体指標の変化量の一例であると言えるため、以下、生体指標の変化量の代表値を単に生体指標の変化量とも呼ぶ。
　具体的には、判定部１３ａは、生体指標１（ここでは、心拍数）の変化量が第１の閾値よりも大きく、かつ、生体指標２（ここでは、心拍揺らぎ）の変化量が第２の閾値よりも大きく、かつ、生体指標３（ここでは、皮膚コンダクタンス又は皮膚温度）の変化量が第３の閾値よりも大きい場合、ストレスの要因は、対人に関する要因であると判定する。また、判定部１３ａは、生体指標１の変化量が第１の閾値よりも大きく、かつ、生体指標２の変化量が第２の閾値よりも小さく、かつ、生体指標３の変化量が第３の閾値よりも大きい場合、ストレスの要因は、痛みであると判定する。また、判定部１３ａは、生体指標１の変化量が第１の閾値よりも小さく、かつ、生体指標２の変化量が第２の閾値よりも大きく、かつ、生体指標３の変化量が第３の閾値よりも小さい場合、ストレスの要因は、思考による疲労であると判定する。
　さらに、判定部１３ａは、生体指標１の変化量と第１の閾値との差、生体指標２の変化量と第２の閾値との差、及び、生体指標３の変化量と第３の閾値との差に応じて、ストレスの強度を判定し、判定結果を当該判定結果に基づく情報として出力する。判定結果に基づく情報は、例えば、ストレスの要因、ストレスの強度及びストレスの低減策の少なくとも１つを含む。ストレスの低減策は、例えば、ストレスの解消方法又はストレスの回避方法などである。ストレスの低減策は、後述する提示情報テーブルに含まれる。判定部１３ａは、記憶部１５ａに格納された提示情報テーブルから適切なストレス低減策を読み出し、提示部１４ａに出力する。
　また、判定部１３ａは、判定結果に基づく情報を記憶部１５ａに格納する。このとき、判定部１３ａは、被測定者がストレスを感じた時間の情報と上記判定結果に基づく情報とを紐づけして記憶部１５ａに格納してもよい。
　提示部１４ａは、判定部１３ａによって出力された上記判定結果に基づく情報を提示する。提示部１４ａは、上記判定結果に基づく情報を音声で提示してもよく、画像で提示してもよい。提示部１４ａが上記情報を音声で提示する場合は、提示部１４ａは、例えば、スピーカーである。また、提示部１４ａが上記情報を画像で提示する場合は、提示部１４ａは、例えば、ディスプレイである。
　記憶部１５ａは、各生体指標の基準値、各生体指標の閾値、及び、提示情報テーブルなどを格納する。提示情報テーブルは、ストレスの要因及び当該ストレスの強度に応じて提示されるストレス低減策などの提示情報のテーブルである。上述したように、各生体指標の基準値及び閾値は、所定の期間で更新されてもよい。なお、提示情報テーブルも同様に、所定の期間で更新されてもよい。
　また、記憶部１５ａは、判定部１３ａが出力したストレスの要因、ストレスの強度及びストレス低減策などの判定結果に基づく情報を格納する。このとき、記憶部１５ａは、判定結果に基づく情報とストレス発生期間とを紐付けて格納してもよい。これにより、被測定者は、所望のタイミングで判定結果に基づく情報を呼び出すことができる。このとき、判定部１３ａは、入力部１６ａにより入力された被測定者の操作に基づいて、判定結果に基づく情報を提示部１４に提示させる。
　入力部１６ａは、被測定者による操作を示す操作信号を判定部１３ａに出力する。入力部１６ａは、例えば、キーボード、マウス、タッチパネル、又は、マイクなどである。操作信号とは、判定結果に基づく情報の抽出方法又は提示部１４ａにおける提示方法などの設定を行う信号である。提示部１４ａには、入力部１６ａに入力された設定に基づき、様々な形式の判定結果が提示される。例えば、所定の期間におけるストレスの変化、被測定者が影響を受けやすいストレスの要因、及び、被測定者に適したストレス低減策などである。これにより、被測定者は、短期的なストレスの傾向を把握できるだけでなく、中期的及び長期的なストレスの傾向を把握することができる。このように、被測定者は、自己に適した効果的なストレス低減策を知ることができるため、中長期的なストレスを制御することができる。
　［ストレス評価方法］
　次に、本実施の形態に係るストレス評価方法について図１３を用いて具体的に説明する。図１３は、実施の形態に係るストレス評価方法を説明するフローチャートである。
　本実施の形態に係るストレス評価方法は、測定された被測定者の（ｉ）心拍数、（ｉｉ）心拍揺らぎ、及び、（ｉｉｉ）皮膚コンダクタンス又は皮膚温度を取得する取得ステップＳ１００と、（ｉ）心拍数の変化量、（ｉｉ）心拍揺らぎの変化量、及び、（ｉｉｉ）皮膚コンダクタンスの変化量、又は、皮膚温度の変化量を算出する算出ステップＳ２００と、（ｉ）心拍数の変化量と、（ｉｉ）心拍揺らぎの変化量と、（ｉｉｉ）皮膚コンダクタンスの変化量又は皮膚温度の変化量の少なくとも一方の変化量と、に基づいて被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定ステップＳ３００と、を含む。心拍数の変化量は、被測定者の安静時の心拍数である基準からの第１センサ部１１ａによって測定された心拍数への変化量であり、心拍揺らぎの変化量は、被測定者の安静時の心拍揺らぎである基準からの第１センサ部１１ａによって測定された心拍揺らぎへの変化量である。また、皮膚コンダクタンスの変化量は、被測定者の安静時の皮膚コンダクタンスである基準からの第２センサ部１１ｂによって測定された皮膚コンダクタンスへの変化量であり、皮膚温度の変化量は、被測定者の安静時の皮膚温度である基準からの第２センサ部１１ｂによって測定された皮膚温度である。判定ステップＳ３００では、（Ｉ）心拍数の変化量と第１の閾値との大小関係を比較し、かつ、（ＩＩ）心拍揺らぎの変化量と第２の閾値との大小関係を比較し、かつ、（ＩＩＩ）皮膚コンダクタンスの変化量又は皮膚温度の変化量と第３の閾値との大小関係を比較することにより、ストレスの要因を判定する。本実施の形態では、さらに、判定ステップＳ３００の判定結果に基づく情報を提示する提示ステップＳ４００を含む。
　以下、各ステップについてより具体的に説明する。
　まず、取得ステップＳ１００では、演算部１２ａは、第１センサ部１１ａ及び第２センサ部１１ｂで測定された被測定者の複数の生体指標を取得する。第１センサ部１１ａでは、第１生体センサ１１１ａで心拍情報（ここでは、心電図）が測定され、第１信号処理部１１２ａで、心拍数の指標及び心拍揺らぎの指標が算出される。また、第２センサ部１１ｂでは、第２生体センサ１１１ｂで温度又は発汗の生体情報が測定され、第２信号処理部１１２ｂで、皮膚温度（ＳＫＴ）又は皮膚コンダクタンス（ＳＣ）が算出される。なお、上述した通り、生体情報は、例えば、心拍、脈拍、呼吸数、血中酸素飽和度、血圧、体温、発汗などのストレスで影響を受ける生理学的な情報であってもよい。特に、心拍情報は、ウエアラブルな生体センサを用いた場合、脈拍、呼吸数、血圧、及び血中酸素飽和度などの他の生体情報よりも被測定者の負担が少ない状態で簡便に、かつ、リアルタイムに測定することができる。そのため、生体情報として被測定者の心拍情報を用いることにより、被測定者のストレスの状態を適切に評価することができる。
　例えば、心拍情報から得られる生体指標は、心拍数の指標であるＲＲＩ、心拍揺らぎの指標であるＣｖＲＲ、ＬＦ、ＨＦ、及び、ＬＦ／ＨＦ等などである。このように１つの生体情報から複数種類の生体指標が得られる。また、上述したように、これらの生体指標の組み合わせにより、比較的高い判定精度でストレスの要因を判定することができるため、信頼性の高い評価が得られる。
　図６を再び参照する。心拍情報は、例えば、心電図であり、図６に示す心電波形となる。心電波形は、心房の電気的興奮を反映するＰ波と、心室の電気的興奮を反映するＱ波、Ｒ波、及びＳ波と、興奮した心室の心筋細胞が再分極する過程を反映するＴ波とから構成されている。これらの心電波形のうち、Ｒ波の波高（電位差）が最も大きく、筋電位などのノイズに対して最も頑健である。このため、これらの心電波形における連続する２つの心拍のＲ波のピークの間隔、つまり、心拍間隔（ＲＲＩ）を算出する。心拍数は、ＲＲＩの逆数に６０を乗じて算出される。
　さらに、第１の知見におけるモニター試験で上述した通り、ＣｖＲＲは、上記式（２）を用いて、ＲＲＩから、任意時間帯におけるＲＲＩの標準偏差ＳＤを心拍間隔の平均値で規格化することにより算出される。
　第１信号処理部１１２ａは、第１生体センサ１１１ａで得られた心拍情報から、左心室が急激に収縮して心臓から血液を送り出す際に発生する電気信号（Ｒ波）を検出し、ＲＲＩを算出する。なお、Ｒ波の検出には、例えば、Ｐａｎ＆Ｔｏｍｐｋｉｎｓ法などの公知の手法が用いられる。
　次に、演算部１２ａにおいて検出されたＲ波から心拍間隔（ＲＲＩ）の変動量を算出する方法について説明する。
　図７を再び参照する。第１信号処理部１１２ａは、得られたＲ波の検出データから、以下のように、ＲＲＩの変動量を算出する。
　図７の（ａ）に示すように、第１信号処理部１１２ａは連続する２つの心拍のＲ波のピークの間隔であるＲＲＩを算出する。第１信号処理部１１２ａは、算出された各ＲＲＩを時間とＲＲＩとの２軸の関係に変換する。変換されたデータは不等間隔の離散的なデータであるため、演算部１２ａは、変換されたＲＲＩの時系列データを、図７の（ｂ）に示す等間隔時系列データに変換する。次いで、演算部１２ａは、この等間隔時系列データに対して、例えば、高速フーリエ変換（ＦＦＴ）を用いて周波数解析することにより、図７の（ｃ）に示す心拍変動の周波数成分を求める。
　心拍変動の周波数成分は、例えば、高周波成分ＨＦと低周波成分ＬＦとに分けることができる。モニター試験で上述した通り、ＨＦは、副交感神経活動量を反映していると考えられる。また、ＬＦは、交感神経及び副交感神経の活動量を反映すると考えられている。そのため、ＬＦとＨＦとの比であるＬＦ／ＨＦは、交感神経活動量を示すと考えられる。
　このように、第１センサ部１１ａでは、心拍情報から複数種類の生体指標を算出される。
　以上のように、取得ステップＳ１００では、演算部１２ａにて、第１センサ部１１ａから出力された２種類の生体指標（ここでは、心拍数及び心拍揺らぎ）及び第２センサ部１１ｂから出力された１種類の生体指標（ここでは、皮膚コンダクタンス）を取得する。
　次いで、算出ステップＳ２００では、演算部１２ａにて、取得ステップＳ１００で取得された各生体指標の変化量を算出する。各生体指標の変化量は、上述した通り、例えば、被測定者に安静時の各生体指標の値を基準値として、各生体指標の基準値と取得された各生体指標の値との比又は差を算出して得られる。演算部１２ａは、記憶部１５ａに格納された各生体指標の基準値を読み出して使用する。
　なお、各生体指標の変化量は、変化量が差で表される場合は、取得ステップＳ１００で取得された各生体指標の値から各生体指標基準値を差し引くことにより算出される。例えば、心拍数の変化量は、取得ステップＳ１００で取得された被測定者の心拍数の値から心拍数の基準値を差し引くことにより算出される。また、変化量が比で表される場合は、取得ステップＳ１００で取得された各生体指標の値を各生体指標の基準値で割ることにより算出される。例えば、心拍数の変化量は、取得ステップＳ１００で取得された被測定者の心拍数の値を心拍数の基準値で割ることにより算出される。
　以上のように、算出ステップＳ２０では、演算部１２ａにて、各生体指標の変化量を算出する。
　次いで、判定ステップＳ３００では、判定部１３ａにて、算出ステップＳ２００で算出された各生体指標の変化量に基づいてストレスの要因を判定する。判定部１３ａは、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、被測定者のストレスの要因を判定する。具体的には、判定ステップＳ３００では、判定部１３ａは、心拍数の変化量が第１の閾値よりも大きく、かつ、心拍揺らぎの変化量が第２の閾値よりも大きく、かつ、皮膚コンダクタンスの変化量又は皮膚温度の変化量が第３の閾値よりも大きい場合、ストレスの要因は、対人に関する要因であると判定する。また、判定部１３ａは、生体指標１の変化量が第１の閾値よりも大きく、かつ、生体指標２の変化量が第２の閾値よりも小さく、かつ、皮膚コンダクタンスの変化量又は皮膚温度の変化量が第３の閾値よりも大きい場合、ストレスの要因は、痛みであると判定する。また、判定部１３ａは、生体指標１の変化量が第１の閾値よりも小さく、かつ、生体指標２の変化量が第２の閾値よりも大きく、かつ、皮膚コンダクタンスの変化量又は皮膚温度の変化量が第３の閾値よりも小さい場合、ストレスの要因は、思考による疲労であると判定する。
　さらに、判定部１３ａは、生体指標１の変化量と第１の閾値との差、心拍揺らぎの変化量と第２の閾値との差、及び、皮膚コンダクタンスの変化量又は皮膚温度の変化量と第３の閾値との差に応じて、ストレスの強度を判定し、判定結果を当該判定結果に基づく情報として出力する。
　なお、第１の閾値は、心拍数の閾値であり、被測定者の安静時の心拍数を基準とする被測定者に対して任意の時間に測定された心拍数である。第２の閾値は、心拍揺らぎの閾値であり、被測定者の安静時の心拍揺らぎを基準とする任意の時間に測定された心拍揺らぎである。第３の閾値は、皮膚コンダクタンス又は皮膚温度の閾値であり、被測定者の安静時の皮膚コンダクタンス又は皮膚温度を基準とする任意の時間に測定された皮膚コンダクタンス又は皮膚温度である。これらの閾値は、演算部１２ａにて算出され、記憶部１５ａに格納される。判定部１３ａは、記憶部１５ａに格納された各生体指標の閾値を読み出して使用する。上述したように、任意の時間とは、例えば、被測定者がストレスを感じる手前の状態にある時を指す。
　各生体指標の閾値は、各生体指標の変化量が正の値である場合の閾値と、各生体指標の変化量が負の値である場合の閾値とが設定される。基準値は変化量のゼロ点である。各生体指標の変化量と閾値との大小関係は、以下のように比較される。生体指標の変化量が正の値である場合、生体指標の変化量と正の閾値との大小関係を比較する。また、生体指標の変化量が負の値である場合、生体指標の変化量の絶対値と負の閾値の絶対値との大小関係を比較する。なお、各生体指標の閾値は、固定値であってもよく、所定の期間で更新されてもよく、日々の測定に基づいて都度更新されてもよい。
　なお、閾値は、線形判別又は決定木等の比較的単純な機械学習によって算出されてもよい。これにより、被測定者に適した判定基準値及び閾値を設定できるため、ストレスの要因をより精度良く判定することができる。
　以上のように、判定ステップＳ３００では、各生体指標の変化量と各生体指標の閾値との大小関係を比較することにより、被測定者のストレスの要因を判定する。
　次いで、提示ステップＳ４００では、提示部１４ａにて、判定部１３ａで判定された判定結果に基づく情報を提示する。提示部１４ａは、判定結果に基づく情報を音声で提示してもよく、画像で提示してもよい。判定結果に基づく情報は、ストレスの要因、ストレスの強度及びストレスの低減策の少なくとも１つを含む。提示部１４ａは、被測定者が入力部１６ａで入力した設定に基づき、様々な形式の判定結果を表示する。
　［ストレス評価装置の使用例］
　次に、本実施の形態にストレス評価装置１００ａの使用例について具体的に説明する。図１４は、本実施の形態に係るストレス評価装置１００ａの使用例を説明する図である。
　図１４に示すように、ストレス評価装置１００ａは、第１センサ部１１ａの一部である第１生体センサ１１１ａ及び第２センサ部１１ｂの一部である第２生体センサ１１１ｂと、第１生体センサ１１１ａ及び第２生体センサ１１１ｂ以外の構成を含む評価端末２０とから構成される。被測定者は、第１生体センサ１１１ａを胸部の肌に接触するように装着し、心電図（ＥＣＧ）を測定する。第１生体センサ１１１ａは、導電性粘着ゲル電極であってもよく、導電性繊維などで構成されるドライ電極であってもよい。第１生体センサ１１１ａは、測定した心拍の電気信号を、通信により評価端末２０に送信する。
　また、第２生体センサ１１１ｂは、腕時計型センサであり、掌に貼り付けて使用するセンサ電極を備える。第２生体センサ１１１ｂは、センサ電極で測定された掌の皮膚電位を測定し、通信により、評価端末２０に送信する。さらに、第２生体センサ１１１ｂは、指先に貼り付けて使用する熱電対型センサを備えてもよい。これにより、第２生体センサ１１１ｂは、熱電対型センサで指先の温度を測定することができる。なお、第１生体センサ１１１ａ及び第２生体センサ１１１ｂと評価端末２０との通信方法は、Ｂｌｕｅｔｏｏｈ（登録商標）などの無線通信であってもよく、有線通信であってもよい。
　評価端末２０は、第１センサ部１１ａの第１信号処理部１１２ａ、第２センサ部１１ｂの第２信号処理部１１２ｂ、演算部１２ａ、判定部１３ａ、提示部１４ａ、記憶部１５ａ及び入力部１６ａを備える。第１信号処理部１１２ａ及び第２信号処理部１１２ｂは、それぞれ第１生体センサ１１１ａ及び第２生体センサ１１１ｂから通信により送信された生体信号を受信する。
　第１信号処理部１１２ａは、受信した心拍の電気信号から心拍数の指標であるＲＲＩ及び心拍揺らぎの指標であるＣｖＲＲを算出しこれらの生体指標を演算部１２ａに出力する。第２信号処理部１１２ｂは、受信した皮膚電位の信号から発汗の指標である皮膚コンダクタンス（ＳＣ）を算出し、ＳＣを演算部１２ａに出力する。なお、第２信号処理部１１２ｂは、第２生体センサ１１１ｂが皮膚温度を測定した場合は、第２生体センサ１１１ｂから皮膚温度の信号を受信し、体温の指標である皮膚温度（ＳＫＴ）を算出し、ＳＫＴを演算部１２ａに出力する。
　演算部１２ａは、第１信号処理部１１２ａが出力したＲＲＩ及びＣｖＲＲを取得し、記憶部１５に格納されたＲＲＩの基準値及びＣｖＲＲの基準値を読み出す。また、演算部１２ａは、第２信号処理部１１２ｂが出力したＳＣを取得し、記憶部１５ａに格納されたＳＣの基準値を読み出す。演算部１２ａは、読み出された基準値を基準とする、これらの生体指標それぞれの変化量を算出する。生体指標の変化量は、差又は比で表される。本実施の形態では、当該変化量は、比で表される。
　また、上述した通り、演算部１２ａは、各生体指標の閾値を算出し、記憶部１５ａに出力する。各生体指標の閾値は、各生体指標の変化量が正の値になる場合の閾値と、各生体指標の変化量が負の値になる場合の閾値とが設定される。基準値は変化量ゼロである。具体的には、各生体指標の変化量が正の値になる場合、正の閾値は、基準値よりも大きい値であり、変化量のグラフ１２０ａ中の第１の閾値１ａ（以下、正の閾値１ａ）、第２の閾値２ａ（以下、正の閾値２ａ）及び第３の閾値３ａ（以下、正の閾値３ａ）である。各生体指標の変化量が負の値になる場合、負の閾値は、基準値よりも小さい値であり、変化量のグラフ１２０中の第１の閾値１ｂ（以下、負の閾値１ｂ）、第２の閾値２ｂ（以下、負の閾値２ｂ）及び第３の閾値３ｂ（以下、負の閾値３ｂ）である。また、演算部１２ａは、各生体指標の基準値を算出し、記憶部１５ａに出力する。各生体指標の基準値は、各生体指標の変化量がゼロである。例えば、変化量のグラフ１２０ａでは、基準値は、正の閾値１ａ及び負の閾値１ｂの間の実線で示される。なお、正の閾値及び負の閾値は、基準値（変化量ゼロ）を挟んで等間隔に設定されてもよく、基準値を挟んで等間隔に設定されなくてもよい。これらの閾値は、各生体指標の変化量の大きさに応じて、適宜設定されてもよい。
　判定部１３ａでは、演算部１２ａが出力した各生体指標の変化量を取得し、記憶部１５ａに格納された各生体指標の閾値を読み出す。判定部１３ａは、各生体指標の変化量と各生体指標の閾値との大小関係を比較し、ストレスの要因を判定する。例えば、各生体指標の変化量が正の値である場合、判定部１３ａは、各生体指標の変化量と正の閾値との大小関係を比較する。また、各生体指標の変化量が負の値である場合、判定部１３ａは、各生体指標の変化量の絶対値と負の閾値の絶対値との大小関係を比較する。以下、変化量のグラフ１２０ａ及び判定表１３０ａを用いて、より具体的に説明する。
　変化量のグラフ１２０ａに示すように、期間Ａ２では、ＲＲＩの変化量の絶対値は、負の閾値１ｂの絶対値よりも大きく、かつ、ＣｖＲＲの変化量は、正の閾値２ａよりも大きく、かつ、皮膚コンダクタンスの変化量は、正の閾値３ａよりも大きい。よって、判定部１３ａは、被測定者が期間Ａ２で感じたストレスの要因は、対人に関する要因であると判定する。また、期間Ｂ２では、ＲＲＩの変化量は、正の閾値１ａよりも大きく、かつ、ＣｖＲＲの変化量の絶対値は、負の閾値２ｂの絶対値よりも小さく、かつ、皮膚コンダクタンスの変化量は、正の閾値３ａよりも大きい。よって、判定部１３ａは、被測定者が期間Ｂ２で感じたストレスの要因は、痛みであると判定する。また、期間Ｃ２では、ＲＲＩの変化量の絶対値は、負の閾値１ｂの絶対値よりも小さく、かつ、ＣｖＲＲの変化量の絶対値は、負の閾値２ｂの絶対値よりも大きく、かつ、皮膚コンダクタンスの変化量の絶対値は、負の閾値３ｂの絶対値よりも小さい。よって、判定部１３ａは、被測定者が期間Ｃ２で感じたストレスの要因は、思考による疲労（思考疲労）であると判定する。
　判定表１３０ａでは、矢印の向き及び本数で、基準値（変化量ゼロ）に基づく各生体指標の変化量の推移を示している。横向きの矢印は、生体指標の変化量が閾値を超える変化を伴わないことを示している。
　さらに、判定部１３ａは、ＲＲＩの変化量の絶対値と第１の閾値の絶対値との差、ＣｖＲＲの変化量の絶対値と第２の閾値の絶対値との差、及び、ＳＣの変化量の絶対値と第３の閾値の絶対値との差に応じて、ストレスの強度を判定する。
　判定部１３ａは、これらの判定結果に基づく情報を提示部１４ａに出力する。提示部１４ａは、例えば、スマートフォンのディスプレイである。また、判定部１３ａは、判定結果に基づく情報を記憶部１５ａに格納する。これにより、被測定者は、所望のタイミングで判定結果に基づく情報を呼び出すことができる。このとき、判定部１３ａは、タッチパネルなどの入力部１６ａにより入力された被測定者の操作に基づいて、判定結果に基づく情報を提示部１４ａに提示させる。例えば、被測定者が評価端末２０の入力部１６ａで必要な情報を抽出する指示を入力すると、判定部１３ａは、被測定者の指示に基づいて提示部１４ａに提示情報１４０ａを提示する。提示情報１４０ａは、被測定者がストレスを感じた時間、ストレスの要因、及びストレスの低減策を含んでいる。ストレスの低減策は、例えば、ストレスの要因に応じたストレス解消方法又はストレス回避方法を提案するメッセージである。当該メッセージは、例えば、ストレスの要因が思考疲労である場合、少し休憩しましょう、又は、ストレッチをしましょう、などであり、対人に関わる要因である場合、少し瞑想しましょう、又は、深呼吸をしましょう、などである。
　以上のように、本実施の形態によれば、被測定者が日常生活を送りながら簡便に、かつ、正確に、ストレスの要因を判定することができる。そのため、被測定者は、従来よりも正確に自身のストレス状態及び適切なストレス低減策を把握することができる。これにより、被測定者は、適切に、かつ、効率良く、自身のストレスの制御を行うことができるため、ストレスの制御を継続して行うことができる。
　以上、本開示に係るストレス評価装置、ストレス評価方法及びプログラムについて、実施の形態に基づいて説明したが、本開示は、これらの実施の形態に限定されるものではない。本開示の主旨を逸脱しない限り、当業者が思いつく各種変形を実施の形態に施したものや、本実施の形態における一部の構成要素を組み合わせて構築される別の形態も、本開示の範囲に含まれる。
　なお、上記実施の形態では、生体情報として心拍情報を用い、心拍情報から得られる複数種類の生体指標として、心拍数の指標及び心拍揺らぎの指標を用いる例を示したが、これに限られない。例えば、自律神経活動度であるエントロピーＥ及び自律神経バランスであるトーンＴを用いてもよい。また、上記実施の形態では、心拍数の指標として、ＲＲＩを用い、心拍揺らぎの指標としてＣｖＲＲ、ＬＦ及びＨＦを用いる例を説明したが、心拍揺らぎを示すこれら以外の指標を用いてもよい。
　また、実施の形態１では、ストレス評価装置１００が生体センサ１１１と評価端末２０とから構成される例を示したが、例えば、第１センサ部１１ａと、第１センサ部１１ａ以外の構成を備える評価端末とから構成されてもよい。
　また、実施の形態２では、ストレス評価装置１００ａが生体センサ１１１ａと評価端末２０とから構成される例を示したが、例えば、第１センサ部１１ａ及び第２センサ部１１ｂと、第１センサ部１１ａ及び第２センサ部１１ｂ以外の構成を備える評価端末とから構成されてもよい。
　また、ストレス評価装置は、全ての構成が１つのデバイスに組み込まれた一体型の装置であってもよい。本実施の形態では、生体センサは心拍センサである例を示したが、生体センサは脈波センサであってもよい。この場合、ストレス評価装置は、ディスプレイを備える腕時計型の装置であってもよい。
　また、実施の形態１では、評価端末２０はスマートフォン又はタブレット端末である例を示したが、スマートフォン又はタブレット端末は提示部１４と入力部１６とを備え、第１信号処理部１１２ａ、演算部１２、判定部１３及び記憶部１５をインターネットなどの通信ネットワークを介して接続されるサーバー上に設けてもよい。
　また、実施の形態２では、評価端末２０はスマートフォン又はタブレット端末である例を示したが、スマートフォン又はタブレット端末は提示部１４ａと入力部１６ａとを備え、第１信号処理部１１２ａ、第２信号処理部１１２ｂ演算部１２ａ、判定部１３ａ及び記憶部１５ａをインターネットなどの通信ネットワークを介して接続されるサーバー上に設けてもよい。
　また、各生体指標の基準値及び閾値が評価端末に設けられた記憶部に格納される形態を一例として示したが、上記基準値及び閾値がインターネット上のサーバーに格納され、評価端末に随時送られてくる形態でもよい。
　また、本開示では、ストレスの要因を判定するための指標の１つとして、皮膚コンダクタンスを挙げたが、精神性発汗を測定できる指標であれば特に限定されない。例えば、皮膚抵抗など皮膚電位又は電流値を測定して得られる指標であってもよく、皮膚表面の湿度など水分量を測定して得られる指標であってもよい。
　また、実施の形態２では、皮膚コンダクタンス又は皮膚温度を掌で測定する例を挙げたが、精神性発汗が生じやすい顔面の一部で測定してもよく、足の裏で測定してもよい。
　また、本開示では、ストレスの要因の一つである対人に関する要因の具体例として、モニター試験における模擬の就職面談を挙げたが、これに限られない。例えば、対人に関する要因は、職場及びプライベートでの人間関係、人前で話をすること、又は、人との交渉など、人と関わる事柄で被測定者が不安又は緊張を感じる要因であればよい。
　また、本開示では、ストレスの要因の一つである痛みの具体例として、電気刺激による痛みを挙げたが、これに限られない。例えば、痛みは、打撲、頭痛、歯痛、裂傷などの身体の痛み、又は、擦る、刺す、切る、打つなどの物理的な刺激に伴う痛みなど、恐怖又は我慢を感じる痛みであればよい。
　また、本開示では、ストレスの要因の一つである思考による疲労の具体例として、思考が必要な作業として暗算及び音声によるじゃんけんの課題を挙げたが、これに限られない。例えば、思考による疲労は、思考が必要な作業として、パソコンでの作業、又は、集中力が必要な実験等の知的な活動など、思考する作業を続けることによって疲労を感じる要因であればよい。
　本開示は、複数種類の生体指標から被測定者のストレスの要因を簡便に、かつ、正確に判定できるストレス評価装置として有用である。
　１１ａ　第１センサ部
　１１ｂ　第２センサ部
　１２、１２ａ　演算部
　１３、１３ａ　判定部
　１４、１４ａ　提示部
　１５、１５ａ　記憶部
　１６、１６ａ　入力部
　２０　評価端末
　１００、１００ａ　ストレス評価装置
　１１１ａ　第１生体センサ
　１１１ｂ　第２生体センサ
　１１２ａ　第１信号処理部
　１１２ｂ　第２信号処理部
　１２０、１２０ａ　変化量のグラフ
　１３０、１３０ａ　判定表
　１４０、１４０ａ　提示情報

Claims (19)

  1. 　被測定者の心拍数及び心拍揺らぎを測定する第１センサ部と、
    　（ｉ）心拍数の変化量、及び、（ｉｉ）心拍揺らぎの変化量を算出する演算部と、
    　（ｉ）前記心拍数の変化量及び（ｉｉ）前記心拍揺らぎの変化量に基づいて前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定部と、
    　を備え、
    　前記心拍数の変化量は、前記被測定者の安静時の心拍数である基準からの前記第１センサ部によって測定された前記心拍数への変化量であり、
    　前記心拍揺らぎの変化量は、前記被測定者の安静時の心拍揺らぎである基準からの前記第１センサ部によって測定された前記心拍揺らぎへの変化量であり、
    　前記判定部は、
    　（Ｉ）前記心拍数の変化量と第１の閾値との大小関係の比較、及び、
    　（ＩＩ）前記心拍揺らぎの変化量と第２の閾値との大小関係の比較を行うことにより、前記ストレスの要因を判定する、
    　ストレス評価装置。
  2. 　前記心拍数の変化量は、第１の時間に測定された前記心拍数の変化量であり、
    　前記心拍揺らぎの変化量は、第２の時間に測定された前記心拍揺らぎの変化量であり、
    　前記第１の閾値は、前記被測定者の安静時の心拍数を基準とする、前記第１及び前記第２の時間とは異なる任意の時間に測定された前記心拍数であり、
    　前記第２の閾値は、前記被測定者の安静時の心拍揺らぎを基準とする、前記任意の時間に測定された前記心拍揺らぎである、
    　請求項１に記載のストレス評価装置。
  3. 　前記判定部は、
    　前記心拍数の変化量が前記第１の閾値よりも大きく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも大きい場合、前記ストレスの要因は、対人に関する要因であると判定する、
    　請求項１又は２に記載のストレス評価装置。
  4. 　前記判定部は、
    　前記心拍数の変化量が前記第１の閾値よりも大きく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも小さい場合、前記ストレスの要因は、痛みであると判定する、
    　請求項１又は２に記載のストレス評価装置。
  5. 　前記判定部は、
    　前記心拍数の変化量が前記第１の閾値よりも小さく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも大きい場合、前記ストレスの要因は、思考による疲労であると判定する、
    　請求項１又は２に記載のストレス評価装置。
  6. 　さらに、前記判定部は、前記心拍数の変化量と前記第１の閾値との差、及び、前記心拍揺らぎの変化量と前記第２の閾値との差に応じて、前記ストレスの強度を判定し、判定結果を前記判定結果に基づく前記情報として出力する、
    　請求項１又は２に記載のストレス評価装置。
  7. 　さらに、前記被測定者の皮膚コンダクタンス又は皮膚温度の少なくとも一方を測定する第２センサ部を備え、
    　前記演算部は、さらに、
    　（ｉｉｉ）皮膚コンダクタンスの変化量、又は、皮膚温度の変化量を算出し、
    　前記皮膚コンダクタンスの変化量は、前記被測定者の安静時の皮膚コンダクタンスである基準からの前記第２センサ部によって測定された前記皮膚コンダクタンスへの変化量であり、
    　前記皮膚温度の変化量は、前記被測定者の安静時の皮膚温度である基準からの前記第２センサ部によって測定された前記皮膚温度への変化量であり、
    　前記判定部は、前記（Ｉ）及び前記（ＩＩ）に加えて、
    　（ＩＩＩ）前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量と第３の閾値との大小関係の比較も行うことにより、前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する、
    　請求項１に記載のストレス評価装置。
  8. 　前記心拍数の変化量は、第１の時間に測定された前記心拍数の変化量であり、
    　前記心拍揺らぎの変化量は、第２の時間に測定された前記心拍揺らぎの変化量であり、
    　前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量は、第３の時間に測定された皮膚コンダクタンス又は前記皮膚温度の変化量であり、
    　前記第１の閾値は、前記被測定者の安静時の心拍数を基準とする、前記第１、前記第２及び前記第３の時間とは異なる任意の時間に測定された前記心拍数であり、
    　前記第２の閾値は、前記被測定者の安静時の心拍揺らぎを基準とする、前記任意の時間に測定された前記心拍揺らぎであり、
    　前記第３の閾値は、前記被測定者の安静時の皮膚コンダクタンスを基準とする、前記任意の時間に測定された前記皮膚コンダクタンス、又は、前記被測定者の安静時の皮膚温度を基準とする、前記任意の時間に測定された前記皮膚温度である、
    　請求項７に記載のストレス評価装置。
  9. 　前記判定部は、
    　前記心拍数の変化量が前記第１の閾値よりも大きく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも大きく、かつ、前記皮膚コンダクタンス又は前記皮膚温度の変化量が前記第３の閾値よりも大きい場合、前記ストレスの要因は、対人に関する要因であると判定する、
    　請求項７又は８に記載のストレス評価装置。
  10. 　前記判定部は、
    　前記心拍数の変化量が前記第１の閾値よりも大きく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも小さく、かつ、前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量が前記第３の閾値よりも大きい場合、前記ストレスの要因は、痛みであると判定する、
    　請求項７～９のいずれか１項に記載のストレス評価装置。
  11. 　前記判定部は、
    　前記心拍数の変化量が前記第１の閾値よりも小さく、かつ、前記心拍揺らぎの変化量が前記第２の閾値よりも大きく、かつ、前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量が前記第３の閾値よりも小さい場合、前記ストレスの要因は、思考による疲労であると判定する、
    　請求項７～９のいずれか１項に記載のストレス評価装置。
  12. 　さらに、前記判定部は、
    　前記心拍数の変化量と前記第１の閾値との差、前記心拍揺らぎの変化量と前記第２の閾値との差、及び、前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量と前記第３の閾値との差に応じて、前記ストレスの強度を判定し、判定結果を前記判定結果に基づく前記情報として出力する、
    　請求項７～９のいずれか１項に記載のストレス評価装置。
  13. 　前記心拍揺らぎは、前記被測定者の心拍間隔を周波数分析して求められる、
    　請求項１～１２のいずれか１項に記載のストレス評価装置。
  14. 　さらに、前記判定部によって出力された前記判定結果に基づく前記情報を提示する提示部を備え、
    　前記情報は、前記ストレスの要因、前記ストレスの強度及び前記ストレスの低減策からなる群から選択される少なくとも１つを含む、
    　請求項１～１３のいずれか１項に記載のストレス評価装置。
  15. 　前記提示部は、音声で提示する、
    　請求項１４に記載のストレス評価装置。
  16. 　前記提示部は、画像で提示する、
    　請求項１４に記載のストレス評価装置。
  17. 　測定された被測定者の心拍数及び心拍揺らぎを取得する取得ステップと、
    　（ｉ）心拍数の変化量、及び、（ｉｉ）心拍揺らぎの変化量を算出する算出ステップと、
    　前記心拍数の変化量及び前記心拍揺らぎの変化量に基づいて前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する判定ステップと、
    　を含み、
    　前記心拍数の変化量は、前記被測定者の安静時の心拍数である基準からの前記第１センサ部によって測定された前記心拍数への変化量であり、
    　前記心拍揺らぎの変化量は、前記被測定者の安静時の心拍揺らぎである基準からの前記第１センサ部によって測定された前記心拍揺らぎへの変化量であり、
    　前記判定ステップでは、
    　（Ｉ）前記心拍数の変化量と第１の閾値との大小関係を比較し、かつ、（ＩＩ）前記心拍揺らぎの変化量と第２の閾値との大小関係を比較することにより、前記ストレスの要因を判定する、
    　ストレス評価方法。
  18. 　前記取得ステップは、さらに、
    　前記被測定者の皮膚コンダクタンス又は皮膚温度の少なくとも一方を取得し、
    　前記算出ステップは、さらに、
    　（ｉｉｉ）皮膚コンダクタンスの変化量、又は、皮膚温度の変化量を算出し、
    　前記皮膚コンダクタンスの変化量は、前記被測定者の安静時の皮膚コンダクタンスである基準からの前記第２センサ部によって測定された前記皮膚コンダクタンスへの変化量であり、
    　前記皮膚温度の変化量は、前記被測定者の安静時の皮膚温度である基準からの前記第２センサ部によって測定された前記皮膚温度であり、
    　前記判定ステップは、
    　前記（Ｉ）かつ前記（ＩＩ）かつ（ＩＩＩ）前記皮膚コンダクタンスの変化量又は前記皮膚温度の変化量と第３の閾値との大小関係を比較することにより、前記被測定者のストレスの要因を判定し、判定結果に基づく情報を出力する、
    　請求項１７に記載のストレス評価方法。
  19. 　請求項１７又は１８に記載のストレス評価方法をコンピュータに実行させるためのプログラム。
PCT/JP2019/016437 2018-05-30 2019-04-17 ストレス評価装置、ストレス評価方法及びプログラム WO2019230235A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201980017146.8A CN111818850A (zh) 2018-05-30 2019-04-17 压力评价装置、压力评价方法以及程序
US17/029,128 US20210000355A1 (en) 2018-05-30 2020-09-23 Stress evaluation device, stress evaluation method, and non-transitory computer-readable medium

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
JP2018103049 2018-05-30
JP2018103066 2018-05-30
JP2018-103049 2018-05-30
JP2018-103066 2018-05-30
JP2019067219A JP7262006B2 (ja) 2018-05-30 2019-03-29 ストレス評価装置、ストレス評価方法及びプログラム
JP2019-067219 2019-03-29

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US17/029,128 Continuation US20210000355A1 (en) 2018-05-30 2020-09-23 Stress evaluation device, stress evaluation method, and non-transitory computer-readable medium

Publications (1)

Publication Number Publication Date
WO2019230235A1 true WO2019230235A1 (ja) 2019-12-05

Family

ID=68697456

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/016437 WO2019230235A1 (ja) 2018-05-30 2019-04-17 ストレス評価装置、ストレス評価方法及びプログラム

Country Status (1)

Country Link
WO (1) WO2019230235A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2021182250A1 (ja) * 2020-03-11 2021-09-16
WO2023162645A1 (ja) * 2022-02-22 2023-08-31 ソニーグループ株式会社 情報処理装置、情報処理方法、プログラム、および情報処理システム

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015532165A (ja) * 2012-10-23 2015-11-09 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． ストレス測定システム
JP2017063966A (ja) * 2015-09-29 2017-04-06 シチズン時計株式会社 疲労度計
JP2018037073A (ja) * 2016-08-29 2018-03-08 パナソニックＩｐマネジメント株式会社 ストレスマネジメントシステム及びストレスマネジメント方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015532165A (ja) * 2012-10-23 2015-11-09 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． ストレス測定システム
JP2017063966A (ja) * 2015-09-29 2017-04-06 シチズン時計株式会社 疲労度計
JP2018037073A (ja) * 2016-08-29 2018-03-08 パナソニックＩｐマネジメント株式会社 ストレスマネジメントシステム及びストレスマネジメント方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2021182250A1 (ja) * 2020-03-11 2021-09-16
WO2021182250A1 (ja) * 2020-03-11 2021-09-16 パナソニックＩｐマネジメント株式会社 信号処理システム、センサシステム、生体管理システム、環境制御システム、信号処理方法、及びプログラム
WO2023162645A1 (ja) * 2022-02-22 2023-08-31 ソニーグループ株式会社 情報処理装置、情報処理方法、プログラム、および情報処理システム

Similar Documents

Publication Publication Date Title
JP7262006B2 (ja) ストレス評価装置、ストレス評価方法及びプログラム
US20190159685A1 (en) Biological information analyzing device, system, and program
JP4243605B2 (ja) 自律神経検査装置
US10034612B2 (en) Biological state eliminating apparatus and method
US20180199893A1 (en) Processing biological data
JP7191159B2 (ja) コンピュータプログラム、及び、被験者の感情状態を提供する方法
JP6579890B2 (ja) 疲労度計
US20160338640A1 (en) Psychological acute stress measurement using a wireless sensor
Choi et al. Removal of respiratory influences from heart rate variability in stress monitoring
Hanson et al. Vagal cardiac control throughout the day: the relative importance of effort–reward imbalance and within-day measurements of mood, demand and satisfaction
JP6513005B2 (ja) 疲労度計
JP5327458B2 (ja) 精神ストレス評価とそれを用いた装置とそのプログラム
JPH0795966A (ja) 脈波解析装置
EP2408358A2 (en) Stress monitor system and method
JP6676499B2 (ja) 疲労度判定装置、疲労度判定方法、疲労度判定プログラム、及び、生体情報測定装置
JP2011115188A (ja) 睡眠状態モニタリング装置、モニタリングシステムおよびコンピュータプログラム
JP2019030389A (ja) 自律神経状態評価装置、自律神経状態評価システム、自律神経状態評価方法及びプログラム
KR100610813B1 (ko) 혈압 측정 회로가 내장된 배터리 팩과 결합 가능한 이동통신 단말기
WO2019230235A1 (ja) ストレス評価装置、ストレス評価方法及びプログラム
CN116636821A (zh) 一种血压管理系统及方法
US11439347B2 (en) Portable dehydration monitoring system
JP2019025253A (ja) 自律神経状態評価装置、自律神経状態評価システム、自律神経状態評価方法及びプログラム
JP7402495B2 (ja) 健康情報検出装置の作動方法及び健康情報検出装置
JP2001204696A (ja) 電子血圧計
US20230139627A1 (en) Method for determining psychological stress

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19811545

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19811545

Country of ref document: EP

Kind code of ref document: A1