JP6927491B2 - Method and device for creating indicators to determine neuropsychiatric status - Google Patents

Description

The present invention relates to a method for creating an index capable of objectively discriminating a neuropsychiatric state and a device for creating the index.

Depression (major depressive disorder), bipolar disorder, schizophrenia, panic disorder, obsessive-compulsive disorder, autonomic imbalance, sleep disorder, etc. due to various factors such as inheritance, trauma, surrounding environment, stress, etc. May cause mental disorders. Therefore, observing the state of psychiatric nerves is effective for early detection of illness and appropriate treatment. The presence, type or severity of mental illness may be determined by a physician interview based on known criteria. For example, for the determination of depression, a depression self-evaluation scale (CES-D), a stress self-evaluation scale (Hoppins Symptom Checklist: HSCL), a depression evaluation scale (HAMD), a Young Mania Rating Scale (YMRS), etc. Criteria are used. However, since the conventional judgment criteria include abstract items and there is room for the subjectivity of the patient to enter, a method that can be judged objectively and mechanically is desired.

For example, US Pat. The high frequency index, low frequency index, and their ratios in the frequency analysis step of the method are measured in each of the three states of the subject in the resting state, the task performing state, and the resting state after the task is performed, and the frequency is measured. Obtained by analysis. However, the evaluation method described in Patent Document 1 assigns a task such as random number generation to the subject, and may impose an excessive burden on the subject. Therefore, the neuropsychiatric state is deteriorated. There is room for improvement from the perspective of easily creating indicators for discrimination.

Japanese Unexamined Patent Publication No. 2013-233256

Therefore, an object of the present invention is to provide a method for creating an index and a device for creating the index, which can easily and objectively determine the neuropsychiatric state.

The present inventors are investigating a method for discriminating a neuropsychiatric state using easily measurable biometric information, and when the subject is in a depressive posture other than the sleep time zone or during awakening. Finding that (LF / HF) / activity is related to the results of interviews with subjects such as the depression self-evaluation scale (CES-D) and stress self-evaluation scale (HSCL), and discriminating the neuropsychiatric status. The present invention was completed with the idea that it is effective as an index for doing so.

That is, the method for creating an index for discriminating the neuropsychiatric state of the present invention (hereinafter, may be simply referred to as "first method") includes a step of measuring the acceleration TA in the height direction of the subject and the following conditions. The steps of calculating the negative acceleration NA from (1) to (4) to calculate the first lying time zone, the second lying time zone, the sleeping time zone, and the third lying time zone, and the constant term. , A variable term for the total time of the third recumbent time zone, and a step of creating an index represented by an equation including the following.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _{m2 : The gap time zone L sm in} which there are two or more of the first recumbent time zones L _m1 and NA <C1 between two adjacent first recumbent time zones L _m11 and L _m12. If it is within the second predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled (3) Calculation of the sleep time zone Sleep time zone: The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time (4) Calculation of the third recumbent time zone Third recumbent time zone: The time zone obtained by excluding the time of the sleep time zone from the total time of the first recumbent time zone L _m1 and the second recumbent time zone L _m2]

Since the first method has the conditions (1) to (4), the sleep time zone can be calculated objectively and accurately. In addition, an index for discriminating the neuropsychiatric state can be easily created from the time of the lying posture of the subject, which is a parameter that can be easily measured. Furthermore, it is possible to objectively determine the neuropsychiatric state by using the created index.

In addition, the method for creating an index for discriminating the neuropsychiatric state of the present invention (hereinafter, may be simply referred to as "second method") is the beat interval of the subject and the acceleration associated with the movement of the subject. Alternatively, the power spectrum obtained by including the step of measuring the amount of activity expressed by the angular velocity and the acceleration TA in the height direction of the subject and the step of converting the beat interval into a frequency spectrum is determined from frequencies Lf1 to Lf2. The first step is to calculate the LF, which is the integrated value, and the HF, which is the value obtained by constantly integrating the power spectrum from frequencies Hf1 to Hf2, and to calculate the negative acceleration NA under the following conditions (1) to (3). An expression that includes a step for calculating the lying time zone, the second lying time zone, the awakening time zone, and the sleeping time zone, a constant term, and a variable term of (LF / HF) / activity amount of the awakening time zone. It has a gist in that it has a step to create an index to be represented. Here, Hf1> Lf1 and Hf2> Lf2.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _m2 : There are two or more first _decubitus time zones L _m1 , and the gap time zone L _{sm where NA <C1 between two adjacent first decubitus time zones L m11} and L m12 is the first. 2 If it is within the predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled (3) Calculation of the awakening time zone and the sleeping time zone Sleeping time zone : The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time Awakening time zone: The time zone excluding the sleep time zone from the predetermined measurement unit time]

Since the second method has the conditions (1) to (3), it is possible to objectively and accurately calculate the awakening time zone and the sleeping time zone. In addition, an index for discriminating the neuropsychiatric state can be easily created from the easily measurable parameters of the awakening time zone (LF / HF) and the amount of activity. Furthermore, it is possible to objectively determine the neuropsychiatric state by using the created index.

In the second method, it is preferable that the index includes the variable term of sleep time zone (LF / HF) × activity amount.

The second method further includes a step of calculating a third recumbent time zone obtained by subtracting the sleep time zone from the total time of the first recumbent time zone and the second recumbent time zone, and the index is , It is preferable to include the variable term of the total time of the third recumbent time zone.

In the second method, it is preferable that the index includes a variable term of pulsation interval × activity amount in the awakening time zone.

In the second method, it is preferable that the index includes the variable term of HF × activity amount in the awakening time zone.

In the second method, it is preferable that the index includes the variable term of the beat interval / activity amount during the sleep time zone.

In the second method, it is preferable that the index includes at least one of the variable term of HF / activity in the sleep time zone and the variable term of the time in the sleep time zone.

The second method may further include the step of time zone is NA ≧ C1 calculates the fourth recumbent hours is within a third predetermined time T _3, the third predetermined time T _3, the first It is preferable that the predetermined time _{is less than T 1} and the index includes the variable term of the total time of the fourth recumbent time zone.

In the first method or the second method, it is preferable to use the RR interval, which is the interval between the R wave and the R wave in the electrocardiographic signal, as the pulsation interval.

In the first method or the second method, it is preferable that the index indicates the presence or absence or degree of depression or autonomic imbalance.

The device for creating an index for discriminating the neuropsychiatric state of the present invention (hereinafter, may be simply referred to as "first device") has a time calculation unit for calculating a sleep time zone and an index for discriminating the neuropsychiatric state. It has an index creation unit to be created, and the time calculation unit calculates the negative acceleration NA from the acceleration TA in the height direction of the subject under the following conditions (1) to (4), and is the first lie. It calculates the position time zone, the second lying time zone, the sleeping time zone, and the third lying time zone, and the index is a constant term and a variable term of the total time of the third lying time zone. It has a gist in the point expressed by the including formula.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _{m2 : The gap time zone L sm in} which there are two or more of the first recumbent time zones L _m1 and NA <C1 between two adjacent first recumbent time zones L _m11 and L _m12. If it is within the second predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled (3) Calculation of the sleep time zone Sleep time zone: The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time (4) Calculation of the third recumbent time zone Third recumbent time zone: The time zone obtained by excluding the time of the sleep time zone from the total time of the first recumbent time zone L _m1 and the second recumbent time zone L _m2]

The time calculation unit of the first device can objectively and accurately calculate the sleep time zone according to the conditions (1) to (4). In addition, the index creation unit can easily create an index for discriminating the neuropsychiatric state from the time of the lying posture of the subject, which is a parameter that can be easily measured. Furthermore, it is possible to objectively determine the neuropsychiatric state by using the created index.

Further, the device for creating an index for discriminating the neuropsychiatric state of the present invention (hereinafter, may be simply referred to as a "second device") obtains a power spectrum obtained by including a step of converting the beat interval into a frequency spectrum. An LF that is a constant integration value from frequencies Lf1 to Lf2, a processing unit that calculates HF that is a constant integration value of a power spectrum from frequencies Hf1 to Hf2, and a time calculation unit that calculates an awakening time zone and a sleep time zone. The time calculation unit has an index creation unit that creates an index for discriminating the neuropsychiatric state, and the time calculation unit has a negative acceleration from the acceleration TA in the height direction of the subject under the following conditions (1) to (3). NA is calculated to calculate the first decubitus time zone, the second decubitus time zone, the awakening time zone, and the sleeping time zone, and the indexes are the constant term and the (LF / HF) of the awakening time zone. / It has a gist in that it is expressed by a variable term of activity amount and an expression including. Here, the amount of activity is an acceleration or an angular velocity accompanying the movement of the subject, and Hf1> Lf1 and Hf2> Lf2.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _m2 : There are two or more first _decubitus time zones L _m1 , and the gap time zone L _{sm where NA <C1 between two adjacent first decubitus time zones L m11} and L m12 is the first. 2 If it is within the predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled (3) Calculation of the awakening time zone and the sleeping time zone Sleeping time zone : The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time Awakening time zone: The time zone excluding the sleep time zone from the predetermined measurement unit time]

The time calculation unit of the second device can objectively and accurately calculate the awakening time zone and the sleeping time zone according to the conditions (1) to (3). In addition, the index creation unit can easily create an index for discriminating the neuropsychiatric state from the easily measurable parameters (LF / HF) and the amount of activity. In addition, it is possible to objectively determine the neuropsychiatric state by using the created index.

