JP6802467B2 - Processing device, judgment device and appropriate range identification method - Google Patents

Processing device, judgment device and appropriate range identification method Download PDF

Info

Publication number
JP6802467B2
JP6802467B2 JP2016161684A JP2016161684A JP6802467B2 JP 6802467 B2 JP6802467 B2 JP 6802467B2 JP 2016161684 A JP2016161684 A JP 2016161684A JP 2016161684 A JP2016161684 A JP 2016161684A JP 6802467 B2 JP6802467 B2 JP 6802467B2
Authority
JP
Japan
Prior art keywords
distribution
value
appropriate range
determination
processing
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
Application number
JP2016161684A
Other languages
Japanese (ja)
Other versions
JP2018029634A (en
Inventor
大森 信行
信行 大森
智啓 村澤
智啓 村澤
淳平 相澤
淳平 相澤
吉人 小山
吉人 小山
浩 栗田
浩 栗田
英哉 百瀬
英哉 百瀬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NAGANO PREFECTURAL GOVERNMENT
Shinshu University NUC
Skinos Co Ltd
Original Assignee
NAGANO PREFECTURAL GOVERNMENT
Shinshu University NUC
Skinos Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NAGANO PREFECTURAL GOVERNMENT, Shinshu University NUC, Skinos Co Ltd filed Critical NAGANO PREFECTURAL GOVERNMENT
Priority to JP2016161684A priority Critical patent/JP6802467B2/en
Publication of JP2018029634A publication Critical patent/JP2018029634A/en
Application granted granted Critical
Publication of JP6802467B2 publication Critical patent/JP6802467B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Description

本発明は、適正範囲内の判定用パラメータを用いて嚥下に関する事項を判定する際の適正範囲を特定する処理装置、およびその処理装置を備えて嚥下に関する事項を判定する判定装置、および適正範囲内の判定用パラメータを用いて嚥下に関する事項を判定する際の適正範囲を特定する適正範囲特定方法に関するものである。 The present invention is a processing device that specifies an appropriate range when determining a matter related to swallowing using a determination parameter within an appropriate range, a determination device that is provided with the processing device and determines a matter related to swallowing, and within the appropriate range. It relates to an appropriate range specifying method for specifying an appropriate range when determining a matter related to swallowing using the determination parameters of.

飲食物を嚥下するときの筋活動に関する被測定量を測定する装置として、下記特許文献1において出願人が開示した嚥下活動測定装置が知られている。この嚥下活動測定装置は、センサシート、レコーダ部、測定データ分析部、筋活動分析部および記録・表示部を備えて構成されている。この嚥下活動測定装置では、被験者の前頚部に貼りつけたセンサシートに配設されている複数の筋電図電極および複数の振動ピックアップから出力される電気信号が記憶装置に保存され、その電気信号に対して測定データ分析部がフィルタ処理を行う。次いで、筋活動分析部が、測定データ分析部によってフィルタ処理された電気信号に基づき、嚥下に関連する筋群の活動継続時間や、筋群の活動開始および活動終了の順序(嚥下に関する事項の判定に用いるパラメータ)を計算する。 As a device for measuring a measured amount related to muscle activity when swallowing food and drink, a swallowing activity measuring device disclosed by the applicant in Patent Document 1 below is known. This swallowing activity measuring device is configured to include a sensor sheet, a recorder unit, a measurement data analysis unit, a muscle activity analysis unit, and a recording / display unit. In this swallowing activity measuring device, electric signals output from a plurality of electromyographic electrodes and a plurality of vibration pickups arranged on a sensor sheet attached to the anterior neck of the subject are stored in a storage device, and the electric signals are stored in the storage device. The measurement data analysis unit performs the filtering process. Next, the muscle activity analysis unit determines the activity duration of the muscle group related to swallowing and the order of the start and end of the activity of the muscle group (determination of matters related to swallowing) based on the electrical signal filtered by the measurement data analysis unit. (Parameter used for) is calculated.

特開2014−168490号公報(第7−8頁、第1図)Japanese Unexamined Patent Publication No. 2014-168490 (Pages 7-8, Fig. 1)

ところが、上記した従来の嚥下活動測定装置には、改善すべき以下の課題がある。すなわち、この嚥下活動測定装置では、測定データ分析部が筋電図電極および振動ピックアップから出力される電気信号をフィルタ処理し、フィルタ処理した後の電気信号に基づいて活動継続時間等を計算している。この場合、フィルタ処理した後の電気信号にもノイズが含まれている。また、咀嚼等の嚥下とは異なる動作に伴う信号成分が含まれていることもある。このため、これらの影響により、嚥下に関する事項の判定に用いる活動継続時間等の判定用パラメータの値が判定用としては不適正な値となることがある。具体的には、ノイズの影響によって実際の活動継続時間よりも短い時間が活動継続時間として計算されることがある。しかしながら、従来の嚥下活動測定装置には、判定用パラメータの値の適正範囲を把握する手段が備えられていないため、不適正な値の判定用パラメータを用いて嚥下に関する事項の判定を行うことによって判定結果が正確に行うことが困難なことがあり、この点の改善が望まれている。 However, the above-mentioned conventional swallowing activity measuring device has the following problems to be improved. That is, in this swallowing activity measuring device, the measurement data analysis unit filters the electric signals output from the EMG electrodes and the vibration pickup, and calculates the activity duration and the like based on the electric signals after the filtering. There is. In this case, the electrical signal after filtering also contains noise. In addition, it may contain signal components associated with movements different from swallowing such as mastication. Therefore, due to these effects, the value of the determination parameter such as the activity duration used for determining the swallowing-related matter may be an inappropriate value for the determination. Specifically, due to the influence of noise, a time shorter than the actual activity duration may be calculated as the activity duration. However, since the conventional swallowing activity measuring device does not have a means for grasping the appropriate range of the value of the determination parameter, the swallowing-related matter is determined by using the determination parameter of the inappropriate value. It may be difficult to make an accurate judgment result, and improvement in this respect is desired.

本発明は、かかる課題に鑑みてなされたものであり、嚥下に関する事項の判定に用いる判定用パラメータの値の適正範囲を特定し得る処理装置、判定装置および適正範囲特定方法を提供することを主目的とする。 The present invention has been made in view of such a problem, and mainly provides a processing device, a determination device, and an appropriate range specifying method capable of specifying an appropriate range of values of a determination parameter used for determining a matter related to swallowing. The purpose.

上記目的を達成すべく請求項1記載の処理装置は、人の頸部にセンサを装着して測定した測定パラメータに基づいて算出した判定用パラメータの値が適正範囲内のときの当該判定用パラメータを用いて嚥下に関する事項を判定する際の当該適正範囲を特定する処理部を備え、前記処理部は、複数回測定した前記測定パラメータに基づいてそれぞれ算出した複数の前記判定用パラメータを統計処理して得た処理結果に基づいて前記適正範囲の最小値とする第1閾値を特定し、前記判定用パラメータの値が前記第1閾値以上の範囲を前記適正範囲として特定する処理装置であって、前記処理部は、前記各判定用パラメータの混合ガウス分布を求める処理を前記統計処理として実行し、前記処理結果としての前記混合ガウス分布から得た結果に基づいて前記第1閾値を特定する。 The processing device according to claim 1 in order to achieve the above object is the determination parameter when the value of the determination parameter calculated based on the measurement parameter measured by attaching a sensor to the neck of a person is within an appropriate range. The processing unit is provided with a processing unit that specifies the appropriate range when determining a matter related to swallowing, and the processing unit statistically processes a plurality of the determination parameters calculated based on the measurement parameters measured a plurality of times. A processing device that specifies a first threshold value as the minimum value of the appropriate range based on the processing result obtained, and specifies a range in which the value of the determination parameter is equal to or greater than the first threshold value as the appropriate range . The processing unit executes a process for obtaining a mixed Gaussian distribution of each determination parameter as the statistical process, and specifies the first threshold value based on the result obtained from the mixed Gaussian distribution as the processing result.

また、請求項記載の処理装置は、請求項記載の処理装置において、前記処理部は、前記混合ガウス分布における混合比が最大の第1分布の確率密度分布と、前記第1分布に隣接する分布であって平均値が第1分布の平均値よりも小さい第2分布の確率密度分布とが等しい前記判定用パラメータの値のうちの当該第1分布の平均値と当該第2分布の平均値との間に位置する当該判定用パラメータの値を前記第1閾値として特定する。 Further, the processing apparatus according to claim 2 is the processing apparatus according to claim 1 , wherein the processing unit has a probability density distribution of the first distribution having the maximum mixing ratio in the mixed Gaussian distribution and an adjacent first distribution. The average value of the first distribution and the average of the second distribution among the values of the determination parameters that are equal to the probability density distribution of the second distribution whose average value is smaller than the average value of the first distribution. The value of the determination parameter located between the value and the value is specified as the first threshold value.

また、請求項記載の処理装置は、請求項1または2記載の処理装置において、前記処理部は、前記処理結果に基づいて前記適正範囲の最大値とする第2閾値を特定し、前記判定用パラメータの値が前記第1閾値以上で前記第2閾値以下の範囲を前記適正範囲として特定する。 Further, the processing apparatus according to claim 3, in the processing apparatus according to claim 1 or 2, wherein the processing unit identifies a second threshold value to the maximum value of the proper range based on the processing result, the determination A range in which the value of the parameter is equal to or higher than the first threshold value and equal to or lower than the second threshold value is specified as the appropriate range.

また、請求項記載の処理装置は、請求項記載の処理装置において、前記処理部は、前記各判定用パラメータの混合ガウス分布を求める処理を前記統計処理として実行し、前記処理結果としての前記混合ガウス分布から得た結果に基づいて前記第2閾値を特定する。 The processing apparatus according to claim 4, wherein, in the processing apparatus according to claim 3, wherein the processing unit executes a process for obtaining a mixed Gaussian distribution of the respective determination parameter as said statistical process, as the process result The second threshold is specified based on the results obtained from the mixed Gaussian distribution.

また、請求項記載の処理装置は、請求項記載の処理装置において、前記処理部は、前記混合ガウス分布における混合比が最大の第1分布の確率密度分布の値と、前記第1分布に隣接する分布であって平均値が第1分布の平均値よりも大きい第3分布の確率密度分布の値とが等しい前記判定用パラメータの値のうちの当該第1分布の平均値と当該第3分布の平均値との間に位置する当該判定用パラメータの値を前記第2閾値として特定する。 Further, the processing apparatus according to claim 5 is the processing apparatus according to claim 4 , wherein the processing unit has a value of a probability density distribution of a first distribution having a maximum mixing ratio in the mixed Gaussian distribution and the first distribution. The average value of the first distribution and the mean value of the first distribution among the values of the determination parameters that are adjacent to and have the same mean value as the value of the probability density distribution of the third distribution that is larger than the average value of the first distribution. The value of the determination parameter located between the average value of the three distributions is specified as the second threshold value.

また、請求項記載の判定装置は、請求項1からのいずれかに記載の処理装置と、当該処理装置によって特定された前記適正範囲内の前記判定用パラメータを用いて前記嚥下に関する事項を判定する判定部とを備えている。 Further, the determination device according to claim 6 uses the processing device according to any one of claims 1 to 5 and the determination parameter within the appropriate range specified by the processing device to perform matters related to swallowing. It is provided with a determination unit for determination.

また、請求項7記載の適正範囲特定方法は、人の頸部にセンサを装着して測定した測定パラメータに基づいて算出した判定用パラメータの値が適正範囲内のときの当該判定用パラメータを用いて嚥下に関する事項を判定する際の当該適正範囲を請求項1から5のいずれかに記載の処理装置における前記処理部に特定させる適正範囲特定方法であって、前記処理部に対して、複数回測定した前記測定パラメータに基づいてそれぞれ算出した複数の前記判定用パラメータを統計処理して得た処理結果に基づいて前記適正範囲の最小値とする第1閾値を特定させると共に前記判定用パラメータの値が前記第1閾値以上の範囲を前記適正範囲として特定させ、かつ前記各判定用パラメータの混合ガウス分布を求める処理を前記統計処理として実行させると共に前記処理結果としての前記混合ガウス分布から得た結果に基づいて前記第1閾値を特定させるFurther, the method for specifying the appropriate range according to claim 7 uses the determination parameter when the value of the determination parameter calculated based on the measurement parameter measured by attaching the sensor to the neck of a person is within the appropriate range. a appropriate range specification method Ru is specific to the processing unit in the processing apparatus according to the appropriate range at the time of determining the matters relating swallowing to any one of claims 1 to 5 Te, to the processing unit, a plurality Based on the processing results obtained by statistically processing a plurality of the determination parameters calculated based on the measurement parameters measured once, the first threshold value to be the minimum value of the appropriate range is specified, and the determination parameters A range having a value equal to or higher than the first threshold value is specified as the appropriate range , and a process for obtaining a mixed Gaussian distribution of each determination parameter is executed as the statistical process, and is obtained from the mixed Gaussian distribution as a result of the process. The first threshold is specified based on the result .

請求項1記載の処理装置、請求項記載の判定装置、および請求項記載の適正範囲特定方法によれば、複数回測定した測定パラメータに基づいてそれぞれ算出した複数の判定用パラメータを統計処理して得た処理結果に基づいて適正範囲の最小値とする第1閾値を特定することにより、嚥下に関する事項の判定に用いる判定用パラメータの適正範囲を確実に特定することができる。このため、この処理装置、判定装置および適正範囲特定方法によれば、算出した判定用パラメータを嚥下に関する事項の判定に用いることが適正であるか否かを判別することが困難な状況、(例えば、適正であるか否かを判別可能な熟練者が存在しない状況)においても、適正範囲と判定用パラメータとを比較することで判定用パラメータを嚥下に関する事項の判定に用いることが適正であるか否かを確実かつ容易に判別することができる結果、嚥下に関する事項の判定を正確に行うことができる。 According to the processing device according to claim 1, the determination device according to claim 6 , and the appropriate range specifying method according to claim 7, a plurality of determination parameters calculated based on measurement parameters measured a plurality of times are statistically processed. By specifying the first threshold value which is the minimum value of the appropriate range based on the processing result obtained in the above procedure, the appropriate range of the determination parameter used for determining the matter related to swallowing can be surely specified. Therefore, according to this processing device, the determination device, and the appropriate range identification method, it is difficult to determine whether or not it is appropriate to use the calculated determination parameter for determining a matter related to swallowing (for example,). Even in a situation where there is no expert who can determine whether or not it is appropriate), is it appropriate to use the determination parameter for determining matters related to swallowing by comparing the appropriate range with the determination parameter? As a result of being able to reliably and easily determine whether or not it is possible, it is possible to accurately determine matters related to swallowing.

