(57)【要約】  医療処置の効果を評価するコンピュータ実行可能なシステムおよび方法であって、その方法は、患者記録データを受信し、評価のための関連する特徴を識別し、第１の処置を識別し、第２の処置を識別し、各患者の症例に重みを割り振り、割り振られた重みを用いて、第２の識別される処置と対比して、識別される処置が識別される効果をもたらす相対的な可能性を決定し、この推定された相対的な可能性を出力することを含む。【選択図】図１
2012530300.tif 000002
【特許請求の範囲】
【請求項１】
医療処置の効果を評価するコンピュータ実行可能な方法であって、
患者記録データを受信し、
評価のための関連する特徴を識別し、
第１の処置を識別し、
第２の処置を識別し、
各患者の症例に重みを割り振り、
割り振られた前記重みを用いて、第２の識別される処置と対比して、識別される処置が識別される効果をもたらす相対的な可能性を決定し、
この推定された相対的な可能性を出力する
ことを含む方法。
【請求項２】
各患者の症例に割り振られた前記重みが、前記患者が被験グループのメンバーであろう可能性に基づいている、請求項１に記載の方法。
【請求項３】
各患者の症例に割り振られた前記重みが、前記患者の特色および特色の組み合わせを目標の母集団に整合させることによって計算される、請求項２に記載の方法。
【請求項４】
前記相対的な可能性を決定する際に回帰が行われる、請求項３に記載の方法。
【請求項５】
ポアソン回帰が行われる、請求項４に記載の方法。
【請求項６】
整合させた前記特色のうちの少なくともいくつかが前記回帰で共変量として用いられる、請求項５に記載の方法。
【請求項７】
医療処置の効果を評価するためのシステムであって、
ネットワーク・インタフェースと、
ネットワークを介してアクセスされるサーバ上にある患者記録データベースと、
患者記録データを受信し、評価のための関連する特徴を識別し、第１の処置を識別し、第２の処置を識別し、前記患者が被験グループのメンバーであろう可能性に基づいて各患者の症例に重みを割り振り、割り振られた前記重みを用いて、第２の識別される処置と対比して、識別される処置が識別される効果をもたらす相対的な可能性を決定し、この推定された相対的な可能性を出力するように構成されている解析サーバと
を備えるシステム。
【請求項８】
医療処置の効果を評価するための解析サーバであって、
プロセッサと、
前記プロセッサに動作可能に連結され、内部にプログラミング命令を蓄積するメモリであって、前記プロセッサがプログラム命令を実行するように動作可能であるメモリとを備え、
前記プログラム命令が、
患者記録データを受信し、
評価のための関連する特徴を識別し、
第１の処置を識別し、
第２の処置を識別し、
前記患者が被験グループのメンバーであろう可能性に基づいて各患者の症例に重みを割り振り、
割り振られた前記重みを用いて、第２の識別される処置と対比して、識別される処置が識別される効果をもたらす相対的な可能性を決定し、
この推定された相対的な可能性を出力する
ことを含む解析サーバ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般に、後向きの医療データを統計的に解析するための方法およびシステムに関し、より詳細には、医療処置の効果を評価するための方法およびシステムに関する。
【背景技術】
【０００２】
コホート研究を用いて医療処置の効果を解析するためのシステムおよび方法がある。これらのシステムおよび方法のほとんどは、特定の処置のランダム化比較試験に基づいている。十分な期間にわたって十分な数の被験者が与えられれば、ランダム化比較試験には、その処置を受けるグループと制御グループの全体に交絡因子（またはグループ間の差異）を一様に分布させることによって、この解析が簡略化されるという利点がある。ランダム化比較試験は、実行するのに費用がかかり、それらの精度はある程度は多数の被験者に依存し、処置の候補の間で典型的でない条件に起因する副次的な交互作用を識別するのにあまり効果的ではないため、これらのシステムおよび方法は制約を受ける。その上、ランダム化比較試験（これは前向きのコホート研究である）では、既存の患者記録データのより広範な集合で利用可能な多量の情報を、後向きのコホート研究で活用できるのと同じようには活かすことができない。
【０００３】
他の既存のシステムおよび方法は、履歴の患者記録データを解析して特定の処置の効果を評価する。これらのシステムは、処置されたグループと未処置グループの全体にわたって履歴データの中に一様に分布していない交絡因子や患者の特徴に関する制御のために主として手操作によるモデリング処理を用いるため、これらのシステムも制約を受ける。これらの手操作による交絡因子の制御処理では偏りが生じる危険性が高いため、既存のシステムによって生み出される結果の質は、オペレータの熟練度に大きく依存する。その上、大きな既存のデータ集合中の多数の潜在的な交絡因子に付随する主要な効果や交互作用の効果により、手操作のモデリングはごく短期間に集中したものになり、人的ミスに影響されやすくなる。
【発明の概要】
【０００４】
本発明は、識別される処置を用いて治療されたときに患者が識別される効果を経験する相対的な可能性を評価する目的で患者の医療記録を統計的に解析するためのソフトウェア・コンポーネントのシステムおよび方法に関する。この解析方法により、どの処置を承認し、使用し、資金提供するかを決定する際に、識別される処置の効力や危険性を評価するために医療分野の意志決定者が依拠することができる結果が生み出される。
【０００５】
いくつかの実施形態では、そのシステムは、２つの処置を評価して、第１の処置を受ける患者が、第２の処置を受ける患者と比較して識別される効果を経験する相対的な可能性を決定する。他の実施形態では、比較される複数の処置は、実際には同じ処置の複数の適用量である。さらに他の実施形態では、そのシステムは、いくつかの個々の比較を行っていくつかの異なる比率比を生成することにより第１の処置をいくつかの異なる第２の処置と比較する。その比率比のそれぞれは、第１の処置と別の処置との間のものとして識別される効果を患者が経験する相対的な可能性を表す。
【０００６】
本発明の一実施形態により、医療処置の効果を評価するコンピュータ実行可能なシステムおよび方法が提供され、その方法は、患者記録データを受信し、評価のための関連する特徴を識別し、第１の処置を識別し、第２の処置を識別し、各患者の症例に重みを割り振り、割り振られた重みを用いて、第２の識別される処置と対比して、識別される処置が識別される効果をもたらす相対的な可能性を決定し、この推定された相対的な可能性を出力することを含む。
【０００７】
別の実施形態では、本発明は、医療処置の効果を評価するためのシステムに関し、そのシステムは、ネットワーク・インタフェースと、ネットワークを介してアクセスされるサーバ上にある患者記録データベースと、患者記録データを受信し、評価のための関連する特徴を識別し、第１の処置を識別し、第２の処置を識別し、患者が被験グループのメンバーであろう可能性に基づいて各患者の症例に重みを割り振り、割り振られた重みを用いて、第２の識別される処置と対比して、識別される処置が識別される効果をもたらす相対的な可能性を決定し、この推定された相対的な可能性を出力するように構成されている解析サーバとを備える。
【０００８】
別の実施形態では、本発明は、医療処置の効果を評価するためのサーバに関し、そのサーバは、プロセッサと、プロセッサに動作可能に連結され、内部にプログラミング命令を蓄積するメモリであって、プロセッサがプログラム命令を実行するように動作可能であるメモリとを備え、プログラム命令が、患者が被験グループのメンバーであろう可能性に基づいて各患者の症例に重みを割り振り、第２の識別される処置と対比して、識別される処置が識別される効果をもたらす相対的な可能性を決定する際に、割り振られた重みを用いる。
【図面の簡単な説明】
【０００９】
【図１】後向き観察のデータに基づいて医療処置の効果を比較するための、本発明の実施形態により構成された医療処置仮説検定システムの概略ブロック図である。
【図２】本発明の実施形態による図１の医療処置仮説検定システムのデータ選択コンポーネントによって実行されるステップの高レベルのフロー図である。
【図３】本発明の実施形態による図１の医療処置仮説検定システムの傾向スコアリング・コンポーネントによって実行されるステップの高レベルのフロー図である。
【図４】本発明の実施形態による図３の傾向スコアリング・コンポーネントの、テーブルのデータを準備するコンポーネントによって実行されるステップのフロー図である。
【図５Ａ】本発明の実施形態によるテーブルのデータを準備する図４の処理の最初における例示的な患者症例データのテーブルの説明図である。
【図５Ｂ】本発明の実施形態によるテーブルのデータを準備する図４の処理のある時点における例示的な患者症例データのテーブルの説明図である。
【図５Ｃ】本発明の実施形態によるテーブルのデータを準備する図４の処理のある時点における例示的な患者症例データのテーブルの説明図である。
【図５Ｄ】本発明の実施形態によるテーブルのデータを準備する図４の処理のある時点における例示的な患者症例データのテーブルの説明図である。
【図５Ｅ】本発明の実施形態によるテーブルのデータを準備する図４の処理のある時点における例示的な患者症例データのテーブルの説明図である。
【図６】本発明の実施形態による、処置グループ内の患者についての重み付けされていない患者データと処置グループ内にない患者についての重み付けされた患者データとを比較するための例示的なテーブルの説明図である。
【図７】図１の回帰モデリング１１８用の例示的な患者症例データのテーブルの説明図であり、この説明図は図３で回帰ステップ用にデータをフォーマット３１２する間にフォーマットされたデータを示す。
【発明を実施するための形態】
【００１０】
図面は、あくまで説明のために本発明の実施形態を描いている。本明細書に記載された本発明の原理から逸脱することなく本明細書で説明される構造および方法の代替実施形態を採用できることが、当業者は以下の記述から容易に理解されよう。
【００１１】
さて、図面に目を向けると、患者データ選択、傾向スコアリング、逆確率重み付け、および比率比の２重にロバストな推定を含む、本発明の実施形態が示されている。調査対象となる２つの識別される処置および識別される効果に基づいて、後向き観察の医療記録データがオペレータにより患者記録データベースから選択される。第１の識別される処置を受けなかった各患者には、その患者の個々の特徴に基づいて、第１の識別される処置を受けたグループにその患者がいたであろう可能性を表す傾向スコアが割り振られる。その割り振られた傾向スコアを用いて各患者のデータを重み付けし、第１の識別される処置を受けなかったグループ（即ち、第２の識別される処置を受けたグループ）の重み付けされたデータが第１の処置を受けたグループの重み付けされていないデータに近似するようにする。重み付けされたデータを用いて回帰を行って、第２の識別される処置を受ける患者と比べて第１の識別される処置を受ける患者が識別される効果を経験するであろう相対的な可能性を推定する。この推定値のことを比率比と呼ぶ。回帰を行う際に重み付けされたデータを用いることによって、比率比の推定は２重にロバストになる。
【００１２】
図１は、後向き観察のデータに基づいて医療処置の効果を比較するための、本発明の実施形態により構成された医療処置仮説検定システムの概略ブロック図である。図１に示すように、そのシステムは、オペレータがネットワーク１０４を介して解析マシン１０６と患者データ・マシン１０８にアクセスするためのオペレータ端末１０２を備える。患者データ・マシン１０８により、オペレータが患者記録データベース１１０と補助データベース１１２にアクセスできる。解析マシン１０６は、プロセッサと、プロセッサに動作可能に連結され、内部にプログラミング命令や他のデータを蓄積するメモリであって、プロセッサがプログラム命令を実行するように動作可能であるメモリと、解析マシンが患者データ・マシン１０８や他の入力元からの入力を受信し、オペレータ端末１０２や他の出力先に結果を出力することを可能にするネットワーク接続とを備えてもよい。
【００１３】
患者記録データベース１１０は、患者の識別番号、誕生日、性別、医療検査の結果といった患者の医療記録からの後向き観察の患者データ、医療提供者により記録された観察データ、および患者により医療提供者に提供された情報を蓄積する。補助データベース１１２は、どの記録が患者記録データベースから検索されたかやそれらがいつ検索されたかについての情報などの、システムの運用に関する情報、システム・データのフォーマット・ルール、および患者記録データの解析に関係する他のデータを蓄積する。
【００１４】
解析マシン１０６により、データ選択コンポーネント１１４、傾向スコアリング・コンポーネント１１６、回帰モデリング・コンポーネント１１８およびデータ出力コンポーネント１２０を含むいくつかのソフトウェア・コンポーネントにオペレータはアクセスできる。これらのソフトウェア・コンポーネントは、コンピュータ・メモリに蓄積されコンピュータ・プロセッサにより実行されるコンピュータ命令の形をとってもよい。データ選択コンポーネント１１４は、関連する患者のデータや属性を患者記録データベース１１０から選択するためのインタフェースをオペレータに提示し、選択されたデータを検索しかつ傾向スコアリング・コンポーネント１１６用にフォーマットする。
【００１５】
傾向スコアリング・コンポーネント１１６は、識別される処置を受ける患者のグループ内に（「処置グループ内に」）その患者がいる可能性を表す傾向スコアを決定し、その傾向スコアを各患者の記録に割り振る。さらに、傾向スコアリング・コンポーネント１１６は、傾向スコアを患者データに適用して患者記録のデータを重み付けし、処置グループ内にない患者のグループについての重み付けされたデータが処置グループ内の患者のグループについての重み付けされていないデータに近似するようにする。
【００１６】
回帰モデリング・コンポーネント１１８により、第２の識別される処置を受ける患者（即ち、処置グループ内にない患者）に比べて第１の識別される処置を受ける患者が識別される効果を経験するであろう相対的な可能性を推定するためのモデルをオペレータが構築し検査するためのインタフェースが与えられる。回帰モデリング・コンポーネント１１８は、傾向スコアリング・コンポーネント１１６によって重み付けされた重み付けデータを受け取る。
【００１７】
データ出力コンポーネント１２０により、解析マシンの結果を提示するためのフォーマットやスタイルをオペレータが選択することができるインタフェースが与えられる。いくつかの実施形態では、データ出力コンポーネント１２０は、解析マシンによって生成されたデータを選択しフォーマットするためのツールを含む。他の実施形態では、データ出力コンポーネント１２０内のツールにより、解析マシンの結果の解釈や理解を助けるための図表やグラフなどの色々な可視化ツールを、オペレータは選択したり操作することができる。
【００１８】
図２は、本発明の実施形態による図１の医療処置仮説検定システムのデータ選択コンポーネント１１４によって実行されるステップの高レベルのフロー図である。図２に示すように、データ選択コンポーネント１１４は、最初に患者記録データベースへの接続を確立２０２して、オペレータが患者データにアクセスできるようにする。本発明の一実施形態によれば、データ選択コンポーネント１１４は、ＯＤＢＣドライバを利用して、患者データ・マシン１０８上の患者記録データベース１１０にネットワーク１０４を介してデータ選択コンポーネント１１４を接続する。次いで、データ選択コンポーネント１１４はデータ選択インタフェースをオペレータに提供２０４して、解析マシン上での後続する解析用の適切な患者記録をオペレータが選択できるようにする。本発明の一実施形態によれば、データ選択インタフェースは、ウェブ・ブラウザ・アプリケーションで閲覧されるＨＴＭＬウェブページとして、オペレータに提示される。オペレータが所望のデータを選択した後に、システムはこの選択を受信し２０６、選択された記録を患者記録データベース１１０から検索２０８する。次いで、システムは、解析マシン用のデータを傾向スコアリング・コンポーネントによってフォーマット２１０し、フォーマットされたデータを解析マシン１０６上のメモリに蓄積２１２する。
【００１９】
図３は、本発明の実施形態による図１の医療処置仮説検定システムの傾向スコアリング・コンポーネント１１６によって実行されるステップの高レベルのフロー図である。図３に描いたように、傾向スコアリング・コンポーネント１１６は初期データ準備ステップ３０２を実行するが、このステップは図４に関してより詳細に説明する。このデータ準備ステップ３０２では、解析マシンは、患者記録をテーブルの形で準備し、「初期値の」（ｉｎｉｔｉａｌ  ”ｓｅｅｄｉｎｇ”）傾向スコアと各患者記録に付随する残差を識別し、選択されたデータ内に記載されている各患者の特徴について、各特性関数や特性関数の積のための列を作成し、それらを評価する。データ準備ステップ３０２の後、解析マシンは、いずれかの特性関数の列または特性関数の積の列と残差の列との間で最大の絶対相関を識別３０４する。以下の公式について、ｔｉは第１の識別される処置を受ける患者ｉの０／１指標を表し、ｐｉは患者ｉが第１の識別される処置を受けた推定確率であり、Ｉｊはｊ番目の特性関数であり、ｎは患者の人数である。最大の絶対相関を識別するために、いずれかの２つの列がどの程度まで相関しているかを決定するには、以下の公式を使用する：ｒｊ＝Ｓｕｍｉ（（ｔｉ−ｐｉ）（Ｉｊｉ−ｍｅａｎ（Ｉｊ））／（（ｎ−１）ｓｄ（ｔ−ｐ）ｓｄ（Ｉｊ））。残差の列に最も相関している列を識別した後で、解析マシンは、識別された列に基づいてすべての患者についての傾向スコアを調整３０６する。傾向スコアを調整する際には、以下の公式を使用する：ｐｉ／（１−ｐｉ）＝ｐｉ／（１−ｐｉ）×ｅｘｐ（（δ×ｓｉｇｎ（ｒｊ２）×（Ｉｊｉ−１））。ここで、δは０．００１などの小さな数に設定されるチューニング・パラメータである。次いで、調整処理の結果得られる新たな傾向スコアを用いて、解析マシンは、未処置の患者記録のそれぞれに適用する重みを決定３０８する。その重みは、以下の公式を用いて傾向スコアから計算される：ｗｉ＝ｐｉ／（１−ｐｉ）。
【００２０】
解析マシンは、処置グループ内にない患者についての重み付けされたデータの総計を、処置グループ内の患者についてのデータの総計と比較３０９して、２つのデータ集合が十分に類似している（「最適に釣り合っている」）かどうかを決定３１０する。それらのデータ集合が最適に釣り合っていない場合には、解析マシンは、調整処理３０６で新たに割り振られた傾向スコアと新たに計算された（傾向スコアの変化によって調整された）残差を含むデータ・テーブルに戻る。次いで、新たな残差で最大の絶対相関を識別３０４し、患者の傾向スコアを調整３０６し、より新しい傾向スコアに基づいて新たな患者データの重みを決定３０８し、重み付けされた未処置データを重み付けされていない処置データと比較３０９して、データ集合が今度は最適に釣り合っているかどうかを決定３１０する処理を、解析マシンは繰り返す。この処理は、データ集合同士が十分に類似するまで繰り返され、その時点で解析マシンは回帰ステップ用にデータをフォーマット３１２する。
【００２１】
図４は、本発明の実施形態による図３の傾向スコアリング・コンポーネントの、テーブルのデータを準備するコンポーネントによって実行されるステップのフロー図である。図４に描いたデータ準備ステップでは、選択された患者データを解釈し、相関に関する傾向スコアリング処理での後続するステップ３０４により解釈できるようにそのデータを整える。図４で詳述される上記のステップで記載されるデータ準備は、図５Ａ〜図５Ｅにも例示的なテーブルで示している。図４に示すように、システムは、最初に、やはり図２で詳しく描いたデータ選択コンポーネント１１４の処理中に蓄積された患者症例データを検索４０２する。５人の例の患者に関する患者症例データを含む例示的なテーブルを、図５Ａに示す。図５Ａの説明図は、少なくとも、患者ＩＤ番号、その患者が処置グループのメンバーだったかどうかを識別するための処置指標、人種などの分類別の特徴に付随するデータ、および血圧などの計測数値を、初期段階での各患者の記録が含み得ることを示す。図４内で、処理中の次のステップは、「初期値の」傾向スコアを各患者に割り振る４０４ことである。本発明の一実施形態では、解析マシンが処置グループ内の患者の人数を患者記録の総数で割ることによって、この初期値の傾向スコアが計算される。データ準備処理での次のステップは、各患者記録についての残差を計算することである。本発明の実施形態では、各患者についての処置指標から傾向スコアを差し引くことによって、その残差が計算される。５人の例の患者のそれぞれに関する傾向スコアおよび残差のデータを含む例示的なテーブルを、図５Ｂおよび図５Ｃにそれぞれ示す。図４に戻って、一実施形態では、解析マシンは次に、データ内に記載された各患者の特徴についての特性関数を定義し、評価する。上記の特性関数が定義され付随する列が生成されていることが図５Ｄの例示的なテーブルで示されており、その図には各条件についての新たな列が示されているが、それは解析マシンが分類別の特徴と計測数値の両方をブール式で記載できるように定義されなければならない。別の実施形態では、公開された簡便なアルゴリズムを用いることによりすべての特性関数を評価することは回避され、この場合には、選択された特性関数だけが定義され、評価される。特性関数は、特性関数に含まれる条件が各患者記録について真かまたは偽かを決定することによって評価され、１または０を列に配置して真または偽を表す。図４の処理での次のステップは、各特性関数についてそれぞれ他の特性関数との積を定義し、評価４１０することである。特性関数の定義と同様に、システムは、各条件についての新たな列を生成し、各条件を評価して、その行で識別される患者が特性関数（この場合、各列の一番上にある１つまたは複数の特性関数の積）で識別される条件を満たすかどうかを示す０または１を付随する列に配置する。このステップで生成されたその列と患者記録の評価を、図５Ｅに示す。
【００２２】
図６は、本発明の実施形態による、処置を受けた患者についての患者データと、処置を受けなかった患者についての重み付けされた患者データとを比較するための例示的なテーブルの説明図である。図６は、データが最適な釣り合いに到達したかどうかを決定３１０する際に図３の比較ステップ３０９で考慮されたデータを示す。図６の例示的なテーブル中のデータは、回帰用にデータを準備３１２するための重み付けの後で十分に釣り合っているデータ集合を示している。
【００２３】
図７は、図１の回帰モデリング１１８用の例示的な患者症例データのテーブルの説明図である。この説明図は図３で回帰ステップ用にデータをフォーマット３１２する間にフォーマットされたデータを示す。この例では、観察された識別される効果についてのデータで、図５Ａからの元々のデータが補われる。図７のテーブルは各患者についての２つの記録を示しており、１つは「処置前の期間」についてのもの、１つは「処置後の期間」についてのものである。図７の「期間」の列は１または０を含み、それぞれ、記録が「処置後」の期間と「処置前」の期間のどちらを記述しているかを示す。図７のテーブルは、観察期間中に患者が識別される効果を経験したかどうかを示す列／指標「Ｙ」、および各患者についての観察期間の長さを記載した列／指標「Ｅ」も含む。図７のテーブルは、傾向スコアの重みについての列／指標「Ｗ」も含み、これは図３に描いて上述したシステム３０８により算出される。
【００２４】
第２の処置を受ける患者に比べて第１の処置を受ける患者が識別される効果を経験する相対的な可能性は、比率比を推定することによって計算される。その比率比（「ＲＲＲ」）は、処置前の期間に生じる識別されるイベントの、処置後の期間に生じる識別されるイベントに対する比率を、処置（第１の処置）グループおよび未処置（第２の処置）グループについて比較したものである。ＲＲＲを計算するための公式は、（比率（後，処置）／比率（前，処置））／（比率（後，未処置）／比率（前，未処置））である。計算後に、解析マシンは、この推定されたＲＲＲを、オペレータ端末、または第１および第２の識別される処置に関係する効力や危険性を評価するのに用いることができる別の場所に出力する。医療分野の意志決定者は、解析マシンからの推定されたＲＲＲの出力を利用して、特定の処置を患者や医療機構のために承認または推薦するかどうかを決定することができる。
【００２５】
潜在的な交絡を調整するために、システムは、図１に示すような重み付けポアソン回帰モデル１１８を推定することによって、２重にロバストに調整された比率比を取得する。その推定値は、数式Ｌ（ｂ，β）＝ＳｕｍｉＷｉ（Ｙｉｆｉ−ｅｘｐ（ｆｉ））を最大化することから導出され、ここでｆｉ＝ｌｏｇ（Ｅｉ）＋ｂ０＋ｂ１Ｔｉ＋ｂ２Ｐｅｒｉｏｄｉ＋ｂ３Ｔｉ＊Ｐｅｒｉｏｄｉ＋β’Ｘである。結果として得られるｅｘｐ（ｂ３）の値は、２重にロバストな調整された比率比である。
【００２６】
計算後に、解析マシンは、２重にロバストな調整された比率比を、オペレータによる点検用にオペレータ端末に、または第１および第２の識別される処置に関係する効力や危険性を評価するのに用いることができる別の場所に出力する。医療分野の意志決定者は、解析マシンから出力されたこの２重にロバストな調整された比率比を利用して、特定の処置を患者や医療機構のために承認または推薦するかどうかを決定することができる。
【００２７】
例示的な実施形態を用いることで本発明を説明してきたが、本明細書に添付された特許請求の範囲およびそれらの均等物により定められる本発明の範囲および精神の中に含まれる色々な修正を説明された実施形態に加えてもよいことが、当業者であれば理解されよう。
【図１】
2012530300.tif 000003
【図２】
2012530300.tif 000004
【図３】
2012530300.tif 000005
【図４】
2012530300.tif 000006
【図５Ａ】
2012530300.tif 000007
【図５Ｂ】
2012530300.tif 000008
【図５Ｃ】
2012530300.tif 000009
【図５Ｄ】
2012530300.tif 000010
【図５Ｅ】
2012530300.tif 000011
【図６】
2012530300.tif 000012
【図７】
2012530300.tif 000013
【国際調査報告】
2012530300.tif 000014