In the second device, it is preferable that the index includes the variable term of sleep time zone (LF / HF) × activity amount.

In the second device, the time calculation unit further calculates the third recumbent time zone, which is the total time of the first recumbent time zone and the second recumbent time zone minus the time of the sleep time zone. Yes, it is preferable that the index includes a variable term of the total time of the third recumbent time zone.

The first method and the first apparatus of the present invention can objectively and accurately calculate the sleep time zone according to the conditions (1) to (4), and are easily measurable parameters. From the time of the examiner's lying posture, an index for discriminating the neuropsychiatric state can be easily created. In addition, the second method and the second device of the present invention can objectively and accurately calculate the awakening time zone and the sleeping time zone according to the conditions (1) to (3), and can easily measure them. From the possible parameters of the awakening time zone (LF / HF) and the amount of activity, it is possible to easily create an index for discriminating the neuropsychiatric state. It is possible to objectively determine the neuropsychiatric state by using the indexes created by these methods and devices.

The explanatory diagram of the power spectrum integral which concerns on this invention is shown. The block diagram which shows the structure of the index making apparatus for discriminating the neuropsychiatric state which concerns on Embodiment 1 of this invention is shown. The block diagram which shows the structure of the index making apparatus which discriminates the neuropsychiatric state which concerns on Embodiment 2 of this invention is shown. The block diagram which shows the structure of the index making apparatus which discriminates the neuropsychiatric state which concerns on Embodiment 3 of this invention is shown.

1. 1. Method of creating an index to determine the neuropsychiatric state [First method]
The method for creating an index for discriminating the neuropsychiatric state (first method) of the present invention is a negative acceleration according to a step of measuring the acceleration TA in the height direction of the subject and the conditions (1) to (4) described later. The step of calculating NA to calculate the 1st decubitus time zone, the 2nd decubitus time zone, the sleeping time zone and the 3rd decubitus time zone, the constant term, and the total time of the 3rd decubitus time zone. It has a variable term and a step to create an index represented by an expression containing. The details of each step will be described.

(A) Measurement Step The first method includes a step of measuring the acceleration TA in the height direction of the subject.

The acceleration TA in the height direction of the subject is used for calculating the first lying time zone, the second lying time zone, the sleeping time zone, and the third lying time zone, which will be described later. The height direction is the direction from the subject's feet to the head. In general, when the accelerometer is adjusted so that the acceleration in the height direction becomes a negative value when standing, the acceleration in the height direction in the lying position tends to be larger than the acceleration in the height direction in the standing or sitting position. be. When the subject is standing and there is no movement, the magnitude of the acceleration in the height direction of the subject is 1. In the case of psychiatric disorders, the activity in the sleep time zone is excessively high, and the activity in the awake time zone is excessively low. Therefore, it can be said that the acceleration TA in the height direction of the subject indicates the state of the neuropsychiatric nerve.

In the measurement step, it is preferable to measure the acceleration TA in the height direction for a predetermined measurement time. The predetermined measurement time refers to the total time for measurement. In order to improve the accuracy of the index, it is preferable that the predetermined measurement time is a length of time during which one or more sleep time zones and one or more awakening time zones can be obtained. Alternatively, the predetermined measurement time may be a length of time during which two or more awakening time zones and two or more sleeping time zones can be obtained. Therefore, the predetermined measurement time is preferably 2 days or more, more preferably 3 days or more, and further preferably 4 days or more. Further, the longer the predetermined measurement time is, the more reliable the data can be acquired. However, in consideration of the burden on the subject, the predetermined measurement time can be set, for example, within 14 days, more preferably within 10 days. .. The sleep time zone is a time zone in which the subject is sleeping, and the awakening time zone is a time zone in which the subject is awake. The sleep time zone will be described in detail later.

(B) Time calculation step In the first method, the negative acceleration NA is calculated under the following conditions (1) to (4), and the first recumbent time zone, the second recumbent time zone, the sleep time zone, and the first method are performed. 3 It has a step of calculating the lying time zone.

