JP6712909B2 - Cognitive change prediction device and program - Google Patents
Cognitive change prediction device and program Download PDFInfo
- Publication number
- JP6712909B2 JP6712909B2 JP2016118918A JP2016118918A JP6712909B2 JP 6712909 B2 JP6712909 B2 JP 6712909B2 JP 2016118918 A JP2016118918 A JP 2016118918A JP 2016118918 A JP2016118918 A JP 2016118918A JP 6712909 B2 JP6712909 B2 JP 6712909B2
- Authority
- JP
- Japan
- Prior art keywords
- change
- feature amount
- heartbeat
- cognitive ability
- amount
- 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
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
この発明は、“Working Memory”や“Selective Attention”等、ユーザの認知的な情報処理能力(以下認知能力と呼称する)の変化を予測するための装置およびプログラムに関する。 The present invention, "Working Memory" Ya "Selective Attention", etc., related equipment Contact and program for predicting changes in cognitive information processing capability of the user (hereinafter referred to as cognitive performance).
近年、例えば従業員の労務管理のために、就業中にユーザの認知能力を推定する技術が種々提案されている。例えば、非特許文献1には、Advanced Trail Making Test(以下ATMTと呼称する)のように、実行中の作業等を中断して数分程度のタスクを実行し、当該タスクの成績からユーザの認知能力を推定する手法が記載されている。また非特許文献2には、ユーザに心電センサを装着し、心電のR波ピーク間隔(以下RRIと呼称する)を測定し、このRRIから算出される複数の特徴量を用いて現時点での認知能力を推定する手法が記載されている。 In recent years, various techniques for estimating the cognitive ability of a user during work have been proposed, for example, for labor management of employees. For example, in Non-Patent Document 1, as in Advanced Trail Making Test (hereinafter referred to as ATMT), a work in progress is interrupted to execute a task for about several minutes, and the user's recognition based on the result of the task. Techniques for estimating capacity are described. Further, in Non-Patent Document 2, an electrocardiographic sensor is attached to a user, an R wave peak interval of the electrocardiogram (hereinafter referred to as RRI) is measured, and a plurality of feature amounts calculated from this RRI are used at the present time. A method for estimating the cognitive ability of humans is described.
ところが、非特許文献1、2に記載された技術は、何れも現時点でのユーザの認知能力を推定するものとなっている。このため、推定結果を労務管理に利用しようとした場合、認知能力が低下した時点でその人を休ませるといった対応は可能であるが、認知能力が低下する前に休ませるという対応は不可能だった。従って、生産性の低下は防げず、労務管理上改善が要望されていた。 However, the technologies described in Non-Patent Documents 1 and 2 all estimate the user's cognitive ability at the present time. Therefore, when trying to use the estimation results for labor management, it is possible to take a rest when the cognitive ability declines, but it is impossible to take a rest before the cognitive ability declines. It was Therefore, reduction in productivity cannot be prevented, and improvement in labor management has been demanded.
この発明は上記事情に着目してなされたもので、その目的とするところは、認知能力の変化を予測できるようにしてユーザの管理をより適切に行えるようにする認知能力変化予測装置およびプログラムを提供しようとするものである。 The present invention has been made in view of the above circumstances, the place of interest, and our cognitive change projection equipment so as to be able to predict changes in cognitive ability allows the user management better It is intended to provide a program.
上記課題を解決するためにこの発明の第1の態様は、予測対象ユーザの心拍間隔を表すデータを取得し、この取得された心拍間隔を表すデータをもとに、時間区間tにおける心拍変動の特徴量と、1つ前の時間区間t−1における心拍変動の特徴量をそれぞれ算出すると共に、当該算出された各心拍変動の特徴量をもとにその変化量を算出する。そして、上記算出された各心拍変動の特徴量およびその変化量と、予め記憶されている学習データとをもとに、上記予測対象ユーザの1つ後の時間区間t+1における認知能力の変化を、前記時間区間tにおける認知能力からの増減として予測し、その予測結果を表す情報を出力するようにしたものである。 In order to solve the above problems, a first aspect of the present invention acquires data representing a heartbeat interval of a prediction target user, and based on the acquired data representing a heartbeat interval, a heartbeat fluctuation of a time section t. The feature amount and the feature amount of the heartbeat variability in the immediately preceding time interval t-1 are calculated, and the amount of change is calculated based on the calculated feature amount of each heartbeat variability. Then, based on the calculated feature amount of each heartbeat variability and its change amount and the learning data stored in advance , the change in the cognitive ability of the prediction target user in the subsequent time interval t+1, The prediction is made as an increase or decrease from the cognitive ability in the time section t, and information indicating the prediction result is output.
この発明の第2の態様は、上記心拍変動の特徴量を算出する手段として、上記取得された心拍間隔を表すデータをもとに時間領域の特徴量を算出する手段と、上記取得された心拍間隔を表すデータをもとに周波数領域の特徴量を算出する手段とを備える。そして、上記時間領域の特徴量を算出する手段により、心拍間隔の平均値と、心拍間隔の標準偏差(SDNN)と、隣接した心拍間隔の差が50msecを超えた比率(pNN50)と、隣接した心拍間隔の差の二乗平均平方和(RMSSD)とのうちの少なくとも1つを算出し、上記周波数領域の特徴量を算出する手段により、心拍間隔をスペクトル解析した際の低周波成分(0.04~0.15Hz)の総パワー(LF)と、心拍間隔をスペクトル解析した際の高周波成分(0.15~0.40Hz)の総パワー(HF)と、上記低周波成分の総パワーLFと上記高周波成分の総パワーHFとの比率(LF/HF)とのうちの少なくとも1つを算出するようにしたものである。 A second aspect of the present invention is, as means for calculating the feature amount of the heartbeat variability, means for calculating a feature amount in a time domain based on the acquired data indicating the heartbeat interval, and the obtained heartbeat. And means for calculating a feature amount in the frequency domain based on the data representing the interval. Then, by the means for calculating the feature amount in the time domain, the average value of the heartbeat intervals, the standard deviation of the heartbeat intervals (SDNN), and the ratio (pNN50) in which the difference between adjacent heartbeat intervals exceeds 50 msec, the The low-frequency component (0.04 to 0.15) when the heartbeat interval is spectrally analyzed by means for calculating at least one of the root mean square difference of the heartbeat intervals (RMSSD) and the feature amount in the frequency domain is calculated. Hz) total power (LF), total power (HF) of high frequency components (0.15 to 0.40 Hz) when spectrum analysis of heartbeat interval, total power LF of the above low frequency components and total power HF of the above high frequency components And at least one of the ratio (LF/HF) and is calculated.
この発明の第3の態様は、上記心拍変動の特徴量の変化量を算出する手段により、上記時間区間tにおける心拍変動の特徴量と、1つ前の時間区間t−1における心拍変動の特徴量との間の変化量または変化率と、上記時間区間tにおける心拍変動の特徴量と、1つ前の時間区間t−1における心拍変動の特徴量との間の相関係数の変化量または変化率とのうちの何れかを算出するようにしたものである。 A third aspect of the invention, the means for calculating a feature amount of variation of the above SL HRV, the feature amount of heart rate variability in the time interval t, the heart rate variability in the previous time interval t-1 a change amount or rate of change between the feature quantity, the amount of change in the correlation coefficient between the feature quantity of heart rate variability in the time interval t, the feature amount of heart rate variability in the previous time interval t-1 Alternatively, one of the change rate and the change rate is calculated.
この発明の第4の態様は、上記心拍変動の特徴量の変化量を算出する際に、上記時間区間tにおける心拍変動の特徴量と、1つ前の時間区間t−1における心拍変動の特徴量に対し、それぞれ主成分の分析を行い、当該主成分の分析結果に基づいて、上記主成分の変化を示すT2 統計量と、次元圧縮の残差を示すQ統計量とをそれぞれ算出するようにしたものである。 According to a fourth aspect of the present invention, when calculating the change amount of the feature amount of the heartbeat variability, the feature amount of the heartbeat variability in the time period t and the feature of the heartbeat variability in the immediately preceding time period t-1. The principal component is analyzed for each quantity, and the T 2 statistic indicating the change of the principal component and the Q statistic indicating the residual of dimension compression are calculated based on the analysis result of the principal component. It was done like this.
この発明の第6の態様は、上記認知能力の変化を予測する手段として、上記学習データをもとに認知能力の変化を表す分類器を作成し、上記算出された心拍変動の特徴量とその変化量を上記分類器に入力することで、上記認知能力の変化を予測する手段と、上記学習データをもとに回帰モデルを作成し、上記算出された心拍変動の特徴量とその変化量を上記回帰モデルに入力することで、上記認知能力の変化量を予測する手段とのうちの何れか一方を備えるようにしたものである。 According to a sixth aspect of the present invention, as a means for predicting a change in the cognitive ability, a classifier representing a change in the cognitive ability is created based on the learning data, and the calculated feature amount of the heartbeat variability and its By inputting the change amount into the classifier, a regression model is created based on the means for predicting the change in the cognitive ability and the learning data, and the calculated feature amount of the heartbeat variability and the change amount thereof. One of the means for predicting the change amount of the cognitive ability by inputting into the regression model is provided.
この発明の第1の態様によれば、予測対象ユーザの例えば現在の時間区間tより1つ後の時間区間t+1における認知能力の変化が、時間区間tにおける認知能力からの増減として予測されてその結果が表示される。このため、例えば認知能力が低下する以前に予測対象ユーザが自主的に休息を取ったり、管理者が予測対象ユーザに対し休息を与える等の対策を適切なタイミングで講じることが可能となる。その結果、タスクを実行中のユーザの作業効率や品質を高く維持することが可能となる。 According to the first aspect of the present invention, a change in the cognitive ability of the prediction target user, for example, in the time section t+1, which is one time after the current time section t, is predicted as an increase or decrease from the cognitive ability in the time section t. The result is displayed. Therefore, for example, it is possible for the prediction target user to take a rest before the cognitive ability deteriorates, or the administrator can take measures such as giving a rest to the prediction target user at an appropriate timing. As a result, it is possible to maintain high work efficiency and quality of the user who is executing the task.
この発明の第2の態様によれば、心拍変動の特徴量として、心拍間隔データをもとに時間領域の特徴量と周波数領域の特徴量がそれぞれ算出される。このため、ユーザの認知能力と関連性を持つ特徴量を確実に得ることが可能となる。 According to the second aspect of the present invention, the feature amount in the time domain and the feature amount in the frequency domain are each calculated as the feature amount of the heartbeat fluctuation based on the heartbeat interval data. Therefore, it is possible to reliably obtain the feature amount that is related to the cognitive ability of the user.
この発明の第3の態様によれば、心拍変動の特徴量の変化量として、時間区間tにおける心拍変動の特徴量と、1つ前の時間区間t−1における心拍変動特徴量との間の変化量または変化率と、上記各心拍変動特徴量間の相関係数の変化量または変化率とのうちの何れかが算出される。このため、心拍変動の特徴量の変化量を時間区間単位で比較的簡易に算出することができる。 According to a third aspect of the present invention, as the feature quantity variation of heart rate variability, and characteristics of heart rate variability in the time interval t, between the HRV feature value of the previous in time interval t-1 a change amount or rate of change, one of a change amount or rate of change of the correlation coefficient between each HRV feature amount is calculated. For this reason, it is possible to relatively easily calculate the amount of change in the characteristic amount of the heartbeat variability on a time interval basis.
この発明の第4の態様によれば、心拍変動の特徴量の変化量として、心拍変動の特徴量の主成分の変化を示すT2 統計量と、次元圧縮の残差を示すQ統計量が算出される。このため、心拍変動の特徴量の変化量をより的確に算出することができる。 According to the fourth aspect of the present invention, as the amount of change in the feature amount of heartbeat variability, there are T 2 statistic indicating the change of the main component of the feature amount of heartbeat variability and Q statistic indicating the residual of dimensional compression. It is calculated. Therefore, the amount of change in the feature amount of heartbeat variability can be calculated more accurately.
この発明の第6の態様は、ユーザの認知能力の変化を、認知能力の変化の有無として予測したり、認知能力の変化量として予測することができる。 In the sixth aspect of the present invention, the change in the cognitive ability of the user can be predicted as the presence or absence of the change in the cognitive ability or the change amount of the cognitive ability.
すなわちこの発明の各態様によれば、認知能力の変化を予測できるようにしてユーザの管理をより適切に行えるようにした認知能力変化予測装置およびプログラムを提供することができる。 That is, according to each aspect of the invention, it is possible to provide a cognitive change projection equipment Contact and programs so as to allow predict changes in cognitive ability allows the user of the management more properly.
[原理]
図1はこの発明の原理説明に使用するための図で、1は予測処理装置、2はデータ入力手段、3は結果出力手段、4は学習データ蓄積手段をそれぞれ示している。
データ入力手段2は、対象ユーザの心拍間隔データを予測処理装置1に入力するもので、例えばウェアラブル心拍センサ、脈波センサ、或いは心拍間隔データが蓄積されたデータベースからなる。結果出力手段3は、上記予測処理装置1により得られた認知能力の予測結果を出力するもので、例えばスマートフォンで使用されるアプリケーション・プログラムからなる。
[principle]
FIG. 1 is a diagram for use in explaining the principle of the present invention, in which 1 is a prediction processing device, 2 is data input means, 3 is result output means, and 4 is learning data storage means.
The data input unit 2 inputs the heartbeat interval data of the target user to the prediction processing device 1, and is composed of, for example, a wearable heartbeat sensor, a pulse wave sensor, or a database in which heartbeat interval data is accumulated. The result output means 3 outputs the prediction result of the cognitive ability obtained by the prediction processing device 1, and is composed of, for example, an application program used in a smartphone.
学習データ蓄積手段4には、予め測定された、予測対象となる複数のユーザの、作業中における心拍変動を示す特徴量と、当該ユーザの作業中における情報処理能力に関する情報とを、時刻が同じもの同士で対応付けた学習データが記憶される。より具体的には、ユーザの認知能力変化を示す値(以下認知能力変化値と呼称する)、HRV特徴量、およびHRV特徴量の変化量を表す情報が、学習データとして記憶される。 In the learning data accumulating means 4, the pre-measured characteristic amounts of the plurality of prediction target users, which indicate heart rate variability during work, and the information regarding the information processing ability of the users during work are stored at the same time. Learning data associated with each other is stored. More specifically, a value indicating a change in cognitive ability of the user (hereinafter referred to as a cognitive ability change value), an HRV feature amount, and information indicating an amount of change in the HRV feature amount are stored as learning data.
認知能力変化値としては、例えばATMT(Advanced Trail Making Test)の成績やその変化量が挙げられる。HRV特徴量としては、例えばRRIの平均値、RRIの標準偏差、またはRRIをスペクトル解析した際の高周波帯(0.15Hz〜0.40Hz)および低周波帯(0.04Hz〜0.15Hz)それぞれのパワースペクトルの比であるLF/HFが挙げられる。HRV特徴量の変動量としては、複数のHRV特徴量に対し主成分分析を適用した後に算出される、主成分得点の変動を示すT2統計量や特徴量間の相関関係の変動を示すQ統計量が挙げられる。 Examples of the cognitive ability change value include, for example, an ATMT (Advanced Trail Making Test) performance and its change amount. As the HRV feature amount, for example, the average value of RRI, the standard deviation of RRI, or the power spectrum of each high frequency band (0.15 Hz to 0.40 Hz) and low frequency band (0.04 Hz to 0.15 Hz) when spectrum analysis of RRI is performed. The ratio is LF/HF. As the variation amount of the HRV feature amount, T 2 statistic indicating the variation of the principal component score calculated after applying the principal component analysis to a plurality of HRV feature amounts and Q indicating the variation of the correlation between the feature amounts. There are statistics.
予測処理装置1は、この発明を実施するための処理機能として、特徴量算出手段11と、特徴量変化算出手段12と、認知能力変化予測手段13とを有する。
特徴量算出手段11は、上記データ入力手段2から入力された心拍間隔データから所定の時間区間ごとにHRV特徴量を算出する。特徴量変化算出手段12は、上記特徴量算出手段11により時間区間ごとに算出されたHRV特徴量をもとに、HRV特徴量の変化量(以下HRV変化量と呼称する)を算出する。認知能力変化予測手段13は、上記特徴量算出手段11により算出されたHRV特徴量および上記特徴量変化算出手段12により算出されたHRV変化量をもとに、上記学習データ蓄積手段4に記憶された学習データを参照して、認知能力の変化を予測する。
The prediction processing device 1 has a feature amount calculation unit 11, a feature amount change calculation unit 12, and a cognitive ability change prediction unit 13 as processing functions for implementing the present invention.
The feature amount calculation means 11 calculates the HRV feature amount for each predetermined time interval from the heartbeat interval data input from the data input means 2. The feature amount change calculation means 12 calculates the change amount of the HRV feature amount (hereinafter referred to as HRV change amount) based on the HRV feature amount calculated for each time section by the feature amount calculation means 11. The cognitive ability change prediction unit 13 is stored in the learning data storage unit 4 based on the HRV feature amount calculated by the feature amount calculation unit 11 and the HRV change amount calculated by the feature amount change calculation unit 12. The change in cognitive ability is predicted by referring to the learned data.
以上のような構成であるから、予測処理装置1は以下のように動作する。図2はその処理手順と処理内容を示すフローチャートである。
予測処理装置1は、先ずステップS11において、データ入力手段2から対象ユーザの心拍間隔データを時系列で受け取り、当該心拍間隔データから所定の時間区間ごとにHRV特徴量を算出する。続いてステップS12において、上記ステップS11により時間区間ごとに算出されたHRV特徴量をもとにHRV変化量を算出する。次にステップS13において、上記ステップS11により算出されたHRV特徴量およびステップS12により算出されたHRV変化量をもとに、学習データ蓄積手段4に記憶された学習データを参照して対象ユーザの認知能力の変化量を予測する。そして、上記予測した認知能力の変化量を表す情報を結果出力手段3へ出力し、対象ユーザまたはその管理者等に提示する。
With the above configuration, the prediction processing device 1 operates as follows. FIG. 2 is a flowchart showing the processing procedure and processing contents.
In step S11, the prediction processing device 1 first receives the heartbeat interval data of the target user from the data input unit 2 in time series, and calculates the HRV feature amount for each predetermined time interval from the heartbeat interval data. Then, in step S12, the HRV change amount is calculated based on the HRV feature amount calculated for each time section in step S11. Next, in step S13, based on the HRV feature amount calculated in step S11 and the HRV change amount calculated in step S12, referring to the learning data stored in the learning data storage means 4, the target user's recognition Predict the amount of change in ability. Then, the information indicating the predicted amount of change in cognitive ability is output to the result output means 3 and presented to the target user or the manager thereof.
従ってこの発明によれば、対象ユーザの認知能力の変化量の予測結果が提示されるので、例えばユーザの認知能力が低下する前に当該ユーザに対し休息を与えるといった、より適切な労務管理を行うことが可能となる。 Therefore, according to the present invention, since the prediction result of the change amount of the cognitive ability of the target user is presented, for example, more appropriate labor management such as giving a rest to the user before the cognitive ability of the user is deteriorated is performed. It becomes possible.
以下、図面を参照してこの発明に係わる実施形態を説明する。
[第1の実施形態]
(構成)
図3は、この発明の第1の実施形態に係る認知能力変化予測装置を含むシステムの構成を示すブロック図である。
このシステムは、認知能力変化予測装置10Aを中核とし、加えて脈波センサ21,22と、ATMT入力部23と、表示部30と、学習データベース(学習DB)40とを備えている。
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
[First Embodiment]
(Constitution)
FIG. 3 is a block diagram showing the configuration of a system including the cognitive ability change prediction device according to the first embodiment of the present invention.
This system has the cognitive ability change prediction device 10A as a core, and additionally includes pulse wave sensors 21 and 22, an ATMT input unit 23, a display unit 30, and a learning database (learning DB) 40.
脈波センサ21は、学習データ測定対象ユーザに装着される。そして、学習データ測定対象ユーザの脈波を計測し、当該脈波の計測データをもとに心拍間隔を表すデータを算出して、その算出データを例えば無線インタフェースを介して認知能力変化予測装置10Aへ送信する機能を備えている。 The pulse wave sensor 21 is attached to the user who measures the learning data. Then, the pulse wave of the learning data measurement target user is measured, data representing a heartbeat interval is calculated based on the measured data of the pulse wave, and the calculated data is used, for example, via a wireless interface to predict the cognitive ability change device 10A. It has a function to send to.
脈波センサ22は、認知能力変化の予測対象ユーザに装着される。そして、認知能力変化の予測対象ユーザの脈波を計測し、当該脈波の計測データをもとに心拍間隔を表すデータを算出して、その算出データを例えば無線インタフェースを介して認知能力変化予測装置10Aへ送信する機能を備える。 The pulse wave sensor 22 is attached to a user who is a target of prediction of a change in cognitive ability. Then, the pulse wave of the user who is the target of prediction of the change in cognitive ability is measured, data representing the heartbeat interval is calculated based on the measured data of the pulse wave, and the calculated data is predicted through, for example, a wireless interface. It has a function of transmitting to the device 10A.
上記各脈波センサ21,22が使用する無線インタフェースとしては、例えばBluetooth(登録商標)等の近距離無線データ通信規格を採用した無線インタフェースが使用される。なお、脈波センサ21,22は個別に用意せず、1台の脈波センサを学習データ作成用と認知能力変化予測用として共用するようにしてもよい。 As the wireless interface used by each of the pulse wave sensors 21 and 22, for example, a wireless interface adopting a short-range wireless data communication standard such as Bluetooth (registered trademark) is used. Note that the pulse wave sensors 21 and 22 may not be individually prepared, and one pulse wave sensor may be shared for learning data creation and cognitive ability change prediction.
ATMT入力部23は、ATMTタスクの実施結果を表すデータを入力するためのインタフェースを備えた携帯端末(例えばノート型のパーソナルコンピュータやタブレット型端末、スマートフォン)からなる。そして、上記入力されたATMTタスクの実施結果を表すデータを、例えば無線インタフェースを介して認知能力変化予測装置10Aに向け送信する。 The ATMT input unit 23 includes a mobile terminal (for example, a notebook personal computer, a tablet terminal, or a smart phone) having an interface for inputting data representing the execution result of the ATMT task. Then, the data representing the execution result of the input ATMT task is transmitted to the cognitive performance change prediction device 10A via, for example, a wireless interface.
表示部30は、認知能力変化予測結果を表示するためのアプリケーション・プログラムを備えた携帯端末(例えばウェアラブル端末やスマートフォン、タブレット型端末)からなる。そして、認知能力変化予測装置10Aから無線インタフェースを介して送られた認知能力変化の予測結果を表す情報を受信し、当該情報をディスプレイに表示する。 The display unit 30 includes a mobile terminal (for example, a wearable terminal, a smartphone, or a tablet terminal) that includes an application program for displaying the cognitive ability change prediction result. Then, the information representing the prediction result of the cognitive ability change sent from the cognitive ability change prediction device 10A via the wireless interface is received, and the information is displayed on the display.
学習DB40は、例えばWeb上のサーバコンピュータまたはクラウドコンピュータに設けられ、認知能力変化予測装置10Aにより作成された学習データをネットワークを介して受信し記憶する。また、認知能力変化予測装置10Aから送られる参照要求に応じて、該当する学習データを要求元の認知能力変化予測装置10Aへ送信する。なお、学習DB40は、認知能力変化予測装置10A内に設けてもよい。 The learning DB 40 is provided in, for example, a server computer or a cloud computer on the Web, and receives the learning data created by the cognitive ability change prediction device 10A via the network and stores it. In addition, in response to the reference request sent from the cognitive ability change prediction device 10A, the corresponding learning data is transmitted to the request source cognitive performance change prediction device 10A. The learning DB 40 may be provided in the cognitive ability change prediction device 10A.
認知能力変化予測装置10Aは、CPU(Central Processing Unit)等のプロセッサ、記憶部、および無線インタフェース等の各種外部インタフェースを備えるパーソナルコンピュータまたはサーバコンピュータからなる。そして、本実施形態を実施するために必要な処理機能として、特徴量算出部101,104と、特徴量変化算出部102,105と、学習データ作成部103と、認知能力変化予測部106とを有している。これらの処理機能は何れも、上記記憶部に格納されたプログラムを上記プロセッサに実行させることにより実現される。 The cognitive ability change prediction device 10A includes a personal computer or a server computer including a processor such as a CPU (Central Processing Unit), a storage unit, and various external interfaces such as a wireless interface. Then, as the processing functions necessary for implementing the present embodiment, the feature amount calculation units 101 and 104, the feature amount change calculation units 102 and 105, the learning data creation unit 103, and the cognitive ability change prediction unit 106 are provided. Have Each of these processing functions is realized by causing the processor to execute the program stored in the storage unit.
上記各処理機能のうち、特徴量算出部101、特徴量変化算出部102および学習データ作成部103は、学習データを作成するために用いられる。 Among the above processing functions, the feature amount calculation unit 101, the feature amount change calculation unit 102, and the learning data creation unit 103 are used to create learning data.
特徴量算出部101は、上記脈波センサ21から時系列で送信される、学習データ測定対象ユーザの心拍間隔データを受信し、この受信された心拍間隔データをもとに所定の時間区間ごとにHRV特徴量を算出する。 The feature amount calculation unit 101 receives the heartbeat interval data of the learning data measurement target user, which is transmitted in time series from the pulse wave sensor 21, and based on the received heartbeat interval data, for each predetermined time interval. HRV feature quantity is calculated.
特徴量変化算出部102は、上記算出されたHRV特徴量をもとに、時間的に隣接する区間のT2統計量およびQ統計量を、HRV特徴量の変化量を表すデータとして算出する。そして、この算出されたT2統計量およびQ統計量を、上記特徴量算出部101により算出されたHRV特徴量と共に学習データ作成部103に渡す。なお、上記特徴量変化算出部102では、上記異なる複数の時間区間において算出された各心拍変動の特徴量間の変化量または変化率を算出する以外に、上記各心拍変動の特徴量間の相関係数の変化量または変化率を算出するようにしてもよい。 Based on the calculated HRV feature amount, the feature amount change calculation unit 102 calculates the T 2 statistic and the Q statistic of the temporally adjacent sections as data representing the change amount of the HRV feature amount. Then, the calculated T 2 statistic and Q statistic are passed to the learning data creation unit 103 together with the HRV feature amount calculated by the feature amount calculation unit 101. In addition, the feature amount change calculation unit 102 calculates the amount of change or rate of change between the feature amounts of each heartbeat variability calculated in the plurality of different time intervals, and also calculates the phase between the feature amounts of each heartbeat variability. You may make it calculate the variation|change_quantity or rate of change of the number of relations.
学習データ作成部103は、上記特徴量変化算出部102から渡された、HRV特徴量、T2統計量およびQ統計量を説明変数、ATMT入力部23において学習データ測定対象ユーザが入力したATMTエラー率の増減を表すデータを被説明変数とする学習データを作成し、この作成した学習データを学習DB40へ送信し記憶させる。 The learning data creation unit 103 uses the HRV feature amount, the T 2 statistic amount, and the Q statistic amount, which are passed from the feature amount change calculation unit 102, as explanatory variables, and the ATMT error input by the learning data measurement target user in the ATMT input unit 23. Learning data in which the data representing the increase/decrease in the rate is the explained variable is created, and the created learning data is transmitted to and stored in the learning DB 40.
一方、上記各処理機能のうち、特徴量算出部104、特徴量変化算出部105および認知能力変化予測部106は、予測対象ユーザの認知能力の変化を予測するために使用される。 On the other hand, among the above processing functions, the feature amount calculation unit 104, the feature amount change calculation unit 105, and the cognitive ability change prediction unit 106 are used to predict the change in the cognitive ability of the prediction target user.
特徴量算出部104は、脈波センサ22から時系列で送信される、予測対象ユーザの心拍間隔データを受信し、この受信された心拍間隔データをもとに所定の時間区間ごとにHRV特徴量を算出する。 The feature amount calculation unit 104 receives the heartbeat interval data of the prediction target user, which is transmitted in time series from the pulse wave sensor 22, and based on the received heartbeat interval data, the HRV feature amount for each predetermined time section. To calculate.
特徴量変化算出部105は、上記算出されたHRV特徴量をもとに、時間的に隣接する区間のT2統計量およびQ統計量を、HRV特徴量の変化量として算出する。そして、この算出されたT2統計量およびQ統計量を、上記特徴量算出部104により算出されたHRV特徴量と共に認知能力変化予測部106へ出力する。なお、上記特徴量変化算出部105では、上記異なる複数の時間区間において算出された各心拍変動の特徴量間の変化量または変化率を算出する以外に、上記各心拍変動の特徴量間の相関係数の変化量または変化率を算出するようにしてもよい。 The feature amount change calculation unit 105 calculates the T 2 statistic and the Q statistic of the temporally adjacent sections as the change amount of the HRV feature amount based on the calculated HRV feature amount. Then, the calculated T 2 statistic and Q statistic are output to the cognitive ability change prediction unit 106 together with the HRV feature amount calculated by the feature amount calculation unit 104. In addition, the feature amount change calculation unit 105 calculates the amount of change or the rate of change between the feature amounts of each heartbeat variability calculated in the plurality of different time intervals, and also calculates the phase between the feature amounts of each heartbeat variability. You may make it calculate the variation|change_quantity or rate of change of the number of relations.
認知能力変化予測部106は、上記特徴量変化算出部105によりHRV変化量として算出されたT2統計量およびQ統計量と、学習DB40に記憶された学習データとをもとに、対象ユーザの認知能力の変化を予測する。そして、当該予測結果の表示データを生成し、表示部30へ送信する。 The cognitive ability change prediction unit 106 uses the T 2 statistic and the Q statistic calculated as the HRV change amount by the feature amount change calculation unit 105, and the learning data stored in the learning DB 40 to target the target user. Predict changes in cognitive ability. Then, the display data of the prediction result is generated and transmitted to the display unit 30.
(動作)
次に、以上のように構成された認知能力変化予測装置10Aの動作を説明する。ここでは、予測する認知能力変化の増減を「増えたか」、「減ったか」のみの2値で予測する場合について説明を行う。
(motion)
Next, the operation of the cognitive ability change prediction device 10A configured as described above will be described. Here, a case will be described in which the predicted increase/decrease in the change in cognitive ability is predicted using only two values, "whether it has increased" and "has decreased".
(1)学習データの作成
予測対象ユーザの認知能力変化の予測処理に先立ち、認知能力変化予測装置10Aは学習データの作成処理を以下のように実行する。図4はその処理手順と処理内容を示すフローチャートである。
(1) Creation of Learning Data Prior to the process of predicting the change in cognitive ability of the prediction target user, the cognitive performance change prediction device 10A executes the process of creating learning data as follows. FIG. 4 is a flowchart showing the processing procedure and processing contents.
学習データの測定対象ユーザには脈波センサ21が装着される。脈波センサ21は測定対象ユーザの脈波を計測し、計測された脈波の波形データをもとに心拍間隔データを算出する。そして、算出された心拍間隔データを無線インタフェースを介して認知能力変化予測装置10Aに向け送信する。 The pulse wave sensor 21 is attached to the user whose learning data is to be measured. The pulse wave sensor 21 measures the pulse wave of the measurement target user, and calculates the heartbeat interval data based on the waveform data of the measured pulse wave. Then, the calculated heartbeat interval data is transmitted to the cognitive performance change prediction device 10A via the wireless interface.
また、上記心拍間隔データの計測動作と並行して学習データの測定対象ユーザは、10分間隔でATMTのタスクを実施し、ATMTのエラー率と実施時刻を自身の認知能力を表すデータとしてATMT入力部23において入力する。ATMT入力部23は、上記ATMTエラー率と実施時刻が入力されるごとに、これらを無線インタフェースを介して認知能力変化予測装置10Aに向け送信する。 Further, in parallel with the operation of measuring the heartbeat interval data, the user to be measured with the learning data performs the task of ATMT at intervals of 10 minutes, and inputs the error rate and execution time of ATMT as data representing his/her own cognitive ability. Input in the section 23. The ATMT input unit 23 transmits each time the ATMT error rate and the implementation time are input to the cognitive performance change prediction device 10A via the wireless interface.
認知能力変化予測装置10Aは、先ず特徴量算出部101の制御の下、ステップS21により上記脈波センサ21から時系列的に送信される心拍間隔データを受信して取り込み、図示しない記憶部に順次保存する。またステップS22により、上記ATMT入力部23から送信されたATMTエラー率および実施時刻を受信し、この受信したATMTエラー率および実施時刻を学習データ作成部103に、また実施時刻を特徴量算出部101に入力する。 Under the control of the feature amount calculation unit 101, the cognitive performance change prediction apparatus 10A first receives and takes in the heartbeat interval data transmitted in time series from the pulse wave sensor 21 in step S21, and sequentially stores the data in a storage unit (not shown). save. In step S22, the ATMT error rate and the execution time transmitted from the ATMT input unit 23 are received, the received ATMT error rate and the execution time are received by the learning data creation unit 103, and the execution time is received by the feature amount calculation unit 101. To enter.
続いて特徴量算出部101は、ステップS23において、上記記憶部に保存されている心拍間隔データを、上記ATMTタスクの実施時刻に同期して10分間隔(区間t (t =1,2,…,n)に分割する。そして、上記区間t ごとに、5分の時間窓を6秒のスライド幅でスライドさせながらそれぞれHRV特徴量を算出し、算出したHRV特徴量を特徴量変化算出部102に渡す。 Then, in step S23, the feature amount calculation unit 101 synchronizes the heartbeat interval data stored in the storage unit with a time interval of 10 minutes (interval t (t=1, 2,... , N), and the HRV feature amount is calculated while sliding the time window of 5 minutes with the slide width of 6 seconds for each section t, and the calculated HRV feature amount is calculated by the feature amount change calculation unit 102. Pass to.
算出されるHRV特徴量には、時間領域の特徴量と周波数領域の特徴量がある。
時間領域の特徴量には、例えば
(1) 心拍間隔の平均値
(2) 心拍間隔の標準偏差(SDNN)
(3) 隣接した心拍間隔の差が50msecを超えた比率(pNN50)
(4) 隣接した心拍間隔の差の二乗平均平方和(RMSSD)
がある。
The calculated HRV feature amount includes a time domain feature amount and a frequency domain feature amount.
The feature quantity in the time domain is, for example,
(1) Average heartbeat interval
(2) Standard deviation of heartbeat interval (SDNN)
(3) Ratio of adjacent heartbeat intervals exceeding 50 msec (pNN50)
(4) Sum of root mean squares of differences between adjacent heartbeat intervals (RMSSD)
There is.
これに対し周波数領域の特徴量には、例えば
(5) 心拍間隔をスペクトル解析した際の低周波成分(0.04~0.15Hz)の総パワー(LF)
(6) 心拍間隔をスペクトル解析した際の高周波成分(0.15~0.40Hz)の総パワー(HF)
(7) 低周波成分の総パワーLFと、高周波成分の総パワーHFとの比率(LF/HF)
がある。
On the other hand, for the feature quantity in the frequency domain, for example,
(5) Total power (LF) of low frequency components (0.04 to 0.15Hz) when spectrum analysis of heartbeat interval
(6) Total power (HF) of high frequency components (0.15 to 0.40Hz) when spectrum analysis of heartbeat interval
(7) Ratio of total power LF of low frequency components to total power HF of high frequency components (LF/HF)
There is.
特徴量算出部101は、HRV特徴量として、例えば上記時間領域の特徴量のうちのSDNNと、周波数領域の特徴量のうちのLF/HFを算出するが、時間領域の特徴量および周波数領域の特徴量のそれぞれにおいて他の特徴量を算出してもよく、また複数の特徴量を同時に算出してもよい。 The feature amount calculation unit 101 calculates, as the HRV feature amount, for example, SDNN of the time domain feature amounts and LF/HF of the frequency domain feature amounts. Other characteristic amounts may be calculated for each of the characteristic amounts, or a plurality of characteristic amounts may be calculated simultaneously.
認知能力変化予測装置10Aは、次に特徴量変化算出部102の制御の下、ステップS24において、上記特徴量算出部101から渡されたHRV特徴量をもとに、区間t (t>1)ごとにHRV変化量を算出する。HRV変化量はT2 統計量およびQ統計量からなる。そして、この算出されたHRV変化量を、上記特徴量算出部101により算出されたHRV特徴量と共に、学習データ作成部103に渡す。 Next, under the control of the feature amount change calculation unit 102, the cognitive ability change prediction device 10A, in step S24, based on the HRV feature amount passed from the feature amount calculation unit 101, a section t (t>1). HRV change amount is calculated for each. The HRV change amount consists of T 2 statistic and Q statistic. Then, the calculated HRV change amount is passed to the learning data creation unit 103 together with the HRV feature amount calculated by the feature amount calculation unit 101.
上記T2 統計量およびQ統計量の算出方法は以下の通りである。
Step1 ;区間t-1におけるHRV特徴量のベクトルX∈RA×B (A=窓数、B=特徴量数)に主成分数Rとして主成分分析を適用し、主成分得点ベクトルS∈RA×Bと主成分の負荷量行列VR を算出する。
The calculation method of the above T 2 statistic and Q statistic is as follows.
Step1; Principal component analysis is applied to the vector of HRV feature quantity X ∈ RA×B (A=number of windows, B=number of feature quantity) as the number of principal component R in the section t−1, and the principal component score vector S∈R The load matrix V R of A×B and the main component is calculated.
Step2 ;区間t におけるHRV特徴量のベクトルX′∈RA×Bに対し、Step1 により算出された負荷量行列VR を用いて主成分得点ベクトルS′∈RA×Bを算出する。 Step2; For the HRV feature vector X'ε RA x B in the section t, the principal component score vector S'ε RA x B is calculated using the load matrix V R calculated in step 1.
Step3 ;主成分得点ベクトルS′における各主成分得点Srと、主成分得点ベクトルS′における各主成分得点S′rより、区間t のT2統計量を、
Step4 ;主成分得点ベクトルS′および主成分の負荷量行列VRより、上記HRV特徴量のベクトルX′を主成分より再構成した再構成特徴量ベクトルX`′を、
X`′=S′VT R
により算出する。
Step4; 'from and the main component of the load matrix V R, the vector X of the HRV feature quantity' principal component score vector S reconfiguration feature vector X ` 'reconstituted from main component,
X'=S'V T R
Calculate by
Step5 ;Q統計量を、次元圧縮の残差を示す値、つまり上記再構成特徴量ベクトルX`′と特徴量ベクトルX′との二乗誤差の総和として、
Q=Σ(X`′−X′)2
により算出する。
Step5; The Q statistic is a value indicating the residual of dimensional compression, that is, the sum of squared errors of the reconstructed feature vector X′′ and the feature vector X′,
Q=Σ(X'-X') 2
Calculate by
認知能力変化予測装置10Aは、最後に学習データ作成部103の制御の下、ステップS25において、区間t (1<t <n)ごとに、区間t ,t-1 におけるHRV特徴量と、区間t-1 をもとに算出された区間t におけるHRV変化量(T2統計量およびQ統計量)を説明変数とし、かつ区間t+1 と区間t のATMTエラー率の増減を被説明変数とする学習データを作成し、当該学習データを学習DB40に保存する。 Finally, under the control of the learning data creation unit 103, the cognitive performance change prediction apparatus 10A, in step S25, for each section t (1<t<n), the HRV feature amount in sections t 1 and t−1 and the section t 1. The HRV variation (T 2 statistic and Q statistic) in the interval t calculated based on -1 is used as the explanatory variable, and the increase or decrease of the ATMT error rate in the interval t+1 and the interval t is used as the explained variable. Learning data is created and the learning data is stored in the learning DB 40.
なお、以上述べた学習データの作成処理は、認知能力変化の予測対象となるユーザが複数の場合にはこれらの各ユーザについてそれぞれ行われる。 Note that the above-described learning data creation process is performed for each of these users when there are a plurality of users who are targets for prediction of changes in cognitive ability.
なお、以上の説明では特定の一人の対象ユーザについて学習データを作成する場合を例にとって述べた。しかし、認知能力変化の予測対象として想定される複数のユーザについても同様に学習データを作成し、学習DB40に記憶する。 In the above description, the case where learning data is created for one specific target user has been described as an example. However, learning data is similarly created and stored in the learning DB 40 for a plurality of users who are assumed to be targets of prediction of changes in cognitive ability.
(2)ユーザの認知能力変化の予測処理
上記学習データの作成が終了すると、認知能力変化予測装置10Aは予測対象ユーザについて認知能力の変化を予測するための処理を以下のように実行する。図5はその処理手順と処理内容を示すフローチャートである。
(2) User's Cognitive Ability Change Prediction Process When the creation of the learning data is completed, the cognitive ability change prediction device 10A executes a process for predicting a change in the cognitive ability of the prediction target user as follows. FIG. 5 is a flowchart showing the processing procedure and processing contents.
対象ユーザには脈波センサ22が装着される。脈波センサ22は対象ユーザの脈波を計測し、計測された脈波の波形データをもとに心拍間隔データを算出する。そして、算出された心拍間隔データを無線インタフェースを介して認知能力変化予測装置10Aに向け送信する。 The pulse wave sensor 22 is attached to the target user. The pulse wave sensor 22 measures the pulse wave of the target user and calculates heartbeat interval data based on the waveform data of the measured pulse wave. Then, the calculated heartbeat interval data is transmitted to the cognitive performance change prediction device 10A via the wireless interface.
認知能力変化予測装置10Aは、先ず特徴量算出部104の制御の下、ステップS31により上記脈波センサ22から時系列的に送信される心拍間隔データを受信して取り込み、図示しない記憶部に一旦保存する。そして、先ず現時点から最新の10分間に上記脈波センサ22から受信される心拍間隔データについて、5分の時間窓を6秒のスライド幅でスライドさせながらそれぞれHRV特徴量を算出する。このとき算出されるHRV特徴量は、先に述べた学習データ作成時に特徴量算出部101で算出された(1) 〜(7) の値と同じである。そして、上記算出された最新の10分間のHRV特徴量を特徴量変化算出部105に渡す。以後、10分ごとに上記HRV特徴量の算出処理を繰り返す。 Under the control of the feature amount calculation unit 104, the cognitive ability change prediction device 10A first receives and takes in the heartbeat interval data transmitted in time series from the pulse wave sensor 22 in step S31, and temporarily stores the data in a storage unit (not shown). save. Then, first, for the heartbeat interval data received from the pulse wave sensor 22 in the latest 10 minutes from the present time point, the HRV characteristic amount is calculated while sliding the time window of 5 minutes with the slide width of 6 seconds. The HRV feature amount calculated at this time is the same as the values (1) to (7) calculated by the feature amount calculation unit 101 at the time of creating the learning data described above. Then, the calculated latest HRV feature amount for 10 minutes is passed to the feature amount change calculation unit 105. After that, the calculation process of the HRV feature amount is repeated every 10 minutes.
認知能力変化予測装置10Aは、次に特徴量変化算出部105の制御の下、ステップS32によりHRV変化量を以下のように算出する。すなわち、上記特徴量算出部104から渡されたHRV特徴量を順次記憶部に保存する。そして、先ず現時点から10分前までの区間におけるHRV特徴量を上記記憶部から読み込み、10分前から20分前の区間に測定・算出されたHRV特徴量をもとに、HRV変化量を算出する。具体的には、上記10分前から20分前の区間に測定・算出されたHRV特徴量に主成分分析を適用した際の負荷量行列VRおよび主成分得点ベクトルS′より、上記現時点から10分前までのHRV特徴量を用いて、T2統計量およびQ統計量を算出する。 Next, the cognitive ability change prediction device 10A calculates the HRV change amount as follows in step S32 under the control of the feature amount change calculation unit 105. That is, the HRV feature amount passed from the feature amount calculation unit 104 is sequentially stored in the storage unit. Then, first, the HRV feature amount in the section from the present time to 10 minutes ago is read from the storage unit, and the HRV change amount is calculated based on the HRV feature amount measured and calculated in the section 10 minutes to 20 minutes ago. To do. Specifically, from the above-mentioned present time, from the load matrix V R and the principal component score vector S′ when the principal component analysis is applied to the HRV feature amount measured and calculated in the section from 10 minutes before to 20 minutes before, The T 2 statistic and the Q statistic are calculated using the HRV features up to 10 minutes ago.
なお、各心拍変動の特徴量間の変化量または変化率を算出する以外に、前記各心拍変動の特徴量間の相関係数の変化量または変化率を算出するようにしてもよい。 Instead of calculating the change amount or change rate between the feature amounts of each heartbeat variability, the change amount or change rate of the correlation coefficient between the feature amounts of each heartbeat change may be calculated.
認知能力変化予測装置10Aは、最後に認知能力変化予測部106の制御の下、ステップS33において、先ず学習DB40からATMTエラー率の増減を表すデータを読み込んで、SVM(Support Vector Machine)により予測用回帰モデルの分類器を作成する。そして、上記特徴量変化算出部105により算出されたHRV特徴量と、その変化量(T2 統計量およびQ統計量)を上記分類器に入力することで、ATMTエラー率の増減、つまり「増えるか」または「減るか」を予測する。そして、当該予測結果を示す表示データを生成し、表示部30へ送信する。この結果、表示部30には、例えば「今後10分間は直前10分間に比べて認知能力が低下すると予測される」といったメッセージが表示される。 Under the control of the cognitive performance change prediction unit 106, the cognitive performance change prediction device 10A finally reads data representing increase or decrease in the ATMT error rate from the learning DB 40 in step S33, and uses the SVM (Support Vector Machine) for prediction. Create a regression model classifier. By inputting the HRV feature amount calculated by the feature amount change calculation unit 105 and the change amount (T 2 statistic amount and Q statistic amount) to the classifier, the ATMT error rate increases or decreases, that is, “increases”. "Or" or "will decrease"? Then, display data indicating the prediction result is generated and transmitted to the display unit 30. As a result, the display unit 30 displays a message, for example, "it is predicted that the cognitive ability will decrease in the next 10 minutes compared to the immediately preceding 10 minutes".
(効果)
以上詳述したように第1の実施形態では、事前に学習データ作成部103において、区間t (1<t <n)ごとに、特徴量算出部101により算出された、区間t ,t-1 におけるHRV特徴量と、特徴量変化算出部102により算出された、区間t-1 をもとに算出された区間t におけるHRV変化量(T2統計量およびQ統計量)を説明変数とし、かつATMT入力部23により入力された、区間t+1 と区間t との間のATMTエラー率の増減を表す2値データを被説明変数とする学習データを作成し、学習DB40に保存する。そして、この状態で予測対象ユーザの心拍間隔データをもとに特徴量算出部104により所定の時間分ごとにHRV特徴量を算出し、当該HRV特徴量をもとに特徴量変化算出部105によりHRV変化量を示すT2統計量およびQ統計量を算出する。そして、認知能力変化予測部106により、学習DB40からATMTエラー率の増減を表すデータを読み込んで分類器を作成し、上記算出されたHRV特徴量と、その変化量(T2 統計量およびQ統計量)を上記分類器に入力することでATMTエラー率の増減を予測して、その結果を表示部30に表示するようにしている。
(effect)
As described above in detail, in the first embodiment, in the learning data creation unit 103, the intervals t 1, t−1 calculated by the feature amount calculation unit 101 for each interval t (1<t<n) in advance. , And the HRV change amount (T 2 statistic and Q statistic) in the section t calculated based on the section t−1 calculated by the feature change calculation unit 102 as an explanatory variable, and Learning data having binary data, which is input by the ATMT input unit 23 and represents the increase or decrease of the ATMT error rate between the sections t+1 and t, as an explained variable, is created and stored in the learning DB 40. Then, in this state, the feature amount calculation unit 104 calculates the HRV feature amount for each predetermined time based on the heartbeat interval data of the prediction target user, and the feature amount change calculation unit 105 calculates the HRV feature amount based on the HRV feature amount. A T 2 statistic and a Q statistic showing the HRV change amount are calculated. Then, the cognitive ability change prediction unit 106 reads data representing the increase or decrease of the ATMT error rate from the learning DB 40 to create a classifier, and calculates the above-mentioned HRV feature amount and its change amount (T 2 statistic and Q statistic). By inputting (amount) into the classifier, increase/decrease in the ATMT error rate is predicted, and the result is displayed on the display unit 30.
従って、予測対象ユーザの認知能力の変化が予測されてその結果が表示されることで、例えば認知能力が低下する前に予測対象ユーザが自主的に休息を取ったり、管理者が予測対象ユーザに対し休息を与える等の対策を適切なタイミングで講じることが可能となる。その結果、例えばパーソナルコンピュータを使用してタスクを実行しているユーザの作業効率や品質を高く維持することが可能となる。 Therefore, by predicting the change in the cognitive ability of the prediction target user and displaying the result, for example, the prediction target user takes a voluntary rest before the cognitive ability declines, or the administrator sets the prediction target user as the prediction target user. It is possible to take measures such as giving rest at appropriate timing. As a result, it is possible to maintain high work efficiency and quality of a user who is performing a task using, for example, a personal computer.
[第2の実施形態]
この発明の第2の実施形態は、予測する認知能力変化の増減を、2値ではなく連続値である増減量として予測するようにしたものである。
[Second Embodiment]
In the second embodiment of the present invention, the predicted increase/decrease in the change in cognitive ability is predicted as an increase/decrease amount that is a continuous value rather than a binary value.
(1)学習データの作成
図6は、認知能力変化予測装置10Aによる学習データ作成処理の手順と処理内容を示すフローチャートである。なお、同図において前記図4と処理内容が同一のステップには同一符号を付して説明を行う。
(1) Creation of Learning Data FIG. 6 is a flowchart showing a procedure and processing contents of learning data creation processing by the cognitive ability change prediction device 10A. In the figure, the steps having the same processing contents as those in FIG. 4 will be described with the same reference numerals.
第2の実施形態において前記第1の実施形態と異なるところは、ATMTエラー率の変化量を被説明変数として学習DB40に保存している点である。すなわち、認知能力変化予測装置10Aの学習データ作成部103は、学習データを作成する際に、ステップS26において、区間t (1<t <n)ごとに、区間t ,t-1 におけるHRV特徴量と、区間t-1 をもとに算出された区間t におけるHRV変化量(T2統計量およびQ統計量)を説明変数とし、かつ区間t+1 と区間t のATMTエラー率の変化量、つまりATMTエラー率そのものを被説明変数とする学習データを作成する。そして、この作成した学習データを学習DB40に保存する。 The second embodiment differs from the first embodiment in that the amount of change in the ATMT error rate is stored in the learning DB 40 as the explained variable. That is, the learning data creation unit 103 of the cognitive ability change prediction device 10A creates the learning data, in step S26, for each interval t (1<t<n), the HRV feature amount in the interval t 1, t-1. And the HRV change amount (T 2 statistic and Q statistic) in the section t calculated based on the section t−1 as an explanatory variable, and the change amount of the ATMT error rate in the section t+1 and the section t, In other words, the learning data is created with the ATMT error rate itself as the explained variable. Then, the created learning data is stored in the learning DB 40.
(2)ユーザの認知能力変化の予測処理
上記学習データの作成が終了すると、認知能力変化予測装置10Aは予測対象ユーザについて認知能力の変化を予測するための処理を以下のように実行する。図7はその処理手順と処理内容を示すフローチャートである。なお、同図において前記図5と処理内容が同一のステップには同一符号を付して説明を行う。
(2) User's Cognitive Ability Change Prediction Process When the creation of the learning data is completed, the cognitive ability change prediction device 10A executes a process for predicting a change in the cognitive ability of the prediction target user as follows. FIG. 7 is a flowchart showing the processing procedure and processing contents. In the figure, steps having the same processing contents as those in FIG. 5 will be described with the same reference numerals.
認知能力変化予測装置10Aの認知能力変化予測部106は、ステップS34において、先ず学習DB40からATMTエラー率の変化量を表すデータを読み込み、SVR(Support Vector Regression)モデルによって推定用回帰モデルを作成する。そして、上記特徴量変化算出部105により算出されたHRV特徴量と、その変化量(T2 統計量およびQ統計量)を上記推定用回帰モデルに入力することでATMTエラー率の変化量を予測し、当該予測結果を示す表示データを生成して表示部30へ送信する。なお、推定に用いる回帰モデルとしては、SVRに限らず、PLS(Partial Least Square)回帰モデルなど他の重回帰モデルを用いてもよい。 In step S34, the cognitive performance change prediction unit 106 of the cognitive performance change prediction device 10A first reads data representing the amount of change in the ATMT error rate from the learning DB 40 and creates an estimation regression model using an SVR (Support Vector Regression) model. .. Then, the change amount of the ATMT error rate is predicted by inputting the HRV feature amount calculated by the feature amount change calculation unit 105 and the change amount (T 2 statistic and Q statistic) into the estimation regression model. Then, the display data indicating the prediction result is generated and transmitted to the display unit 30. The regression model used for estimation is not limited to SVR, and other multiple regression models such as PLS (Partial Least Square) regression model may be used.
以上のような処理を実行することで、ATMTエラー率の増減に加え、その増加量または減少量を予測することができる。そして、当該予測結果を示す表示データを生成し、表示部30へ送信することで、表示部30に例えば「今後10分間は直前10分間に比べて認知能力が20%低下すると予測される」といったメッセージを表示することができる。従って、予測対象ユーザに対し、休息させるか否かに止まらず、どの程度の時間休息させればよいか、さらには適当な休息の種類(仮眠、体操、軽食等)を指示することが可能となる。 By executing the processing as described above, it is possible to predict the increase or decrease of the ATMT error rate and the increase or decrease thereof. Then, by generating display data indicating the prediction result and transmitting the display data to the display unit 30, the display unit 30 displays, for example, that “the cognitive ability is predicted to decrease by 20% for the next 10 minutes compared to the immediately preceding 10 minutes”. A message can be displayed. Therefore, it is possible to instruct the prediction target user not only whether or not to rest, but also for how long the rest should be performed and the type of proper rest (nap, gymnastics, snacks, etc.). Become.
[第3の実施形態]
この発明の第3の実施形態は、対象ユーザの認知能力変化の予測処理後に、当該予測対象期間において実際に計測・算出されたHRV特徴量をもとに上記対象ユーザの認知能力変化を推定し、この認知能力変化の推定結果をもとに新たな学習データを作成して学習DBに追加するようにしたものである。
[Third Embodiment]
The third embodiment of the present invention estimates the change in cognitive ability of the target user based on the HRV feature amount actually measured and calculated in the prediction target period after the process of predicting the change in cognitive ability of the target user. , Based on the estimation result of the change in cognitive ability, new learning data is created and added to the learning DB.
図8は、この発明の第3の実施形態に係る認知能力変化予測装置の構成を示すブロック図である。なお、同図において前記図3と同一部分には同一符号を付して詳しい説明は省略する。 FIG. 8 is a block diagram showing the configuration of the cognitive ability change prediction device according to the third embodiment of the present invention. In the figure, the same parts as those in FIG. 3 are designated by the same reference numerals, and detailed description thereof will be omitted.
本実施形態の認知能力変化予測装置10Bは、認知能力変化推定部107をさらに備えている。この認知能力変化推定部107は、対象ユーザの認知能力変化の予測処理後に、当該予測対象期間において特徴量算出部104により実際に算出されたHRV特徴量を読み込む。そして、この算出されたHRV特徴量をもとに上記予測対象期間における対象ユーザの認知能力を推定し、その推定結果を学習データ作成部103に渡して新たな学習データを作成させる。 The cognitive ability change prediction device 10B of the present embodiment further includes a cognitive ability change estimation unit 107. The cognitive ability change estimation unit 107 reads the HRV feature amount actually calculated by the feature amount calculation unit 104 during the prediction target period after the target user's cognitive ability change prediction process. Then, the cognitive ability of the target user in the prediction target period is estimated based on the calculated HRV feature amount, and the estimation result is passed to the learning data creating unit 103 to create new learning data.
このような構成であるから、認知能力変化予測装置10Bでは、認知能力変化予測部106により例えば今後10分間における対象ユーザの認知能力変化の予測処理が終了すると、続いて認知能力変化推定部107により上記今後10分間における実際の認知能力の推定が行われる。 With such a configuration, in the cognitive performance change prediction device 10B, when the cognitive performance change prediction unit 106 finishes the prediction process of the cognitive performance change of the target user for the next 10 minutes, for example, the cognitive performance change estimation unit 107 then The actual cognitive ability is estimated in the next 10 minutes.
この認知能力の推定処理は以下のように行われる。すなわち、先ず特徴量算出部により、脈波センサ22から受信した最新の10分間の心拍間隔データについて、5分の時間窓を6秒のスライド幅でスライドさせながらそれぞれHRV特徴量を算出する。次に認知能力変化推定部107により、上記算出されたHRV特徴量をもとに、上記10分間における対象ユーザの実際の認知能力変化を推定する。具体的には、認知能力変化の推定結果に基づくATMTエラー率の変化量を推定する。なお、上記HRV特徴量を用いた認知能力変化の推定処理法方については、非特許文献2に詳しく記載されている。 The process of estimating the cognitive ability is performed as follows. That is, first, the feature amount calculation unit calculates the HRV feature amount for each of the latest 10-minute heartbeat interval data received from the pulse wave sensor 22 while sliding a 5-minute time window with a slide width of 6 seconds. Next, the cognitive ability change estimation unit 107 estimates the actual cognitive ability change of the target user in the 10 minutes based on the calculated HRV feature amount. Specifically, the change amount of the ATMT error rate is estimated based on the estimation result of the change in cognitive ability. Note that Non-Patent Document 2 describes in detail how to estimate the change in cognitive ability using the HRV feature amount.
認知能力変化推定部107は、上記算出されたATMTエラー率の変化量の推定結果を、認知能力変化の予測処理過程において特徴量算出部104および特徴量変化算出部105によりそれぞれ算出されたHRV特徴量およびHRV変化量(T2 統計量とQ統計量)と共に、学習データ作成部103に渡す。学習データ作成部103は、上記認知能力変化推定部107から渡された、HRV特徴量およびHRV変化量(T2 統計量とQ統計量)を説明変数とし、かつATMTエラー率の変化量の推定結果を被説明変数とする新たな学習データを作成し、この作成した学習データを学習DB40に書き込む。 The cognitive ability change estimation unit 107 calculates the estimation result of the change amount of the ATMT error rate, which is calculated by the feature amount calculation unit 104 and the feature amount change calculation unit 105 in the process of predicting the change in cognitive ability. The amount and the HRV change amount (T 2 statistic and Q statistic) are passed to the learning data creation unit 103. The learning data creation unit 103 uses the HRV feature amount and the HRV change amount (T 2 statistic amount and Q statistic amount) passed from the cognitive ability change estimation unit 107 as explanatory variables, and estimates the change amount of the ATMT error rate. New learning data with the result as the explained variable is created, and the created learning data is written in the learning DB 40.
以上述べたように第3の実施形態では、対象ユーザの認知能力変化の予測処理後に、当該予測対象期間において実際に計測・算出されたHRV特徴量をもとに上記対象ユーザの認知能力変化を推定し、この認知能力変化の推定結果をもとに新たな学習データを作成して学習DBに追加するようにしている。従って、対象ユーザの最新の期間における認知能力変化の推定結果に基づく学習データを記憶することができ、これによりそれ以降の認知能力変化の予測をより一層精度良く行うことが可能となる。 As described above, in the third embodiment, after the target user's cognitive ability change prediction processing, the target user's cognitive ability change is calculated based on the HRV feature amount actually measured and calculated in the prediction target period. The learning data is estimated and new learning data is created based on the estimation result of the change in cognitive ability and is added to the learning DB. Therefore, the learning data based on the estimation result of the change in cognitive ability of the target user in the latest period can be stored, and thereby, the change in cognitive ability thereafter can be predicted more accurately.
[その他の実施形態]
認知能力の指標としては、ATMTのエラー率以外に、ATMTの正解数やエラー数を用いてもよく、さらにはタスク中のキーボードの打鍵数やアクティブウィンドウの切換数等のパーソナルコンピュータの操作ログを用いてもよい。
[Other Embodiments]
In addition to the error rate of ATMT, the number of correct answers and the number of errors of ATMT may be used as an index of cognitive ability. You may use.
また、前記各実施形態では認知能力変化予測装置をサーバコンピュータ等に設けた場合を例にとって説明したが、対象ユーザが所持するスマートフォンやタブレット型端末等の携帯端末に設けてもよい。
その他、認知能力変化予測装置の構成や予測処理の手順、処理内容等については、この発明の要旨を逸脱しない範囲で種々変形して実施可能である。
Further, in each of the above-described embodiments, the case in which the cognitive ability change prediction device is provided in the server computer or the like has been described as an example, but the device may be provided in a mobile terminal such as a smartphone or a tablet type terminal owned by the target user.
In addition, the configuration of the cognitive ability change prediction device, the procedure of the prediction processing, the processing content, and the like can be variously modified and implemented without departing from the scope of the present invention.
要するにこの発明は、上記各実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記各実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、各実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態に亘る構成要素を適宜組み合せてもよい。 In short, the present invention is not limited to the above-mentioned embodiments as they are, and can be embodied by modifying the constituent elements within a range not departing from the gist of the invention in an implementation stage. Further, various inventions can be formed by appropriately combining a plurality of constituent elements disclosed in each of the above embodiments. For example, some constituent elements may be deleted from all the constituent elements shown in each embodiment. Furthermore, the constituent elements of different embodiments may be combined appropriately.
1…処理装置、2…データ入力手段、3…結果出力手段と、4…学習データ蓄積手段、11…特徴量算出手段と、12…特徴量変化算出手段、13…認知能力変化予測手段と、21,22…脈波センサ、23…ATMT入力部、30…表示部、40…学習データベース、101,104…特徴量算出部、102,105…特徴量変化算出部、103…学習データ作成部、106…認知能力変化予測部、107…認知能力変化推定部。 DESCRIPTION OF SYMBOLS 1... Processing device, 2... Data input means, 3... Result output means, 4... Learning data storage means, 11... Feature amount calculation means, 12... Feature amount change calculation means, 13... Cognitive ability change prediction means, 21, 22... Pulse wave sensor, 23... ATMT input unit, 30... Display unit, 40... Learning database, 101, 104... Feature amount calculating unit, 102, 105... Feature amount change calculating unit, 103... Learning data creating unit, 106... Cognitive ability change prediction unit, 107... Cognitive ability change estimation unit.
Claims (6)
前記取得された心拍間隔を表すデータをもとに、時間区間tにおける心拍変動の特徴量と、1つ前の時間区間t−1における心拍変動の特徴量をそれぞれ算出する手段と、
前記算出された各心拍変動の特徴量をもとに、当該各心拍変動の特徴量の変化量を算出する手段と、
前記算出された各心拍変動の特徴量およびその変化量と、予め記憶されている学習データとをもとに、前記予測対象ユーザの1つ後の時間区間t+1における認知能力の変化を、前記時間区間tにおける認知能力からの増減として予測する手段と、
前記予測対象ユーザの認知能力の変化の予測結果を表す情報を出力する手段と
を具備することを特徴とする認知能力変化予測装置。 Means for acquiring data representing the heartbeat interval of the prediction target user,
Means for calculating the characteristic amount of the heartbeat fluctuation in the time interval t and the characteristic amount of the heartbeat fluctuation in the immediately preceding time interval t-1, based on the acquired data representing the heartbeat interval,
Based on the calculated feature amount of each heartbeat variability, means for calculating the amount of change in the feature amount of each heartbeat variability,
Based on the calculated feature amount of each heartbeat variability and its change amount, and the learning data stored in advance, the change in the cognitive ability of the prediction target user in the subsequent time interval t+1 is calculated as the time. Means for predicting increase/decrease from cognitive ability in section t,
A unit for outputting information representing a prediction result of a change in cognitive ability of the prediction target user.
前記取得された心拍間隔を表すデータをもとに時間領域の特徴量を算出する手段と、
前記取得された心拍間隔を表すデータをもとに周波数領域の特徴量を算出する手段と
を備え、
前記時間領域の特徴量を算出する手段は、心拍間隔の平均値と、心拍間隔の標準偏差(SDNN)と、隣接した心拍間隔の差が50msecを超えた比率(pNN50)と、隣接した心拍間隔の差の二乗平均平方和(RMSSD)とのうちの少なくとも1つを算出し、
前記周波数領域の特徴量を算出する手段は、心拍間隔をスペクトル解析した際の低周波成分(0.04~0.15Hz)の総パワー(LF)と、心拍間隔をスペクトル解析した際の高周波成分(0.15~0.40Hz)の総パワー(HF)と、前記低周波成分の総パワーLFと前記高周波成分の総パワーHFとの比率(LF/HF)とのうちの少なくとも1つを算出する
ことを特徴とする請求項1に記載の認知能力変化予測装置。 The means for calculating the feature amount of the heartbeat variability,
A means for calculating a characteristic amount in the time domain based on the data representing the obtained heartbeat interval;
A means for calculating a feature amount in the frequency domain based on the acquired data representing the heartbeat interval,
The means for calculating the feature amount in the time domain includes an average value of heartbeat intervals, a standard deviation of heartbeat intervals (SDNN), a ratio (pNN50) in which the difference between adjacent heartbeat intervals exceeds 50 msec, and an adjacent heartbeat interval. And calculating at least one of the sum of root mean square (RMSSD) of the differences,
The means for calculating the feature amount in the frequency domain is the total power (LF) of the low frequency component (0.04 to 0.15 Hz) when the heartbeat interval is spectrally analyzed and the high frequency component (0.15 to 0.40 Hz) total power (HF) and at least one of a ratio (LF/HF) of the total power LF of the low-frequency component and the total power HF of the high-frequency component. The apparatus for predicting change in cognitive ability according to claim 1.
ことを特徴とする請求項1に記載の認知能力変化予測装置。 The means for calculating the change amount of the feature amount of the heartbeat variability is the change amount or the change between the feature amount of the heartbeat change in the time period t and the feature amount of the heartbeat change in the immediately preceding time period t-1. Rate, the change amount or change rate of the correlation coefficient between the feature amount of the heartbeat variability in the time section t and the feature amount of the heartbeat change in the immediately preceding time section t-1. The cognitive ability change prediction device according to claim 1, wherein
前記学習データをもとに認知能力の変化を表す分類器を作成し、前記算出された各心拍変動の特徴量とその変化量を前記分類器に入力することで、前記認知能力の変化を予測する手段と、
前記学習データをもとに回帰モデルを作成し、前記算出された各心拍変動の特徴量とその変化量を前記回帰モデルに入力することで、前記認知能力の変化量を予測する手段と
のうちの何れか一方を備えることを特徴とする請求項1に記載の認知能力変化予測装置。 The means for predicting a change in the cognitive ability,
A classifier that represents a change in cognitive ability is created based on the learning data, and the calculated feature amount of each heartbeat variability and its change amount are input to the classifier to predict the change in the cognitive ability. Means to do
A means for predicting a change amount of the cognitive ability by creating a regression model based on the learning data and inputting the calculated feature amount of each heartbeat variability and its change amount into the regression model. The cognitive ability change prediction apparatus according to claim 1, further comprising:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016118918A JP6712909B2 (en) | 2016-06-15 | 2016-06-15 | Cognitive change prediction device and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016118918A JP6712909B2 (en) | 2016-06-15 | 2016-06-15 | Cognitive change prediction device and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2017221416A JP2017221416A (en) | 2017-12-21 |
| JP6712909B2 true JP6712909B2 (en) | 2020-06-24 |
Family
ID=60687428
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016118918A Active JP6712909B2 (en) | 2016-06-15 | 2016-06-15 | Cognitive change prediction device and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6712909B2 (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2021248744A1 (en) * | 2020-02-07 | 2022-08-25 | Eyetracking Llc | Evaluation of a person or system through measurement of physiological data |
| JP7561365B2 (en) * | 2020-06-18 | 2024-10-04 | パナソニックIpマネジメント株式会社 | Information processing method, information processing system, and program |
| CN113468077B (en) * | 2021-09-06 | 2021-12-10 | 北京无疆脑智科技有限公司 | Cognitive ability testing method, device, electronic device and storage medium |
| EP4541270A4 (en) * | 2022-07-14 | 2025-07-02 | Medical Ai Co Ltd | Method, program and device for generating long-term heart rate variability based on short-term measurement signals |
| CN121100277A (en) * | 2023-05-15 | 2025-12-09 | 株式会社岛津制作所 | Waveform analysis method, waveform analysis device, and analysis device |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110230732A1 (en) * | 2009-09-14 | 2011-09-22 | Philometron, Inc. | System utilizing physiological monitoring and electronic media for health improvement |
| EP2477538A4 (en) * | 2009-09-14 | 2013-01-09 | Carl Frederick Edman | System utilizing physiological monitoring and electronic media for health improvement |
| JP6344912B2 (en) * | 2013-12-13 | 2018-06-20 | 国立大学法人京都大学 | Epileptic seizure sign detection apparatus, epileptic seizure sign detection model generation apparatus, epileptic seizure sign detection method, epileptic seizure sign detection model generation method, epileptic seizure sign detection program, and epileptic seizure sign detection model generation program |
| JP6216294B2 (en) * | 2014-07-30 | 2017-10-18 | 日本電信電話株式会社 | Multiple regression analysis apparatus and multiple regression analysis method |
-
2016
- 2016-06-15 JP JP2016118918A patent/JP6712909B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2017221416A (en) | 2017-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6712909B2 (en) | Cognitive change prediction device and program | |
| CA3070980C (en) | Methods and systems for forecasting seizures | |
| CN106202850B (en) | Exercise management methods and equipment | |
| US10360806B2 (en) | Method and apparatus for providing exercise program based on feedback | |
| US20130260357A1 (en) | Skill Screening | |
| CN111629667A (en) | Exercise load estimation method, exercise load estimation device, and recording medium | |
| JP2019534088A (en) | System, method, and computer program for quantifying physical fatigue of a subject | |
| CN105144174A (en) | Organism age prediction device and method | |
| Ali-Eldin et al. | Measuring cloud workload burstiness | |
| JP2018011721A (en) | Stress estimation device, stress estimation method, and stress estimation program | |
| JP2016147010A (en) | Work alertness estimation device, method and program | |
| Sonne et al. | A modified version of the three-compartment model to predict fatigue during submaximal tasks with complex force-time histories | |
| JP2022098608A (en) | Heartbeat fluctuation index estimation device, heartbeat fluctuation index estimation method, and heartbeat fluctuation index estimation program | |
| US20220406439A1 (en) | System and method for adaptive configuration of computerized cognitive training programs | |
| US8694635B1 (en) | Time series technique for analyzing performance in an online professional network | |
| US20190114934A1 (en) | System and method for evaluating actions performed to achieve communications | |
| JP5913493B2 (en) | Deficient biological signal estimation method | |
| JP5993968B2 (en) | Sleep stage estimation apparatus, method and program | |
| JP2021101846A (en) | Information processing device and program | |
| JP7347822B2 (en) | Mental and physical state estimation device, mental and physical state estimation method, and program | |
| US20240206825A1 (en) | Estimation device, estimation system, estimation method, and storage medium | |
| JP2022049584A (en) | Information processing device and program | |
| Fakhfakh et al. | QoS aggregation for service orchestrations based on workflow pattern rules and MCDM method: evaluation at design time and runtime | |
| US20230317235A1 (en) | Effective mobility | |
| US20230010067A1 (en) | Systems and methods for measuring performance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180620 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20190328 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190409 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190605 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20191126 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200110 |
|
| 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: 20200526 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200602 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6712909 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |