JP2005052164A - 疾病リスクシミュレーションシステム - Google Patents
疾病リスクシミュレーションシステム Download PDFInfo
- Publication number
- JP2005052164A JP2005052164A JP2003180975A JP2003180975A JP2005052164A JP 2005052164 A JP2005052164 A JP 2005052164A JP 2003180975 A JP2003180975 A JP 2003180975A JP 2003180975 A JP2003180975 A JP 2003180975A JP 2005052164 A JP2005052164 A JP 2005052164A
- Authority
- JP
- Japan
- Prior art keywords
- simulation
- disease
- screen
- values
- simulation system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Images
Landscapes
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Abstract
【課題】様々な疫学的調査・解析を基にして多数年間に渡る総合的なシミュレーションシステムを構築する。
【解決手段】所定の属性値及び検査値を入力する要素データ入力手段と、
所定の疾病に係る複数のステージ間の推移と、所定のイベント発生とを確率的予測手法により予測し、一定期間のイベント、生存年数、及びそれらに係る費用をシミュレーションするシミュレーション演算手段と、
将来における費用に関する基本的データとその効果に関する基本的データを設定する詳細条件設定入力手段と、
上記のシミュレーション演算手段の作出する演算結果を、夫々に対応して出力するシミュレーション結果出力手段とを含む疾病リスクシミュレーションシステムを構築する。
【選択図】図1
【解決手段】所定の属性値及び検査値を入力する要素データ入力手段と、
所定の疾病に係る複数のステージ間の推移と、所定のイベント発生とを確率的予測手法により予測し、一定期間のイベント、生存年数、及びそれらに係る費用をシミュレーションするシミュレーション演算手段と、
将来における費用に関する基本的データとその効果に関する基本的データを設定する詳細条件設定入力手段と、
上記のシミュレーション演算手段の作出する演算結果を、夫々に対応して出力するシミュレーション結果出力手段とを含む疾病リスクシミュレーションシステムを構築する。
【選択図】図1
Description
【0001】
【発明の属する技術分野】
本発明は、疾病に関わる種々のイベントの発生をシミュレーションするシステムに関する。
【0002】
【従来の技術】
疾病に係る疫学的調査及びその解析は、医療及び保健関係者のみならず、保険関係者や医療経済関係者にとっても関心事の一つである。
【0003】
昨今では特に、生活習慣病と冠動脈疾患発症率の関係を検討した「Framingham Heart Study」がよく知られている。また、「UKPDS」、「熊本スタディ」等も、同様の疫学的調査・解析の成果である。
【0004】
それら疫学的調査・解析を基にした、コンピュータシミュレーションシステムが数多く構築され更に発表されつつある。
【0005】
特許文献1では、主として「Framingham Heart Study」を基礎として、将来の心血管疾患発症リスクおよび任意の治療目標を達成した際のリスクをシミュレーションする計算機システムを提供する。該計算機システムとは、被検査者に係る検査値入力値と既往歴入力値を含む危険因子入力値を入力する入力表示画面部と、その入力表示画面部からのデータを受け取り、疫学調査により導出された方程式に従って被検査者の疾患発症危険シミュレーション処理を行う処理部と、それにより生成された結果を出力する出力表示画面部とを有する、疾患治療指導計算機システムである。
【0006】
【特許文献1】
特開2002−24401号公報
【0007】
【発明が解決しようとする課題】
本発明は、様々な疫学的調査・解析を基にして多数年間に渡る総合的なシミュレーションシステムを構築することを目的とするものである。その際、人種差を反映した調査結果の加味、仮想値の設定、費用のパラメータの設定、費用対効果の算出などを考慮する。
【0008】
【課題を解決するための手段】
本発明は、上記の目的を達成するために為されたものである。本発明に係る請求項1に記載の疾病リスクシミュレーションシステムは、
所定の属性値及び検査値を入力する要素データ入力手段と、
所定の疾病に係る複数のステージ間の推移と、所定のイベント発生とを確率的予測手法により予測し、一定期間のイベント、及び生存年数をシミュレーションするシミュレーション演算手段と、
上記のシミュレーション演算手段の作出する演算結果を、夫々に対応して出力するシミュレーション結果出力手段とを含む疾病リスクシミュレーションシステムである。
【0009】
本発明に係る請求項2に記載の疾病リスクシミュレーションシステムは、
シミュレーション演算手段において、ステージ間推移の推移確率として、標本対象の人種が限定された統計に基づく確率値が利用される、
請求項1に記載の疾病リスクシミュレーションシステムである。
【0010】
本発明に係る請求項3に記載の疾病リスクシミュレーションシステムは、
シミュレーション演算手段において、所定のイベント発生のイベント発生率が、標本対象の人種が限定された統計により補正される、
請求項1又は請求項2に記載の疾病リスクシミュレーションシステムである。
【0011】
本発明に係る請求項4に記載の疾病リスクシミュレーションシステムは、
シミュレーション演算手段において、死亡率を他の疾病に係るステージの状態に依存させる、
請求項1乃至請求項3のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステムである。
【0012】
本発明に係る請求項5に記載の疾病リスクシミュレーションシステムは、
シミュレーション演算手段において、所定の疾病に係る複数のステージ間の推移と、所定のイベント発生とに係る費用がシミュレーションされ、
将来における費用に関する基本的データとその効果に関する基本的データを設定する詳細条件設定入力手段を、更に含む、
請求項1乃至請求項4のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステムである。
【0013】
本発明に係る請求項6に記載の疾病リスクシミュレーションシステムは、
上記の詳細条件設定入力手段において、夫々の疾病のイベント又は夫々の疾病の個別のステージに関し、
費用に関する基本的データとして、少なくとも、基本値、追加値及び追加仮想値が設定され、更に、効果に関する基本的データとして、少なくとも、効用値が設定され、
それら設定値がシミュレーション演算手段での将来における費用対効果に係る演算で利用され、
シミュレーション結果出力手段に出力される、
請求項5に記載の疾病リスクシミュレーションシステムである。
【0014】
本発明に係る請求項7に記載の疾病リスクシミュレーションシステムは、
所定の疾病に係る複数のステージが、マルコフモデルに従い形成されており、
確率的統計手法として、モンテカルロシミュレーションの手法が利用される、請求項1乃至請求項6のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステムである。
【0015】
【発明の実施の形態】
以下において、図面を参照しつつ本発明に係る好適な実施の形態を説明する。
【0016】
(1)稼動装置について
図3は、本発明に係る疾患リスクシミュレーションシステム(装置)として稼動し得る、計算機システム2の1つの形態である。本発明として動作する計算機は、構成や能力において特に限定されない。図3の計算機システム2は、表示部4、入力部6、本体部8を備える。さらに、図示していないが、本体部8には、例えばCPUを備えるマザーボードなどからなる処理部12、ハードディスクなどからなるメモリ部10、及びCD−ROM等の外部記憶媒体14に格納されるデータを読み取る読み出し部(図示せず。)が、含まれる。
【0017】
本発明に係る疾患リスクシミュレーションシステム(装置)の処理プログラムは、計算機システム2のメモリ部10に格納されている。若しくは、CD−ROM等の外部記憶媒体14に格納されていてもよい。その場合、上記処理プログラムは、(CD−ROM等の)外部記憶媒体14から読み出し部を介して読み出され、必要に応じてメモリ部10若しくは処理部12に保持される。計算機システム2が外部ネットワーク(図示せず。)を介して外部サーバと接続し、該サーバに用意される本発明に係る疾患リスクシミュレーションの処理プログラムを利用して(ダウンロードして)、本発明に係る疾患リスクシミュレーションを行なう、という形態であってもよい。
【0018】
(2)ソフトウエアの全体概要
以下の説明における本発明に係る疾患リスクシミュレーションシステムは、糖尿病リスクシミュレーションに関するものである。言うまでも無く、本発明に係るリスクシミュレーションの概念は、他の疾患・疾病にも利用され得るものである。
【0019】
本発明に係る疾患リスクシミュレーションシステム、即ち糖尿病リスクシミュレーションシステムは、糖尿病患者の年齢、性別などの基本的な属性データ、及び血糖値、血圧などの検査値データを基にして、当該患者の網膜症、腎症、CHD(coronary heart disease;虚血性心疾患)及び死亡の、糖尿病に関わる主要なイベントの発生を長期間に渡りシミュレーションする。
【0020】
具体的にはシミュレーションにより、生存年数、一定期間におけるイベント(例えば、網膜症における失明や、腎症における透析)の数(ないしは発生率)、それらに伴う費用を推計し、適宜グラフやテーブルに示す。シミュレーションの基本原則は、各合併症のステージ間の推移と、CHD及び死亡の発生とを、確率的に予想するモンテカルロシミュレーションによるものである。
【0021】
(3)表示画面の推移・構成(画面フロー)
図2に示すように、本発明に係る糖尿病リスクシミュレーションシステムでは大きく分類すると、表示部4にて3類型の画面が表示される。つまり、
・初期データ入力画面100
・シミュレーション結果画面200
・詳細設定画面300
の3類型である。後述するように、夫々の類型は更に数種類の画面を含む。
【0022】
図2の画面の全体構成を簡単に説明する。初期データ入力画面100では、糖尿病患者の年齢、性別などの基本的な属性データ、及び血糖値、血圧などの検査値データが入力・設定される。そこでシミュレーションの実施が指示されると、シミュレーションの実施後、シミュレーション結果画面200が表示される。シミュレーション結果画面200は複数種類のシミュレーション結果画面から構成される。
【0023】
糖尿病患者の年齢、性別などの基本的な属性データ、及び血糖値、血圧などの検査値データの値が変更されるときには、通常、初期データ入力画面100へ戻される。ただし、シミュレーション結果画面200も、入力・設定値の変更を為し得る機能を併有する。よって、ある種類のシミュレーション結果画面200から別の種類のシミュレーショ結果画面200に直接移行したり、一つの種類のシミュレーション結果画面200で(入力・設定値を変更しつつ)シミュレーションを繰り返したりすることが可能である。
【0024】
詳細設定画面300では、将来における費用及びその効果に関する基本的データが設定される。初期データ入力画面100では、糖尿病患者個人に係るデータが入力設定されるのに対して、当該詳細設定画面300では患者に拠って変動しない費用及び効果に係るデータが設定・修正される。図2では基本的に初期データ入力画面100と詳細設定画面300との間で画面推移することを示しているが、シミュレーション結果画面200と詳細設定画面300との間で画面推移できるようにも構成されている。
【0025】
(4)ソフトウエア主要機能(4つのサブモデル)
本発明に係る糖尿病リスクシミュレーションシステムを構築するソフトウエアの主要な機能は、
・(糖尿病に関連する合併症である)網膜症の推移と腎症の推移を確率的に推計すること、
・(糖尿病患者でリスクが増大する)CHDの発生と死亡の発生を確率的に推計すること、
である。これらの推計のために、本発明に係る糖尿病リスクシミュレーションシステムは、
(i)網膜症の推移
(ii)腎症の推移
(iii)CHDの発生
(iv)死亡の発生
に関する4つのサブモデルを含む。
【0026】
上記の4つのサブモデルは、疫学データを利用する。利用する疫学データとしては、様々なものが存在し得る。以下の説明では、日本国における熊本スタディ、英国におけるUKPDS、米国におけるFramingham Heart Study及びHonolulu Heart Programを利用するものを挙げている。ここで、Honolulu Heart Programは日系人を対象とするものである(“Ralph B. D’Agostino etal. JAMA. 286. 180−187.2001”参照)。
【0027】
(i)網膜症の推移(網膜症サブモデル)
網膜症サブモデルは、
「網膜症なし」、
「単純網膜症」、
「前増殖網膜症」、
「増殖網膜症」、
「失明」
の5つのステージによって構成されるマルコフモデルである。これらのステージ間で一定の推移確率で、より重症度の高いステージへの推移が発生するとしており、逆戻りは設定していない。推移は、年次の推移を検討する。よって推移確率は年次の推移確率である。
【0028】
上記の推移確率は、熊本スタディのデータを基に設定されている。但し、「増殖網膜症」から「失明」への推移確率は、熊本スタディでは報告されていないため、米国における報告値(Scottら, 1988)が用いられている。
【0029】
「網膜症なし」から「増殖網膜症」までの各区間の推移確率は、熊本スタディのデータを基に、HbA1Cによるリスク調整の対象としている。即ち、HbA1C値が悪化する程、「網膜症なし→単純網膜症」、「単純網膜症→前増殖網膜症」、「前増殖網膜症→増殖網膜症」の各区間における網膜症の進展が早くなることになる。リスク調整については、後で説明する。
【0030】
更に、「単純網膜症→前増殖網膜症」への推移確率は、UKPDSのデータを基に収縮期血圧(SBP)によるリスク調整も行う。
【0031】
(ii)腎症の推移(腎症サブモデル)
腎症のサブモデルは、
「腎症なし」、
「微量アルブミン尿期」、
「顕性蛋白尿期」、
「腎不全保存期」、
「血液透析」
の5つのステージによって構成されるマルコフモデルである。網膜症サブモデルと同様、これらのステージ間で一定の確率で、より重症度の高いステージへの推移が発生するとしており、逆戻りは設定していない。推移はやはり年次の推移を検討するため、推移確率は年次の推移確率である。それら推移確率は、熊本スタディのデータを基に設定されている。
【0032】
「腎症なし」から「顕性蛋白尿期」までの各区間の推移確率は、熊本スタディのデータを基に、HbA1Cによるリスク調整の対象としている。更に、「腎症なし」から「顕性蛋白尿期」までの各区間の推移確率は、UKPDSのデータを基に、収縮期血圧(SBP)によるリスク調整も行う。
【0033】
リスク調整について説明する。上記(i)(ii)で示したリスク調整は、本発明では例えば、熊本スタディ(血糖コントロール)及びUKPDS(血圧コントロール)の結果を下記の式に代入して実施する。なお、下記の式は“Jonathan Betz Brown et al. Diabetes research and clinical practice.50(Suppl.3).s15−s46,2000”に示されるものである。
【数1】
RRi=RRt^(ln(vi/vm)/ln(vc/vt))
RRi:i番目にシミュレーションされる患者の調整係数
RRt:比較試験で治療の結果観察された平均相対危険
vi:i番目の患者のリスク因子の数値
vm:推移確率の基本値を計算した集団におけるリスク因子の平均値
vt:比較試験で治療群が達成したリスク因子の平均値
vc:比較試験でのコントロール群のリスク因子の平均値
vmについては熊本スタディにおける平均値(HbA1C:9.125、SBP:121)を使用する。vt、vcについては、HbA1Cによる調整時は7.2、9.4を、SBPによる調整時は144、154をそれぞれ使用する。相対リスク(RRt)は、下の表1と表2のように、熊本スタディにおける従来インスリン療法群と強化インスリン療法群の相対リスクを、またSBP修正時の相対リスク(RRt)はUKPDSにおける従来療法群と介入群との相対リスクを使用する。
【0034】
【表1】
【表2】
【0035】
例えば、HbA1C:8%のときの「網膜症なし」から「単純網膜症」への推移確率の調整は以下のようになる。
【数2】
0.052=3.72^(ln(8/9.125)/ln(9.4/7.2))
従って、推移確率は基本値の9.3%をかけて
【数3】
9.3%×0.052=0.48%
となる。この例のような「網膜症なし」から「単純網膜症」への場合はSBPによる調整は必要ないが、SBPによる調整が必要な場合にはHbA1Cによる調整のあとに更にSBPによる調整を行なう。
【0036】
(iii)CHDの発生(虚血性心疾患サブモデル)
CHD(虚血性心疾患)発生率の計算では、「Framingham Heart Study」の解析結果に基づくリスク計算式を基本としている。このリスク計算式では、年齢、性別、総コレステロール、HDLコレステロール、SBP、喫煙習慣の有無、耐糖能異常の有無及び左室肥大の有無が、リスクを計算するための因子とされる。
【0037】
但し、人種差を考慮し、日系人を対象とした「Honolulu HeartProgram」の結果を参考に補正を行っている。
【0038】
例えば、以下のAndersonの計算式により求められる発症確率pに一律「1/2.3」を掛けるという補正である。
【0039】
なお、「Framingham Heart Study」の解析結果に基づくリスク計算式については、いくつか提示されている。本発明では、例えば、Andersonらにより提示されている式を用いる。Andersonらの計算式は以下の通りである。
【0040】
{smoking:喫煙習慣ありの場合1、それ以外0}
{LVH:左室肥大ありの場合1、それ以外0}
{DM:耐糖能異常ありの場合1、それ以外0}
と定義する。ここで、本発明に係る糖尿病リスクシミュレーションシステムは、“耐糖能異常あり”の患者のデータを取り扱うことが前提となっている、として、“DM=1”と固定されていてもかまわない。
続いてまず、中間指標の「a」、「m」、「μ」、「σ」、「u」を算出する。
【0041】
aの計算;
【数4】
a=11.1122−0.9119*log(SBP)−0.2767*smoking−0.7181*log(TC/HDL)−0.5865*LVH
【0042】
mの計算;
男の場合
【数5】
m=a−1.4792*log(age)−0.1759*DM
女の場合
【数6】
m=a−5.8549+1.8515*[log(age/74)]^2−0.3758*DM
【0043】
μの計算;
【数7】
μ=4.4181+m
【0044】
σの計算;
【数8】
σ=exp(−0.3155−0.2784*m)
【0045】
uの計算;
【数9】
u=(log(t)−μ)/σ
【0046】
以上の5つの中間指標a、m、μ、σ、uによって、t年後のCHD発症確率pは下記の式で求められる。
【数10】
p=1−exp(−exp(u))
【0047】
(iv)死亡の発生(死亡サブモデル)
死亡率は、厚生労働省発表の簡易生命表における年齢・性別死亡率を基本値としている。それに加えて、「顕性蛋白尿期」では自然死亡率の1.5倍、「腎不全保存期」は2倍とする。更に、「血液透析」では一律14.1%/年の死亡率とし、虚血性心疾患発症後は16.2%/年の死亡率とする。
【0048】
上記の簡易生命表以外の倍率は、“日本内科学会雑誌 87(7):81−86、1998”に基づく。
【0049】
(5)モンテカルロシミュレーションフローチャート
図1にて、本発明の糖尿病リスクシミュレーションシステムに係るモンテカルロシミュレーションの概略のフローチャートを示す。シミュレーション回数(n)は操作者が自由に設定し得る。図1では「n=10000」としている。図1により、モンテカルロシミュレーションの概略の処理の流れを説明する。
【0050】
(ステップS00):シミュレーションを開始する。
【0051】
(ステップS02):n回目の(n番目の)シミュレーションを開始する。まず必要パラメータを初期化する。例えば、
死亡フラグ=0
網膜症ステージ=R1(「網膜症なし」)
腎症ステージ=N1(「腎症なし」)
CHDフラグ=0(「無し」)
等である。
【0052】
(ステップS04):網膜症サブモデル計算を行う。網膜症ステージ毎に設定されている推移確率trと、プログラムで発生させる乱数ranとを比較し、
tr>ran
ならば網膜症ステージをひとつ進展させる。
【0053】
(ステップS06):腎症サブモデル計算を行う。腎症ステージ毎に設定されている推移確率trと、プログラムで発生させる乱数ranを比較し、
tr>ran
ならば腎症ステージをひとつ進展させる。
【0054】
(ステップS08):CHDサブモデル計算を行う。フラミンガムリスク推計式により、CHD発生率を計算する。但し前述のように、「Honolulu Heart Program」の結果による補正をかけている。発生があれば、CHDフラグを“1”にする。更に、CHDフラグ=1のとき、すなわちCHD発生時には、死亡率「16.2%」と、プログラムで発生した乱数ranを比較し、
16.2>ran
なら死亡フラグを1にする。
【0055】
(ステップS10):死亡サブモデル計算を行う。各年の死亡率は簡易生命表による年齢・性別死亡率(以下、“自然死亡率”という。)を基本値として、腎症サブモデルで計算された腎症ステージに応じて死亡率trを以下の表3のとおり求める。
【表3】
死亡率trと、プログラムで発生した乱数ranを比較し、
tr>ran
なら死亡フラグを1にする。
【0056】
(ステップS12):死亡判定を行う。死亡フラグが1なら死亡と判定し、ステップS14に進む。死亡フラグが0ならば、ステップS04に戻り、上記の4つのサブモデルを繰り返す。
【0057】
即ち、n番目のシミュレーションについては、死亡フラグが1となるまで繰り返されることになる。
【0058】
(ステップS14):シミュレーション回数を判定する。シミュレーション回数が設定値(n=10000)に達したら、シミュレーションを終了するためステップS18に進む。シミュレーション回数が設定値に未達であれば、シミュレーション回数をインクリメントして(ステップS16)、次回数のシミュレーションを行う(ステップS02)。
【0059】
(ステップS18):シミュレーションを終了する。
【0060】
(6)画面について
次に、表示画面について説明する。
【0061】
図4は、初期データ入力画面100である。この画面では、糖尿病患者の年齢、性別などの基本的な属性データ、血糖値、血圧、コレステロールなどの検査値データ、及び喫煙習慣・左室肥大などのデータが入力・設定される。糖尿病の合併症のうち、網膜症と腎症の状態についても入力する。シミュレーション開始ボタン108でシミュレーションを実施する。また、詳細設定ボタンで、医療費などのパラメータ設定を行なう詳細設定画面300に進む。
【0062】
図5は、シミュレーション結果画面200のうちの網膜症シミュレーション画面202である。初期データ入力画面で入力した値(現在値)によってシミュレーションされた今後30年間にわたる網膜症の進展予測が表示される。ここで、「腎症」ボタン212を選択すれば、腎症進展についての同様のグラフ(腎症シミュレーション画面(図示せず。))が表示される。
【0063】
上記の網膜症シミュレーション画面202は、「網膜症」ボタン210を表示する画面にて該「網膜症」ボタン210を選択することによっても表示される。同様に、上記の腎症シミュレーション画面(図示せず。)も、「腎症」ボタン212を表示する画面にて該「腎症」ボタン212を選択することによっても表示される。
【0064】
図6は、シミュレーション結果画面200のうちのイベント発生率シミュレーション画面222である。各画面の「イベント発生率」ボタン214が選択されると、初期データ入力画面100で入力された値(現在値)によってシミュレーションされた生涯にわたり発生が予測される各種イベント発生率が表示される(左部226)。
【0065】
図7は、シミュレーション結果画面200のうちの医療費シミュレーション画面232である。各画面の「医療費」ボタン216が選択されると、初期データ入力画面で入力した値(現在値)によってシミュレーションされた生涯にわたり発生が予測される糖尿病合併症に係る医療費が表示される。この医療費の計算において基礎のデータとなるのは、
・初期データ入力画面100で入力された値(現在値)、
・図1に示すシミュレーションフローチャートにより計算され、図5の網膜症シミュレーション画面202、図6のイベント発生率シミュレーション画面222、及び腎症シミュレーション画面(図示せず。)などに表示される確率(発生率)、
・詳細設定画面300で入力される合併症別のステージ毎の年間費用及び合併症別のイベント発生時の年間費用、
である。
【0066】
図8は、シミュレーション結果画面200のうちのシミュレーション一覧画面242である。各画面の「一覧」ボタン218が選択されると、初期データ入力画面100で入力された値(現在値)によってシミュレーションされたイベント発生率、費用対効果に関する情報が一覧表示される。費用対効果については、後で説明する。
【0067】
図9も、シミュレーション結果画面200のうちの網膜症シミュレーション画面202である。ところで、「仮想値」の入力設定欄には、現在値と異なるデータが入力・設定され得る。そうすると、現在値によるシミュレーション結果(画面左部)を表示したまま、仮想値シミュレーションを実施することができる。図9の表示例では、HbA1c、SBP、TCを現在値よりも少なくしてシミュレーションを実施しようとしている。「シミュレーション開始」ボタンで仮想値によるシミュレーションが実施される。
【0068】
図10も、シミュレーション結果画面200のうちの網膜症シミュレーション画面202であるが、入力された仮想値によってシミュレーションされた30年間にわたる網膜症の進展予測が、すでに実施された現在値によるシミュレーション結果の隣に表示される。ここで、「腎症」ボタン212を選択すれば、腎症進展についての同様のグラフ(腎症シミュレーション画面(図示せず。))が表示される。
【0069】
図11も、図6と同様のシミュレーション結果画面200のうちのイベント発生率シミュレーション画面222である。但し図6と異なり、仮想値によりシミュレーションされた生涯にわたる各種イベント発生率が表示されている。このように、「仮想値」の入力設定欄に現在値と異なるデータが入力・設定された状態で各画面の「イベント発生率」214が選択されると、現在値シミュレーション結果の隣に仮想値によりシミュレーションされた生涯にわたる各種イベント発生率が表示される。
【0070】
図12も、図7と同様のシミュレーション結果画面200のうちの医療費シミュレーション画面232である。但し図7と異なり、仮想値によりシミュレーションされた生涯にわたる各種医療費が表示されている。このように、「仮想値」の入力設定欄に現在値と異なるデータが入力・設定された状態で各画面の「医療費」ボタン216が選択されると、現在値シミュレーション結果の隣に仮想値によりシミュレーションされた生涯にわたる各種医療費が表示される。
【0071】
図13も、図8と同様のシミュレーション結果画面200のうちのシミュレーション一覧画面242である。但し図8と異なり、特に仮想値を利用してシミュレーションされたイベント発生率、費用対効果に関する情報が一覧表示されている。このように、「仮想値」の入力設定欄に現在値と異なるデータが入力・設定された状態で各画面の「一覧」ボタン218が選択されると、現在値及び仮想値によってシミュレーションされたイベント発生率、費用対効果に関する情報が一覧表示される。
【0072】
図13において、ARR(%)とRRR(%)とは、以下の式で定義される。
ARR(%)=仮想値(%)−現在値(%)
RRR(%)=100×{仮想値(%)−現在値(%)}/現在値(%)
また、「生存年」は、シミュレーションにおいて計算される余命年の平均値である。「質調整生存年」は、上記「生存年」において、糖尿病合併症別の各ステージ推移やイベント発生が生じたときに、図14に示される「効用値」を生存年として累積したものである。
【0073】
更に図13において、費用対効果(生存年)と費用対効果(質調整生存年)とは、以下の式で定義される。
費用対効果(生存年)=総医療費/生存年
費用対効果(質調整生存年)=総医療費/質調整生存年
増分費用対効果(生存年)と増分費用対効果(質調整生存年)とは、以下の式で定義される。
増分費用対効果(生存年)=増分総医療費/増分生存年
増分費用対効果(質調整生存年)=増分総医療費/増分質調整生存年
ここで、「増分総医療費」は、図13に示す通り、仮想値の総医療費から現在値の総医療費を引いた値である。
【0074】
図14は、詳細設定画面300である。各画面の「詳細設定」ボタンが選択されると、表示・設定される。シミュレーションで使用する費用パラメータ(304、306、308)、効用値310、又はシミュレーションにおけるシミュレーション回数312を変更することができる。更に、割引率設定ボタン314により、割引率を設定する画面(図示せず。)に移行することもできる。
【0075】
ここで「割引率」とは、費用と質調整生存年とを、一定の年率で現在価値に割り戻すための数値である。
【0076】
詳細設定画面300では、費用パラメータとして合併症基本診療費304と追加診療費とを区別して入力・設定可能なように構成されている。更に、追加診療費として現在値306と仮想値308とを区別して入力設定可能なように構成されている。
【0077】
例えば、糖尿病合併症の各ステージでの基本的な医療費を「合併症基本診療費」に設定し、加えて、初期データ入力画面100での現在値を維持したりシミュレーション結果画面200での仮想値を達成したりするために必要な診療費を「追加診療費・現在値」「追加診療費・仮想値」に適宜設定する、というようなことができる。
【0078】
以上、本発明の具体的な実施の形態を説明したが、本発明はこれらに限られず更に種々改変可能である。
【0079】
上記の本発明に係る実施の形態では、糖尿病合併症に関するリスクシミュレーションのためのシステムを取り上げてきた。このシステムは、その他の疾患・疾病に関するリスクシミュレーションのためにも応用可能である。
【0080】
既述のように、本発明の実施の形態は、計算機システム2が外部ネットワーク(図示せず。)を介して外部サーバと接続し、該外部サーバに用意される本発明に係る疾患リスクシミュレーションの処理プログラムを利用して(ダウンロードして)、本発明に係る疾患リスクシミュレーションを行なう、というようなものであってもよい。そうすると、該外部サーバはシミュレーション手段を提供するとともに、フィードバックデータを大規模に収集することも可能になる。
【0081】
【発明の効果】
本発明を利用することにより、以下のような効果を得ることができる。
【0082】
様々な疫学的調査・解析を基にして、多数年間(数十年間)に渡る累積的なシミュレーションを行なうことができる。その際、人種差を考慮したシミュレーションを実施できる。
【0083】
また、費用のパラメータの設定、及び費用対効果の算出などにより、経済的効用の視点からのシミュレーションを可能にする。
【図面の簡単な説明】
【図1】糖尿病リスクシミュレーションシステムのモンテカルロシミュレーションの概略のフローチャートである。
【図2】糖尿病リスクシミュレーションシステムの画面の分類である。
【図3】疾患リスクシミュレーションシステム(装置)として稼動し得る計算機システムの1つの形態である。
【図4】初期データ入力画面である。
【図5】シミュレーション結果画面のうちの網膜症シミュレーション画面である。
【図6】シミュレーション結果画面のうちのイベント発生率シミュレーション画面である。
【図7】シミュレーション結果画面のうちの医療費シミュレーション画面である。
【図8】シミュレーション結果画面のうちのシミュレーション一覧画面である。
【図9】「仮想値」の入力設定欄に仮想値が入力された、シミュレーション結果画面のうちの網膜症シミュレーション画面である。
【図10】仮想値によってシミュレーションされた30年間にわたる網膜症の進展予測を伴う、シミュレーション結果画面のうちの網膜症シミュレーション画面である。
【図11】仮想値によりシミュレーションされた生涯にわたる各種イベント発生率の表示を伴う、シミュレーション結果画面のうちのイベント発生率シミュレーション画面である。
【図12】仮想値によりシミュレーションされた生涯にわたる各種医療費の表示を伴う、シミュレーション結果画面のうちの医療費シミュレーション画面である。
【図13】現在値及び仮想値によってシミュレーションされたイベント発生率、費用対効果に関する情報が一覧表示された、シミュレーション結果画面のうちのシミュレーション一覧画面である。
【図14】詳細設定画面である。
【符号の説明】
100・・・初期データ入力画面、200・・・シミュレーション結果画面、300・・・詳細設定画面。
【発明の属する技術分野】
本発明は、疾病に関わる種々のイベントの発生をシミュレーションするシステムに関する。
【0002】
【従来の技術】
疾病に係る疫学的調査及びその解析は、医療及び保健関係者のみならず、保険関係者や医療経済関係者にとっても関心事の一つである。
【0003】
昨今では特に、生活習慣病と冠動脈疾患発症率の関係を検討した「Framingham Heart Study」がよく知られている。また、「UKPDS」、「熊本スタディ」等も、同様の疫学的調査・解析の成果である。
【0004】
それら疫学的調査・解析を基にした、コンピュータシミュレーションシステムが数多く構築され更に発表されつつある。
【0005】
特許文献1では、主として「Framingham Heart Study」を基礎として、将来の心血管疾患発症リスクおよび任意の治療目標を達成した際のリスクをシミュレーションする計算機システムを提供する。該計算機システムとは、被検査者に係る検査値入力値と既往歴入力値を含む危険因子入力値を入力する入力表示画面部と、その入力表示画面部からのデータを受け取り、疫学調査により導出された方程式に従って被検査者の疾患発症危険シミュレーション処理を行う処理部と、それにより生成された結果を出力する出力表示画面部とを有する、疾患治療指導計算機システムである。
【0006】
【特許文献1】
特開2002−24401号公報
【0007】
【発明が解決しようとする課題】
本発明は、様々な疫学的調査・解析を基にして多数年間に渡る総合的なシミュレーションシステムを構築することを目的とするものである。その際、人種差を反映した調査結果の加味、仮想値の設定、費用のパラメータの設定、費用対効果の算出などを考慮する。
【0008】
【課題を解決するための手段】
本発明は、上記の目的を達成するために為されたものである。本発明に係る請求項1に記載の疾病リスクシミュレーションシステムは、
所定の属性値及び検査値を入力する要素データ入力手段と、
所定の疾病に係る複数のステージ間の推移と、所定のイベント発生とを確率的予測手法により予測し、一定期間のイベント、及び生存年数をシミュレーションするシミュレーション演算手段と、
上記のシミュレーション演算手段の作出する演算結果を、夫々に対応して出力するシミュレーション結果出力手段とを含む疾病リスクシミュレーションシステムである。
【0009】
本発明に係る請求項2に記載の疾病リスクシミュレーションシステムは、
シミュレーション演算手段において、ステージ間推移の推移確率として、標本対象の人種が限定された統計に基づく確率値が利用される、
請求項1に記載の疾病リスクシミュレーションシステムである。
【0010】
本発明に係る請求項3に記載の疾病リスクシミュレーションシステムは、
シミュレーション演算手段において、所定のイベント発生のイベント発生率が、標本対象の人種が限定された統計により補正される、
請求項1又は請求項2に記載の疾病リスクシミュレーションシステムである。
【0011】
本発明に係る請求項4に記載の疾病リスクシミュレーションシステムは、
シミュレーション演算手段において、死亡率を他の疾病に係るステージの状態に依存させる、
請求項1乃至請求項3のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステムである。
【0012】
本発明に係る請求項5に記載の疾病リスクシミュレーションシステムは、
シミュレーション演算手段において、所定の疾病に係る複数のステージ間の推移と、所定のイベント発生とに係る費用がシミュレーションされ、
将来における費用に関する基本的データとその効果に関する基本的データを設定する詳細条件設定入力手段を、更に含む、
請求項1乃至請求項4のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステムである。
【0013】
本発明に係る請求項6に記載の疾病リスクシミュレーションシステムは、
上記の詳細条件設定入力手段において、夫々の疾病のイベント又は夫々の疾病の個別のステージに関し、
費用に関する基本的データとして、少なくとも、基本値、追加値及び追加仮想値が設定され、更に、効果に関する基本的データとして、少なくとも、効用値が設定され、
それら設定値がシミュレーション演算手段での将来における費用対効果に係る演算で利用され、
シミュレーション結果出力手段に出力される、
請求項5に記載の疾病リスクシミュレーションシステムである。
【0014】
本発明に係る請求項7に記載の疾病リスクシミュレーションシステムは、
所定の疾病に係る複数のステージが、マルコフモデルに従い形成されており、
確率的統計手法として、モンテカルロシミュレーションの手法が利用される、請求項1乃至請求項6のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステムである。
【0015】
【発明の実施の形態】
以下において、図面を参照しつつ本発明に係る好適な実施の形態を説明する。
【0016】
(1)稼動装置について
図3は、本発明に係る疾患リスクシミュレーションシステム(装置)として稼動し得る、計算機システム2の1つの形態である。本発明として動作する計算機は、構成や能力において特に限定されない。図3の計算機システム2は、表示部4、入力部6、本体部8を備える。さらに、図示していないが、本体部8には、例えばCPUを備えるマザーボードなどからなる処理部12、ハードディスクなどからなるメモリ部10、及びCD−ROM等の外部記憶媒体14に格納されるデータを読み取る読み出し部(図示せず。)が、含まれる。
【0017】
本発明に係る疾患リスクシミュレーションシステム(装置)の処理プログラムは、計算機システム2のメモリ部10に格納されている。若しくは、CD−ROM等の外部記憶媒体14に格納されていてもよい。その場合、上記処理プログラムは、(CD−ROM等の)外部記憶媒体14から読み出し部を介して読み出され、必要に応じてメモリ部10若しくは処理部12に保持される。計算機システム2が外部ネットワーク(図示せず。)を介して外部サーバと接続し、該サーバに用意される本発明に係る疾患リスクシミュレーションの処理プログラムを利用して(ダウンロードして)、本発明に係る疾患リスクシミュレーションを行なう、という形態であってもよい。
【0018】
(2)ソフトウエアの全体概要
以下の説明における本発明に係る疾患リスクシミュレーションシステムは、糖尿病リスクシミュレーションに関するものである。言うまでも無く、本発明に係るリスクシミュレーションの概念は、他の疾患・疾病にも利用され得るものである。
【0019】
本発明に係る疾患リスクシミュレーションシステム、即ち糖尿病リスクシミュレーションシステムは、糖尿病患者の年齢、性別などの基本的な属性データ、及び血糖値、血圧などの検査値データを基にして、当該患者の網膜症、腎症、CHD(coronary heart disease;虚血性心疾患)及び死亡の、糖尿病に関わる主要なイベントの発生を長期間に渡りシミュレーションする。
【0020】
具体的にはシミュレーションにより、生存年数、一定期間におけるイベント(例えば、網膜症における失明や、腎症における透析)の数(ないしは発生率)、それらに伴う費用を推計し、適宜グラフやテーブルに示す。シミュレーションの基本原則は、各合併症のステージ間の推移と、CHD及び死亡の発生とを、確率的に予想するモンテカルロシミュレーションによるものである。
【0021】
(3)表示画面の推移・構成(画面フロー)
図2に示すように、本発明に係る糖尿病リスクシミュレーションシステムでは大きく分類すると、表示部4にて3類型の画面が表示される。つまり、
・初期データ入力画面100
・シミュレーション結果画面200
・詳細設定画面300
の3類型である。後述するように、夫々の類型は更に数種類の画面を含む。
【0022】
図2の画面の全体構成を簡単に説明する。初期データ入力画面100では、糖尿病患者の年齢、性別などの基本的な属性データ、及び血糖値、血圧などの検査値データが入力・設定される。そこでシミュレーションの実施が指示されると、シミュレーションの実施後、シミュレーション結果画面200が表示される。シミュレーション結果画面200は複数種類のシミュレーション結果画面から構成される。
【0023】
糖尿病患者の年齢、性別などの基本的な属性データ、及び血糖値、血圧などの検査値データの値が変更されるときには、通常、初期データ入力画面100へ戻される。ただし、シミュレーション結果画面200も、入力・設定値の変更を為し得る機能を併有する。よって、ある種類のシミュレーション結果画面200から別の種類のシミュレーショ結果画面200に直接移行したり、一つの種類のシミュレーション結果画面200で(入力・設定値を変更しつつ)シミュレーションを繰り返したりすることが可能である。
【0024】
詳細設定画面300では、将来における費用及びその効果に関する基本的データが設定される。初期データ入力画面100では、糖尿病患者個人に係るデータが入力設定されるのに対して、当該詳細設定画面300では患者に拠って変動しない費用及び効果に係るデータが設定・修正される。図2では基本的に初期データ入力画面100と詳細設定画面300との間で画面推移することを示しているが、シミュレーション結果画面200と詳細設定画面300との間で画面推移できるようにも構成されている。
【0025】
(4)ソフトウエア主要機能(4つのサブモデル)
本発明に係る糖尿病リスクシミュレーションシステムを構築するソフトウエアの主要な機能は、
・(糖尿病に関連する合併症である)網膜症の推移と腎症の推移を確率的に推計すること、
・(糖尿病患者でリスクが増大する)CHDの発生と死亡の発生を確率的に推計すること、
である。これらの推計のために、本発明に係る糖尿病リスクシミュレーションシステムは、
(i)網膜症の推移
(ii)腎症の推移
(iii)CHDの発生
(iv)死亡の発生
に関する4つのサブモデルを含む。
【0026】
上記の4つのサブモデルは、疫学データを利用する。利用する疫学データとしては、様々なものが存在し得る。以下の説明では、日本国における熊本スタディ、英国におけるUKPDS、米国におけるFramingham Heart Study及びHonolulu Heart Programを利用するものを挙げている。ここで、Honolulu Heart Programは日系人を対象とするものである(“Ralph B. D’Agostino etal. JAMA. 286. 180−187.2001”参照)。
【0027】
(i)網膜症の推移(網膜症サブモデル)
網膜症サブモデルは、
「網膜症なし」、
「単純網膜症」、
「前増殖網膜症」、
「増殖網膜症」、
「失明」
の5つのステージによって構成されるマルコフモデルである。これらのステージ間で一定の推移確率で、より重症度の高いステージへの推移が発生するとしており、逆戻りは設定していない。推移は、年次の推移を検討する。よって推移確率は年次の推移確率である。
【0028】
上記の推移確率は、熊本スタディのデータを基に設定されている。但し、「増殖網膜症」から「失明」への推移確率は、熊本スタディでは報告されていないため、米国における報告値(Scottら, 1988)が用いられている。
【0029】
「網膜症なし」から「増殖網膜症」までの各区間の推移確率は、熊本スタディのデータを基に、HbA1Cによるリスク調整の対象としている。即ち、HbA1C値が悪化する程、「網膜症なし→単純網膜症」、「単純網膜症→前増殖網膜症」、「前増殖網膜症→増殖網膜症」の各区間における網膜症の進展が早くなることになる。リスク調整については、後で説明する。
【0030】
更に、「単純網膜症→前増殖網膜症」への推移確率は、UKPDSのデータを基に収縮期血圧(SBP)によるリスク調整も行う。
【0031】
(ii)腎症の推移(腎症サブモデル)
腎症のサブモデルは、
「腎症なし」、
「微量アルブミン尿期」、
「顕性蛋白尿期」、
「腎不全保存期」、
「血液透析」
の5つのステージによって構成されるマルコフモデルである。網膜症サブモデルと同様、これらのステージ間で一定の確率で、より重症度の高いステージへの推移が発生するとしており、逆戻りは設定していない。推移はやはり年次の推移を検討するため、推移確率は年次の推移確率である。それら推移確率は、熊本スタディのデータを基に設定されている。
【0032】
「腎症なし」から「顕性蛋白尿期」までの各区間の推移確率は、熊本スタディのデータを基に、HbA1Cによるリスク調整の対象としている。更に、「腎症なし」から「顕性蛋白尿期」までの各区間の推移確率は、UKPDSのデータを基に、収縮期血圧(SBP)によるリスク調整も行う。
【0033】
リスク調整について説明する。上記(i)(ii)で示したリスク調整は、本発明では例えば、熊本スタディ(血糖コントロール)及びUKPDS(血圧コントロール)の結果を下記の式に代入して実施する。なお、下記の式は“Jonathan Betz Brown et al. Diabetes research and clinical practice.50(Suppl.3).s15−s46,2000”に示されるものである。
【数1】
RRi=RRt^(ln(vi/vm)/ln(vc/vt))
RRi:i番目にシミュレーションされる患者の調整係数
RRt:比較試験で治療の結果観察された平均相対危険
vi:i番目の患者のリスク因子の数値
vm:推移確率の基本値を計算した集団におけるリスク因子の平均値
vt:比較試験で治療群が達成したリスク因子の平均値
vc:比較試験でのコントロール群のリスク因子の平均値
vmについては熊本スタディにおける平均値(HbA1C:9.125、SBP:121)を使用する。vt、vcについては、HbA1Cによる調整時は7.2、9.4を、SBPによる調整時は144、154をそれぞれ使用する。相対リスク(RRt)は、下の表1と表2のように、熊本スタディにおける従来インスリン療法群と強化インスリン療法群の相対リスクを、またSBP修正時の相対リスク(RRt)はUKPDSにおける従来療法群と介入群との相対リスクを使用する。
【0034】
【表1】
【表2】
【0035】
例えば、HbA1C:8%のときの「網膜症なし」から「単純網膜症」への推移確率の調整は以下のようになる。
【数2】
0.052=3.72^(ln(8/9.125)/ln(9.4/7.2))
従って、推移確率は基本値の9.3%をかけて
【数3】
9.3%×0.052=0.48%
となる。この例のような「網膜症なし」から「単純網膜症」への場合はSBPによる調整は必要ないが、SBPによる調整が必要な場合にはHbA1Cによる調整のあとに更にSBPによる調整を行なう。
【0036】
(iii)CHDの発生(虚血性心疾患サブモデル)
CHD(虚血性心疾患)発生率の計算では、「Framingham Heart Study」の解析結果に基づくリスク計算式を基本としている。このリスク計算式では、年齢、性別、総コレステロール、HDLコレステロール、SBP、喫煙習慣の有無、耐糖能異常の有無及び左室肥大の有無が、リスクを計算するための因子とされる。
【0037】
但し、人種差を考慮し、日系人を対象とした「Honolulu HeartProgram」の結果を参考に補正を行っている。
【0038】
例えば、以下のAndersonの計算式により求められる発症確率pに一律「1/2.3」を掛けるという補正である。
【0039】
なお、「Framingham Heart Study」の解析結果に基づくリスク計算式については、いくつか提示されている。本発明では、例えば、Andersonらにより提示されている式を用いる。Andersonらの計算式は以下の通りである。
【0040】
{smoking:喫煙習慣ありの場合1、それ以外0}
{LVH:左室肥大ありの場合1、それ以外0}
{DM:耐糖能異常ありの場合1、それ以外0}
と定義する。ここで、本発明に係る糖尿病リスクシミュレーションシステムは、“耐糖能異常あり”の患者のデータを取り扱うことが前提となっている、として、“DM=1”と固定されていてもかまわない。
続いてまず、中間指標の「a」、「m」、「μ」、「σ」、「u」を算出する。
【0041】
aの計算;
【数4】
a=11.1122−0.9119*log(SBP)−0.2767*smoking−0.7181*log(TC/HDL)−0.5865*LVH
【0042】
mの計算;
男の場合
【数5】
m=a−1.4792*log(age)−0.1759*DM
女の場合
【数6】
m=a−5.8549+1.8515*[log(age/74)]^2−0.3758*DM
【0043】
μの計算;
【数7】
μ=4.4181+m
【0044】
σの計算;
【数8】
σ=exp(−0.3155−0.2784*m)
【0045】
uの計算;
【数9】
u=(log(t)−μ)/σ
【0046】
以上の5つの中間指標a、m、μ、σ、uによって、t年後のCHD発症確率pは下記の式で求められる。
【数10】
p=1−exp(−exp(u))
【0047】
(iv)死亡の発生(死亡サブモデル)
死亡率は、厚生労働省発表の簡易生命表における年齢・性別死亡率を基本値としている。それに加えて、「顕性蛋白尿期」では自然死亡率の1.5倍、「腎不全保存期」は2倍とする。更に、「血液透析」では一律14.1%/年の死亡率とし、虚血性心疾患発症後は16.2%/年の死亡率とする。
【0048】
上記の簡易生命表以外の倍率は、“日本内科学会雑誌 87(7):81−86、1998”に基づく。
【0049】
(5)モンテカルロシミュレーションフローチャート
図1にて、本発明の糖尿病リスクシミュレーションシステムに係るモンテカルロシミュレーションの概略のフローチャートを示す。シミュレーション回数(n)は操作者が自由に設定し得る。図1では「n=10000」としている。図1により、モンテカルロシミュレーションの概略の処理の流れを説明する。
【0050】
(ステップS00):シミュレーションを開始する。
【0051】
(ステップS02):n回目の(n番目の)シミュレーションを開始する。まず必要パラメータを初期化する。例えば、
死亡フラグ=0
網膜症ステージ=R1(「網膜症なし」)
腎症ステージ=N1(「腎症なし」)
CHDフラグ=0(「無し」)
等である。
【0052】
(ステップS04):網膜症サブモデル計算を行う。網膜症ステージ毎に設定されている推移確率trと、プログラムで発生させる乱数ranとを比較し、
tr>ran
ならば網膜症ステージをひとつ進展させる。
【0053】
(ステップS06):腎症サブモデル計算を行う。腎症ステージ毎に設定されている推移確率trと、プログラムで発生させる乱数ranを比較し、
tr>ran
ならば腎症ステージをひとつ進展させる。
【0054】
(ステップS08):CHDサブモデル計算を行う。フラミンガムリスク推計式により、CHD発生率を計算する。但し前述のように、「Honolulu Heart Program」の結果による補正をかけている。発生があれば、CHDフラグを“1”にする。更に、CHDフラグ=1のとき、すなわちCHD発生時には、死亡率「16.2%」と、プログラムで発生した乱数ranを比較し、
16.2>ran
なら死亡フラグを1にする。
【0055】
(ステップS10):死亡サブモデル計算を行う。各年の死亡率は簡易生命表による年齢・性別死亡率(以下、“自然死亡率”という。)を基本値として、腎症サブモデルで計算された腎症ステージに応じて死亡率trを以下の表3のとおり求める。
【表3】
死亡率trと、プログラムで発生した乱数ranを比較し、
tr>ran
なら死亡フラグを1にする。
【0056】
(ステップS12):死亡判定を行う。死亡フラグが1なら死亡と判定し、ステップS14に進む。死亡フラグが0ならば、ステップS04に戻り、上記の4つのサブモデルを繰り返す。
【0057】
即ち、n番目のシミュレーションについては、死亡フラグが1となるまで繰り返されることになる。
【0058】
(ステップS14):シミュレーション回数を判定する。シミュレーション回数が設定値(n=10000)に達したら、シミュレーションを終了するためステップS18に進む。シミュレーション回数が設定値に未達であれば、シミュレーション回数をインクリメントして(ステップS16)、次回数のシミュレーションを行う(ステップS02)。
【0059】
(ステップS18):シミュレーションを終了する。
【0060】
(6)画面について
次に、表示画面について説明する。
【0061】
図4は、初期データ入力画面100である。この画面では、糖尿病患者の年齢、性別などの基本的な属性データ、血糖値、血圧、コレステロールなどの検査値データ、及び喫煙習慣・左室肥大などのデータが入力・設定される。糖尿病の合併症のうち、網膜症と腎症の状態についても入力する。シミュレーション開始ボタン108でシミュレーションを実施する。また、詳細設定ボタンで、医療費などのパラメータ設定を行なう詳細設定画面300に進む。
【0062】
図5は、シミュレーション結果画面200のうちの網膜症シミュレーション画面202である。初期データ入力画面で入力した値(現在値)によってシミュレーションされた今後30年間にわたる網膜症の進展予測が表示される。ここで、「腎症」ボタン212を選択すれば、腎症進展についての同様のグラフ(腎症シミュレーション画面(図示せず。))が表示される。
【0063】
上記の網膜症シミュレーション画面202は、「網膜症」ボタン210を表示する画面にて該「網膜症」ボタン210を選択することによっても表示される。同様に、上記の腎症シミュレーション画面(図示せず。)も、「腎症」ボタン212を表示する画面にて該「腎症」ボタン212を選択することによっても表示される。
【0064】
図6は、シミュレーション結果画面200のうちのイベント発生率シミュレーション画面222である。各画面の「イベント発生率」ボタン214が選択されると、初期データ入力画面100で入力された値(現在値)によってシミュレーションされた生涯にわたり発生が予測される各種イベント発生率が表示される(左部226)。
【0065】
図7は、シミュレーション結果画面200のうちの医療費シミュレーション画面232である。各画面の「医療費」ボタン216が選択されると、初期データ入力画面で入力した値(現在値)によってシミュレーションされた生涯にわたり発生が予測される糖尿病合併症に係る医療費が表示される。この医療費の計算において基礎のデータとなるのは、
・初期データ入力画面100で入力された値(現在値)、
・図1に示すシミュレーションフローチャートにより計算され、図5の網膜症シミュレーション画面202、図6のイベント発生率シミュレーション画面222、及び腎症シミュレーション画面(図示せず。)などに表示される確率(発生率)、
・詳細設定画面300で入力される合併症別のステージ毎の年間費用及び合併症別のイベント発生時の年間費用、
である。
【0066】
図8は、シミュレーション結果画面200のうちのシミュレーション一覧画面242である。各画面の「一覧」ボタン218が選択されると、初期データ入力画面100で入力された値(現在値)によってシミュレーションされたイベント発生率、費用対効果に関する情報が一覧表示される。費用対効果については、後で説明する。
【0067】
図9も、シミュレーション結果画面200のうちの網膜症シミュレーション画面202である。ところで、「仮想値」の入力設定欄には、現在値と異なるデータが入力・設定され得る。そうすると、現在値によるシミュレーション結果(画面左部)を表示したまま、仮想値シミュレーションを実施することができる。図9の表示例では、HbA1c、SBP、TCを現在値よりも少なくしてシミュレーションを実施しようとしている。「シミュレーション開始」ボタンで仮想値によるシミュレーションが実施される。
【0068】
図10も、シミュレーション結果画面200のうちの網膜症シミュレーション画面202であるが、入力された仮想値によってシミュレーションされた30年間にわたる網膜症の進展予測が、すでに実施された現在値によるシミュレーション結果の隣に表示される。ここで、「腎症」ボタン212を選択すれば、腎症進展についての同様のグラフ(腎症シミュレーション画面(図示せず。))が表示される。
【0069】
図11も、図6と同様のシミュレーション結果画面200のうちのイベント発生率シミュレーション画面222である。但し図6と異なり、仮想値によりシミュレーションされた生涯にわたる各種イベント発生率が表示されている。このように、「仮想値」の入力設定欄に現在値と異なるデータが入力・設定された状態で各画面の「イベント発生率」214が選択されると、現在値シミュレーション結果の隣に仮想値によりシミュレーションされた生涯にわたる各種イベント発生率が表示される。
【0070】
図12も、図7と同様のシミュレーション結果画面200のうちの医療費シミュレーション画面232である。但し図7と異なり、仮想値によりシミュレーションされた生涯にわたる各種医療費が表示されている。このように、「仮想値」の入力設定欄に現在値と異なるデータが入力・設定された状態で各画面の「医療費」ボタン216が選択されると、現在値シミュレーション結果の隣に仮想値によりシミュレーションされた生涯にわたる各種医療費が表示される。
【0071】
図13も、図8と同様のシミュレーション結果画面200のうちのシミュレーション一覧画面242である。但し図8と異なり、特に仮想値を利用してシミュレーションされたイベント発生率、費用対効果に関する情報が一覧表示されている。このように、「仮想値」の入力設定欄に現在値と異なるデータが入力・設定された状態で各画面の「一覧」ボタン218が選択されると、現在値及び仮想値によってシミュレーションされたイベント発生率、費用対効果に関する情報が一覧表示される。
【0072】
図13において、ARR(%)とRRR(%)とは、以下の式で定義される。
ARR(%)=仮想値(%)−現在値(%)
RRR(%)=100×{仮想値(%)−現在値(%)}/現在値(%)
また、「生存年」は、シミュレーションにおいて計算される余命年の平均値である。「質調整生存年」は、上記「生存年」において、糖尿病合併症別の各ステージ推移やイベント発生が生じたときに、図14に示される「効用値」を生存年として累積したものである。
【0073】
更に図13において、費用対効果(生存年)と費用対効果(質調整生存年)とは、以下の式で定義される。
費用対効果(生存年)=総医療費/生存年
費用対効果(質調整生存年)=総医療費/質調整生存年
増分費用対効果(生存年)と増分費用対効果(質調整生存年)とは、以下の式で定義される。
増分費用対効果(生存年)=増分総医療費/増分生存年
増分費用対効果(質調整生存年)=増分総医療費/増分質調整生存年
ここで、「増分総医療費」は、図13に示す通り、仮想値の総医療費から現在値の総医療費を引いた値である。
【0074】
図14は、詳細設定画面300である。各画面の「詳細設定」ボタンが選択されると、表示・設定される。シミュレーションで使用する費用パラメータ(304、306、308)、効用値310、又はシミュレーションにおけるシミュレーション回数312を変更することができる。更に、割引率設定ボタン314により、割引率を設定する画面(図示せず。)に移行することもできる。
【0075】
ここで「割引率」とは、費用と質調整生存年とを、一定の年率で現在価値に割り戻すための数値である。
【0076】
詳細設定画面300では、費用パラメータとして合併症基本診療費304と追加診療費とを区別して入力・設定可能なように構成されている。更に、追加診療費として現在値306と仮想値308とを区別して入力設定可能なように構成されている。
【0077】
例えば、糖尿病合併症の各ステージでの基本的な医療費を「合併症基本診療費」に設定し、加えて、初期データ入力画面100での現在値を維持したりシミュレーション結果画面200での仮想値を達成したりするために必要な診療費を「追加診療費・現在値」「追加診療費・仮想値」に適宜設定する、というようなことができる。
【0078】
以上、本発明の具体的な実施の形態を説明したが、本発明はこれらに限られず更に種々改変可能である。
【0079】
上記の本発明に係る実施の形態では、糖尿病合併症に関するリスクシミュレーションのためのシステムを取り上げてきた。このシステムは、その他の疾患・疾病に関するリスクシミュレーションのためにも応用可能である。
【0080】
既述のように、本発明の実施の形態は、計算機システム2が外部ネットワーク(図示せず。)を介して外部サーバと接続し、該外部サーバに用意される本発明に係る疾患リスクシミュレーションの処理プログラムを利用して(ダウンロードして)、本発明に係る疾患リスクシミュレーションを行なう、というようなものであってもよい。そうすると、該外部サーバはシミュレーション手段を提供するとともに、フィードバックデータを大規模に収集することも可能になる。
【0081】
【発明の効果】
本発明を利用することにより、以下のような効果を得ることができる。
【0082】
様々な疫学的調査・解析を基にして、多数年間(数十年間)に渡る累積的なシミュレーションを行なうことができる。その際、人種差を考慮したシミュレーションを実施できる。
【0083】
また、費用のパラメータの設定、及び費用対効果の算出などにより、経済的効用の視点からのシミュレーションを可能にする。
【図面の簡単な説明】
【図1】糖尿病リスクシミュレーションシステムのモンテカルロシミュレーションの概略のフローチャートである。
【図2】糖尿病リスクシミュレーションシステムの画面の分類である。
【図3】疾患リスクシミュレーションシステム(装置)として稼動し得る計算機システムの1つの形態である。
【図4】初期データ入力画面である。
【図5】シミュレーション結果画面のうちの網膜症シミュレーション画面である。
【図6】シミュレーション結果画面のうちのイベント発生率シミュレーション画面である。
【図7】シミュレーション結果画面のうちの医療費シミュレーション画面である。
【図8】シミュレーション結果画面のうちのシミュレーション一覧画面である。
【図9】「仮想値」の入力設定欄に仮想値が入力された、シミュレーション結果画面のうちの網膜症シミュレーション画面である。
【図10】仮想値によってシミュレーションされた30年間にわたる網膜症の進展予測を伴う、シミュレーション結果画面のうちの網膜症シミュレーション画面である。
【図11】仮想値によりシミュレーションされた生涯にわたる各種イベント発生率の表示を伴う、シミュレーション結果画面のうちのイベント発生率シミュレーション画面である。
【図12】仮想値によりシミュレーションされた生涯にわたる各種医療費の表示を伴う、シミュレーション結果画面のうちの医療費シミュレーション画面である。
【図13】現在値及び仮想値によってシミュレーションされたイベント発生率、費用対効果に関する情報が一覧表示された、シミュレーション結果画面のうちのシミュレーション一覧画面である。
【図14】詳細設定画面である。
【符号の説明】
100・・・初期データ入力画面、200・・・シミュレーション結果画面、300・・・詳細設定画面。
Claims (7)
- 所定の属性値及び検査値を入力する要素データ入力手段と、
所定の疾病に係る複数のステージ間の推移と、所定のイベント発生とを確率的予測手法により予測し、一定期間のイベント、及び生存年数をシミュレーションするシミュレーション演算手段と、
上記のシミュレーション演算手段の作出する演算結果を、夫々に対応して出力するシミュレーション結果出力手段とを含む疾病リスクシミュレーションシステム。 - シミュレーション演算手段において、ステージ間推移の推移確率として、標本対象の人種が限定された統計に基づく確率値が利用される、
請求項1に記載の疾病リスクシミュレーションシステム。 - シミュレーション演算手段において、所定のイベント発生のイベント発生率が、標本対象の人種が限定された統計により補正される、
請求項1又は請求項2に記載の疾病リスクシミュレーションシステム。 - シミュレーション演算手段において、死亡率を他の疾病に係るステージの状態に依存させる、
請求項1乃至請求項3のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステム。 - シミュレーション演算手段において、所定の疾病に係る複数のステージ間の推移と、所定のイベント発生とに係る費用がシミュレーションされ、
将来における費用に関する基本的データとその効果に関する基本的データを設定する詳細条件設定入力手段を、更に含む、
請求項1乃至請求項4のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステム。 - 上記の詳細条件設定入力手段において、夫々の疾病のイベント又は夫々の疾病の個別のステージに関し、
費用に関する基本的データとして、少なくとも、基本値、追加値及び追加仮想値が設定され、更に、効果に関する基本的データとして、少なくとも、効用値が設定され、
それら設定値がシミュレーション演算手段での将来における費用対効果に係る演算で利用され、
シミュレーション結果出力手段に出力される、
請求項5に記載の疾病リスクシミュレーションシステム。 - 所定の疾病に係る複数のステージが、マルコフモデルに従い形成されており、
確率的統計手法として、モンテカルロシミュレーションの手法が利用される、請求項1乃至請求項6のうちのいずれか一つの請求項に記載の疾病リスクシミュレーションシステム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003180975A JP2005052164A (ja) | 2003-06-11 | 2003-06-25 | 疾病リスクシミュレーションシステム |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003166794 | 2003-06-11 | ||
JP2003180975A JP2005052164A (ja) | 2003-06-11 | 2003-06-25 | 疾病リスクシミュレーションシステム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2005052164A true JP2005052164A (ja) | 2005-03-03 |
Family
ID=34379701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003180975A Withdrawn JP2005052164A (ja) | 2003-06-11 | 2003-06-25 | 疾病リスクシミュレーションシステム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2005052164A (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013061785A (ja) * | 2011-09-13 | 2013-04-04 | Kurume Univ | 生死予測装置とその方法、生死予測プログラム及び記録媒体 |
JP2018005726A (ja) * | 2016-07-06 | 2018-01-11 | オムロンヘルスケア株式会社 | リスク分析システム及びリスク分析方法 |
JP2018055728A (ja) * | 2016-07-06 | 2018-04-05 | オムロンヘルスケア株式会社 | リスク分析システム及びリスク分析方法 |
JP2018077896A (ja) * | 2018-01-10 | 2018-05-17 | オムロンヘルスケア株式会社 | リスク分析システム及びリスク分析方法 |
JP2018088263A (ja) * | 2018-01-10 | 2018-06-07 | オムロンヘルスケア株式会社 | 血圧リスク分析装置及び血圧リスク分析方法 |
-
2003
- 2003-06-25 JP JP2003180975A patent/JP2005052164A/ja not_active Withdrawn
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013061785A (ja) * | 2011-09-13 | 2013-04-04 | Kurume Univ | 生死予測装置とその方法、生死予測プログラム及び記録媒体 |
JP2018005726A (ja) * | 2016-07-06 | 2018-01-11 | オムロンヘルスケア株式会社 | リスク分析システム及びリスク分析方法 |
JP2018055728A (ja) * | 2016-07-06 | 2018-04-05 | オムロンヘルスケア株式会社 | リスク分析システム及びリスク分析方法 |
US11087406B2 (en) | 2016-07-06 | 2021-08-10 | Omron Healthcare Co., Ltd. | Risk analysis system and risk analysis method |
JP2018077896A (ja) * | 2018-01-10 | 2018-05-17 | オムロンヘルスケア株式会社 | リスク分析システム及びリスク分析方法 |
JP2018088263A (ja) * | 2018-01-10 | 2018-06-07 | オムロンヘルスケア株式会社 | 血圧リスク分析装置及び血圧リスク分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7487134B2 (en) | Medical risk stratifying method and system | |
JP6075973B2 (ja) | 健康状態判定装置およびその作動方法 | |
US20090177453A1 (en) | Medical diagnosis support computer system, computer program, and server computer | |
US20060265136A1 (en) | Treatment effect prediction system, a treatment effect prediction method, and a computer program product thereof | |
JP6282783B2 (ja) | 分析システム及び分析方法 | |
JP2020537232A (ja) | 全母集団から任意に選択された部分母集団内の被験者における非健康状態のリスク、発生又は進行を予測する医療装置及びコンピュータ実装方法 | |
JP2007265347A (ja) | 健康指導支援装置、システム、及びプログラム | |
WO2020027213A1 (ja) | 認知症リスクの提示システムおよび方法 | |
Kroll | Biological variation of glucose and insulin includes a deterministic chaotic component | |
Zhao et al. | Interindividual glucose dynamics in different frequency bands for online prediction of subcutaneous glucose concentration in type 1 diabetic subjects | |
Yau et al. | Reinforcement learning models and algorithms for diabetes management | |
EP3718116B1 (en) | Apparatus for patient data availability analysis | |
Van Herpe et al. | An adaptive input–output modeling approach for predicting the glycemia of critically ill patients | |
Cescon et al. | Subspace-based linear multi-step predictors in type 1 diabetes mellitus | |
JP2005052164A (ja) | 疾病リスクシミュレーションシステム | |
JP7238705B2 (ja) | 診療支援方法、診療支援システム、学習モデルの生成方法、および、診療支援プログラム | |
US9202009B2 (en) | Diagnostic support apparatus for diabetes and computer program product | |
Enright et al. | Modelling and monitoring the individual patient in real time | |
Gruber | Targeted learning in healthcare research | |
Horvat et al. | To improve a prediction model, give it time | |
Hettiarachchi et al. | G2P2C—A modular reinforcement learning algorithm for glucose control by glucose prediction and planning in Type 1 Diabetes | |
US20210043327A1 (en) | Computing device, portable device and computer-implemented method for predicting major adverse cardiovascular events | |
JP7479317B2 (ja) | 情報処理装置、方法及びプログラム | |
US20220068485A1 (en) | Computing system that generates patient-specific outcome predictions | |
Sun et al. | Package ‘CopulaCenR’ |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Withdrawal of application because of no request for examination |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20060905 |