JP6658077B2 - How to assess stumbling risk - Google Patents
How to assess stumbling risk Download PDFInfo
- Publication number
- JP6658077B2 JP6658077B2 JP2016034126A JP2016034126A JP6658077B2 JP 6658077 B2 JP6658077 B2 JP 6658077B2 JP 2016034126 A JP2016034126 A JP 2016034126A JP 2016034126 A JP2016034126 A JP 2016034126A JP 6658077 B2 JP6658077 B2 JP 6658077B2
- Authority
- JP
- Japan
- Prior art keywords
- acceleration
- toe clearance
- minimum toe
- minimum
- clearance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 206010017577 Gait disturbance Diseases 0.000 title claims description 44
- 230000001133 acceleration Effects 0.000 claims description 259
- 238000000034 method Methods 0.000 claims description 27
- 238000000611 regression analysis Methods 0.000 claims description 9
- 238000000513 principal component analysis Methods 0.000 claims description 6
- 230000037396 body weight Effects 0.000 claims description 5
- 230000002596 correlated effect Effects 0.000 claims 1
- 210000003371 toe Anatomy 0.000 description 146
- 210000001624 hip Anatomy 0.000 description 27
- 210000002683 foot Anatomy 0.000 description 18
- 239000003550 marker Substances 0.000 description 10
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 9
- 238000011156 evaluation Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 230000005484 gravity Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 125000002066 L-histidyl group Chemical group [H]N1C([H])=NC(C([H])([H])[C@](C(=O)[*])([H])N([H])[H])=C1[H] 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000035622 drinking Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 210000001872 metatarsal bone Anatomy 0.000 description 1
- 210000000452 mid-foot Anatomy 0.000 description 1
- 238000000491 multivariate analysis Methods 0.000 description 1
- 210000004197 pelvis Anatomy 0.000 description 1
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
本発明は、歩行時の最小つま先クリアランスからつまずきリスクを評価する方法とそのための装置に関する。 The present invention relates to a method for evaluating a risk of stumbling from a minimum toe clearance during walking and an apparatus therefor.
歩行において、片方の足の踵着地からその足の次の踵着地までは１歩行周期と言われ、つま先離地からその足の次の踵着地までは遊脚期と言われる。この遊脚期におけるつま先と地面との距離がつま先クリアランスであり、遊脚中期におけるつま先クリアランスの最小値が最小つま先クリアランスである。最小つま先クリアランスを測定することにより歩行中のつまずきリスクを評価することができる。 In walking, the period from the heel landing of one foot to the next heel landing of the foot is referred to as one walking cycle, and the period from toe landing to the next heel landing of the foot is referred to as a free leg period. The distance between the toe and the ground in the swing phase is the toe clearance, and the minimum value of the toe clearance in the middle swing phase is the minimum toe clearance. Measuring the minimum toe clearance can assess the risk of stumbling during walking.
従来、つま先クリアランスの測定はモーションキャプチャシステムを用いて行なわれているが、モーションキャプチャシステムは高価であり、測定方法が煩雑である。これに対し、床反力計で測定される左右方向、前後方向、鉛直方向の床反力データを説明変数とし、身長で正規化した最小つま先クリアランスを目的変数として重回帰分析することにより回帰式を得、その回帰式を用いて任意の被験者について床反力データから最小つま先クリアランスを推定し、その推定値を所定の閾値と対比することにより、あるいは推定値のばらつきの大小により、つまずきリスクを評価する方法が提案されている（特許文献１）。 Conventionally, the measurement of toe clearance has been performed using a motion capture system, but the motion capture system is expensive and the measurement method is complicated. On the other hand, the regression equation was obtained by performing multiple regression analysis using the floor reaction force data measured by the floor reaction force meter in the horizontal direction, the front-rear direction, and the vertical direction as explanatory variables, and the minimum toe clearance normalized by height as the objective variable. , And using the regression equation to estimate the minimum toe clearance from the floor reaction force data for any subject and compare the estimated value to a predetermined threshold, or the magnitude of the variance in the estimated value, to reduce the risk of stumbling. An evaluation method has been proposed (Patent Document 1).
一方、背屈力が低下すると、遊脚期につま先が垂れてつまずきが誘発されることから、腰に加速度計を装着し、一方の足の蹴り出しから他方の足のみで立脚する立脚中期までの腰部における前後方向の平均加速度を求め、その平均加速度と背屈力との関係を得ておき、その関係を用いて任意の被験者について腰部の加速度から背屈力を推定し、つまずきに対して対策を講じることが提案されている（特許文献２） On the other hand, when the dorsiflexion force decreases, the toes hang down during the swing phase and stumbling is induced, so an accelerometer is attached to the waist and from the kicking out of one leg to the middle stage of standing with only the other leg Find the average acceleration in the back and forth direction of the waist, obtain the relationship between the average acceleration and the dorsiflexion force, estimate the dorsiflexion force from the waist acceleration for any subject using that relationship, It is proposed to take measures (Patent Document 2)
特許文献１に記載の方法では、被験者の歩行の床反力データを得るために、歩行特徴の検査会場などで被験者に床反力計上を歩行してもらうことが必要となる。しかしながら、そうして得た床反力データが、被験者の日常的な歩行特徴を表しているとは限らない。そのため、この方法では、被験者の日常のつまずきリスクを評価することが難しい。 In the method described in Patent Literature 1, in order to obtain floor reaction force data of the subject's walking, it is necessary to have the subject walk the floor reaction force meter at an examination hall for a walking characteristic or the like. However, the floor reaction force data obtained in this manner does not always represent the subject's daily walking characteristics. Therefore, it is difficult for this method to evaluate a subject's daily trip risk.
特許文献２に記載の方法では、背屈力を推定するにあたり、腰に装着した加速度計で腰の加速度を測定するので、被験者の日常的な歩行における加速度を測定することができる。しかしながら、この方法は、最小つま先クリアランスを直接的に推定するものはなく、低下によりつまずきが誘発される背屈力を推定する方法であるため、つまずきリスクの評価精度が劣る。
In the method described in
これに対し、本発明の課題は、腰に装着した加速度計を用いて最小つま先クリアランスを正確に推定し、つまずきリスクを評価できるようにすることにある。 On the other hand, an object of the present invention is to accurately estimate a minimum toe clearance using an accelerometer attached to the waist and to evaluate a risk of tripping.
本発明者は、腰に装着した加速度計で測定される加速度から導出される加速度パラメータを説明変数として最小つま先クリアランスを重回帰分析により推定するにあたり、１歩行周期内の特定の期間の加速度が最小つま先クリアランスとの関連が強いのではないかとの着想をもとに、測定された加速度の１歩行周期を時間分割し、時間分割した加速度から最小つま先クリアランスと相関の高いものを加速度パラメータとして使用すると、最小つま先クリアランスの推定精度が顕著に向上することを見出し、本発明を想到した。 When estimating the minimum toe clearance by multiple regression analysis using the acceleration parameter derived from the acceleration measured by the accelerometer mounted on the waist as an explanatory variable, the acceleration in a specific period within one walking cycle is minimized. Based on the idea that there is a strong relationship with toe clearance, one walking cycle of the measured acceleration is time-divided, and a time-divided acceleration that has a high correlation with the minimum toe clearance is used as an acceleration parameter. The present inventors have found that the accuracy of estimating the minimum toe clearance is remarkably improved, and arrived at the present invention.
即ち、本発明は、最小つま先クリアランスを目的変数とし、歩行における腰部の前後方向、左右方向及び鉛直方向の加速度から導出される加速度パラメータを説明変数に含む回帰式を使用して、被験者の右足又は左足について最小つま先クリアランスを算出し、つまずきリスクを評価するつまずきリスクの評価方法であって、
加速度パラメータとして、最小つま先クリアランスを推定する足について、１歩行周期における立脚期中期の加速度、中期両足支持期の加速度、遊脚期前半の加速度、及び遊脚期後半の加速度から選ばれる加速度を使用するつまずきリスクの評価方法を提供する。
That is, the present invention, the minimum toe clearance purposes variable, the longitudinal direction of the waist in walking, using a regression equation including the explanatory variable acceleration parameters derived from the left and right direction and the vertical direction of the acceleration, the subjects' right leg Or, calculating the minimum toe clearance for the left foot, a stumbling risk evaluation method to evaluate the stumbling risk,
As the acceleration parameter, for the foot for which the minimum toe clearance is to be estimated, an acceleration selected from the acceleration in the middle of the stance phase in one walking cycle, the acceleration in the mid-leg support phase, the acceleration in the first half of the swing phase, and the acceleration in the second half of the swing phase is used. Provides a method of assessing stumbling risk.
特に本発明は、上述のつまずきリスクの評価方法において、加速度パラメータとして、１歩行周期における立脚期中期の左右方向の加速度、中期両足支持期の前後方向の加速度、遊脚期前半の前後方向の加速度及び左右方向の加速度、並びに遊脚期後半の前後方向の加速度及び左右方向の加速度を使用する方法や、加速度パラメータとして、さらに、最小つま先クリアランスを推定する足についての踵着地から次の踵着地までの１歩行周期における所定時間間隔の前後方向、左右方向及び鉛直方向の加速度の主成分分析により算出される主成分得点であって、実測された最小つま先クリアランスと相関の高い主成分得点を使用する方法を提供する。 In particular, the present invention provides the method for evaluating a stumbling risk as described above, wherein the acceleration parameter includes, as an acceleration parameter, a lateral acceleration in a middle stance phase in one walking cycle, a longitudinal acceleration in a mid-legs support phase, and a longitudinal acceleration in a first half of a swing phase. And the method of using the acceleration in the left-right direction, and the acceleration in the front-rear direction and the left-right direction in the latter half of the swing period, and as an acceleration parameter, further, from the heel landing to the next heel landing for the foot for which the minimum toe clearance is estimated. Use the principal component score calculated by the principal component analysis of the acceleration in the front-rear direction, the left-right direction, and the vertical direction at a predetermined time interval in one walking cycle of the walking cycle, and having a high correlation with the actually measured minimum toe clearance. Provide a way.
また、本発明は、歩行中の被験者の腰部の加速度から最小つま先クリアランスを推定する最小つま先クリアランスの推定装置であって、
被験者に携帯される加速度センサ、及び加速度センサで計測された加速度を用いて最小つま先クリアランスを算出し出力する演算装置を有し、
加速度センサは、歩行中の被験者の前後方向、左右方向及び鉛直方向の加速度を測定可能であり、
演算装置は、上述の回帰式を記憶する機能、
最小つま先クリアランスを推定する足について、踵着地から次の踵着地までの１歩行周期の加速度を時間分割して所定の時間間隔の加速度を取得する機能、
該所定時間間隔の加速度から、前記回帰式の加速度パラメータである立脚期中期の加速度、中期両足支持期の加速度、遊脚期前半の加速度、及び遊脚期後半の加速度から選ばれる加速度を抽出する機能、
抽出した加速度を前記回帰式で使用して最小つま先クリアランスを算出する機能、
を有する最小つま先クリアランスの推定装置を提供する。
Further, the present invention is a minimum toe clearance estimating apparatus for estimating the minimum toe clearance from the acceleration of the waist of the subject during walking,
An acceleration sensor carried by the subject, and an arithmetic unit that calculates and outputs the minimum toe clearance using the acceleration measured by the acceleration sensor,
The acceleration sensor can measure the acceleration in the front-rear direction, the left-right direction, and the vertical direction of the subject while walking,
The arithmetic unit has a function of storing the above-described regression equation,
For foot to estimate the minimum toe clearance, it functions the acceleration of 1 walking period from the heel landing until the next heel wear land by dividing time to get the acceleration of a predetermined time interval,
From the acceleration at the predetermined time interval, an acceleration selected from the acceleration parameters of the regression equation, ie, the acceleration in the middle stance phase, the acceleration in the mid-leg and foot support phase, the acceleration in the first half of the swing phase, and the acceleration in the second half of the swing phase is extracted. function,
A function of calculating the minimum toe clearance using the extracted acceleration in the regression equation,
And an apparatus for estimating a minimum toe clearance.
本発明のつまずきリスクの評価方法によれば、歩行における腰の加速度を用いて、簡便に且つ正確に最小つま先クリアランスを推定することができ、推定した最小つま先クリアランスに基づいてつまずきリスクを的確に評価することができる。 According to the stumbling risk evaluation method of the present invention, the minimum toe clearance can be easily and accurately estimated using the waist acceleration during walking, and the stumbling risk is accurately evaluated based on the estimated minimum toe clearance. can do.
また、本発明のつま先クリアランスの推定装置によれば、腰に装着する加速度センサを用いて、本発明のつまずきリスクの評価方法が使用する最小つま先クリアランスの算出方法により最小つま先クリアランスを算出し、出力することができる。 According to the toe clearance estimating apparatus of the present invention, the minimum toe clearance is calculated by the method of calculating the minimum toe clearance used by the stumbling risk evaluation method of the present invention using an acceleration sensor attached to the waist, and output. can do.
したがって、本発明は、日常生活におけるつまずきリスクを知り、また、つまずき防止策を講じる上で有用となる。 Therefore, the present invention is useful for knowing a stumbling risk in daily life and taking a stumbling prevention measure.
以下、図面を参照しつつ、本発明を詳細に説明する。 Hereinafter, the present invention will be described in detail with reference to the drawings.
＜概要＞
本発明のつまずきリスクの評価方法では、概略、最小つま先クリアランスを目的変数とし、歩行における腰部の前後方向、左右方向及び鉛直方向の加速度から導出される加速度パラメータを説明変数に含む回帰式を使用して、任意の被験者の最小つま先クリアランスを算出し、つまずきリスクを評価する。
<Overview>
In the method for evaluating a stumbling risk of the present invention, a regression equation including, as an explanatory variable, an acceleration parameter derived from the acceleration of the waist in the front-rear direction, the left-right direction, and the vertical direction during walking with the minimum toe clearance as an objective variable is used. Then, the minimum toe clearance of any subject is calculated, and the risk of stumbling is evaluated.
ここで、最小つま先クリアランスとは、遊脚中期において最もつま先と地面との距離が近づく際の距離をいう。この最小つま先クリアランスには、その物理量のみならず、身長若しくは体重により、又はその双方の値を用いて標準化したものも含まれる。最小つま先クリアランスは、遊脚中期におけるつま先の高さであるため、身体の質量の影響を受けることが考えられる。そこで最小つま先クリアランスを体重で標準化すると、回帰式による最小つま先クリアランスの推定精度が高めることができる。 Here, the minimum toe clearance refers to the distance when the distance between the toe and the ground is shortest in the middle stage of the swing. This minimum toe clearance includes not only its physical quantity, but also standardized by height or weight, or both. Since the minimum toe clearance is the height of the toes in the middle stage of swing, it is considered that the minimum toe clearance is affected by the body mass. Therefore, if the minimum toe clearance is standardized by body weight, the accuracy of estimating the minimum toe clearance by the regression equation can be improved.
＜回帰式＞
本発明で使用する回帰式としては最小つま先クリアランスを目的変数とし、加速度パラメータを説明変数に使用して重回帰分析することにより得たものを使用することが好ましい。
<Regression equation>
As the regression equation used in the present invention, it is preferable to use a regression equation obtained by performing a multiple regression analysis using the minimum toe clearance as an objective variable and the acceleration parameter as an explanatory variable.
この回帰式で使用する加速度パラメータは、複数の被験者について測定した、歩行における腰部の前後方向、左右方向及び鉛直方向の加速度から導出される。加速度の測定は、各被験者の腰部に装着した３軸加速度センサにより測定することができる。一方、最小つま先クリアランスはモーションキャプチャにより測定することができ、左右の１方向周期ごとに測定される。 The acceleration parameters used in this regression equation are derived from the acceleration of the waist in walking in the front-rear direction, the left-right direction, and the vertical direction measured for a plurality of subjects. The acceleration can be measured by a three-axis acceleration sensor attached to the waist of each subject. On the other hand, the minimum toe clearance can be measured by motion capture, and is measured every left and right one-direction cycle.
加速度パラメータには、少なくとも、（ａ）最小つま先クリアランスを推定する足について、１歩行周期における立脚期中期の加速度、中期両足支持期の加速度、遊脚期前半の前後方向の加速度、及び遊脚期後半の加速度から選ばれる加速度を含める。 The acceleration parameters include at least (a) the acceleration in the middle stance phase in one walking cycle, the acceleration in the mid-legs support phase, the acceleration in the front-rear direction in the first half of the swing phase, and the swing phase in the walking phase for the foot for which the minimum toe clearance is estimated. Include the acceleration selected from the second half acceleration.
この他、加速度パラメータとしては、（ｂ）測定された加速度を時間分割して得られる所定時間間隔の加速度を主成分分析して得られる主成分得点であって、最小つま先クリアランスとの相関性が高いものを併用することが好ましい。（ａ）及び（ｂ）の２種の加速度パラメータを使用することにより、最小つまずきリスクの推定精度を顕著に高めることができる。 In addition, as the acceleration parameter, (b) a principal component score obtained by performing principal component analysis on acceleration at predetermined time intervals obtained by time-dividing the measured acceleration, and a correlation with the minimum toe clearance is obtained. It is preferable to use a higher one in combination. By using the two types of acceleration parameters (a) and (b), the accuracy of estimating the minimum trip risk can be significantly increased.
このうち（ａ）の加速度パラメータとしては、次の（ａ1）〜（ａ6）から選択される加速度を使用することが好ましく、特に（ａ1）〜（ａ6）の全てを使用することが好ましい。
(ａ１)踵着地から次の踵着地までの１歩行周期における立脚期中期の左右方向の加速度、好ましくは、図１に示すように１歩行周期を１００％としたときの時間率（以下、単に時間率という）が２０〜３０％の左右方向の加速度。
（ａ2）１歩行周期における中期両足支持期の前後方向の加速度、好ましくは時間率５０〜６０％の前後方向の加速度。
（ａ3）１歩行周期における遊脚期前半の前後方向の加速度、好ましくは時間率７０〜７５％の前後方向の加速度。
（ａ4）１歩行周期における遊脚期前半の左右方向の加速度、好ましくは時間率７５〜８０％の左右方向の加速度。
（ａ5）１歩行周期における遊脚期後半の前後方向の加速度、好ましくは時間率９０〜９５％の前後方向の加速度。
（ａ6）１歩行周期における遊脚期後半の前後方向の加速度、好ましくは時間率９２〜１００％の左右方向の加速度。
Among them, as the acceleration parameter (a), it is preferable to use acceleration selected from the following (a1) to (a6), and it is particularly preferable to use all of (a1) to (a6).
(a1) The acceleration in the left-right direction in the middle of the stance phase in one walking cycle from the heel landing to the next heel landing, preferably, the time rate when one walking cycle is 100% as shown in FIG. (Referred to as time rate) is a lateral acceleration of 20 to 30%.
(A2) Acceleration in the front-rear direction during the mid-foot support period in one walking cycle, preferably in the front-rear direction at a time rate of 50-60%.
(A3) Front-back acceleration in the first half of the swing phase in one walking cycle, preferably front-back acceleration at a time rate of 70 to 75%.
(A4) Horizontal acceleration in the first half of the swing phase in one walking cycle, preferably horizontal acceleration at a time rate of 75 to 80%.
(A5) The longitudinal acceleration in the latter half of the swing phase in one walking cycle, preferably the longitudinal acceleration at a time rate of 90 to 95%.
(A6) Front-back acceleration in the latter half of the swing period in one walking cycle, preferably a left-right acceleration at a time rate of 92 to 100%.
これら（ａ1）〜（ａ6）の加速度は、例えば、右足の最小つま先クリアランスと加速度とを関連づける場合、加速度センサで測定された腰部の前後方向、左右方向又は鉛直方向の加速度において右足の着床に伴う鉛直方向の加速度の極値に基づいて１歩行周期を切り出し、その右足の１歩行周期の加速度を時間分割して所定時間間隔の加速度とし、所定時間間隔の加速度とその１歩行周期における最小つま先クリアランスとの相関係数を多数の被験者から求めた場合に、相関係数が好ましくは０．２５以上、より好ましくは０．４以上となる時間率の加速度である。 These accelerations (a1) to (a6) are, for example, when associating the minimum toe clearance of the right foot with the acceleration, the landing of the right foot in the front-back, left-right, or vertical acceleration of the waist measured by the acceleration sensor. One walking cycle is cut out based on the extreme value of the accompanying vertical acceleration, the acceleration of one walking cycle of the right foot is time-divided into acceleration at a predetermined time interval, and the acceleration at a predetermined time interval and the minimum toe in the one walking cycle are calculated. When the correlation coefficient with the clearance is obtained from a large number of subjects, the acceleration is a time rate at which the correlation coefficient becomes preferably 0.25 or more, more preferably 0.4 or more.
なお、（ａ1）〜（ａ6）の加速度を定めるにあたり、１歩行周期を時間分割した所定時間間隔の加速度とその１歩行周期における最小つま先クリアランスとの相関性を調べるときに、各被験者の歩行から１歩行周期を所定数切り出して各歩行周期の所定時間間隔の加速度と最小つま先クリアランスとの相関係数を算出してもよく、切り出した複数の１歩行周期の加速度から平均化した１歩行周期の加速度と平均化した最小つま先クリアランスを算出し、平均化した１歩行周期の加速度を時間分割した所定時間間隔の加速度と平均化した最小つま先クリアランスとの相関係数を算出してもよい。最小つま先クリアランスのばらつきを重視する点からは前者が好ましい。 In determining the accelerations (a1) to (a6), when examining the correlation between the acceleration at a predetermined time interval obtained by time-dividing one walking cycle and the minimum toe clearance in the one walking cycle, it is necessary to determine the acceleration of each walking cycle. A predetermined number of one walking cycle may be cut out and a correlation coefficient between the acceleration at a predetermined time interval of each walking cycle and the minimum toe clearance may be calculated. The acceleration and the averaged minimum toe clearance may be calculated, and the correlation coefficient between the averaged acceleration of one walking cycle and the acceleration at a predetermined time interval obtained by time division and the averaged minimum toe clearance may be calculated. The former is preferred from the viewpoint of emphasizing the dispersion of the minimum toe clearance.
１歩行周期の加速度データの時間分割においては、その分割の間隔が長すぎると最小つま先クリアランスの推定精度が低下し、短すぎると計算量が増大するので５ｍ秒以上１５ｍ秒以下の所定間隔で分割することが好ましい。 In the time division of the acceleration data in one walking cycle, if the interval is too long, the accuracy of estimating the minimum toe clearance decreases. If the interval is too short, the amount of calculation increases. Is preferred.
また、（ａ1）〜（ａ6）の加速度としては、それぞれの時間率の範囲の加速度のいずれかを代表値として使用してもよく、その時間率の範囲の加速度の平均を使用してもよい。 As the accelerations (a1) to (a6), any of the accelerations in the respective time rate ranges may be used as representative values, or the average of the accelerations in the time rate ranges may be used. .
（ａ）の加速度パラメータとして、（ａ1）〜（ａ6）の時間率の範囲の、左右方向の加速度、前後方向の加速度、鉛直方向の加速度を合成したものを使用してもよい。 As the acceleration parameter (a), a combination of the acceleration in the left-right direction, the acceleration in the front-rear direction, and the acceleration in the vertical direction in the range of the time rates (a1) to (a6) may be used.
（ａ）の加速度パラメータとして、（ａ1）〜（ａ6）の時間率の範囲の左右方向の加速度、前後方向の加速度、鉛直方向の加速度又はこれらを合成したの加速度を、身長、体重又は身長と体重の双方により標準化した加速度を使用してもよい。 As the acceleration parameter of (a), the acceleration in the left-right direction, the acceleration in the front-rear direction, the acceleration in the vertical direction, or the acceleration obtained by combining these in the range of the time rate of (a1) to (a6) is defined as height, weight, or height. Acceleration normalized by both body weights may be used.
一方、（ｂ）の主成分得点からなる加速度パラメータとしては、前記複数の被験者について測定した前後方向、左右方向及び鉛直方向の加速度を時間分割し、それにより得られる前後方向、左右方向及び鉛直方向の所定時間間隔の加速度を主成分分析して主成分得点を算出し、各主成分得点と、最小つま先クリアランスを推定する側の足の最小つま先クリアランスの測定値との相関性を調べ、相関係数の高い主成分得点を使用する。好ましくは、最小つま先クリアランスの測定値に対して相関係数が０．２５以上の主成分得点を１又は複数個使用する。 On the other hand, as the acceleration parameter consisting of the principal component score of (b), the acceleration in the front-rear direction, the left-right direction, and the vertical direction measured for the plurality of subjects are time-divided, and the front-rear direction, the left-right direction, and the vertical direction obtained thereby are obtained. Principal component analysis is performed on the acceleration at a predetermined time interval to calculate the principal component score, and the correlation between each principal component score and the measured value of the minimum toe clearance of the foot on which the minimum toe clearance is to be estimated is examined, and the phase relationship is determined. Use the high number of principal component scores. Preferably, one or more principal component scores having a correlation coefficient of 0.25 or more with respect to the measured value of the minimum toe clearance are used.
（ｂ）の主成分得点からなる加速度パラメータを得るにあたり、所定時間間隔の加速度は、（ａ）の加速度パラメータを得る場合と同様に１歩行周期の加速度データを切り出し、それを時間分割することにより得ることができる。 In obtaining the acceleration parameter consisting of the principal component score of (b), the acceleration at a predetermined time interval is obtained by cutting out the acceleration data of one walking cycle and dividing it by time similarly to the case of obtaining the acceleration parameter of (a). Obtainable.
本発明においては、最小つま先クリアランスを推定する回帰式の説明変数として、上述の（ａ）の加速度パラメータ、好ましくは（ａ）及び（ｂ）の加速度パラメータに加えて、被験者の年齢、性別、身長、体重、ＢＭＩ、等から選ばれる属性パラメータを使用することがより好ましい。言い換えると、本発明で使用する回帰式は、上述の（ａ）及び（ｂ）の加速度パラメータと、必要に応じて属性パラメータを説明変数とし、最小つま先クリアランスを目的変数として重回帰分析することにより得たものが好ましい。重回帰分析では、説明変数とする（ａ）及び（ｂ）の加速度パラメータ、並びに属性パラメータをステップワイズ法により選択することが好ましい。 In the present invention, as the explanatory variables of the regression equation for estimating the minimum toe clearance, in addition to the acceleration parameter of (a) described above, preferably the acceleration parameters of (a) and (b), the age, sex, and height of the subject It is more preferable to use an attribute parameter selected from the group consisting of, weight, BMI, and the like. In other words, the regression equation used in the present invention is obtained by performing multiple regression analysis using the above acceleration parameters (a) and (b) and, if necessary, attribute parameters as explanatory variables and minimum toe clearance as an objective variable. Those obtained are preferred. In the multiple regression analysis, it is preferable to select the acceleration parameters (a) and (b) as the explanatory variables and the attribute parameters by a stepwise method.
＜回帰式を用いた最小つま先クリアランスの算出＞
任意の被験者の最小つま先クリアランスは、その被験者の歩行における腰部の前後方向、左右方向及び鉛直方向の加速度を測定し、回帰式の説明変数である加速度パラメータを得、その加速度パラメータを上述の回帰式で用いることにより算出する。こうして算出した最小つま先クリアランスを被験者の最小つま先クリアランスの推定値とする。
<Calculation of minimum toe clearance using regression equation>
The minimum toe clearance of any subject is determined by measuring the acceleration of the subject's waist in the front-rear direction, left-right direction, and vertical direction during walking, and obtaining an acceleration parameter that is an explanatory variable of the regression equation. It is calculated by using. The calculated minimum toe clearance is used as an estimated value of the minimum toe clearance of the subject.
より具体的には、例えば右足の最小つま先クリアランスを算出する場合、測定された腰部の加速度から１歩行周期を、右足の着床に伴う鉛直方向の加速度の極値により切り出し、１歩行周期毎にそれを時間分割して所定時間間隔の加速度を求め、所定時間間隔の加速度から少なくとも（ａ）の加速度パラメータ、好ましくは（ａ）及び（ｂ）の加速度パラメータを求め、１歩行周期毎の最小つま先クリアランスを算出することが好ましい。これにより最小つま先クリアランスのばらつきをみることができる。 More specifically, for example, when calculating the minimum toe clearance of the right foot, one walking cycle is cut out from the measured waist acceleration by the extreme value of the vertical acceleration accompanying the landing of the right foot, and every one walking cycle It is time-divided to determine acceleration at predetermined time intervals, and at least acceleration parameters (a), preferably (a) and (b), are determined from accelerations at predetermined time intervals. Preferably, the clearance is calculated. As a result, variations in the minimum toe clearance can be seen.
一方、（ａ）及び（ｂ）の加速度パラメータを求めるにあたり、切り出した複数の１歩行周期の加速度から平均化した１歩行周期の加速度を求め、その平均化した１歩行周囲の加速度を時間分割して所定時間間隔の加速度を求め、所定時間間隔の加速度から（ａ）及び（ｂ）の加速度パラメータを求めてもよい。こうして得た加速度パラメータを用いて回帰式から算出した最小つま先クリアランスは、複数の歩行周期における最小つま先クリアランスの平均といえる。最小つま先クリアランスのばらつきを重視する点からは、切り出した１歩行周期ごとに最小つま先クリアランスを求めることが好ましい。 On the other hand, in obtaining the acceleration parameters of (a) and (b), an acceleration of one walking cycle averaged from a plurality of cut-out accelerations of one walking cycle is obtained, and the averaged acceleration of one walking cycle is time-divided. Thus, the acceleration at predetermined time intervals may be obtained from the acceleration, and the acceleration parameters (a) and (b) may be obtained from the acceleration at predetermined time intervals. The minimum toe clearance calculated from the regression equation using the acceleration parameters thus obtained can be said to be an average of the minimum toe clearances in a plurality of walking cycles. From the viewpoint that emphasis is placed on the variation of the minimum toe clearance, it is preferable to obtain the minimum toe clearance for each cutout walking cycle.
こうして得た最小つま先クリアランスの推定値は、それを算出する回帰式の説明変数として少なくとも（ａ）の加速度パラメータが使用されているので最小つま先クリアランスの推定精度が高く、特に（ａ）及び（ｂ）の２種の加速度パラメータを使用した場合には、（ａ）、（ｂ）いずれかの加速度パラメータが説明変数である場合に比して、最小つま先クリアランスの推定精度が顕著に向上し、最小つま先クリアランスの推定値と実測値との相関の程度を表す決定係数を０．４以上とすることができる。
したがって、本発明によれば、最小つま先クリアランスの推定値に基づいてつまずきリスクを評価することができる。
The estimated value of the minimum toe clearance obtained in this way has high estimation accuracy of the minimum toe clearance because at least the acceleration parameter (a) is used as an explanatory variable of the regression equation for calculating the minimum toe clearance, and in particular, (a) and (b) When the two types of acceleration parameters are used, the accuracy of estimating the minimum toe clearance is remarkably improved as compared with the case where one of the acceleration parameters is an explanatory variable. The determination coefficient indicating the degree of correlation between the estimated value of the toe clearance and the actually measured value can be 0.4 or more.
Therefore, according to the present invention, the stumbling risk can be evaluated based on the estimated value of the minimum toe clearance.
＜つまずきリスクの評価＞
最小つま先クリアランスに基づいてつまずきリスクを評価する方法の具体的な態様としては、例えば、最小つま先クリアランスとつまずき易さとの関係を予め調べ、つまずき易さを段階的に評価する場合の最小つま先クリアランスの閾値を設定し、その閾値に対して任意の被験者の最小つま先クリアランスの推定値を比較することによりその被験者のつまずき易さを段階的に評価することができる。
<Evaluation of stumbling risk>
As a specific embodiment of the method of evaluating the stumbling risk based on the minimum toe clearance, for example, a relationship between the minimum toe clearance and the ease of stumbling is checked in advance, and the minimum toe clearance in the case of evaluating the ease of stumbling stepwise is evaluated. By setting a threshold value and comparing the estimated value of the minimum toe clearance of any subject to the threshold value, the subject's ease of tripping can be evaluated stepwise.
所定の時間帯ごとに、最小つま先クリアランスの平均を算出し、時間帯毎につまずきリスクを対比できるようにしてもよい。最小つま先クリアランスの平均として、当該時間帯の歩行の１歩行周期ごとの最小つま先クリアランスの積算値を当該時間帯の歩行のストライド数（即ち、歩数）で割ったものを使用することができる。 The average of the minimum toe clearance may be calculated for each predetermined time period, and the risk of stumbling may be compared for each time period. As the average of the minimum toe clearance, a value obtained by dividing the integrated value of the minimum toe clearance for each walking cycle of the walking in the time period by the number of stride of walking in the time period (that is, the number of steps) can be used.
また、最小つま先クリアランスのばらつきが大きいとつまずき易いことが知られていることから、所定の時間帯ごとに最小つま先クリアランスのばらつきを算出し、そのばらつきに基づいてつまずきリスクを評価してもよい。このばらつきとしては、標準偏差、分散などを使用することができる。 Further, since it is known that if the variation in the minimum toe clearance is large, it is easy to trip, the variation in the minimum toe clearance may be calculated for each predetermined time period, and the risk of tripping may be evaluated based on the variation. As the variation, a standard deviation, a variance, or the like can be used.
所定の時間帯ごとに最小つま先クリアランス又はそのばらつきがわかることにより、例えば、被験者に歩き方のアドバイスをする者が、被験者に１日のうちの時間帯に応じたアドバイスを提供することが可能となる。 By knowing the minimum toe clearance or its variation for each predetermined time period, for example, a person who advises the subject on how to walk can provide the subject with advice according to the time period of the day. Become.
＜最小つま先クリアランスの推定装置＞
任意の被験者について、腰部の加速度の測定から最小つま先クリアランスを推定する最小つま先クリアランスの推定装置としては、日常の歩行における最小つま先クリアランスを推定できるように、日常生活において携帯可能な装置であることが好ましい。
<Minimum toe clearance estimation device>
For any subject, the minimum toe clearance estimating device for estimating the minimum toe clearance from the measurement of the acceleration of the waist may be a device that is portable in daily life so that the minimum toe clearance in daily walking can be estimated. preferable.
そこで、本発明の最小つま先クリアランスの推定装置は、被験者に携帯される加速度センサ、及び加速度センサで計測された加速度を用いて最小つま先クリアランスを算出し出力する演算装置を備え、さらに必要に応じて算出したつま先クリアランスを表示するディスプレイを備える。なお、演算装置が算出した最小つま先クリアランスは、携帯電話、プリンタ等に出力できるようにしてもよい。 Therefore, the minimum toe clearance estimating device of the present invention includes an acceleration sensor carried by the subject, and a calculation device that calculates and outputs the minimum toe clearance using acceleration measured by the acceleration sensor, and further includes an as needed. A display for displaying the calculated toe clearance is provided. Note that the minimum toe clearance calculated by the arithmetic device may be output to a mobile phone, a printer, or the like.
この装置において加速度センサとしては、歩行中の被験者の前後方向、左右方向及び鉛直方向の加速度を測定可能なものを使用する。このような加速度センサとしては、ＸＹＺの３軸の加速度波形を抽出することのできる携帯端末をあげることができる。 In this device, an acceleration sensor that can measure the acceleration of the walking subject in the front-back direction, the left-right direction, and the vertical direction is used. As such an acceleration sensor, a portable terminal capable of extracting acceleration waveforms of three axes of XYZ can be cited.
演算装置は、最小つま先クリアランスを算出する前述の回帰式を記憶する機能、
左右の足のうち最小つま先クリアランスを推定する足について、踵着地から次の踵着地までの１歩行周期の加速度を時間分割して所定の時間間隔の加速度を取得する機能、
該所定時間間隔の加速度から、前記回帰式の加速度パラメータである所定の時間率の範囲
の加速度を抽出する機能、
抽出した加速度と算出した主成分得点を前記回帰式で使用して最小つま先クリアランスを算出し、最小つま先クリアランスの推定値として出力する機能を有する。
The arithmetic unit has a function of storing the above-described regression equation for calculating the minimum toe clearance,
The legs of estimating the minimum toe clearance of the left and right feet, the ability acceleration of 1 walking period from the heel landing to the next heel deposition locations by dividing the time to obtain the acceleration of the predetermined time interval,
A function of extracting, from the acceleration at the predetermined time interval, an acceleration in a predetermined time rate range, which is an acceleration parameter of the regression equation;
It has a function of calculating a minimum toe clearance by using the extracted acceleration and the calculated principal component score in the regression equation, and outputting it as an estimated value of the minimum toe clearance.
これらの機能に加えて、演算装置は、所定の時間率の願意の所定方向の加速度を抽出する機能と、所定時間間隔の加速度を用いて回帰式の加速度パラメータである主成分得点を算出する機能も有することが好ましい。 In addition to these functions, the arithmetic device has a function of extracting acceleration in a desired direction at a predetermined time rate and a function of calculating a principal component score, which is a regression-type acceleration parameter, using acceleration at predetermined time intervals. It is also preferred to have
このうち、所定の時間間隔の加速度を取得する機能の具体的態様としては、加速度センサで測定された前後方向、左右方向又は鉛直方向のそれぞれの加速度から１歩行周期の加速度を複数切り出し、切り出した各歩行周期ごとに加速度を時間分割して所定時間間隔の加速度を取得するものが好ましい。また、複数の歩行周期の加速度から平均化した１歩行周期の加速度を算出し、平均化した１歩行周期の加速度から所定時間間隔の加速度を取得できるようにしてもよい。 Among these, as a specific mode of the function of acquiring the acceleration at a predetermined time interval, a plurality of accelerations in one walking cycle are cut out and cut out from the respective accelerations in the front-rear direction, the left-right direction, or the vertical direction measured by the acceleration sensor. It is preferable that the acceleration is time-divided for each walking cycle to obtain acceleration at predetermined time intervals. Further, the acceleration of one walking cycle averaged from the accelerations of a plurality of walking cycles may be calculated, and the acceleration at a predetermined time interval may be obtained from the averaged acceleration of one walking cycle.
また、前後方向、左右方向及び鉛直方向という成分ごとの加速度ではなく、３方向を一括してベクトルとし、そのベクトルの大きさを加速度パラメータとして使用できるようにしてもよい Instead of the acceleration for each component of the front-rear direction, the left-right direction, and the vertical direction, three directions may be collectively set as a vector, and the magnitude of the vector may be used as an acceleration parameter.
演算装置に、所定の時間帯における歩行の最小つま先クリアランスの平均やばらつきを算出する機能をもうけ、それらが出力されるようにしてもよい。 The arithmetic unit may be provided with a function of calculating the average and the variation of the minimum toe clearance of walking in a predetermined time zone, and these may be output.
また、演算装置に、回帰式を用いて算出した最小つま先クリアランス、最小つま先クリアランスの平均、又は最小つま先クリアランスのばらつきに基づいてつまずきリスクを評価する機能をもたせ、評価結果が出力されるようにしてもよい。即ち、この機能により、前述のつまずきリスクの評価方法が実施されるようにしてもよい。 Further, the arithmetic unit has a function of evaluating a stumbling risk based on the minimum toe clearance calculated using the regression equation, the average of the minimum toe clearance, or the variation of the minimum toe clearance, so that the evaluation result is output. Is also good. That is, the above-described method of evaluating the stumbling risk may be performed by this function.
以下、実施例に基づいて本発明を具体的に説明する。
実施例１
（１）回帰式の取得
歩行時に痛みがない自力で歩行可能な２０歳から７３歳の男女１３３名（男性７４名、女性５９名）を被験者とした。被験者には、前日の過度な運動と飲酒を控えるように指示した。
Hereinafter, the present invention will be specifically described based on examples.
Example 1
(1) Acquisition of
被験者の腰の前後方向、左右方向、及び鉛直方向の加速度を測定するために、被験者の腰に反射マーカを取り付け、各時間帯の歩行におけるマーカの位置をモーションキャプチャシステム（Vicon社製のVicon MXシステム、Vicon Nexus）で計測することにより、腰の位置情報をもとに加速度を算出した。 A reflection marker is attached to the subject's waist to measure the acceleration of the subject's waist in the front-rear direction, the left-right direction, and the vertical direction, and the position of the marker during walking in each time zone is measured by a motion capture system (Vicon MX manufactured by Vicon Corporation). Acceleration was calculated based on waist position information by measuring with the system, Vicon Nexus).
回帰式を得るモデル構築のための前処理は以下の手順で行った。
1）計測された全被験者のマーカ座標の生データに４次のButterworthローパスフィルタをかけノイズを除去した。カットオフ周波数は１０Ｈｚとした。
2）腰（仙骨部）の３軸の加速度を抽出した。加速度は、マーカ座標を２回微分することで取得した。
Pre-processing for constructing a model for obtaining a regression equation was performed in the following procedure.
1) A 4th-order Butterworth low-pass filter was applied to the measured raw data of the marker coordinates of all subjects to remove noise. The cutoff frequency was 10 Hz.
2) Three-axis acceleration of the hip (sacrum) was extracted. The acceleration was obtained by differentiating the marker coordinates twice.
測定した３軸方向の加速度から１歩行周期を右足の着床に伴う鉛直方向の加速度の極値により切り出し、その１歩行周期の加速度を１０１等分（約１０ｍ秒間隔）に分割することで合計３０３(時間率０〜１００％）の加速度成分を得た。３軸方向の加速度成分について多変量解析手法の一つである分散共分散行列を用いた主成分分析を行うことで、合計２６の加速度主成分を抽出し、加速度主成分ごとに加速度の数値と主成分得点係数を掛け合わせ、それらの和をとることにより主成分得点を算出した。
One walking cycle is cut out from the measured three-axis direction acceleration by the extreme value of the vertical acceleration associated with the landing of the right foot, and the acceleration of one walking cycle is divided into 101 equal parts (at intervals of about 10 msec) to obtain a total. An acceleration component of 303 (
また、同じ被験者の右足の第１中足骨に反射マーカを取り付け、各時間帯の歩行におけるマーカの位置をモーションキャプチャシステム（Vicon社製のVicon MXシステム、Vicon Nexus）で計測することにより、１歩行周期の最小つま先クリアランスを測定した。 In addition, by attaching a reflective marker to the first metatarsal of the right leg of the same subject and measuring the position of the marker during walking in each time zone with a motion capture system (Vicon MX system, Vicon Nexus manufactured by Vicon), 1 The minimum toe clearance of the walking cycle was measured.
各加速度主成分の主成分得点と、最小つま先クリアランスとの相関性を調べ、相関性の高い第２１主成分と第２４主成分を抽出した。図２に、第２１主成分得点と最小つま先クリアランスとの関係を示し、図３に第２４主成分得点と最小つま先クリアランスとの関係を示す。 The correlation between the principal component score of each acceleration principal component and the minimum toe clearance was examined, and the 21st and 24th principal components having high correlation were extracted. FIG. 2 shows the relationship between the 21st principal component score and the minimum toe clearance, and FIG. 3 shows the relationship between the 24th principal component score and the minimum toe clearance.
一方、３０３の加速度成分のそれぞれについて、最小つま先クリアランスの測定値との相関性を調べ、相関係数が０．２５以上となる時間率の加速度成分を選択した。 On the other hand, for each of the acceleration components 303, the correlation with the measured value of the minimum toe clearance was examined, and the acceleration component of the time rate at which the correlation coefficient was 0.25 or more was selected.
また、各被験者の性別に、平均０、分散１となる標準化した数値を与えた。 In addition, standardized numerical values of 0 mean and 1 variance were given for each subject's gender.
最小つま先クリアランスを目的変数とし、各被験者の性別、上述の相関係数が０．２５以上の加速度成分、第２１主成分得点及び第２４主成分得点を説明変数として重回帰分析（ステップワイズ法）を行い、次の回帰式を得た。式中、Ｘt%は左右方向の時間率ｔ％の加速度を表し、Ｙt%は前後方向の時間率ｔ％の加速度を表す。 Multiple regression analysis (stepwise method) using the minimum toe clearance as the objective variable, the gender of each subject, the acceleration component having the above-mentioned correlation coefficient of 0.25 or more, the 21st principal component score, and the 24th principal component score as explanatory variables And the following regression equation was obtained. In the equation, Xt% represents acceleration at a time rate t% in the left-right direction, and Yt% represents acceleration at a time rate t% in the front-rear direction.
最小つま先クリアランス（ｃｍ）
＝4.464＋(性別)×0.367＋(Ｘ24％)×(-0.018)＋(Ｙ53％)×(-0.018)＋(Ｙ71％)×(-0.017)＋(Ｘ77％)×(-0.042)＋(Ｙ94％)×(0.007)＋(Ｘ95％)×(-0.019)＋(第２１主成分得点)×(0.056)＋（第２４主成分得点）×(-0.054)
Minimum toe clearance (cm)
= 4.464 + (sex) x 0.367 + (X24%) x (-0.018) + (Y53%) x (-0.018) + (Y71%) x (-0.017) + (X77%) x (-0.042) + ( Y94%) x (0.007) + (X95%) x (-0.019) + (21st principal component score) x (0.056) + (24th principal component score) x (-0.054)
この回帰式より算出された最小つま先クリアランスの推定値と最小つま先クリアランスの測定値との関係を図４に示す。これらの相関の程度を示す決定係数Ｒ2は０．４２１であり、この回帰式による最小つま先クリアランスの推定精度が高いことがわかる。 FIG. 4 shows the relationship between the estimated value of the minimum toe clearance calculated from the regression equation and the measured value of the minimum toe clearance. Determination coefficient R 2 indicating the degree of correlation is 0.421, it is understood that the estimation accuracy of the minimum toe clearance by the regression equation is high.
（２）時間帯毎の最小つま先クリアランスの推定
回帰式の取得に関わらなかった自力で歩行が可能な２８歳の男性を被験者とし、この被験者の日常生活での０時〜３時、３時〜６時、６時〜９時、９時〜１２時、１２時〜１５時、１５時〜１８時、１８時〜２１時、２１時〜２４時の時間帯ごとに加速度の計測が可能であったストライド数の歩行について、Ｘ（左右方向）、Ｙ（前後方向）、Ｚ（鉛直方向）の加速度を測定した。なお、０時〜３時は、０時以降３時前までを意味する。他の時間帯についても同様である。
(2) Estimation of minimum toe clearance for each time zone A 28-year-old man who can walk on his / her own irrespective of acquisition of the regression equation was set as a subject, and the subject's daily life was from 0:00 to 3:00, 3:00 to 3:00 Acceleration can be measured every time period from 6:00, 6:00 to 9:00, 9:00 to 12:00, 12:00 to 15:00, 15:00 to 18:00, 18:00 to 21:00, and 21:00 to 24:00. With respect to walking with the number of strides, accelerations in X (left and right directions), Y (front and rear directions), and Z (vertical direction) were measured. Note that 0:00 to 3:00 means from 0:00 to 3:00. The same applies to other time zones.
時間帯ごとに、測定した加速度から１歩行周期を切り出し、各歩行周期を１０１等分して時間率０〜１００％の加速度成分を求め、その加速度成分を使用して上述の回帰式から最小つま先クリアランス（ＭＴＣ）を算出し、時間帯ごとのストライド数における積算値を求め、積算値を時間帯ごとのストライド数で割ることにより最小つま先クリアランス（ＭＴＣ）の平均値を求めた。また、時間帯ごとの最小つま先クリアランスのばらつきをみるため、最小つま先クリアランス（ＭＴＣ）の標準偏差を求めた。結果を表１に示す。 For each time zone, one walking cycle is cut out from the measured acceleration, each walking cycle is divided into 101 equal parts, an acceleration component of a time rate of 0 to 100% is obtained, and the minimum toe is obtained from the regression equation using the acceleration component. The clearance (MTC) was calculated, the integrated value of the number of strides for each time zone was obtained, and the average value of the minimum toe clearance (MTC) was obtained by dividing the integrated value by the number of strides for each time zone. In addition, the standard deviation of the minimum toe clearance (MTC) was obtained in order to see the variation of the minimum toe clearance for each time zone. Table 1 shows the results.
表１の結果から、時間帯ごとに最小つま先クリアランスの推定値が異なり、つまずきリスクが異なることから、時間帯毎につまずき予防のためのアドバイスを変えることの有用性がわかる。例えば、今回加速度計測を行った被験者は６時から９時の間と、１８時から２１時の間で最小つま先クリアランスの推定値が日中よりも高い値となった。これは、出勤、帰宅時において通常よりも速歩で歩いているためであると考えられる。一方で、帰宅時である１８時から２１時の間では最も最小つま先クリアランスの標準偏差が高い。これは帰宅時に最もつまずきリスクが増加していることを意味しており、この時間帯でのつまずきリスクを予防する施策を提示することが重要である。 From the results in Table 1, the estimated value of the minimum toe clearance differs for each time zone, and the stumbling risk differs, so that the usefulness of changing advice for stumbling prevention for each time zone can be understood. For example, for the subject who performed the acceleration measurement this time, the estimated value of the minimum toe clearance became higher than that during the day between 6:00 and 9:00 and between 18:00 and 21:00. This is considered to be because the user walks faster than usual when going to work or returning home. On the other hand, the standard deviation of the smallest toe clearance is high between 18:00 and 21:00 when returning home. This means that the risk of stumbling increases most when returning home, and it is important to present measures to prevent the risk of stumbling during this time.
（３）個別の被験者における最小つま先クリアランスの推定精度
（２）の被験者の最小つま先クリアランスを（１）と同様にして１２時から１５時に測定した。その結果、（２）の推定値は３．７０ｃｍであったのに対し、日常生活における最小つま先クリアランスの実測値は３．７３ｃｍであり、推定精度の高いことが確認できた。
(3) Estimation accuracy of minimum toe clearance in individual subject The minimum toe clearance of the subject in (2) was measured from 12:00 to 15:00 in the same manner as in (1). As a result, the estimated value of (2) was 3.70 cm, whereas the measured value of the minimum toe clearance in daily life was 3.73 cm, confirming that the estimation accuracy was high.
比較例１
（１）の重回帰分析において、説明変数として、主成分得点を使用せず、各被験者の性別及び相関係数が０．２５以上の加速度成分だけを使用する以外は実施例１と同様にして次の回帰式を求めた。式中、Ｘt%は左右方向の時間率ｔ％の加速度を表し、Ｙt%は前後方向の時間率ｔ％の加速度を表し、Ｚt%は鉛直方向の時間率ｔ％の加速度を表す。
Comparative Example 1
In the multiple regression analysis of (1), the same as in Example 1 except that the principal component score is not used as an explanatory variable and only the acceleration component having a gender and a correlation coefficient of 0.25 or more for each subject is used. The following regression equation was determined. In the formula, Xt% represents the acceleration at the time rate t% in the left-right direction, Yt% represents the acceleration at the time rate t% in the front-rear direction, and Zt% represents the acceleration at the time rate t% in the vertical direction.
最小つま先クリアランス（ｃｍ）
＝4.617＋(性別)×0.331＋(Ｚ58％)×(0.009)＋(Ｙ64％)×(-0.067)＋(Ｙ65％)×(0.173)＋(Ｙ66％)×(-0.104)＋(Ｘ84％)×(-0.021)＋(Ｘ85％)×(0.04)＋(Ｙ92％)×(-0.092)＋(Ｙ93％)×(0.205)＋(Ｙ94％)×(-0.1)＋(Ｘ95％)×(-0.29)＋(Ｘ96％)×(0.521)＋(Ｘ97％)×(-0.272)
Minimum toe clearance (cm)
= 4.617 + (sex) x 0.331 + (Z58%) x (0.009) + (Y64%) x (-0.067) + (Y65%) x (0.173) + (Y66%) x (-0.104) + (X84% ) X (-0.021) + (X85%) x (0.04) + (Y92%) x (-0.092) + (Y93%) x (0.205) + (Y94%) x (-0.1) + (X95%) x (-0.29) + (X96%) x (0.521) + (X97%) x (-0.272)
この回帰式より算出された最小つま先クリアランスの推定値と最小つま先クリアランスの測定値との関係を図５に示す。これらの相関の程度を表す決定係数Ｒ2は０．３３６であり、実施例１に比して比較例１の回帰式による最小つま先クリアランスの推定精度が低いことがわかる。 FIG. 5 shows the relationship between the estimated value of the minimum toe clearance and the measured value of the minimum toe clearance calculated from the regression equation. Determination coefficient R 2 indicating the degree of correlation is 0.336, it can be seen that low estimation accuracy of the minimum toe clearance by regression of Comparative Example 1 than in Example 1.
実施例２
（１）回帰式の取得
歩行時に痛みがない自力で歩行可能な２０歳から７７歳の健常成人１２０名分（男性６３名、女性５７名）を被験者とした。被験者らは歩行動作に影響を与えるような傷害歴や疾病歴などを負っていない健常成人であり、計測時には痛みなどを訴えることなく歩行補助具や装具などを用いずに２足で歩行できる者であった。
Example 2
(1) Acquisition of regression equation The subjects were 120 healthy adults (63 males and 57 females) from 20 to 77 years old who can walk on their own without pain when walking. Subjects are healthy adults who do not have a history of injury or illness that may affect their walking behavior, and who can walk on two feet without complaining of pain at the time of measurement and without using walking aids or equipment. Met.
１０ｍほどの直線歩行が可能な、タイルカーペットが敷かれた実験室で次のように歩行計測を行った。 The walking was measured as follows in a laboratory where a tile carpet was laid, which was capable of walking straight about 10 m.
モーションキャプチャシステム（Vicon社製のVicon MXシステム、Vicon Nexus）を用いて、被験者体表面に貼付した計５７点のマーカ座標を２００Ｈｚで計測した。マーカの貼付位置はHelen Hayesマーカセットを基準とした。被験者らは、歩行路を裸足で端から端までまっすぐ歩行した。歩行時の速度や歩幅、視線などに関しては特に規定せず、普段通り歩行するよう指示した。各試行の開始は実験者が口頭で合図した。計測は、被験者の疲労を考慮しつつ、最低５試行計測した。 Using a motion capture system (Vicon MX system, Vicon Nexus manufactured by Vicon), a total of 57 marker coordinates attached to the subject's body surface were measured at 200 Hz. The marker attachment position was based on the Helen Hayes marker set. Subjects walked straight across the walking path barefoot. The walking speed, stride length, line of sight, etc. were not specified, and instructions were given to walk as usual. The start of each trial was verbally signaled by the experimenter. The measurement was performed for at least 5 trials while considering the fatigue of the subject.
回帰式を得るモデル構築のための前処理として、ノイズ除去と、腰（仙骨部）の３軸の加速度の抽出を実施例１と同様に行った。また、マーカ座標から取得した腰の３軸の加速度に重力の成分を加えることにより加速度ベクトルの大きさを算出した。 As preprocessing for constructing a model for obtaining a regression equation, noise removal and extraction of three-axis acceleration of the hip (sacral portion) were performed in the same manner as in the first embodiment. Further, the magnitude of the acceleration vector was calculated by adding a gravity component to the acceleration of the three axes of the waist acquired from the marker coordinates.
算出された加速度ベクトルの大きさの波形から一歩行周期分のデータを切り出し、一歩行周期を１００等分（時間率０〜１００％）することで時間正規化した。また、加速度に関しては、正規化していない加速度成分１０１個、体重で正規化した加速度成分１０１個、身長で正規化した加速度成分１０１個の計３０３個の成分をあらかじめ算出した。以上のことから１試行分のデータとして、最小つま先クリアランスの測定値一個と、腰の加速度ベクトルの大きさ３０３個を得た。
Data for one walking cycle was cut out from the waveform of the magnitude of the calculated acceleration vector, and the time was normalized by dividing one walking cycle into 100 equal parts (
健常成人１２０名×５試行分、計６００試行分のデータを用いて腰の加速度ベクトルの大きさから最小つま先クリアランスを推定するモデルを構築した。その際１２０名の被験者を、年齢や性別が等質となるように２群（Ａ群６０名とＢ群６０名）とし、「モデル構築群」と「モデル評価群」として分類した。最小つま先クリアランスは、遊脚中期におけるつま先の高さであるため、身体の質量の影響を受けることが考えられる。一方、骨盤は身体の重心が近いこともあり、腰の加速度は身体の質量の影響を受けにくいと考えられる。そこで最小つま先クリアランスについては体重で正規化し、体重の影響を取り除いた。 A model for estimating the minimum toe clearance from the magnitude of the waist acceleration vector was constructed using data for a total of 600 trials for 120 healthy adults × 5 trials. At that time, 120 subjects were divided into two groups (60 in group A and 60 in B group) so that their ages and genders were homogeneous, and were classified as a “model construction group” and a “model evaluation group”. Since the minimum toe clearance is the height of the toes in the middle stage of swing, it is considered that the minimum toe clearance is affected by the body mass. On the other hand, the pelvis may be close to the center of gravity of the body, and it is considered that the acceleration of the waist is hardly affected by the mass of the body. Therefore, the minimum toe clearance was normalized by body weight to remove the effect of body weight.
「モデル構築群」における最小つま先クリアランスを目的変数とし、説明変数として、上述したベクトル３０３成分を用いて重回帰分析（ステップワイズ法）を利用し、次の回帰式を得た。式中、ｔ％は時間率ｔ％の加速度ベクトルの大きさを表す。 The following regression equation was obtained by using multiple regression analysis (stepwise method) using the above-described vector 303 component as the explanatory variable and the minimum toe clearance in the “model building group” as the objective variable. In the equation, t% represents the magnitude of the acceleration vector at the time rate t%.
最小つま先クリアランス（cm）
＝0.055＋(24％:身長で正規化)×(-0.009)＋(51％：正規化せず)×(-0.007)＋(99％：正規化せず)×(0.002)＋ (77％:体重で正規化)×(0.421)＋(94％:正規化せず)×(0.005)
Minimum toe clearance (cm)
= 0.055 + (24%: normalized by height) x (-0.009) + (51%: not normalized) x (-0.007) + (99%: not normalized) x (0.002) + (77% : Normalized by weight) x (0.421) + (94%: not normalized) x (0.005)
この回帰式より算出されたモデル構築群における最小つま先クリアランスの推定値と最小つま先クリアランスの測定値との関係を図６に示す。これらの相関の程度を示す決定係数Ｒ2は0.561であり、この回帰式による最小つま先クリアランスの推定精度が高いことがわかる。 FIG. 6 shows the relationship between the estimated value of the minimum toe clearance and the measured value of the minimum toe clearance in the model building group calculated from the regression equation. Determination coefficient R 2 indicating the degree of correlation is 0.561, it is understood that the estimation accuracy of the minimum toe clearance by the regression equation is high.
一方、この回帰式により算出されたモデル評価群における最小つま先クリアランスの推定値と最小つま先クリアランスの測定値との関係を図７に示す。相関係数ｒは0.512となり、回帰式の推定精度の高さがわかる。 On the other hand, FIG. 7 shows the relationship between the estimated value of the minimum toe clearance and the measured value of the minimum toe clearance in the model evaluation group calculated by the regression equation. The correlation coefficient r is 0.512, which indicates that the estimation accuracy of the regression equation is high.
Claims (17)
加速度パラメータとして、最小つま先クリアランスを推定する足について、１歩行周期における立脚期中期の加速度、中期両足支持期の加速度、遊脚期前半の加速度、及び遊脚期後半の加速度から選ばれる加速度を使用するつまずきリスクの評価方法。 Minimum and toe clearance purposes variable, the longitudinal direction of the waist in walking, using a regression equation which includes the left-right direction and the vertical direction of the explanatory variable acceleration parameters derived from the acceleration, the minimum toe clearance for the right foot or left foot of the subjects' Is calculated, and the stumbling risk is evaluated.
As the acceleration parameter, for the foot for which the minimum toe clearance is to be estimated, an acceleration selected from the acceleration in the middle of the stance phase in one walking cycle, the acceleration in the mid-leg support phase, the acceleration in the first half of the swing phase, and the acceleration in the second half of the swing phase is used. How to evaluate stumbling risk.
被験者に携帯される加速度センサ、及び加速度センサで計測された加速度を用いて最小つま先クリアランスを算出し出力する演算装置を有し、
加速度センサは、歩行中の被験者の前後方向、左右方向及び鉛直方向の加速度を測定可能であり、
演算装置は、請求項１記載の回帰式を記憶する機能、
最小つま先クリアランスを推定する足について、踵着地から次の踵着地までの１歩行周期の加速度を時間分割して所定の時間間隔の加速度を取得する機能、
該所定時間間隔の加速度から、前記回帰式の加速度パラメータである立脚期中期の加速度、中期両足支持期の加速度、遊脚期前半の加速度、及び遊脚期後半の加速度から選ばれる加速度を抽出する機能、
抽出した加速度を前記回帰式で使用して最小つま先クリアランスを算出する機能、
を有する最小つま先クリアランスの推定装置。 A minimum toe clearance estimating device for estimating the minimum toe clearance from the acceleration of the waist of the subject during walking,
An acceleration sensor carried by the subject, and an arithmetic unit that calculates and outputs the minimum toe clearance using the acceleration measured by the acceleration sensor,
The acceleration sensor can measure the acceleration in the front-rear direction, the left-right direction, and the vertical direction of the subject while walking,
An arithmetic unit configured to store the regression equation according to claim 1;
For foot to estimate the minimum toe clearance, it functions the acceleration of 1 walking period from the heel landing until the next heel wear land by dividing time to get the acceleration of a predetermined time interval,
From the acceleration at the predetermined time interval, an acceleration selected from the acceleration parameters of the regression equation, ie, the acceleration in the middle stance phase, the acceleration in the mid-leg and foot support phase, the acceleration in the first half of the swing phase, and the acceleration in the second half of the swing phase is extracted. function,
A function of calculating the minimum toe clearance using the extracted acceleration in the regression equation,
Device for estimating the minimum toe clearance.
演算装置が、１歩行周期における所定時間間隔の加速度から、加速度パラメータである立脚期中期の左右方向の加速度、中期両足支持期の前後方向の加速度、遊脚期前半の前後方向の加速度及び左右方向の加速度、並びに遊脚期後半の前後方向の加速度及び左右方向の加速度を抽出する機能、及び
抽出した加速度を前記回帰式で使用して最小つま先クリアランスを算出する機能、
を有する請求項１１記載の最小つま先クリアランスの推定装置。 The regression equation stored in the arithmetic unit is defined as acceleration parameters as acceleration parameters in the left-right direction in the middle of the stance phase in one walking cycle, the acceleration in the front-rear direction during the mid-leg support period, the acceleration in the front-rear direction in the first half of the swing phase, and the left-right direction. Including the acceleration in the longitudinal direction and the lateral acceleration in the latter half of the swing phase,
The arithmetic unit calculates the acceleration parameters from the acceleration at a predetermined time interval in one walking cycle to the acceleration in the left-right direction during the middle of the stance phase, the acceleration in the front-rear direction during the mid-leg support period, the acceleration in the front-rear direction in the first half of the swing phase, and the left-right direction. The function of extracting the acceleration in the longitudinal direction and the lateral acceleration in the latter half of the swing period, and the function of calculating the minimum toe clearance using the extracted acceleration in the regression equation,
The minimum toe clearance estimating device according to claim 11, comprising:
演算装置が、被験者の１歩行周期における所定時間間隔の加速度から、加速度パラメータである主成分得点を算出する機能、及び
抽出した加速度と算出した主成分得点を回帰式で使用して最小つま先クリアランスを算出する機能、
を有する請求項１１又は１２記載の最小つま先クリアランスの推定装置。 The regression equation stored in the arithmetic unit is a principal component score calculated by principal component analysis of the acceleration in the front-rear direction, the left-right direction, and the vertical direction at predetermined time intervals in one walking cycle as the acceleration parameter, and the minimum toe Includes principal component scores that are highly correlated with clearance,
The arithmetic unit calculates a principal component score, which is an acceleration parameter, from the acceleration of the subject at a predetermined time interval in one walking cycle, and calculates a minimum toe clearance by using the extracted acceleration and the calculated principal component score in a regression equation. Function to calculate,
The apparatus for estimating a minimum toe clearance according to claim 11 or 12, comprising:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016034126A JP6658077B2 (en) | 2016-02-25 | 2016-02-25 | How to assess stumbling risk |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016034126A JP6658077B2 (en) | 2016-02-25 | 2016-02-25 | How to assess stumbling risk |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2017148287A JP2017148287A (en) | 2017-08-31 |
| JP6658077B2 true JP6658077B2 (en) | 2020-03-04 |
Family
ID=59740007
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016034126A Active JP6658077B2 (en) | 2016-02-25 | 2016-02-25 | How to assess stumbling risk |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6658077B2 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3530253B1 (en) | 2018-02-27 | 2021-06-02 | Jtekt Corporation | Walking assist device |
| US20230414130A1 (en) * | 2020-12-01 | 2023-12-28 | Nec Corporation | Walking index calculation device, walking index calculation system, walking index calculation method, and program recording medium |
| US20240350034A1 (en) * | 2021-09-27 | 2024-10-24 | Nec Corporation | Estimation device, information presentation system, estimation method, and recording medium |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4962765B2 (en) * | 2005-09-26 | 2012-06-27 | アイシン精機株式会社 | Walking analysis device and walking analysis method |
| US8647287B2 (en) * | 2008-12-07 | 2014-02-11 | Andrew Greenberg | Wireless synchronized movement monitoring apparatus and system |
| JP5115673B2 (en) * | 2010-03-16 | 2013-01-09 | 株式会社村田製作所 | Walking shoes |
| US9307932B2 (en) * | 2010-07-14 | 2016-04-12 | Ecole Polytechnique Federale De Lausanne (Epfl) | System and method for 3D gait assessment |
| US9524424B2 (en) * | 2011-09-01 | 2016-12-20 | Care Innovations, Llc | Calculation of minimum ground clearance using body worn sensors |
| US10945641B2 (en) * | 2011-10-09 | 2021-03-16 | The Medical Research, Infrastructure and Health Services Fund of the Tel Aviv Medical Center | Virtual reality for movement disorder diagnosis and/or treatment |
| JP5915990B2 (en) * | 2012-01-04 | 2016-05-11 | 国立研究開発法人産業技術総合研究所 | Stumbling risk assessment device, Stumbling risk assessment system |
-
2016
- 2016-02-25 JP JP2016034126A patent/JP6658077B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2017148287A (en) | 2017-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6183906B2 (en) | Gait estimation device and program, fall risk calculation device and program | |
| US7912672B2 (en) | Method and device for evaluating displacement signals | |
| JP6664746B2 (en) | Walking evaluation system and operation method of walking evaluation system | |
| JP6127873B2 (en) | Analysis method of walking characteristics | |
| JP6127455B2 (en) | Walking age evaluation method | |
| KR102235926B1 (en) | Gait Analysis System Using Smart Insole | |
| JP6829988B2 (en) | Momentum estimation device, momentum estimation program, and momentum estimation system | |
| JP6603055B2 (en) | Physical fitness estimation method | |
| JP6658077B2 (en) | How to assess stumbling risk | |
| JP6649323B2 (en) | Gait analysis system and method | |
| CN110495890A (en) | Detection method, detection device and detection system | |
| KR101746619B1 (en) | Apparatus for analysing gait and method for analysing gait using the same | |
| JP2012161402A (en) | Exercise characteristics evaluation system and exercise characteristics evaluation method | |
| KR102039236B1 (en) | Indicator Derivatives, Wearables, and Mobile Devices | |
| JP2020120807A (en) | Fall risk evaluation device, fall risk evaluation method and fall risk evaluation program | |
| JP6303746B2 (en) | Setting the walking target | |
| JP6127454B2 (en) | Walking age display system | |
| JP5915990B2 (en) | Stumbling risk assessment device, Stumbling risk assessment system | |
| JP6738249B2 (en) | Gait analysis method and gait analysis device | |
| WO2022249608A1 (en) | Ground-reaction-force-index estimation system, ground-reaction-force-index estimation method, and ground-reaction-force-index estimation program | |
| JP4723078B2 (en) | Body fat scale | |
| JP6552875B2 (en) | Moving motion analysis apparatus, method and program | |
| CN112839569B (en) | Method and system for assessing human movement | |
| JP6917047B2 (en) | Motor function evaluation system with multiple wearable sensors | |
| CN113040752A (en) | Exercise amount monitoring method and system based on heart rate |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20181203 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20190814 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190827 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20191024 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20200107 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200120 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 6658077 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |