JP2008128781A - 装着式温度測定装置および体温推定方法 - Google Patents

Abstract

【課題】 時系列の体表温度の検出データに基づいて精度良く体温を推定することができる装着式温度測定装置および体温推定方法を提供する。
【解決手段】 センサ装置１には体表温度検出部３、補助体表温度検出部４および外気温度検出部５を設ける。そして、センサ装置１のコントロールユニット７は、各温度検出部３〜５を用いて被測定者Ｈの就寝中に検出した温度検出データＤ1〜Ｄ3を外部の処理装置１３に転送する。そして、処理装置１３は、温度検出データＤ1〜Ｄ3を時系列に並べた温度データ群Ｄg1〜Ｄg3と共に、一緒に測定した体温データＤ4を用いて、ＰＬＳ回帰分析を行い、体温を推定する逆演算モデルＢを構築する。これにより、処理装置１３は、逆演算モデルＢを用いて時系列の温度データ群Ｄg1〜Ｄg3から体温を推定することができる。
【選択図】 図６

Description

本発明は、例えば就寝中に測定した被測定者の体表温度に基づいて口中温度等の体温を推定する装着式温度測定装置および体温推定方法に関する。

従来、妊娠や避妊の目安となる周期を知る目的で女性のホルモンバランス状況を把握するために、基礎体温の日毎の変動を数ヶ月間に渡って追跡記録し、その変動周期や変動幅を数値的に分析する方法が行われてきた。基礎体温は一日のうちで最も新陳代謝の少ない就寝中の体温を計測することが望ましいが、就寝中に体温測定を行うことは一般的に困難である。このため、例えば起床時の口中温度を体温として計測し、この計測した体温を基礎体温の代わりに用いてきた（例えば特許文献１参照）。しかし、女性の生活の多様化を鑑みると、毎朝、計測時刻や計測方法等の条件を一定とすることは難しく、長期間に亘る計測の継続には苦痛を伴うと共に、短時間の口中温度の計測では計測結果に揺らぎが生じ易かった。このため、就寝している被測定者の腹部等に身体表面の温度（体表温度）を測定するセンサ装置を取り付け、このセンサ装置によって長時間に亘って自動測定した体表温度に基づいて基礎体温に対応した当日の代表温度を決めるものが知られている（例えば特許文献２参照）。このとき、代表温度は、例えば測定した体表温度の最高温度や平均温度としていた。

特開平１１−３１６１６１号公報 特開２００４−１６３３９１号公報

ところで、従来技術では、被測定者の寝返り等によって体表に対するセンサ装置の密着状態が変化し、測定結果にばらつきが生じることがある。例えば、被測定者の腹部にセンサ装置を取付けた状態で、仰向けに寝たときにはセンサ装置の密着度が低いのに対し、うつぶせに寝たときにはセンサ装置の密着度が高くなる。さらに、うつぶせに寝た場合には、腹部の圧迫によって血流が増加し、この血流の増加に伴って体表温度が上昇する可能性や、電気敷き毛布の利用に伴って体表温度が上昇する可能性もある。この結果、被測定者が一時的にうつぶせに寝たときでも、このときの体表温度の測定結果が上昇してしまい、体表温度および代表温度の測定精度が低下する傾向がある。

また、従来技術では、測定した体表温度の範囲で代表温度を決定している。このとき、従来技術でも、代表温度を用いて基礎体温の変動を把握することはある程度は可能である。しかし、一般的に体表温度は外気や血流等の影響を受け易いため、実測値の最高温度や平均温度を代表温度として決定すると、代表温度には外気等の影響による揺らぎが生じる可能性がある。ここで、高温期と低温期からなる基礎体温の二相性の把握においては、低温期から高温期への＋０．３℃の上昇を的確に捉える必要がある。これに対し、従来技術では、高温期での体温の上昇か、計測の揺らぎによる変化かを区別することが難しいという問題がある。

本発明は上述した従来技術の問題に鑑みなされたもので、本発明の目的は、時系列の体表温度の検出データに基づいて精度良く体温を推定することができる装着式温度測定装置および体温推定方法を提供することにある。

上述した課題を解決するために請求項１の発明は、被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定する装着式温度測定装置であって、既知の被測定者について開始時刻から終了時刻までの間に亘って体表温度データ群を測定すると共にこのときの体温を測定し、この測定した体表温度データ群と体温とを用いて、ＰＬＳ回帰分析によって、得られた体表温度データ群から体温を逆演算することができる逆演算モデルを構築し、未知の被測定者について開始時刻から終了時刻までの間に亘って体表温度データ群を測定し、未知の被測定者について得られた体表温度データ群と前記逆演算モデルとに基づいて体温を予測する構成としたことにある。

請求項２の発明では、被測定者の身体に常時装着可能なケーシングと、該ケーシングのうち被測定者の体表面側に設けられ体表の温度を検出する体表温度検出手段と、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で該体表温度検出手段による体表温度検出データを読込む読込み手段と、該読込み手段による複数の体表温度検出データを体表温度データ群として記憶する記憶手段と、該記憶手段に記憶した体表温度データ群と前記ＰＬＳ回帰分析によって構築した逆演算モデルとに基づいて、被測定者の体温を予測する体温予測手段とを備える構成としている。

請求項３の発明では、前記記憶手段は、前記読込み手段による体表温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として体表温度データ群を記憶し、該体表温度データ群のうち異常な値の体表温度検出データに対して、主成分分析による推定モデルを用いて正常な値の体表温度検出データを推定し、前記逆演算モデルは、該推定した体表温度検出データを含む体表温度データ群を用いて、ＰＬＳ回帰分析によって構築する構成としている。

請求項４の発明では、前記ケーシングには、被測定者の体表面とは反対側に位置して外気の温度を検出する外気温度検出手段を設け、前記読込み手段は、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で前記体表温度検出データと一緒に該外気温度検出手段による外気温度検出データを読込み、前記記憶手段は、体表温度データ群に加えて、前記読込み手段による複数の外気温度検出データを外気温度データ群として記憶し、前記逆演算モデルは、既知の被測定者から測定した体表温度データ群と体温とに加えて外気温度データ群を用いて、ＰＬＳ回帰分析によって、得られた体表温度データ群および外気温度データ群から体温を逆演算することができる構成とし、前記体温予測手段は、該記憶手段に記憶した体表温度データ群および外気温度データ群と前記ＰＬＳ回帰分析によって構築した逆演算モデルとに基づいて、被測定者の体温を予測する構成としている。

請求項５の発明では、前記記憶手段は、前記読込み手段による外気温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として外気温度データ群を記憶し、該外気温度データ群のうち異常な値の外気温度検出データに対して、主成分分析による推定モデルを用いて正常な値の外気温度検出データを推定し、前記逆演算モデルは、該推定した外気温度検出データを含む外気温度データ群を用いて、ＰＬＳ回帰分析によって構築する構成としている。

請求項６の発明では、前記ケーシングには、被測定者の体表面側のうち外周縁側に位置して被測定者の体表の温度を検出する補助体表温度検出手段を設け、前記読込み手段は、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で前記体表温度検出データと一緒に該補助体表温度検出手段による補助体表温度検出データを読込み、前記記憶手段は、体表温度データ群に加えて、前記読込み手段による複数の補助体表温度検出データを補助体表温度データ群として記憶し、前記逆演算モデルは、既知の被測定者から測定した体表温度データ群と体温とに加えて補助体表温度データ群を用いて、ＰＬＳ回帰分析によって、得られた体表温度データ群および補助体表温度データ群から体温を逆演算することができる構成とし、前記体温予測手段は、該記憶手段に記憶した体表温度データ群および補助体表温度データ群と前記ＰＬＳ回帰分析によって構築した逆演算モデルとに基づいて、被測定者の体温を予測する構成としている。

請求項７の発明では、前記記憶手段は、前記読込み手段による補助体表温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として補助体表温度データ群を記憶し、該補助体表温度データ群のうち異常な値の補助体表温度検出データに対して、主成分分析による推定モデルを用いて正常な値の補助体表温度検出データを推定し、前記逆演算モデルは、該推定した補助体表温度検出データを含む補助体表温度データ群を用いて、ＰＬＳ回帰分析によって構築する構成としている。

請求項８の発明では、前記ケーシングには、前記記憶手段に記憶したデータを外部の処理装置に転送する転送手段を設け、前記逆演算モデルは、外部の処理装置を用いて構築する構成としている。

請求項９の発明では、前記逆演算モデルを用いて予測した被測定者の体温を複数日に亘って記憶する体温記憶手段を設け、該体温記憶手段に記憶した複数日の体温に対して高周波な変動成分を除去するローパスフィルタ手段を設ける構成としている。

請求項１０の発明では、被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定する体温推定方法であって、既知の被測定者について開始時刻から終了時刻までの間に亘って体表温度データ群を測定すると共にこのときの体温を測定し、この測定した体表温度データ群と体温とを用いて、ＰＬＳ回帰分析によって、得られた体表温度データ群から体温を逆演算することができる逆演算モデルを構築し、未知の被測定者について開始時刻から終了時刻までの間に亘って体表温度データ群を測定し、未知の被測定者について得られた体表温度データ群と前記逆演算モデルとに基づいて体温を予測する構成としている。

請求項１，１０の発明によれば、予め測定した体表温度データ群と体温とを用いてＰＬＳ（Partial Least Square）回帰分析によって逆演算モデルを構築する。ここで、ＰＬＳ回帰分析では、ｎ個の体表温度データからなる体表温度データ群（例えばｎ次元）と例えば口中温度等の体温との間の共分散（Convariance）を最大化するように各空間内で新軸の引き換えを行いながら、その低次元化（例えばｎ次元よりも低次元（ｎ＞ｍ）となるｍ次元）された新軸を用いて体表温度データ群と体温との間のモデル化を行い、逆演算モデルを構築する。このとき、低次元化された新軸は隠れ軸（Latent Variable）と呼ばれる。