Condition (1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA

Negative acceleration NA has a sign of minus and is ^{represented by a ratio to gravity acceleration g (= 9.8 m / s 2} ) (unit: dimensionless quantity). To calculate the negative acceleration, the acceleration TA in the height direction measured in the measurement step (A) is used.

Condition (2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Lying time zone L _m2 : There are two or more first lying time zones L _m1 , and the gap time zone L _sm where NA <C1 between two adjacent first lying time zones L _m11 and L _m12 If it is within the second predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled.

The recumbent time zone indicates a lying posture, for example, a supine posture in which the patient is lying on his / her back, a lateral lying posture in which he / she is lying on his / her side, and a prone position in which he / she is lying on his / her stomach. Includes time. In the first method, the decubitus time zone is divided into a first decubitus time zone and a second decubitus time zone for calculation.

The first recumbent hours L _m1 NA ≧ C1 (C1 is a constant) is the time period is a first predetermined time above T _1. A negative acceleration NA of C1 or higher indicates that the subject is in the recumbent position. In the calculation of the first recumbent time zone L _m1 , a threshold value based _{on the first predetermined time T 1} is set in order to improve the calculation accuracy of the sleep time zone.

The first predetermined time T ₁ can be set, for example, preferably 30 minutes or more, more preferably 45 minutes or more, still more preferably 1 hour or more, but set to 2 hours or less, or 1.5 hours or less. You can also. It can be said that the time zone in which NA ≥ C1 is less than the first predetermined time T ₁ is actually in the lying position, but it can be evaluated as the original sleep from the viewpoint of activity rhythm such as napping and nap for a relatively short time. In order to prevent the time that cannot be calculated as the sleep time zone, the time zone in which NA ≧ C1 is _{less than the first predetermined time T 1} is not regarded as the first recumbent time zone.

The second recumbent time zone L _m2 is a gap time zone in which there are two or more first recumbent time zones L _m1 and NA <C1 between two adjacent first recumbent time zones L _m11 and L _m12. When L _sm is within the second predetermined time T ₂ , it is the total time zone of the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm. NA <C1 indicates that the patient is in a posture other than the recumbent position. The reason why the second decubitus time zone is calculated in this way is that when the first decubitus time zone L _m1 is fragmented, the gap between the two first decubitus time zones L _m11 and L _m12. This is because it is calculated as one recumbent time zone (second recumbent time zone L _m2 ) including the time zone L _sm. This makes it easier to calculate the sleep time zone even if the subject frequently takes a posture other than the recumbent position during the sleep time zone.

The second predetermined time T ₂ can be set, for example, preferably 30 minutes or more, more preferably 45 minutes or more, still more preferably 1 hour or more, and may be set to 2 hours or less, or 1.5 hours or less. can. Further, the first predetermined time T ₁ and the second predetermined time T ₂ may be the same or different.

The value of the constant C1 (unit: dimensionless quantity) is not particularly limited, but is preferably -0.75 or more, more preferably -0.62 or more, and preferably -0.5 or more, for example. More preferred. The constant C1 is allowed to be -0.4 or less, or -0.45 or less.

Condition (3) Calculation of sleep time zone Sleep time zone: The longest time zone of _{the first decubitus time zone L m1} and the second decubitus time zone L _{m2 during a predetermined measurement unit time.}

The sleep time zone is a time zone in which the subject is sleeping, and refers to the longest time zone of _{the first recumbent time zone L m1} and the second recumbent time zone L _{m2 during a predetermined measurement unit time.} When the subject self-reports, there is room for assumptions and misunderstandings, but in the present invention, since the measured acceleration TA in the height direction of the subject is used, the calculation of the sleep time zone is objective. It can be done accurately and accurately.

The predetermined measurement unit time is a time length set for estimating the daily sleep time zone. The predetermined measurement unit time is preferably 12 hours or more, more preferably 18 hours or more, and preferably 24 hours or less. That is, the number of sleep time zones is preferably one per day. Since it is common for both day shift workers and night shift workers to be awake around 18:00, the starting point of the predetermined measurement unit time is preferably set from 17:00 to 19:00.

Condition (4) Calculation of the 3rd decubitus time zone 3rd decubitus time zone: The time zone obtained by excluding the sleep time zone from the total time of _{the 1st decubitus time zone L m1} and the 2nd decubitus time zone L _m2.

The third recumbent time zone is the time during which the patient is in the recumbent posture other than the sleep time zone. Here, the total time of the third decubitus time zone includes the first decubitus time zone constituting the second decubitus time zone. If you have an irregular life or symptoms such as chronic fatigue, depression, or insomnia, the total time in the 3rd decubitus time zone tends to be longer, so the total time in the 3rd decubitus time zone is also neuropsychiatric. It can be said that it shows the state of.

(C) Index Creation Step The first method includes a step of creating an index represented by an equation including a constant term and a variable term of the total time in the third recumbent time zone. Since it is necessary to process a large amount of data, the index is preferably created by a machine such as a computer.

The index is an expression containing a constant term and a variable term. The constant term a ₀ is an arbitrary real number, and its value is not particularly limited and may be 0. The variable term is represented by the product of the nth-order variable x _{i and the coefficient a i} , which is an arbitrary real number. The variable term of the total time of the third recumbent time zone means that the variable term _{includes the variable x i} of the nth order (the total time of the third recumbent time zone). The degree n of the variable x _i is preferably 1 or more or 2 or more in order to simplify the expression, and 5 or less or 4 or less in order to prevent the expression from becoming complicated.

In the first method, it is possible to easily create an index for discriminating the neuropsychiatric state from the acceleration TA in the height direction of the subject, which is a parameter that can be easily measured. In addition, it can be said that the third decubitus time zone indicates the state of the psychiatric nerve. By setting the variable term as described above, an index with high discrimination accuracy can be obtained. Furthermore, it is possible to objectively determine the neuropsychiatric state by using the created index. In the above index, the total time of the third recumbent time zone is preferably the total time of the third recumbent time zone per predetermined measurement unit time.

The index is preferably expressed by an equation consisting of a constant term and a variable term of the total time in the third recumbent time zone. By setting the index in this way, the accuracy of discriminating the neuropsychiatric state can be further improved.

In order to improve the accuracy of discriminating the neuropsychiatric state using the index, a variable term other than the variable term of the total time in the third recumbent time zone may be included. In the formula of the index, if the coefficient of each variable term is an arbitrary real number, its magnitude relation is not particularly specified.

The first method includes a step of measuring the beat interval of the subject and an activity amount represented by an acceleration or an angular velocity accompanying the movement of the subject, and a step of converting the beat interval into a frequency spectrum. It may have a step of calculating LF which is a value obtained by definitely integrating the power spectrum obtained in 1 to frequencies Lf1 to Lf2 and HF which is a value obtained by definitely integrating the power spectrum from frequencies Hf1 to Hf2. In that case, the indicators are (LF / HF) x activity variable term of sleep time zone, beat interval / activity variable term of sleep time zone, HF / activity variable term of sleep time zone, sleep time. It may contain at least one of the variable terms of the time of the band.

(C) In the index creation step, (LF / HF) × activity amount in the sleep time zone, pulsation interval / activity amount in the sleep time zone, and HF / activity amount in the sleep time zone are average values in the sleep time zone, respectively. It is preferable that the value is, and more preferably, it is an average value in all the measured sleep time zones. By setting the variable term in this way, it becomes easy to create an index. Similarly, in other variable terms described later, it is preferable to use the average value of the variables in the awakening time zone or the sleeping time zone.

The first method may further include a step of calculating the awakening time zone under the following condition (I). In that case, the index is at least one of the variable term of (LF / HF) / activity amount in the awakening time zone, the beat interval of the awakening time zone × the constant term of the activity amount, and the variable term of HF of the awakening time zone × the variable amount of the activity amount. Or one may be included.
[conditions:
(I) Calculation of awakening time zone Awakening time zone: Time zone obtained by excluding sleep time zone from predetermined measurement unit time]

The method of measuring the pulsation interval and the amount of activity, and the method of calculating the LF, HF, and the awakening time zone will be described in detail in the second method described later.

The first method may have the step of the time zone is NA ≧ C1 calculates the fourth recumbent hours is within a third predetermined time T _3. In that case, it is preferable that the third predetermined time T ₃ is less than the first predetermined time T ₁ and the index includes the variable term of the total time of the fourth recumbent time zone. The band 4 lying position time, while indicating the period during which taking a supine posture other than band first supine time, relatively and its time (third predetermined time) is less than T ₁ first predetermined time It can be said that it is a short time. However, if the patient is in the lying position during awakening even for a short time, for example, the sympathetic nerve and the parasympathetic nerve may be out of balance, so the total time of the 4th lying time zone is also the psychiatric nerve. It can be said that it shows the state of. Therefore, by setting the variable term as described above, an index with high discrimination accuracy can be obtained. The total time of the 4th decubitus time zone is preferably the total time of the 4th decubitus time zone per predetermined measurement unit time.

The index created by the first method may be used to determine the neuropsychiatric status of both male and female subjects, but according to the verification results described later, the first method is the mentality of female subjects. It can be preferably used for discriminating the neurological state.

[Second method]
The method for creating an index for discriminating the neuropsychiatric state (second method) of the present invention includes the beat interval of the subject, the amount of activity expressed by the acceleration or angular velocity accompanying the movement of the subject, and the test. The LF, which is a constant integral of the power spectrum obtained by including the step of measuring the acceleration TA in the height direction of the person and the step of converting the beat interval into the frequency spectrum, and the power spectrum. Negative acceleration NA is calculated according to the step of calculating HF, which is a constant integrated value from frequencies Hf1 to Hf2, and the following conditions (1) to (3), and the first decubitus time zone and the second decubitus time zone. , A step of calculating the awakening time zone and the sleeping time zone, and a step of creating an index represented by an equation including a constant term and a variable term of (LF / HF) / activity amount of the awakening time zone. Have. Here, Hf1> Lf1 and Hf2> Lf2. The details of each step will be described.

(A) Measurement step In the second method, the beat interval of the subject, the amount of activity expressed by the acceleration or angular velocity accompanying the movement of the subject, the acceleration TA in the height direction of the subject, and the acceleration TA in the height direction of the subject, Has a step to measure.

The beat interval refers to the interval between heartbeats or pulses (unit: ms). The heartbeat interval is acquired by reading the interval between R waves from the electrocardiogram or measuring the interval between adjacent heartbeats. The pulse interval is acquired by measuring the interval between adjacent pulses. The beat interval or its pulsation is said to be an index of physical and psychological stress, and reflects the balance of the neuropsychiatric states of the sympathetic and parasympathetic nerves, which are the autonomic nervous system.

As the pulsation interval, it is preferable to use the RR interval (hereinafter, referred to as “RRI”), which is the interval between the R wave and the R wave in the electrocardiographic signal. In RRI, since the peak of the signal appears clearly, erroneous recognition of the peak position is unlikely to occur, so that the accuracy of the beat interval can be improved.

The amount of activity is the acceleration or angular velocity associated with the movement of the subject. Acceleration A is a composite value of the acceleration accompanying the movement of the subject and the gravitational acceleration acting on the subject, and is ^{expressed as a ratio to the gravitational acceleration g (= 9.8 m / s 2} ) (unit: dimensionless). amount). Specifically, as expressed by the following equation (1), the acceleration A is the square of the accelerations x, y, and z in the X-axis, Y-axis, and Z-axis directions, which are the accelerations associated with the movement of the subject. It is the value obtained by subtracting (g / g) = 1 as ^{the gravitational acceleration g (= 9.8 m / s 2} ) acting on the subject from the square root of the sum (where the unit g is the magnitude of the gravitational acceleration). show). When the subject is standing and there is no movement, the values of x, y, and z are almost 0, so that the magnitude of the acceleration in the height direction of the subject is almost 1.

_{The angular velocity Ω is the square root of the sum of squares of the angular velocities ω x} , ω _y , and ω _z around the X-axis, Y-axis, and Z-axis of the subject, and the unit is rad / s or 1 / s. The angular velocity Ω is expressed by the following equation (2).

Since the angular velocity detects rotation, it is suitable for detecting, for example, the frequency of turning over of a subject during sleep time. In order to make it easier to detect the posture of the subject, the amount of activity is preferably acceleration.

In the case of psychiatric disorders, the amount of activity during sleep is excessively high, and the amount of activity during wakefulness is excessively low. Therefore, it can be said that the amount of activity also indicates the state of neuropsychiatry.

In the second method, the acceleration TA in the height direction of the subject is used for calculating the first recumbent time zone, the second recumbent time zone, the awakening time zone, and the sleep time zone, which will be described later. The acceleration TA in the height direction can be measured as described in the first method.

In the measurement step, it is preferable to measure the pulsation interval, the amount of activity, and the acceleration TA in the height direction for a predetermined measurement time. The predetermined measurement time can be set in the same manner as in the case of measuring the acceleration TA in the height direction by the first method.

(B) LF, HF calculation step In the second method, the power spectrum obtained by including the step of converting the beat interval into a frequency spectrum is LF, which is a constant integration value from frequencies Lf1 to Lf2, and the power spectrum is frequency. It has a step of calculating HF, which is a value obtained by constant integration from Hf1 to Hf2.

LF is a value obtained by constant integration of a power spectrum obtained by including a step of frequency spectrum conversion of a beat interval, which is a time signal f, from frequencies Lf1 to Lf2, and HF defines the power spectrum from frequencies Hf1 to Hf2. It is an integrated value, and Hf1> Lf1 and Hf2> Lf2. ^{For example, LF is a constant integration of the power spectrum F 2} (first power spectrum) obtained by squaring the beat interval, which is the time signal f, obtained by frequency spectrum conversion (frequency spectrum F) from frequencies Lf1 to Lf2. The HF can be a value obtained by definitely integrating the power spectrum F ² (first power spectrum) from the frequencies Hf1 (> Lf1) to Hf2 (> Lf2). The unit of LF and HF calculated using the first power spectrum F ^{2 is ms 2} . As a method of frequency spectrum transform, for example, a fast Fourier transform (FFT), a wavelet analysis, a maximum entropy method, or the like can be used. In this specification, the case where FFT is used will be described as an example, but of course, other methods can also be used.

_{In the present specification, the discrete Fourier transform G of the beat interval RRI k obtained} by spline-interlacing the beat interval and re-sampling at the sampling interval Δt is expressed by the following equation (3), and the power spectrum F ² (first power). The spectrum) (unit: ms ² / Hz) is expressed by the following equation (4). Here, k represents a time series, N represents the number of data, S is an arbitrary scale, and generally S = 1 in the power spectrum.