また、の処理装置、この判定装置およびこの適正範囲特定方法によれば、各判定用パラメータの混合ガウス分布を求める処理を統計処理として実行し、処理結果としての混合ガウス分布から得た結果に基づいて第1閾値を特定することにより、適正範囲の最小値として信頼性の高い第1閾値を特定することができる。 Further, this processing apparatus, a result according to the determination device and the appropriate range specification method this, which executes the processing for obtaining a mixed Gaussian distribution of each determination parameter as statistical processing, were obtained from a mixed Gaussian distribution as the processing result By specifying the first threshold value based on the above, a highly reliable first threshold value can be specified as the minimum value in the appropriate range.

また、請求項記載の処理装置、および請求項記載の判定装置によれば、混合ガウス分布における混合比が最大の第1分布の確率密度分布と、第1分布に隣接する分布であって平均値が第1分布の平均値よりも小さい第2分布の確率密度分布とが等しい判定用パラメータの値のうちの第1分布の平均値と第2分布の平均値との間に位置する判定用パラメータの値を第1閾値として特定することにより、適正範囲の最小値としてより信頼性の高い第1閾値を特定することができる。 Further, according to the processing apparatus according to claim 2 and the determination apparatus according to claim 6 , the probability density distribution of the first distribution having the maximum mixing ratio in the mixed Gaussian distribution and the distribution adjacent to the first distribution. Judgment that the mean value is smaller than the mean value of the first distribution and the probability density distribution of the second distribution is equal. Judgment located between the mean value of the first distribution and the mean value of the second distribution among the values of the judgment parameters. By specifying the value of the parameter as the first threshold value, a more reliable first threshold value can be specified as the minimum value in the appropriate range.

また、請求項記載の処理装置、および請求項記載の判定装置によれば、複数の判定用パラメータを統計処理して得た処理結果に基づいて適正範囲の最大値とする第2閾値を特定することにより、極端に短い判定用パラメータおよび極端に長い判定用パラメータの双方を、判定に用いる判定用パラメータから除外することが可能な適正範囲を特定することができるため、判定用パラメータを嚥下に関する事項の判定に用いることが適正であるか否かをより確実に判別することができる。 The processing apparatus according to claim 3, and according to the determination apparatus of claim 6, the second threshold value to the maximum value of the proper range based on the processing result obtained by statistically processing a plurality of determination parameter By specifying, it is possible to specify an appropriate range in which both extremely short judgment parameters and extremely long judgment parameters can be excluded from the judgment parameters used for judgment, so that the judgment parameters are swallowed. It is possible to more reliably determine whether or not it is appropriate to use it for determining matters related to.

また、請求項記載の処理装置、および請求項記載の判定装置によれば、各判定用パラメータの混合ガウス分布を求める処理を統計処理として実行し、処理結果としての混合ガウス分布から得た結果に基づいて第2閾値を特定することにより、適正範囲の最大値として信頼性の高い第2閾値を特定することができる。 Further, according to the processing apparatus according to claim 4 and the determination apparatus according to claim 6 , a process for obtaining a mixed Gaussian distribution of each determination parameter is executed as statistical processing, and the mixed Gaussian distribution as a processing result is obtained. By specifying the second threshold value based on the result, it is possible to specify a highly reliable second threshold value as the maximum value in the appropriate range.