ＰＬＳ回帰分析には低次元化された隠れ軸を用いるから、体表温度データ群の空間内の共線性（Multi-Collinearity）を効果的に処理することができ、ノイズに対してロバストな（強い）モデルパラメータを用いて逆演算モデルを構築することができる。このため、例えばＰＣＲ（Principal Components Regression）回帰分析等の他のモデリング技術に比べて体温の推定能力が高い。

また、ＰＬＳ回帰分析では、入力データである体表温度データ群の空間内のモデル化も同時に行う。このため、モデリング哲学として、出力となる体温の予測誤差だけでなく、キー（key）情報が一番存在している入力データの空間内のモデル化も同時に行っており、より真に近い逆演算モデルを構築することができる。

さらに、低次元化された隠れ軸を検証することによって、簡易的に変数関係およびサンプル関係を把握することができる。即ち、低次元化された隠れ軸によって、体表温度データ群に含まれる複数の体表温度検出データの中で体温の推定に対する影響が大きいものと、小さいものとを把握することができる。このため、影響の大きい時間帯の体表温度データを測定することによって、正確な体温を推定することができる。これにより、例えば体表温度検出データを記憶するメモリの容量が制限される場合でも、影響の大きい時間帯の体表温度検出データを測定するだけで、正確な体温を推定することができる。

また、逆演算モデルを用いて体表温度データ群から体温を推定するから、例えば毎日の体温に基づいて被測定者の健康管理を行う場合でも容易に利用することができる。さらに、逆演算モデルを用いて、数時間に亘って計測した体表温度検出データから体温を推定するから、口中温度等の体温を直接計測した場合に比べて、計測の揺らぎの影響を抑制することができる。

なお、体温は、口中温度に限らず、例えば脇下の温度でもよく、外気と隔絶された理想的な状態であれば皮膚の温度でもよい。さらに、体温は、計測が可能であれば体内温度であることが望ましい。

請求項２の発明によれば、ケーシングのうち被測定者の体表面側に体表温度検出手段を設けたから、読込み手段は、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で該体表温度検出手段を用いて体表温度検出データを読込む。そして、記憶手段は該読込み手段による複数の体表温度検出データを時系列の体表温度データ群として記憶するから、体温予測手段は、該記憶手段に記憶した体表温度データ群とＰＬＳ回帰分析によって構築した逆演算モデルとに基づいて、被測定者の体温を予測することができる。

請求項３の発明によれば、記憶手段は、読込み手段による体表温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として体表温度データ群を記憶する。一般に人体の体表温度は３４℃程度であり、この温度付近まで上昇しないとケーシング等が十分に温められた熱平衡状態にならない。また、人体の体表温度は４０℃を超えることはない。このため、熱平衡状態で体表温度検出手段によって検出される温度が例えば３２℃〜４０℃の許容範囲を逸脱することはないものと考えられる。体表温度検出データがこの温度の許容範囲以外の値となった場合には、例えば被測定者が身体からケーシングを取外したとき等が考えられる。

また、就寝中は被測定者の体表温度は比較的安定しているから、例えば１０分程度の測定間隔で±１℃の許容範囲を超えて変化することはないものと考えられる。体表温度検出データがこの温度変化の許容範囲以外の値となった場合には、例えば被測定者が一時的に起きたときや、被測定者が寝具を跳ね除けたとき等が考えられる。

このように温度の許容範囲以外または温度変化の許容範囲以外の体表温度検出データは、体温の推定値に対する寄与率は低いものと考えられる。また、温度および温度変化の許容範囲以外の体表温度検出データを用いて逆演算モデルと構築すると、体温の推定値の精度が低下する可能性がある。このため、記憶手段は、体表温度検出データのうち許容範囲以外の値を異常な値として体表温度データ群を記憶する。

一方、ＰＬＳ回帰分析を行うためには、体表温度データ群に含まれる体表温度検出データは全て揃う必要がある。このため、本発明では、既知の被測定者について、測定した体表温度データ群のうち異常な値の体表温度検出データに対して、主成分分析（Principal Component Analysis）による推定モデルを用いて正常な値の体表温度検出データを推定する。そして、この推定した体表温度検出データを含む体表温度データ群を用いて、ＰＬＳ回帰分析によって逆演算モデルを構築する。この結果、体表温度データ群に異常な値の体表温度検出データを含むときでも、異常な値の影響を受けることなく逆演算モデルを構築することができ、体温の推定値の精度を高めることができる。

また、異常な値の体表温度検出データは、容量の小さなデータ（例えばnull data等）で代用することや削除することが可能である。このため、体表温度データ群の容量を小さくすることができるから、記憶手段で記憶可能な期間を延長することができると共に、記憶したデータを外部に転送する際には、転送するデータの容量を削減でき、データの転送時間を短縮することができる。

請求項４の発明によれば、ケーシングのうち被測定者の体表面とは反対側に外気温度検出手段を設けたから、読込み手段は、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で体表温度検出データと一緒に外気温度検出手段による外気温度検出データを読込む。そして、記憶手段は体表温度データ群に加えて該読込み手段による複数の外気温度検出データを時系列の外気温度データ群として記憶するから、既知の被測定者についての体表温度データ群、外気温度データ群および体温を用いて、ＰＬＳ回帰分析によって、逆演算モデルを構築する。

このとき、逆演算モデルは、体表温度データ群に加えて、外気温度データ群を用いて構築されるから、例えば体表温度検出データが外気温度の影響を受けるときでも、この外気温度の影響を除去することができる。

ここで、体表温度検出手段および外気温度検出手段はいずれもケーシングに設けられているから、これらの温度検出データは相互に大きな相関性を有する。このように大きな相関性を有する２つの温度検出データに対して、各温度検出データが独立していることを前提とする分析（例えば重回帰分析（Multilinear Regression））を行ってモデルを構築すると、ノイズの影響を受け易くなり、信頼性のあるモデルが構築できない。

これに対し、本発明では、ＰＬＳ回帰分析を用いて逆演算モデルを構築するから、体表温度データ群、外気温度データ群と体温との間で共分散が最大化し、かつ低次元化しながら、逆演算モデルの構築を行う。このため、体表温度データ群と外気温度データ群との間で大きな相関性を有する場合でも、ノイズに対してロバストな逆演算モデルを構築することができる。

請求項５の発明によれば、記憶手段は、読込み手段による外気温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として外気温度データ群を記憶する。このとき、外気温度検出手段もケーシングに設けられているから、外気温度検出手段は、例えば布団内で被測定者の体温によって温められた空気の温度を検出すると共に、ケーシングを介して被測定者の体表温度が熱伝導して被測定者の体表温度付近まで温められる。このため、外気温度検出手段は、体表温度検出手段と比較的近い温度を検出するものと考えられる。

このように外気温度検出手段も体表温度検出手段に近い温度を検出するから、例えば３２℃〜４０℃程度の温度の許容範囲以外または±１℃／１０分程度の温度変化の許容範囲以外の外気温度検出データは、体温の推定値に対する寄与率が低いと考えられる。また、温度および温度変化の許容範囲以外の外気温度検出データを用いて逆演算モデルと構築すると、体温の推定値の精度が低下する可能性がある。このため、記憶手段は、外気温度検出データのうち許容範囲以外の値を異常な値として外気温度データ群を記憶する。

一方、ＰＬＳ回帰分析を行うためには、外気温度データ群に含まれる外気温度検出データは全て揃う必要がある。このため、本発明では、既知の被測定者について、測定した外気温度データ群のうち異常な値の外気温度検出データに対して、主成分分析（Principal Component Analysis）による推定モデルを用いて正常な値の外気温度検出データを推定する。そして、体表温度データ群、体温に加えて、この推定した外気温度検出データを含む外気温度データ群を用いて、ＰＬＳ回帰分析によって逆演算モデルを構築する。この結果、外気温度データ群に異常な値の外気温度検出データを含むときでも、異常な値の影響を受けることなく逆演算モデルを構築することができ、体温の推定値の精度を高めることができる。

また、異常な値の外気温度検出データは、容量の小さなデータで代用することや削除することが可能である。このため、例えば記憶手段で長期間に亘るデータを保存する場合やデータを外部に転送する場合には、外気温度データ群の容量を削減することができる。

請求項６の発明によれば、ケーシングには被測定者の体表面側のうち外周縁側に位置して被測定者の体表の温度を検出する補助体表温度検出手段を設けたから、読込み手段は、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で体表温度検出データと一緒に補助体表温度検出手段による補助体表温度検出データを読込む。そして、記憶手段は体表温度データ群に加えて該読込み手段による複数の補助体表温度検出データを時系列の補助体表温度データ群として記憶するから、既知の被測定者についての体表温度データ群、補助体表温度データ群および体温を用いて、ＰＬＳ回帰分析によって、逆演算モデルを構築する。

このとき、逆演算モデルは、体表温度データ群に加えて、補助体表温度データ群を用いて構築されるから、例えば体表温度検出データが外気温度の影響を受けるときでも、この外気温度の影響を除去することができる。即ち、補助体表温度検出手段は、ケーシングの外周縁側に位置しているから、例えばケーシングの中央側に位置する体表温度検出手段に比べて外気温度の影響を受け易い。このため、補助体表温度検出データは体表温度検出データに比べて外気を反映した値となるから、逆演算モデルは、補助体表温度検出データによって、外気温度の影響を相殺することができる。

ここで、体表温度検出手段および補助体表温度検出手段はいずれもケーシングに設けられているから、これらの温度検出データは相互に大きな相関性を有する。このように大きな相関性を有する２つの温度検出データに対して、各温度検出データが独立していることを前提とする分析（例えば重回帰分析）を行ってモデルを構築すると、ノイズの影響を受け易くなり、信頼性のあるモデルが構築できない。

