(57)【要約】【課題】  過去の花粉飛散量など−の環境に関する時系列データおよび花粉症状などの症状の時系列データと、将来の環境に関する予測データとから、将来の症状を予測し、対策支援を可能とする疾病対策支援方法を提供することにある。また、通信ネットワークに接続され、かつ、前記の疾病対策支援方法を有するサーバおよび端末から構成される疾病対策支援システムを提供する。【解決手段】  コンピュータプログラム１Ａとして疾病対策支援方法を有し、かつ、個人に関する疾病症状時系列データ１Ｃと、環境に関する時系列データ１Ｂとを有し、症状に応じた対策情報を提供するサーバ１と、通信回線１０を介してサーバ１から症状に応じた対策情報を受信する端末２とを備える。【選択図】  図６
2005063218.tif 000002
【特許請求の範囲】
【請求項１】
将来の疾病症状を、過去の疾病症状時系列データと、過去の環境に関する時系列データと、将来の環境予測データとの関数とみなし、
前記データが得られるごとに、前記関数に含まれる係数を推定し、疾病症状の予測式を算出し、
最新の前記データを前記疾病症状の予測式に入力し、
前記疾病症状の予測式の値により将来の疾病症状を予測し、症状に応じた対策情報を提供する、
ことを特徴とする疾病対策支援方法。
【請求項２】
予測した疾病症状と実際の疾病症状との誤差を用いて、前記疾病症状の予測式の各係数を更新することを特徴とする請求項１に記載の疾病対策支援方法。
【請求項３】
過去の疾病症状時系列データと、過去の環境に関する時系列データと、将来の環境予測データとからなるベクトルを構成し、
前記ベクトルが得られるごとに、実際の疾病症状に応じて分類しカテゴリ付けしたものをデータベースとし、
最新の前記ベクトルを前記データベースと比較することにより、この最新のベクトルを分類し、
得られたカテゴリに対応する疾病症状により、疾病症状を予測し、症状に応じた対策情報を提供する、
ことを特徴とする疾病対策支援方法。
【請求項４】
将来の疾病症状を、過去の疾病症状時系列データと、過去の環境に関する時系列データと、将来の環境予測データとの関数とみなし、
予測する疾病症状の特性がそれぞれ異なり、かつ、指標付けされた前記関数を複数用意し、
ユーザの疾病症状データおよび環境データを前記複数の関数に入力し、
実際の疾病症状に最も適合する値を出力する関数を選択し、
その関数に付けられた指標にユーザの症状特性を分類する、
ことを特徴とする疾病対策支援方法。
【請求項５】
過去の疾病症状時系列データと、過去の環境に関する時系列データと、将来の環境予測データとからなるベクトルを構成し、
複数の前記ベクトルを実際の疾病症状に応じて分類し、カテゴリ付けしたものをデータベースとし、
予測する疾病症状の特性がそれぞれ異なり、かつ、指標付けされた前記データベースを複数用意し、
ユーザの疾病症データおよび環境データからなるベクトルを、前記複数のデータベースにて分類し、
指標付けされた複数のデータベースの中で正しい疾病症状予測を行ったデータベースを選択し、そのデータベースに付けられた指標にユーザの症状特性を分類する、
ことを特徴とする疾病対策支援方法。
【請求項６】
疾病症状と環境状態をセンサによって自動的に記録することを特徴とする請求項１〜５のいずれか１項に記載の疾病対策支援方法。
【請求項７】
疾病症状レベルをユーザが設定可能であることを特徴とする請求項１〜６のいずれか１項に記載の疾病対策支援方法。
【請求項８】
対象とする疾病は花粉症であることを特徴とする請求項１〜７のいずれか１項に記載の疾病対策支援方法。
【請求項９】
請求項１から７のいずれかの疾病対策支援方法を用いた疾病対策支援システムであって、
コンピュータプログラム（１Ａ）として、請求項１から７のいずれかの疾病対策支援方法を有し、かつ、個人に関する疾病症状時系列データ（１Ｃ）と環境に関する時系列データ（１Ｂ）とを有し、症状に応じた対策情報を提供するサーバ（１）と、
通信ネットワークを介して前記サーバ（１）から症状に応じた対策情報を受信する端末（２）と、
を備えることを特徴とする疾病対策支援システム。
【請求項１０】
請求項１から７のいずれかの疾病対策支援方法を用いた疾病対策支援システムであって、
環境に関する時系列データ（２１Ａ）を有するサーバ（２１）と、
個人に関する疾病症状時系列データ（２２Ｂ）を有し、通信ネットワーク（３０）を介して前記サーバ（２１）から環境に関する時系列データ（２１Ａ）を受信し、かつ、コンピュータプログラム（２２Ａ）として請求項１から７のいずれかの疾病対策支援方法を有し、症状に応じた対策情報を提供する端末（２２）と、
を備えることを特徴とする疾病対策支援システム。
【請求項１１】
対象とする疾病は花粉症であることを特徴とする請求項９または１０に記載の疾病対策支援システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、環境に関する時系列データと、個人に関する疾病症状時系列データとから将来の疾病症状を予測し、症状に応じた対策情報を提供する疾病対策支援方法に関する。また、通信ネットワークに接続され、かつ、疾病対策支援方法をプログラムとして有するサーバおよび端末から構成される疾病対策支援システムに関する。
【背景技術】
【０００２】
従来、医療機関における診療を補う目的で、例えば在宅健康管理システムがある。この在宅健康管理システムは、個人が日常生活の中で手軽に健康状態を診断したり、あるいは、疾病の兆候を検出して、疾病の予防をするためのものである。
【０００３】
このように、在宅にて環境状態に起因する疾病の予防や、健康状態の診断を行うものとして、日々の健康データに基づいた個人健康履歴のデータベースと、気象の変化とより、多変量解析を用いて、疾病発生予測をするものがある（例えば、特許文献１）。また、環境センサによる環境データと、健康センサによる健康データとを用いて、個人の疾病発生や症状変化を予測するものがある（例えば、特許文献２）。さらに、遺伝情報や健康情報などの個人情報を端末に保持し、環境情報および環境と健康の相関情報をサーバから受信することにより、個人の健康状態を診断するものがある（例えば、特許文献３）。
【特許文献１】特開２００１−０６７４０３（特願平１１−２３９０１６号）  在宅健康管理システム
【特許文献２】特開２００２−３１１１５８（特願２００１−１２１６９０）  医学気象予報配信システム、医学気象予報配信方法、医学気象予報配信プログラム
【特許文献３】特開２００３−０６７４９９（特願２００１−２５１９１６）  健康診断方法および健康診断システム、健康診断システム用サーバおよび端末、ならびに健康診断システム用プログラム
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかし、前記の従来技術には、以下のような問題があった。すなわち、いずれも環境データと個人の健康データとの相関情報から、個人の健康状態を診断したり、予測したりする。しかし、その相関情報といったもの自体、どのように作成するかが問題となる。環境に起因する疾病、例えば花粉症に関しても、花粉の飛散量と花粉症状との相関に関してはまだ不明な点が多く、花粉の飛散量が花粉症状に直接結びつくわけではなく、また個人差も大きい。現在の花粉飛散量の多寡だけでなく、過去の飛散量の履歴に症状が大きく依存する例も多い。
【０００５】
さらに、花粉症を抑える薬は症状の程度に応じた複数種類のものがあり、どれを利用するかは患者の判断に任されている。あるいは、マスクをした方が良い、花粉の飛散量が多いところには行かないほうが良いなどの様々な判断を患者は自分で下す必要がある。このような環境要因に影響を受ける花粉症やアトピーなどのアレルギー疾患や、日焼けなどの症状に応じた対策や判断を支援するようなシステムが必要とされている。
【０００６】
本発明は、前記の課題を解決し、過去の花粉飛散量などの環境に関する時系列データおよび花粉症状などの症状の時系列データと、将来の環境に関する予測データとから、将来の症状を予測し、対策支援を可能とする疾病対策支援方法を提供することにある。また、通信ネットワークに接続され、かつ、前記の疾病対策支援方法を有するサーバおよび端末から構成される疾病対策支援システムを提供することにある。
【課題を解決するための手段】
【０００７】
前記課題を解決するために、請求項１の発明は、将来の疾病症状を、過去の疾病症状時系列データと、過去の環境に関する時系列データと、将来の環境予測データとの関数とみなし、前記データが得られるごとに、前記関数に含まれる係数を推定し、疾病症状の予測式を算出し、最新の前記データを前記疾病症状の予測式に入力し、前記疾病症状の予測式の値により将来の疾病症状を予測し、症状に応じた対策情報を提供することを特徴とする疾病対策支援方法である。
請求項２の発明は、請求項１に記載の疾病対策支援方法において、予測した疾病症状と実際の疾病症状との誤差を用いて、前記疾病症状の予測式の各係数を更新することを特徴とする。
請求項３の発明は、過去の疾病症状時系列データと、過去の環境に関する時系列データと、将来の環境予測データとからなるベクトルを構成し、前記ベクトルが得られるごとに、実際の疾病症状に応じて分類しカテゴリ付けしたものをデータベースとし、最新の前記ベクトルを前記データベースと比較することにより、この最新のベクトルを分類し、得られたカテゴリに対応する疾病症状により、疾病症状を予測し、症状に応じた対策情報を提供することを特徴とする疾病対策支援方法である。
請求項４の発明は、将来の疾病症状を、過去の疾病症状時系列データと、過去の環境に関する時系列データと、将来の環境予測データとの関数とみなし、予測する疾病症状の特性がそれぞれ異なり、かつ、指標付けされた前記関数を複数用意し、ユーザの疾病症状データおよび環境データを前記複数の関数に入力し、実際の疾病症状に最も適合する値を出力する関数を選択し、その関数に付けられた指標にユーザの症状特性を分類することを特徴とする疾病対策支援方法である。
請求項５の発明は、過去の疾病症状時系列データと、過去の環境に関する時系列データと、将来の環境予測データとからなるベクトルを構成し、複数の前記ベクトルを実際の疾病症状に応じて分類し、カテゴリ付けしたものをデータベースとし、予測する疾病症状の特性がそれぞれ異なり、かつ、指標付けされた前記データベースを複数用意し、ユーザの疾病症データおよび環境データからなるベクトルを、前記複数のデータベースにて分類し、指標付けされた複数のデータベースの中で正しい疾病症状予測を行ったデータベースを選択し、そのデータベースに付けられた指標にユーザの症状特性を分類することを特徴とする疾病対策支援方法である。
請求項６の発明は、請求項１〜５のいずれか１項に記載の疾病対策支援方法において、疾病症状と環境状態をセンサによって自動的に記録することを特徴とする。
請求項７の発明は、請求項１〜６のいずれか１項に記載の疾病対策支援方法において、疾病症状レベルをユーザが設定可能であることを特徴とする。
請求項８の発明は、請求項１〜７のいずれか１項に記載の疾病対策支援方法において、対象とする疾病は花粉症であることを特徴とする。
請求項９の発明は、請求項１から７のいずれかの疾病対策支援方法を用いた疾病対策支援システムであって、コンピュータプログラムとして、請求項１から７のいずれかの疾病対策支援方法を有し、かつ、個人に関する疾病症状時系列データと環境に関する時系列データとを有し、症状に応じた対策情報を提供するサーバと、通信ネットワークを介して前記サーバから症状に応じた対策情報を受信する端末とを備えることを特徴とする疾病対策支援システムである。
請求項１０の発明は、請求項１から７のいずれかの疾病対策支援方法を用いた疾病対策支援システムであって、環境に関する時系列データを有するサーバと、個人に関する疾病症状時系列データを有し、通信ネットワークを介して前記サーバから環境に関する時系列データを受信し、かつ、コンピュータプログラムとして請求項１から７のいずれかの疾病対策支援方法を有し、症状に応じた対策情報を提供する端末とを備えることを特徴とする疾病対策支援システムである。
請求項１１の発明は、請求項９または１０に記載の疾病対策支援システムにおいて、対象とする疾病は花粉症であることを特徴とする。
【０００８】
本発明によれば、本発明を利用する患者、例えば花粉症患者に信頼性が高く確実な疾病対策支援方法および疾病対策支援システムの提供を可能とする。これにより、利用するユーザは将来の症状の程度を予測することができ、服薬の種類や量、あるいは、旅行目的地に応じた症状対策をすることができる。具体的には、センサによって計測した過去の花粉飛散量といった、環境に関する時系列データおよび症状の時系列データと、将来の環境に関する時系列データとから、将来の症状の予測式を算出し、あるいは、データを分類することにより、症状を予測することを特徴とする。
【０００９】
また、本発明の疾病対策支援方法および疾病対策支援システムは、予測症状と実際の症状との誤差を用いて、症状の予測式の各係数を更新することを特徴とする。実際の症状は、センサで測定したもの、あるいは、自分が感じた症状の程度を手入力にてシステムに入力したものである。症状レベルには個人差があるため、ユーザが設定可能である。これにより、花粉などに対する感受性の経時変化にシステムを追従させ、常にユーザの特性に合った症状の予測性能を高めることができる。
【００１０】
また、本発明の疾病対策支援方法および疾病対策支援システムは、ユーザの症状をよく表す予測式あるいはデータベースを選択し、ユーザの症状特性を分類する。これにより、ユーザの症状特性をいくつかのパターンに分類し、典型的な症状対策情報を提供することができる。
【００１１】
また、本発明の疾病対策支援方法および疾病対策支援システムは、通信ネットワークによって接続されたサーバと端末とから構成され、コンピュータプログラムとして疾病対策支援方法を有し、かつ、個人に関する疾病症状時系列データと環境に関する時系列データとを有し、症状に応じた対策情報を提供する。これにより、ユーザは、携帯電話などにより容易に疾病対策支援サービスを利用できたり、あるいは端末内に疾病症状データを留めておくことにより、プライバシーを確保した疾病対策支援サービスを利用できるようになる。
【発明の効果】
【００１２】
本発明は、将来の症状を、過去の環境に関する時系列データおよび症状の時系列データと、将来の環境に関する予測データとの関数とみなし、データが得られるごとに、関数に含まれる各係数を推定し、将来の症状の予測式を算出し、最新のデータを予測式に入力し、その値により、症状の予測を可能とする。
【００１３】
また、本発明は、予測した症状と実際の症状との誤差により、症状の予測式を更新することができる。
【００１４】
また、本発明は、過去の症状時系列データと、過去の環境に関する時系列データと、将来の環境に関する予測データとからなるベクトルを構成し、ベクトルが得られるごとに分類し実際の疾病症状に応じてカテゴリ付けしたものをデータベースとし、データベースと比較することにより最新のベクトルを分類し、得られたカテゴリに対応する疾病症状により疾病症状を予測し、症状に応じた対策情報を提供するデータが得られるごとにベクトルを分類することにより、症状の予測を可能とする。
【００１５】
また、本発明は、ユーザの症状をよく表す予測式あるいはデータベースを選択し、ユーザの症状特性を分類する。また、本発明は、疾病症状をセンサによって自動的に記録することにより、症状の予測式を自動的に更新することを可能とする。また、本発明は、疾病症状レベルをユーザが設定することにより、ユーザの特性に合わせた症状予測を可能とする。
【００１６】
さらに、本発明は、コンピュータプログラムとして疾病対策支援方法を有し、かつ個人に関する疾病症状時系列データと環境に関する時系列データを有し、症状に応じた対策情報を提供するサーバと、通信ネットワークを介してサーバから症状に応じた対策情報を受信する端末とを備えることにより、ユーザの症状に応じた高精度なアドバイスを提供することが可能となり、また、ユーザの特性に応じて自動的に予測式を更新することが可能となる。
【００１７】
また、本発明は、環境に関する時系列データを有するサーバと、個人に関する疾病症状時系列データを有し、通信ネットワークを介してサーバから環境に関する時系列データを受信し、かつコンピュータプログラムとして疾病対策支援方法を有し、症状に応じた対策情報を提供する端末からなることにより、疾病歴などユーザのプライバシーを保護した高精度なアドバイスを提供することが可能となる。
【発明を実施するための最良の形態】
【００１８】
本発明の疾病対策支援方法および疾病対策支援システムについて実施例を挙げ、図に基づいて説明する。図１は本疾病対策支援方法のフローチャートの例である。以下の説明では、各符号を次のとおりとする。
a[t]を、時刻ｔにおける症状の程度を数値化した症状データ、
b[t]を、時刻ｔにおける大気浮遊物濃度など数値化した環境データ、
a〜[t+1]を、時刻t+1における症状予測データ、
b〜[t+1]を、時刻t+1における環境予測データ、
Ｃを、実数からなる係数ベクトル、
とする。症状データa[t]については個人差があるため、ユーザが自分に最も適合するレベルパターンを選択できる。
【００１９】
疾病対策支援方法は、前記の各データとベクトルとから、次のような処理を行う。つまり、症状データa[t]に関しては、時刻t-n+1まで遡ったｎ個のデータ、環境データb[t]に関しては、時刻t-m+1まで遡ったｍ個のデータ、および、時刻ｔ＋１における環境予測データb〜[t+1]１個からなるベクトルをＸと表す（ステップＳ１１）。
Ｘ=(a[t],a[t-1],…,a[t-n+1],b〜[t+1],b[t],b[t-1],…,b[t-m+1])Ｔ  …式１
ここで、Ｔは転置を表す。症状の予測式を線形関数で構成する場合、
Ｃ=(c1,c2,…,cn+m+1)Ｔ    …式２
を用いて、
a〜[t+1]=ＸＴＣ    …式３
と表すことができる。
【００２０】
式３において、係数ベクトルＣは次のように求める。ある期間のデータから行列Ａとベクトルｙとを、
Ａ1=(a[t-1],a[t-2],…,a[t-n],b[t],b[t-1],…,b[t-m])Ｔ
Ａ2=(a[t-2],a[t-3],…,a[t-n-1],b[t-1],b[t-2],…,b[t-m-1])Ｔ
・
・
Ａｋ=(a[t-k],a[t-k-1],…,a[t-k-n+1],b[t-k+1],b[t-k],…,b[t-k-m+1])Ｔ
Ａ=(Ａ1,Ａ2,…,Ａk)Ｔ    …式４
ｙ=(a[t],a[t-1],…,a[t-k+1])Ｔ    …式５
と定める（ステップＳ１２、ステップＳ１３）。これらのトレーニングデータＡとｙとを用いて、係数ベクトルＣは、
Ｃ=(ＡＴＡ)-１ＡＴｙ    …式６
のように計算することができる（ステップＳ１４）。ここで、−１は逆行列を表す。
【００２１】
新たに環境データおよび症状データが得られ、ベクトルＸを構成すると（ステップＳ１５）、得られた係数ベクトルＣより式３を求め、症状を予測する（ステップＳ１６）。次の時刻に、再びステップＳ１１から繰り返す。
【００２２】
ところで、症状の予測式を非線形関数で構成する場合、時刻t+1における症状予測データを、
a〜[t+1]=f(Ｘ,Ｗ), Ｗ=(w1,w2,…,wj)Ｔ    …式７
と定める。ここで、f(Ｘ,Ｗ)は、ＸとＷの非線形関数、Ｗは係数ベクトルである。この係数ベクトルＷを、ある期間のトレーニングデータＸと、時刻ｔにおける症状の程度を数値化した症状データa[t]と、時刻t+1における症状予測データa〜[t]とを用いて更新する。
【００２３】
図２は、本疾病対策支援方法のフローチャートの例である。まず、環境データb[t]と症状データa[t]とよりベクトルＸを構成する（ステップＳ２１）。予測式による病状予測データa〜[t]と実際の病状データa[t]との差の自乗値をある期間総和した、
Ｅ=1/2Σ(a〜[t]-a[t])2    …式８
を用いて、係数ベクトルＷを以下のように更新する（ステップＳ２２、Ｓ２３、Ｓ２４）。
dＥ/dＷ=Σ(a〜[t]-a[t]), ∂f(Ｘ,Ｗ)/∂Ｗ, ΔＷ=-αdＥ/dＷ    …式９
ここで、αは正の実数、ΔＷはＷの更新量である。
【００２４】
このように、予測値と実際の値との誤差(a〜[t]-a[t])を用いて、係数ベクトルＷ、すなわち関数ｆを更新できるので、関数ｆを毎時刻更新することにより、常に個人の特性に合わせた最新の予測式に更新できる。予測式の値が実際の値に等しい場合は誤差が０となり、予測式の更新は行われない。関数ｆとしては、ニューラルネット（例えば参考文献１）などが考えられ、その場合、係数ベクトルＷはニューラルネットの結合荷重値および閾値である。新たに環境データおよび症状データが得られ、ベクトルＸを構成すると（ステップＳ２５）、得られた係数ベクトルＷより式７を求め、症状を予測する（ステップＳ２６）。次の時刻に再びステップＳ２１から繰り返す。
参考文献１：パターン識別  R. O. Duda, P. E. Hart, and D. G. Stork John Wiley & Sons 2001
【００２５】
さらに、ベクトルＸの値そのものを分類することにより、症状を予測することも可能である。この場合、分類するカテゴリは、t+1での症状a[t+1]とする。図３は、本疾病対策支援方法のフローチャートの例である。まず、環境データb[t]と症状データa[t]とより、ベクトルＸを構成する（ステップＳ３１）。ベクトルＸにカテゴリを対応付け、ある期間のデータを収集して、トレーニングデータとする（ステップＳ３２、ステップＳ３３）。
【００２６】
新たに環境データおよび症状データが得られ、ベクトルＸを構成すると（ステップＳ３４）、それをカテゴリが未知の新たなＸとして属するカテゴリを識別し（ステップＳ３５）、症状を予測する（ステップＳ３６）。トレーニングデータは、ベクトルＸを分類するためのデータベースとも言える。ステップＳ３５での識別方法としては、判別分析法、サポートベクトルマシン法（例えば前記の参考文献１）などを適用することができる。次の時刻に、再びステップＳ３１から繰り返す。
【００２７】
前記のいずれにおいても、初めてトレーニングを行う場合、ある期間は、トレーニングデータを取得するだけで、予測を行うことができない。しかし、この間は例えば、過去のデータより自分に最も適合する予測式あるいはトレーニングデータのテンプレートを選択し、予測式やトレーニングデータが得られるまで、一時的に利用することなどが考えられる。
【００２８】
ユーザの疾病特性を典型的ないくつかのタイプに分類することにより、疾病対策を行うことも有益である。たとえば、花粉症の場合、花粉の飛散シーズンになると、すぐに症状が悪化する人もいるし、徐々に症状が悪化する人もいる。前もって自分がどのようなタイプか分かっていれば、花粉飛散シーズン前に薬を利用するなどの対策を講じることも可能となる。このように分類する方法としては、症状の予測式を利用するものと、Ｘの値を分類するデータベースを利用するものとがある。
【００２９】
予測式を利用する場合について、図４の本疾病対策支援方法のフローチャートの例を用いて説明する。症状の予測式を利用する場合、症状のタイプにより複数の予測式を用意し、それぞれＡ、Ｂ、Ｃなどの指標を付けておく。たとえば、花粉症の場合、
指標Ａは、花粉飛散シーズンになるとすぐに症状が現れるが、症状は弱いタイプ、
指標Ｂは、花粉飛散シーズンになってもすぐに症状は現れないが、少しすると症状が悪化するタイプ、
指標Ｃは、初めからひどい症状のタイプ、
などである。予測式は、線形でも非線形でも構わないが、典型的な症状の何人かのユーザのデータより作成したものが好ましい。
【００３０】
以下、ユーザの症状特性の分類方法について述べる。まず、環境データb[t]と症状データa[t]とより、ベクトルＸを構成する（ステップＳ４１）。ベクトルＸを指標付けした複数の予測式に入力し（ステップＳ４２）、実際の症状に一番近い値を出力した予測式を選択する（ステップＳ４３）。選択された予測式の指標にユーザを分類する（ステップＳ４４）。
【００３１】
データベースを利用する場合について、図５の本疾病対策支援方法のフローチャートの例を用いて説明する。ベクトルＸの値を分類するデータベースを利用する場合は、症状のタイプにより複数のデータベースを用意し、それぞれＡ、Ｂ、Ｃなどの指標を付けておく。データベースは、典型的な症状の何人かのユーザのデータより作成したものが好ましい。
【００３２】
以下、ユーザの症状特性の分類方法について述べる。まず、環境データb[t]と症状データa[t]とよりベクトルＸを構成する（ステップＳ５１）。ベクトルＸを指標付けした複数のデータベースで分類し（ステップＳ５２）、分類されたカテゴリが実際の症状と適合したデータベースを選択する（ステップＳ５３）。選択されたデータベースの指標にユーザを分類する（ステップ５４）。
【００３３】
前記のいずれの分類法において、使用するデータは１日だけのデータに限らず、まとまったある期間のデータがあれば、トータルで最も良く適合した予測式やデータベースを選択する。これによって、ユーザの症状特性タイプをより正確に推定することが可能となる。
【００３４】
また、前記いずれにおいても、実際の症状の値は、ユーザが感じた症状の程度を手入力にてシステムに投入したり、あるいは症状を計測できるセンサ、たとえば、花粉症であればくしゃみの回数を計測できるマイク、加速度センサ、振動センサなどをシステムに接続し、自動的に入力できる。症状については個人差があるため、ユーザは自分に適合した症状レベルパターンを選択できるようにする。環境データについては、ユーザの所在位置の近くに設置された環境センサや、ユーザが携帯できる環境センサから得られるものである。
【００３５】
前記の疾病対策支援方法を用いた疾病対策支援システムについて説明する。この疾病対策支援システムの例を図６に示す。図６の疾病対策支援システムは、サーバ１と、通信回線１０を介してサーバ１に接続可能な端末２と、端末２に接続されているセンサ３とで構成されている。
【００３６】
この疾病対策支援システムでは、端末２は個人が所持する。センサ３は、症状や環境を計測する。センサ３として、たとえば、くしゃみの回数を測定するマイクロホンや加速度センサや振動センサ、あるいは、花粉飛散量を計測するセンサがある。センサ３による計測結果である個人の症状データや環境データは、端末２により、通信回線１０を経由してサーバ１に送信される。
【００３７】
センサ３を用いずに、手入力にてユーザが症状の程度を入力することも、もちろん可能であるし、環境センサについても、端末に接続せずに、ユーザの所在位置の近くに設置されたサーバ１にデータ送信可能なものを利用することもできる。
【００３８】
サーバ１は、過去の環境データと環境予測データとからなる環境データ１Ｂと、個人の症状データ１Ｃとを保持し、疾病対策支援プログラム１Ａによって症状の予測を行う。サーバ１は、その予測に基づいた症状対策情報を、通信回線１０によって個人の所持する端末２に送信する。端末２は、データ通信と表示とのできる簡易なもの、たとえば、携帯電話のようなものでよい。
【００３９】
疾病対策支援システムのもう一つの例を図７に示す。図７の疾病対策支援システムは、サーバ２１と、通信回線３０を介してサーバ２１に接続可能な端末２２と、端末２２に接続されているセンサ２３とで構成されている。サーバ２１は、過去の環境データおよび環境予測データからなる環境データ２１Ａのみ保持し、通信回線３０によって、環境予測データを個人の端末２２に送信する。
【００４０】
個人の所持する端末２２には、症状を測定するセンサ２３、たとえば、花粉症によるくしゃみの回数を測定するマイクロホンや加速度センサや振動センサが接続されている。センサ２３を用いずに、手入力にてユーザが症状の程度を入力することも、もちろん可能である。端末２２は、疾病対策支援プログラム２２Ａを内蔵し、受信した環境データ２１Ａと、個人の症状データ２２Ｂとから症状の予測を行い、それに基づいた症状対策情報を個人に提供する。個人の症状データは、個人の所持する端末２２以外からは流出しないので、個人のプライバシーは保護される。端末２２は、データ通信と表示との他に、メモリと計算機能とを有する簡易なもの、例えばＰＤＡ（Personal Digital Assistance：携帯情報端末）のようなものがある。
【００４１】
以下では、本実施形態による計算例について説明する。
（計算例１）
今、花粉の飛散量を、
０：少
１：中
２：多
とし、症状を、
０：軽
１：中
２：重
とする。花粉の飛散量に関して、次のような時系列データが得られた。
0,0,1,0,0,1,1,1,2,1,2,2,1,2,1,1,1,0,1,0,1,1,2,2,2,1,1,1,0,1    …式１０
ここで、データは１日ごと全３０日のデータで、式１０の左側ほど新しいデータである。
【００４２】
式１０のデータに対応する症状に関して、次のような時系列データが得られた。
0,0,0,1,1,1,1,1,2,2,2,1,1,1,1,1,1,1,1,1,2,2,1,1,0,0,0,0,0,0    …式１１
【００４３】
もちろん、飛散量や症状のレベルのパターンは一例に過ぎず、レベル数を変えてもよいし、飛散量と症状とのレベル数を同じにする必要もない。特に、症状のレベルについては個人差があるため、個人の症状に最も適合したレベルのパターンを選択できる。ここでは軽、中、重の３レベルであるが、これを軽、中、やや重、重、最悪といったような５レベルにすることも考えられる。このデータを用いて、症状の予測モデルを構築する。なお、飛散量は正しく予測できるものとする。式１において、例えば、ｎ＝３、ｍ＝２、すなわち
Ｘ=(a[t],a[t-1],a[t-2],b〜[t+1],b[t],b[t-1])Ｔ    …式１２
とし、式４において、例えばｋ＝１０、すなわち１０日間のデータを用いて、
Ｃ=(c1,c2,…,c6)Ｔ    …式１３
を推定した。
【００４４】
まず、左から１８番目以降のデータを用いて、係数ベクトルＣを推定した。その結果、
Ｃ=(-0.037, -0.023, 0.370, -0.454, 0.385, 0.725)Ｔ    …式１４
となった。この場合、１８番目以降１０日間の症状（式１１）に関しては、
1,1,1,2,2,1,1,0,0,0    …式１５
のところを、
1,1,1,2,2,1,1,0,1,0    …式１６
と推定している。ここで、予測値は、式３の値を小数点以下四捨五入して、整数化している。得られた係数ベクトルＣを用いて、翌日の症状（１７番目のデータ）を予測すると
１    …式１７
となり、実際の症状を予測できたことになる。翌日、ｋを１つ増やし、すなわちｋ＝１１として、同様に係数ベクトルＣを更新し、症状を予測する。ｋ＝１０のまま新しいデータを追加し、一番古いデータを棄却して、計算することもできる。これを毎日繰り返す。
【００４５】
このようにして得られた１番目から１７番目までの症状の予測は、
0,0,1,1,1,1,1,2,2,2,1,1,1,1,1,0,1    …式１８
となり、実際の症状
0,0,0,1,1,1,1,1,2,2,2,1,1,1,1,1,1    …式１９
と比較すると、ほぼ症状を予測できることが分かる。
【００４６】
（計算例２）
同様に、ニューラルネットにて症状の予測モデルを構築した。ニューラルネットは３層とし、入力層はＸの次元数６に対応した６素子で、その関数は、
y(x)=ｘ
中間層は３素子で、その関数は、
y(x)=1/(1+e-(x-z))
ここで、ｚは閾値である。出力層は１素子で、その関数は、
y(x)=ｘ
である。出力層の出力が症状の予測値となる。
【００４７】
まず、１８番目以降１０日間の症状に関して５００００回学習を行い、式７、式９の係数ベクトルＷに対応する結合荷重値および閾値を決定した。入力層のｉ番目の素子と中間層のｊ番目の素子との間の結合加重をｗ１ｉｊとし、中間層のｉ晩目の素子と出力層のｊ番目の素子との間の結合加重をｗ２ｉｊとし、中間層のｉ番目の素子の閾値をｚｉとし、
Ｗ=(w111,w112,w113,w121,w122,w123,…,w163,w211,w221,w231,z1,z2,z3)    …式２０
とすると、得られた係数ベクトルＷは、
Ｗ=(3.848, -0.817, -0.581, 7.343, 1.192, 6.138, -4.437, -0.161, -0.810, 2.167, -0.468, -6.107, 1.984, 0.306, 5.383, -6.569, -0.8126, 1.351, 1.120, -3.027, 1.277, 0.0020, 0.00029, 0.00041)    …式２１
となった。この期間のデータに関しては、実際の症状が、
1,1、1,2,2,1,1,0,0,0    …式２２
であるところを、
1,1,1,2,2,1,1,0,0,0    …式２３
とし、全て正しく学習している。このモデルで翌日の予測を行うと、
１    …式２４
となり、実際の症状を予測できたことになる。
【００４８】
新たなデータを学習データに追加して、２５０００回学習した。毎日このような学習を繰り返す。このようにして得られた１番目から１７番目までの症状の予測は、
0,0,0,1,2,1,1,1,2,2,1,1,1,0,1,1,1    …式２５
となり、実際の症状である
0,0,0,1,1,1,1,1,2,2,2,1,1,1,1,1,1    …式２６
と比較すると、ほぼ症状を予測できることが分かる。
【００４９】
（計算例３）
次に、ベクトルＸを分類することにより、症状を予測する。判別分析法にてベクトルＸの値を分類した。まず、１８番目以降１０日間のデータを用いて、１７番目のデータのカテゴリを識別し、症状を予測する。翌日、１７番目以降１１日間のデータを用いて、１６番目のデータのカテゴリを識別し、同様に症状を予測する。このようにして得られた１番目から１７番目までの症状の予測は、
0,0,1,1,1,1,1,2,2,2,1,1,1,1,1,1,0    …式２７
となる。実際の症状である、
0,0,0,1,1,1,1,1,2,2,2,1,1,1,1,1,1    …式２８
と比較すると、ほぼ症状を予測できることが分かる。
【００５０】
（計算例４）
同様に、サポートベクトルマシン法にてベクトルＸの値を分類した。判別分析法と同様、得られた１番目から１７番目までの症状の予測は、
0,1,1,1,1,1,1,1,2,2,1,1,1,1,1,1,0    …式２９
となり、実際の症状である、
0,0,0,1,1,1,1,1,2,2,2,1,1,1,1,1,1    …式３０
と比較すると、ほぼ症状を予測できることが分かる。
【００５１】
このようにして予測した症状に基づいて、行動アドバイスを生成する。たとえば、症状１と予測したならば、効力が弱めの薬Ａを、症状２と予測したならば、強力な薬Ｂを飲むようなアドバイスを生成する。あるいは、症状１と予測したならば、花粉の予測飛散量が多いＣ地域には行かないほうがよい、症状２と予測したならば、花粉の予測飛散量が少ないＤ地域ならば行っても大丈夫だ、などのアドバイスを生成する。
【００５２】
服薬やマスクなどによって症状が緩和する効果については、新たなパラメタとして、予測式や分類するベクトルに含ませたり、服薬したりマスクをした日の見かけの花粉飛散量を低減させるといった方法により、予測精度を高めることが可能である。たとえば、
１：服薬した、
０：服薬なし
といった、d[t]というパラメタを新たに作成したり、服薬した日は、
b[t]->b[t]-1
とするといったルールを構築することが考えられる。
【００５３】
ベクトルＸを構成する症状時系列データの数ｎと、環境時系列データの数ｍとは、症状予測に影響を及ぼす症状と環境の履歴の度合いと考えることができる。たとえば、花粉症の場合、現時点までの花粉飛散量の総量が花粉症の程度に影響を及ぼすと考えられ、そのような場合は、環境時系列データの数ｍの値を大きくしたり、あるいは新たなパラメタとして、現時点までの花粉飛散量の総量をモデルに付加したりすることが考えられる。
【００５４】
（計算例５）
環境データと症状データとが式１０、式１１で与えられるユーザの症状特性タイプを、予測式あるいはデータベースを利用することにより分類する。Ａ、Ｂ、Ｃという指標を付けた３つの予測式あるいはデータベースを用意する。これらの予測式あるいはデータベースは、典型的な症状パターンを示す複数のユーザのデータを基に、既に述べた方法により作成される。これら３つの予測式あるいはデータベースに、式１０、式１１を入力して予測された症状が、
Ａ:0,0,0,1,1,1,1,1,2,1,2,1,1,1,1,0,1,1,1,1,2,2,1,1,0,0,0,0,0,0
（２箇所誤り）    …式３１
Ｂ:0,0,0,1,1,1,1,1,2,2,1,1,1,0,1,1,0,1,1,1,2,2,1,1,0,1,0,0,0,0
（４箇所誤り）    …式３２
Ｃ:0,1,0,1,1,1,2,2,2,2,2,1,1,1,1,1,1,1,2,2,2,2,1,0,0,0,0,0,0,0
（６箇所誤り）    …式３３
となったとする。この場合は、ユーザの症状特性タイプを、誤り箇所の最も少ない指標Ａとすることができる。このように、ユーザの症状特性タイプが分かると、どのような治療薬を利用するかという判断の支援が可能となったり、長期的な症状の予測が可能となったり、あるいは、同じタイプのユーザの症状対策を参考にし自分も試してみる、などの利点が生じる。
【図面の簡単な説明】
【００５５】
【図１】本発明による疾病対策支援方法を示すフローチャートである。
【図２】本発明による疾病対策支援方法を示すフローチャートである。
【図３】本発明による疾病対策支援方法を示すフローチャートである。
【図４】本発明による疾病対策支援方法を示すフローチャートである。
【図５】本発明による疾病対策支援方法を示すフローチャートである。
【図６】本発明の疾病対策支援方法コンピュータプログラムおよび症状時系列データをサーバが有する疾病対策支援システムの一例を示すブロックである。
【図７】本発明の疾病対策支援方法コンピュータプログラムおよび症状時系列データを端末が有する疾病対策支援システムを示すブロックである。
【符号の説明】
【００５６】
１、２１    サーバ
１Ａ、２２Ａ    疾病対策支援プログラム
１Ｂ、２１Ａ    環境データ
１Ｃ、２２Ｂ    症状データ
２、２２    端末
３、２３    センサ
１０、３０    通信回線
【図１】
2005063218.tif 000003
【図２】
2005063218.tif 000004
【図３】
2005063218.tif 000005
【図４】
2005063218.tif 000006
【図５】
2005063218.tif 000007
【図６】
2005063218.tif 000008
【図７】
2005063218.tif 000009