また、請求項記載の処理装置、および請求項記載の判定装置によれば、混合ガウス分布における混合比が最大の第1分布の確率密度分布の値と、第1分布に隣接する分布であって平均値が第1分布の平均値よりも大きい第3分布の確率密度分布の値とが等しい判定用パラメータの値のうちの第1分布の平均値と第3分布の平均値との間に位置する判定用パラメータの値を第2閾値として特定することにより、適正範囲の最大値としてより信頼性の高い第2閾値を特定することができる。 Further, according to the processing apparatus according to claim 5 and the determination apparatus according to claim 6 , the value of the probability density distribution of the first distribution having the maximum mixing ratio in the mixed Gaussian distribution and the distribution adjacent to the first distribution. The mean value is larger than the mean value of the first distribution. The value of the probability density distribution of the third distribution is equal. Among the values of the judgment parameters, between the mean value of the first distribution and the mean value of the third distribution. By specifying the value of the determination parameter located in, as the second threshold value, a more reliable second threshold value can be specified as the maximum value in the appropriate range.

判定装置1の構成を示す構成図である。It is a block diagram which shows the structure of the determination device 1. センサシート11の構成を示す平面図である。It is a top view which shows the structure of the sensor sheet 11. センサシート11を装着した判定対象者100(被験者300)の装着状態を示す説明図である。It is explanatory drawing which shows the wearing state of the determination target person 100 (subject 300) wearing a sensor sheet 11. 活動時間Taの算出方法を説明する第1の説明図である。It is 1st explanatory drawing explaining the calculation method of activity time Ta. 活動時間Taの算出方法を説明する第2の説明図である。It is a 2nd explanatory diagram explaining the calculation method of the activity time Ta. 適正範囲Tfの特定方法を説明する説明図である。It is explanatory drawing explaining the method of specifying an appropriate range Tf.

以下、処理装置、判定装置および適正範囲特定方法の形態について、添付図面を参照して説明する。 Hereinafter, the forms of the processing device, the determination device, and the appropriate range specifying method will be described with reference to the attached drawings.

最初に、判定装置の一例としての図1に示す判定装置1の構成について説明する。判定装置1は、嚥下に関する事項を判定可能に構成されている。具体的には、判定装置1では、嚥下に関する事項として、判定対象の人(以下「判定対象者100」ともいう)が嚥下機能低下者(嚥下の機能が低下している人)であるか嚥下機能正常者(嚥下の機能が低下していない人)であるかを嚥下に関する事項として判定(以下「嚥下機能判定」ともいう)することが可能となっている。 First, the configuration of the determination device 1 shown in FIG. 1 as an example of the determination device will be described. The determination device 1 is configured to be able to determine matters related to swallowing. Specifically, in the determination device 1, as a matter related to swallowing, whether the person to be determined (hereinafter, also referred to as “determination target person 100”) is a person with reduced swallowing function (a person with reduced swallowing function) or swallowing. It is possible to determine whether a person has normal function (a person whose swallowing function has not deteriorated) as a matter related to swallowing (hereinafter, also referred to as "swallowing function determination").

また、判定装置1は、図1に示すように、センサシート11、マイクロフォン12、操作部13、記憶部14、表示部15および処理部16を備えて構成されている。この場合、記憶部14および処理部16によって処理装置が構成される。 Further, as shown in FIG. 1, the determination device 1 includes a sensor sheet 11, a microphone 12, an operation unit 13, a storage unit 14, a display unit 15, and a processing unit 16. In this case, the processing device is composed of the storage unit 14 and the processing unit 16.

センサシート11は、図2に示すように、粘着シート21、複数(一例として8つ)の筋電図電極22a〜22h(センサの一例であって、以下、区別しないときには「筋電図電極22」ともいう)、および取り付け指標23を備えて構成されている。粘着シート21は、一面が粘着面となっている。なお、同図では、粘着シート21の粘着面とは逆側の表面から見た状態でセンサシート11を図示している。 As shown in FIG. 2, the sensor sheet 11 includes an adhesive sheet 21, and a plurality of (eight as an example) electromyographic electrodes 22a to 22h (an example of a sensor, hereinafter, when not distinguished, “electromyographic electrodes 22”. (Also referred to as), and a mounting index 23. One side of the adhesive sheet 21 is an adhesive surface. In the figure, the sensor sheet 11 is shown as viewed from the surface opposite to the adhesive surface of the adhesive sheet 21.

筋電図電極22は、飲食物を嚥下するときの筋活動に伴って生じる電圧を検出して電圧の経時変化(測定パラメータ)を示す電気信号S1を出力する。このセンサシート11では、図2に示すように、一例として、2つの筋電図電極22を1組とする4組(合計8つ:筋電図電極22a,22bの組、筋電図電極22c,22dの組、筋電図電極22e,22fの組、筋電図電極22g,22hの組の4組)の筋電図電極22を備えている。また、各筋電図電極22は、粘着シート21の粘着面に固定されている。 The electromyogram electrode 22 detects a voltage generated by muscle activity when swallowing food and drink, and outputs an electric signal S1 indicating a change with time (measurement parameter) of the voltage. In this sensor sheet 11, as shown in FIG. 2, as an example, four sets of two EMG electrodes 22 as one set (a total of eight: a set of EMG electrodes 22a and 22b, an EMG electrode 22c). , 22d, EMG electrodes 22e, 22f, EMG electrodes 22g, 22h) EMG electrodes 22 are provided. Further, each EMG electrode 22 is fixed to the adhesive surface of the adhesive sheet 21.

指標23は、図3に示すように、人(判定対象者100や後述する被験者300)の頸部Neにセンサシート11を貼付(装着)する際に位置合わせをするための指標(目印)であって、図2に示すように、粘着シート21の縁部に形成された切り欠きによって構成されている。 As shown in FIG. 3, the index 23 is an index (mark) for aligning the sensor sheet 11 when the sensor sheet 11 is attached (attached) to the neck Ne of a person (determination target person 100 or subject 300 described later). Therefore, as shown in FIG. 2, it is composed of notches formed at the edges of the adhesive sheet 21.

この場合、このセンサシート11では、図3に示すように、人の頸部Neにおける喉頭隆起Lp(喉仏)の近傍(外側)に指標23が位置するように粘着シート21を頸部Neに貼付したときに(以下、この状態を「規定状態」ともいう)、各筋電図電極22が測定対象筋(嚥下に関与する筋)に対向するように、各筋電図電極22および指標23の位置が規定されている。具体的には、規定状態において、測定対象筋の一例としての胸骨甲状筋、胸骨舌骨筋、甲状舌骨筋およびオトガイ筋(いずれも図示を省略する)の各近傍に各筋電図電極22が対向するように構成されている。より具体的には、規定状態において、筋電図電極22a,22bが胸骨甲状筋の近傍に対向し、筋電図電極22c,22dが胸骨舌骨筋の近傍に対向し、筋電図電極22e,22fが甲状舌骨筋の近傍に対向し、筋電図電極22g,22hがオトガイ筋の近傍に対向するように構成されている。 In this case, in this sensor sheet 11, as shown in FIG. 3, the adhesive sheet 21 is attached to the neck Ne so that the index 23 is located near (outside) the Adam's apple Lp (Adam's apple) in the human neck Ne. (Hereinafter, this state is also referred to as a "specified state"), so that each EMG electrode 22 faces the measurement target muscle (muscle involved in swallowing), each EMG electrode 22 and the index 23 The position is specified. Specifically, in the specified state, each electromyogram electrode 22 is located in the vicinity of the sternothyroid muscle, the sternohyoid muscle, the thyrohyoid muscle, and the mentalis muscle (all not shown) as examples of the muscles to be measured. Are configured to face each other. More specifically, in the specified state, the EMG electrodes 22a and 22b face the vicinity of the sternothyroid muscle, the EMG electrodes 22c and 22d face the vicinity of the sternohyoid muscle, and the EMG electrodes 22e , 22f are configured to face the vicinity of the sternohyoid muscle, and the electromyogram electrodes 22g, 22h are configured to face the vicinity of the mentalis muscle.

マイクロフォン12は、センサの他の一例であって、図3に示すように、人の頸部Neにおける喉頭隆起Lpの近傍に粘着テープによって貼付され、嚥下の際に飲食物が咽頭を通過するときの音(嚥下音)を検出して、嚥下音の大きさに応じた電圧の経時変化(測定パラメータ)を示す電気信号S2を出力する。 The microphone 12 is another example of the sensor, and as shown in FIG. 3, when the food or drink passes through the pharynx when swallowing, it is attached by an adhesive tape near the laryngeal ridge Lp in the human neck Ne. (Swallowing sound) is detected, and an electric signal S2 indicating a time-dependent change (measurement parameter) of the voltage according to the loudness of the swallowing sound is output.

操作部13は、各種の操作が可能に構成され、操作に応じた操作信号を出力する。 The operation unit 13 is configured to enable various operations and outputs an operation signal according to the operation.