これに対し、本発明では、ＰＬＳ回帰分析を用いて逆演算モデルを構築するから、体表温度データ群、補助体表温度データ群と体温との間で共分散が最大化し、かつ低次元化しながら、逆演算モデルの構築を行う。このため、体表温度データ群と補助体表温度データ群との間で大きな相関性を有する場合でも、ノイズに対してロバストな逆演算モデルを構築することができる。

請求項７の発明によれば、記憶手段は、読込み手段による補助体表温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として補助体表温度データ群を記憶する。このとき、補助体表温度検出手段もケーシングに設けられているから、補助体表温度検出手段は、体表温度検出手段と同様に、ケーシングのうち被測定者の体表面側に設けられているから、補助体表温度検出手段は、体表温度検出手段と比較的近い温度を検出するものと考えられる。

このように補助体表温度検出手段も体表温度検出手段に近い温度を検出するから、例えば３２℃〜４０℃程度の温度の許容範囲以外または±１℃／１０分程度の温度変化の許容範囲以外の補助体表温度検出データは、体温の推定値に対する寄与率が低いと考えられる。また、温度および温度変化の許容範囲以外の補助体表温度検出データを用いて逆演算モデルと構築すると、体温の推定値の精度が低下する可能性がある。このため、記憶手段は、補助体表温度検出データのうち許容範囲以外の値を異常な値として補助体表温度データ群を記憶する。

一方、ＰＬＳ回帰分析を行うためには、補助体表温度データ群に含まれる補助体表温度検出データは全て揃う必要がある。このため、本発明では、既知の被測定者について、測定した補助体表温度データ群のうち異常な値の補助体表温度検出データに対して、主成分分析による推定モデルを用いて正常な値の補助体表温度検出データを推定する。そして、体表温度データ群、体温に加えて、この推定した補助体表温度検出データを含む補助体表温度データ群を用いて、ＰＬＳ回帰分析によって逆演算モデルを構築する。この結果、補助体表温度データ群に異常な値の補助体表温度検出データを含むときでも、異常な値の影響を受けることなく逆演算モデルを構築することができ、体温の推定値の精度を高めることができる。

また、異常な値の補助体表温度検出データは、容量の小さなデータで代用することや削除することが可能である。このため、例えば記憶手段で長期間に亘るデータを保存する場合やデータを外部に転送する場合には、補助体表温度データ群の容量を削減することができる。

請求項８の発明によれば、ケーシングには記憶手段に記憶したデータを外部の処理装置に転送する転送手段を設け、逆演算モデルは、外部の処理装置を用いて構築する構成としている。ここで、例えば女性の基礎体温の変動周期を把握するために、就寝中の数時間の体表温度を１０分間隔で測定した体表温度データ群を用いて体温を精度良く推定する場合には、少なくとも３周期分の温度検出データを用いて逆演算モデルを構築するのが好ましい。また、逆演算モデルを特定の個人ではなく、不特定の人に適用するためには、複数名の温度検出データを用いて逆演算モデルを構築するのが好ましい。この場合、被測定者に挿着する小型のケーシング内に設けられた制御回路をなすマイクロコンピュータ等ではＰＬＳ回帰分析を用いて逆演算モデルを構築するのは難しい。

これに対し、本発明では、有線方式、無線方式等の転送手段を用いて記憶手段に記憶した体表温度検出データ等をサーバコンピュータ等の外部の処理装置に転送する。このため、処理装置は、大量の温度検出データを演算処理し、逆演算モデルを速やかに構築することができる。

なお、構築した逆演算モデルは、そのまま外部の処理装置に設けてもよい。この場合、例えば未知の被測定者は、起床後に転送手段を用いて記憶手段に記憶した体表温度データ等を外部の処理装置に転送する。これにより、外部の処理装置は、逆演算モデルを用いて当日の体温を推定し、この推定結果を被測定者がアクセス可能なホームページ内に蓄積する。この結果、被測定者は、外部の処理装置に作成されたホームページにアクセスすることによって、当日の体温を確認することができる。また、外部の処理装置は、予め登録された被測定者の電子メールのアドレスに向けて、推定した体温を送信してもよい。

さらに、構築した逆演算モデルは、ケーシング内の制御回路等に設けてもよい。この場合、ケーシング内の制御回路では、単純な行列計算を行うだけで足りるため、小型の制御回路でも体温の推定演算を行うことができる。この結果、制御回路は、ケーシングに設けられた表示器に当日の体温の推定結果を表示することができるから、被測定者は体温を速やかに確認することができる。

請求項９の発明によれば、体温記憶手段に記憶した複数日の体温に対して高周波な変動成分を除去するローパスフィルタ手段を設ける構成としている。このため、逆演算モデルを用いて推定した体温に誤差が含まれるときでも、ローパスフィルタ手段は、この誤差に基づく体温の高周波な変動成分を除去することができる。この結果、例えば女性の基礎体温は略１ヶ月周期で変動するのに対し、このような低周波な体温の変動を明確化することができ、高温期と低温期とが入れ替わる時期を精度良く把握することができる。

なお、体温記憶手段およびローパスフィルタ手段は、逆演算モデルと同様に、例えば外部の処理装置に設ける構成としてもよく、ケーシング内の制御回路等に設ける構成としてもよい。

以下、本発明の実施の形態による装着式温度測定装置を添付図面に従って詳細に説明する。

まず、図１ないし図１１は第１の実施の形態を示し、図において、１は後述する外部の処理装置１３と一緒に装着式温度測定装置を構成するウエアラブルセンサ装置（以下、センサ装置１という）で、該センサ装置１は、後述のケーシング２、体表温度検出部３、補助体表温度検出部４、外気温度検出部５、コントロールユニット７、表示部１１等によって大略構成されている。

２はセンサ装置１の本体（主支持体）をなすケーシングで、該ケーシング２は、図１ないし図４に示すように例えば樹脂材料を用いて略長円形（小判形）の箱形状に形成され、その内部には後述の回路基板６等が収容されている。また、ケーシング２のうち被測定者Ｈの体表面側（正面側）には、２個の略円形状に貫通した開口２Ａ，２Ｂが設けられている。ここで、開口２Ａはケーシング２の中央側に位置し、開口２Ｂはケーシング２の外周縁側に位置している。

３は体表の温度を検出する体表温度検出部（体表温度検出手段）で、該体表温度検出部３は、図２に示すように、ケーシング２の中央側に位置する開口２Ａに取付けられている。また、体表温度検出部３は、例えばサーミスタ等からなる測温素子と、該測温素子を覆う金属製のカバーとによって構成されている。そして、体表温度検出部３はケーシング２のうち被測定者Ｈの体表面側に配置されるから、被測定者Ｈの体表面にカバーが接触し、体表温度がカバーを通じて測温素子に熱伝導する。これにより、体表温度検出部３は、測温素子によって被測定者Ｈの体表温度に応じた信号を出力する。

なお、体表温度検出部３は、被測定者Ｈの体表面に接触して体表温度を直接的に検出する構成とした。しかし、本発明はこれに限らず、例えば体表温度検出部３と被測定者Ｈの体表面との間に下着等を挟む場合には、体表温度検出部３は、下着等を介して被測定者Ｈの体表温度を間接的に検出する構成としてもよい。即ち、体表温度検出部３は、被測定者Ｈの皮膚表面の絶対温度を検出する必要はなく、皮膚表面の温度に対応して変化する相対的な温度を検出すれば足りる。従って、毎日の計測条件がほぼ一定であれば、例えば下着１枚を挟んで体表温度検出部３を被測定者Ｈの体表面に密着させる構成としてもよい。

４は外気に影響を受けた体表の温度を補助的に検出する補助体表温度検出部（補助体表温度検出手段）で、該補助体表温度検出部４は、図２に示すように、ケーシング２の外周縁側に位置する開口２Ｂに取付けられている。また、補助体表温度検出部４は、体表温度検出部３とほぼ同様に、測温素子とカバーとによって構成されている。そして、被測定者Ｈの体表面にカバーが接触することによって、体表温度がカバーを通じて測温素子に熱伝導する。また、補助体表温度検出部４は、ケーシング２の外周縁側に位置しているから、ケーシング２の中央側に位置する体表温度検出部３に比べて外気温度の影響を受け易い。このため、補助体表温度検出部４は、体表温度検出部３に比べて外気温度を反映した温度を検出し、例えば外気によって体表温度よりも低下した温度に応じた信号を測温素子から出力する。なお、補助体表温度検出部４も、体表温度検出部３と同様に、下着等を挟んだ状態で被測定者Ｈの体表温度を間接的に検出する構成としてもよい。

５は外気温度を検出する外気温度検出部（外気温度検出手段）で、該外気温度検出部５は、図３および図４に示すように、ケーシング２のうち被測定者Ｈの体表面とは反対側（背面側）に位置して、後述する表示部１１の裏面に取付けられている。また、外気温度検出部５は、体表温度検出部３とほぼ同様に、測温素子とカバーとによって構成されている。但し、外気温度検出部５は、体表温度検出部３とは異なり、ケーシング２の内部に配置されている。また、外気温度検出部５と体表温度検出部３、補助体表温度検出部４との間には、後述の回路基板６が配置されている。そして、外気温度検出部５は、ケーシング２の背面側に位置しているから、ケーシング２を介して外気温度を検出し、外気温度に応じた信号を測温素子から出力する。

６はケーシング２の内部に収容された回路基板で、該回路基板６には、図５に示すようにマイクロコンピュータ等からなる読込み手段としてのコントロールユニット７が搭載されている。そして、コントロールユニット７は、その入力側が体表温度検出部３、補助体表温度検出部４および外気温度検出部５に接続されると共に、その出力側が後述の表示部１１に接続されている。また、コントロールユニット７には記憶手段として例えばＲＯＭ、ＲＡＭ等からなる記憶部８が設けられている。

