JP4933227B2 - 装着式温度測定装置および身体状態推定方法 - Google Patents

Description

本発明は、例えば就寝中に測定した被測定者の体表温度に基づいて被測定者が高温期と低温期のうちいずれの状態にあるかを推定する装着式温度測定装置および身体状態推定方法に関する。

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

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

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

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

本発明は上述した従来技術の問題に鑑みなされたもので、本発明の目的は、複数日に亘る体表温度の検出データに基づいて被測定者が高温期と低温期のうちいずれの状態にあるかを推定することができる装着式温度測定装置および身体状態推定方法を提供することにある。

上述した課題を解決するために請求項１の発明は、被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定し、該体温の推定値に基づいて被測定者が低温期と高温期とのうちいずれの状態にあるのか判断する装着式温度測定装置であって、前記推定した被測定者の体温の推定値を複数日に亘って記憶する体温記憶手段を設け、被測定者からの月経開始連絡に基づいて月経開始日を決定し、該月経開始日から数えて第４日目から第１１日目までの体温の推定値を平均した平均値を低温期の基準温度とし、該低温期の基準温度よりも高い温度として、低温期から高温期への上昇を捉えることが可能な所定の温度差を該低温期の基準温度に加えた温度を高温期の基準温度とし、前記体温記憶手段に記憶した体温の推定値と前記高温期の基準温度とを比較し、前記月経開始日以降で前記体温の推定値が前記高温期の基準温度よりも低いときには、前記被測定者は低温期にあるものと判断し、前記体温の推定値が前記高温期の基準温度よりも上昇したときには、前記被測定者は低温期から高温期に移行したものと判断する構成としたことにある。

請求項２の発明では、基礎体温の１周期を所定の割り当て日数をもった８個のステージに区分し、第１〜第４ステージは低温期に対応し、残余の第５〜第８ステージは高温期に対応するものとし、前記月経開始日を前記第１ステージの開始日として、各ステージの割り当て日数が経過すると、次のステージに順次移行するものとし、前記第４ステージでは、前記被測定者が低温期から高温期に移行したと判断したときに、第５ステージに移行し、それ以外のときには第４ステージに留まる構成としている。

請求項３の発明では、被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定し、該体温の推定値に基づいて被測定者が低温期と高温期とのうちいずれの状態にあるのか判断する身体状態推定方法であって、前記推定した被測定者の体温の推定値を複数日に亘って記憶し、被測定者からの月経開始連絡に基づいて月経開始日を決定し、該月経開始日から数えて第４日目から第１１日目までの体温の推定値を平均した平均値を低温期の基準温度とし、該低温期の基準温度よりも高い温度として、低温期から高温期への上昇を捉えることが可能な所定の温度差を該低温期の基準温度に加えた温度を高温期の基準温度とし、前記体温記憶手段に記憶した体温の推定値と前記高温期の基準温度とを比較し、前記月経開始日以降で前記体温の推定値が前記高温期の基準温度よりも低いときには、前記被測定者は低温期にあるものと判断し、前記体温の推定値が前記高温期の基準温度よりも上昇したときには、前記被測定者は低温期から高温期に移行したものと判断する構成としている。
請求項４の発明では、被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定し、該体温の推定値に基づいて被測定者の状態を推定する身体状態推定方法であって、前記推定した被測定者の体温の推定値を複数日に亘って記憶し、月経開始日から数えて第４日目から第１１日目までの体温の推定値を平均した平均値を低温期の基準温度とし、該低温期の基準温度よりも高い温度として、低温期から高温期への上昇を捉えることが可能な所定の温度差を該低温期の基準温度に加えた温度を高温期の基準温度とする構成としている。

請求項１，３の発明によれば、月経開始日から数えて第４日目から第１１日目までの体温の推定値を平均した平均値を低温期の基準温度とし、該低温期の基準温度よりも高い温度として、低温期から高温期への上昇を捉えることが可能な所定の温度差を該低温期の基準温度に加えた温度を高温期の基準温度とした。このため、体温記憶手段に記憶した体温の推定値と高温期の基準温度とを比較し、月経開始日以降で体温の推定値が高温期の基準温度よりも低いときには、被測定者は低温期にあるものと判断し、体温の推定値が高温期の基準温度よりも上昇したときには、被測定者は低温期から高温期に移行したものと判断することができる。