記憶部14は、処理部16によって実行される後述する算出処理によって算出される活動時間Ta(判定用パラメータに相当する)を記憶する。また、記憶部14は、処理部16によって実行される後述する特定処理によって特定される閾値Tt1(第1閾値に相当する)を記憶する。 The storage unit 14 stores the activity time Ta (corresponding to the determination parameter) calculated by the calculation process described later executed by the processing unit 16. Further, the storage unit 14 stores the threshold value Tt1 (corresponding to the first threshold value) specified by the specific processing described later executed by the processing unit 16.

表示部15は、処理部16の制御に従い、センサシート11から出力される電気信号S1,S2の波形の画像や、処理部16によって実行される判定処理の結果を表示する。 The display unit 15 displays images of the waveforms of the electric signals S1 and S2 output from the sensor sheet 11 and the result of the determination process executed by the processing unit 16 under the control of the processing unit 16.

処理部16は、操作部13から出力される操作信号に従って各種の処理を行う。具体的には、処理部16は、センサシート11の筋電図電極22およびマイクロフォン12からそれぞれ出力される電気信号S1,S2によって特定される電圧の経時変化(測定パラメータ)を測定する。また、処理部16は、測定した電圧の経時変化に基づき、嚥下に伴う筋活動の活動時間Ta(判定用パラメータ)を算出する算出処理を実行する。さらに、処理部16は、活動時間Taの適正範囲Tfの最小値としての閾値Tt1を特定する特定処理を実行する。また、処理部16は、判定部として機能し、適正範囲Tf内の活動時間Taを用いて嚥下機能判定(嚥下に関する事項の判定)を行う。 The processing unit 16 performs various processes according to the operation signal output from the operation unit 13. Specifically, the processing unit 16 measures the time-dependent change (measurement parameter) of the voltage specified by the electric signals S1 and S2 output from the electromyographic electrode 22 and the microphone 12 of the sensor sheet 11, respectively. Further, the processing unit 16 executes a calculation process for calculating the activity time Ta (determination parameter) of the muscle activity associated with swallowing based on the time-dependent change of the measured voltage. Further, the processing unit 16 executes a specific process for specifying the threshold value Tt1 as the minimum value of the appropriate range Tf of the activity time Ta. In addition, the processing unit 16 functions as a determination unit, and determines the swallowing function (determination of matters related to swallowing) using the activity time Ta within the appropriate range Tf.

この判定装置1では、上記したように、人の頸部Neに装着した筋電図電極22およびマイクロフォン12から出力される電気信号S1,S2によって特定される電圧の経時変化(測定パラメータ)に基づいて算出した筋活動の活動時間Ta(判定用パラメータ)を用いて判定対象者100の嚥下機能判定(嚥下に関する事項の判定)を行うことが可能となっている。また、この判定装置1では、図4に示すように、電気信号S1によって特定される電圧の電圧値Vmが基準値Vsとクロスした時刻(立ち上がり時刻、および立ち下がり時刻)を開始時刻Tsおよび終了時刻Teとして活動時間Taを算出する。 As described above, the determination device 1 is based on the time-dependent change (measurement parameter) of the voltage specified by the electric signals S1 and S2 output from the electromyographic electrode 22 and the microphone 12 attached to the human neck Ne. It is possible to determine the swallowing function (determination of matters related to swallowing) of the determination target person 100 by using the activity time Ta (determination parameter) of the muscle activity calculated as described above. Further, in this determination device 1, as shown in FIG. 4, the start time Ts and the end time (rise time and fall time) when the voltage value Vm of the voltage specified by the electric signal S1 crosses the reference value Vs are set. The activity time Ta is calculated as the time Te.

ここで、電気信号S1には、一般的にノイズ含まれているため、図5に示すように、ノイズによって短時間の間に電圧値Vmが基準値Vsをクロスして、嚥下の際に筋が活動する実際の活動時間よりも短い時間が活動時間Taとして算出されることがある。したがって、正しい判定を行うためには、活動時間Taが極端に短いとき(活動時間Taが適正範囲Tf外のとき)には、その活動時間Taを判定に用いる活動時間Taから除外する(フィルタリングを行う)必要がある。この判定装置1では、判定に用いる活動時間Taの適正範囲Tfを特定することが可能となっている。以下、この活動時間Taの適正範囲Tfの特定方法について説明する。 Here, since the electric signal S1 generally contains noise, as shown in FIG. 5, the voltage value Vm crosses the reference value Vs in a short time due to the noise, and the muscle occurs during swallowing. The time shorter than the actual activity time during which is active may be calculated as the activity time Ta. Therefore, in order to make a correct judgment, when the activity time Ta is extremely short (when the activity time Ta is outside the appropriate range Tf), the activity time Ta is excluded from the activity time Ta used for the judgment (filtering). There is a need to do. In this determination device 1, it is possible to specify the appropriate range Tf of the activity time Ta used for the determination. Hereinafter, a method for specifying the appropriate range Tf of the activity time Ta will be described.