ここで、記憶部８には、コントロールユニット７を作動させるプログラムと該プログラムで使用する開始時刻ｔ1、終了時刻ｔ2、時間間隔Δｔが予め格納されると共に、コントロールユニット７の作動によって後述の体表温度Ｔ1、補助体表温度Ｔ2、外気温度Ｔ3が記憶される構成となっている。

このとき、開始時刻ｔ1と終了時刻ｔ2は就寝中の時刻として例えば午前０時（ｔ1＝０：００ａｍ）と午前６時（ｔ2＝６：００ａｍ）とにそれぞれ設定され、時間間隔Δｔは例えば１０分程度の値に設定される。なお、開始時刻ｔ1、終了時刻ｔ2は、被測定者Ｈが例えば夜間勤務者であれば、その就寝の開始時刻および終了時刻に設定されるものである。また、時間間隔Δｔも、１０分に限らず、例えば３〜２０分程度の間で測定条件等に応じて適宜設定されるものである。

また、コントロールユニット７は、時刻を計時するタイマ９を有すると共に、例えば２個のボタンスイッチ１０Ａ，１０Ｂからなるスイッチ部１０が接続されている。そして、コントロールユニット７は、ケーシング２に搭載されたコイン型リチウム電池等の電源７Ａによって駆動し、スイッチ部１０を操作することによって記憶部８からプログラムを読出して作動する。これにより、コントロールユニット７は、タイマ９による時刻が開始時刻ｔ1に達すると、開始時刻ｔ1から終了時刻ｔ2までの間に亘って一定の時間間隔Δｔ毎に体表温度検出部３、補助体表温度検出部４、外気温度検出部５から検出温度に応じた体表温度検出データＤ1、補助体表温度検出データＤ2、外気温度検出データＤ3を読込む。

そして、コントロールユニット７は、これらの温度検出データＤ1〜Ｄ3を記憶部８に順次記憶する。これにより、記憶部８には、開始時刻ｔ1から終了時刻ｔ2までの複数の体表温度検出データＤ1からなる時系列の体表温度データ群Ｄg1、複数の補助体表温度検出データＤ2からなる時系列の補助体表温度データ群Ｄg2、および複数の外気温度検出データＤ3からなる時系列の外気温度データ群Ｄg3が記憶される。

また、記憶部８は、コントロールユニット７による温度検出データＤ1〜Ｄ3のうち予め設定された温度の許容範囲（例えば３２℃〜４０℃）以内の値で、かつ温度変化の許容範囲（例えば±１℃／１０分）以内の値を正常な値とし、許容範囲以外の値を異常な値として各温度データ群Ｄg1〜Ｄg3を記憶する。具体的には、正常な値は温度検出データＤ1〜Ｄ3をそのまま記憶し、異常な値は欠損データ（例えば正常な値と識別可能なエラーデータ）が記憶される。

ここで、一般に人体の体表温度は３４℃程度であり、この温度付近まで上昇しないとケーシング２等が十分に温められた熱平衡状態にならない。また、人体の体表温度は４０℃を超えることはない。このため、温度の許容範囲は、熱平衡状態で各温度検出部３〜５によって検出される温度が例えば３２℃〜４０℃の範囲を逸脱することはないとの考えに基づいて決められている。また、就寝中は被測定者の体表温度は比較的安定している。このため、温度変化の許容範囲は、例えば１０分程度の時間間隔Δｔの前，後で±１℃の許容範囲を超えて変化することはないとの考えに基づいて決められている。

なお、温度の許容範囲および温度変化の許容範囲は、例示した値に限らず、測定条件等を考慮して適宜設定してもよい。

また、温度検出データＤ1〜Ｄ3は、検出した温度Ｔ1〜Ｔ3そのものでもよく、例えば予め設定された基準温度数値Ｔ0と温度Ｔ1〜Ｔ3との差分値でもよい。この場合、基準温度数値Ｔ0は、データ容量を削減する意味から、例えば体表温度Ｔ1の平均値（例えば３４℃）等に設定するのが好ましい。

１１はケーシング２の背面側に設けられた液晶画面等からなる表示部で、該表示部１１はコントロールユニット７に接続され、例えばコントロールユニット７の駆動状態を表示する。また、表示部１１は、記憶部８に記憶した温度検出データＤ1〜Ｄ3を後述する外部の処理装置１３に転送する転送手段を構成し、スイッチ部１０を操作することによって、例えばＱＲコード（登録商標）等の２次元コードを表示する。ここで、この２次元コードは、記憶部８に記憶した１日分（１晩分）の温度検出データＤ1〜Ｄ3およびデータのＵＲＬ（Uniform Resource Locator）等の送信先アドレスの情報を含む。このため、被測定者Ｈは、２次元コードの読取り機能を有する携帯電話ＰＴ等を使用することによって、携帯電話ＰＴ内に２次元コード内の情報を読み込む。これにより、被測定者Ｈは、図６に示すように、インターネット等を介して携帯電話ＰＴを外部の処理装置１３にアクセスし、処理装置１３に向けて温度検出データＤ1〜Ｄ3を転送することができる。

なお、転送手段は、表示部１１による２次元コード表示を用いる構成とした。しかし、本発明はこれに限らず、コントロールユニット７と携帯電話ＰＴとの間を例えば有線方式または無線方式を用いて接続するコネクタ等の接続部を設け、該接続部によって転送手段を構成してもよい。

１２はケーシング２を被測定者Ｈに取付ける装着手段としてのクリップで、該クリップ１２は、ケーシング２の中央側に位置して可撓性をもって帯状をなして上部（外部）に向けて延びると共に、その先端が衣服等を挟持可能な構成となっている。このため、クリップ１２は、図１に示すように、ショーツ等の下着を挟むことによって、センサ装置１を被測定者Ｈの腹部に固定する。

なお、センサ装置１は被測定者Ｈの腹部に限らず、例えば被測定者Ｈの胸部に固定する構成としてもよい。この場合、クリップ１２は、就寝時に着用する下着等に取付けるものである。

１３は図６に示すようにセンサ装置１によって計測した温度検出データＤ1〜Ｄ3を用いて被測定者Ｈの体温Ｔ4を予測する体温予測手段としての処理装置で、該処理装置１３は、例えばサーバコンピュータ等によって構成されている。

また、処理装置１３は、既知の被測定者Ｈについて測定した温度データ群Ｄg1〜Ｄg3および体温データＤ4（口中温度）を用いて、ＰＬＳ回帰分析によって、得られた温度データ群Ｄg1〜Ｄg3から体温Ｔ4を逆演算することができる逆演算モデルＢを構築する。

具体的には、例えば複数名の既知の被測定者Ｈについて、センサ装置１を用いて、温度データ群Ｄg1〜Ｄg3を測定すると共に、起床時に当日の体温Ｔ4として口中温度を体温計等を用いて測定し、体温データＤ4として記録する。処理装置１３は、このようにして得られた温度データ群Ｄg1〜Ｄg3および体温データＤ4に対して、後述するＰＬＳ回帰分析を行い、温度データ群Ｄg1〜Ｄg3から体温Ｔ4を逆演算する逆演算モデルＢを構築する。

そして、逆演算モデルＢの構築後は、処理装置１３は、逆演算モデルＢを用いて未知の被測定者Ｈから送られてくる温度データ群Ｄg1〜Ｄg3から体温Ｔ4を予測し、体温の推定値Ｔ4′を算出する。また、処理装置１３は、算出した推定値Ｔ4′を数日〜数ヶ月間に亘って記憶（蓄積）する体温記憶手段としての体温記憶部１３Ａを備えている。さらに、処理装置１３は、体温記憶部１３Ａに記憶した体温の推定値Ｔ4′を被測定者Ｈ毎に表示可能なホームページを備えている。

このため、処理装置１３のホームページに対して被測定者Ｈが例えば携帯電話ＰＴを用いてアクセスすると、被測定者Ｈの要求に応じて、体温記憶部１３Ａに記憶した体温の推定値Ｔ4′を、携帯電話ＰＴ、コンピュータに表示する。これにより、携帯電話ＰＴ等の画面には、例えば温度データ群Ｄg1〜Ｄg3の測定日に対応した体温の推定値Ｔ4′が表示されると共に、１ヶ月間の体温の推定値Ｔ4′が表示される。

なお、被測定者Ｈは携帯電話ＰＴに限らず、各種の携帯端末、コンピュータ等を用いて処理装置１３にアクセスする構成としてもよい。また、処理装置１３は、予め登録されたメールアドレスに向けて、体温の推定値Ｔ4′を送信する構成としてもよい。

本実施の形態による装着式温度測定装置は上述のような構成を有するもので、次に、ＰＬＳ回帰分析を用いた逆演算モデルＢの構築方法について説明する。

まず、複数名の既知の被測定者Ｈ1〜Ｈnについて、センサ装置１を用いて体表温度Ｔ1、補助体表温度Ｔ2および外気温度Ｔ3を測定し、それぞれの時系列の温度データ群Ｄg1〜Ｄg3を記憶する。また、これらの温度データ群Ｄg1〜Ｄg3を測定した日の体温Ｔ4として、起床時の口中温度を体温計等を用いて測定し、記録する。そして、このような温度データ群Ｄg1〜Ｄg3および体温Ｔ4の測定を例えば３ヶ月に亘って行う。

次に、図７に示すように、被測定者Ｈ1について、各測定日毎に温度データ群Ｄg1〜Ｄg3を列方向に並べると共に、測定日順に行方向に並べる。被測定者Ｈ2〜Ｈnについても同様な作業を行い、全ての被測定者Ｈ1〜Ｈnについての温度データ群Ｄg1〜Ｄg3を用いて、時系列のデータポイント数×測定サンプル数の温度データ群行列Ｘを構成する。このとき、測定サンプル数は、被測定者Ｈ1〜Ｈnの人数と測定日との積となる。一方、体温データＤ4は、温度データ群Ｄg1〜Ｄg3の測定日毎に並べて、１（測定日の体温の数）×測定サンプル数の体温データ行列Ｙを構成する。