On the other hand, as the values of LF and HF, the method of the present invention is also obtained by definitely integrating the power spectrum F (second power spectrum) (unit: ms) obtained from the value obtained by converting the beat interval into a frequency spectrum in a predetermined interval. include. As described above, if the value obtained by converting the beat interval into the frequency spectrum is used as the power spectrum, the values of LF and HF can be calculated more easily. The units of LF and HF calculated using the second power spectrum F are dimensionless quantities. The power spectrum F (second power spectrum) is represented by the following equation (5).

A detailed calculation method of LF and HF will be described with reference to FIG. FIG. 1 is an explanatory diagram of a power spectrum integral according to the present invention. The vertical axis of FIG. 1 is the power spectral density (unit: ms ² / Hz), and the horizontal axis is the frequency (unit: Hz). LF is a value obtained by definitely integrating the power spectrum (for example, the first power spectrum F ² ) from, for example, 0.04 Hz (Lf1) to 0.15 Hz (Lf2), and is the portion hatched by diagonal lines in FIG. The area. Here, Lf1 <Lf2. On the other hand, HF is a value obtained by definitely integrating the power spectrum (for example, the first power spectrum F ² ) from, for example, 0.15 Hz (Hf1) to 0.4 Hz (Hf2), and is hatched by vertical lines in FIG. It is the area of the part. Here, Hf1 <Hf2. In FIG. 1, the integration range is set so that both Lf2 and Hf1 are equal to 0.15 Hz, but if the relationship of Lf1 <Hf1 and Lf2 <Hf2 is satisfied, Lf2 and Hf1 have the same value. May be different. Here, the method of power spectrum integration ^{has been described using the first power spectrum F 2} , but the definite integral by the second power spectrum F can also be performed in the same manner.

The power spectrum obtained by frequency spectrum conversion is divided into LF, which is a component derived from fluctuations in blood pressure and is also called a Mayer-Wave-related component, and HF, which is a component derived from respiration. The blood pressure fluctuation component LF is a power spectrum around 0.1 Hz and is related to both sympathetic nerve activity and parasympathetic nerve activity. On the other hand, the respiratory component HF has a power spectrum around 0.3 Hz and is considered to be related to parasympathetic nerve activity. From the above, the integration range of LF indicating sympathetic nerve activity and parasympathetic nerve activity preferably includes at least 0.1 Hz and Lf1 <0.1 <Lf2. Further, Lf1 is more preferably 0.03 Hz or higher, and further preferably 0.04 Hz or higher. Further, Lf1 is preferably 0.05 Hz or less, and more preferably 0.045 Hz or less. Lf2 is preferably 0.13 Hz or higher, more preferably 0.14 Hz or higher, preferably 0.16 Hz or lower, and more preferably 0.15 Hz or lower. The integration range of HF indicating parasympathetic nerve activity preferably includes at least 0.3 Hz and Hf1 <0.3 <Hf2. Hf1 is more preferably 0.14 Hz or higher, further preferably 0.15 Hz or higher, and may be 0.17 Hz or lower, or 0.16 Hz or lower. Hf2 is preferably 0.38 Hz or higher, more preferably 0.39 Hz or higher, more preferably 0.41 Hz or lower, and even more preferably 0.4 Hz or lower.

(C) Time calculation step In the second method, the negative acceleration NA is calculated under the following conditions (1) to (3), and the first decubitus time zone, the second decubitus time zone, the awakening time zone, and sleep are calculated. It has a step of calculating the time zone.