適正範囲Tfの特定にあたって複数(この例では、19人)の被験者300を選定する。次に、各被験者300が飲食物を嚥下するときの測定パラメータを測定して記憶させる。具体的には、図3に示すように、被験者300の目標位置(この例では、喉頭隆起Lpの左側の位置)とセンサシート11の取り付け指標23とが一致するようにセンサシート11の粘着シート21を被験者300の頸部Neにおける前頸部Anに貼付する。この際に、センサシート11の筋電図電極22が被験者300の測定対象筋に対向する。また、被験者300の頸部Neにおける喉頭隆起Lpの近傍に粘着テープを用いてマイクロフォン12を貼付する。 A plurality of subjects 300 (19 in this example) are selected to specify the appropriate range Tf. Next, the measurement parameters when each subject 300 swallows food and drink are measured and stored. Specifically, as shown in FIG. 3, the adhesive sheet of the sensor sheet 11 so that the target position of the subject 300 (in this example, the position on the left side of the adam's apple Lp) and the attachment index 23 of the sensor sheet 11 coincide with each other. 21 is attached to the anterior neck An in the neck Ne of the subject 300. At this time, the electromyogram electrode 22 of the sensor sheet 11 faces the measurement target muscle of the subject 300. In addition, the microphone 12 is attached to the vicinity of the laryngeal ridge Lp in the neck Ne of the subject 300 using an adhesive tape.

次いで、操作部13を操作して、処理部16に対して処理の開始を指示し、続いて、被験者300に対して、飲食物を嚥下するよう指示する。指示に従って被験者300が飲食物を嚥下したときには、センサシート11の筋電図電極22が、嚥下の際の筋活動に伴って生じる電圧を検出して、測定パラメータとしての電圧の経時変化を示す電気信号S1を出力する。また、マイクロフォン12が、嚥下の際の嚥下音を検出して、測定パラメータとしての嚥下音の大きさに応じた電圧の経時変化を示す電気信号S2を出力する。 Next, the operation unit 13 is operated to instruct the processing unit 16 to start the processing, and then instruct the subject 300 to swallow food and drink. When the subject 300 swallows food or drink according to the instruction, the electromyographic electrode 22 of the sensor sheet 11 detects the voltage generated by the muscle activity during swallowing, and indicates the time-dependent change of the voltage as a measurement parameter. The signal S1 is output. Further, the microphone 12 detects the swallowing sound during swallowing and outputs an electric signal S2 indicating a time-dependent change in voltage according to the loudness of the swallowing sound as a measurement parameter.

次いで、処理部16は、筋電図電極22およびマイクロフォン12からそれぞれ出力された電気信号S1,S2に基づいて判定用パラメータとしての活動時間Taを算出する算出処理を実行する。この算出処理では、処理部16は、まず、電気信号S2に基づき、嚥下音が最大となった時刻(嚥下音ピーク時刻)を測定する。 Next, the processing unit 16 executes a calculation process for calculating the activity time Ta as a determination parameter based on the electric signals S1 and S2 output from the EMG electrode 22 and the microphone 12, respectively. In this calculation process, the processing unit 16 first measures the time when the swallowing sound is maximized (the peak time of the swallowing sound) based on the electric signal S2.

続いて、処理部16は、電気信号S1に基づき、電圧値Vm(筋活動に伴う電圧)を測定し、電圧値Vmが予め決められた基準値Vs以上に立ち上がる立ち上がり時刻を特定する。この場合、立ち上がり時刻が複数存在するときには、各立ち上がり時刻前において電圧値Vmが基準値Vs未満の時間が予め規定された規定時間Tr1以上連続する立ち上がり時刻を特定し、その中で最も嚥下音ピーク時刻に近い立ち上がり時刻を活動時間Taの開始時刻Tsとする。 Subsequently, the processing unit 16 measures the voltage value Vm (voltage associated with muscle activity) based on the electric signal S1 and specifies the rise time at which the voltage value Vm rises above a predetermined reference value Vs. In this case, when there are a plurality of rise times, a rise time in which the voltage value Vm is less than the reference value Vs is specified for a predetermined time Tr1 or more consecutively before each rise time, and the swallowing sound peak is the most among them. The rise time close to the time is defined as the start time Ts of the activity time Ta.

また、開始時刻Ts後において電圧値Vmが基準値Vs以上から基準値Vs未満に立ち下がって基準値Vs未満となっている時間が予め規定された規定時間Tr2以上連続する立ち下がり時刻を特定し、その中で開始時刻Tsに最も近い立ち下がり時刻を活動時間Taの終了時刻Teとする。なお、一例として、嚥下に関する反応のない時間帯における電圧値Vmの平均値の2倍の値を基準値Vsとすることができる。 Further, after the start time Ts, the time when the voltage value Vm falls from the reference value Vs or more to less than the reference value Vs and becomes less than the reference value Vs is specified as a continuous fall time of Tr2 or more for a predetermined time. , The fall time closest to the start time Ts is set as the end time Te of the activity time Ta. As an example, the reference value Vs can be twice the average value of the voltage value Vm in the time zone when there is no reaction related to swallowing.

次いで、処理部16は、開始時刻Tsおよび終了時刻Teから活動時間Taを算出する。続いて、処理部16は、活動時間Taを記憶部14に記憶させる。以上により、1人の被験者300についての1つの活動時間Taを算出する1回の算出処理が終了する。次いで、同じ被験者300について、2回目の算出処理を実行する。なお、この例では、1人の被験者300について、固さ(嚥下難易性)が異なる2種類の飲食物を嚥下させての算出処理を、飲食物毎に6回(1人の被験者300について12回)実行するものとする。したがって、この例では、19人の被験者300について合計で228回の算出処理を実行するものとする。 Next, the processing unit 16 calculates the activity time Ta from the start time Ts and the end time Te. Subsequently, the processing unit 16 stores the activity time Ta in the storage unit 14. As described above, one calculation process for calculating one activity time Ta for one subject 300 is completed. Next, the second calculation process is executed for the same subject 300. In this example, for one subject 300, the calculation process of swallowing two kinds of foods and drinks having different hardness (difficulty of swallowing) is performed 6 times for each food and drink (12 for one subject 300). (Times) It shall be executed. Therefore, in this example, it is assumed that the calculation process is executed a total of 228 times for 19 subjects 300.

続いて、全ての被験者300についての算出処理が終了したときには、操作部13を操作して、嚥下機能判定に用いるのに適した活動時間Taの適正範囲Tfの最小値である閾値Tt1を特定する特定処理の実行を処理部16に対して指示する。この特定処理では、処理部16は、複数の被験者300について算出した複数の活動時間Taを統計処理して得た処理結果に基づいて閾値Tt1を特定する。この場合、処理部16は、各活動時間Taの混合ガウス分布を求める処理を統計処理として実行する。 Subsequently, when the calculation processing for all the subjects 300 is completed, the operation unit 13 is operated to specify the threshold value Tt1 which is the minimum value of the appropriate range Tf of the activity time Ta suitable for use in determining the swallowing function. Instructs the processing unit 16 to execute the specific processing. In this specific processing, the processing unit 16 specifies the threshold value Tt1 based on the processing result obtained by statistically processing the plurality of activity time Tas calculated for the plurality of subjects 300. In this case, the processing unit 16 executes a process of obtaining a mixed Gaussian distribution of each activity time Ta as a statistical process.

ここで、混合ガウス分布を求める処理によって閾値Tt1を特定するのに際して、次の3つを仮定する。第1に、あるクラスタ中の活動時間Taの分布はあるパラメータの確率密度分布に従い、異なるクラスタはそれとは異なるパラメータの確率密度分布に由来すると仮定する。第2に、算出した全ての活動時間Taについての全体の確率密度分布は、複数の確率密度分布の重み付き線形和で表現される混合分布に従うと仮定する。第3に、全体の分布を構成する個々の分布がいずれも正規分布であると仮定する。 Here, the following three are assumed when the threshold value Tt1 is specified by the process of obtaining the mixed Gaussian distribution. First, it is assumed that the distribution of activity time Ta in one cluster follows the probability density distribution of one parameter, and that different clusters derive from the probability density distribution of different parameters. Second, it is assumed that the overall probability density distribution for all the calculated activity time Ta follows a mixture distribution represented by a weighted linear sum of the plurality of probability density distributions. Third, it is assumed that the individual distributions that make up the overall distribution are all normal distributions.

上記の3つを仮定すると、混合分布による活動時間Taの確率密度分布p(x)は、次に示す式(1)となる。
p(x)=Σ[i=1,n]ξ(x:μσ )・・・(1)
ここでxは、計算された活動時間Taに対応する確率変数であり、nは仮定する分布数(クラスタ数)である。また、Ni(x:μσ )は、平均値がμで分散がσ のi番目の分布(正規分布)である。また、ξはi番目の分布における混合比であって、確率密度分布p(x)におけるi番目の分布の寄与度合いに相当する。
また、n個の分布における混合比ξiの合計は、次に示す式(2)のように1となる。
Σ[n,i=1]ξ=1・・・(2)
Assuming the above three assumptions, the probability density distribution p (x) of the activity time Ta based on the mixture distribution is given by the following equation (1).
p (x) = Σ [i = 1, n] ξ i N i (x: μ i σ i 2 ) ... (1)
Here, x is a random variable corresponding to the calculated activity time Ta, and n is the assumed number of distributions (number of clusters). Further, Ni (x: μ i σ i 2 ) is the i-th distribution (normal distribution) having an average value of μ i and a variance of σ i 2 . Further, ξ i is a mixing ratio in the i-th distribution, and corresponds to the degree of contribution of the i-th distribution in the probability density distribution p (x).
Further, the total of the mixing ratios ξi in the n distributions is 1 as shown in the following equation (2).
Σ [n, i = 1] ξ i = 1 ... (2)

この判定装置1では、処理部16は、算出した全ての活動時間Taの全体の確率密度分布p(x)を混合ガウス分布によって表現する上記の式(1)の各パラメータ、すなわち各正規分布の平均値μ、分散σ 、混合比ξをEMアルゴリズムにより特定する。この場合、図6に示すように、上記した算出処理によって算出した各活動時間Taの確率密度分布p(x)は、式(1)の各パラメータによって同図に示す各曲線で表される。 In this determination device 1, the processing unit 16 expresses the overall probability density distribution p (x) of all the calculated activity time Tas by a mixed Gaussian distribution, that is, each parameter of the above equation (1), that is, each normal distribution. The mean value μ i , the variance σ i 2 , and the mixing ratio ξ i are specified by the EM algorithm. In this case, as shown in FIG. 6, the probability density distribution p (x) of each activity time Ta calculated by the above calculation process is represented by each curve shown in the figure by each parameter of the equation (1).

また、処理部16は、算出した全ての活動時間Taの平均値よりも十分に小さな平均値μを有する分布N(1≦k≦nth)に由来する活動時間Taを適正範囲Tf外の活動時間Taとし、それ以外の分布N(nth+1≦k≦n)に由来する活動時間Taを適正範囲Tf内の活動時間Taと判断する。つまり、処理部16は、分布N内のいずれかの活動時間Taを閾値Tt1とする。 The processing unit 16, the appropriate range Tf outside the active time Ta derived from the distribution N k than the average value of all active time Ta calculated with a sufficiently small average value μ i (1 ≦ k ≦ n th) of the active time Ta, determines activity time Ta derived other distribution n k (n th + 1 ≦ k ≦ n) and activity time Ta in the proper range Tf. That is, the processing unit 16 sets any activity time Ta in the distribution Nk as the threshold value Tt1.

具体的には、処理部16は、混合比ξが最大の分布(第1分布に相当し、以下、この分布を「分布k」ともいう)の確率密度分布p(x)と、分布kに隣接する分布であって平均値μk−1が分布kの平均値μよりも小さい分布(第2分布に相当し、以下、この分布を「分布k−1」ともいう)の確率密度分布p(x)とが等しくなる活動時間Ta(つまり次に示す式(3)を満たすx)のうちの、分布k-1の平均値μk−1と分布kの平均値μとの間のx(xに対応する活動時間Ta)を閾値Tt1とする。
ξk−1k−1(x:μk−1,σ k−1)=ξ(x:μ,σ )・・・(3)
図6に示すように、この例では、約0.25秒が閾値Tt1として特定される。この場合、この例では、適正範囲Tfの最小値としての閾値Tt1を特定し、適正範囲Tfの最大値を規定しない。したがって、この例では、閾値Tt1以上の範囲が適正範囲Tfとなる。
Specifically, the processing unit 16 has a probability density distribution p (x) and a distribution k of a distribution having the maximum mixing ratio ξ i (corresponding to the first distribution, hereinafter this distribution is also referred to as “distribution k”). Probability density of a distribution adjacent to and whose mean value μ k-1 is smaller than the mean value μ k of the distribution k (corresponding to the second distribution, hereinafter this distribution is also referred to as “distribution k-1”). Of the activity time Ta (that is, x satisfying the following equation (3)) in which the distribution p (x) is equal, the mean value μ k-1 of the distribution k -1 and the mean value μ k of the distribution k ) Let x (activity time Ta corresponding to x) between them be the threshold Tt1.
ξ k-1 N k-1 (x: μ k-1 , σ 2 k-1 ) = ξ k N k (x: μ k , σ 2 k ) ... (3)
As shown in FIG. 6, in this example, about 0.25 seconds is specified as the threshold Tt1. In this case, in this example, the threshold value Tt1 as the minimum value of the appropriate range Tf is specified, and the maximum value of the appropriate range Tf is not specified. Therefore, in this example, the range of the threshold value Tt1 or more is the appropriate range Tf.