次に、温度データ群行列Ｘの各項目（各検出時刻毎の温度検出データＤ1〜Ｄ3）および体温データ行列Ｙの各項目（体温データＤ4）について、以下に示す数１の式を用いてスケール化（Mean Centering and Unit variance scaling）を行う。

このとき、各項目の平均値および標準偏差の具体例として、例えば午前１時４０分の体表温度検出データＤ1では、複数の測定日と複数名の被測定者に対しての平均値（例えば３５．８４５℃）および標準偏差（例えば０．６０３）となる。同様に、他の測定時刻の体表温度検出データＤ1に対する平均値および標準偏差も求めることができ、補助体表温度検出データＤ2、外気温度検出データＤ3、体温データＤ4についても、平均値および標準偏差を求めることができる。このため、これらの平均値および標準偏差を用いて、行列Ｘ，Ｙの各項目についてスケール前（Scale前）の値からスケール後（Scale後）の値を演算する。

次に、ＰＬＳ回帰分析では、スケール化した後の温度データ群行列Ｘおよび体温データ行列Ｙが以下の数２〜数４に示す関係を満足するものとする。

ここで、Ｔはスコア行列、Ｗ^*はローディングウエイト行列（Loading Weight Matrix）、ＰはＸローディング行列、ＱはＹローディング行列をそれぞれ示している。また、Ｅは温度データ群の残差行列（時系列のデータポイント数×測定サンプル数の行列）、Ｆは体温データの残差行列（１×測定サンプル数の行列）を示している。また、上付き記号（superscript）Ｔは転置行列を示している。

温度データ群行列Ｘと体温データ行列Ｙとの共分散が最大で、各隠れ軸（スコア行列Ｔの列）が互いに直交すると共に、残差行列Ｅ，Ｆが最小となるとの条件に基づいて、ローディング行列Ｐ，Ｑ，Ｗ^*を求めることができる。この結果、数４の式からスコア行列Ｔを求めることができるから、未知の被測定者Ｈの温度データ群Ｄg1〜Ｄg3からなる温度データ群行列Ｘnew（時系列のデータポイント数×１の行列）に対して、体温の推定値Ｙ′（Ｙ′＝Ｔ4′）は以下の数５の式に基づいて算出することができる。

この結果、未知の被測定者Ｈについて、温度データ群Ｄg1〜Ｄg3からなる温度データ群行列Ｘnewから体温Ｔ4を予測する逆演算モデルＢ（Ｂ＝Ｗ^*・Ｑ^T）を構築することができる。

なお、隠れ軸の数は、例えばリーブ・ワン・アウト法（leave-one-out validation）を用いて決定する。具体的には、逆演算モデルＢの構築に用いたトレーニングセットとなる温度データ群行列Ｘおよび体温データ行列Ｙのうち、１つのサンプル（例えば被測定者Ｈ1の任意の測定日の温度データ群Ｄg1〜Ｄg3および体温データＤ4）を一時的に除外して、残余のサンプルから一時的な逆演算モデルを構築する。この一時的な逆演算モデルから、除外された温度データ群Ｄg1〜Ｄg3を用いて体温を推定し、この推定値と既知の体温データＤ4との残差（推定誤差）を求める。次に、除外したサンプルをトレーニングセットに戻し、他のサンプルを除外して同様の作業を行う。そして、この作業を全てのサンプルについて繰り返し、全ての推定誤差の総和を求める。これらの推定誤差の総和の演算を、隠れ軸の数を変化させる毎に行うと、このような推定誤差の総和が最小となる隠れ軸の数が分かる。これにより、隠れ軸の数を決定することができる。

以上のように、ＰＬＳ回帰分析によって逆演算モデルＢを構築する際には、既知の被測定者Ｈ1〜Ｈnについて温度データ群Ｄg1〜Ｄg3からなる温度データ群行列Ｘは全ての温度検出データＤ1〜Ｄ3が揃う必要がある。一方、センサ装置１は、温度検出データＤ1〜Ｄ3のうち予め設定された温度の許容範囲（例えば３２℃〜４０℃）以内の値で、かつ温度変化の許容範囲（例えば±１℃／１０分）以内の値を正常な値とし、許容範囲以外の値を異常な値として各温度データ群Ｄg1〜Ｄg3を記憶する。このため、温度データ群Ｄg1〜Ｄg3には欠損データが含まれており、このままでは逆演算モデルＢを構築することができない。

そこで、本実施の形態では、温度データ群Ｄg1〜Ｄg3には異常な値の温度検出データＤ1〜Ｄ3に対して、主成分分析による推定モデルを用いて正常な値を推定し、この推定値を用いて欠損データを補完している。

このため、次に、温度データ群Ｄg1〜Ｄg3中の欠損データの補完処理について説明する。

まず、正常な値の温度検出データＤ1〜Ｄ3に対して、欠損データの割合が多過ぎると、正常な値を推定することができなくなる。このため、欠損データの割合が６０％を超えるサンプル（温度データ群Ｄg1〜Ｄg3およびこのときの体温データＤ4）は削除する。

また、各温度データ群Ｄg1〜Ｄg3のうち開始時刻ｔ1付近の時間帯および終了時刻ｔ2付近の時間帯は、欠損データが多い。この理由は、開始時刻ｔ1付近の時間帯については、センサ装置１の温度検出部３〜５が十分に温められていないと考えられる。また、終了時刻ｔ2付近の時間帯については、起床直前に伴って寝返り等が多くなるためと考えられる。このため、各温度データ群Ｄg1〜Ｄg3のうち開始時刻ｔ1付近の時間帯および終了時刻ｔ2付近の時間帯の温度検出データＤ1〜Ｄ3は削除した。

この結果、図７に示すように、逆演算モデルＢの構築および体温Ｔ4の推定値Ｙ′の算出には、体表温度データ群Ｄg1は例えば午前１時４０分（１：４０ａｍ）から午前５時３０分（５：３０ａｍ）までの体表温度検出データＤ1を用い、補助体表温度データ群Ｄg2は例えば午前１時３０分（１：３０ａｍ）から午前４時５０分（４：５０ａｍ）までの補助体表温度検出データＤ2を用い、外気温度データ群Ｄg1は例えば午前２時１０分（２：１０ａｍ）から午前３時２０分（３：２０ａｍ）までの外気温度検出データＤ3を用いた。

以上により、温度データ群Ｄg1〜Ｄg3中の欠損データの割合は５％程度となった。この状態で、図８に示す温度データ群Ｄg1〜Ｄg3中の欠損データの補完処理を行った。

まず、ステップ１では、欠損データ部分に初期値を代入し、欠損データを補った温度データ群行列Ｘ（温度データ群Ｄg1〜Ｄg3）を構成する。このとき、初期値としては、欠損データに対応する温度検出データＤ1〜Ｄ3の平均値またはＮＩＰＡＬＳ（繰り返しアルゴリズム）のＳＣＰ（Single Component Projection）による推定値を用いる。

次に、ステップ２では、欠損データを補った温度データ群行列Ｘを用いて主成分分析を行い、以下の数６に示すＰＣＡ（Principal Component Analysis）モデルを構築する。

ここで、数２の式と同様に、Ｔはスコア行列、ＰはＸローディング行列をそれぞれ示している。このとき、スコア行列Ｔ、Ｘローディング行列Ｐは、温度データ群行列Ｘの共分散が最大となると共に、各隠れ軸（スコア行列Ｔの列）が互いに直交するとの条件に基づいて、ローディング行列Ｐおよびスコア行列Ｔを求めることができる。

次に、ステップ３では、ＰＣＡモデルが構築できたら、数６に示すように、該ＰＣＡモデルを用いて温度データ群行列の推定値Ｘ′を算出する。

次に、ステップ４では、温度データ群行列Ｘと推定値Ｘ′とのうち欠損データ箇所の収束性を検討する。そして、各欠損データ箇所の推定値の変動が所定のしきい値よりも小さい場合には、欠損データの推定値を記憶して処理を終了する。

一方、各欠損データ箇所の推定値の変動が所定のしきい値よりも大きい場合には、ステップ５に移って欠損データ箇所をステップ３で算出した推定値Ｘ′のデータに置き換えて、ステップ２以降の処理を繰り返す。これにより、欠損データを適切な推定値Ｘ′で補完することができる。

以上のような欠損データの補完処理を行った後に、ＰＬＳ回帰分析を行うことによって、逆演算モデルＢを構築することができる。次に、実際に逆演算モデルＢを構築して体温の推定値Ｔ4′を演算した。

ここで、既知の６５名の女性の被測定者Ｈ1〜Ｈ65にセンサ装置１を装着してもらい、６ヶ月間に亘って温度データ群Ｄg1〜Ｄg3および体温データＤ4（口中温度）を収集した。

最初に、これら全てのサンプル（温度データ群Ｄg1〜Ｄg3および体温データＤ4）を用いて逆演算モデルＢを構築した。しかし、逆演算モデルＢの構築に用いたデータ（トレーニングデータ）を用いて推定誤差を検討したところ、約±０．３〜０．４℃の推定誤差をもつ結果となり、精度が悪くなった。これは被測定者Ｈ1〜Ｈ65毎にセンサ装置１の装着状況や口中温度の測定精度が異なり、欠損データの多い被測性者や比較的精度の低い被測性者のデータ等が外乱として作用するためと考えられる。

このため、例えば１ヶ月程度の間に高温期と低温期とを周期的に繰り返す再現性が高い体温データＤ4、即ち基礎体温の周期性が明確な体温データＤ4が得られた３名の被測定者Ｈ2，Ｈ13，Ｈ20を特定し、この３名の被測定者Ｈ2，Ｈ13，Ｈ20についての温度データ群Ｄg1〜Ｄg3および体温データＤ4を用いて逆演算モデルＢを構築した。他の被測定者のデータについては、この逆演算モデルＢを用いて体温Ｔ4を推定演算し、逆演算モデルＢの正当性のチェックに用いた。

