(57)【要約】【課題】妊婦の血圧の測定値において、正しく測定できている測定値を誤って外れ値として検出してしまうことを低減する外れ値検出装置を提供する。【解決手段】外れ値検出装置は、所定の条件で順次に測定した妊婦の収縮期血圧値と拡張期血圧値とを含む測定値を取得する取得部と、前記妊婦の妊娠日数の情報及び／又は測定値の変動に基づいて、前記測定値の外れ値を判定する閾値を経時的に変化させる閾値決定部と、前記収縮期血圧値の時間的変化、前記拡張期血圧値の時間的変化、及び前記閾値に基づいて、前記測定値の外れ値を検出する検出部と、を有する。【選択図】図３
2020110486.tif 000002
【特許請求の範囲】
【請求項１】
所定の条件で順次に測定した妊婦の収縮期血圧値と拡張期血圧値とを含む測定値を取得する取得部と、
前記妊婦の妊娠日数の情報及び／又は測定値の変動に基づいて、前記測定値の外れ値を判定する閾値を経時的に変化させる閾値決定部と、
前記収縮期血圧値の時間的変化、前記拡張期血圧値の時間的変化、及び前記閾値に基づいて、前記測定値の外れ値を検出する検出部と、
を有する、外れ値検出装置。
【請求項２】
前記閾値決定部は、妊娠日数が分娩時妊娠日数に近いほど、前記閾値が大きくなるように前記閾値を決定する、請求項１に記載の外れ値検出装置。
【請求項３】
前記閾値決定部は、前記妊婦の前記妊娠日数の前後における前記測定値の変動から、前記妊婦の前記妊娠日数に対応する閾値を決定する、請求項１に記載の外れ値検出装置。
【請求項４】
前記閾値決定部は、複数の妊婦の前記妊娠日数の前後における前記測定値の変動から、前記妊婦の前記妊娠日数に対応する閾値を決定する、請求項１に記載の外れ値検出装置。
【請求項５】
前記閾値は、前記収縮期血圧値の閾値である第１の閾値、及び前記拡張期血圧値の閾値である第２の閾値を含む、請求項２乃至４のいずれか一項に記載の外れ値検出装置。
【請求項６】
前記収縮期血圧値の時間的変化を表す第１の平滑化曲線、及び前記拡張期血圧値の時間的変化を表す第２の平滑化曲線を算出する算出部を有し、
前記検出部は、前記収縮期血圧値と前記第１の平滑化曲線との残差が前記閾値以上であり、かつ前記拡張期血圧値と前記第２の平滑化曲線との残差が前記閾値以上である場合、前記収縮期血圧値又は前記拡張期血圧値を外れ値と判断する、請求項５に記載の外れ値検出装置。
【請求項７】
前記算出部は、前記取得部が取得した前記測定値のうち、予め定められた除外条件を満たす測定値を除外して、前記第１の平滑化曲線、及び前記第２の平滑化曲線を算出する、請求項６に記載の外れ値検出装置。
【請求項８】
前記検出部は、前記測定値が測定エラーの属性を有する場合、測定エラーの場合に用いる閾値を用いて、前記測定値の外れ値を検出する、請求項１乃至７のいずれか一項に記載の外れ値検出装置。
【請求項９】
前記検出部が検出した外れ値に、外れ値であることを示すフラグを付加して出力する出力部を有する、請求項１乃至８のいずれか一項に記載の外れ値検出装置。
【請求項１０】
所定の条件で順次に測定した妊婦の収縮期血圧値と拡張期血圧値とを含む測定値を取得する処理と、
前記妊婦の妊娠日数の情報及び／又は測定値の変動に基づいて、前記測定値の外れ値を判定する閾値を経時的に変化させる処理と、
前記収縮期血圧値の時間的変化、前記拡張期血圧値の時間的変化、及び前記閾値に基づいて、前記測定値の外れ値を検出する処理と、
をコンピュータが実行する、外れ値検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、外れ値検出装置、外れ値検出システム、及び外れ値検出方法
に関する。
【背景技術】
【０００２】
血圧に関する測定には、何かしらの測定の誤りや人的な誤りが含まれる場合がある。これらの誤りには、例えば、カフ巻きが不十分で正しく測定できなかった場合、安静にせずに測定した場合、測定値を手動でシステムに入力するときに値を間違って入力した場合等が含まれる。従って、血圧に関する測定値を用いて何らかの判断を行う場合、異常な測定値を除外、又は補正する必要がある。このような異常な測定値は、統計的に通常取り得ない外れ値となる場合が多く、これらを検出することが重要となる。
【０００３】
これに関連する技術として、生体指標の測定値から生体指標の変動パターンを作成し、新たな測定値が得られた場合に、新たな測定値と変動パターンを比較することにより、新たな測定値が異常値か否かを判断する技術が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−５７５５２
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし、例えば、妊娠期間中の妊婦の血圧測定においては、短い妊娠期間中に血圧のトレンドが大きく変化し、特に出産前後では大きく変化する。そのため、従来の手法をそのまま適用しただけでは、正しく測定できている測定値を誤って外れ値として検出してしまう場合がある。
【０００６】
本発明の一実施形態は、上記の問題点に鑑みてなされたものであって、妊婦の血圧の測定値において、正しく測定できている測定値を誤って外れ値として検出してしまうことを低減する外れ値検出装置を提供する。
【課題を解決するための手段】
【０００７】
上記の課題を解決するため、本発明の一実施形態に係る外れ値検出装置は、所定の条件で順次に測定した妊婦の収縮期血圧値と拡張期血圧値とを含む測定値を取得する取得部と、前記妊婦の妊娠日数の情報及び／又は測定値の変動に基づいて、前記測定値の外れ値を判定する閾値を経時的に変化させる閾値決定部と、前記収縮期血圧値の時間的変化、前記拡張期血圧値の時間的変化、及び前記閾値に基づいて、前記測定値の外れ値を検出する検出部と、を有する。
【発明の効果】
【０００８】
本発明の一実施形態によれば、妊婦の血圧の測定値において、正しく測定できている測定値を誤って外れ値として検出してしまうことを低減する外れ値検出装置を提供することができる。
【図面の簡単な説明】
【０００９】
【図１】一実施形態に係る外れ値検出システムのシステム構成の例を示す図である。
【図２】一実施形態に係る外れ値検出装置のハードウェア構成の例を示す図である。
【図３】一実施形態に係る外れ値検出装置の機能構成の例を示す図である。
【図４】一実施形態に係る記憶部が記憶する測定データの例を示す図である。
【図５】一実施形態に係る外れ値検出装置の処理の概要を示すフローチャートである。
【図６】一実施形態に係る外れ値検出装置の処理の具体的な一例を示すフローチャートである。
【図７】一実施形態に係る外れ値を検出する処理の具体的な一例を示すフローチャートである。
【図８】一実施形態に係る外れ値検出装置が出力する出力データの例を示す図である。
【発明を実施するための形態】
【００１０】
以下、図面を参照して本発明の実施の形態を説明する。なお、以下で説明する実施の形態は一例であり、本発明が適用される実施の形態は、以下の実施の形態に限られない。
【００１１】
＜システム構成＞
図１は、一実施形態に係る外れ値検出システムのシステム構成の例を示す図である。外れ値検出システム１００は、一例として、血圧を測定する血圧計１０２と、血圧計１０２が測定した血圧の測定値を順次に記憶した測定データから外れ値を検出する外れ値検出装置１０１とを含む。
【００１２】
別の一例として、外れ値検出装置１０１は、通信ネットワーク１０６等を介して通信可能なサーバ装置１０５等から、血圧の測定値を順次に記憶した測定データを取得し、取得した測定データから外れ値を検出するものであっても良い。
【００１３】
血圧計１０２等により測定する血圧の測定値には、何かしらの測定誤りや人的誤りが含まれる場合がある。この測定値の誤りには、例えば、カフ巻きが不十分で正しく測定できなかった場合、安静にせずに測定した場合、測定値を手動でシステムに入力するときに値を間違って入力した場合、別の人が誤って計測した場合等が含まれ得る。
【００１４】
これらの測定値を用いて何らかの判断を行う場合、異常な測定値を除外又は補正することが望ましい。異常な測定値は、統計的に通常は取り得ない外れ値となる場合が多く、外れ値検出装置１０１は、血圧の測定値からこのような外れ値を検出して、検出結果を、例えば、表示装置１０３、記憶装置１０４、又はサーバ装置１０５等に出力する。
【００１５】
一般的に、血圧の測定値のように繰り返し測定する測定値は一定ではなく、例えば、正規分布等のように所定の分布に従うことが多い。このように、測定値が所定の分布に従う場合、統計的に外れ値を検出する方法が用いられる。例えば、測定値が正規分布に従う場合、標準偏差をσとして、平均値等から３σ以上離れた測定値を外れ値として検出することが行われている。
【００１６】
しかし、例えば、妊娠期間中の妊婦の血圧測定においては、比較的短い妊娠期間中に血圧のトレンド（傾向）が大きく変化し、特に出産前後では大きく変化する。このような場合、従来の技術をそのまま適用してしまうと、正しく測定できている測定値を誤って外れ値として検出してしまう場合がある。
【００１７】
そこで、本実施形態に係る外れ値検出装置１０１は、妊婦の妊娠日数の情報及び／又は測定値の変動に基づいて、妊婦の血圧の測定値における外れ値を判定する閾値を経時的に変化させる機能を有している。例えば、外れ値検出装置１０１は、妊婦本人の経時的な血圧値もしくは他の妊婦群の統計的かつ経時的な血圧値から、当該妊婦の各妊娠日数時点での血圧変化の許容できる閾値を算出し、その閾値により外れ値を判定する基準を緩める（例えば、閾値を大きくする）。また、外れ値検出装置１０１は、出産日前後の所定の期間（例えば、前後２０〜２５日間程度）において、段階的に閾値を変化させることで、外れ値を判定する基準を緩める（例えば、閾値を大きくする）。
【００１８】
また、外れ値検出装置１０１は、収縮期血圧値（ＳＢＰ：Systolic Blood Pressure）の時間的変化と、拡張期血圧値（ＤＢＰ：Diastolic Blood Pressure）の時間的変化と、前述した閾値と、に基づいて測定値の外れ値を検出する。
【００１９】
これにより、本実施形態に係る外れ値検出装置１０１は、例えば、出産日の前後等、血圧値のトレンドが大きく変化する場合でも、正しく測定できている測定値を誤って外れ値として検出してしまうことを低減することができる。また、血圧値のトレンドには個人差があるが、本実施形態に係る外れ値検出装置１０１は、妊婦の血圧値のトレンドに基づいて外れ値を検出するので、より正しく外れ値を検出することができるようになる。
【００２０】
好ましくは、外れ値検出装置１０１は、出産日前後の所定の期間において、収縮期血圧値（以下、ＳＢＰと呼ぶ）と、拡張期血圧値（以下、ＤＢＰと呼ぶ）との両方が、妊婦の血圧値のトレンドから乖離しているときに、測定値を外れ値と判断する。これは、例えば、カフ巻きが不十分で正しく測定できなかった場合、安静にせずに測定した場合等、測定自体が適切に行われていない場合には、ＳＢＰ及びＤＢＰの両方の測定値が、血圧値のトレンドから乖離すると考えられるためである。
【００２１】
以上、本実施形態によれば、妊婦の血圧の測定値において、正しく測定できている測定値を誤って外れ値として検出してしまうことを低減する外れ値検出装置１０１を提供することができる。
【００２２】
なお、図１に示す外れ値検出システム１００のシステム構成は一例である。例えば、外れ値検出装置１０１の機能のうち、少なくとも一部は、血圧計１０２が有していても良いし、サーバ装置１０５が有していても良い。また、表示装置１０３、記憶装置１０４等は、外れ値検出装置１０１が有していても良い。
【００２３】
＜ハードウェア構成＞
（外れ値検出装置のハードウェア構成）
図２は、一実施形態に係る外れ値検出装置１０１のハードウェア構成の例を示す図である。外れ値検出装置１０１は、物理的には、プロセッサ２０１、メモリ２０２、ストレージ２０３、入力装置２０４、出力装置２０５、通信装置２０６、及びバス２０７等を含むコンピュータとして構成されても良い。なお、以下の説明では、「装置」という文言は、回路、デバイス、ユニット等に読み替えることができる。
【００２４】
プロセッサ２０１は、例えば、オペレーティングシステムを動作させてコンピュータ全体を制御する。プロセッサ２０１は、周辺装置とのインタフェース、制御装置、演算装置、レジスタ等を含む中央処理装置（ＣＰＵ：Central Processing Unit）で構成されても良い。
【００２５】
また、プロセッサ２０１は、プログラム（プログラムコード）、ソフトウェアモジュールやデータを、ストレージ２０３及び／又は通信装置２０６からメモリ２０２に読み出し、これらに従って各種の処理を実行する。プログラムとしては、外れ値検出装置１０１の動作の少なくとも一部をコンピュータに実行させるプログラムが用いられる。外れ値検出装置１０１において実行される各種処理は、１つのプロセッサ２０１で実行されても良いし、２以上のプロセッサ２０１により同時又は逐次に実行されても良い。プロセッサ２０１は、１以上のチップで実装されても良い。なお、プログラムは、電気通信回線を介してネットワークから送信されても良い。
【００２６】
メモリ２０２は、コンピュータ読み取り可能な記憶媒体であり、例えば、ＲＯＭ（Read Only Memory）、ＥＰＲＯＭ（Erasable Programmable ＲＯＭ）、ＥＥＰＲＯＭ（Electrically Erasable Programmable ＲＯＭ）、ＲＡＭ（Random Access Memory）等の少なくとも１つで構成されても良い。メモリ２０２は、レジスタ、キャッシュ、メインメモリ（主記憶装置）等と呼ばれても良い。メモリ２０２は、本発明の一実施の形態に係る外れ値検出方法を実施するために実行可能なプログラム（プログラムコード）、ソフトウェアモジュール等を保存することができる。
【００２７】
ストレージ２０３は、コンピュータ読み取り可能な記憶媒体であり、例えば、ＣＤ−ＲＯＭ（Compact Disc ＲＯＭ）等の光ディスク、ハードディスクドライブ、フレキシブルディスク、光磁気ディスク（例えば、コンパクトディスク、デジタル多用途ディスク、Ｂｌｕ−ｒａｙ（登録商標）ディスク）、スマートカード、フラッシュメモリ（例えば、カード、スティック、キードライブ）、フロッピー（登録商標）ディスク、磁気ストリップ、ＳＳＤ（Solid State Drive）、その他不揮発性メモリ等の少なくとも１つで構成されても良い。ストレージ２０３は、補助記憶装置と呼ばれても良い。上述の記憶媒体は、例えば、メモリ２０２及び／又はストレージ２０３を含むデータベース、サーバその他の適切な媒体であっても良い。
【００２８】
入力装置２０４は、外部からの入力を受け付ける入力デバイス（例えば、キーボード、マウス、マイクロフォン、スイッチ、ボタン、センサ等）である。出力装置２０５は、外部への出力を実施する出力デバイス（例えば、ディスプレイ、スピーカ、ＬＥＤランプ等）である。なお、入力装置２０４及び出力装置２０５は、一体となった構成（例えば、タッチパネルディスプレイ）であっても良い。
【００２９】
通信装置２０６は、有線及び／又は無線ネットワークを介してコンピュータ間の通信を行うためのハードウェア（送受信デバイス）であり、例えばネットワークデバイス、ネットワークコントローラ、ネットワークカード、通信モジュール等ともいう。また、通信装置２０６は、近距離無線通信により、外部装置と直接通信を行う機能を有していても良い。
【００３０】
上記のプロセッサ２０１やメモリ２０２等の各装置は、情報を通信するためのバス２０７で接続される。バス２０７は、単一のバスで構成されても良いし、装置間で異なるバスで構成されても良い。また、バスは有線で接続されても良いし、無線で接続されても良い。
【００３１】
＜機能構成＞
（外れ値検出装置の機能構成）
図３は、一実施形態に係る外れ値検出装置の機能構成の例を示す図である。外れ値検出装置１０１は、例えば、図２のプロセッサ２０１で所定のプログラムを実行することにより、測定データ取得部３０１、記憶部３０２、情報取得部３０３、データ抽出部３０４、算出部３０５、閾値決定部３０６、外れ値検出部３０７、及び出力部３０８等を実現している。なお、上記の各機能構成のうち、少なくとも一部は、ハードウェアによって実現されるものであっても良い。
【００３２】
測定データ取得部（取得部）３０１は、例えば、図２のプロセッサ２０１で実行されるプログラムによって実現され、所定の条件で順次に測定した妊婦のＳＢＰ（収縮期血圧）とＤＢＰ（拡張期血圧）とを含む測定値を取得する。例えば、測定データ取得部３０１は、図１の血圧計１０２で測定された妊婦の血圧の測定値を取得し、記憶部３０２に順次に記憶する。
【００３３】
或いは、測定データ取得部３０１は、図１のサーバ装置１０５、記憶装置１０４等から、妊婦の血圧の測定値を順次に記憶した測定データを取得し、記憶部３０２に記憶するものであっても良い。
【００３４】
記憶部３０２は、例えば、図２のプロセッサ２０１で実行されるプログラム、及びストレージ２０３、メモリ２０２等によって実現され、測定データ取得部３０１が取得した利用者の血圧の測定値等を記憶する。
【００３５】
図４は、一実施形態に係る記憶部が記憶する測定データの例を示す図である。図４（ａ）は、出産前に測定した妊婦の血圧の測定値を順次に記憶した、産前の測定データ４０１の例を示している。図４（ａ）の例では、産前の測定データ４０１には、項目として、「ＩＤ」、「妊娠日数ｔ」、「日付」、「時刻」、「ＳＢＰ」、「ＤＢＰ」、「入力方法」、「測定エラー」等の情報が記憶されている。
【００３６】
「ＩＤ」の情報は、妊婦を識別する識別情報である。なお、ＩＤは、妊婦の氏名等の情報であっても良い。好ましくは、記憶部３０２には、図４（ａ）に示すように、妊婦毎に各情報が記憶されている。
【００３７】
「妊娠日数ｔ」の情報は、基本的には、妊娠前の最終月経開始日からの経過日数を示す情報であり、２８０日目を出産予定日とする。「日付」の情報は、測定値を測定した日付を示す情報である。「時刻」の情報は、測定値を測定した時刻を示す情報である。「ＳＢＰ」の情報は、収縮期血圧の測定値である。「ＤＢＰ」の情報は、拡張期血圧の測定値である。
【００３８】
「入力方法」の情報は、測定値の入力方法が、自動であるか、手動であるかを示す情報である。「測定エラー」の情報は、入力方法が自動である場合に、血圧計１０２等で検出されたエラーの情報が記憶される。
【００３９】
図４（ｂ）は、出産後に測定した妊婦の血圧の測定値を順次に測定した、産後の測定データ４０２の例を示している。産後の測定データ４０２には、産前の測定データ４０１と同様に、項目として、「ＩＤ」、「産後日数ｔ」、「日付」、「時刻」、「ＳＢＰ」、「ＤＢＰ」、「入力方法」、「測定エラー」等の情報を含む。ただし、産後の測定データ４０２では、「妊娠日数ｔ」に代わり、出産日からの経過日数を示す「産後日数ｔ」になっている。なお、他の情報については、産前の測定データ４０１と同様である。
【００４０】
ここで、図３に戻り、外れ値検出装置１０１の機能構成の説明を続ける。
【００４１】
情報取得部３０３は、例えば、図２のプロセッサ２０１で実行されるプログラムによって実現され、妊婦の妊娠日数の情報（例えば、妊娠日数ｔ、出産予定日、出産日等）を取得する。
【００４２】
一例として、情報取得部３０３は、図１のサーバ装置１０５から出産予定日、出産日等の妊娠日数に関する情報を取得する。または、情報取得部３０３は、図２の出力装置２０５等に、例えば、出産予定日、出産日等の妊娠日数に関する情報の入力を受付する入力画面を表示し、入力装置２０４等に入力された妊娠日数に関する情報を取得する。
【００４３】
また、情報取得部３０３は、例えば、図４（ａ）に示すような産前の測定データ４０１に記憶された測定日である「日付」の情報と、取得した妊娠日数に関する情報（例えば、出産予定日の情報）とに基づいて、「妊娠日数ｔ」の情報を算出し、産前の測定データ４０１に記憶する。
【００４４】
さらに、情報取得部３０３は、例えば、図４（ｂ）に示すような産後の測定データ４０２に記憶された「日付」の情報と、取得した妊娠日数に関する情報（例えば、出産日の情報）とに基づいて、「産後日数ｔ」の情報を算出し、産後の測定データ４０２に記憶する。
【００４５】
別の一例として、情報取得部３０３は、取得した妊娠日数に関する情報（例えば、妊娠日数ｔ、出産予定日、出産日等）を、例えば、そのまま記憶部３０２に記憶しておくものであっても良い。この場合、妊娠日数ｔの情報は、必要なときに、「日付」の情報と、記憶部３０２に記憶した妊娠日数に関する情報（例えば、出産予定日の情報）とを用いて算出する。この場合、産前の測定データ４０１には、「妊娠日数ｔ」の情報が含まれていなくても良く、また産後の測定データ４０２には、「産後日数ｔ」の情報が含まれていなくても良い。
【００４６】
データ抽出部３０４は、例えば、図２のプロセッサ２０１で実行されるプログラムによって実現される。データ抽出部３０４は、記憶部３０２に記憶された、産前の測定データ４０１、産後の測定データ４０２から、算出部３０５が妊婦の血圧値のトレンドを算出する際に用いる測定データ（ＳＢＰ、及びＤＢＰ）を抽出する。
【００４７】
算出部３０５は、例えば、図２のプロセッサ２０１で実行されるプログラムによって実現され、対象となる、ある妊婦の血圧値のトレンドを算出する。例えば、算出部３０５は、ＳＢＰの測定値の時間的変化（トレンド）を表す第１の平滑化曲線、及びＤＢＰの測定値の時間的変化を表す第２の平滑化曲線を算出する。
【００４８】
例えば、算出部３０５は、ＬＯＷＥＳＳ（Locally Weighted Scatterplot Smoother）等の平滑化曲線アルゴリズムを用いて、ユーザ（妊婦）の血圧値の傾向を表す平滑化曲線を算出する。なお、平滑化の方式は、ＬＯＷＥＳＳに限られず、ＬＯＥＳＳ（Locally estimated Scatterplot Smoothing）、平均移動、基底関数等を用いた方式を用いても良い。
【００４９】
なお、本実施形態では、図４（ａ）、（ｂ）で説明したように、出産前と出産後とで、妊娠日数ｔと産後日数ｔの定義が異なる。そのため、算出部３０５は、妊婦の血圧値のトレンドを表す平滑化曲線を算出する場合、出産前と出産後とに分けて、別々に処理を行う。
【００５０】
また、血圧値は、例えば、朝、晩等、測定する時間帯によって血圧の測定値が変化する。そのため、算出部３０５は、例えば、朝、晩等の時間帯に分けて、別々に妊婦の血圧値のトレンドを表す平滑化曲線を算出することが望ましい。
【００５１】
例えば、データ抽出部３０４は、算出部３０５が、出産前の朝の平滑化曲線を算出する場合、図４（ａ）に示すような産前の測定データ４０１から、朝の時間帯に該当する測定データを抽出する。また、データ抽出部３０４は、算出部３０５が、出産後の晩の平滑化曲線を算出する場合、図４（ｂ）に示すような産後の測定データ４０２から、晩の時間帯に該当する測定データを抽出する。
【００５２】
好ましくは、データ抽出部３０４は、抽出する測定データの中から、明らかな外れ値を除外する。例えば、データ抽出部３０４は、図４（ａ）、（ｂ）に示すような測定データにおいて、入力方法が「手動」となっており、かつ「ＳＢＰ」、及び「ＤＢＰ」の値がデフォルト値である場合、第１種の外れ値として、当該測定値を、平滑化曲線を算出するための測定データから除外する。
【００５３】
デフォルト値は、例えば、血圧の測定値を入力するアプリケーションにおいて、測定値の入力する前に、利用者の入力の手間を最小限にするために一般利用者の平均的な値などとして予め設定されている値であり、血圧値を手動で入力した場合に入力間違いとして高頻度に出現するため、当該測定値を第１の外れ値としてフラグを付加して管理する。例えば、ある血圧入力アプリケーションにおいてＳＢＰのデフォルト値が１２０ｍｍＨｇに設定され、ＤＢＰのデフォルト値が８０ｍｍＨｇに設定されている場合がある。
【００５４】
閾値決定部３０６は、例えば、図２のプロセッサ２０１で実行されるプログラムによって実現され、外れ値検出部３０７が、外れ値を判定するための閾値を経時的に変化させる。例えば、閾値決定部３０６は、当該妊婦の経時的な血圧値もしくは他の妊婦群の統計的かつ経時的な血圧値から、当該妊婦の各妊娠日数時点での血圧変化の許容できる閾値を変動させる。また、閾値決定部３０６は、出産日前後の所定の期間において、妊婦の出産予定日、又は出産日が近い程、閾値が大きくなるように（外れ値を判定する基準が緩くなるように）閾値を決定する。なお、閾値決定部３０６による閾値の決定方法の具体的な例については後述（〔００７８〕）する。
【００５５】
外れ値検出部（検出部）３０７は、例えば、図２のプロセッサ２０１で実行されるプログラムによって実現され、ＳＢＰの時間的変化、ＤＢＰの時間的変化、及び閾値決定部３０６が決定した閾値に基づいて、血圧の測定値から外れ値を検出する。
【００５６】
好ましくは、外れ値検出部３０７は、出産日前後の所定の期間において、ＳＢＰとＤＢＰとの両方の測定値が、算出部３０５が算出した血圧値のトレンドから乖離している場合、当該測定値を外れ値と判断する。なお、外れ値検出部３０７による外れ値の検出方法の具体的な例については後述（〔００７８〕）する。
【００５７】
出力部３０８は、図２のプロセッサ２０１で実行されるプログラムによって実現される。出力部３０８は、外れ値検出部３０７が検出した外れ値に、外れ値であることを示すフラグを付加した血圧値の測定結果を、例えば、図１の表示装置１０３、記憶装置１０４、又はサーバ装置１０５等に出力する。
【００５８】
なお、図３に示す外れ値検出装置１０１の機能構成は一例である。例えば、図３の情報取得部３０３は、図１の血圧計１０２等が備えていても良いし、図３の記憶部３０２は、図１の記憶装置１０４、又はサーバ装置１０５が備えていても良い。このように、図３に示す外れ値検出装置１０１の各機能構成は、外れ値検出システム１００に含まれる各装置のいずれかが備えていれば良い。
【００５９】
＜処理の流れ＞
続いて、本実施形態に係る外れ値検出方法の処理の流れについて説明する。
【００６０】
（処理の概要）
図５は、一実施形態に係る外れ値検出装置の処理の概要を示すフローチャートである。この処理は、外れ値検出装置１０１が実行する処理の概要を示している。なお、図５に示す処理の開始時点において、例えば、図４（ａ）に示すような産前の測定データ４０１、及び図４（ｂ）に示すような産後の測定データ４０２が、外れ値検出装置１０１の記憶部３０２に記憶されているものとする。
【００６１】
ステップＳ５０１において、外れ値検出装置１０１のデータ抽出部３０４は、所定の条件で順次に測定した血圧の測定値（ＳＢＰ、及びＤＢＰ）を、記憶部３０２から取得する。
【００６２】
外れ値検出装置１０１は、前述したように、産前の測定データ４０１と、産後の測定データ４０２のように、条件が異なる測定データを別々に処理する。また、外れ値検出装置１０１は、同じ産前の測定データ４０１であっても、例えば、朝の時間帯に測定した測定データと、晩の時間帯に測定して測定データのように、条件が異なる測定データを別々に処理する。ここでは、一例として、データ抽出部３０４が、図４（ａ）に示すような産前の測定データ４０１から、朝の時間帯における血圧の測定値（ＳＢＰ、及びＤＢＰ）を取得するものとして、以下の説明を行う。
【００６３】
好ましくは、データ抽出部３０４は、取得したＳＢＰ、及びＤＢＰの測定値のうち、明らかな外れ値を除外する。例えば、データ抽出部３０４は、図４（ａ）に示すような産前の測定データ４０１において、入力方法が「手動」となっており、かつ「ＳＢＰ」、及び「ＤＢＰ」の値がデフォルト値である測定値（ＳＢＰ、及びＤＢＰ）を除外する。
【００６４】
ステップＳ５０２において、外れ値検出装置１０１の算出部３０５は、データ抽出部３０４が取得したＳＢＰ、及びＤＢＰを用いて、ＳＢＰの時間的変化を示す第１の平滑化曲線ｆ＿ＳＢＰ、及びＤＢＰの時間的変化を示す第２の平滑化曲線ｆ＿ＤＢＰを算出する。例えば、算出部３０５は、ＬＯＷＥＳＳ等の平滑化曲線アルゴリズムを用いて、第１の平滑化曲線ｆ＿ＳＢＰ、及び第２の平滑化曲線ｆ＿ＤＢＰを算出する。
【００６５】
ステップＳ５０３において、外れ値検出装置１０１の外れ値検出部３０７は、算出部３０５が算出した第１の平滑化曲線ｆ＿ＳＢＰ、第２の平滑化曲線ｆ＿ＤＢＰ、及び閾値決定部３０６が決定する閾値を用いて、ＳＢＰ、及びＤＢＰの外れ値を検出する。なお、閾値決定部３０６が閾値を決定する処理、及び、外れ値検出部３０７が、外れ値を検出する処理については後述（〔００７８〕）する。
【００６６】
ステップＳ５０４において、外れ値検出部３０７によって検出された外れ値に、外れ値であることを示すフラグを付加して、例えば、図１の表示装置１０３、記憶装置１０４、又はサーバ装置１０５等に測定結果を出力する。
【００６７】
（外れ値検出装置の処理の具体例）
図６は、一実施形態に係る外れ値検出装置の処理の具体的な一例を示すフローチャートである。この処理は、外れ値検出処理の具体的な一例を示している。
【００６８】
なお、図６に示す処理の開始時点において、妊婦Ａの産前の測定データ４０１、及び産後の測定データ４０２が、記憶部３０２に記憶されているものとする。また、産前の測定データ４０１、及び産後の測定データ４０２には、１日２回、朝晩に測定した血圧の測定値が記憶されているものとする。
【００６９】
ステップＳ６０１において、外れ値検出装置１０１の測定データ取得部３０１は、記憶部３０２に記憶されている妊婦Ａの産前の測定データ４０１、及び産後の測定データ４０２を取得し、取得した測定データを朝と晩の血圧値に分類する。なお、朝晩の分類は、例えば、図４（ａ）、（ｂ）に示す測定データの「時刻」の情報に基づいて行う。
【００７０】
ステップＳ６０２〜Ｓ６０８において、外れ値検出装置１０１は、朝、晩それぞれのデータ群毎に、各妊娠日数ｔの測定値について、ステップＳ６０３〜Ｓ６０７の処理をループして実行する。なお、ここでは、データ群が朝である場合の例について説明するが、外れ値検出装置１０１は、データ群が晩である場合についても同様の処理を行う。
【００７１】
ステップＳ６０３において、外れ値検出装置１０１のデータ抽出部３０４は、ＳＢＰ、ＤＢＰのトレンド（平滑化曲線）を算出する前に、影響のある外れ値を予め除外するために、第１の外れ値検出を行う。具体的な一例として、朝のデータ群において、入力方法が「手動」となっており、かつＳＢＰ及びＤＢＰの値がデフォルト値である測定値を第１種外れ値としてフラグを付加する。なお、入力方法が「手動」となっており、かつＳＢＰ及びＤＢＰの値がデフォルト値である測定値は、予め定められた除外条件の一例である。除外条件は、他の条件であっても良い。
【００７２】
ここで、第１の外れ値と判断された測定値は、ＳＢＰ、ＤＢＰのトレンド算出用データから除外されるが、外れ値検出装置１０１が最終的に出力する測定結果には、第１種外れ値のフラグを付加して出力される。
【００７３】
ステップＳ６０４において、外れ値検出装置１０１の算出部３０５は、産前の測定データ４０１、産後の測定データ４０２のそれぞれについて、ＳＢＰの平滑化曲線ｆ＿ＳＢＰ、及びＤＢＰの平滑化曲線ｆ＿ＤＢＰを算出する。例えば、算出部３０５は、ＬＯＷＥＳＳ等の平滑化曲線のアルゴリズムを用いて、妊婦ＡのＳＢＰの平滑化曲線ｆ＿ＳＢＰ、及びＤＢＰの平滑化曲線ｆ＿ＤＢＰを算出する。なお、算出部３０５は、ＬＯＷＥＳＳに限られず、ＬＯＥＳＳ、移動平均、基底関数を用いたフィッティング等により、ｆ＿ＳＢＰ、ｆ＿ＤＢＰを算出しても良い。
【００７４】
ステップＳ６０５では、外れ値検出装置１０１の外れ値検出部３０７は、産前の測定データ４０１のＳＢＰ、ＤＢＰについて、ＳＢＰの平滑化曲線ｆ＿ＳＢＰ、ＤＢＰの平滑化曲線ｆ＿ＤＢＰ、及び閾値決定部３０６が決定する閾値を用いて、外れ値を検出する。例えば、外れ値検出部３０７は、測定値が、平滑化曲線より閾値以上離れていれば、第２種の外れ値としてフラグを付加する。なお、閾値決定部３０６及び外れ値検出部３０７の具体的な処理の例については後述（〔００７８〕）する。
【００７５】
ステップＳ６０６において、外れ値検出部３０７は、産後の測定データ４０１のＳＢＰ、ＤＢＰについて、ＳＢＰの平滑化曲線ｆ＿ＳＢＰ、ＤＢＰの平滑化曲線ｆ＿ＤＢＰ、及び閾値決定部３０６が決定する閾値を用いて、外れ値を検出する。
【００７６】
ステップＳ６０７において、外れ値検出装置１０１の出力部３０８は、第１種外れ値のフラグを付加した測定値、及び第２種外れ値のフラグを付加した測定値を含む測定結果を、例えば、表示装置１０３、記憶装置１０４、又はサーバ装置１０５等に出力する。
【００７７】
外れ値検出装置１０１は、上記の処理を、各測定値に対して実行することにより、各測定値が外れ値であるか否かを判断し、外れ値と判断した測定値に、外れ値であることを示すフラグを付加することができる。
【００７８】
（外れ値を検出する処理の具体例）
図７は、一実施形態に係る外れ値を検出する処理の具体的な一例を示すフローチャートである。この処理は、図６のステップＳ６０５において、外れ値検出装置１０１が実行する第２種外れ値の検出処理の具体的な一例を示している。
【００７９】
ステップＳ７０１において、外れ値検出装置１０１の外れ値検出部３０７は、産前の各時点において、ＳＢＰとｆ＿ＳＢＰとの残差ε＿ＳＢＰ、ＳＢＰの標準偏差σ＿ＳＢＰ、ＤＢＰとｆ＿ＤＢＰとの残差ε＿ＤＢＰ、及びＤＢＰの標準偏差σ＿ＤＢＰを算出する。なお、標準偏差の算出は、例えば、産前の測定データ４０１の全体から計算しても良いし、一定の区間に区切ってその中で計算するものであっても良い。
【００８０】
ステップＳ７０２〜Ｓ７１５において、外れ値検出装置１０１は、各妊娠日数ｔでの測定値の外れ値を判定するために、各妊娠日数ｔの測定値について、ループして処理を行う。例えば、外れ値検出装置１０１は、産前の測定値であれば、妊娠日数ｔが０〜２９０のうち、実際のデータが存在する範囲を選択し、各妊娠日数ｔについて、ステップＳ７０３〜Ｓ７１４の処理を実行する。また、外れ値検出装置１０１は、産後の測定値であれば、産後０日以降のデータのうち、実際のデータが存在する範囲を選択し、各産後日数ｔについて、ステップＳ７０３〜Ｓ７１４の処理を実行する。
【００８１】
ステップＳ７０３において、外れ値検出装置１０１の外れ値検出部３０７は、
ａｂｓ（ｔ−ｔ＿ｄ）＞ｔ＿１、かつ
ａｂｓ（ε＿ＳＢＰ（ｔ））≧ｂ＿１
であるか否かを判断する。
【００８２】
ここで、ａｂｓ（ｔ−ｔ＿ｄ）は、妊娠日数ｔと、分娩時妊娠日数ｔ＿ｄとの差の大きさであり、ｔ＿１は、予め定められた出産前後期間１（例えば、２４日等）である。これは、出産日の前後ｔ＿１より遠い期間（第２の期間という）において、予め定められた閾値ｂ＿１よりも大きい残差ε＿ＳＢＰであった場合に、それを外れ値と判断する。後段のステップＳ９０５やステップＳ９０７では標準偏差σ＿ＳＢＰなどから〔００８６〕で説明する閾値１，閾値２を利用して外れ値判断を行うが、血圧計測点数が少ない場合などでσ＿ＳＢＰが大きくなると閾値１，閾値２も大きくなり、本来外れ値とするべきものを見逃す危険性がある。その場合にも、本来外れ値とする測定値を外れ値と判定するために、ステップＳ９０３にて固定値のｂ＿１を用いて外れ値をまず判定する。なお、ステップＳ９０３の処理は補助的であるため、ステップＳ９０３を省略しても良い。または、ｔ＿１＝０として、全ての妊娠日数ｔをステップＳ９０３の対象としても良い。
【００８３】
また、ε＿ＳＢＰ（ｔ）は、妊娠日数ｔにおけるＳＢＰとｆ＿ＳＢＰとの残差ε＿ＳＢＰであり、ｂ＿１は、妊娠日数ｔが、出産前後でない場合の血圧外れ値の基準（閾値）である。ｔ＿１やｂ＿１などは、従来の経験に則った範囲で設定可能で、例えば、発明した閾値決定方法Ａのモデル（図６の一連の流れ）に対して血圧データとその外れ値の正例負例から学習して、もっとも外れ値の識別性能が高くなるように計算することができ、一例として、ｔ＿１＝２４［日］、ｂ＿１＝４０［ｍｍＨｇ］を適用することができるが、これに限られない。
【００８４】
ステップＳ７０３の判断結果が「ＹＥＳ」である場合、外れ値検出部３０７は、処理をステップＳ７０４に移行させる。一方、ステップＳ７０３の判断結果が「ＮＯ」である場合、外れ値検出部３０７は、処理をステップＳ７０５に移行させる。
【００８５】
ステップＳ７０４に移行すると、外れ値検出部３０７は、妊娠日数ｔのＳＢＰの測定値と共に、ＳＢＰの第２種外れ値情報として外れ値理由１のフラグを付加して、処理をステップＳ９０９に移行させる。なお、外れ値理由１のフラグが付加された測定値は、妊娠日数ｔが、分娩時妊娠日数ｔ＿ｄからｔ＿１日より離れた第２の期間であり、残差ε＿ＳＢＰがｂ＿１以上離れた、明らかな外れ値であることを示している。
【００８６】
ステップＳ７０３からステップＳ７０５に移行すると、外れ値検出部３０７は、
測定値の属性が測定エラーでない、かつ
ａｂｓ（ε＿ＳＢＰ（ｔ））＞閾値１、かつ
ａｂｓ（ε＿ＤＢＰ（ｔ））＞閾値２
であるか否かを判断する。
【００８７】
ここで、閾値１は、ＳＢＰに関する閾値（第１の閾値）であり、閾値２は、ＤＢＰに関する閾値（第２の閾値）である。
【００８８】
閾値決定部３０６及び外れ値検出部３０７は、妊娠の経過に応じて、外れ値と判定するための閾値１〜２を妊娠日数ごとに変化させることで、急激な変化に対応した外れ値判定を行う。この閾値決定の方法として、以下の３つをあげる。
閾値決定方法Ａ：妊娠日数ｔが分娩時妊娠日数ｔ＿ｄに近いほど、閾値が大きくなるように閾値を決定する方法
閾値決定方法Ｂ：当該妊婦の妊娠日数ｔの血圧値の前後の変動から妊娠日数ｔに対応する閾値を決定する方法
閾値決定方法Ｃ：他の妊婦を含めて妊娠日数ｔの血圧値の前後の変動から妊娠日数ｔに対応する閾値を決定する方法
閾値決定部３０６は、これらの何れの方法を用いても良い。
【００８９】
まず、閾値決定方法Ａについて説明する。これは、妊娠日数ｔが分娩時妊娠日数ｔ＿ｄに近いほど、閾値が大きくなるように閾値１、閾値２を決定する方法である。
【００９０】
閾値決定方法Ａの閾値１は、例えば、σ＿ＳＢＰの定数ｓ倍、又は固定値ｂ＿２の大きい方に、分娩時妊娠日数ｔ＿ｄに近づくほど大きくなる値を足したものである。例えば、閾値決定方法Ａの閾値１は、次の式で表すことができる。
閾値１＝ｍａｘ（ｓ＊σ＿ＳＢＰ，ｂ＿２）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
【００９１】
ここで、ｍａｘは関数であり、ｍａｘ（Ａ，Ｂ）はＡ，Ｂのうち大きい方を表し、ｍａｘ（Ａ，Ｂ，Ｃ）はＡ，Ｂ，Ｃのうち最大のものを表す、プログラミング言語などで一般的な表記法である。ａｂｓは絶対値を表す関数であり、例えばａｂｓ（−１０）は１０を表す。ｓは、ε＿ＳＢＰが、σ＿ＳＢＰの何倍乖離している場合に外れ値とするかを決める定数であり、例えば、ｓ＝３を適用することができる。ＳＢＰの各妊娠日数ｔでの測定値が正規分布をしていた場合、ｓ＝３では測定値の分布の外側約０．３％を外れ値と見なすような基準となる。ｓには０以上の任意の値が取れるが、主に１〜４の値を使うのが望ましい。
【００９２】
ｂ＿２は、ステップＳ９０３のｂ＿１と同様の目的で用いられ、血圧計測点数が少ない場合などでσ＿ＳＢＰが極端に大きい場合にも、明らかに異常な外れ値を判定するための基準である。ｂ＿２やｂ＿４は、発明した閾値決定方法Ａのモデル（図６の一連の流れ）に対して血圧データとその外れ値の正例負例から学習して、もっとも外れ値の識別性能が高くなるように計算することができる。事前に入手した妊婦数３００名のデータに基づいて計算した場合、例えば、ｂ＿２＝２９．５［ｍｍＨｇ］を適用することができる。ｂ＿４は、妊娠日数ｔが、分娩時妊娠日数ｔ＿ｄに近づくにつれて閾値１を増加させる量（基準緩和度）であり、同様の計算により、例えば、ｂ＿４＝３［ｍｍＨｇ］を適用することができる。
【００９３】
また、「ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４」は、分娩時妊娠日数ｔ＿ｄに近づくほど大きくなる値の一例である。なお、上記のｓ、ｂ＿２、ｂ＿４の値は一例であり、他の値であっても良い。
【００９４】
同様に、閾値決定方法Ａの閾値２は次の式で表すことができる。
閾値２＝ｍａｘ（ｓ＊σ＿ＤＢＰ，ｂ＿２）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
閾値２はＤＢＰに対する基準であるが、ステップＳ９０５では、ＳＢＰの外れ値判定にＳＢＰの変動だけを見るのではなく、ＤＢＰの変動も見ることが特長である。これは、測定に異常があった場合などはＳＢＰとＤＢＰの両方が外れ値となりやすいという特徴を反映したものである。つまり、妊娠日数ｔにおいてＤＢＰが正常の範囲である場合で残差ε＿ＳＢＰが大きくなった場合は、ＳＢＰは外れ値ではなく、正しく測定できているものの極端に高値又は低値であったと見なす。
【００９５】
次に、閾値決定方法Ｂについて説明する。これは、当該妊婦の妊娠日数ｔの血圧値の前後ｕ日の変動から妊娠日数ｔに対応する閾値１〜２を決定する方法である。閾値決定方法Ｂの閾値１は、次の式で表すことができる。
閾値１＝ｍａｘ（ｓ＊σ'＿ＳＢＰ，ｂ＿２）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
【００９６】
σ＿ＳＢＰが当該妊婦の対象期間（例えば産前）全てのＳＢＰ測定値から計算された、静的な（妊娠日数ｔによって変動しない）標準偏差であったのに対し、σ'＿ＳＢＰは、当該妊婦の妊娠日数ｔの前後ｕ日の期間内におけるＳＢＰ測定値から計算された動的な（妊娠日数ｔによって変動する）標準偏差である。これにより、当該妊婦本人の血圧が大きく変動する場合にはσ'＿ＳＢＰが大きくなり、変動に合わせて閾値を変化させることができる。例えばｕ＝７日とすると、妊娠日数ｔの前後７日の計１５日におけるＳＢＰ測定値の標準偏差からσ'＿ＳＢＰが計算される。この１５日のＳＢＰ測定値が急激に変動していた場合、σ'＿ＳＢＰが大きくなり閾値１を大きく変化させ、逆に１５日のＳＢＰ測定値の変動が少ない場合、σ'＿ＳＢＰが小さくなり閾値１を小さく抑えることで、ＳＢＰの変動に追従することができる。さらに、閾値決定方法Ｂではσ'＿ＳＢＰによって出産日前後の変動に追従できることから、閾値１の第２項であるｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４を省略しても良い。つまりｂ＿４＝０としても良い。
【００９７】
同様に、閾値決定方法Ｂの閾値２は次の式で表すことができる。
閾値２＝ｍａｘ（ｓ＊σ'＿ＤＢＰ，ｂ＿２）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
ここでσ'＿ＤＢＰは、当該妊婦の妊娠日数ｔの前後ｕ日の期間内におけるＤＢＰ測定値から計算された動的な（妊娠日数ｔによって変動する）標準偏差である。また、閾値決定方法Ｂの閾値１と同様、第２項であるｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４を省略しても良い。つまりｂ＿４＝０としても良い。
【００９８】
次に、閾値決定方法Ｃについて説明する。これは、他の妊婦を含めて妊娠日数ｔの血圧値の前後ｕ日の変動から妊娠日数ｔに対応する閾値１〜２を決定する方法である。閾値決定方法Ｃの閾値１は、次の式で表すことができる。
閾値１＝ｍａｘ（ｓ'＊σ'＿ＳＢＰ＋ｓ''＊σ''＿ＳＢＰ，ｂ＿２）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
【００９９】
閾値決定方法Ｂと同様に、σ'＿ＳＢＰが、当該妊婦の妊娠日数ｔの前後ｕ日の期間内におけるＳＢＰ測定値から計算された標準偏差である。これに加え、当該妊婦だけでない多数の妊婦での妊娠日数ｔの前後ｕ日の期間内におけるＳＢＰ測定値から標準偏差を計算したものがσ''＿ＳＢＰである。閾値決定方法Ｃでは、これら２つの標準偏差をｓ'およびｓ''で重み付けして足し合わせた値と、ｂ＿２の大きいほうが閾値１の第１項として採用される。望ましくはｓ≒ｓ'＋ｓ''であり、例えばｓ'＝２，Ｓ''＝１などの値が取り得る。σ'＿ＳＢＰとσ''＿ＳＢＰを組み合わせることによって、当該妊婦の妊娠日数ｔ近辺の変動だけでなく、他の妊婦の妊娠日数ｔの前後ｕ日の期間内の変動を加味することができ、例えば当該妊婦の測定時点数が少ないなどの理由でσ'＿ＳＢＰが過小に算出された場合でも、他の妊婦のσ''＿ＳＢＰでその影響を軽減することができる。
さらに、閾値決定方法Ｃではσ'＿ＳＢＰやσ''＿ＳＢＰによって出産日前後の変動に追従できることから、閾値１の第２項であるｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４を省略しても良い。つまりｂ＿４＝０としても良い。
【０１００】
同様に、閾値決定方法Ｃの閾値２は次の式で表すことができる。
閾値２＝ｍａｘ（ｓ'＊σ'＿ＤＢＰ＋ｓ''＊σ''＿ＤＢＰ，ｂ＿２）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
ここでσ'＿ＤＢＰは、当該妊婦の妊娠日数ｔの前後ｕ日の期間内におけるＤＢＰ測定値から計算された動的な標準偏差であり、当該妊婦だけでない多数の妊婦での妊娠日数ｔの前後ｕ日の期間内におけるＤＢＰ測定値から標準偏差を計算したものがσ''＿ＤＢＰである。また、閾値決定方法Ｃの閾値１と同様、第２項であるｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４を省略しても良い。つまりｂ＿４＝０としても良い。
【０１０１】
ステップＳ７０５の判断結果が「ＹＥＳ」である場合、外れ値検出部３０７は、処理をステップＳ７０６に移行させる。一方、ステップＳ７０５の判断結果が「ＮＯ」である場合、外れ値検出部３０７は、処理をステップＳ７０７に移行させる。
【０１０２】
ステップＳ７０６に移行すると、外れ値検出部３０７は、妊娠日数ｔのＳＢＰの測定値と共に、ＳＢＰの第２種外れ値情報として外れ値理由２のフラグを付加して、処理をステップＳ９０９に移行させる。なお、外れ値理由２のフラグが付加された測定値は、測定値の属性が測定エラーでない場合の外れ値であることを示している。
【０１０３】
ステップＳ７０５からステップＳ７０７に移行すると、外れ値検出部３０７は、
測定値の属性が測定エラーであり、かつ
ａｂｓ（ε＿ＳＢＰ（ｔ））＞閾値３
であるか否かを判断する。
【０１０４】
閾値３は、測定値が測定エラーである場合の閾値である。外れ値検出部３０７は、妊娠の経過に応じて、外れ値と判定するための閾値３を妊娠日数ごとに変化させることで、急激な変化に対応した外れ値判定を行う。この閾値決定の方法として、ステップＳ９０５と同様に閾値決定方法Ａ、Ｂ、Ｃの３つをあげる。これらの何れの方法を用いても良い。
【０１０５】
まず、閾値決定方法Ａについて説明する。これは、妊娠日数ｔが分娩時妊娠日数ｔ＿ｄに近いほど、閾値が大きくなるように閾値３を決定する方法である。
閾値３は、例えば、次の式で表すことができる。
閾値３＝ｍａｘ（ｓ＊σ＿ＳＢＰ，ｂ＿３）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
【０１０６】
ここで、ｂ＿３は、測定値の属性が測定エラーである場合の血圧外れ値の基準であり、例えば、ｂ＿３＝２７．０［ｍｍＨｇ］を適用することができる。これは、発明した閾値決定方法Ａのモデル（図６の一連の流れ）に対して血圧データとその外れ値の正例負例から学習して、もっとも外れ値の識別性能が高くなるように計算した場合である。なお、このｂ＿３の値は一例であり、他の値であっても良い。
【０１０７】
次に、閾値決定方法Ｂについて説明する。これは、当該妊婦の妊娠日数ｔの血圧値の前後ｕ日の変動から妊娠日数ｔに対応する閾値３を決定する方法である。閾値決定方法Ｂの閾値３は、次の式で表すことができる。
閾値３＝ｍａｘ（ｓ＊σ'＿ＳＢＰ，ｂ＿３）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
【０１０８】
閾値決定方法Ｂではσ'＿ＳＢＰによって出産日前後の変動に追従できることから、閾値３の第２項であるｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４を省略しても良い。つまりｂ＿４＝０としても良い。
【０１０９】
次に、閾値決定方法Ｃについて説明する。これは、他の妊婦を含めて妊娠日数ｔの血圧値の前後ｕ日の変動から妊娠日数ｔに対応する閾値３を決定する方法である。閾値決定方法Ｃの閾値３は、次の式で表すことができる。
閾値３＝ｍａｘ（ｓ'＊σ'＿ＳＢＰ＋ｓ''＊σ''＿ＳＢＰ，ｂ＿３）
＋ｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４
【０１１０】
閾値決定方法Ｃではσ'＿ＳＢＰやσ''＿ＳＢＰによって出産日前後の変動に追従できることから、閾値３の第２項であるｍａｘ（ｔ＿２−ａｂｓ（ｔ−ｔ＿ｄ），０）＊ｂ＿４を省略しても良い。つまりｂ＿４＝０としても良い。
【０１１１】
ステップＳ７０７の判断結果が「ＹＥＳ」である場合、外れ値検出部３０７は、処理をステップＳ７０８に移行させる。一方、ステップＳ７０７の判断結果が「ＮＯ」である場合、外れ値検出部３０７は、処理をステップＳ７０９に移行させる。
【０１１２】
ステップＳ７０８に移行すると、外れ値検出部３０７は、妊娠日数ｔのＳＢＰの測定値と共に、ＳＢＰの第２種外れ値情報として外れ値理由３のフラグを付加して、処理をステップＳ９０９に移行させる。なお、外れ値理由３のフラグが付加された測定値は、測定値の属性が測定エラーである場合の外れ値であることを示している。
【０１１３】
ステップＳ７０９〜Ｓ７１４において、外れ値検出装置１０１は、ステップＳ７０３〜Ｓ７０８でＳＢＰに対して行った処理と同様の処理を、ＤＢＰに対して行う。このとき、ステップＳ７０３〜Ｓ７０８におけるＳＢＰをＤＢＰと読み替え、ＤＢＰをＳＢＰと読み替える。つまり、閾値１および３はＤＢＰに対するものであり、閾値２はＳＢＰに対するものとなる。
【０１１４】
上記の処理により、産前の測定データ４０１のＳＢＰ、及びＤＢＰの外れ値に、外れ値であることを示すフラグを付加することができる。
【０１１５】
また、外れ値検出装置１０１は、図９に示す処理と同様にして、産後の測定データ４０２のＳＢＰ、及びＤＢＰの外れ値に、外れ値であることを示すフラグを付加することができる。
【０１１６】
（出力データの例）
図１０は、一実施形態に係る外れ値検出装置が出力する出力データの例を示す図である。図１０（ａ）は、産前の出力データ１００１の例を示している。産前の出力データ１００１は、一例として、図４（ａ）に示す産前の測定データ４０１の項目に加えて、「ＳＢＰの第１種外れ値」、「ＤＢＰの第１種外れ値」、「ｆ＿ＳＢＰ」、「ε＿ＳＢＰ」、「σ＿ＳＢＰ」、「ｆ＿ＤＢＰ」、「ε＿ＤＢＰ」、「σ＿ＤＢＰ」、「ＳＢＰの第２種外れ値」、「ＤＢＰの第２種外れ値」等の情報が追加されている。これらは一例であって、閾値決定方法によっては、それに使用した量としてσ'＿ＳＢＰ、σ'＿ＤＢＰ、σ''＿ＳＢＰ、σ''＿ＤＢＰ、閾値１、閾値２、閾値３などを追記しても良い。
【０１１７】
「第１種外れ値」の情報は、例えば、図８のステップＳ８０３において、外れ値検出装置１０１が、入力方法が「手動」となっており、かつＳＢＰ及びＤＢＰの値がデフォルト値である測定値に付加する第１種外れ値フラグの情報である。図１０（ａ）の例では、第１外れ値と判断された測定値にフラグ（例えば、「○」）が付加されている。
【０１１８】
「ｆ＿ＳＢＰ」、「ｆ＿ＤＢＰ」の情報は、例えば、図８のステップＳ８０４で算出された、ＳＤＰの平滑化曲線ｆ＿ＳＢＰ、及びＤＢＰの平滑化曲線ｆ＿ＤＢＰの妊娠日数ｔに対応する値を示している。
【０１１９】
「ε＿ＳＢＰ」、「ε＿ＤＢＰ」の情報は、例えば、図９のステップＳ９０１で算出された、ＳＢＰとｆ＿ＳＢＰとの残差ε＿ＳＢＰ、及びＤＢＰとｆ＿ＤＢＰとの残差ε＿ＤＢＰの値を示している。
【０１２０】
「σ＿ＳＢＰ」、「σ＿ＤＢＰ」の情報は、例えば、図９のステップＳ９０１で算出された、ＳＢＰの標準偏差σ＿ＳＢＰ、及びＳＢＰの標準偏差σ＿ＤＢＰの値を示している。
【０１２１】
「第２種外れ値」の情報は、図９のステップＳ９０４、Ｓ９０６、又はＳ９０８で付加された、外れ値理由１〜３の情報を示している。
【０１２２】
図１０（ｂ）は、産後の出力データ１００２の例を示している。産後の出力データ１００２は、一例として、図４（ｂ）に示す産後の測定データ４０２の項目に加えて、「第１種外れ値」、「ｆ＿ＳＢＰ」、「ε＿ＳＢＰ」、「σ＿ＳＢＰ」、「ｆ＿ＤＢＰ」、「ε＿ＤＢＰ」、「σ＿ＤＢＰ」、「第２種外れ値」等の情報が追加されている。なお、産後の出力データ１００２に追加されている各情報は、産前の出力データ１００１に追加されている各情報と同様なので、ここでは説明を省略する。
【０１２３】
このように、外れ値検出装置１０１は、例えば、図４（ａ）、（ｂ）に示すような血圧の測定値から外れ値を検出し、検出された外れ値に、外れ値であることを示すフラグを付加して、図１０（ａ）、（ｂ）に示すような出力データ（測定結果）を出力する。
【０１２４】
＜実施形態の効果＞
本実施形態に係る外れ値検出装置１０１は、血圧の測定値の平滑化曲線を元に外れ値判定の閾値を経時的に変動させ、かつＳＢＰとＤＢＰの両方の変化を組み合わせて判定する。これにより、外れ値検出装置１０１は、例えば、妊婦の血圧測定のように、短い妊娠期間中に血圧のトレンドが大きく変化する場合であっても、正しく測定できている測定値を誤って外れ値として検出してしまうことを抑制することができる。また、ＳＢＰとＤＢＰの両方の変化を利用した検出法は、妊婦以外にも一般的な血圧測定での外れ値検出でも効果がある。
【０１２５】
＜補足＞
なお、図３の機能構成図は、機能単位のブロックを示している。これらの機能ブロックは、ハードウェア及び／又はソフトウェアの任意の組み合わせによって実現される。また、各機能ブロックの実現手段は特に限定されない。すなわち、各機能ブロックは、物理的及び／又は論理的に結合した１つの装置により実現されても良いし、物理的及び／又は論理的に分離した２つ以上の装置を直接的及び／又は間接的に（例えば、有線及び／又は無線）で接続し、これら複数の装置により実現されても良い。
【０１２６】
また、図２に示す外れ値検出装置１０１のハードウェア構成は、図に示した各装置を１つ又は複数含むように構成されても良いし、一部の装置を含まずに構成されても良い。また、外れ値検出装置１０１は、マイクロプロセッサ、デジタル信号プロセッサ（ＤＳＰ：Digital Signal Processor）、ＡＳＩＣ（Application Specific Integrated Circuit）、ＰＬＤ（Programmable Logic Device）、ＦＰＧＡ（Field Programmable Gate Array）等のハードウェアを含んで構成されても良く、当該ハードウェアにより、各機能ブロックの一部又は全てが実現されても良い。例えば、プロセッサ２０１は、これらのハードウェアの少なくとも１つで実装されても良い。
【０１２７】
本明細書で説明した各態様／実施形態の処理手順、シーケンス、フローチャート等は、矛盾のない限り、順序を入れ替えても良い。例えば、本明細書で説明した方法については、例示的な順序で様々なステップの要素を提示しており、提示した特定の順序に限定されない。
【０１２８】
入出力された情報等は特定の場所（例えば、メモリ）に保存されても良いし、管理テーブルで管理しても良い。入出力される情報等は、上書き、更新、又は追記され得る。出力された情報等は削除されても良い。入力された情報等は他の装置へ送信されても良い。
【０１２９】
判定は、１ビットで表される値（０か１か）によって行われても良いし、真偽値（Boolean：true又はfalse）によって行われても良いし、数値の比較（例えば、所定の値との比較）によって行われても良い。
【０１３０】
本明細書で説明した各態様／実施形態は単独で用いても良いし、組み合わせて用いても良いし、実行に伴って切り替えて用いても良い。また、所定の情報の通知（例えば、「Ｘであること」の通知）は、明示的に行うものに限られず、暗黙的（例えば、当該所定の情報の通知を行わない）ことによって行われても良い。
【０１３１】
ソフトウェアは、ソフトウェア、ファームウェア、ミドルウェア、マイクロコード、ハードウェア記述言語と呼ばれるか、他の名称で呼ばれるかを問わず、命令、命令セット、コード、コードセグメント、プログラムコード、プログラム、サブプログラム、ソフトウェアモジュール、アプリケーション、ソフトウェアアプリケーション、ソフトウェアパッケージ、ルーチン、サブルーチン、オブジェクト、実行可能ファイル、実行スレッド、手順、機能などを意味するよう広く解釈されるべきである。
【０１３２】
また、ソフトウェア、命令などは、伝送媒体を介して送受信されても良い。例えば、ソフトウェアが、同軸ケーブル、光ファイバケーブル、ツイストペア及びデジタル加入者回線（ＤＳＬ）などの有線技術及び／又は赤外線、無線及びマイクロ波などの無線技術を使用してウェブサイト、サーバ、又は他のリモートソースから送信される場合、これらの有線技術及び／又は無線技術は、伝送媒体の定義内に含まれる。
【０１３３】
本明細書で説明した情報、信号等は、様々な異なる技術のいずれかを使用して表されても良い。例えば、上記の説明全体に渡って言及され得るデータ、命令、コマンド、情報、信号、ビット、シンボル、チップなどは、電圧、電流、電磁波、磁界若しくは磁性粒子、光場若しくは光子、又はこれらの任意の組み合わせによって表されても良い。
【０１３４】
なお、本明細書で説明した用語及び／又は本明細書の理解に必要な用語については、同一の又は類似する意味を有する用語と置き換えても良い。
【０１３５】
また、本明細書で説明した情報、パラメータなどは、絶対値で表されても良いし、所定の値からの相対値で表されても良いし、対応する別の情報で表されてもよい。
【０１３６】
本明細書で使用する「に基づいて」という記載は、別段に明記されていない限り、「のみに基づいて」を意味しない。言い換えれば、「に基づいて」という記載は、「のみに基づいて」と「に少なくとも基づいて」の両方を意味する。
【０１３７】
「含む(ｉｎｃｌｕｄｉｎｇ)」、「含んでいる（ｃｏｍｐｒｉｓｉｎｇ）」、及びそれらの変形が、本明細書、或いは特許請求の範囲で使用されている限り、これら用語は、用語「備える」と同様に、包括的であることが意図される。さらに、本明細書、或いは特許請求の範囲において使用されている用語「または（or）」は、排他的論理和ではないことが意図される。
【０１３８】
本開示の全体において、例えば、英語でのa, an, 及びtheのように、翻訳により冠詞が追加された場合、これらの冠詞は、文脈から明らかにそうではないことが示されていなければ、複数のものを含むものとする。
【０１３９】
以上、本発明について詳細に説明したが、当業者にとっては、本発明が本明細書中に説明した実施形態に限定されるものではないということは明らかである。本発明は、特許請求の範囲の記載により定まる本発明の趣旨及び範囲を逸脱することなく修正及び変更態様として実施することができる。従って、本明細書の記載は、例示説明を目的とするものであり、本発明に対して何ら制限的な意味を有するものではない。
【符号の説明】
【０１４０】
１００      外れ値検出システム
１０１      外れ値検出装置
３０１      測定データ取得部（取得部）
３０５      算出部
３０６      閾値決定部
３０７      外れ値検出部（検出部）
３０８      出力部
【図１】
2020110486.tif 000003
【図２】
2020110486.tif 000004
【図３】
2020110486.tif 000005
【図４】
2020110486.tif 000006
【図５】
2020110486.tif 000007
【図６】
2020110486.tif 000008
【図７】
2020110486.tif 000009
【図８】
2020110486.tif 000010