次いで、処理部16は、特定した閾値Tt1を記憶部14に記憶させて、特定処理を終了する。 Next, the processing unit 16 stores the specified threshold value Tt1 in the storage unit 14, and ends the specific processing.

次に、判定装置1を用いて判定対象者100についての嚥下機能判定を行う方法について説明する。 Next, a method of determining the swallowing function of the determination target person 100 using the determination device 1 will be described.

まず、被験者300を対象とした測定パラメータの測定と同じ手順で、判定対象者100についての測定パラメータを測定する。次いで、処理部16に算出処理を実行させる。続いて、操作部13を操作して、処理部16に嚥下機能判定を実行させる。この嚥下機能判定では、処理部16は、記憶部14から閾値Tt1を読み出し、次いで、算出処理で算出した活動時間Taが、閾値Tt1以上であるか否か(適正範囲Tf内であるか否か)を判別する。この場合、処理部16は、活動時間Taが、閾値Tt1未満である(適正範囲Tf内ではない)と判別したときには、嚥下機能判定を行うことなく、閾値Tt1未満であること(活動時間Taを判別に用いることが不適切であること)を表示部15に表示させる。 First, the measurement parameter for the determination target person 100 is measured by the same procedure as the measurement of the measurement parameter for the subject 300. Next, the processing unit 16 is made to execute the calculation process. Subsequently, the operation unit 13 is operated to cause the processing unit 16 to execute the swallowing function determination. In this swallowing function determination, the processing unit 16 reads the threshold value Tt1 from the storage unit 14, and then whether or not the activity time Ta calculated by the calculation process is equal to or greater than the threshold value Tt1 (whether or not it is within the appropriate range Tf). ) Is determined. In this case, when the processing unit 16 determines that the activity time Ta is less than the threshold value Tt1 (not within the appropriate range Tf), the processing unit 16 does not determine the swallowing function and determines that the activity time Ta is less than the threshold value Tt1 (activity time Ta). (It is inappropriate to use for discrimination) is displayed on the display unit 15.

一方、処理部16は、活動時間Taが、閾値Tt1以上である(適正範囲Tf内である)と判別したときには、活動時間Taと予め決められた判定基準時間とを比較して、判定対象者100が嚥下機能低下者(嚥下の機能が低下している人)であるか嚥下機能正常者(嚥下の機能が低下していない人)であるかを判定する。次いで、処理部16は、判定結果を表示部15に表示させて、嚥下機能判定を終了する。 On the other hand, when the processing unit 16 determines that the activity time Ta is equal to or greater than the threshold value Tt1 (within the appropriate range Tf), the processing unit 16 compares the activity time Ta with a predetermined determination reference time to determine the person to be determined. It is determined whether 100 is a person with decreased swallowing function (a person with decreased swallowing function) or a person with normal swallowing function (a person with not decreased swallowing function). Next, the processing unit 16 displays the determination result on the display unit 15 and ends the swallowing function determination.

このように、この処理装置、判定装置1および適正範囲特定方法によれば、複数回測定した電圧の経時変化(測定パラメータ)に基づいてそれぞれ算出した複数の活動時間Ta(判定用パラメータ)を統計処理して得た処理結果に基づいて適正範囲Tfの最小値とする閾値Tt1を特定することにより、嚥下に関する事項の判定に用いる活動時間Taの適正範囲Tfを確実に特定することができる。このため、この処理装置、判定装置1および適正範囲特定方法によれば、算出した活動時間Taを嚥下に関する事項の判定に用いることが適正であるか否かを判別することが困難な状況、(例えば、適正であるか否かを判別可能な熟練者が存在しない状況)においても、適正範囲Tfと活動時間Taとを比較することで活動時間Taを嚥下に関する事項の判定に用いることが適正であるか否かを確実かつ容易に判別することができる結果、嚥下に関する事項の判定を正確に行うことができる。 As described above, according to the processing device, the determination device 1, and the appropriate range specifying method, a plurality of activity time Tas (determination parameters) calculated based on the time-dependent changes (measurement parameters) of the voltage measured a plurality of times are statistically measured. By specifying the threshold value Tt1 which is the minimum value of the appropriate range Tf based on the processing result obtained by the processing, the appropriate range Tf of the activity time Ta used for determining the matter related to swallowing can be surely specified. Therefore, according to this processing device, the determination device 1, and the appropriate range specifying method, it is difficult to determine whether or not it is appropriate to use the calculated activity time Ta for determining items related to swallowing. For example, even in a situation where there is no expert who can determine whether or not it is appropriate), it is appropriate to use the activity time Ta for determining matters related to swallowing by comparing the appropriate range Tf and the activity time Ta. As a result of being able to reliably and easily determine whether or not there is, it is possible to accurately determine matters related to swallowing.

また、この処理装置、判定装置1および適正範囲特定方法によれば、各活動時間Taの混合ガウス分布を求める処理を統計処理として実行し、処理結果としての混合ガウス分布から得た結果に基づいて閾値Tt1を特定することにより、適正範囲Tfの最小値として信頼性の高い閾値Tt1を特定することができる。 Further, according to this processing device, the determination device 1, and the appropriate range specifying method, a process of obtaining a mixed Gaussian distribution of each activity time Ta is executed as a statistical process, and based on the result obtained from the mixed Gaussian distribution as the processing result. By specifying the threshold value Tt1, it is possible to specify the highly reliable threshold value Tt1 as the minimum value of the appropriate range Tf.

また、この処理装置、判定装置1および適正範囲特定方法によれば、混合ガウス分布における混合比ξが最大の分布kの確率密度分布p(x)と、分布kに隣接する分布であって平均値μk−1が分布kの平均値μよりも小さい分布k−1の確率密度分布p(x)とが等しくなる活動時間Taのうちの、分布k−1の平均値μk−1と分布kの平均値μとの間に位置する活動時間Taを閾値Tt1として特定することにより、適正範囲Tfの最小値としてより信頼性の高い閾値Tt1を特定することができる。 Further, according to this processing device, the determination device 1, and the appropriate range specifying method, the probability density distribution p (x) having the maximum mixing ratio ξ i in the mixed Gaussian distribution and the distribution adjacent to the distribution k. The mean value μ k-1 of the distribution k-1 is equal to the probability density distribution p (x) of the distribution k-1 whose mean value μ k-1 is smaller than the mean value μ k of the distribution k. By specifying the activity time Ta located between 1 and the mean value μ k of the distribution k as the threshold value Tt1, a more reliable threshold value Tt1 can be specified as the minimum value of the appropriate range Tf.

なお、判定装置および判定方法は、上記の構成および方法に限定されず、適宜変更することができる。例えば、混合ガウス分布における混合比ξが最大の分布kの確率密度分布p(x)と、分布kに隣接する分布であって平均値μk−1が分布kの平均値μよりも小さい分布k−1の確率密度分布p(x)とが等しくなる活動時間Taのうちの、分布k-1の平均値μk−1と分布kの平均値μとの間の活動時間Taを閾値Tt1として特定する例について上記したが、混合ガウス分布における混合比ξが最大の分布kに隣接する分布であって平均値μk−1が分布kの平均値μよりも小さい分布k−1の平均値μk−1を閾値Tt1として特定する構成および方法を採用することもできる。 The determination device and determination method are not limited to the above configurations and methods, and can be changed as appropriate. For example, the probability density distribution p (x) having the maximum mixing ratio ξ i in the mixed Gaussian distribution and the distribution adjacent to the distribution k whose mean value μ k-1 is larger than the mean value μ k of the distribution k. Of the activity time Ta at which the probability density distribution p (x) of the small distribution k-1 is equal, the activity time Ta between the mean value μ k-1 of the distribution k -1 and the mean value μ k of the distribution k Was described above as an example of specifying as the threshold value Tt1, but the distribution in which the mixing ratio ξ i in the mixed Gaussian distribution is adjacent to the maximum distribution k and the mean value μ k-1 is smaller than the mean value μ k of the distribution k. It is also possible to adopt a configuration and method in which the mean value μ k-1 of k-1 is specified as the threshold value Tt1.