このとき、温度データ群行列Ｘ（温度データ群Ｄg1〜Ｄg3）は、１０分間隔の時系列の温度検出データＤ1〜Ｄ3によって構成されるため、例えば５３次元のデータである。これに対し、本実施の形態では、この５３次元の温度データ群行列ＸをＰＬＳ回帰分析によって１０本の隠れ軸による新次元軸（１０次元）に引き換えることによって、温度データ群行列ＸによるＸ空間の約８８％を説明することができ、体温データ行列ＹによるＹ空間の約５０％を説明することができた。この結果、３名の被測定者Ｈ2，Ｈ13，Ｈ20によるデータを用いた逆演算モデルＢでは、トレーニングデータに対する推定誤差が約±０．２℃程度となり、推定精度を高めることができた。

図９に逆演算モデルＢを用いてトレーニングデータ（被測定者Ｈ2，Ｈ13，Ｈ20の温度データ群Ｄg1〜Ｄg3）から算出した体温の推定値Ｔ4′を示す。この結果から、逆演算モデルＢを用いて算出した体温の推定値Ｔ4′も、実際の体温Ｔ4とほぼ同様に高温期と低温期とを有する周期性をもって変化することが分かる。また、体温の推定値Ｔ4′の誤差は、略±０．２℃程度となることも分かる。

また、図１０に逆演算モデルＢを用いてバリデーションデータ（被測定者Ｈ2，Ｈ13，Ｈ20以外の被測定者の温度データ群Ｄg1〜Ｄg3）から算出した体温の推定値Ｔ4′を示す。この結果から、体温の推定値Ｔ4′の誤差は少し増加するものの、トレーニングデータ以外に対しても、逆演算モデルＢを用いて算出した体温の推定値Ｔ4′は、実際の体温Ｔ4とほぼ同様に高温期と低温期とを有する周期性をもって変化することが分かる。

図１１に逆演算モデルＢのモデル係数（スコア行列Ｔの係数）を示す。この係数は、体温の推定値Ｔ4′に対する各温度検出データ群Ｄg1〜Ｄg3の影響の度合いに対応している。この結果から、体表温度検出データＤ1は午前２時４０分と午前５時の係数が大きく、補助体表温度検出データＤ2は午前２時３０分の係数が大きく、外気温度検出データＤ3は午前２時２０分の係数が大きくなっている。これにより、午前２時２０分から午前２時４０分の時間帯における各温度検出データＤ1〜Ｄ3が支配的に寄与していることが分かる。このため、センサ装置１の記憶部８の容量が制限される場合等には、午前２〜午前３時頃の温度検出データＤ1〜Ｄ3を記憶すれば、推定精度を維持しつつ有効な逆演算モデルＢが構築できることが分かる。

かくして、本実施の形態によれば、予め測定した温度データ群Ｄg1〜Ｄg3と体温データＤ4とを用いてＰＬＳ回帰分析によって逆演算モデルＢを構築したから、温度データ群Ｄg1〜Ｄg3からなる例えば５３次元の温度データ群行列Ｘを１０次元に低次元化された隠れ軸を用いて表すことができる。このため、温度データ群行列Ｘの空間内の共線性を効果的に処理することができ、ノイズに対してロバストなモデルパラメータを用いて逆演算モデルＢを構築することができる。この結果、例えばＰＣＲ回帰分析等の他のモデリング技術に比べて体温の推定能力を高めることができる。

また、ＰＬＳ回帰分析では、入力データである温度データ群行列Ｘ（温度データ群Ｄg1〜Ｄg3）の空間内のモデル化も同時に行う。このため、モデリング哲学として、出力となる体温の予測誤差だけでなく、重要なキー情報が一番存在している温度データ群行列Ｘの空間内のモデル化も同時に行っており、より真に近い逆演算モデルＢを構築することができる。

さらに、低次元化された隠れ軸を検証することによって、簡易的に変数関係およびサンプル関係を把握することができる。即ち、低次元化された隠れ軸によって、温度データ群行列Ｘに含まれる複数の温度検出データＤ1〜Ｄ3の中で体温の推定に対する影響が大きいものと、小さいものとを把握することができる。このため、影響の大きい時間帯の温度検出データＤ1〜Ｄ3を測定することによって、正確な体温を推定することができる。これにより、例えば温度検出データＤ1〜Ｄ3を記憶する記憶部８の容量が制限される場合でも、影響の大きい時間帯の温度検出データＤ1〜Ｄ3を測定するだけで、正確な体温を推定することができる。

また、逆演算モデルＢを用いて温度データ群Ｄg1〜Ｄg3から体温Ｔ4を推定するから、例えば毎日の体温に基づいて被測定者の健康管理を行う場合でも容易に利用することができる。さらに、逆演算モデルＢを用いて、数時間に亘って計測した温度検出データＤ1〜Ｄ3から体温Ｔ4を推定するから、口中温度等の体温を直接計測した場合に比べて、計測の揺らぎの影響を抑制することができる。

また、本実施の形態では、記憶部８は、各温度検出データＤ1〜Ｄ3のうち予め設定された温度の許容範囲以内の値で、かつ温度変化の許容範囲以内の値を正常な値とし、温度または温度変化のうち少なくともいずれか一方で許容範囲以外の値を異常な値として体表温度データ群を記憶する。このため、逆演算モデルＢを構築するときに、温度の許容範囲以外または温度変化の許容範囲以外の温度検出データＤ1〜Ｄ3を用いることがなくなる。この結果、温度検出データＤ1〜Ｄ3の異常な値が逆演算モデルＢに影響することがなくなり、逆演算モデルＢによる体温の推定値Ｔ4′の精度を高めることができる。

一方、ＰＬＳ回帰分析を行うためには、温度データ群Ｄg1〜Ｄg3に含まれる温度検出データＤ1〜Ｄ3は全て揃う必要がある。このため、本実施の形態では、既知の被測定者について、測定した温度データ群Ｄg1〜Ｄg3のうち異常な値の温度検出データＤ1〜Ｄ3に対して、主成分分析による推定モデル（ＰＣＡモデル）を用いて正常な値の温度検出データＤ1〜Ｄ3を推定する。そして、この推定した温度検出データＤ1〜Ｄ3を含む温度データ群Ｄg1〜Ｄg3を用いて、ＰＬＳ回帰分析によって逆演算モデルＢを構築する。この結果、温度データ群Ｄg1〜Ｄg3に異常な値の温度検出データＤ1〜Ｄ3を含むときでも、異常な値の影響を受けることなく逆演算モデルＢを構築することができ、体温の推定値Ｔ4′の精度を高めることができる。

また、本実施の形態では、センサ装置１には、体表温度検出部３に加えて、補助体表温度検出部４および外気温度検出部５を設けたから、記憶部８には、体表温度データ群Ｄg1に加えて、補助体表温度データ群Ｄg2および外気温度データ群Ｄg3を記憶することができる。そして、既知の被測定者についての温度データ群Ｄg1〜Ｄg3および体温データＤ4を用いて、ＰＬＳ回帰分析によって、逆演算モデルＢを構築する。

このとき、各温度データ群Ｄg1〜Ｄg3に対する外気温度の影響が相互に異なるのに加え、逆演算モデルＢは、体表温度データ群Ｄg1に加えて、補助体表温度データ群Ｄg2および外気温度データ群Ｄg3を用いて構築する。このため、例えば体表温度検出データＤ1が外気温度の影響を受けるときでも、逆演算モデルＢは、この外気温度の影響を除去して、体温の推定値Ｔ4′を演算することができる。

ここで、体表温度検出部３、補助体表温度検出部４および外気温度検出部５はいずれもケーシング２に設けられているから、これらの温度検出データＤ1〜Ｄ3は相互に大きな相関性を有する。これに対し、本実施の形態では、ＰＬＳ回帰分析を用いて逆演算モデルＢを構築するから、温度データ群Ｄg1〜Ｄg3と体温データＤ4との間で共分散が最大化し、かつ低次元化しながら、逆演算モデルＢの構築を行う。このため、温度データ群Ｄg1〜Ｄg3の間で相互に大きな相関性を有する場合でも、ノイズに対してロバストな逆演算モデルＢを構築することができる。

また、本実施の形態では、センサ装置１のケーシング２には記憶部８に記憶した温度検出データＤ1〜Ｄ3を外部の処理装置１３に転送する表示部１１を設ける構成とした。このため、表示部１１に表示した２次元コードを用いて、記憶部８内の温度検出データＤ1〜Ｄ3を外部の処理装置１３に転送し、該処理装置１３を用いて、逆演算モデルＢを構築することができる。

このため、例えば複数名の被測定者について数ヶ月間に亘る温度検出データＤ1〜Ｄ3を用いて逆演算モデルＢを構築するときでも、処理装置１３は、これらの大量の温度検出データＤ1〜Ｄ3を大容量の体温記憶部１３Ａ内に蓄積し、高速に演算処理することができ、逆演算モデルＢを速やかに構築することができる。

なお、前記第１の実施の形態では、逆演算モデルＢは外部の処理装置１３内に設ける構成とした。しかし、本発明はこれに限らず、例えば逆演算モデルＢを処理装置１３で構築した後に、センサ装置１のコントロールユニット７に移植する構成としてもよい。この場合、ケーシング２内のコントロールユニット７では、単純な行列計算を行うだけで足りるため、小型のコントロールユニット７でも体温Ｔ4の推定演算を行うことができる。この結果、コントロールユニット７は、ケーシング２に設けられた表示部１１に温度検出データＤ1〜Ｄ3を測定した当日の体温Ｔ4の推定結果を表示することができるから、被測定者Ｈは体温を速やかに確認することができる。

