(57)【要約】【課題】ユーザに分かり易い指標を用いて歩行能力を評価する装置を提供する。【解決手段】歩行評価装置２は、被験者から複数の歩行データを取得する歩行データ取得手段２２１と、歩行データから、被験者の歩行状態に関連する複数の歩行指標を算出する歩行指標算出手段２２２と、予め蓄積された複数のサンプル者の蓄積歩行データ２３１に基づき、歩行指標の各々に対応する歩行指標スコアを算出する歩行指標スコア算出手段２２３と、蓄積歩行データ２３１に基づき、分類された複数の関連グループに属する歩行指標スコアから、関連グループ毎に、被験者のサブスコアを算出するサブスコア算出手段２２４と、関連グループ毎に算出される被験者のサブスコアに基づいて、被験者の歩行スコアを算出する歩行スコア算出手段２２５と、を備える。【選択図】図３
2020151470.tif 000002
【特許請求の範囲】
【請求項１】
被験者から複数の歩行データを取得する歩行データ取得手段と、
前記歩行データから、前記被験者の歩行状態に関連する複数の歩行指標を算出する歩行指標算出手段と、
予め蓄積された複数のサンプル者の蓄積歩行データに基づき、前記歩行指標の各々に対応する歩行指標スコアを算出する歩行指標スコア算出手段と、
前記蓄積歩行データに基づき、分類された複数の関連グループに属する前記歩行指標スコアから、前記関連グループ毎に、前記被験者のサブスコアを算出するサブスコア算出手段と、
前記関連グループ毎に算出される前記被験者の前記サブスコアに基づいて、前記被験者の歩行スコアを算出する歩行スコア算出手段と、
を備える歩行評価装置。
【請求項２】
前記歩行スコア算出手段は、前記関連グループ毎に算出される前記被験者の前記サブスコアの総和を、前記被験者の前記歩行スコアとして算出する、請求項１に記載の歩行評価装置。
【請求項３】
前記歩行スコア算出手段は、前記関連グループ毎に重み付けを行い、重み付けされた前記関連グループ毎のサブスコアの総和を、前記歩行スコアとして算出する、請求項２に記載の歩行評価装置。
【請求項４】
前記歩行指標スコア算出手段は、前記複数のサンプル者の前記蓄積歩行データに基づく前記歩行指標の分布から、前記被験者の歩行指標をスコア化して歩行指標スコアを算出する、請求項１から３のいずれかに記載の歩行評価装置。
【請求項５】
前記複数の関連グループの分類は、因子分析によりなされている請求項１から４のいずれかに記載の歩行評価装置。
【請求項６】
被験者から複数の歩行データを取得する歩行データ取得ステップと、
前記歩行データから、前記被験者の歩行状態に関連する複数の歩行指標を算出する歩行指標算出ステップと、
予め蓄積された複数のサンプル者の蓄積歩行データに基づき、前記歩行指標の各々に対応する歩行指標スコアを算出する歩行指標スコア算出ステップと、
前記蓄積歩行データに基づき、分類された複数の関連グループに属する前記歩行指標スコアから、前記関連グループ毎に、前記被験者のサブスコアを算出するサブスコア算出ステップと、
前記関連グループ毎に算出される前記被験者の前記サブスコアに基づいて、前記被験者の歩行スコアを算出する歩行スコア算出ステップと、
を含む歩行評価方法。
【請求項７】
請求項１から５のいずれかに記載の歩行評価装置の各手段として、コンピュータを機能させるためのプログラム。 
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被験者の歩行能力を評価する歩行評価装置、歩行評価方法およびプログラムに関する。
【背景技術】
【０００２】
従来より、歩行能力を評価する指標（以下、歩行指標と呼ぶ）には、据え置き型の歩行計測装置から得られる歩行指標と、ポータブル型の小型センサを用いた歩行計測装置から得られる歩行指標とが用いられている。
【０００３】
据え置き型の歩行計測装置の一例を特許文献１に示す。据え置き型の歩行計測装置には、床面に圧力センサを敷き詰めた歩行路や、負荷センサを内蔵した床反力計がある。このような歩行路または床反力計により得られるデータは、ステップ時間と位置とに依存した時間・距離データ（例えば、ステップ時間やステップ幅など）や、関節の周りに生じる回転力（関節トルクと呼ばれる）に関するデータであり、歩行指標としては非常にスタンダードであることから、一般的な研究ではよく用いられている。
【０００４】
小型センサを用いた歩行計測装置の一例を特許文献２および３に示す。歩行計測に用いられている小型センサには、加速度センサや角速度センサがある。近年では、多チャンネルで同時に複数方向の計測が可能なセンサが主流となっている。センサからの電波が届く範囲であれば何歩でも計測することが可能であることや、センサから得られる歩行データに適切な波形処理を施すことにより歩行動作の規則性や円滑性等に関する指標も抽出することができることから、近年では、小型センサは高齢者の歩行動作に関する研究にも用いられている。
【０００５】
センサを用いてパーキンソン病等の脳神経性疾患の自動診断を行う装置の一例を特許文献４に示す。特許文献４の装置によると、患者に取り付けられたセンサにより患者の運動を測定し、測定により得られたデータの特徴量に基づいて診断データを生成している。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００１−２１８７５４号公報
【特許文献２】特開２０１８−１４３４１２号公報
【特許文献３】特開２０１７−２２１５０２号公報
【特許文献４】国際公開第２０１７／０６５２４１号
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、特許文献１に例示するような据え置き型の歩行計測装置では、被験者の歩行指標を計測することが可能な距離が、例えば最大で５メートル程度と非常に短いため、歩行動作として近年注目されている歩行動作のばらつき指標（例えば、ストライド時間の変動係数など）はそもそも計算することができない。また、据え置き型の歩行計測装置は、研究所等を中心とした実験室ベースの機器であることから、高価であり操作性にも欠けており、家庭用の一般ユーザへの普及が期待できない。
【０００８】
小型センサを用いた歩行計測装置であれば、歩行指標を抽出することができる。しかしながら、小型センサから得られる歩行指標の意味を理解するためには専門的な知識が必要であるため、一般ユーザにとっては歩行指標の理解が困難であるという問題が依然として存在している。そのため、例えば特許文献２および３に例示するシステムにおいても、ユーザは、自身の歩行状態が望ましくない歩行状態であることを、システムから発せられる光、音、振動またはメッセージ等により把握しているのが実情である。
【０００９】
特許文献４の装置においても、生成される診断データはパーキンソン病に関連するHoehn-Yahr分類の度数（０度、Ｉ度、ないしＶ度）であり、患者が診断データの意味を理解するのは困難である。
【００１０】
上記した以外の、例えばビデオカメラ等の撮影装置による歩行計測により得られる歩行指標についても、上述した課題と同様の課題を有しており、健康産業やリハビリテーションにおける歩行解析の一般化を阻んでいる原因となっている。歩行能力が低下している例えば高齢者や術後のリハビリテーション中の患者に対して、具体的にどのような歩行動作を行えば望ましい歩行状態となりリハビリテーションの効果が向上しているのかを、ユーザに分かり易く提示することが求められている。
【００１１】
本発明は、上記課題を解決するためになされたものであり、その目的は、ユーザに分かり易い指標を用いて歩行能力を評価する装置、方法およびプログラムを提供することにある。
【課題を解決するための手段】
【００１２】
上記目的を達成するための本発明は、例えば以下に示す態様を含む。
（項１）
被験者から複数の歩行データを取得する歩行データ取得手段と、
前記歩行データから、前記被験者の歩行状態に関連する複数の歩行指標を算出する歩行指標算出手段と、
予め蓄積された複数のサンプル者の蓄積歩行データに基づき、前記歩行指標の各々に対応する歩行指標スコアを算出する歩行指標スコア算出手段と、
前記蓄積歩行データに基づき、分類された複数の関連グループに属する前記歩行指標スコアから、前記関連グループ毎に、前記被験者のサブスコアを算出するサブスコア算出手段と、
前記関連グループ毎に算出される前記被験者の前記サブスコアに基づいて、前記被験者の歩行スコアを算出する歩行スコア算出手段と、
を備える歩行評価装置。
（項２）
前記歩行スコア算出手段は、前記関連グループ毎に算出される前記被験者の前記サブスコアの総和を、前記被験者の前記歩行スコアとして算出する、項１に記載の歩行評価装置。
（項３）
前記歩行スコア算出手段は、前記関連グループ毎に重み付けを行い、重み付けされた前記関連グループ毎のサブスコアの総和を、前記歩行スコアとして算出する、項２に記載の歩行評価装置。
（項４）
前記歩行指標スコア算出手段は、前記複数のサンプル者の前記蓄積歩行データに基づく前記歩行指標の分布から、前記被験者の歩行指標をスコア化して歩行指標スコアを算出する、項１から３のいずれかに記載の歩行評価装置。
（項５）
前記複数の関連グループの分類は、因子分析によりなされている項１から４のいずれかに記載の歩行評価装置。
（項６）
被験者から複数の歩行データを取得する歩行データ取得ステップと、
前記歩行データから、前記被験者の歩行状態に関連する複数の歩行指標を算出する歩行指標算出ステップと、
予め蓄積された複数のサンプル者の蓄積歩行データに基づき、前記歩行指標の各々に対応する歩行指標スコアを算出する歩行指標スコア算出ステップと、
前記蓄積歩行データに基づき、分類された複数の関連グループに属する前記歩行指標スコアから、前記関連グループ毎に、前記被験者のサブスコアを算出するサブスコア算出ステップと、
前記関連グループ毎に算出される前記被験者の前記サブスコアに基づいて、前記被験者の歩行スコアを算出する歩行スコア算出ステップと、
を含む歩行評価方法。
（項７）
項１から５のいずれかに記載の歩行評価装置の各手段として、コンピュータを機能させるためのプログラム。
【発明の効果】
【００１３】
本発明によると、ユーザに分かり易い指標を用いて歩行能力を評価する装置、方法およびプログラムを提供することができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の一実施形態に係る歩行評価システムの概略的な構成を模式的に示す図である。
【図２】センサにより取得される歩行データの一例を示す。
【図３】本発明の一実施形態に係る歩行評価装置の機能を説明するためのブロック図である。
【図４】本発明の一実施形態に係る歩行評価装置が行うデータ処理の手順を示すフローチャートである。
【図５】歩行指標値の分布を示す模式図である。
【図６】歩行スコアおよびサブスコアの表示例である。
【図７】実施例２に係る歩行スコアのヒストグラム及び正規分位点プロットである。
【発明を実施するための形態】
【００１５】
以下、本発明の実施形態を、添付の図面を参照して詳細に説明する。なお、以下の説明および図面において、同じ符号は同じまたは類似の構成要素を示すこととし、よって、同じまたは類似の構成要素に関する重複した説明を省略する。
［システム構成］
【００１６】
図１は、本発明の一実施形態に係る歩行評価システムの概略的な構成を模式的に示す図である。
【００１７】
一実施形態に係る歩行評価システム１０は、被験者９に取り付けられて被験者９の歩行データ（第１の歩行データ）を測定するセンサ１（１ａ〜１ｃ）と、センサ１から送信される歩行データに基づいて、被験者９の歩行能力を評価する歩行評価装置２と、を備える。
【００１８】
本実施形態では、センサ１（１ａ〜１ｃ）と歩行評価装置２とは、bluetooth（登録商標）やWiFi（登録商標）等の無線によって接続されているが、センサ１と歩行評価装置２との接続方法は特に無線に限定されず、有線接続であってもよい。何れの場合にも、歩行評価装置２には、汎用のパーソナルコンピュータや、スマートフォン等を用いることができる。また、歩行評価システム１０は、歩行評価装置２にネットワーク３で接続されたサーバ４をさらに備えることができる。サーバ４には、歩行評価装置２が備える機能ブロックの一部をクラウド化することができる。
【００１９】
センサ１は、加速度または角速度を感知し、感知した加速度または角速度の時系列の波形データである歩行データとして、歩行評価装置２に送信する。本実施形態では、３つのセンサ１（１ａ〜１ｃ）が被験者９に取り付けられている。センサ１ａとしては、右足の踵付近に取り付けられた三軸角速度センサを用い、センサ１ｂとしては、左足の踵付近に取り付けられた三軸角速度センサを用いる。また、センサ１ｃとしては、被験者９の体幹付近、例えば重心位置に対応する第三腰椎棘突起付近に取り付けられた三軸加速度センサを用いる。なお、以下においてセンサ１が測定する加速度または加速度の三軸とは、被験者９の垂直方向、側方および前後方向を意味する。
【００２０】
図２に、センサにより取得される歩行データの一例を示す。（Ａ）は体幹位置センサ１ｃにより取得された加速度に関する歩行データであり、測定方向は垂直方向である。（Ｂ）は、右足踵位置センサ１ａにより取得された角速度に関する歩行データであり、（Ｃ）は、左足踵位置センサ１ｂにより取得された角速度に関する歩行データである。なお、図２には例示していないが、例えば体幹位置センサ１ｃは、垂直方向に限らず、側方および前後方向の測定方向についても、加速度に関する歩行データを取得することができる。
［装置の構成］
【００２１】
図３は、本発明の一実施形態に係る歩行評価装置の機能を説明するためのブロック図である。
【００２２】
一実施形態に係る歩行評価装置２は、ハードウェアの構成として、通信インタフェース部（通信Ｉ／Ｆ部）２１と、データ処理部２２と、補助記憶装置２３と、表示部２４とを備えている。なお、図示していないが、歩行評価装置２は、他のハードウェアの構成として、データ処理を行うＣＰＵ等のプロセッサと、プロセッサがデータ処理の作業領域に使用するメモリとをさらに備えている。
【００２３】
通信インタフェース部２１は、センサ１やサーバ４等の外部機器とのデータの送受信を行う。本実施形態では、センサ１から送信される歩行データの受信を行う。
【００２４】
補助記憶装置２３は、オペレーティングシステム（ＯＳ）、各種制御プログラム、および、プログラムによって生成されたデータなどを記憶する不揮発性の記憶装置であり、例えば、ｅＭＭＣ（embedded Multi Media Card）やＳＳＤ（Solid State Drive）等によって構成される。
【００２５】
本実施形態では、補助記憶装置２３には、サンプル歩行データ２３１（以下、第２の歩行データまたは蓄積歩行データとも記載する）と、歩行評価プログラム２３２とが記憶されている。サンプル歩行データ２３１は、センサ１を用いて事前に所定数のサンプルに対して予め測定された歩行データであり、本実施形態では、補助記憶装置２３に予め記憶されている。例示的には、サンプルとは、健康上の大きな問題を有さない高齢者であり、本実施形態では、サンプル歩行データ２３１は、このようなサンプルの３７８人分の歩行データである。
【００２６】
歩行評価プログラム２３２は、ソフトウェアによる機能ブロックである後述するデータ処理部２２を実現するためのコンピュータプログラムである。歩行評価プログラム２３２は、インターネット等のネットワーク３を介して歩行評価装置２にインストールしてもよい。あるいは、歩行評価プログラム２３２を記録したメモリカード等のコンピュータ読み取り可能な非一時的な有体の記録媒体を歩行評価装置２に読み取らせることにより、歩行評価プログラム２３２を歩行評価装置２にインストールしてもよい。サンプル歩行データ２３１も歩行評価プログラム２３２と同様に、ネットワーク３または記録媒体等を介して、歩行評価装置２にインストールしてもよい。
【００２７】
表示部２４は、例えば液晶ディスプレイおよび有機ＥＬディスプレイ等で構成することができる。なお、表示部２４を、歩行評価装置２とは別個の装置として構成してもよい。
【００２８】
また、歩行評価装置２は、ソフトウェアの構成としてデータ処理部２２を備えている。データ処理部２２は、プロセッサが歩行評価プログラム２３２を実行することにより実現される機能ブロックである。
【００２９】
データ処理部２２は、被験者９の歩行データに基づいて、歩行状態に関連する複数の歩行指標のそれぞれに対応する歩行指標スコアを算出し、算出した歩行指標スコアに基づいて、被験者９の歩行状態を示す歩行スコアを算出する機能を有している。この機能を実現するために、データ処理部２２は、歩行データ取得部２２１と、歩行指標算出部２２２と、歩行指標スコア算出部２２３と、サブスコア算出部２２４と、歩行スコア算出部２２５と、歩行指標分類部２２６と、を備えている。なお、データ処理部２２を、集積回路上に形成された論理回路によってハードウェア的に実現してもよい。
【００３０】
歩行データ取得部２２１は２種類の歩行データを取得する。後述する歩行指標を予め分類する際には、歩行データ取得部２２１は、複数のサンプルから予め取得されたサンプル歩行データ２３１（第２の歩行データ）を取得する。被験者９の歩行能力を評価する際には、歩行データ取得部２２１は、通信Ｉ／Ｆ部２１を介して、被験者９に取り付けられたセンサ１（１ａ〜１ｃ）から歩行データ（第１の歩行データ）を取得する。歩行データ取得部２２１は、取得した歩行データに対して所定の整流処理を施すことができる。整流処理としては、カットオフ周波数が２０Ｈｚのローパスフィルタ（バターワースフィルタ）を用いることができる。
【００３１】
歩行指標算出部２２２は、取得した歩行データに基づいて、歩行状態に関連する複数の歩行指標を算出する。歩行指標スコア算出部２２３は、算出された複数の歩行指標のそれぞれに対応する歩行指標スコアを算出する。なお、歩行指標算出部２２２および歩行指標スコア算出部２２３はいずれも、サンプル歩行データ２３１と、被験者９から取得された歩行データとの両方を取り扱うことができる。
【００３２】
サブスコア算出部２２４および歩行スコア算出部２２５は、算出された複数の歩行指標スコアに基づいて、被験者９の歩行状態を示す歩行スコアを算出する。サブスコア算出部２２４は、関連グループごとに歩行指標スコアに基づいてサブスコアを算出する。複数の歩行指標は、複数のサンプルから予め取得されたサンプル歩行データ２３１に基づいて、歩行指標分類部２２６により複数の関連グループに分類されており、サブスコア算出部２２４は、関連グループごとに歩行指標スコアに基づいてサブスコアを算出する。歩行スコア算出部２２５は、サブスコア算出部２２４によって関連グループごとに算出されたサブスコアに基づいて、歩行スコアを算出する。本実施形態では、サブスコア算出部２２４は、関連グループごとに歩行指標スコアの和をサブスコアとして算出し、歩行スコア算出部２２５は、各サブスコアの総和を歩行スコアとして算出する。
【００３３】
歩行指標分類部２２６は、複数の歩行指標のそれぞれに対応する歩行指標スコアに基づいて、複数の歩行指標を複数の関連グループに分類する。複数の歩行指標のそれぞれに対応する歩行指標スコアは、歩行指標スコア算出部２２３により、複数のサンプルから予め取得されたサンプル歩行データ２３１に基づいて算出されている。本実施形態では、歩行指標分類部２２６は、歩行指標スコアに対して因子分析を行うことにより、複数の歩行指標を複数の関連グループに分類する。
［歩行指標］
【００３４】
歩行評価装置２は、被験者９の左右の踵付近および体幹付近に取り付けられたセンサ１（１ａ〜１ｃ）を用いて測定された歩行データから、以下に例示する種々の歩行指標を算出することができる。
【００３５】
本実施形態では、以下に例示する、歩行速度［歩行指標No.1］、ストライド時間（ST）［歩行指標No.2］、ストライド時間のばらつき（STV）［歩行指標No.3］、遊脚期時間のばらつき（SwTV）［歩行指標No.4］、三軸各方向の規則性（AC-VT, AC-ML, AC-AP）［歩行指標No.5〜No.7］、および三軸各方向の円滑性（HR-VT, HR-ML, HR-AP）［歩行指標No.8〜No.10］、の１０種類の歩行指標のそれぞれに対応する歩行指標スコアを算出し、算出した歩行指標スコアに基づいて、被験者９の歩行状態を示す歩行スコアを算出する。
【００３６】
なお、歩行指標に関する以下の説明において、足部とは例えば左右の踵を意味し、体幹とは例えば腰を意味する。
【００３７】
・歩行速度［歩行指標No.1］
踵の前後方向の連続した加速度波形から重力加速度をキャンセルし、ドリフトを取り除いた後に得られたデータを積分することにより算出する。
【００３８】
歩行速度の詳細については、文献「Mariani B, Hoskovec C, Rochat S, Bula C, Penders J, Aminian K. 3D gait assessment in young and elderly subjects using foot-worn inertial sensors. Journal of Biomechanics. 2010; 16;43(15):2999-3006」を参照されたい。
【００３９】
・ストライド時間（ST: Stride Time）［歩行指標No.2］
・ストライド時間のばらつき（STV: Stride Time Variability: STV）［歩行指標No.3］
左右それぞれの踵の角速度波形からＭ字型の特徴的なピークを同定し、歩行中のヒールコンタクトの時間を抽出する。ヒールコンタクトとは、歩行のステップにおいて踵が地面と接地することを意味する。Ｍ字型を構成する左のピークがヒールコンタクトであるので、その時間データの差分をストライド時間として、例えば５ストライド分のストライド時間を左右それぞれの踵について計算する。左右で５ストライド分の平均値を計算しストライド時間（ST）とする。また、算出したストライド時間の標準偏差と平均値とから変動係数を求めて１００倍したものをストライド時間のばらつき（STV）とする。
【００４０】
・遊脚期時間（SwT: Swing Time）および遊脚期時間のばらつき（SwTV: Swing Time Variability）［歩行指標No.4］
足部の角速度波形から特徴的なピークを同定することにより算出する。波形において隣り合うピーク間の差分はそれぞれ立脚時間と遊脚時間となり、いずれか一方が決まればもう一方が決まる。そのため、どちらか一方だけを歩行指標とすればよい。本実施形態では遊脚時間にフォーカスして、その平均値およびそのばらつきの程度を算出する。計算式は上記したストライド時間のばらつき（STV）と同様である。
【００４１】
・規則性（AC: AutocorrelationまたはRegularity）［歩行指標No.5〜No.7］
ストライド間の規則性を意味する。自己相関係数とも呼ばれる。自己相関係数は、腰の三方向（垂直方向、側方、前後方向）の加速度波形と、左右両踵の三面（矢状面、水平面、前額面）の角速度波形とから算出する。自己相関係数は、複製した波形と例えば１ストライド分波形データをずらした際の相関係数であり、０か１までの値をとる。
【数１】
2020151470.tif 000003
【００４２】
ここで、Ｎは解析対象となったデータの数（加速度データまたは加速度データの数）であり、ｍ＝０，±１，±２，…，±Ｎ−１である。
【００４３】
自己相関係数の詳細については、文献「Rolf Moe-Nilssen, Jorunn L. Helbostad. Estimation of gait cycle characteristics by trunk accelerometry. Journal of Biomechanics 2004; 37: 121-126」を参照されたい。
【００４４】
・円滑性（HR: Harmonic RatioまたはSmoothness）［歩行指標No.8〜No.10］
ストライド内の流暢性またはステップの対称性を意味する。腰の三方向（垂直方向、側方、前後方向）の加速度波形からHarmonic Ratioを求める。各加速度波形をストライドごとの例えば５つの波形に分け、各波形に対して離散フーリエ変換を行う。１ストライドは２ステップであることから、ノイズを含まないストライド内の波形は２の倍数で構成された波形が含まれていることになる。よって、２，４，６，…と続く偶数のharmonicsのパワーを合算したものが、ノイズのない波形のパワーを表し、それ以外の１，３，５，…と続く奇数のharmonicsのパワーを合算したものが、ノイズ成分のパワー、すなわち本来なら含まれるはずがない動きの成分のパワーを表す。よって、奇数のパワーの合算値と偶数のパワーの合算値との比率がHarmonic Ratioであり、歩行中の体幹における衝撃の流暢性もしくはステップの対称性を表す。
【００４５】
円滑性の詳細については、文献「Menz, Hylton B., Stephen R. Lord, and Richard C. Fitzpatrick. Acceleration patterns of the head and pelvis when walking on level and irregular surfaces. Gait posture 2003; 18: 35-46.」を参照されたい。
［手順］
【００４６】
図４は、本発明の一実施形態に係る歩行評価装置が行うデータ処理の手順を示すフローチャートである。
【００４７】
ステップＳ１〜Ｓ４の処理は、複数のサンプルから取得されたサンプル歩行データ２３１（第２の歩行データ）に基づいて、複数の歩行指標を関連グループに分類する処理である。ステップＳ５〜Ｓ９の処理は、複数の関連グループに分類された歩行指標の情報と、被験者９の歩行データ（第１の歩行データ）とに基づいて、被験者９の歩行スコアを算出する処理である。
【００４８】
ステップＳ１において、サンプルの歩行データを取得する。本実施形態では、事前に測定されている３７８人分のサンプルの歩行データであるサンプル歩行データ２３１は、補助記憶装置２３に予め記憶されている。歩行データ取得部２２１は、補助記憶装置２３からサンプル歩行データ２３１を取得する。
【００４９】
ステップＳ２において、サンプルの歩行データに基づいて、各サンプルの歩行指標No.1〜No.10を算出する。歩行指標算出部２２２は、サンプル歩行データ２３１に含まれている３７８人分のそれぞれのサンプルについて、歩行データに基づいて１０種類の歩行指標を算出する。
【００５０】
ステップＳ３において、サンプルの歩行指標に対応する歩行指標スコアを算出する。歩行指標スコア算出部２２３は、算出した１０種類の歩行指標のそれぞれの分布に基づいて、各サンプルの歩行指標No.1〜No.10に対応する歩行指標スコアを算出する。
【００５１】
図５は、歩行指標値の分布を示す模式図である。図５を参照して、歩行指標の値と歩行指標スコアとの対応関係を説明する。
【００５２】
１０種類の歩行指標のうち、歩行速度［歩行指標No.1］を一例として説明する。事前に測定されている３７８人分のサンプル歩行データ２３１から算出される歩行速度について統計処理が行われ、歩行速度について代表値（２５％四分位、中央値、および７５％四分位）が求められる。この代表値に従って、図５に示すように歩行速度の歩行指標値の分布を例えば４つの区分に分割し、それぞれの区分と歩行指標スコアとを対応させる。これにより、サンプル毎に歩行速度は０〜３の範囲の整数値をとる歩行指標スコアとして表現される。例えば、或るサンプルの歩行速度が、中央値から７５％四分位の区分内の値である場合、そのサンプルの歩行速度の歩行指標スコアは「２」と表現される。
【００５３】
歩行速度以外の他の９種類の歩行指標［歩行指標No.2〜No.10］についてもこの歩行速度の例と同様に、０〜３の範囲の整数値をとる歩行指標スコアとして表現される。一例として、サンプル歩行データ２３１について１０種類の各歩行指標の歩行指標スコアを算出した結果の一例を表１に示す。
【表１】
2020151470.tif 000004
【００５４】
ステップＳ４において、サンプルの歩行指標に対応する歩行指標スコアに基づいて、歩行指標を関連グループに分類する。歩行指標分類部２２６は、３７８人分のそれぞれのサンプルについてそれぞれの歩行指標について算出した歩行指標スコアに基づいて、歩行指標を複数の関連グループに分類する。
【００５５】
本実施形態では、歩行指標分類部２２６は、表１に例示するデータ行列の歩行指標スコアに対して、プロマックス回転を伴う因子分析を行う。プロマックス回転は、因子分析における斜交回転法の一つであり、プロマックス回転によると、因子間に相関が生じることが許容される。なお、プロマックス回転に代えて、バリマックス回転を伴う因子分析を行ってもよい。バリマックス回転は、因子分析における直交回転法の一つである。
【００５６】
因子分析を行ううえで、本実施形態では、因子数を決定する際の固有値の下限は１とし、因子負荷量が０．５未満の変数（歩行指標）は歩行指標スコアより除外する。すなわち、どの歩行指標の関連グループにも分類されず、歩行動作の特徴を表すパラメータとして必要ないと判断される歩行指標は、この段階において除外される。因子分析を行った結果である因子負荷行列を表２に示す。
【表２】
2020151470.tif 000005
【００５７】
表２を参照して、歩行指標を関連グループに分類する手順について具体的に説明する。表２に示されているそれぞれの負荷の値から、因子１については、垂直方向規則性（AC-VT）、側方規則性（AC-ML）、および前後方向規則性（AC-AP）が、関連グループを構成していると判断される。関連グループ内に含まれる歩行指標の特徴から、因子１の歩行指標の関連グループは、歩行の規則性に関連する関連グループであると判断される。すなわち、垂直方向規則性（AC-VT）、側方規則性（AC-ML）、および前後方向規則性（AC-AP）のそれぞれの歩行指標スコアは、歩行の規則性に関連するスコアであると判断される。本実施形態では、各歩行指標は０〜３の範囲の整数値を取る歩行指標スコアとして表現されるため、歩行の規則性に関連するスコアは、３種類の歩行指標により９点満点で表される。
【００５８】
同様に、因子２については、歩行速度とストライド時間とが、関連グループを構成していると判断される。関連グループ内に含まれる歩行指標の特徴から、これらは歩調に関連する関連グループを構成していると判断される。歩調に関連するスコアは６点満点で表される。また、因子３については、ストライド時間のばらつきと遊脚期時間のばらつきとが、関連グループを構成していると判断される。関連グループ内に含まれる歩行指標の特徴から、これらは歩行のばらつきに関連する関連グループを構成していると判断される。歩行のばらつきに関連するスコアは６点満点で表される。更には、因子４については、垂直方向円滑性（HR-VT）、側方円滑性（HR-ML）、および前後方向円滑性（HR-AP）が、関連グループを構成していると判断される。関連グループ内に含まれる歩行指標の特徴から、これらは歩行のなめらかさに関連する関連グループを構成していると判断される。歩行のなめらかさに関連するスコアは、９点満点で表される。
【００５９】
以上により、サンプルの歩行データに基づいた統計的な分析から、１０種類の歩行指標は、本実施形態では、歩行の規則性（RegularityまたはAutocorrelation）、歩調（Pace）、歩行のばらつき（Variability）、および歩行のなめらかさ（SmoothnessまたはHarmonic Ratio）の４つの関連グループに分類される。これら４つの関連グループのスコアはそれぞれ、９点満点、６点満点、６点満点、および９点満点で表されることから、以下のステップＳ５〜Ｓ９の処理により算出する、１０種類の歩行指標に関連する被験者９の歩行スコアは、３０点満点で表される。
【００６０】
ステップＳ５において、被験者９から歩行データを取得する。歩行データ取得部２２１は、通信Ｉ／Ｆ部２１を介して、被験者９に取り付けられたセンサ１（１ａ〜１ｃ）から歩行データを取得する。
【００６１】
ステップＳ６において、被験者９の歩行データに基づいて、被験者の歩行指標No.1〜No.10を算出する。歩行指標算出部２２２は、被験者９から取得した歩行データから、被験者９についての１０種類の歩行指標を算出する。
【００６２】
ステップＳ７において、被験者９の歩行指標に対応する歩行指標スコアを算出する。歩行指標スコア算出部２２３は、算出されたそれぞれの歩行指標の値を、図５を参照して説明した歩行指標の値と歩行指標スコアとの対応関係に基づいて、歩行指標スコアとして表現する。
【００６３】
ステップＳ８において、ステップＳ３と同様なプロセスを経て、被験者９のサブスコアを関連グループ毎に算出する。サブスコア算出部２２４は、被験者９の各歩行指標に関する歩行指標スコアと、ステップＳ４において複数の関連グループに分類された歩行指標の情報とに基づいて、歩行指標の関連グループ毎にサブスコアを算出する。
【００６４】
本実施形態では、歩行指標スコア算出部２２３は、被験者９の垂直方向規則性（AC-VT）、側方規則性（AC-ML）、および前後方向規則性（AC-AP）のそれぞれの歩行指標の値から、それぞれに対応する歩行指標スコアを算出する。サブスコア算出部２２４は、算出された３つの歩行指標スコアの和を算出することにより、歩行の規則性に関連するサブスコアを算出する。
【００６５】
同様に、歩行指標スコア算出部２２３は、被験者９の歩行速度およびストライド時間のそれぞれの歩行指標の値から、それぞれに対応する歩行指標スコアを算出する。サブスコア算出部２２４は、算出された２つの歩行指標スコアの和を算出することにより、歩調に関連するサブスコアを算出する。
【００６６】
歩行指標スコア算出部２２３は、被験者９のストライド時間のばらつきおよび遊脚期時間のばらつきのそれぞれの歩行指標の値から、それぞれに対応する歩行指標スコアを算出する。サブスコア算出部２２４は、算出された２つの歩行指標スコアの和を算出することにより、歩行のばらつきに関連するサブスコアを算出する。
【００６７】
歩行指標スコア算出部２２３は、被験者９の垂直方向円滑性（HR-VT）、側方円滑性（HR-ML）、および前後方向円滑性（HR-AP）のそれぞれの歩行指標の値から、それぞれに対応する歩行指標スコアを算出する。サブスコア算出部２２４は、算出された３つの歩行指標スコアの和を算出することにより、歩行のなめらかさに関連するサブスコアを算出する。
【００６８】
ステップＳ９において、被験者９の歩行スコアを算出する。歩行スコア算出部２２５は、ステップＳ８において算出された４つのサブスコアに基づいて、歩行スコアを算出する。本実施形態では、図６に例示するように、歩行スコア算出部２２５は、歩行の規則性に関連するサブスコア（９点満点）と、歩調に関連するサブスコア（６点満点）と、歩行のばらつきに関連するサブスコア（６点満点）と、歩行のなめらかさに関連するサブスコア（９点満点）との総和を算出して、被験者９の歩行スコア（３０点満点）とする。算出された被験者９の歩行スコア２４１は表示部２４に表示される。またこの際、歩行スコア２４１とともに、歩行スコア２４１の内訳である４つのサブスコア２４２を、関連グループ名２４３とともに表示部２４に表示してもよい。図６の例示においては、被験者９のサブスコアは、規則性グループで６点、歩調グループで４点、ばらつきグループで５点、なめらかさグループで８点と示され、歩行スコアは２３点となっている。
【００６９】
好ましくは、歩行スコア２４１は、被験者９が取り得る満点のスコア表示（図示する例では３０点）とともに表示される。好ましくは、各サブスコア２４２は、各関連グループが取り得る満点のスコア表示とともに表示される。
［効果］
【００７０】
以上、本発明によると、被験者９の歩行状態に関連する歩行スコアを算出することにより、ユーザに分かり易い指標を用いて歩行能力を評価することができる。本発明によると、歩行能力が低下している例えば高齢者や術後のリハビリテーション中の患者は、歩行スコアの増加または減少を通じて、経時的なリハビリテーションの効果の有無を容易に把握することができる。また、運動指導をしているセラピストは、今後実施するべき適切な運動プログラムを作成することができる。歩行スコアは、複数の歩行指標が分類されることにより、歩行能力が多面的に評価された包括的なスコアの数値として算出され表示されるので、ユーザは、小型センサから得られる歩行指標の意味を理解するための専門的な知識を得る必要はない。
【００７１】
また、歩行スコアの内訳であるサブスコアを表示することにより、ユーザは、自身の歩行動作のどのような観点について歩行スコアが増加または減少しているのかを容易に把握することできる。これにより、具体的にどのような歩行動作をすれば良いのかといった、リハビリテーションの方向性もユーザは容易に把握することができる。
【００７２】
本発明によると、複数の歩行指標は、サンプルの歩行データに基づいた統計的な分析から複数の関連グループに分類されており、ユーザに提示される歩行スコアも、関連グループ毎のスコアであるサブスコアに基づいて算出されている。これにより、歩行スコアは歩行能力が多面的に評価されたスコアとなっている。このような歩行スコアは、例えば高齢者の転倒リスク、フレイル、認知症等の健康リスクをより早期に捉えるための指標として用いることができる。特に、例えば歩行指標として左右の足の動きの規則性を含む場合には、身体機能の左右差が著しい脳卒中患者のリハビリテーションにおいて、下肢機能の回復の程度を示す指標としても用いることができる。
［その他の形態］
【００７３】
以上、本発明を特定の実施形態によって説明したが、本発明は上記した実施形態に限定されるものではない。
【００７４】
データ処理部２２は一体の構成である必要はない。例えば、データ処理部２２を構成する歩行データ取得部２２１、歩行指標算出部２２２、歩行指標スコア算出部２２３、サブスコア算出部２２４、歩行スコア算出部２２５、および歩行指標分類部２２６が別所に配置され、これらが互いにネットワークで接続されていてもよい。
【００７５】
また、データ処理部２２は全てが歩行評価装置２内に配置される必要はなく、一部が別所に配置されていてもよいし、一部が重複して別所に配置されていてもよい。例えば図１に示すように、歩行評価システム１０は、歩行評価装置２にネットワーク３で接続されたサーバ４をさらに備えることができ、歩行評価装置２が備える機能ブロックの一部を、サーバ４上にクラウド化して構成してもよい。同様に、サンプル歩行データ２３１も歩行評価装置２内に記録される必要はなく、サーバ４上の記憶装置に記録されていてもよい。
【００７６】
例えば、歩行指標算出部２２２、歩行指標スコア算出部２２３、サブスコア算出部２２４、歩行スコア算出部２２５および歩行指標分類部２２６をサーバ４上に構成し、サンプル歩行データ２３１をサーバ４に記録しておき、複数の歩行指標を関連グループに分類するための上記した実施形態のステップＳ１〜Ｓ４における処理を、サーバ４側で実行してもよい。この場合、歩行評価装置２のデータ処理部２２は、歩行データ取得部２２１、歩行指標算出部２２２、歩行指標スコア算出部２２３、サブスコア算出部２２４および歩行スコア算出部２２５を備え、被験者９の歩行スコアを算出するための上記した実施形態のステップＳ５〜Ｓ９における処理を実行する。
【００７７】
サブスコアを算出する方法は、上記実施形態の方法に限定されない。サブスコア算出部２２４は、関連グループ内に含まれている複数の歩行指標の歩行指標スコアの平均を算出することにより、サブスコアを算出してもよいし、各歩行指標スコアに重み付けを行い、重み付けされた各歩行指標スコアの和を算出することにより、サブスコアを算出してもよい。すなわち、サブスコア算出部２２４は、関連グループ毎にサブスコアを算出すればよく、各サブスコアは、関連グループ内に含まれている複数の歩行指標の歩行指標スコアに基づいて算出すればよい。
【００７８】
歩行スコアを算出する方法は、上記実施形態の方法に限定されない。歩行スコア算出部２２５は、各サブスコアに重み付けを行い、重み付けされた各サブスコアの総和を被験者９の歩行スコアとして算出してもよい。すなわち、表２に示す因子間に所定の重みを考慮して、歩行スコアを算出してもよい。
【００７９】
歩行指標は上記実施形態にて用いた１０種類に限定されず、以下に例示するような、歩行状態に関連する種々の歩行指標を用いることができる。なお、歩行評価に用いる歩行指標が、被験者９の足部から得られる歩行データにより算出される場合、センサ１は被験者９の体幹に取り付けられる必要は無い。
【００８０】
・ステップ時間の比率（Step time ratio）
左右の足の運びの時間に着目した全体的な対称性の指標であり、解析対象区間の左右ステップ時間の平均の比率である。
【００８１】
・ケイデンス
ストライド時間（ST）の平均値の逆数を求め、さらに２倍することにより算出する。ケイデンスはいわゆるピッチに相当し、単位時間当たりの歩数を表す。
【００８２】
・体幹加速度RMS
解析区間における体幹の加速度波形データの二乗平均根（RMS: root mean square）であり、三軸のそれぞれについて求める。以下の式により算出される。以下において言及するRMSも同じ計算式により算出される。
【数２】
2020151470.tif 000006
【００８３】
ここで、Ｘｉは、加速度データまたは角速度データであり、Ｎは解析対象となったデータの数（加速度データまたは加速度データの数）である。
【００８４】
・足部角速度RMS
解析区間における足部の角速度波形データの二乗平均根（RMS: root mean square）であり、三面のそれぞれについて求める。上記RMSを求める計算式により算出される。
【００８５】
・下肢RMS 比率
三面のそれぞれについて求められた、上記した足部角速度RMSの左右比である。
【００８６】
・交互相関（CC: Cross Correlation）
左右の足部の角速度波形の一致度を表す指標である。解析区間の両足部の波形のうち、一方を一歩時間分ずらし、他方の波形との相関係数を求めることにより算出する。両足のスイングが同じようなスピードで動いていると、交互相関の値は大きくなる。
【数３】
2020151470.tif 000007
【００８７】
式中、各指標の計算は以下の通りである。
【数４】
2020151470.tif 000008
【００８８】
ここで、ＸＲｉは右足の角速度であり、ＸＬｉは左足の角速度であり、ｍ＝０，±１，±２，…，±Ｎ−１である。この式で得られた数のベクトルの最大値がＣＣｍａｘとなる。
【００８９】
ＡＣＲ（ｍ）は、ＸＲｉの自己相関係数で次の式で示される。
【数５】
2020151470.tif 000009
【００９０】
同様に、ＡＣＬ（ｍ）は次の式で示される。
【数６】
2020151470.tif 000010
【００９１】
交互相関の詳細については、文献「Gouwanda, Darwin, and SMN Arosha Senanayake. "Identifying gait asymmetry using gyroscopes-A cross-correlation and Normalized Symmetry Index approach." Journal of Biomechanics 2011; 44.5: 972-978」を参照されたい。
【実施例１】
【００９２】
１．計測における事前準備
三軸の加速度および角速度を測定することが可能な小型センサを、被験者の左右の踵にそれぞれ取り付ける。さらに、もう１つの同じ小型センサを、被験者の重心位置に近似する第三腰椎棘突起付近に、伸縮性のベルトを巻き付けた上から取り付ける。約１５メートルの整地を歩行路とする。歩行路のうち、前後のそれぞれ約２．５メートルは加速路および減速路として、中間の約１０メートルの区間を計測区間とする。センサを取り付けた状態で被験者に静止立位をとってもらい、センサのキャリブレーションを行う。
【００９３】
２．歩行計測
前方の目線の高さに設置された目標物を見て歩くように被験者に指示を行い、「普段歩いているように、一番快適なスピードで歩いて下さい」と指示を行ったのち、練習試行を１度行う。その後本番試行を行い、加速度および角速度の歩行データの波形を取得する。加速度および角速度はサンプリング周期５００Ｈｚで取り込まれ、整流処理を行った後、歩行指標の計算が行われる。
【００９４】
３．歩行指標の分類
歩行指標を分類するためのサンプル歩行データを取得するために、上記した歩行計測を３７８人の高齢者に対して行った。３７８人分のサンプルの概要および用いた歩行指標の概要を表３に示す。
【表３】
2020151470.tif 000011
【００９５】
３７８人分のサンプルの歩行指標のデータについて統計処理を行い、それぞれの歩行指標について、歩行指標の値と歩行指標スコアとの対応関係を決定した。決定された対応関係を表４に示す。
【表４】
2020151470.tif 000012
【００９６】
３７８人分のサンプルの各歩行指標について歩行指標スコアを算出し、算出した歩行指標スコアに基づいて、歩行指標を複数の関連グループに分類した。分類はプロマックス回転を伴う因子分析により行った。得られた因子負荷行列は表２に示す通りである。
【００９７】
４．歩行スコアの妥当性に関する検証
３７８人のサンプルに対し、5-chair stand (5CS) テスト、タンデム立位、過去１年の転倒経験、転倒恐怖感（FoF: Fear of Falling）のそれぞれと、歩行スコアとの関連性を、ｔ検定を用いて検討し、歩行スコアの妥当性を評価した。結果を表５に示す。
【表５】
2020151470.tif 000013
【００９８】
表５において、ｐ値が０．０５以下の場合、相関性があることを意味する。すなわち、高齢者であることは、歩行スコア、歩調指標、ばらつき指標、およびなめらかさ指標のそれぞれと関連性があることを示している。同様に、筋力は、歩行スコア、歩調指標、およびなめらかさ指標のそれぞれと関連性があることを示している。バランス機能は、歩行スコア、規則性指標、歩調指標、およびばらつき指標のそれぞれと関連性があることを示している。転倒経験は、歩行スコアおよび規則性指標のそれぞれと関連性があることを示している。転倒恐怖感は、歩行スコア、規則性指標、歩調指標、およびなめらかさ指標のそれぞれと関連性があることを示している。
【００９９】
この結果より、歩行スコアは、加齢変化、筋力、バランス機能、転倒関連指標（転倒恐怖感、転倒歴）など、高齢者の健康状態を示す指標として妥当であることが示唆された。また、サブスコアは項目ごとに異なる結果となり、それぞれ高齢者の健康の異なる側面を示していると考えられた。このことから、歩行指標を分類して使用することの妥当性も確認された。
【実施例２】
【０１００】
実施例２では、実施例１とは別の地域に在住する高齢者に対して歩行スコアと年齢や転倒歴などとの関連を調べ、歩行スコアの妥当性を検証した。
【０１０１】
［対象および方法］
・対象
実施例２では、高齢者向け体力測定会に参加した地域高齢者１２５名のうち、包含基準と除外基準を満たし、質問紙への回答および歩行計測で加速度データにノイズが少なく、正常に計測できたものを対象とした。包含基準は、年齢が６５歳以上であって、歩行に際して歩行補助具を必要としないものとした。また、除外基準は、脳卒中の既往があるもの、またパーキンソン様症状を示すものとした。神経症状を有する高齢者では、歩行指標が健常者とは大きく異なる値となり、歩行指標の分布にゆがみが生じるためである。対象となった高齢者には事前に自記式のアンケートを記入してもらい、年齢、既往歴（高血圧症、糖尿病、呼吸器疾患、脳卒中、パーキンソン病）、過去一年間の転倒歴（はい・いいえ）、転倒恐怖感（はい・いいえ）について回答して貰った。また、身長、体重を自記式のアンケート回収後に計測した。
【０１０２】
・歩行計測
歩行計測には、三軸角速度センサと三軸加速度センサとを内蔵する多チャンネル慣性センサ（MVP-RF8-HC、マイクロストーン株式会社）を用いた。センサの取り付けにはサージカルテープを使用し、重心位置に近似する第三腰椎棘突起部および右踵にしっかりと固定した。センサが対象者の身体への動きを妨げないことを確認した後、静止立位での較正を実施した。続いて、対象者に歩行計測の手順を説明し歩行の練習をしてもらい、その後計測を実施した。対象者には普段通りの速度で１５メートルの歩行路を歩くように指示し、その際の踵（右）の矢状面の角速度と第三腰椎棘突起部における直行する３方向（鉛直、側方、前後）の加速度とを計測した。
【０１０３】
得られた波形データに、カットオフ周波数を２０Ｈｚとしたローパスフィルターを用いてフィルタリング処理を行った後、全波形データから最初の２ストライドを除いた安定歩行中の連続した５ストライドを抽出し、歩行指標を算出した。歩行指標は、ストライド時間（ST）、ストライド時間のバラつき（STV）、歩行の規則性（Regularity）、歩行のなめらかさ（Smoothness）を用いた。歩行指標の計算方法については上記した実施形態の方法に従った。
【０１０４】
続いて、各歩行指標の値は、上記した実施形態の結果を参照して１点刻みで０点から３点までの歩行スコアに変換した。歩行スコア算出のアルゴリズムは上記した実施形態に示す通りである。次に、歩行スコアは、事前に分類された４つの歩容カテゴリ（歩行の規則性（Regularity）、歩調（Pace）、歩行のばらつき（Variability）、歩行のなめらかさ（Smoothness））のうち、該当するカテゴリ内で合計し（規則性（スコア範囲：０〜９）、歩調（スコア範囲：０〜６）、ばらつき（スコア範囲：０〜６）、なめらかさ（スコア範囲：０〜９））、各サブスコアとした。最後にサブスコアを合計し歩行スコア（３０点満点）とした。
【０１０５】
・身体機能検査
身体機能測定には、高齢者の運動機能評価として広く用いられている五回立ち上がり（5-chair stand (5CS)）テストと片脚立ちテストを用いた。五回立ち上がりテストの方法には上記した実施形態を参考に次のような方法とした。被検者は椅子坐位をとり、起坐位から最大努力で五回立ち座りする動作を行ってもらった。検者は、動作の開始から終了までの時間をストップウォッチで計測した。測定は２回行い、小さい方の値を解析に用いた。片脚立ちテストは、両手を腰に添えた状態で、片脚立位を保持できる時間を計測した。両手のどちらかが腰から離れたり、支持脚とは反対の足が着地したりしたら計測を終了とした。最大時間は１２０秒とし、両側の計測を実施した。測定は対象者が申し出た得意なほうの脚から実施し、解析には１２０秒までの最大値を利用した。
【０１０６】
・統計解析
統計解析は次の手順で行った。まず歩行スコアとサブスコアの分布の正規性をグラフおよびShapiro-Wilkの正規性の検定を用いて確認した。次に歩行スコアおよびサブスコアと、年齢（前期高齢者および後期高齢者に分類）、過去一年間の転倒歴、及び転倒恐怖感との関連を対応のないt検定を用いて調べた。続いて、五回立ち上がりテスト及び片脚立ちテストの結果から、対象者を第一四分位以下に該当する低スコア群とそれ以外の高スコア群に分け、対応のないt検定を用いて群間比較を行った。統計学的検討には市販の統計ソフトウェア（JMP 14.0, SAS Institute Japan）を使用し、統計学的有意水準は５％とした。
【０１０７】
［結果］
最終的な解析対象となった参加者は８７名であった（包含基準：６５歳以上１１０名、除外基準：脳卒中の既往３名、質問紙の回答不完全１名、歩行計測の欠測１９名）。対象者の属性は表６に示す。図７の（Ａ）には歩行スコアのヒストグラムを示し（Ｂ）には正規分位点プロットを示す。図７の結果およびShapiro-Wilkの正規性の検定の結果から歩行スコアは正規分布していると判断された（Ｗ値＝０．９８９、Ｐ値＝０．７１４）。一方、各サブスコアは正規分布しているとは認められなかった（規則性：Ｗ値＝０．９２８、Ｐ値＜０．００１、歩調：Ｗ値＝０．９２１、Ｐ値＜０．００１、ばらつき：Ｗ値＝０．９２０、Ｐ値＜０．００１、なめらかさ：Ｗ値＝０．９５４、Ｐ値＝０．００４）。
【表６】
2020151470.tif 000014
【０１０８】
歩行指標および歩行スコアの結果を表７に示す。また、群間比較の結果を表８に示す。
【表７】
2020151470.tif 000015
【表８】
2020151470.tif 000016
【０１０９】
対応のないt検定の結果、歩行スコアは、前期高齢者と比較して後期高齢者が、転倒歴のない高齢者と比較して転倒歴のある高齢者が有意に低い値を示した（年齢：前期高齢者  ｖｓ  後期高齢者＝１７．８±６．１  ｖｓ  １３．８±６．０、ｔ値＝３．０、Ｐ値＝０．００４、過去一年間の転倒歴：非転倒者ｖｓ転倒者＝１２．９±５．４  ｖｓ  １７．０±６．４、ｔ値＝２．６、Ｐ値＝０．０１２）。しかし、第一四分位以下に該当する低スコア群とそれ以外の高スコア群で実施した運動機能の群間比較では有意な差は認められなかった（五回立ち上がりテストの第一四分位：７．７１（秒）、片脚立ちテストの第一四分位：２１．１４（秒））。サブスコアに分けて検討した結果では、規則性が前期高齢者と比較して後期高齢者が有意に低い値を示し（前期高齢者  ｖｓ  後期高齢者＝５．０±２．８  ｖｓ  ３．２±２．８、ｔ値＝２．８、Ｐ値＝０．００６）、なめらかさが転倒歴のない高齢者と比較して転倒歴のある高齢者で有意に低い値を示した（非転倒者ｖｓ転倒者＝５．６±２．３  ｖｓ  ３．５±２．２、ｔ値＝３．４、Ｐ値＝０．００１）。
【０１１０】
［考察］
実施例２では、実施例１とは別地域に在住する高齢者を対象に検証を行った。多チャンネル慣性センサを用いた歩行計測を実施し、センサ由来の歩行指標に基づいて我々が考案した歩行スコアを算出した。統計的手法によって歩行スコアと転倒リスク等との関連を調べ歩行スコアの妥当性について検討を行った。その結果、歩行スコアは、前期高齢者と比較して後期高齢者で、また非転倒者と比較して転倒者で有意に低値を示した。しかし、運動機能および転倒恐怖感とは有意な差は認められなかった。
【０１１１】
実施例２の結果から、歩行スコアは地域高齢者の健康リスクの指標として妥当だと判断した。実施例１では、歩行スコアは年齢、下肢筋力の低下、バランス機能の低下、過去一年間の転倒歴、転倒恐怖感と有意に関連していた。このうち、年齢および過去一年間の転倒歴については実施例２の結果と一致していた。このことから、加齢変化や転倒リスクに関して言えば、歩行スコアは地域高齢者の健康リスクの指標として妥当だと言える。このような結果が得られた理由として、歩行スコアが、それぞれが異なる特徴を有する複数の下位項目のサブスコアの合算値であることが挙げられる。すなわち、サブスコアの一つ一つが抽出できる異常歩行の特徴は小さいけれども、その差を合算することで差が明瞭になり、歩行スコアの異常歩行を抽出する能力が高まったからだと推測される。
【０１１２】
実施例２のサブスコアの結果では、規則性が年齢と、なめらかさが過去一年間の転倒歴と関連していた。しかし、実施例１では別のサブスコアが年齢および過去一年間の転倒歴と関連していた（年齢は、歩調、ばらつき、なめらかさが関連し、過去一年間の転倒歴は、規則性が関連していた）。このような相反する結果になった理由として、前述したように実施例２の対象者の運動機能が実施例１の対象者よりも高かったことが挙げられる。運動能力と転倒の関係は直線的ではなくＵ字形となることが報告されており、どのくらいの速度で歩けるかによって転倒リスクの高い歩容は異なると考えられている。これらのことから、実施例２において実施例１とは異なるサブスコアが年齢や転倒歴と有意な関連性を示したのは、実施例２の対象者の運動機能が高く、転倒に結び付く歩容が運動機能の低い高齢者の歩容とは異なっていたからではないか、と考えられる。また、実用的な視点からは、運動機能の異なる対象を扱ったいずれの実施例でもサブスコアが健康リスクを捉えていたことから、歩行スコアは、たとえ対象者の運動機能が違っても、いずれかのサブスコアが健康リスクを捉えている可能性が高く、地域高齢者の運動機能評価として幅広く用いることができるスコアだと考えられる。
【０１１３】
歩行スコアの平均値は１６．１で範囲は１から３０点であった。さらに、分布の正規性の検討の結果、歩行スコアの値は正規分布することが確認された。この結果は実施例１と同様であり、本スコアは連続変数として使用できる歩行指標であると考えられた。予防理学療法の介入研究においてメインアウトとして使用される指標のなかには、天井効果や床効果のために分布が歪んでしまい、連続変数として扱うことが難しい場合がある（例えばtimed up and go testのスコア）。その点、歩行スコアは、運動機能が異なっても正規分布を示しており、地域高齢者の運動機能の変化を捉えるような予防理学療法研究に適していると考えられる。
【０１１４】
［結論］
歩行スコアは、対象者の運動機能によらず地域高齢者における加齢変化や転倒リスクを示す指標として妥当である。また、歩行スコアは連続尺度の指標として使用が可能である。これらのことから、歩行スコアは予防理学療法研究への適用が期待される。
【符号の説明】
【０１１５】
１（１ａ〜１ｃ）  センサ
２      歩行評価装置
３      ネットワーク
４      サーバ
９      被験者
１０    歩行評価システム
２１    通信インタフェース部
２２    データ処理部
２２１  歩行データ取得部
２２２  歩行指標算出部
２２３  歩行指標スコア算出部
２２４  サブスコア算出部
２２５  歩行スコア算出部
２２６  歩行指標分類部
２３    補助記憶装置
２３１  サンプル歩行データ（第２の歩行データ）
２３２  歩行評価プログラム
２４    表示部
２４１  歩行スコア
２４２  サブスコア
２４３  関連グループ名 
【図１】
2020151470.tif 000017
【図２】
2020151470.tif 000018
【図３】
2020151470.tif 000019
【図４】
2020151470.tif 000020
【図５】
2020151470.tif 000021
【図６】
2020151470.tif 000022
【図７】
2020151470.tif 000023