また、混合ガウス分布を求める処理以外の他の統計処理によって閾値Tt1を特定する構成および方法を採用することもできる。例えば、算出した全ての活動時間Taの平均値μと標準偏差σから計算したμ−σを閾値Tt1とする構成および方法を採用することができる。また、データ群の中のある値を対象に外れ値であるかの判断を行う手法であるSmirnov‐Grubbs(SG)Testによって閾値Tt1を特定する構成および方法を採用することもできる。 Further, a configuration and a method for specifying the threshold value Tt1 by statistical processing other than the processing for obtaining a mixed Gaussian distribution can also be adopted. For example, a configuration and method in which the threshold value Tt1 is set to μ−σ calculated from the average value μ of all the calculated activity time Tas and the standard deviation σ can be adopted. Further, it is also possible to adopt a configuration and a method of specifying the threshold value Tt1 by the Smirnov-Grubbs (SG) Test, which is a method of determining whether or not a certain value in the data group is an outlier.

また、上記の構成および方法では、適正範囲Tfの最小値としての閾値Tt1だけを特定しているが、複数の活動時間Taを統計処理して得た処理結果に基づいて適正範囲Tfの最大値としての閾値Tt2(第2閾値に相当する)を特定して、活動時間Taの値が閾値Tt1以上で閾値Tt2以下の範囲を適正範囲Tfとする構成および方法を採用することができる。この構成および方法によれば、極端に短い活動時間Taおよび極端に長い活動時間Taの双方を、判定に用いる活動時間Taから除外する(フィルタリングを行う)ことが可能な適正範囲Tfを特定することができる(適正範囲Tfを厳格に規定することができる)ため、活動時間Taを嚥下に関する事項の判定に用いることが適正であるか否かをより確実に判別することができる。 Further, in the above configuration and method, only the threshold value Tt1 as the minimum value of the appropriate range Tf is specified, but the maximum value of the appropriate range Tf is based on the processing result obtained by statistically processing a plurality of activity time Tas. A configuration and method can be adopted in which the threshold value Tt2 (corresponding to the second threshold value) is specified, and the range in which the value of the activity time Ta is equal to or greater than the threshold value Tt1 and equal to or less than the threshold value Tt2 is set as the appropriate range Tf. According to this configuration and method, an appropriate range Tf that can exclude (filter) both extremely short activity time Ta and extremely long activity time Ta from the activity time Ta used for determination is specified. (The appropriate range Tf can be strictly defined), so that it is possible to more reliably determine whether or not it is appropriate to use the activity time Ta for determining matters related to swallowing.

また、上記した閾値Tt2を特定する際に、各活動時間Taの混合ガウス分布を求める処理を統計処理として実行し、処理結果としての混合ガウス分布から得た結果に基づいて閾値Tt2を特定することができる。この構成および方法によれば、適正範囲Tfの最大値として信頼性の高い閾値Tt2を特定することができる。 Further, when specifying the threshold value Tt2 described above, a process of obtaining a mixed Gaussian distribution of each activity time Ta is executed as a statistical process, and the threshold value Tt2 is specified based on the result obtained from the mixed Gaussian distribution as the processing result. Can be done. According to this configuration and method, a highly reliable threshold value Tt2 can be specified as the maximum value of the appropriate range Tf.

また、混合ガウス分布から得た結果に基づいて閾値Tt2を特定する際に、混合ガウス分布における混合比ξが最大の分布kの確率密度分布p(x)と、分布kに隣接する分布であって平均値μk+1が分布kの平均値μよりも大きい分布k+1(第3分布に相当する)の確率密度分布p(x)とが等しくなる活動時間Taのうちの、分布kの平均値μと分布k+1の平均値μk+1との間の活動時間Taを閾値Tt2として特定することができる。この構成および方法によれば、適正範囲Tfの最大値としてより信頼性の高い閾値Tt2を特定することができる。 Also, when specifying the threshold Tt2 based on the results obtained from the Gaussian mixture, a probability density distribution p of the largest distribution k mixing ratio xi] i is the Gaussian mixture (x), a distributed adjacent to distribution k The mean of the distribution k of the activity time Ta at which the mean value μ k + 1 is equal to the probability density distribution p (x) of the distribution k + 1 (corresponding to the third distribution) that is larger than the mean value μ k of the distribution k. The activity time Ta between the value μ k and the mean value μ k + 1 of the distribution k + 1 can be specified as the threshold Tt2. According to this configuration and method, a more reliable threshold value Tt2 can be specified as the maximum value of the appropriate range Tf.

また、閾値Tt2を特定する際に、分布k+1の平均値μk+1を閾値Tt2として特定する構成および方法を採用することもできる。また、上記した閾値Tt1の特定と同様にして、混合ガウス分布を求める処理以外の他の統計処理によって閾値Tt2を特定する構成および方法を採用することもできる。例えば、算出した全ての活動時間Taの平均値μと標準偏差σから計算したμ+σを閾値Tt2とする構成および方法を採用することができる。また、Smirnov‐Grubbs(SG)Testによって閾値Tt2を特定する構成および方法を採用することもできる。 Further, when specifying the threshold value Tt2, a configuration and a method of specifying the average value μ k + 1 of the distribution k + 1 as the threshold value Tt2 can also be adopted. Further, similarly to the above-mentioned specification of the threshold value Tt1, a configuration and a method for specifying the threshold value Tt2 by statistical processing other than the processing for obtaining the mixed Gaussian distribution can also be adopted. For example, it is possible to adopt a configuration and method in which the threshold value Tt2 is set to μ + σ calculated from the average value μ of all the calculated activity time Tas and the standard deviation σ. It is also possible to adopt a configuration and method for specifying the threshold value Tt2 by the Smirnov-Grubbs (SG) Test.

また、センサシート11に配設する筋電図電極22(センサ)の数は、上記した構成で例示した8つに限定されず、任意に規定することができる。また、測定パラメータとして測定した電圧の経時変化に基づいて判定用パラメータとしての活動時間Taを算出する例について上記したが、センサとしての振動ピックアップを用いて測定パラメータとしての振動に基づいて判定用パラメータとしての振幅や加速度を算出する構成および方法を採用することもできる。 Further, the number of electromyographic electrodes 22 (sensors) arranged on the sensor sheet 11 is not limited to the eight illustrated in the above configuration, and can be arbitrarily specified. Further, the example of calculating the activity time Ta as a judgment parameter based on the time course of the voltage measured as a measurement parameter has been described above, but the judgment parameter is based on the vibration as a measurement parameter using a vibration pickup as a sensor. It is also possible to adopt a configuration and a method for calculating the amplitude and acceleration of the above.

また、判定装置1を用いて嚥下に関する事項としての嚥下機能判定を行う例について上記したが、嚥下に関する事項として、判定対象の飲食物がが嚥下容易飲食物(嚥下が容易な飲食物)であるか嚥下困難飲食物(嚥下が困難な飲食物)であるかの判定を判定装置1を用いて行うこともできる。 Further, the example of determining the swallowing function as a matter related to swallowing using the determination device 1 has been described above, but as a matter related to swallowing, the food or drink to be determined is a food or drink that is easy to swallow (food or drink that is easy to swallow). It is also possible to use the determination device 1 to determine whether the food or drink is difficult to swallow (food or drink that is difficult to swallow).

1 判定装置
12 マイクロフォン
14 記憶部
16 処理部
22 筋電図電極
100 判定対象者
300 被験者
k 分布
k−1 分布
k+1 分布
Ne 頸部
p(x) 確率密度分布
Ta 活動時間
Tf 適正範囲
Tt1 閾値
Tt2 閾値
μ平均値
ξ混合比
1 Judgment device 12 Microphone 14 Storage unit 16 Processing unit 22 EMG electrode 100 Judgment target person 300 Subject k distribution k-1 distribution k + 1 distribution Ne Neck p (x) Probability density distribution Ta Activity time Tf Appropriate range Tt1 Threshold Tt2 Threshold μ i mean value ξ i mixing ratio

Claims (7)