また、前記第１の実施の形態では、センサ装置１には体表温度検出部３に加えて補助体表温度検出部４および外気温度検出部５を設け、これらによる温度検出データＤ1〜Ｄ3を用いて逆演算モデルＢを構築する構成とした。しかし、本発明はこれに限らず、例えばセンサ装置１の外気温度検出部５を省き、温度検出データＤ1，Ｄ2を用いて逆演算モデルを構築する構成としてもよい。この場合、第１の実施の形態と同じ３名の被測定者Ｈ2，Ｈ13，Ｈ20の温度検出データＤ1，Ｄ2を用いて逆演算モデルを構築したところ、温度データ群行列ＸによるＸ空間の約９１％を説明することができたものの、体温データ行列ＹによるＹ空間の約４４％しか説明できず、Ｙ空間の説明率が低下した。これに伴い、推定精度も約６％低下した。

また、センサ装置１の補助体表温度検出部４を省き、温度検出データＤ1，Ｄ3を用いて逆演算モデルを構築する構成としてもよい。この場合、第１の実施の形態と同じ３名の被測定者Ｈ2，Ｈ13，Ｈ20の温度検出データＤ1，Ｄ3を用いて逆演算モデルを構築したところ、温度データ群行列ＸによるＸ空間の約９１％を説明することができたものの、体温データ行列ＹによるＹ空間の約４３％しか説明できず、Ｙ空間の説明率が低下した。これに伴い、推定精度も約７％低下した。

さらに、センサ装置１の体表温度検出部３を省き、補助体表温度検出データＤ2と外気温度検出データＤ3とを用いて逆演算モデルを構築する構成としてもよい。この構成では、補助体表温度検出部４は、体表温度検出手段として使用する。この場合、第１の実施の形態と同じ３名の被測定者Ｈ2，Ｈ13，Ｈ20の温度検出データＤ1，Ｄ3を用いて逆演算モデルを構築したところ、温度データ群行列ＸによるＸ空間の約９２％を説明することができたものの、体温データ行列ＹによるＹ空間の約４４％しか説明できず、Ｙ空間の説明率が低下した。これに伴い、推定精度も約６％低下した。

このように、３個の温度検出部３〜５を用いた場合に比べて、２個の温度検出部を用いた場合には、いずれも推定精度が約６％程度低下する。しかし、例えば体表温度検出データＤ1のみを用いた場合と比べたときには、推定精度を高めることができ、外気温度等の影響を除去するには有効である。

また、本発明では、例えばセンサ装置１から補助体表温度検出部４および外気温度検出部５を省き、温度検出データＤ1のみを用いて逆演算モデルを構築する構成としてもよい。この場合、逆演算モデルによる推定精度は低下するものの、製造コストを低くすることができる。

また、前記第１の実施の形態では、例えば３名の被測定者の温度検出データＤ1〜Ｄ3を用いて逆演算モデルＢを構築する構成としたが、４名以上の被測定者の温度検出データを用いて逆演算モデルを構築してもよく、１名または２名の被測定者の温度検出データを用いて逆演算モデルを構築してもよい。複数名の被測定者の温度検出データを用いた場合には、逆演算モデルの汎用性が高くなる。

また、逆演算モデルＢを構築するための既知の被測定者と、逆演算モデルＢを用いて体温を推定する未知の被測定者とは、異なる人物でもよく、同一人物でもよい。即ち、既知の被測定者とは、逆演算モデルを構築するために、体表温度検出データ等に加えて体温（口中温度）を測定した被測定者を示すものである。一方、未知の被測定者とは、体温を除いて体表温度検出データ等を測定した被測定者を示すものである。

特に、逆演算モデルＢを構築するための既知の被測定者と、逆演算モデルＢを用いて体温を推定する未知の被測定者とが、１名の同一人物である場合には、汎用性は低下するものの、当該被測定者については、逆演算モデルの推定精度が向上する。このため、例えば夜間勤務者のように、他の人物と生活パターンが著しく異なる人には、その人に特有な逆演算モデルを構築するのが有効である。

次に、図１２および図１３は本発明の第２の実施の形態を示し、本実施の形態の特徴は、処理装置には、体温記憶部に記憶した複数日の体温に対して高周波な変動成分を除去するローパスフィルタ部を設ける構成としたことにある。なお、本実施の形態では、前述した第１の実施の形態と同一の構成要素に同一の符号を付し、その説明を省略するものとする。

２１は処理装置１３に設けられたローパスフィルタ部（ローパスフィルタ手段）で、該ローパスフィルタ部２１は、処理装置１３に搭載されたソフトウエア、または処理装置１３に接続されたデジタルフィルタ回路等のハードウエアによって構成されている。そして、ローパスフィルタ部２１は、体温記憶部１３Ａ内に記憶した複数日の体温の推定値Ｔ4′に対して、高温期と低温期からなる基礎体温の周期性を残しつつ、高周波な温度の変動成分を除去する処理を施す。即ち、例えば１周期が１０〜２０日以上となるような低周波（長周期）な温度の変動成分は残し、これよりも高周波（短周期）な温度の変動成分を除去する。これにより、図１３に示すように、逆演算モデルＢが１〜３日程度の間の急激に変化する体温の推定値Ｔ4′を算出するときでも、高温期と低温期とが略１ヶ月単位で繰り返す基礎体温の周期性を明確化することができる。

かくして、このように構成される第２の実施の形態でも、第１の実施の形態とほぼ同様の作用効果を得ることができる。特に、第２の実施の形態では、ローパスフィルタ部２１を用いて、体温記憶部１３Ａに記憶した複数日の体温に対して高周波な変動成分を除去する構成としたから、逆演算モデルＢを用いて推定した体温に誤差が含まれるときでも、ローパスフィルタ部２１は、この誤差に基づく体温の高周波な変動成分を除去することができる。この結果、例えば女性の基礎体温は略１ヶ月周期で変動するのに対し、このような低周波な体温の変動を明確化することができ、高温期と低温期とが入れ替わる時期を精度良く把握することができる。

次に、図１４および図１５は本発明の第３の実施の形態を示し、本実施の形態の特徴は、処理装置は、被測定者の身体状態を円グラフを用いて表示すると共に、被測定者の体温の推定値Ｔ4′を文字記号を用いた棒グラフで表示する構成としたことにある。なお、本実施の形態では、前述した第１の実施の形態と同一の構成要素に同一の符号を付し、その説明を省略するものとする。

図１４は処理装置１３によってホームページ内に表示された円グラフ３１を示している。ここで、円グラフ３１は８等分の扇形をなす第１〜第８ステージ３１Ａ〜３１Ｈに分割され、半分の第１〜第４ステージ３１Ａ〜３１Ｄは低温期に対応し、残余の第５〜第８ステージ３１Ｅ〜３１Ｈは高温期に対応している。これにより、円グラフ３１は、被測定者が低温期および高温期（第１〜第８ステージ３１Ａ〜３１Ｈ）のいずれのステージに該当するのかを表示するものである。

そして、円グラフ３１の１周分が基礎体温の１周期に該当するものとし、１サイクル（１周）は基本的に２８日としている。このとき、第１，第２，第５，第６ステージ３１Ａ，３１Ｂ，３１Ｅ，３１Ｆはそれぞれ４日が割り当てられ、第３，第４，第７，第８ステージ３１Ｃ，３１Ｄ，３１Ｇ，３１Ｈはそれぞれ３日が割り当てられる。また、第１ステージ３１Ａの初日は、被測定者が申告する月経開始連絡に基づいて決定する。

また、低温期の基準温度は、月経開始日から４〜１１日目までの平均値とする。一方、高温期の基準温度は、低温期の基準温度に対して所定の温度差（例えば０．３℃）を加えた温度とする。

そして、処理装置１３は、原則として第１ステージ３１Ａの初日から経過した日数に応じて第１〜第４ステージ３１Ａ〜３１Ｄのいずれに該当するかを判定し、該当するステージ中に例えば文字記号（例えば星印）を表示する。なお、該当するステージは、文字記号に限らず、例えば点滅等によって表示する構成としてもよい。

但し、低温期の期間には個人差があるので、１２日目以降の体温の推定値Ｔ4′が低温期の基準温度に比べて０．３℃よりも上昇した日を第５ステージ３１Ｅに突入したものと判定し、第５ステージ３１Ｅ中に文字記号を表示する。逆に、１２日目以降であっても低温期の基準温度に比べて０．３℃まで上昇しないときには、第４ステージ３１Ｄに留まっているものとして、第４ステージ３１Ｄ中に文字記号を表示する。特に、１５日目以降でも低温期の基準温度に比べて０．３℃の上昇が確認できないときには、低温期が長引いていると判断し、その旨を表示する。

また、第５ステージ３１Ｅに突入した後は、処理装置１３は、再び経過した日数に応じて第５〜第８ステージ３１Ｅ〜３１Ｈのいずれに該当するかを判定し、該当するステージ中に例えば文字記号（例えば星印）を表示する。

また、処理装置１３は、第５ステージ３１Ｅの初日から１５日以上が経過しても次回の月経連絡が無い場合は、体温の推定値Ｔ4′が低温期の基準温度よりも低下したか否かを確認する。そして、体温の推定値Ｔ4′が低温期の基準温度よりも低下した場合には、処理装置１３は、第１ステージ３１Ａに突入したものと判定し、第１ステージ３１Ａ中に文字記号を表示すると共に、月経開始の連絡を忘れている可能性がある旨を表示する。

一方、体温の推定値Ｔ4′が低温期の基準温度まで低下しておらず、高温期にあると考えられる場合には、処理装置１３は、第８ステージ３１Ｈ中に文字記号を表示すると共に、妊娠の可能性がある旨を表示する。

図１５は処理装置１３によって円グラフ３１とは別個にホームページ内に表示された棒グラフ３２を示している。ここで、棒グラフ３２は被測定者が申告する月経開始連絡に基づいて月経開始日を決定し、その日付けと一緒に月経開始を示す記号（例えば丸記号）で表示する。