Condition (1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA

The condition (1) can be performed in the same manner as the condition (1) of the (B) time calculation step of the first method.

Condition (2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Lying time zone L _m2 : There are two or more first lying time zones L _m1 , and the gap time zone L _sm where NA <C1 between two adjacent first lying time zones L _m11 and L _m12 If it is within the second predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled.

Also in the second method, the decubitus time zone is divided into the first decubitus time zone and the second decubitus time zone for calculation. The condition (2) can be performed in the same manner as the condition (2) of the time calculation step (B) of the first method.

Condition (3) Calculation of awakening time zone and sleeping time zone Sleeping time zone: The longest time zone of _{the first decubitus time zone L m1} and the second decubitus time zone L _{m2 during a predetermined measurement unit time.} : Time zone excluding sleep time zone from predetermined measurement unit time

As described in the first method, the sleep time zone is the time zone in which the subject is sleeping, and the first recumbent time zone L _m1 and the second recumbent time zone L _{m2 during the predetermined measurement unit time.} Refers to the longest time zone. The awakening time zone refers to the time zone in which the subject is awake, that is, the time zone in which the subject is awake, and is a time zone other than the sleeping time zone. When the subject self-reports, there is room for assumptions and misunderstandings, but in the present invention, since the measured acceleration TA in the height direction of the subject is used, the awakening time zone and the sleeping time zone are used. The calculation can be performed objectively and accurately.

The predetermined measurement unit time can be set in the same manner as in the first method.

Since the second method has the conditions (1) to (3), it is possible to objectively and accurately calculate the awakening time zone and the sleeping time zone.

(D) Index Creation Step The second method includes a step of creating an index represented by an equation including a constant term and a variable term of (LF / HF) / activity amount in the awakening time zone. Since it is necessary to process a large amount of data, the index is preferably created by a machine such as a computer.

Similar to the first method, the index is an expression containing a constant term and a variable term. The constant term a ₀ is an arbitrary real number, and its value is not particularly limited and may be 0. The variable term is represented by the product of the nth-order variable x _{i and the coefficient a i} , which is an arbitrary real number. The variable term of (LF / HF) / activity amount in the awakening time zone includes _{the variable x i} of {(LF / HF)} / (activity amount in the awakening time zone) in which the variable term is nth order. It means that you are. The degree n of the variable x _i is preferably 1 or more or 2 or more in order to simplify the expression, and 5 or less or 4 or less in order to prevent the expression from becoming complicated.

In the second method of the present invention, it is possible to easily create an index for discriminating the neuropsychiatric state from the parameters (LF / HF) that can be easily measured and calculated and the amount of activity. In addition, it is possible to objectively determine the neuropsychiatric state by using the created index.

(D) In the index creation step, the (LF / HF) and the amount of activity in the awakening time zone are preferably the average value of the (LF / HF) / the amount of activity in the awakening time zone, and all the measured awakening times. More preferably, it is the average value of (LF / HF) / activity in the band. By setting the variable term in this way, it becomes easy to create an index. Similarly, in other variable terms described later, it is preferable to use the average value of the variables in the awakening time zone or the sleeping time zone.

In the second method, it is preferable that the index includes the variable term of sleep time zone (LF / HF) × activity amount. As described above, it can be said that LF, HF and the amount of activity indicate the state of the psychiatric nerve. Therefore, by setting the variable term as described above, an index with high discrimination accuracy can be obtained. Above all, the index may be expressed by an equation consisting of a constant term, a variable term of (LF / HF) / activity amount in the awakening time zone, and a variable term of (LF / HF) × activity amount in the sleeping time zone. preferable. By setting the index in this way, the accuracy of discriminating the neuropsychiatric state can be further improved.

It is preferable that the second method further includes a step of calculating a third recumbent time zone obtained by subtracting the sleep time zone from the total time of the first recumbent time zone and the second recumbent time zone. In that case, it is preferable that the index includes the variable term of the total time of the third recumbent time zone. As explained in the time calculation step (B) of the first method, it can be said that the third recumbent time zone indicates the state of the psychiatric nerve. By setting the variable term in this way, an index with high discrimination accuracy can be obtained.

In the second method, it is preferable that the index includes a variable term of pulsation interval × activity amount in the awakening time zone. As described above, since the pulsation interval and the amount of activity in the awakening time zone indicate the state of the psychiatric nerve, by setting the variable term in this way, an index with high discrimination accuracy can be obtained.

In the second method, it is preferable that the index includes the variable term of HF × activity amount in the awakening time zone. Since the HF and the amount of activity in the awakening time zone indicate the state of the psychiatric nerve, by setting the variable term in this way, an index with high discrimination accuracy can be obtained.

In the second method, it is preferable that the index includes the variable term of the beat interval / activity amount during the sleep time zone. Since the beat interval and the amount of activity in the sleep time zone indicate the state of the psychiatric nerve, setting the variable term in this way provides an index with high discrimination accuracy.

In the second method, it is preferable that the index includes at least one of the variable term of HF / activity in the sleep time zone and the variable term of the time in the sleep time zone. Since the HF and the amount of activity in the sleep time zone indicate the state of the psychiatric nerve, by setting the variable term in this way, an index with high discrimination accuracy can be obtained. Among them, the indicators are the constant term, the variable term of (LF / HF) / activity amount in the awakening time zone, the variable term of (LF / HF) × activity amount in the sleeping time zone, and the variable term of the third recumbent time zone. The variable term of total time, the variable term of beat interval x activity amount in the awakening time zone, the variable term of HF x activity amount in the awakening time zone, the variable term of the beat interval / activity amount in the sleeping time zone, It is preferably expressed by an equation consisting of a variable term of HF / activity during sleep time. By setting the index in this way, the accuracy of discriminating the neuropsychiatric state can be further improved.

The second method may further comprise the step of time zone is NA ≧ C1 calculates the fourth recumbent hours is within a third predetermined time T _3. In that case, it is preferable that the third predetermined time T ₃ is less than the first predetermined time T ₁ and the index includes the variable term of the total time of the fourth recumbent time zone. The band 4 lying position time, while indicating the period during which taking a supine posture other than band first supine time, relatively and its time (third predetermined time) is less than T ₁ first predetermined time It can be said that it is a short time. However, if the patient is in the lying position during the awakening time even for a short time, for example, the sympathetic nerve and the parasympathetic nerve may be out of balance, so the total time in the fourth lying time is also It can be said that it indicates the state of the psychiatric nerve. Therefore, by setting the variable term as described above, an index with high discrimination accuracy can be obtained.

The third predetermined time T ₃ can be set, for example, preferably 30 minutes or more, more preferably 45 minutes or more, still more preferably 1 hour or more, but set to 2 hours or less, or 1.5 hours or less. You can also.

In the second method, the indicators are a constant term, a variable term of (LF / HF) / activity amount in the awakening time zone, a variable term of (LF / HF) × activity amount in the sleeping time zone, and a third. Variable term of total time in lying time zone, variable term of beat interval x activity amount in awake time zone, variable term of HF x activity amount in awake time zone, beat interval / activity amount in sleep time zone It can be expressed by an expression consisting of the variable term of HF / activity amount in the sleep time zone, the variable term of the time in the sleep time zone, and the variable term of the total time in the 4th lying time zone. preferable. By setting the index in this way, the accuracy of discriminating the neuropsychiatric state can be further improved.

The index created by the second method may be used to determine the neuropsychiatric status of both male and female subjects, but according to the verification results described later, the second method is the mental state of male subjects. It can be preferably used for discriminating the neurological state.

The first method or the second method further includes a step of calculating VLF, which is a value obtained by definitely integrating the power spectrum obtained including the step of converting the beat interval into a frequency spectrum from frequencies Lf3 to Lf4. May be. In that case, the index may include at least one of the VLF variable term of the awakening time zone and the VLF variable term of the sleeping time zone. Here, Lf3 <Lf4. Since VLF is considered to be related to sympathetic nerve activity, the neuropsychiatric state can also be determined by setting the index in this way. Lf3 is preferably 0.0025 Hz or higher, more preferably 0.003 Hz or higher, and even more preferably 0.0033 Hz or higher. Lf3 is preferably 0.005 Hz or less, and more preferably 0.004 Hz or less. Lf4 is preferably 0.025 Hz or higher, more preferably 0.03 Hz or higher, preferably 0.06 Hz or lower, more preferably 0.05 Hz or lower, and 0.04 Hz. The following is more preferable. Further, Lf4 may be less than Lf1 or more than Lf1, but it is preferable that Lf4 = Lf1.

The first method or the second method further includes a step of calculating ULF, which is a value obtained by definitely integrating the power spectrum obtained including the step of converting the beat interval into a frequency spectrum from frequencies Lf5 to Lf6. You may be. In that case, the index may include at least one of the variable term of ULF in the wake time zone and the variable term of ULF in the sleep time zone. Here, Lf5 <Lf6. Since ULF is considered to be related to sympathetic nerve activity, the neuropsychiatric state can also be determined by setting the index in this way. Lf5 is preferably 0.0015 Hz or less, more preferably 0.001 Hz or less, and may be 0 Hz or more, but more preferably 0 Hz. Lf6 is preferably 0.0025 Hz or higher, more preferably 0.003 Hz or higher, further preferably 0.0033 Hz or higher, and preferably 0.0045 Hz or lower, preferably 0.004 Hz. The following is more preferable. Further, Lf6 may be less than Lf3 or more than Lf3, but it is preferable that Lf3 = Lf6.

For the calculation of VLF and ULF, for example, the first power spectrum F ² may be used or the second power spectrum F may be used as the power spectrum.

In the first method or the second method, the indicators are the SDNN variable term of the awakening time zone, the SDNN variable term of the sleeping time zone, the CVRR variable term of the awakening time zone, and the CVRR variable term of the sleeping time zone. , RMSSD variable term of awakening time zone, RMSSD variable term of sleeping time zone, NN50 variable term of awakening time zone, NN50 variable term of sleeping time zone, pNN50 variable term of awakening time zone, sleeping time zone It may contain at least one of the variable terms of pNN50. SDNN (Standard Deviation of the Normal to Normal Interval) is the standard deviation of the beat interval. CVRR (Coefficient of Variation of R-R intervals) is a value obtained by dividing the SDNN by the average value of the beat interval (preferably the heartbeat interval) and multiplying by 100. RMSDD (Root Mean Square of the Successive Differences) is the square root of the mean square of the difference between adjacent beat intervals. NN50 is the total number of beats with an adjacent beat interval difference of more than 50 ms, and pNN50 is the percentage of beats with an adjacent beat interval difference of more than 50 ms. Since SDNN, CVRR, RMSSD, NN50, and pNN50 all indicate the state of the psychiatric nerve, it is possible to create an index with even higher discrimination accuracy.

In the first method or the second method, at least one of the calculated beat interval, activity amount, LF, HF, and acceleration in the height direction should be regarded as an abnormal value before the index creation step. It is preferable to determine whether or not it is included, and to exclude the value determined to be an abnormal value from the target of index creation. As a result, it is possible to prevent the abnormal value from affecting the index and the discrimination result using the index.

When calculating the first decubitus time zone, the second decubitus time zone, the awakening time zone, or the sleeping time zone in the first method or the second method, for example, with an acceleration TA or a negative acceleration NA in the height direction. The acceleration represented is preferably a value after performing a morphology calculation on the acceleration-time waveform. Morphology operations are used in image processing for noise removal. Therefore, if the value after performing the morphology calculation on the acceleration-time waveform is applied to each conditional expression, the obtained acceleration is shorter than the predetermined measurement time (for example, the predetermined measurement time). Values that change within 1/150 hours are removed. Therefore, the entire outline of the acceleration-time waveform is extracted, and it becomes easy to calculate the awakening time zone and the sleep time zone. The morphology calculation may be performed on the acceleration TA in the height direction or on the negative acceleration NA.

In the first method or the second method, binarization processing may be performed on the acceleration TA in the height direction prior to the morphology calculation. In the binarization process, for example, if the negative acceleration NA is equal to or higher than the threshold value C2 (unit: dimensionless quantity), the negative acceleration NA is regarded as 0, and if the negative acceleration NA is less than C2, it is regarded as 1. The binarization process can shorten the processing time required for the morphology operation. The value of the threshold value C2 is not particularly limited, but is preferably -0.75 or more, more preferably -0.62 or more, and further preferably -0.5 or more, for example. The threshold value C2 is allowed to be -0.4 or less, or -0.45 or less. Although the example of performing the morphology calculation on the data of the negative acceleration NA has been described, the negative acceleration NA may be calculated after performing the morphology calculation on the acceleration TA in the height direction.

The morphology operation includes, for example, an expansion operation for thickening a line, a contraction operation for thinning a line, an opening process for performing an expansion operation after a contraction operation, and a closing process for performing a contraction operation after an expansion operation. When calculating the awakening time zone and the sleeping time zone using the acceleration after the morphology calculation, it is preferable that the morphology calculation is at least one of an opening process and a closing process performed in a predetermined time width. Further, it is more preferable to perform both the opening process and the closing process as the morphology operation. By combining the expansion calculation and the contraction calculation, it becomes easy to extract the entire outline of the acceleration-time waveform, so that it becomes easier to calculate the awakening time zone and the sleep time zone.

The number of times the opening process and the closing process are performed is not particularly limited, but it is preferable that the opening process and the closing process are performed once or more, and it is more preferable that the opening process and the closing process are performed twice or more each.

The time width for performing the expansion calculation and the contraction calculation may be appropriately set, but it is preferable to increase the time width during processing each time the number of processing is repeated. By gradually increasing the time width of the opening process and the closing process in this way, it is possible to prevent the data of the acceleration changed in a short time as compared with the predetermined measurement time from being removed.

In the first method or the second method, at least one of the value of the constant term and the value of the coefficient of the variable term included in the formula representing the index is calculated by multivariate analysis prior to the index creation step. Is preferable. This makes it possible to quantitatively estimate the relationship between the neuropsychiatric state as the objective variable and the total time of the third recumbent time zone or the (LF / HF) / activity amount of the awakening time zone as the explanatory variable. .. In calculating the coefficients of the constant term and the variable term, it is preferable to use data on the presence / absence, type or severity of mental disorders of a plurality of subjects. The number of subjects required for multivariate analysis is preferably the number of explanatory variables (variable terms) × 10 or more, more preferably the number of explanatory variables × 50 or more, and further preferably the number of explanatory variables × 100 or more. Is.

The neuropsychiatric state, which is the objective variable, is preferably quantified according to the presence, type, or severity of mental disorders. The quantification of the neuropsychiatric state includes a method of binarizing the presence or absence of a psychiatric disorder and a method of classifying the psychiatric disorder according to the type or severity of the psychiatric disorder based on known criteria. Known criteria include DSM (Diagnostic and Statistical Manual of Mental Disorders), which is a criterion for mental disorders, HSCL, which is a criterion for stress, and CES-, which is a criterion for depression. D, HAMD, YMRS, PANSS (Positive and Negative Syndrome Scale), which is a criterion for schizophrenia, Toho Medical Index, which is a criterion for autonomic dysfunction, and the like can be used.

For multivariate analysis, linear regression analysis such as simple regression analysis and multiple regression analysis, binary logistic regression analysis, polynomial logistic regression analysis, cumulative logistic regression analysis and other logistic regression analysis can be used. Let Y be an index indicating the psychiatric state. In the case of linear regression analysis, the index Y is represented by Y = a ₀ + a ₁ x ₁ + a ₂ x ₂ + ... + a _i x _i . It is preferable to use binary logistic regression analysis to create an index for determining the presence or absence of mental disorders. In that case, the index Y is preferably represented by Y = 1 / {1 + exp (−b)}, b = a ₀ + a ₁ x ₁ + a ₂ x ₂ + ... + a _i x _i . Note that a ₁ , a ₂ , ... a _i are coefficients of explanatory variables and are arbitrary real numbers, but a ₁ ≠ 0. x ₁ , x ₂ , ... x _i is a variable, and i is an integer of 1 or more. When a plurality of variable terms are included in the expression of the index, if the coefficient of each variable term is an arbitrary real number, the magnitude relationship is not particularly specified. The index created by binary logistic regression analysis represents the neuropsychiatric state as a binary of 0 and 1, or a number between 0 and 1. Multiple regression analysis is preferred to create indicators for determining the type or severity of psychiatric disorders.

In the first method or the second method, the indicators are depression, depression (major depressive disorder), bipolar disorder, schizophrenia, panic disorder, obsessive-compulsive disorder, autonomic imbalance, sleep disorder, etc. It is preferable to indicate the presence or absence and degree of mental disorder. By using the created index, it is possible to determine the presence or absence and degree of specific mental disorders.

2. Method for Discriminating Psychoneurological State The present invention further includes a method for discriminating a neuropsychiatric state having a step of discriminating the neuropsychiatric state using the created index. In that case, the index standard (for example, if the index is a numerical value X or more, it is determined to be depression, etc., and the index ranking standard) is a conventionally known judgment standard (for example, in the case of depression, CES-D). Etc.) are preferably the same or correlated. As the index standard, the above-mentioned judgment standard can be used as a preferable objective variable. As a result, when the neuropsychiatric state is quantified using an index, the meaning of the number becomes clear. A specific example of a method for determining a neuropsychiatric state will be described. Prepare the index created by the first method or the second method. Here, an example is shown in which the index is represented by an expression including a constant term and a variable term of the total time of the third lying time zone, but the variable term "total time of the third lying time zone" is shown. As the value of, for example, the total time of the third recumbent time zone per predetermined measurement unit time of the measured subject is substituted. By comparing the obtained index with the standard, the neuropsychiatric state can be determined. It is preferable that the determination of the psychiatric state using the index is performed by a machine such as a computer.

In the step of determining the neuropsychiatric state, a threshold value may be set for the index to determine the neuropsychiatric state. If the calculated index is above the threshold, you have (or may have) a mental illness, and if it is below the threshold, you are not (or may have) a mental illness. It is unlikely that this is the case). The value of the threshold value is not particularly limited, but for example, in the case of an index obtained by binary logistic regression analysis, the threshold value may be set to 0.4 or more, 0.5 or more, or 0.6 or less. For example, the threshold value can be set to 0.5, and less than 0.5 can be determined as a healthy state, and 0.5 or more can be determined as a depressed state. Further, a plurality of threshold values may be provided. By setting a plurality of threshold values, it becomes possible to determine the type and severity of mental disorders.