請求項２の発明によれば、基礎体温の１周期を所定の割り当て日数をもった８個のステージに区分し、第１〜第４ステージは低温期に対応し、残余の第５〜第８ステージは高温期に対応するものとした。そして、第４ステージでは、被測定者が低温期から高温期に移行したと判断したときに、第５ステージに移行し、それ以外のときには第４ステージに留まる。これにより、被測定者が高温期と低温期のうちいずれの状態にあるか容易に把握することができる。
請求項４の発明によれば、月経開始日から数えて第４日目から第１１日目までの体温の推定値を平均した平均値を低温期の基準温度とし、所定の温度差を該低温期の基準温度に加えた温度を高温期の基準温度とした。このため、体温の推定値と高温期の基準温度とを比較し、月経開始日以降で体温の推定値が高温期の基準温度よりも低いときには、被測定者は低温期にあるものと判断し、体温の推定値が高温期の基準温度よりも上昇したときには、被測定者は低温期から高温期に移行したものと判断することができる。

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

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

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

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

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

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

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

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

ここで、記憶部８には、コントロールユニット７を作動させるプログラムと該プログラムで使用する開始時刻ｔ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 (4)

1. 被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定し、該体温の推定値に基づいて被測定者が低温期と高温期とのうちいずれの状態にあるのか判断する装着式温度測定装置であって、
   前記推定した被測定者の体温の推定値を複数日に亘って記憶する体温記憶手段を設け、
   被測定者からの月経開始連絡に基づいて月経開始日を決定し、
   該月経開始日から数えて第４日目から第１１日目までの体温の推定値を平均した平均値を低温期の基準温度とし、
   該低温期の基準温度よりも高い温度として、低温期から高温期への上昇を捉えることが可能な所定の温度差を該低温期の基準温度に加えた温度を高温期の基準温度とし、
   前記体温記憶手段に記憶した体温の推定値と前記高温期の基準温度とを比較し、前記月経開始日以降で前記体温の推定値が前記高温期の基準温度よりも低いときには、前記被測定者は低温期にあるものと判断し、前記体温の推定値が前記高温期の基準温度よりも上昇したときには、前記被測定者は低温期から高温期に移行したものと判断する構成としてなる装着式温度測定装置。
2. 基礎体温の１周期を所定の割り当て日数をもった８個のステージに区分し、第１〜第４ステージは低温期に対応し、残余の第５〜第８ステージは高温期に対応するものとし、
   前記月経開始日を前記第１ステージの開始日として、各ステージの割り当て日数が経過すると、次のステージに順次移行するものとし、
   前記第４ステージでは、前記被測定者が低温期から高温期に移行したと判断したときに、第５ステージに移行し、それ以外のときには第４ステージに留まる構成としてなる請求項１に記載の装着式温度測定装置。
3. 被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定し、該体温の推定値に基づいて被測定者が低温期と高温期とのうちいずれの状態にあるのか判断する身体状態推定方法であって、
   前記推定した被測定者の体温の推定値を複数日に亘って記憶し、
   被測定者からの月経開始連絡に基づいて月経開始日を決定し、
   該月経開始日から数えて第４日目から第１１日目までの体温の推定値を平均した平均値を低温期の基準温度とし、
   該低温期の基準温度よりも高い温度として、低温期から高温期への上昇を捉えることが可能な所定の温度差を該低温期の基準温度に加えた温度を高温期の基準温度とし、
   前記体温記憶手段に記憶した体温の推定値と前記高温期の基準温度とを比較し、前記月経開始日以降で前記体温の推定値が前記高温期の基準温度よりも低いときには、前記被測定者は低温期にあるものと判断し、前記体温の推定値が前記高温期の基準温度よりも上昇したときには、前記被測定者は低温期から高温期に移行したものと判断する構成としてなる身体状態推定方法。
4. 被測定者の体表温度を開始時刻から終了時刻までの間に亘って測定した時系列の体表温度データ群を用いて被測定者の体温を推定し、該体温の推定値に基づいて被測定者の状態を推定する身体状態推定方法であって、
   前記推定した被測定者の体温の推定値を複数日に亘って記憶し、
   月経開始日から数えて第４日目から第１１日目までの体温の推定値を平均した平均値を低温期の基準温度とし、
   該低温期の基準温度よりも高い温度として、低温期から高温期への上昇を捉えることが可能な所定の温度差を該低温期の基準温度に加えた温度を高温期の基準温度とする構成としてなる身体状態推定方法。

Images