また、処理装置１３は、過去３周期分の体温の推定値Ｔ4′の平均値Ｍを演算する。そして、体温の推定値Ｔ4′が平均値Ｍよりも例えば０．１℃未満の範囲で高いときには、高温期であることを示す文字記号（例えば黒星印）を１個表示する。また、体温の推定値Ｔ4′が平均値Ｍよりも例えば０．１℃以上０．２℃未満の範囲で高いときには、高温期であることを示す文字記号を２個表示する。さらに、体温の推定値Ｔ4′が平均値Ｍよりも例えば０．２℃以上高いときには、高温期であることを示す文字記号を３個表示する。

一方、体温の推定値Ｔ4′が平均値Ｍと同一または平均値Ｍよりも例えば０．１℃未満の範囲で低いときには、低温期であることを示す文字記号（例えば白星印）を１個表示する。また、体温の推定値Ｔ4′が平均値Ｍよりも例えば０．１℃以上０．２℃未満の範囲で低いときには、低温期であることを示す文字記号を２個表示する。さらに、体温の推定値Ｔ4′が平均値Ｍよりも例えば０．２℃以上低いときには、低温期であることを示す文字記号を３個表示する。

これにより、折れ線グラフによる表示を用いた場合に比べて、例えば携帯電話ＰＴの画面でも見易く表示することができる。

なお、高温期を示す文字記号を赤色、低温期を示す文字記号を青色等のように、色分け表示する構成としてもよい。

また、被測定者に当日の体調（良、悪）を申告してもらい、棒グラフ３２には、日付に対応した位置に、申告された当日の体調を簡単な文字記号を用いて表示する構成としてもよい。

かくして、このように構成される第３の実施の形態でも、第１の実施の形態とほぼ同様の作用効果を得ることができる。特に、処理装置１３は、被測定者の身体状態を示す円グラフ３１と、被測定者の体温の推定値Ｔ4′を示す棒グラフ３２とを表示する構成としたから、被測定者が円グラフ３１、棒グラフ３２を目視することによって、当日の身体状態等を容易に把握することができる。

なお、本実施の形態では、第１の実施の形態による逆演算モデルＢを用いて演算した体温の推定値Ｔ4′を用いて、円グラフ３１、棒グラフ３２を作成する構成とした。しかし、本発明はこれに限らず、例えば第２の実施の形態のように、被測定者の体温の推定値に対してローパスフィルタ処理を行った後に、円グラフ、棒グラフを作成する構成としてもよい。

本発明の第１の実施の形態によるウエアラブルセンサ装置を被測定者が装着した状態を示す外観図である。 図１中のウエアラブルセンサ装置を単体で示す正面図である。 図２中のウエアラブルセンサ装置を示す背面図である。 ウエアラブルセンサ装置を図２中の矢示IV−IV方向からみた断面図である。 ウエアラブルセンサ装置に設けられたコントロールユニットの構成を示すブロック図である。 ウエアラブルセンサ装置と外部の処理装置とから構成した装着式温度測定装置を示す説明図である。 処理装置による逆演算モデルの構築に用いる温度データ群行列Ｘおよび体温データ行列Ｙを示す説明図である。 欠損データの補完処理を示す流れ図である。 逆演算モデルを用いてトレーニングデータから算出した体温の推定値を示す特性線図である。 逆演算モデルを用いてバリデーションデータから算出した体温の推定値を示す特性線図である。 逆演算モデルのモデル係数を示す説明図である。 第２の実施の形態による装着式温度測定装置を示す説明図である。 ローパスフィルタ処理を行った体温の推定値を示す特性線図である。 第３の実施の形態による円グラフを示す説明図である。 第３の実施の形態による棒グラフを示す説明図である。

符号の説明

１ ウエアラブルセンサ装置
２ ケーシング
３ 体表温度検出部（体表温度検出手段）
４ 補助体表温度検出部（補助体表温度検出手段）
５ 外気温度検出部（外気温度検出手段）
７ コントロールユニット（データ読込み手段）
８ 記憶部（記憶手段）
１１ 表示部（転送手段）
１３ 処理装置

Claims (10)

1. 被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定する装着式温度測定装置であって、
   既知の被測定者について開始時刻から終了時刻までの間に亘って体表温度データ群を測定すると共にこのときの体温を測定し、
   この測定した体表温度データ群と体温とを用いて、ＰＬＳ回帰分析によって、得られた体表温度データ群から体温を逆演算することができる逆演算モデルを構築し、
   未知の被測定者について開始時刻から終了時刻までの間に亘って体表温度データ群を測定し、
   未知の被測定者について得られた体表温度データ群と前記逆演算モデルとに基づいて体温を予測する構成としてなる装着式温度測定装置。
2. 被測定者の身体に常時装着可能なケーシングと、該ケーシングのうち被測定者の体表面側に設けられ体表の温度を検出する体表温度検出手段と、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で該体表温度検出手段による体表温度検出データを読込む読込み手段と、該読込み手段による複数の体表温度検出データを体表温度データ群として記憶する記憶手段と、該記憶手段に記憶した体表温度データ群と前記ＰＬＳ回帰分析によって構築した逆演算モデルとに基づいて、被測定者の体温を予測する体温予測手段とを備える構成としてなる請求項１に記載の装着式温度測定装置。
3. 前記記憶手段は、前記読込み手段による体表温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として体表温度データ群を記憶し、
   該体表温度データ群のうち異常な値の体表温度検出データに対して、主成分分析による推定モデルを用いて正常な値の体表温度検出データを推定し、
   前記逆演算モデルは、該推定した体表温度検出データを含む体表温度データ群を用いて、ＰＬＳ回帰分析によって構築する構成としてなる請求項２に記載の装着式温度測定装置。
4. 前記ケーシングには、被測定者の体表面とは反対側に位置して外気の温度を検出する外気温度検出手段を設け、
   前記読込み手段は、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で前記体表温度検出データと一緒に該外気温度検出手段による外気温度検出データを読込み、
   前記記憶手段は、体表温度データ群に加えて、前記読込み手段による複数の外気温度検出データを外気温度データ群として記憶し、
   前記逆演算モデルは、既知の被測定者から測定した体表温度データ群と体温とに加えて外気温度データ群を用いて、ＰＬＳ回帰分析によって、得られた体表温度データ群および外気温度データ群から体温を逆演算することができる構成とし、
   前記体温予測手段は、該記憶手段に記憶した体表温度データ群および外気温度データ群と前記ＰＬＳ回帰分析によって構築した逆演算モデルとに基づいて、被測定者の体温を予測する構成としてなる請求項２または３に記載の装着式温度測定装置。
5. 前記記憶手段は、前記読込み手段による外気温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として外気温度データ群を記憶し、
   該外気温度データ群のうち異常な値の外気温度検出データに対して、主成分分析による推定モデルを用いて正常な値の外気温度検出データを推定し、
   前記逆演算モデルは、該推定した外気温度検出データを含む外気温度データ群を用いて、ＰＬＳ回帰分析によって構築する構成としてなる請求項４に記載の装着式温度測定装置。
6. 前記ケーシングには、被測定者の体表面側のうち外周縁側に位置して被測定者の体表の温度を検出する補助体表温度検出手段を設け、
   前記読込み手段は、開始時刻から終了時刻までの間に亘って予め決められた時間間隔で前記体表温度検出データと一緒に該補助体表温度検出手段による補助体表温度検出データを読込み、
   前記記憶手段は、体表温度データ群に加えて、前記読込み手段による複数の補助体表温度検出データを補助体表温度データ群として記憶し、
   前記逆演算モデルは、既知の被測定者から測定した体表温度データ群と体温とに加えて補助体表温度データ群を用いて、ＰＬＳ回帰分析によって、得られた体表温度データ群および補助体表温度データ群から体温を逆演算することができる構成とし、
   前記体温予測手段は、該記憶手段に記憶した体表温度データ群および補助体表温度データ群と前記ＰＬＳ回帰分析によって構築した逆演算モデルとに基づいて、被測定者の体温を予測する構成としてなる請求項２，３，４または５に記載の装着式温度測定装置。
7. 前記記憶手段は、前記読込み手段による補助体表温度検出データのうち予め設定された温度および温度変化の許容範囲以内の値を正常な値とし、許容範囲以外の値を異常な値として補助体表温度データ群を記憶し、
   該補助体表温度データ群のうち異常な値の補助体表温度検出データに対して、主成分分析による推定モデルを用いて正常な値の補助体表温度検出データを推定し、
   前記逆演算モデルは、該推定した補助体表温度検出データを含む補助体表温度データ群を用いて、ＰＬＳ回帰分析によって構築する構成としてなる請求項６に記載の装着式温度測定装置。
8. 前記ケーシングには、前記記憶手段に記憶したデータを外部の処理装置に転送する転送手段を設け、
   前記逆演算モデルは、外部の処理装置を用いて構築する構成としてなる請求項２，３，４，５，６または７に記載の装着式温度測定装置。
9. 前記逆演算モデルを用いて予測した被測定者の体温を複数日に亘って記憶する体温記憶手段を設け、
   該体温記憶手段に記憶した複数日の体温に対して高周波な変動成分を除去するローパスフィルタ手段を設ける構成としてなる請求項１，２，３，４，５，６，７または８に記載の装着式温度測定装置。
10. 被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定する体温推定方法であって、
    既知の被測定者について開始時刻から終了時刻までの間に亘って体表温度データ群を測定すると共にこのときの体温を測定し、
    この測定した体表温度データ群と体温とを用いて、ＰＬＳ回帰分析によって、得られた体表温度データ群から体温を逆演算することができる逆演算モデルを構築し、
    未知の被測定者について開始時刻から終了時刻までの間に亘って体表温度データ群を測定し、
    未知の被測定者について得られた体表温度データ群と前記逆演算モデルとに基づいて体温を予測する体温推定方法。

Images