According to the index creation method of the present invention, it is possible to create and discriminate an index that can be objectively discriminated while reducing the burden on the subject. Such indicators are useful for checking and screening for the presence, type or severity of mental illness.

3. 3. Device for creating an index for discriminating the neuropsychiatric state The device for creating an index for discriminating the neuropsychiatric state (first device) of the present invention has a time calculation unit for calculating a sleep time zone and an index for discriminating the neuropsychiatric state. It has an index creation unit to create. Further, the device for creating an index (second device) for discriminating the neuropsychiatric state of the present invention includes a processing unit for calculating LF and HF, a time calculation unit for calculating an awakening time zone and a sleep time zone, and a psychiatric nerve. It has an index creation unit that creates an index for determining the state. In the present specification, the "device for creating an index for discriminating the neuropsychiatric state" may be simply referred to as a "device". Examples of the device include a personal computer capable of transmitting and receiving various data and performing various arithmetic processes, a computer such as a microcomputer, a tablet terminal, and a smartphone. The first device will be described in the first embodiment, and the second device will be described in the second and third embodiments.

(Embodiment 1)
FIG. 2 is a block diagram showing the configuration of the device 10 (10A) according to the first embodiment of the present invention. As shown in FIG. 2, the device 10A may be provided with a receiving unit 11 for receiving data transmitted from the sensor 50 capable of measuring the acceleration TA in the height direction of the subject.

The sensor 50 includes a measuring unit 51 that detects the acceleration TA in the height direction. The measurement unit 51 of the sensor 50 measures the acceleration in at least one of the X-axis, Y-axis, and Z-axis directions of the subject corresponding to the acceleration TA in the height direction and transmits it to the time calculation unit 12. The type of the sensor 50 that measures the acceleration is not particularly limited, and for example, a piezo resistor type acceleration sensor, a piezoelectric type acceleration sensor, a capacitance type acceleration sensor, or the like can be used. Since the piezoresistive acceleration sensor uses a semiconductor, it is small and easy to mass-produce. Piezoelectric accelerometers can easily detect relatively high accelerations. Capacitive accelerometers are more sensitive than piezoresistive accelerometers, have a wider range of detectable accelerations, and are less temperature dependent.

The receiving unit 11 may be provided in the time calculation unit 12 or the index creating unit 13. As a method of transmitting / receiving data from the sensor 50 to the device 10A, wireless communication may be used or wired communication may be used. In particular, when data is transmitted / received by wireless communication, it is preferable to reduce the frequency of transmission / reception by transmitting a plurality of data at once in order to improve the life of the built-in battery.

The time calculation unit 12 receives the data of the acceleration TA in the height direction of the subject transmitted from the sensor 50 via the reception unit 11. The time calculation unit 12 calculates the negative acceleration NA from the acceleration TA in the height direction of the subject under the following conditions (1) to (4), and calculates the first decubitus time zone, the second decubitus time zone, and sleep. Calculate the time zone and the third recumbent time zone. The details of the calculation method of the 1st decubitus time zone, the 2nd decubitus time zone, the sleep time zone and the 3rd decubitus time zone are as described in "1. Method of creating an index for discriminating the neuropsychiatric state". be.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _{m2 : The gap time zone L sm in} which there are two or more of the first recumbent time zones L _m1 and NA <C1 between two adjacent first recumbent time zones L _m11 and L _m12. If it is within the second predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled (3) Calculation of the sleep time zone Sleep time zone: The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time (4) Calculation of the third recumbent time zone Third recumbent time zone: The time zone obtained by excluding the time of the sleep time zone from the total time of the first recumbent time zone L _m1 and the second recumbent time zone L _m2]

The time calculation unit 12 calculates the total time of the third recumbent time zone. As a result, the index creation unit 13 can create an index including the variable term of the total time in the third recumbent time zone.

Time calculating unit 12 further, NA ≧ C1 at a time zone may calculate the fourth recumbent hours is within a third predetermined time T _3. As a result, the index creation unit 13 can create an index including the variable term of the total time in the fourth recumbent time zone.

Although not shown, the time calculation unit 12 may have a morphology calculation unit. The morphology calculation unit performs the morphology calculation in order to remove the noise of the negative acceleration-time waveform. Further, the time calculation unit 12 may have a binarization processing unit that binarizes the magnitude of the negative acceleration value with a predetermined value as a threshold value before the processing by the morphology calculation unit. The morphology operation and the binarization process can be performed by the method described above.

The index creation unit 13 creates an index for discriminating the neuropsychiatric state using the data of the third recumbent time zone obtained by the time calculation unit 12. The index is represented by an expression including a constant term and a variable term of the total time in the third recumbent time zone. In the index creation unit 13, it is preferable that the coefficients of the constant term and the variable term included in the index are created by multivariate analysis (more preferably regression analysis). The index creation unit 13 may have an average calculation unit that calculates an average value of data in each time zone.

The index created by the index creation unit 13 may include a variable term of the total time in the fourth recumbent time zone. Not limited to this, the index created by the index creating unit 13 may include at least one of the various variable terms mentioned in "1. Method for creating an index for discriminating the neuropsychiatric state".

Although not shown, the device 10A may have a function of removing abnormal values of data. That is, the device 10A may include an abnormal value detecting unit and an abnormal value removing unit. The abnormal value detection unit and the abnormal value removal unit can be preferably provided in the front stage of the time calculation unit 12 or the front stage of the index creation unit 13. As a result, the data from which the abnormal values have been removed can be transmitted to the index creation unit 13, so that an index with high accuracy can be created.

Further, the computer 1 may be provided with a discriminating unit 20 for discriminating the neuropsychiatric state based on the index created by the index creating unit 13. By comparing the created index with the standard, the discriminating unit 20 can objectively and mechanically discriminate the presence or absence and degree of mental disorder.