人の頸部にセンサを装着して測定した測定パラメータに基づいて算出した判定用パラメータの値が適正範囲内のときの当該判定用パラメータを用いて嚥下に関する事項を判定する際の当該適正範囲を特定する処理部を備え、
前記処理部は、複数回測定した前記測定パラメータに基づいてそれぞれ算出した複数の前記判定用パラメータを統計処理して得た処理結果に基づいて前記適正範囲の最小値とする第1閾値を特定し、前記判定用パラメータの値が前記第1閾値以上の範囲を前記適正範囲として特定する処理装置であって、
前記処理部は、前記各判定用パラメータの混合ガウス分布を求める処理を前記統計処理として実行し、前記処理結果としての前記混合ガウス分布から得た結果に基づいて前記第1閾値を特定する処理装置。
When the value of the judgment parameter calculated based on the measurement parameter measured by attaching the sensor to the human neck is within the proper range, the appropriate range when judging matters related to swallowing using the judgment parameter Equipped with a processing unit to identify
The processing unit specifies a first threshold value as the minimum value of the appropriate range based on the processing result obtained by statistically processing a plurality of the determination parameters calculated based on the measurement parameters measured a plurality of times. A processing device that specifies a range in which the value of the determination parameter is equal to or greater than the first threshold value as the appropriate range .
The processing unit executes a process of obtaining a mixed Gaussian distribution of each determination parameter as the statistical process, and specifies the first threshold value based on the result obtained from the mixed Gaussian distribution as the processing result. ..
前記処理部は、前記混合ガウス分布における混合比が最大の第1分布の確率密度分布と、前記第1分布に隣接する分布であって平均値が第1分布の平均値よりも小さい第2分布の確率密度分布とが等しい前記判定用パラメータの値のうちの当該第1分布の平均値と当該第2分布の平均値との間に位置する当該判定用パラメータの値を前記第1閾値として特定する請求項記載の処理装置。 In the processing unit, the probability density distribution of the first distribution having the maximum mixing ratio in the mixed Gaussian distribution and the second distribution which is adjacent to the first distribution and whose mean value is smaller than the mean value of the first distribution. Of the values of the determination parameters having the same probability density distribution of, the value of the determination parameter located between the average value of the first distribution and the average value of the second distribution is specified as the first threshold value. The processing apparatus according to claim 1 . 前記処理部は、前記処理結果に基づいて前記適正範囲の最大値とする第2閾値を特定し、前記判定用パラメータの値が前記第1閾値以上で前記第2閾値以下の範囲を前記適正範囲として特定する請求項1または2記載の処理装置。 The processing unit specifies a second threshold value as the maximum value of the appropriate range based on the processing result, and the appropriate range is a range in which the value of the determination parameter is equal to or greater than the first threshold value and equal to or less than the second threshold value. The processing apparatus according to claim 1 or 2, which is specified as. 前記処理部は、前記各判定用パラメータの混合ガウス分布を求める処理を前記統計処理として実行し、前記処理結果としての前記混合ガウス分布から得た結果に基づいて前記第2閾値を特定する請求項記載の処理装置。 A claim in which the processing unit executes a process of obtaining a mixed Gaussian distribution of each determination parameter as the statistical process, and specifies the second threshold value based on the result obtained from the mixed Gaussian distribution as the processing result. 3. The processing apparatus according to 3 . 前記処理部は、前記混合ガウス分布における混合比が最大の第1分布の確率密度分布の値と、前記第1分布に隣接する分布であって平均値が第1分布の平均値よりも大きい第3分布の確率密度分布の値とが等しい前記判定用パラメータの値のうちの当該第1分布の平均値と当該第3分布の平均値との間に位置する当該判定用パラメータの値を前記第2閾値として特定する請求項記載の処理装置。 In the processing unit, the value of the probability density distribution of the first distribution having the maximum mixing ratio in the mixed Gaussian distribution and the second distribution adjacent to the first distribution whose average value is larger than the average value of the first distribution. Probability of 3 distributions The value of the judgment parameter located between the mean value of the first distribution and the mean value of the third distribution among the values of the judgment parameters that are equal to the value of the density distribution is the first. 2. The processing apparatus according to claim 4, which is specified as a threshold value. 請求項1からのいずれかに記載の処理装置と、当該処理装置によって特定された前記適正範囲内の前記判定用パラメータを用いて前記嚥下に関する事項を判定する判定部とを備えている判定装置。 A determination device including the processing device according to any one of claims 1 to 5 and a determination unit that determines a matter related to swallowing using the determination parameter within the appropriate range specified by the processing device. .. 人の頸部にセンサを装着して測定した測定パラメータに基づいて算出した判定用パラメータの値が適正範囲内のときの当該判定用パラメータを用いて嚥下に関する事項を判定する際の当該適正範囲を請求項1から5のいずれかに記載の処理装置における前記処理部に特定させる適正範囲特定方法であって、
前記処理部に対して、複数回測定した前記測定パラメータに基づいてそれぞれ算出した複数の前記判定用パラメータを統計処理して得た処理結果に基づいて前記適正範囲の最小値とする第1閾値を特定させると共に前記判定用パラメータの値が前記第1閾値以上の範囲を前記適正範囲として特定させ、かつ前記各判定用パラメータの混合ガウス分布を求める処理を前記統計処理として実行させると共に前記処理結果としての前記混合ガウス分布から得た結果に基づいて前記第1閾値を特定させる適正範囲特定方法。
When the value of the judgment parameter calculated based on the measurement parameter measured by attaching the sensor to the human neck is within the proper range, the appropriate range when judging the matter related to swallowing using the judgment parameter a appropriate range specification method Ru is specific to the processing unit in the processing apparatus according to any one of claims 1 to 5,
A first threshold value, which is the minimum value in the appropriate range, is set for the processing unit based on the processing results obtained by statistically processing a plurality of the determination parameters calculated based on the measurement parameters measured a plurality of times. A process of specifying and specifying a range in which the value of the determination parameter is equal to or greater than the first threshold value as the appropriate range and obtaining a mixed Gaussian distribution of each determination parameter is executed as the statistical process and as the processing result. An appropriate range specifying method for specifying the first threshold value based on the result obtained from the mixed Gaussian distribution .
JP2016161684A 2016-08-22 2016-08-22 Processing device, judgment device and appropriate range identification method Active JP6802467B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016161684A JP6802467B2 (en) 2016-08-22 2016-08-22 Processing device, judgment device and appropriate range identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016161684A JP6802467B2 (en) 2016-08-22 2016-08-22 Processing device, judgment device and appropriate range identification method

Publications (2)

Publication Number Publication Date
JP2018029634A JP2018029634A (en) 2018-03-01
JP6802467B2 true JP6802467B2 (en) 2020-12-16

Family

ID=61303864

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016161684A Active JP6802467B2 (en) 2016-08-22 2016-08-22 Processing device, judgment device and appropriate range identification method

Country Status (1)

Country Link
JP (1) JP6802467B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2019225243A1 (en) * 2018-05-23 2021-07-01 パナソニックIpマネジメント株式会社 Eating and swallowing function evaluation method, program, eating and swallowing function evaluation device and eating and swallowing function evaluation system
JP6802470B1 (en) * 2019-08-28 2020-12-16 国立大学法人信州大学 Judgment device
JP6996722B1 (en) 2020-10-08 2022-01-17 株式会社スキノス Detection device and detection method
WO2023176112A1 (en) * 2022-03-18 2023-09-21 国立大学法人大阪大学 Program, information processing method, and information processing device

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3928463B2 (en) * 2002-04-12 2007-06-13 松下電工株式会社 Sleep motion detector
JP5022782B2 (en) * 2007-06-05 2012-09-12 株式会社日立製作所 Biopsy device
DK2667777T3 (en) * 2011-01-28 2016-10-03 Nestec Sa APPARATUS FOR DIAGNOSIS OF SAG dysfunction
JP5292505B1 (en) * 2012-11-20 2013-09-18 株式会社東芝 Behavior estimation device, threshold calculation device, behavior estimation method, behavior estimation program
JP6073709B2 (en) * 2013-03-01 2017-02-01 長野県 Sensor sheet, swallowing activity measuring device, and swallowing activity measuring method
EP3284407B1 (en) * 2013-08-26 2018-12-05 Hyogo College Of Medicine Swallowing estimation device, information terminal device, and program

Also Published As

Publication number Publication date
JP2018029634A (en) 2018-03-01

Similar Documents

Publication Publication Date Title
JP6802467B2 (en) Processing device, judgment device and appropriate range identification method
US7427270B2 (en) Sleep stage determination apparatus
US6011990A (en) Method and device for evaluating an EEG carried out in the context of anaesthesia or intensive care
US10039465B2 (en) Method for automatically evaluating an absence EEG, computer program and evaluating device therefor
Churchland et al. Shifts in the population response in the middle temporal visual area parallel perceptual and motor illusions produced by apparent motion
US8271077B1 (en) Learning optimization using biofeedback
JP7322227B2 (en) detector
US20170007147A1 (en) Cognitive function evaluation apparatus, method, system and program
US8280498B2 (en) Method, ventilator and control unit for determining a position of an esophageal catheter
EP2734104B1 (en) System for detecting the onset of seizures
JP3658580B2 (en) Sleep depth estimating device and futon provided with the same
Charupanit et al. A simple statistical method for the automatic detection of ripples in human intracranial EEG
JP2012139442A (en) Physiological evaluation apparatus for food texture
Dalebout et al. Reliability of the mismatch negativity in the responses of individual listeners
JP2011167398A (en) Biological state determining apparatus
JP5831622B2 (en) Arousal level determination device, arousal level determination program, and arousal level determination method
JP6073709B2 (en) Sensor sheet, swallowing activity measuring device, and swallowing activity measuring method
EP2989979B1 (en) Respiratory waveform recognition method and system
KR20210009708A (en) Method, apparatus and system for mornitoring eating
TW202108077A (en) Focus automatic detection method and system
JP6692351B2 (en) Method for determining the perceptual ability of a subject
JP5609985B2 (en) Arousal level determination device, arousal level determination method, and arousal level determination program
JP2018134125A (en) Processing device and measuring device
JP2018064759A (en) Arousal state determination device, arousal state determination method, and program
JP2021105878A (en) Mastication support system

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20190204

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20190204

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190306

RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7426

Effective date: 20190306

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190805

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200721

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200804

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201002

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: 20201020

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20201102

R150 Certificate of patent or registration of utility model

Ref document number: 6802467

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250