Although not shown, the computer 1 may be provided with a notification unit that notifies the subject or the like of the determination result of the neuropsychiatric state from the determination unit 20. The notification method is not particularly limited to voice, still image, video and the like. It is also possible to change the content of the notification according to the level of expertise of the person to be notified, such as a specialist such as a doctor or counselor, the subject or his / her family. The discrimination result may be transmitted to a notification device other than the computer 1 and the result may be notified to the subject or the like. Examples of the notification device include an external monitor, a mobile phone, a smartphone, a tablet terminal, a speaker, an earphone, and the like.

According to the first device, the index described in the first method of "1. Method for creating an index for discriminating a psychiatric state" can be created.

(Embodiment 2)
FIG. 3 shows a block diagram showing the configuration of the device 10 (10B) according to the second embodiment of the present invention. The same components as those of the device 10A of the first embodiment are designated by the same numbers, and the description thereof will be omitted. The device 10B includes a time calculation unit 12, an index creation unit 13, and a processing unit 14 (hereinafter, referred to as “first processing unit 14A”).

The sensor 50 preferably includes a measuring unit 51 that detects the pulsation interval, the amount of activity, and the acceleration TA in the height direction.

The beat interval refers to a heartbeat or pulse interval, and it is preferable to use the RR interval (RRI) as the beat interval. In this embodiment, an example in which RRI is measured as the beat interval and acceleration is measured as the amount of activity is shown.

When the heartbeat interval is used as the beat interval, it is preferable that the sensor 50 is small and lightweight, and the electrode on the back surface of the main body is attached to the skin of the chest of the subject together with the main body. Further, the sensor 50 does not need to have the electrode and the main body integrated, and clothes, underwear, a belt or the like having an electrode made of conductive fibers or a conductive sheet / film may be used as the electrode. The measuring unit 51 of the sensor 50 measures the electrocardiographic signal with the electrodes in close contact with the chest of the subject, calculates the RRI based on the electrocardiographic signal, and transmits it to the device 10B. Electrodes can be placed on the abdomen, back, and lumbar region as well as on the chest. The RRI may be calculated by the receiving unit 11 or the time calculation unit 12.

The measuring unit 51 of the sensor 50 may measure the pulse wave. The pulse wave can be measured by irradiating a person's fingertip, earlobe, or the like with near-infrared rays having a wavelength of 700 nm to 1200 nm, and measuring the amount of near-infrared rays reflected in contact or non-contact. The sensor that measures the pulse wave has an advantage that it is easy to attach to the body, and in particular, the non-contact type that measures the pulse wave has a possibility of becoming widespread because it eliminates the trouble of attaching the sensor to the body.

The measuring unit 51 of the sensor 50 measures the amount of activity represented by the acceleration or the angular velocity accompanying the movement of the subject. Above all, the amount of activity is preferably acceleration. This eliminates the need to measure the amount of activity and the acceleration TA in the height direction separately. The type of the sensor for measuring the acceleration is not particularly limited, and examples thereof include the sensor used for measuring the acceleration TA in the height direction described in the first embodiment. The type of the sensor for measuring the angular velocity is not particularly limited, and for example, a rotary type, a vibration type, a gas type, an optical fiber type, and a ring laser type angular velocity sensor can be used.

The first processing unit 14A definitely integrates the power spectrum obtained by including the step of converting the beat interval into the frequency spectrum, LF, which is a value obtained by definitely integrating the frequencies Lf1 to Lf2, and the power spectrum from the frequencies Hf1 to Hf2. Calculate the value HF. In addition, it is preferable that the first processing unit 14A also calculates LF / HF.

Although not shown, the first processing unit 14A may include a frequency spectrum conversion unit and a power spectrum integration calculation unit. The frequency spectrum conversion unit converts the RRI, which is a time signal transmitted from the reception unit 11, into a frequency spectrum by using a frequency spectrum conversion method such as FFT. Next, the power spectrum integration calculation unit calculates the power spectrum from the spectrum obtained by the frequency spectrum conversion unit and performs integration in a predetermined frequency range to obtain LF and HF. Specifically, LF and HF can be calculated as described in "1. Method for creating an index for discriminating a neuropsychiatric state". As the power spectrum, for example, the first power spectrum F ² may be used, or the second power spectrum F may be used.

The first processing unit 14A includes a VLF which is a value obtained by definitely integrating the power spectrum obtained by including the step of converting the beat interval into a frequency spectrum from frequencies Lf3 to Lf4, and a step of converting the beat interval into a frequency spectrum. At least one of ULF, which is a value obtained by constantly integrating the obtained power spectrum from frequencies Lf5 to Lf6, may be calculated. By calculating the VLF in the first processing unit 14A, the index creating unit 13 can create an index including the variable term of the VLF. Further, by calculating the ULF in the first processing unit 14A, the index creating unit 13 can create an index including the variable term of the ULF.

It is preferable that the time calculation unit 12 classifies the data transmitted from the sensor 50 into the data of the awakening time zone and the data of the sleeping time zone. The time calculation unit 12 of the device 10B calculates (LF / HF) in the awakening time zone and the value of the activity amount in the awakening time zone.

The time calculation unit 12 of the device 10B calculates the negative acceleration NA from the acceleration TA in the height direction of the subject under the following conditions (1) to (3), and calculates the first decubitus time zone and the second decubitus time. Calculate the zone, wake time zone and sleep time zone. Specifically, among the LF and HF calculated by the first processing unit 14A and the amount of activity transmitted from the receiving unit 11, the data corresponding to the calculation condition of the awakening time zone is extracted. As a result, the value of the awakening time zone (LF / HF) and the amount of activity in the awakening time zone are calculated.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _m2 : There are two or more first _decubitus time zones L _m1 , and the gap time zone L _{sm where NA <C1 between two adjacent first decubitus time zones L m11} and L m12 is the first. 2 If it is within the predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled (3) Calculation of the awakening time zone and the sleeping time zone Sleeping time zone : The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time Awakening time zone: The time zone excluding the sleep time zone from the predetermined measurement unit time]

The time calculation unit 12 may further calculate the third recumbent time zone obtained by subtracting the sleep time zone from the total time of the first recumbent time zone and the second recumbent time zone. As a result, the index creation unit 13 can create an index including the variable term of the total time in the third recumbent time zone.

The index creation unit 13 creates an index using the data obtained by the time calculation unit 12 and the first processing unit 14A. The index is expressed by an expression including a constant term and a variable term of (LF / HF) / activity amount in the awakening time zone.

In the second embodiment, an example in which the index includes the variable term of (LF / HF) / activity amount in the awakening time zone has been described, but the index created by the device 10B is not limited to this, and the index is sleep. It may include a variable term of time zone (LF / HF) × activity amount. The index may also include a variable term for the total time of the third recumbent time zone. Furthermore, the variable term of the total time in the third recumbent time zone, the variable term of the beat interval of the awakening time zone x the amount of activity, the variable term of the variable term of HF x the amount of activity in the awakening time zone, and the beat of the sleep time zone. Interval / activity variable term, sleep time zone HF / activity variable term, sleep time zone time variable term, 4th lying time zone total time variable term, awakening time zone VLF It may contain at least one of the variable term of VLF in the sleeping time zone, the variable term of ULF in the awakening time zone, and the variable term of ULF in the sleeping time zone. By setting the index in this way, the accuracy of discriminating the neuropsychiatric state can be further improved.

(Embodiment 3)
Embodiment 3 is a configuration example of an apparatus having a processing unit that further processes beat interval data into data such as SDNN. FIG. 4 shows a block diagram showing the configuration of the device 10 (10C) according to the third embodiment of the present invention. The same components as those of the devices 10A and 10B of the first and second embodiments are designated by the same numbers, and the description thereof will be omitted. The device 10C includes a time calculation unit 12, an index creation unit 13, and a processing unit 14 (hereinafter, referred to as “second processing unit 14B”).

In the second processing unit 14B, it is preferable that at least one of SDNN, CVRR, RMSSD, NN50, and pNN50 is calculated using the beat interval data. As a result, the index creation unit 13 can create an index containing at least one of SDNN, CVRR, RMSSD, NN50, and pNN50 as a variable term.

Although not shown, the apparatus 10C may be provided with both the first processing unit 14A and the second processing unit 14B. This makes it possible to create an index with high discrimination accuracy. The configurations described in the first to third embodiments may be used in combination as appropriate.

(Verification 1)
According to the first method of the present invention, the following objective variables and explanatory variables were set for the data of 20 female subjects including depressed patients and healthy subjects, and a regression equation was created by multiple regression analysis. SPSS Statistics (manufactured by IBM) was used as the statistical analysis software.
Objective variable: Degree of depression Explanatory variable: Total time in the 3rd decubitus time zone The regression equation as the obtained index is Y = a ₀ + a ₁ × (total time in the 1st decubitus time zone), a ₀ = 8, a ₁ = 0.039.

The index created for 20 female subjects was applied. Regarding the score calculated by the index, it was determined that 14 points or less was healthy, 14 points or more and less than 20 points was mild depression, 20 points or more and less than 25 points was moderate depression, and 25 points or more was severe depression. .. The detection rate of subjects judged to be depressed based on the CES-D interview results was 55%, and the detection rate of subjects determined to be healthy based on the CES-D interview results was 67%. It was found that the neuropsychiatric state can be discriminated by using the index prepared by the first method of the present invention.

(Verification 2)
According to the second method of the present invention, the following objective variables and explanatory variables were set for the data of 17 male subjects including depressed patients and healthy subjects, and a regression equation was created by multiple regression analysis. SPSS Statistics (manufactured by IBM) was used as the statistical analysis software.
Objective variable: Degree of depression Explanatory variable: Sleep time zone (LF / HF) x activity amount, activity time zone (LF / HF) / activity amount The regression equation as the obtained index is Y = a ₀ + a. ₁ x {(LF / HF) x activity during sleep time) + a ₂ x {(LF / HF) / activity during activity time}, a ₀ = 2, a ₁ = 266.175, a ₂ =- It was 0.947.

The index created was applied to 17 male subjects. Regarding the score calculated by the index, it was determined that 15 points or less was healthy, 15 points or more and less than 25 points was mild depression, 25 points or more and less than 35 points was moderate depression, and 35 points or more was severe depression. .. The detection rate of subjects judged to be depressed based on the CES-D interview results was 75%, and the detection rate of subjects determined to be healthy based on the CES-D interview results was 67%. It was found that the neuropsychiatric state can be discriminated by using the index prepared by the second method of the present invention.

1: Computers 10, 10A, 10B, 10C: Index creating device for discriminating psychiatric state 11: Receiving unit 12: Time calculation unit 13: Index creating unit 14: Processing unit 14A: First processing unit 14B: Second processing Unit 20: Discrimination unit 50: Sensor 51: Measurement unit

Claims (14)

Steps to measure the acceleration TA in the height direction of the subject,
Negative acceleration NA is calculated according to the following conditions (1) to ( 5 ), and the first decubitus time zone, the second decubitus time zone, the sleep time zone , the third decubitus time zone, and the fourth decubitus time zone are calculated. And the steps to calculate
A step of creating an index represented by an expression including a constant term, a variable term of the total time of the third recumbent time zone, and a variable term of the total time of the fourth recumbent time zone.
A method for creating an index for discriminating a psychiatric state, which is characterized by having.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _m2 : There are two or more of the first recumbent time zones L _m1 , and the gap time zone L _sm where NA <C1 between two adjacent first recumbent time zones L _m11 and L _{m12 is} If within a second predetermined time _{T 2,} the adjacent two first supine time zone _{L m11, L m12} and the gap time period _{L sm} total time period of (3) calculating sleeping hours of sleep time zone: The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time (4) Calculation of the third recumbent time zone Third recumbent time zone: the first recumbent hours L _m1 and the second supine hours L _{m 2} or we time zone excluding the sleep time period time
(5) Calculation of the 4th lying time zone
Fourth recumbent hours: NA ≧ C1 at which time zone a time zone is within a third predetermined time T _3, where the third predetermined time T ₃ is the first less than the predetermined time T _1] Steps to measure the pulsatile interval of the subject, the amount of activity expressed by the acceleration or angular velocity associated with the movement of the subject, and the acceleration TA in the height direction of the subject.
LF, which is a definite integral value of the power spectrum obtained by including the step of converting the beat interval into a frequency spectrum from frequencies Lf1 to Lf2, and HF, which is a definite integral value of the power spectrum from frequencies Hf1 to Hf2, are calculated. Steps to do and
Negative acceleration NA is calculated according to the following conditions (1) to ( 4 ), and the first decubitus time zone, the second decubitus time zone, the awakening time zone , the sleep time zone, and the fourth decubitus time zone are calculated. Steps and
A step of creating an index represented by an equation including a constant term, a variable term of (LF / HF) / activity in the awakening time zone, and a variable term of the total time in the fourth recumbent time zone. ,
A method for creating an index for discriminating a psychiatric state, which is characterized by having.
Here, Hf1> Lf1 and Hf2> Lf2.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _m2 : There are two or more of the first recumbent time zones L _m1 , and the gap time zone L _sm where NA <C1 between two adjacent first recumbent time zones L _m11 and L _{m12 is} If it is within the second predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled (3) Calculation of the awakening time zone and the sleeping time zone Sleep time zone: The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time Awakening time zone: The sleep time zone from the predetermined measurement unit time Excluded time zone
(4) Calculation of the 4th lying time zone
Fourth recumbent hours: NA ≧ C1 at which time zone a time zone is within a third predetermined time T _3, where the third predetermined time T ₃ is the first less than the predetermined time T _1] The method according to claim 2, wherein the index includes a variable term of (LF / HF) × activity amount of the sleep time zone. Further comprising the step of calculating a third supine time period excluding the first supine time zone and the second supine time zone or al the sleeping time zone time,
The method according to claim 2 or 3, wherein the index includes a variable term of the total time of the third recumbent time zone. The method according to any one of claims 2 to 4, wherein the index includes a variable term of beat interval × activity amount in the awakening time zone. The method according to any one of claims 2 to 5, wherein the index includes a variable term of HF × activity amount in the awakening time zone. The method according to any one of claims 2 to 6, wherein the index includes a variable term of beat interval / activity amount in the sleep time zone. The method according to any one of claims 2 to 7, wherein the index includes at least one of the variable term of HF / activity amount in the sleep time zone and the variable term of the time in the sleep time zone. The method according to any one of claims 2 to 8 , wherein the RR interval, which is the interval between the R wave and the R wave in the electrocardiographic signal, is used as the pulsation interval. The method according to any one of claims 1 to 9 , wherein the index indicates the presence or absence or degree of depression or autonomic imbalance. The time calculation unit that calculates the sleep time zone and
It has an index creation unit that creates an index to determine the psychiatric state, and
The time calculation unit calculates the negative acceleration NA from the acceleration TA in the height direction of the subject under the following conditions (1) to (5 ), and calculates the first decubitus time zone, the second decubitus time zone, and sleep. It calculates the time zone , the 3rd decubitus time zone, and the 4th decubitus time zone.
The index is represented by an equation including a constant term, a variable term of the total time of the third lying time zone, and a variable term of the total time of the fourth lying time zone. A device for creating an index to determine the neuropsychiatric state.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _m2 : There are two or more of the first recumbent time zones L _m1 , and the gap time zone L _sm where NA <C1 between two adjacent first recumbent time zones L _m11 and L _{m12 is} If within a second predetermined time _{T 2,} the adjacent two first supine time zone _{L m11, L m12} and the gap time period _{L sm} total time period of (3) calculating sleeping hours of sleep time zone: The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time (4) Calculation of the third recumbent time zone Third recumbent time zone: the first recumbent hours L _m1 and the second supine hours L _{m 2} or we time zone excluding the sleep time period time
(5) Calculation of the 4th lying time zone
Fourth recumbent hours: NA ≧ C1 at which time zone a time zone is within a third predetermined time T _3, where the third predetermined time T ₃ is the first less than the predetermined time T _1] LF, which is a value obtained by definitely integrating the power spectrum obtained by including the step of converting the beat interval into a frequency spectrum from frequencies Lf1 to Lf2, and HF, which is a value obtained by definitely integrating the power spectrum from frequencies Hf1 to Hf2, are calculated. Processing unit and
A time calculation unit that calculates the awakening time zone and the sleeping time zone,
It has an index creation unit that creates an index to determine the psychiatric state, and
The time calculation unit calculates the negative acceleration NA from the acceleration TA in the height direction of the subject under the following conditions (1) to (4 ), and calculates the first decubitus time zone, the second decubitus time zone, and awakening. It calculates the time zone , sleep time zone, and 4th decubitus time zone.
The index is expressed by an equation including a constant term, a variable term of (LF / HF) / activity amount in the awakening time zone, and a variable term of the total time in the fourth lying time zone. A device for creating an index for discriminating a characteristic neuropsychiatric state.
Here, the amount of activity is an acceleration or an angular velocity accompanying the movement of the subject, and Hf1> Lf1 and Hf2> Lf2.
[conditions:
(1) Calculation of negative acceleration NA When TA ≧ 0 when the subject is standing, NA = (-1) × (TA)
If TA <0 when the subject is standing, NA = TA
(2) Calculation of 1st decubitus time zone and 2nd decubitus time zone 1st decubitus time zone L _m1 : NA ≧ C1 (C1 is a constant) is the first predetermined time T ₁ or more. Position time zone L _m2 : There are two or more of the first recumbent time zones L _m1 , and the gap time zone L _sm where NA <C1 between two adjacent first recumbent time zones L _m11 and L _{m12 is} If it is within the second predetermined time T _{2, the} time zone in which the two adjacent first recumbent time zones L _m11 and L _m12 and the gap time zone L _sm are totaled (3) Calculation of the awakening time zone and the sleeping time zone Sleep time zone: The longest time zone of the first recumbent time zone L _m1 and the second recumbent time zone L _m2 during the predetermined measurement unit time Awakening time zone: The sleep time zone from the predetermined measurement unit time Excluded time zone
(4) Calculation of the 4th lying time zone
Fourth recumbent hours: NA ≧ C1 at which time zone a time zone is within a third predetermined time T _3, where the third predetermined time T ₃ is the first less than the predetermined time T _1] The device according to claim 12 , wherein the index includes a variable term of (LF / HF) × activity amount of the sleep time zone. The time calculation unit is further provided for calculating a third supine time period excluding the first supine time zone and the second supine time zone or al the sleeping time zone time,
The device according to claim 12 or 13 , wherein the index includes a variable term of the total time of the third recumbent time zone.

Images