臨床検査値の管理装置、管理方法及び管理プログラム
【課題】検査過誤を高精度に検出することができる臨床検査値の管理装置、管理方法及び管理プログラムを提供する。
【解決手段】処理済及び未処理検査値が被験者、検査項目、及び採取年月日と対応させて記録部に記録されており、演算処理部は、第１又は第２検査項目且つ第１又は第２採取年月日に対応する処理済検査値を読み出して検査項目及び被験者ごとの処理済検査値差を計算し、相関係数の絶対値が第１のしきい値よりも大きい場合、第１及び第２検査項目を検査過誤の判断に使用可能な組として決定して第２主成分の標準偏差を計算し、対応する処理済及び未処理検査値を読み出して検査項目及び被験者ごとの第３及び第４採取年月日の処理済及び未処理検査値差を計算し、第２主成分を標準偏差で除した値の絶対値が第２のしきい値よりも大きい場合、その処理済及び未処理検査値のいずれかに過誤があると判断する。
【解決手段】処理済及び未処理検査値が被験者、検査項目、及び採取年月日と対応させて記録部に記録されており、演算処理部は、第１又は第２検査項目且つ第１又は第２採取年月日に対応する処理済検査値を読み出して検査項目及び被験者ごとの処理済検査値差を計算し、相関係数の絶対値が第１のしきい値よりも大きい場合、第１及び第２検査項目を検査過誤の判断に使用可能な組として決定して第２主成分の標準偏差を計算し、対応する処理済及び未処理検査値を読み出して検査項目及び被験者ごとの第３及び第４採取年月日の処理済及び未処理検査値差を計算し、第２主成分を標準偏差で除した値の絶対値が第２のしきい値よりも大きい場合、その処理済及び未処理検査値のいずれかに過誤があると判断する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、統計処理を用いた臨床検査値の管理に関し、特に、主成分分析を用いて臨床検査値を処理し、その結果から検査過誤の有無を判断する管理装置、管理方法及び管理プログラムに関する。
【背景技術】
【０００２】
医療機関においては、日々様々な臨床検査が行われ、その結果が診療に利用されている。測定結果である臨床検査値（以下、単に「検査値」とも記す）には測定の不備や検査の対象者（以下「被験者」とも記す）の取り違いにより発生した、間違った検査値（以下「検査過誤」と記す）が一定の割合で含まれている。これに対し医療機関では検査値の信頼性を確保するため検査過誤が疑われた検査値に対し再検査を施行する。沢山の再検査を施行すれば検査値の信頼性は向上する。しかし、再検査のための時間的、経済的コスト（以下「再検査コスト」と記す。）が増大する。従って、診療に貢献するには検査過誤を正確に検出することができる高精度の精度管理法を用いて検査値の信頼性の確保と再検査コストの低減の相反する要求を両立させることが必要である。
【０００３】
検査値を管理する方法は、例えば下記特許文献１〜４に記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第２８２８６０９号明細書
【特許文献２】特開平８−１４７３９６号公報
【特許文献３】特開平１１−４５３０２号公報
【特許文献４】特開２００６−３１２６４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記したように、診療に貢献するには、検査結果の信頼性確保と再検査コスト低減の相反する要求を両立させることが必要である。しかし、上記の特許文献１〜４に記載された管理方法では、検査過誤の見逃しが多く間違った検査値を高精度に検出することができない。
【０００６】
本発明は、上記した課題を解決すべく、検査値の集合を統計処理することによって、検査過誤を高精度に検出することができる臨床検査値の管理装置、管理方法及び管理プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明者は、上記した課題を解決すべく鋭意研究した結果、ある採取年月日における検査値と別の採取年月日における検査値とから求めた差の間に高い相関を有する検査項目の組合せが存在することに着目した。これに主成分分析を組み合わせることにより検査過誤を高精度に検出できることを見出し、これに基づき本願発明をするに至った。なお、下記において、括弧を付して図面の符号を記載しているが、これは本発明の理解を容易にするためのものであり、本発明を限定するためのものではない。ここで用いる主成分分析により第１主成分、第２主成分が計算される。
【０００８】
本発明に係る第１の臨床検査値の管理装置は、臨床検査値の検査過誤を検出可能な装置であって、
記録部（１３）と演算処理部（１１）とを備え、
前記記録部（１３）に、複数の処理済検査値及び未処理検査値が被験者、検査項目、及び採取年月日と対応させて記録されており、
前記演算処理部（１１）は
第１又は第２検査項目（ｊ、ｏ）且つ第１又は第２採取年月日（ｋ、ｍ）に対応する前記複数の処理済検査値を、前記記録部（１３）から読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目（ｊ、ｏ）各々における前記第１及び第２採取年月日（ｋ、ｍ）の前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を前記被験者（ｉ）ごとに計算し、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）の相関係数（ｒ）を計算し、前記相関係数（ｒ）の絶対値が第１のしきい値（Ｔｈ１）よりも大きい場合、前記第１及び第２検査項目（ｊ、ｏ）を検査過誤の判断に使用可能な組として決定し、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）の第２主成分（ａｉｊｏ）を、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を用いて計算し、前記第２主成分（ａｉｊｏ）の標準偏差（σｊｏ）を計算し、前記組として決定された第１及び第２検査項目（ｊ、ｏ）に該当する前記処理済検査値及び未処理検査値を前記記録部（１３）から読み出し、前記第１及び第２検査項目（ｊ、ｏ）各々における第３採取年月日（ｒ）の前記処理済検査値及び第４採取年月日（ｓ）の前記未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）を前記被験者（ｐ）ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を計算し、前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を前記標準偏差（σｊｏ）で除した値の絶対値が、第２のしきい値（Ｔｈ２）よりも大きい場合、前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）に対応する前記処理済検査値及び未処理検査値のいずれかに検査過誤があると判断し、
前記第１のしきい値（Ｔｈ１）が０．５以上であり、
前記第２のしきい値（Ｔｈ２）が２以上である
ことを特徴としている。なお、本発明の「処理済検査値」とは、既に本発明が適用され、あるいは既に医師などに報告されて診療に利用されて、再検査対象外となっている検査値を意味し、「未処理検査値」とは、医師などに報告される前、即ち診療に利用される前の検査値であり、本発明が適用されておらず、検査過誤の検出が必要なものを意味する。
【０００９】
また、本発明に係る第２の臨床検査値の管理装置は、上記第１の臨床検査値の管理装置において、
前記演算処理部（１１）が、読み出された前記複数の処理済検査値を前記検査項目ごとにトランケーションして、前記各処理済検査値をパラメータ作成用検査値とするか否かを判断すること、及び、
歪度を基に前記検査項目ごとの前記パラメータ作成用検査値群の分布形を判断し、該判断結果に応じて前記各パラメータ作成用検査を変換することを含むことを特徴としている。
【００１０】
また、本発明に係る第３の臨床検査値の管理装置は、前記第１採取年月日（ｋ）と前記第２採取年月日（ｍ）との差、及び前記第３採取年月日（ｒ）と前記第４採取年月日（ｓ）との差がともに許容日差の範囲内にあることを特徴としている。なお、本発明において、「許容日差」とは同一被験者及び同一検査項目で２つの検査値それぞれに対応する採取年月日の差を示す整数値であり、「許容日差関係にある」とは許容日差が整数値の範囲内であることを示す。
【００１１】
また、本発明に係る臨床検査値の管理方法は、臨床検査値の検査過誤を検出可能な方法であって、
複数の処理済検査値及び未処理検査値が、被験者、検査項目、及び採取年月日と対応させて記録された記録部（１３）から、第１又は第２検査項目（ｊ、ｏ）且つ第１又は第２採取年月日（ｋ、ｍ）に対応する前記複数の処理済検査値を読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目（ｊ、ｏ）各々における前記第１及び第２採取年月日（ｋ、ｍ）の前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を前記被験者（ｉ）ごとに計算し、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）の相関係数（ｒ）を計算する第１ステップと、
前記相関係数（ｒ）の絶対値が第１のしきい値（Ｔｈ１）よりも大きいか否かを判断し、大きい場合、前記第１及び第２検査項目（ｊ、ｏ）を検査過誤の判断に使用可能な組として決定する第２ステップと、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）の第２主成分（ａｉｊｏ）を、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を用いて計算する第３ステップと、
前記第２主成分（ａｉｊｏ）の標準偏差（σｊｏ）を計算する第４ステップと、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）に該当する前記処理済検査値及び未処理検査値を前記記録部（１３）から読み出し、前記第１及び第２検査項目（ｊ、ｏ）各々における第３採取年月日（ｒ）の前記処理済検査値及び第４採取年月日（ｓ）の前記未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）を前記被験者（ｐ）ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を計算する第５ステップと、
前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を前記標準偏差（σｊｏ）で除した値の絶対値が、第２のしきい値（Ｔｈ２）よりも大きい場合、前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）に対応する処理済検査値及び未処理検査値のいずれかに検査過誤があると判断する第６ステップと、を含み、
前記第１のしきい値（Ｔｈ１）が０．５以上であり、
前記第２のしきい値（Ｔｈ２）が２以上である
ことを特徴としている。
【００１２】
また、本発明に係る臨床検査値の管理プログラムは、臨床検査値の検査過誤を検出可能なコンピュータプログラムであって、
コンピュータに、
複数の処理済検査値及び未処理検査値が、被験者、検査項目、及び採取年月日と対応させて記録された記録部（１３）から、第１又は第２検査項目（ｊ、ｏ）且つ第１又は第２採取年月日（ｋ、ｍ）に対応する前記複数の処理済検査値を読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目（ｊ、ｏ）各々における前記第１及び第２採取年月日（ｋ、ｍ）の前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を前記被験者（ｉ）ごとに計算し、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）の相関係数（ｒ）を計算する第１の機能と、
前記相関係数（ｒ）の絶対値が第１のしきい値（Ｔｈ１）よりも大きいか否かを判断し、大きい場合、前記第１及び第２検査項目（ｊ、ｏ）を検査過誤の判断に使用可能な組として決定する第２の機能と、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）の第２主成分（ａｉｊｏ）を、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を用いて計算する第３機能と、
前記第２主成分（ａｉｊｏ）の標準偏差（σｊｏ）を計算する第４機能と、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）に該当する前記処理済検査値及び未処理検査値を前記記録部（１３）から読み出し、前記第１及び第２検査項目（ｊ、ｏ）各々における第３採取年月日（ｒ）の前記処理済検査値及び第４採取年月日（ｓ）の前記未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）を前記被験者（ｐ）ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を計算する第５の機能と、
前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を前記標準偏差（σｊｏ）で除した値の絶対値が、第２のしきい値（Ｔｈ２）よりも大きい場合、前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）に対応する処理済検査値及び未処理検査値のいずれかに検査過誤があると判断する第６の機能と、を実現させ、
前記第１のしきい値（Ｔｈ１）が０．５以上であり、
前記第２のしきい値（Ｔｈ２）が２以上である
ことを特徴としている。
【発明の効果】
【００１３】
本発明によれば、臨床検査値群の中から、検査過誤を高精度に検出することができる。検出可能な検査過誤は、医師などに報告されて診療に利用された過去の検査値の検査過誤に限定されず、医師などに報告される以前の、未だ診療に利用されていない検査値の検査過誤も含まれる。この後者の検査過誤を検出できる点が、臨床検査値の管理において特に重要である。従って、臨床検査値の管理を高精度且つ効率的に行うことができ、検査結果の信頼性確保と再検査コスト低減の相反する要求を両立させることが可能になる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施の形態に係る臨床検査値の管理装置の概略構成を示すブロック図である。
【図２】図１に示した管理装置による、検査過誤の検出機能を示すフローチャートである。
【図３】図２に示した処理を説明する相関図である。
【発明を実施するための形態】
【００１５】
以下、本発明に係る実施の形態を、添付した図面に基づいて説明する。
【００１６】
図１は、本発明の実施の形態に係る臨床検査値の管理装置（以下、単に「管理装置」とも記す）の概略構成を示すブロック図である。管理装置（１）は、演算処理部（１１）と、一時記憶部（１２）と、記録部（１３）と、操作部（１４）と、表示部（１５）と、入出力インタフェース部（１６）（以下、入出力ＩＦ部（１６）と記す）と、通信インタフェース部（１７）（以下、通信ＩＦ部（１７）と記す）と、これら各部間でデータを交換するための内部バス（１８）とを備えて構成されている。
【００１７】
操作部（１４）は、演算処理部（１１）に対する指示やデータを入力するための手段である。表示部（１５）は、演算処理部（１１）による処理結果などを表示する。入出力ＩＦ部（１６）は、操作部（１４）及び表示部（１５）との、インタフェースを担う。通信ＩＦ部（１７）は、ネットワーク（ＬＡＮ、イントラネット、インターネットなど）を介して、外部の測定装置（２）やデータベース（３）（以下、ＤＢ（３）と記す）とデータ交換を行う。測定装置（２）は、臨床検査用の装置であり、例えば、血液分析装置、血液凝固測定装置、染色装置など（東芝製TBA-200FR NEO、Sysmex製 XE-5000など）である。また、測定装置（２）は、複数備えていてもよい。
【００１８】
管理装置（１）を、例えばコンピュータを利用して構成する場合、演算処理部（１１）、一時記憶部（１２）、及び記録部（１３）には、それぞれＣＰＵ、ＲＡＭ、及びハードディスクドライブを用いることができる。また、操作部（１４）には、コンピュータ用のキーボード、マウス、タッチパネルなどを使用することができる。入出力ＩＦ部（１６）には、操作部（１４）に応じたシリアル若しくはパラレルインタフェースを採用すればよい。また、入出力ＩＦ部（１６）は、ビデオメモリ及びＤＡ変換部を備え、表示部（１５）のビデオ方式に応じたアナログ信号を出力することによって、表示部（１５）に情報を提示するための画像が表示される。
【００１９】
図２は、管理装置（１）の機能、即ち、検査過誤の検出機能を示すフローチャートである。以下においては、特に断らない限り、管理装置（１）の演算処理部（１１）が行う処理として説明する。演算処理部（１１）は、操作部（１４）が操作されて通信ＩＦ部（１７）を介してデータを取得して記録部（１３）に記録し、適宜記録部（１３）からデータを一時記憶部（１２）に読み出し、所定の処理を行った後、その結果を記録部（１３）に記録する。また、演算処理部（１１）は、操作部（１４）の操作を促す画面データや処理結果を表示する画面データを生成し、入出力ＩＦ部（１６）のビデオＲＡＭを介して、これらの画像を表示部（１５）に表示する。また、後述するしきい値などは、予め初期値として記録部（１３）に記録されているとする。
【００２０】
以下の説明で、処理対象とする検査項目は、例えば次の合計１１８項目である。但し、これらに限定される訳ではなく、これらの一部が別の検査項目で置き換えられても、これらに別の検査項目が追加されても、これらの一部が削除されてもよい。
・血清学検査（１９項目）：TPAb、HbeAg、HbeAg、HbsAg、HbsAb、HTLV-1Ab、HIVAb、HCVAb、HCVcAg、BNP尿HCG、血清HCG、FT3、FT4、TSH、AFP、CEA、CA19-9、IRI
・生化学検査（４２項目）：TP、K、ALB、Cl、TBIL、S-OSM、DBIL、HDLc、AST、UIBC、ALT、Fe、CHE、U-OSM、ALP、24H-CCR、LAP、1H-CCR、γ-GT、アンモニア、CK、髄液蛋白定量、LD、髄液糖定量、CHOL、髄液Na、TGL、髄液K、AMY、髄液Cl、IP、FBS、Ca、RBP、UN、TTR、CRE、TRF、UA、LDLc、Na、CRP
・血液学検査（５７項目）：PT-T、PT-%、APTT-T、APTT-%、フィブリノーゲン、FDP-E、AT3、α2PI、PLG、出血時間、PT比、PT(INR)、第13因子、プロテインC、D-ダイマー、プラスミンα2PI複合体、トロンビンATIII複合体、可溶性トロンボモジュリン、WBC、RBC、HGB、HCT、MCV、MCH、MCHC、RDW、PLT、MPV、PDW、PCT、好中球、リンパ球、単球、好酸球、好塩基球、網赤血球数、網赤血球数%、ヘモグロビンA1C、IRF、好塩基球（用手）、好酸球（用手）、骨髄球（用手）、後骨髄球（用手）、かん状該球（用手）、分葉核球（用手）、リンパ球（用手）、単球（用手）、異形リンパ球（用手）、異常リンパ球（用手）、異常細胞（用手）、芽球（用手）、前骨髄球（用手）、赤芽球（用手）、NNC（用手）、MEG（用手）、plasma cell（用手）顆粒リンパ球（用手）。
【００２１】
記録部（１３）には、予め、検体コード、検査項目コード、検査値、検査日、及び採取年月日の情報がそれぞれの対応が分かるように記録されている。例えば、検体コード、検査項目コード、検査値、検査日、及び採取年月日がテーブルとして記録されているとする。また、記録部（１３）には、予め、検体コード及び被験者コードの対応も分かるように記録されており、例えば、検体コード及び被験者コードがテーブルとして記録されているとする。
【００２２】
以下、図２を参照して、本実施形態に係る検査過誤を検出する機能について説明する。ステップＳ１からステップＳ１５はパラメータ作成のための過程（パラメータ計算過程）であり、ステップＳ１６からステップＳ２１は管理対象検査値に対し検査過誤を検出する過程（管理過程）である。
【００２３】
ステップＳ１において、通信ＩＦ部（１７）を介して、測定装置（２）またはＤＢ（３）から処理済検査値及び未処理検査値を取得し、記録部（１３）に記録する。取得されたデータは、複数の検体に関する検査データである。従って、複数の検体についての検査値が、検体を特定する情報（例えば検体コード）、検査項目を特定する情報（例えば検査項目コード）、検査日（例えば年月日）、及び採取年月日と対応させて、記録部（１３）に記録される。すなわち、複数の検体、検査項目、検査結果、検査日、及び採取年月日が、それらの対応関係が分かるように記録される。また、検体及び被験者についても、それらの対応関係が分かるように記録部（１３）に記録される。
【００２４】
ステップＳ２において、検査項目（上記の１１８項目）の中から１つの検査項目を指定し、処理済検査値を採取年月日とともに記録部（１３）から読み出す。ステップＳ２〜Ｓ４の処理は、全ての検査項目について行うこととする。
【００２５】
ステップＳ３において、ステップＳ２で読み出された処理済検査値に対してトランケーションを行う。具体的には、検査項目ごとに集計し平均と標準偏差を算出する。これを基に平均＋３倍の標準偏差以上を示す処理済検査値、及び平均−３倍の標準偏差以下を示す処理済検査値を除外する。これにより異常な検査値が含まれることが多い極端な外れ値を除外することができる。この過程がトランケーションである。なお、検査項目によっては、トランケーション処理を行わなくてもよい場合もある。トランケーション処理を完了した処理済検査値をパラメータ作成用検査値とする。
【００２６】
ステップＳ４において、ステップＳ３のパラメータ作成用検査値（ステップＳ４〜Ｓ１５の処理において、検査値とはパラメータ作成用検査値を意味することとする）に対し、検査項目ごとの歪度を公知の次式１で計算する。このとき、検査値をそのまま用いて歪度を計算した結果（以下、歪度１と記す）と検査値を自然対数の底を用いて対数変換した変換値から歪度を計算した結果（以下、歪度２と記す）を比較する。具体的には、歪度１の絶対値と歪度２の絶対値のうち、どちらが０に近いかを判断する。歪度１の絶対値の方が０に近ければ、その検査項目の検査値の分布は正規分布に近似しているとして、記録部（１３）に記録されている検査値をそのままステップＳ６からステップＳ２１の計算に用いる。これに対し、歪度２の絶対値の方が０に近ければ、その検査項目の検査値の分布は対数正規分布に近似しているとして検査値を対数変換した値をステップＳ６からステップＳ２１の計算に用いる。但し、対数変換以前の検査値が０．０００である場合、対数変換不可能であるため近似値として０．００１を対数変換した値を用いる。なお、最初から正規分布する検査項目のみを処理対象とする場合には、ステップＳ４の処理は不要である。
【００２７】
これは、検査値の分布が正規分布に近似する検査項目と、対数正規分布に近似する検査項目があり、本法では正規分布に近似することを前提としているため、対数正規分布に近似する検査項目の検査値を対数変換により正規分布に近似させる措置である。
ＳｑＸ＝［Σｉ（ｘｉ−ｘａｖ）３／ｎ０］／ｘσ３ （式１）
ここで、Ｘは検査項目、ｘｉ（ｉ＝１〜ｎ０）は多数の検体を測定して得られた検査項目（Ｘ）の検査値を示す。ｘａｖは検査項目（Ｘ）の平均値である。ｎ０は検査値数である。ｘσは検査項目（Ｘ）の標準偏差である。Σｉは、添え字（ｉ）についての和を意味する。
【００２８】
ステップＳ５において、検査項目（上記の１１８項目）の中から１つの検査項目（ｊとする）を選択する。
【００２９】
ステップＳ６において、ステップＳ５で指定した検査項目（ｊ）の第１採取年月日の検査値「ｘｉｊｋ（ｉ＝１〜ｎ１）」と第２採取年月日の検査値「ｙｉｊｍ（ｉ＝１〜ｎ１）」との差である検査値差「ｚｉｊｋｍ（ｉ＝１〜ｎ１）」を算出し、このｚｉｊｋｍを記録部（１３）に記録する。ここで、ｘｉｊｋは被験者（ｉ）、検査項目（ｊ）、及び第１採取年月日（ｋ）の検査値、ｙｉｊｍは被験者（ｉ）、検査項目（ｊ）、第２採取年月日（ｍ）の検査値を意味しており、ｋ＜ｍ、即ち第１採取年月日は第２採取年月日より過去の採取年月日である。である。ｚｉｊｋｍは、同一被験者（ｉ）及び同一検査項目（ｊ）における第１採取年月日（ｋ）の検査値と第２採取年月日（ｍ）の検査値との差、すなわち、ｘｉｊｋとｙｉｊｍとの差である。ｎ１は被験者数である。採取年月日が同日（ｋ＝ｍ）の場合は、時差、即ちｘｉｊｋとｙｉｊｍとの検査時刻間差で区別する。ステップＳ６の処理は、検査項目ごとに、重複しない全ての第１及び第２採取年月日の組み合わせについて行う。なお、第１の採取年月日及び第２の採取年月日は許容日差の範囲内にあることが好ましく、この許容日差は、特に限定されるものではないが、通常０日〜７日程度である。
【００３０】
ステップＳ７において、全ての検査項目についてステップＳ５及びステップＳ６の処理が終了したか否かを判断し、終了した場合、ステップＳ８に移行する。未処理の検査項目があればステップＳ５に戻る。ステップＳ５に戻った場合、未処理の検査項目の中から１つの検査項目を選択し、ステップＳ６の処理を実行する。
【００３１】
ステップＳ８において、検査項目（上記の１１８項目）の中から異なる２つの検査項目（ｊ、ｏ）を選択し、それらの検査値差「ｚｉｊｋｍ（ｉ＝１〜ｎ1）」及び「ｚｉｏｋｍ（ｉ＝１〜ｎ1）」を記録部（１３）から読み出す。ここで、ｚｉｊｋｍは、同一被験者（ｉ）及び同一検査項目（ｊ）における第１採取年月日（ｋ）の検査値と第２採取年月日（ｍ）の検査値との差であり、同様に、ｚｉｏｋｍは、同一被験者（ｉ）及び同一検査項目（ｏ）における第１採取年月日（ｋ）の検査値と第２採取年月日（ｍ）の検査値との差である。第１採取年月日（ｋ）及び第２採取年月日（ｍ）は、被験者ｉ、検査項目ｊ、及び検査項目ｏについて同一とする。なお、検査項目（ｊ、ｏ）は本発明の第１及び第２検査項目に相当する。
【００３２】
ステップＳ９において、ステップＳ８で読み出したｎ１個の検査値差（ｚｉｊｋｍ）と検査値差（ｚｉｏｋｍ）との検査値差対（ペア）「ｚｉｊｋｍ、ｚｉｏｋｍ（ｉ＝１〜ｎ1）」を作成する。なお、ステップＳ９の処理は、選択した検査項目（ｊ、ｏ）ごとに、重複しない全ての第１及び第２採取年月日の組み合わせについて行う。
【００３３】
ステップＳ１０において、ステップＳ９で作成された検査値差対「ｚｉｊｋｍ、ｚｉｏｋｍ（ｉ＝１〜ｎ1）」を用いて、相関係数（ｒ）を公知の次式２で計算する。
ｒ＝Σｉ（ｚｉｊｋｍ−ｚｊ）（ｚｉｏｋｍ−ｚｏ）／｛Σｉ（ｚｉｊｋｍ−ｚｊ）２Σｉ（ｚｉｏｋｍ−ｚｏ）２｝１／２ （式２）
ここで、ｚｊ、ｚｏは各検査値差の平均値であり、Σｉは、添え字（ｉ）についての和を意味する。
【００３４】
ステップＳ１１において、ステップＳ１０で計算された相関係数（ｒ）の絶対値が所定のしきい値Ｔｈ１（例えば０．５）よりも大きく、且つ、２つの検査項目の検査値差対の数が所定数以上（例えば５０以上、望ましくは１００以上）である場合、ステップＳ１２に移行し、そうでなければステップＳ１３に移行する。検査値差対の数についての条件は、統計的に有効な結果が得られるようにするためである。
【００３５】
ステップＳ１２において、ステップ１１で選択された２つの検査項目（ｊ、ｏ）を後述する検査過誤の可能性の判定の対象として決定し、例えば、２つの検査項目（ｊ、ｏ）の情報を対応付けて記録部（１３）に記録する。
【００３６】
ステップＳ１３において、重複せずに選択された全ての検査項目の組み合わせについて、ステップＳ８〜ステップＳ１２の処理が終了したか否かを判断し、終了した場合、ステップＳ１４に移行する。残っている組み合わせがあればステップＳ８に戻る。ステップＳ８に戻った場合、再び検査項目の中から異なる２つの検査項目を選択し、ステップＳ９〜ステップＳ１３の処理を実行する。このとき、既にステップＳ８で選択された項目と重複しないように、２つの検査項目を選択する。
【００３７】
以上のように、ステップＳ１〜ステップＳ１３によって、相関の高い２つの検査項目の組が複数決定される。以下のステップＳ１４〜ステップＳ１６の処理においては、ここで決定された組が使用される。
【００３８】
ステップＳ１４において、ステップＳ１２で決定された組の検査値差対の第２主成分を計算する。これをａｉｊｏとする。ここで、ａｉｊｏは、２つの検査項目（ｊ、ｏ）の対応する検査値差「ｚｉｊｋｍ、ｚｉｏｋｍ（ｉ＝１〜ｎ１）」に対して、公知の主成分分析を適用して決定された第２主成分（短軸成分）である。ここで、対応する検査値差（ｚｉｊｋｍ、ｚｉｏｋｍ）とは、被験者（ｉ）、第１採取年月日（ｋ）、及び第２採取年月日（ｍ）が同一の検査値を意味する。検査値差（ｚｉｊｋｍ、ｚｉｏｋｍ）のそれぞれの平均をｚｊ、ｚｏ、標準偏差をｚｊσ、ｚｏσとして、ａｉｊｏは、
相関係数（ｒ）＞０．５の場合、
ａｉｊｏ＝｛（ｚｉｊｋｍ−ｚｊ）／ｚｊσ−（ｚｉｏｋｍ−ｚｏ）／ｚｏσ｝／２１／２ で、
相関係数（ｒ）＜−０．５の場合、
ａｉｊｏ＝｛（ｚｉｊｋｍ−ｚｊ）／ｚｊσ＋（ｚｉｏｋｍ−ｚｏ）／ｚｏσ｝／２１／２ で計算される。
なお、ステップＳ１４の処理は、組として決定された検査項目（ｊ、ｏ）ごとに、重複しない全ての第１及び第２採取年月日の組み合わせについて行う。
【００３９】
ステップＳ１５において、ステップＳ１４で計算された第２主成分「ａｉｊｏ（ｉ＝１〜ｎ１）」の標準偏差を計算する。これを２つの検査項目（ｊ、ｏ）の第２主成分の標準偏差（σｊｏ）としてパラメータに設定し、記録部（１３）に記録する。
【００４０】
ステップＳ１６において、ステップＳ１２で決定された複数の組のうち、１つの組（検査項目（ｊ、ｏ））を選択して記録部（１３）から読み出し、さらに、この組に対応する、被験者（ｐ）の処理済検査値及び未処理検査値を記録部（１３）から読み出す。読み出した処理済検査値及び未処理検査値について、検査項目（ｊ、ｏ）ごとに任意に選択した第３採取年月日（ｒ）と第４採取年月日（ｓ）との検査値差（ｚｐｊｒｓ 、ｚｐｏｒｓ）を計算する。このとき、第３採取年月日は処理済検査値の採取年月日であり、第４採取年月日は未処理検査値の採取年月日であり、且つ、ｒ＜ｓである。即ち第３採取年月日は第４採取年月日より過去の採取年月日である。なお、第３の採取年月日及び第４の採取年月日は、特に限定されるものではないが、上述した第１の採取年月日及び第２の採取年月日と同じ許容日差の範囲内にあることが好ましい。この結果対（ｚｐｊｒｓ 、ｚｐｏｒｓ）に対応する検査値は、パラメータ作成用検査値からパラメータであるσｊｏが設定された後、本法によって検査過誤検出の対象となる管理対象検査値である。ステップＳ4において、正規分布に近似していると判断された検査項目は、そのままの検査値を用い、対数正規分布に近似していると判断された検査項目は対数変換された検査値を用いる。
【００４１】
ステップＳ１７において、次式３によってＳＤＩｐｊｏを計算する。
ＳＤＩｐｊｏ＝ａｐｊｏ／σｊｏ （式３）
ここで、ａｐｊｏは、ステップＳ１６において算出されたｚｐｊｒｓ 、ｚｐｏｒｓに対し、以下のように計算される。すなわち、ｚｐｊｒｓ 、ｚｐｏｒｓに対しステップＳ１４で求めた平均ｚｊ、ｚｏ、標準偏差ｚｊσ、ｚｏσを用いて、ａｐｊｏは、
相関係数（ｒ）＞０．５の場合、
ａｐｊｏ＝｛（ｚｐｊｒｓ−ｚｊ）／ｚｊσ−（ｚｐｏｒｓ−ｚｏ）／ｚｏσ｝／２１／２ で、
相関係数（ｒ）＜−０．５の場合、
ａｐｊｏ＝｛（ｚｐｊｒｓ−ｚｊ）／ｚｊσ＋（ｚｐｏｒｓ−ｚｏ）／ｚｏσ｝／２１／２ で計算される。σｊｏは、ステップＳ１５で計算された第２主成分の標準偏差である。主成分分析において、第１主成分（長軸成分）は系統誤差成分、第２主成分（短軸成分）は偶発誤差成分を表すと考えられる。なお、主成分分析は公知であるので説明を省略する。
【００４２】
ステップＳ１８において、ステップＳ１７で計算されたＳＤＩｐｊｏの絶対値が、所定のしきい値Ｔｈ２よりも大きいか否かを判断する。この条件を満たす場合、ステップＳ１９に移行し、そうでなければステップＳ２０に移行する。しきい値Ｔｈ２は、標準偏差（σｊｏ）を基準として設定することができる。しきい値Ｔｈ２は、例えば２以上、望ましくは４以上である。
【００４３】
ステップＳ１９において、ステップＳ１８の条件を満たす検査値差対（ｚｐｊｒｓ、ｚｐｏｒｓ）に対応する処理済検査値及び未処理検査値のいずれかに検査過誤の可能性があると判断し、その情報を記録する。例えば、この検査値差対（ｚｐｊｒｓ 、ｚｐｏｒｓ）に対応する２つの検査項目コード（ｊ、ｏ）、２つの採取年月日（ｒ、ｓ）、及びそれぞれの検査値を対応付けて記録部（１３）に記録する。
【００４４】
図３は、２つの検査項目（ｊ、ｏ）について、上記の処理を概念的に説明するための相関図である。図３では省略しているが、検査項目（ｊ）における第３及び第４採取年月日の検査値差と検査項目（ｏ）における第３及び第４採取年月日の検査値差との検査値差対（ペア）が複数プロットされている。楕円は、プロットされた検査値差対の分布範囲を示している。検査項目（ｊ）の検査値差と検査項目（ｏ）の検査値差は正の相関を示している。ベクトル（α）は分布範囲内の検査値差対を表すベクトルであり、ベクトル（β）は分布範囲を逸脱している検査値差対を表すベクトルである。また、ベクトル（α）、（β）は、短軸方向に平行であり、それぞれの長さは、長軸から検査値差対までの垂直距離（短軸方向の距離）である。このベクトル（α）、（β）の「長さ」が偶発誤差成分（第２主成分）である。ベクトル（β）のように偶発誤差成分が任意の管理限界（しきい値Ｔｈ２）を逸脱する場合、２つの検査値差に対応する処理済検査値と未処理検査値のいずれかに過誤があると判断できる。
【００４５】
ステップＳ２０において、全ての組について処理を終了したか否かを判断し、残っている組があればステップＳ１６に戻り、全て終了していればステップＳ２１に移行する。ステップＳ１６に戻った場合、再び１つの組を選択し、ステップＳ１７〜ステップＳ２０の処理を実行する。このとき、既にステップＳ１６で選択された組を除外して選択する。
【００４６】
ステップＳ２１において、ステップＳ１９で記録された検査項目、検査日（測定年月日）、検査値、及び関連する情報を表示部（１５）に表示する。なお、ステップＳ４において対数変換した検査値は変換前の検査値に戻して表示する。表１に、表示の一例を示す。ここで、第３採取年月日（前回日）と第４採取年月日（今回日）は西暦の年月日が８桁で表されている。
【００４７】
【表１】
【００４８】
以上によって、検査過誤検出の対象となる実検体の検査値対のうち、検査過誤に該当する可能性の高い検査値を提示することができる。
【００４９】
なお、上記処理において、検査値群が異なれば、即ち、医療機関や検査値の測定装置が異なれば、ステップＳ１２で決定される相関の高い検査項目や第２主成分の標準偏差（σｊｏ）が異なり得る。従って、ステップＳ１〜Ｓ１５の処理は、検査値群が追加、変更される毎に行うことが望ましい。その一方、特定の医療機関で同じ測定装置を使用していれば、既に決定された相関の高い２つの検査項目を用いて、ステップＳ１６〜Ｓ２１の処理を行っても、高い信頼性で検査過誤の可能性を提示することができる。また、相関係数（ｒ）の絶対値が０．６以上である組を使用すれば、医療機関によらず高い信頼性で検査過誤の可能性を提示することができる。
【００５０】
以上、本発明の実施の形態に係る臨床検査値の管理装置、及びそれを用いた管理方法について説明したが、本発明は上記に限定されない。例えば、管理装置（１）は、図１の構成から変更して構成され得る。また、図２に示したフローチャートは、種々変更されて実行され得る。
【００５１】
例えば、図１では、管理装置（１）が処理済検査値及び未処理検査値を、ネットワークを介して取得する場合を説明したが、可搬性の記録媒体（光磁気記録媒体、フラッシュメモリなど）を介して取得してもよい。
【００５２】
また、専用の管理装置として実現する代わりに、上記した機能をコンピュータに実現させるプログラムを、既存の臨床検査装置に搭載してもよい。
【００５３】
また、しきい値は任意である。例えば、Ｔｈ１として０．６以上の値を使用すれば、医療機関などに依存せずに一般的に使用できる検査項目及び採取年月日を決定できる。
【００５４】
また、Ｔｈ２として、標準偏差の４倍よりも大きい値を使用することも、標準偏差の４倍よりも小さい値を使用することもできる。但し、Ｔｈ２は、小さ過ぎると無駄な再検査が増加し、大きすぎると結果の信頼性が低下するので好ましくない。少なくとも標準偏差の２倍であることが必要である。
【００５５】
また、図２の判断ステップ（Ｓ１１、Ｓ１８）において、値を計算する度に判断する代わりに、計算結果を記録し、全ての計算が完了した後に、計算結果の全てについて、条件を満たすか否かを判断してもよい。
【００５６】
また、測定された検査値が、正規分布または対数正規分布する場合に限定されず、これらの分布に対応しない場合には、変換後の検査値が正規分布となるように変換すればよい。
【００５７】
また、上記では、未処理検査値（医師に報告されていない検査値）の検査過誤を検出する場合を説明したが、これに限定されない。将来の検査値（将来において測定された場合の検査値）の検査過誤を検出することも可能である。また、処理済検査値（医師に報告されて診療に利用された検査値）についても、同様に検査過誤を検出することが可能である。
【実施例】
【００５８】
以下に実施例を示し、本発明の特徴とするところをより一層明確にする。
【００５９】
２０１０年１月１日〜２０１０年６月３０日の間に、大分大学医学部附属病院で検査された検体の検査値に、本発明を適用した。検体の総数は２０８８２１であり、上記で示した生化学検査、血液学検査、血清学検査の合計１１８項目の検査値を用いた。
【００６０】
表２に、図２のステップＳ１〜Ｓ１３に対応する処理を行った結果、得られた相関の高い２つの検査項目（｜ｒ｜＞０．５）の組の一部を示す。全部でなく一部を示すのは相関の高い組合せ数が膨大であるためである。これらの検査項目の組の検査値を用いて、検査過誤の可能性を判断することが望ましい。なお、｜ｒ｜＞０．６の検査項目の組には、右端列に“１”が表示されている。これらの検査項目の組は、別の医療機関においても、検査過誤の判断に使用することができる。
【００６１】
【表２】
【符号の説明】
【００６２】
１ 管理装置
１１ 演算処理部
１２ 一時記憶部
１３ 記録部
１４ 操作部
１５ 表示部
１６ 入出力インタフェース部（入出力ＩＦ部）
１７ 通信インタフェース部（通信ＩＦ部）
１８ 内部バス
２ 測定装置
３ データベース
４ ネットワーク
【技術分野】
【０００１】
本発明は、統計処理を用いた臨床検査値の管理に関し、特に、主成分分析を用いて臨床検査値を処理し、その結果から検査過誤の有無を判断する管理装置、管理方法及び管理プログラムに関する。
【背景技術】
【０００２】
医療機関においては、日々様々な臨床検査が行われ、その結果が診療に利用されている。測定結果である臨床検査値（以下、単に「検査値」とも記す）には測定の不備や検査の対象者（以下「被験者」とも記す）の取り違いにより発生した、間違った検査値（以下「検査過誤」と記す）が一定の割合で含まれている。これに対し医療機関では検査値の信頼性を確保するため検査過誤が疑われた検査値に対し再検査を施行する。沢山の再検査を施行すれば検査値の信頼性は向上する。しかし、再検査のための時間的、経済的コスト（以下「再検査コスト」と記す。）が増大する。従って、診療に貢献するには検査過誤を正確に検出することができる高精度の精度管理法を用いて検査値の信頼性の確保と再検査コストの低減の相反する要求を両立させることが必要である。
【０００３】
検査値を管理する方法は、例えば下記特許文献１〜４に記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第２８２８６０９号明細書
【特許文献２】特開平８−１４７３９６号公報
【特許文献３】特開平１１−４５３０２号公報
【特許文献４】特開２００６−３１２６４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記したように、診療に貢献するには、検査結果の信頼性確保と再検査コスト低減の相反する要求を両立させることが必要である。しかし、上記の特許文献１〜４に記載された管理方法では、検査過誤の見逃しが多く間違った検査値を高精度に検出することができない。
【０００６】
本発明は、上記した課題を解決すべく、検査値の集合を統計処理することによって、検査過誤を高精度に検出することができる臨床検査値の管理装置、管理方法及び管理プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明者は、上記した課題を解決すべく鋭意研究した結果、ある採取年月日における検査値と別の採取年月日における検査値とから求めた差の間に高い相関を有する検査項目の組合せが存在することに着目した。これに主成分分析を組み合わせることにより検査過誤を高精度に検出できることを見出し、これに基づき本願発明をするに至った。なお、下記において、括弧を付して図面の符号を記載しているが、これは本発明の理解を容易にするためのものであり、本発明を限定するためのものではない。ここで用いる主成分分析により第１主成分、第２主成分が計算される。
【０００８】
本発明に係る第１の臨床検査値の管理装置は、臨床検査値の検査過誤を検出可能な装置であって、
記録部（１３）と演算処理部（１１）とを備え、
前記記録部（１３）に、複数の処理済検査値及び未処理検査値が被験者、検査項目、及び採取年月日と対応させて記録されており、
前記演算処理部（１１）は
第１又は第２検査項目（ｊ、ｏ）且つ第１又は第２採取年月日（ｋ、ｍ）に対応する前記複数の処理済検査値を、前記記録部（１３）から読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目（ｊ、ｏ）各々における前記第１及び第２採取年月日（ｋ、ｍ）の前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を前記被験者（ｉ）ごとに計算し、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）の相関係数（ｒ）を計算し、前記相関係数（ｒ）の絶対値が第１のしきい値（Ｔｈ１）よりも大きい場合、前記第１及び第２検査項目（ｊ、ｏ）を検査過誤の判断に使用可能な組として決定し、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）の第２主成分（ａｉｊｏ）を、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を用いて計算し、前記第２主成分（ａｉｊｏ）の標準偏差（σｊｏ）を計算し、前記組として決定された第１及び第２検査項目（ｊ、ｏ）に該当する前記処理済検査値及び未処理検査値を前記記録部（１３）から読み出し、前記第１及び第２検査項目（ｊ、ｏ）各々における第３採取年月日（ｒ）の前記処理済検査値及び第４採取年月日（ｓ）の前記未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）を前記被験者（ｐ）ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を計算し、前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を前記標準偏差（σｊｏ）で除した値の絶対値が、第２のしきい値（Ｔｈ２）よりも大きい場合、前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）に対応する前記処理済検査値及び未処理検査値のいずれかに検査過誤があると判断し、
前記第１のしきい値（Ｔｈ１）が０．５以上であり、
前記第２のしきい値（Ｔｈ２）が２以上である
ことを特徴としている。なお、本発明の「処理済検査値」とは、既に本発明が適用され、あるいは既に医師などに報告されて診療に利用されて、再検査対象外となっている検査値を意味し、「未処理検査値」とは、医師などに報告される前、即ち診療に利用される前の検査値であり、本発明が適用されておらず、検査過誤の検出が必要なものを意味する。
【０００９】
また、本発明に係る第２の臨床検査値の管理装置は、上記第１の臨床検査値の管理装置において、
前記演算処理部（１１）が、読み出された前記複数の処理済検査値を前記検査項目ごとにトランケーションして、前記各処理済検査値をパラメータ作成用検査値とするか否かを判断すること、及び、
歪度を基に前記検査項目ごとの前記パラメータ作成用検査値群の分布形を判断し、該判断結果に応じて前記各パラメータ作成用検査を変換することを含むことを特徴としている。
【００１０】
また、本発明に係る第３の臨床検査値の管理装置は、前記第１採取年月日（ｋ）と前記第２採取年月日（ｍ）との差、及び前記第３採取年月日（ｒ）と前記第４採取年月日（ｓ）との差がともに許容日差の範囲内にあることを特徴としている。なお、本発明において、「許容日差」とは同一被験者及び同一検査項目で２つの検査値それぞれに対応する採取年月日の差を示す整数値であり、「許容日差関係にある」とは許容日差が整数値の範囲内であることを示す。
【００１１】
また、本発明に係る臨床検査値の管理方法は、臨床検査値の検査過誤を検出可能な方法であって、
複数の処理済検査値及び未処理検査値が、被験者、検査項目、及び採取年月日と対応させて記録された記録部（１３）から、第１又は第２検査項目（ｊ、ｏ）且つ第１又は第２採取年月日（ｋ、ｍ）に対応する前記複数の処理済検査値を読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目（ｊ、ｏ）各々における前記第１及び第２採取年月日（ｋ、ｍ）の前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を前記被験者（ｉ）ごとに計算し、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）の相関係数（ｒ）を計算する第１ステップと、
前記相関係数（ｒ）の絶対値が第１のしきい値（Ｔｈ１）よりも大きいか否かを判断し、大きい場合、前記第１及び第２検査項目（ｊ、ｏ）を検査過誤の判断に使用可能な組として決定する第２ステップと、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）の第２主成分（ａｉｊｏ）を、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を用いて計算する第３ステップと、
前記第２主成分（ａｉｊｏ）の標準偏差（σｊｏ）を計算する第４ステップと、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）に該当する前記処理済検査値及び未処理検査値を前記記録部（１３）から読み出し、前記第１及び第２検査項目（ｊ、ｏ）各々における第３採取年月日（ｒ）の前記処理済検査値及び第４採取年月日（ｓ）の前記未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）を前記被験者（ｐ）ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を計算する第５ステップと、
前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を前記標準偏差（σｊｏ）で除した値の絶対値が、第２のしきい値（Ｔｈ２）よりも大きい場合、前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）に対応する処理済検査値及び未処理検査値のいずれかに検査過誤があると判断する第６ステップと、を含み、
前記第１のしきい値（Ｔｈ１）が０．５以上であり、
前記第２のしきい値（Ｔｈ２）が２以上である
ことを特徴としている。
【００１２】
また、本発明に係る臨床検査値の管理プログラムは、臨床検査値の検査過誤を検出可能なコンピュータプログラムであって、
コンピュータに、
複数の処理済検査値及び未処理検査値が、被験者、検査項目、及び採取年月日と対応させて記録された記録部（１３）から、第１又は第２検査項目（ｊ、ｏ）且つ第１又は第２採取年月日（ｋ、ｍ）に対応する前記複数の処理済検査値を読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目（ｊ、ｏ）各々における前記第１及び第２採取年月日（ｋ、ｍ）の前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を前記被験者（ｉ）ごとに計算し、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）の相関係数（ｒ）を計算する第１の機能と、
前記相関係数（ｒ）の絶対値が第１のしきい値（Ｔｈ１）よりも大きいか否かを判断し、大きい場合、前記第１及び第２検査項目（ｊ、ｏ）を検査過誤の判断に使用可能な組として決定する第２の機能と、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）の第２主成分（ａｉｊｏ）を、前記パラメータ作成用検査値の差（ｚｉｊｋｍ、ｚｉｏｋｍ）を用いて計算する第３機能と、
前記第２主成分（ａｉｊｏ）の標準偏差（σｊｏ）を計算する第４機能と、
前記組として決定された第１及び第２検査項目（ｊ、ｏ）に該当する前記処理済検査値及び未処理検査値を前記記録部（１３）から読み出し、前記第１及び第２検査項目（ｊ、ｏ）各々における第３採取年月日（ｒ）の前記処理済検査値及び第４採取年月日（ｓ）の前記未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）を前記被験者（ｐ）ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を計算する第５の機能と、
前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）を前記標準偏差（σｊｏ）で除した値の絶対値が、第２のしきい値（Ｔｈ２）よりも大きい場合、前記処理済検査値及び未処理検査値の差（ｚｐｊｒｓ、ｚｐｏｒｓ）の第２主成分（ａｐｊｏ）に対応する処理済検査値及び未処理検査値のいずれかに検査過誤があると判断する第６の機能と、を実現させ、
前記第１のしきい値（Ｔｈ１）が０．５以上であり、
前記第２のしきい値（Ｔｈ２）が２以上である
ことを特徴としている。
【発明の効果】
【００１３】
本発明によれば、臨床検査値群の中から、検査過誤を高精度に検出することができる。検出可能な検査過誤は、医師などに報告されて診療に利用された過去の検査値の検査過誤に限定されず、医師などに報告される以前の、未だ診療に利用されていない検査値の検査過誤も含まれる。この後者の検査過誤を検出できる点が、臨床検査値の管理において特に重要である。従って、臨床検査値の管理を高精度且つ効率的に行うことができ、検査結果の信頼性確保と再検査コスト低減の相反する要求を両立させることが可能になる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施の形態に係る臨床検査値の管理装置の概略構成を示すブロック図である。
【図２】図１に示した管理装置による、検査過誤の検出機能を示すフローチャートである。
【図３】図２に示した処理を説明する相関図である。
【発明を実施するための形態】
【００１５】
以下、本発明に係る実施の形態を、添付した図面に基づいて説明する。
【００１６】
図１は、本発明の実施の形態に係る臨床検査値の管理装置（以下、単に「管理装置」とも記す）の概略構成を示すブロック図である。管理装置（１）は、演算処理部（１１）と、一時記憶部（１２）と、記録部（１３）と、操作部（１４）と、表示部（１５）と、入出力インタフェース部（１６）（以下、入出力ＩＦ部（１６）と記す）と、通信インタフェース部（１７）（以下、通信ＩＦ部（１７）と記す）と、これら各部間でデータを交換するための内部バス（１８）とを備えて構成されている。
【００１７】
操作部（１４）は、演算処理部（１１）に対する指示やデータを入力するための手段である。表示部（１５）は、演算処理部（１１）による処理結果などを表示する。入出力ＩＦ部（１６）は、操作部（１４）及び表示部（１５）との、インタフェースを担う。通信ＩＦ部（１７）は、ネットワーク（ＬＡＮ、イントラネット、インターネットなど）を介して、外部の測定装置（２）やデータベース（３）（以下、ＤＢ（３）と記す）とデータ交換を行う。測定装置（２）は、臨床検査用の装置であり、例えば、血液分析装置、血液凝固測定装置、染色装置など（東芝製TBA-200FR NEO、Sysmex製 XE-5000など）である。また、測定装置（２）は、複数備えていてもよい。
【００１８】
管理装置（１）を、例えばコンピュータを利用して構成する場合、演算処理部（１１）、一時記憶部（１２）、及び記録部（１３）には、それぞれＣＰＵ、ＲＡＭ、及びハードディスクドライブを用いることができる。また、操作部（１４）には、コンピュータ用のキーボード、マウス、タッチパネルなどを使用することができる。入出力ＩＦ部（１６）には、操作部（１４）に応じたシリアル若しくはパラレルインタフェースを採用すればよい。また、入出力ＩＦ部（１６）は、ビデオメモリ及びＤＡ変換部を備え、表示部（１５）のビデオ方式に応じたアナログ信号を出力することによって、表示部（１５）に情報を提示するための画像が表示される。
【００１９】
図２は、管理装置（１）の機能、即ち、検査過誤の検出機能を示すフローチャートである。以下においては、特に断らない限り、管理装置（１）の演算処理部（１１）が行う処理として説明する。演算処理部（１１）は、操作部（１４）が操作されて通信ＩＦ部（１７）を介してデータを取得して記録部（１３）に記録し、適宜記録部（１３）からデータを一時記憶部（１２）に読み出し、所定の処理を行った後、その結果を記録部（１３）に記録する。また、演算処理部（１１）は、操作部（１４）の操作を促す画面データや処理結果を表示する画面データを生成し、入出力ＩＦ部（１６）のビデオＲＡＭを介して、これらの画像を表示部（１５）に表示する。また、後述するしきい値などは、予め初期値として記録部（１３）に記録されているとする。
【００２０】
以下の説明で、処理対象とする検査項目は、例えば次の合計１１８項目である。但し、これらに限定される訳ではなく、これらの一部が別の検査項目で置き換えられても、これらに別の検査項目が追加されても、これらの一部が削除されてもよい。
・血清学検査（１９項目）：TPAb、HbeAg、HbeAg、HbsAg、HbsAb、HTLV-1Ab、HIVAb、HCVAb、HCVcAg、BNP尿HCG、血清HCG、FT3、FT4、TSH、AFP、CEA、CA19-9、IRI
・生化学検査（４２項目）：TP、K、ALB、Cl、TBIL、S-OSM、DBIL、HDLc、AST、UIBC、ALT、Fe、CHE、U-OSM、ALP、24H-CCR、LAP、1H-CCR、γ-GT、アンモニア、CK、髄液蛋白定量、LD、髄液糖定量、CHOL、髄液Na、TGL、髄液K、AMY、髄液Cl、IP、FBS、Ca、RBP、UN、TTR、CRE、TRF、UA、LDLc、Na、CRP
・血液学検査（５７項目）：PT-T、PT-%、APTT-T、APTT-%、フィブリノーゲン、FDP-E、AT3、α2PI、PLG、出血時間、PT比、PT(INR)、第13因子、プロテインC、D-ダイマー、プラスミンα2PI複合体、トロンビンATIII複合体、可溶性トロンボモジュリン、WBC、RBC、HGB、HCT、MCV、MCH、MCHC、RDW、PLT、MPV、PDW、PCT、好中球、リンパ球、単球、好酸球、好塩基球、網赤血球数、網赤血球数%、ヘモグロビンA1C、IRF、好塩基球（用手）、好酸球（用手）、骨髄球（用手）、後骨髄球（用手）、かん状該球（用手）、分葉核球（用手）、リンパ球（用手）、単球（用手）、異形リンパ球（用手）、異常リンパ球（用手）、異常細胞（用手）、芽球（用手）、前骨髄球（用手）、赤芽球（用手）、NNC（用手）、MEG（用手）、plasma cell（用手）顆粒リンパ球（用手）。
【００２１】
記録部（１３）には、予め、検体コード、検査項目コード、検査値、検査日、及び採取年月日の情報がそれぞれの対応が分かるように記録されている。例えば、検体コード、検査項目コード、検査値、検査日、及び採取年月日がテーブルとして記録されているとする。また、記録部（１３）には、予め、検体コード及び被験者コードの対応も分かるように記録されており、例えば、検体コード及び被験者コードがテーブルとして記録されているとする。
【００２２】
以下、図２を参照して、本実施形態に係る検査過誤を検出する機能について説明する。ステップＳ１からステップＳ１５はパラメータ作成のための過程（パラメータ計算過程）であり、ステップＳ１６からステップＳ２１は管理対象検査値に対し検査過誤を検出する過程（管理過程）である。
【００２３】
ステップＳ１において、通信ＩＦ部（１７）を介して、測定装置（２）またはＤＢ（３）から処理済検査値及び未処理検査値を取得し、記録部（１３）に記録する。取得されたデータは、複数の検体に関する検査データである。従って、複数の検体についての検査値が、検体を特定する情報（例えば検体コード）、検査項目を特定する情報（例えば検査項目コード）、検査日（例えば年月日）、及び採取年月日と対応させて、記録部（１３）に記録される。すなわち、複数の検体、検査項目、検査結果、検査日、及び採取年月日が、それらの対応関係が分かるように記録される。また、検体及び被験者についても、それらの対応関係が分かるように記録部（１３）に記録される。
【００２４】
ステップＳ２において、検査項目（上記の１１８項目）の中から１つの検査項目を指定し、処理済検査値を採取年月日とともに記録部（１３）から読み出す。ステップＳ２〜Ｓ４の処理は、全ての検査項目について行うこととする。
【００２５】
ステップＳ３において、ステップＳ２で読み出された処理済検査値に対してトランケーションを行う。具体的には、検査項目ごとに集計し平均と標準偏差を算出する。これを基に平均＋３倍の標準偏差以上を示す処理済検査値、及び平均−３倍の標準偏差以下を示す処理済検査値を除外する。これにより異常な検査値が含まれることが多い極端な外れ値を除外することができる。この過程がトランケーションである。なお、検査項目によっては、トランケーション処理を行わなくてもよい場合もある。トランケーション処理を完了した処理済検査値をパラメータ作成用検査値とする。
【００２６】
ステップＳ４において、ステップＳ３のパラメータ作成用検査値（ステップＳ４〜Ｓ１５の処理において、検査値とはパラメータ作成用検査値を意味することとする）に対し、検査項目ごとの歪度を公知の次式１で計算する。このとき、検査値をそのまま用いて歪度を計算した結果（以下、歪度１と記す）と検査値を自然対数の底を用いて対数変換した変換値から歪度を計算した結果（以下、歪度２と記す）を比較する。具体的には、歪度１の絶対値と歪度２の絶対値のうち、どちらが０に近いかを判断する。歪度１の絶対値の方が０に近ければ、その検査項目の検査値の分布は正規分布に近似しているとして、記録部（１３）に記録されている検査値をそのままステップＳ６からステップＳ２１の計算に用いる。これに対し、歪度２の絶対値の方が０に近ければ、その検査項目の検査値の分布は対数正規分布に近似しているとして検査値を対数変換した値をステップＳ６からステップＳ２１の計算に用いる。但し、対数変換以前の検査値が０．０００である場合、対数変換不可能であるため近似値として０．００１を対数変換した値を用いる。なお、最初から正規分布する検査項目のみを処理対象とする場合には、ステップＳ４の処理は不要である。
【００２７】
これは、検査値の分布が正規分布に近似する検査項目と、対数正規分布に近似する検査項目があり、本法では正規分布に近似することを前提としているため、対数正規分布に近似する検査項目の検査値を対数変換により正規分布に近似させる措置である。
ＳｑＸ＝［Σｉ（ｘｉ−ｘａｖ）３／ｎ０］／ｘσ３ （式１）
ここで、Ｘは検査項目、ｘｉ（ｉ＝１〜ｎ０）は多数の検体を測定して得られた検査項目（Ｘ）の検査値を示す。ｘａｖは検査項目（Ｘ）の平均値である。ｎ０は検査値数である。ｘσは検査項目（Ｘ）の標準偏差である。Σｉは、添え字（ｉ）についての和を意味する。
【００２８】
ステップＳ５において、検査項目（上記の１１８項目）の中から１つの検査項目（ｊとする）を選択する。
【００２９】
ステップＳ６において、ステップＳ５で指定した検査項目（ｊ）の第１採取年月日の検査値「ｘｉｊｋ（ｉ＝１〜ｎ１）」と第２採取年月日の検査値「ｙｉｊｍ（ｉ＝１〜ｎ１）」との差である検査値差「ｚｉｊｋｍ（ｉ＝１〜ｎ１）」を算出し、このｚｉｊｋｍを記録部（１３）に記録する。ここで、ｘｉｊｋは被験者（ｉ）、検査項目（ｊ）、及び第１採取年月日（ｋ）の検査値、ｙｉｊｍは被験者（ｉ）、検査項目（ｊ）、第２採取年月日（ｍ）の検査値を意味しており、ｋ＜ｍ、即ち第１採取年月日は第２採取年月日より過去の採取年月日である。である。ｚｉｊｋｍは、同一被験者（ｉ）及び同一検査項目（ｊ）における第１採取年月日（ｋ）の検査値と第２採取年月日（ｍ）の検査値との差、すなわち、ｘｉｊｋとｙｉｊｍとの差である。ｎ１は被験者数である。採取年月日が同日（ｋ＝ｍ）の場合は、時差、即ちｘｉｊｋとｙｉｊｍとの検査時刻間差で区別する。ステップＳ６の処理は、検査項目ごとに、重複しない全ての第１及び第２採取年月日の組み合わせについて行う。なお、第１の採取年月日及び第２の採取年月日は許容日差の範囲内にあることが好ましく、この許容日差は、特に限定されるものではないが、通常０日〜７日程度である。
【００３０】
ステップＳ７において、全ての検査項目についてステップＳ５及びステップＳ６の処理が終了したか否かを判断し、終了した場合、ステップＳ８に移行する。未処理の検査項目があればステップＳ５に戻る。ステップＳ５に戻った場合、未処理の検査項目の中から１つの検査項目を選択し、ステップＳ６の処理を実行する。
【００３１】
ステップＳ８において、検査項目（上記の１１８項目）の中から異なる２つの検査項目（ｊ、ｏ）を選択し、それらの検査値差「ｚｉｊｋｍ（ｉ＝１〜ｎ1）」及び「ｚｉｏｋｍ（ｉ＝１〜ｎ1）」を記録部（１３）から読み出す。ここで、ｚｉｊｋｍは、同一被験者（ｉ）及び同一検査項目（ｊ）における第１採取年月日（ｋ）の検査値と第２採取年月日（ｍ）の検査値との差であり、同様に、ｚｉｏｋｍは、同一被験者（ｉ）及び同一検査項目（ｏ）における第１採取年月日（ｋ）の検査値と第２採取年月日（ｍ）の検査値との差である。第１採取年月日（ｋ）及び第２採取年月日（ｍ）は、被験者ｉ、検査項目ｊ、及び検査項目ｏについて同一とする。なお、検査項目（ｊ、ｏ）は本発明の第１及び第２検査項目に相当する。
【００３２】
ステップＳ９において、ステップＳ８で読み出したｎ１個の検査値差（ｚｉｊｋｍ）と検査値差（ｚｉｏｋｍ）との検査値差対（ペア）「ｚｉｊｋｍ、ｚｉｏｋｍ（ｉ＝１〜ｎ1）」を作成する。なお、ステップＳ９の処理は、選択した検査項目（ｊ、ｏ）ごとに、重複しない全ての第１及び第２採取年月日の組み合わせについて行う。
【００３３】
ステップＳ１０において、ステップＳ９で作成された検査値差対「ｚｉｊｋｍ、ｚｉｏｋｍ（ｉ＝１〜ｎ1）」を用いて、相関係数（ｒ）を公知の次式２で計算する。
ｒ＝Σｉ（ｚｉｊｋｍ−ｚｊ）（ｚｉｏｋｍ−ｚｏ）／｛Σｉ（ｚｉｊｋｍ−ｚｊ）２Σｉ（ｚｉｏｋｍ−ｚｏ）２｝１／２ （式２）
ここで、ｚｊ、ｚｏは各検査値差の平均値であり、Σｉは、添え字（ｉ）についての和を意味する。
【００３４】
ステップＳ１１において、ステップＳ１０で計算された相関係数（ｒ）の絶対値が所定のしきい値Ｔｈ１（例えば０．５）よりも大きく、且つ、２つの検査項目の検査値差対の数が所定数以上（例えば５０以上、望ましくは１００以上）である場合、ステップＳ１２に移行し、そうでなければステップＳ１３に移行する。検査値差対の数についての条件は、統計的に有効な結果が得られるようにするためである。
【００３５】
ステップＳ１２において、ステップ１１で選択された２つの検査項目（ｊ、ｏ）を後述する検査過誤の可能性の判定の対象として決定し、例えば、２つの検査項目（ｊ、ｏ）の情報を対応付けて記録部（１３）に記録する。
【００３６】
ステップＳ１３において、重複せずに選択された全ての検査項目の組み合わせについて、ステップＳ８〜ステップＳ１２の処理が終了したか否かを判断し、終了した場合、ステップＳ１４に移行する。残っている組み合わせがあればステップＳ８に戻る。ステップＳ８に戻った場合、再び検査項目の中から異なる２つの検査項目を選択し、ステップＳ９〜ステップＳ１３の処理を実行する。このとき、既にステップＳ８で選択された項目と重複しないように、２つの検査項目を選択する。
【００３７】
以上のように、ステップＳ１〜ステップＳ１３によって、相関の高い２つの検査項目の組が複数決定される。以下のステップＳ１４〜ステップＳ１６の処理においては、ここで決定された組が使用される。
【００３８】
ステップＳ１４において、ステップＳ１２で決定された組の検査値差対の第２主成分を計算する。これをａｉｊｏとする。ここで、ａｉｊｏは、２つの検査項目（ｊ、ｏ）の対応する検査値差「ｚｉｊｋｍ、ｚｉｏｋｍ（ｉ＝１〜ｎ１）」に対して、公知の主成分分析を適用して決定された第２主成分（短軸成分）である。ここで、対応する検査値差（ｚｉｊｋｍ、ｚｉｏｋｍ）とは、被験者（ｉ）、第１採取年月日（ｋ）、及び第２採取年月日（ｍ）が同一の検査値を意味する。検査値差（ｚｉｊｋｍ、ｚｉｏｋｍ）のそれぞれの平均をｚｊ、ｚｏ、標準偏差をｚｊσ、ｚｏσとして、ａｉｊｏは、
相関係数（ｒ）＞０．５の場合、
ａｉｊｏ＝｛（ｚｉｊｋｍ−ｚｊ）／ｚｊσ−（ｚｉｏｋｍ−ｚｏ）／ｚｏσ｝／２１／２ で、
相関係数（ｒ）＜−０．５の場合、
ａｉｊｏ＝｛（ｚｉｊｋｍ−ｚｊ）／ｚｊσ＋（ｚｉｏｋｍ−ｚｏ）／ｚｏσ｝／２１／２ で計算される。
なお、ステップＳ１４の処理は、組として決定された検査項目（ｊ、ｏ）ごとに、重複しない全ての第１及び第２採取年月日の組み合わせについて行う。
【００３９】
ステップＳ１５において、ステップＳ１４で計算された第２主成分「ａｉｊｏ（ｉ＝１〜ｎ１）」の標準偏差を計算する。これを２つの検査項目（ｊ、ｏ）の第２主成分の標準偏差（σｊｏ）としてパラメータに設定し、記録部（１３）に記録する。
【００４０】
ステップＳ１６において、ステップＳ１２で決定された複数の組のうち、１つの組（検査項目（ｊ、ｏ））を選択して記録部（１３）から読み出し、さらに、この組に対応する、被験者（ｐ）の処理済検査値及び未処理検査値を記録部（１３）から読み出す。読み出した処理済検査値及び未処理検査値について、検査項目（ｊ、ｏ）ごとに任意に選択した第３採取年月日（ｒ）と第４採取年月日（ｓ）との検査値差（ｚｐｊｒｓ 、ｚｐｏｒｓ）を計算する。このとき、第３採取年月日は処理済検査値の採取年月日であり、第４採取年月日は未処理検査値の採取年月日であり、且つ、ｒ＜ｓである。即ち第３採取年月日は第４採取年月日より過去の採取年月日である。なお、第３の採取年月日及び第４の採取年月日は、特に限定されるものではないが、上述した第１の採取年月日及び第２の採取年月日と同じ許容日差の範囲内にあることが好ましい。この結果対（ｚｐｊｒｓ 、ｚｐｏｒｓ）に対応する検査値は、パラメータ作成用検査値からパラメータであるσｊｏが設定された後、本法によって検査過誤検出の対象となる管理対象検査値である。ステップＳ4において、正規分布に近似していると判断された検査項目は、そのままの検査値を用い、対数正規分布に近似していると判断された検査項目は対数変換された検査値を用いる。
【００４１】
ステップＳ１７において、次式３によってＳＤＩｐｊｏを計算する。
ＳＤＩｐｊｏ＝ａｐｊｏ／σｊｏ （式３）
ここで、ａｐｊｏは、ステップＳ１６において算出されたｚｐｊｒｓ 、ｚｐｏｒｓに対し、以下のように計算される。すなわち、ｚｐｊｒｓ 、ｚｐｏｒｓに対しステップＳ１４で求めた平均ｚｊ、ｚｏ、標準偏差ｚｊσ、ｚｏσを用いて、ａｐｊｏは、
相関係数（ｒ）＞０．５の場合、
ａｐｊｏ＝｛（ｚｐｊｒｓ−ｚｊ）／ｚｊσ−（ｚｐｏｒｓ−ｚｏ）／ｚｏσ｝／２１／２ で、
相関係数（ｒ）＜−０．５の場合、
ａｐｊｏ＝｛（ｚｐｊｒｓ−ｚｊ）／ｚｊσ＋（ｚｐｏｒｓ−ｚｏ）／ｚｏσ｝／２１／２ で計算される。σｊｏは、ステップＳ１５で計算された第２主成分の標準偏差である。主成分分析において、第１主成分（長軸成分）は系統誤差成分、第２主成分（短軸成分）は偶発誤差成分を表すと考えられる。なお、主成分分析は公知であるので説明を省略する。
【００４２】
ステップＳ１８において、ステップＳ１７で計算されたＳＤＩｐｊｏの絶対値が、所定のしきい値Ｔｈ２よりも大きいか否かを判断する。この条件を満たす場合、ステップＳ１９に移行し、そうでなければステップＳ２０に移行する。しきい値Ｔｈ２は、標準偏差（σｊｏ）を基準として設定することができる。しきい値Ｔｈ２は、例えば２以上、望ましくは４以上である。
【００４３】
ステップＳ１９において、ステップＳ１８の条件を満たす検査値差対（ｚｐｊｒｓ、ｚｐｏｒｓ）に対応する処理済検査値及び未処理検査値のいずれかに検査過誤の可能性があると判断し、その情報を記録する。例えば、この検査値差対（ｚｐｊｒｓ 、ｚｐｏｒｓ）に対応する２つの検査項目コード（ｊ、ｏ）、２つの採取年月日（ｒ、ｓ）、及びそれぞれの検査値を対応付けて記録部（１３）に記録する。
【００４４】
図３は、２つの検査項目（ｊ、ｏ）について、上記の処理を概念的に説明するための相関図である。図３では省略しているが、検査項目（ｊ）における第３及び第４採取年月日の検査値差と検査項目（ｏ）における第３及び第４採取年月日の検査値差との検査値差対（ペア）が複数プロットされている。楕円は、プロットされた検査値差対の分布範囲を示している。検査項目（ｊ）の検査値差と検査項目（ｏ）の検査値差は正の相関を示している。ベクトル（α）は分布範囲内の検査値差対を表すベクトルであり、ベクトル（β）は分布範囲を逸脱している検査値差対を表すベクトルである。また、ベクトル（α）、（β）は、短軸方向に平行であり、それぞれの長さは、長軸から検査値差対までの垂直距離（短軸方向の距離）である。このベクトル（α）、（β）の「長さ」が偶発誤差成分（第２主成分）である。ベクトル（β）のように偶発誤差成分が任意の管理限界（しきい値Ｔｈ２）を逸脱する場合、２つの検査値差に対応する処理済検査値と未処理検査値のいずれかに過誤があると判断できる。
【００４５】
ステップＳ２０において、全ての組について処理を終了したか否かを判断し、残っている組があればステップＳ１６に戻り、全て終了していればステップＳ２１に移行する。ステップＳ１６に戻った場合、再び１つの組を選択し、ステップＳ１７〜ステップＳ２０の処理を実行する。このとき、既にステップＳ１６で選択された組を除外して選択する。
【００４６】
ステップＳ２１において、ステップＳ１９で記録された検査項目、検査日（測定年月日）、検査値、及び関連する情報を表示部（１５）に表示する。なお、ステップＳ４において対数変換した検査値は変換前の検査値に戻して表示する。表１に、表示の一例を示す。ここで、第３採取年月日（前回日）と第４採取年月日（今回日）は西暦の年月日が８桁で表されている。
【００４７】
【表１】
【００４８】
以上によって、検査過誤検出の対象となる実検体の検査値対のうち、検査過誤に該当する可能性の高い検査値を提示することができる。
【００４９】
なお、上記処理において、検査値群が異なれば、即ち、医療機関や検査値の測定装置が異なれば、ステップＳ１２で決定される相関の高い検査項目や第２主成分の標準偏差（σｊｏ）が異なり得る。従って、ステップＳ１〜Ｓ１５の処理は、検査値群が追加、変更される毎に行うことが望ましい。その一方、特定の医療機関で同じ測定装置を使用していれば、既に決定された相関の高い２つの検査項目を用いて、ステップＳ１６〜Ｓ２１の処理を行っても、高い信頼性で検査過誤の可能性を提示することができる。また、相関係数（ｒ）の絶対値が０．６以上である組を使用すれば、医療機関によらず高い信頼性で検査過誤の可能性を提示することができる。
【００５０】
以上、本発明の実施の形態に係る臨床検査値の管理装置、及びそれを用いた管理方法について説明したが、本発明は上記に限定されない。例えば、管理装置（１）は、図１の構成から変更して構成され得る。また、図２に示したフローチャートは、種々変更されて実行され得る。
【００５１】
例えば、図１では、管理装置（１）が処理済検査値及び未処理検査値を、ネットワークを介して取得する場合を説明したが、可搬性の記録媒体（光磁気記録媒体、フラッシュメモリなど）を介して取得してもよい。
【００５２】
また、専用の管理装置として実現する代わりに、上記した機能をコンピュータに実現させるプログラムを、既存の臨床検査装置に搭載してもよい。
【００５３】
また、しきい値は任意である。例えば、Ｔｈ１として０．６以上の値を使用すれば、医療機関などに依存せずに一般的に使用できる検査項目及び採取年月日を決定できる。
【００５４】
また、Ｔｈ２として、標準偏差の４倍よりも大きい値を使用することも、標準偏差の４倍よりも小さい値を使用することもできる。但し、Ｔｈ２は、小さ過ぎると無駄な再検査が増加し、大きすぎると結果の信頼性が低下するので好ましくない。少なくとも標準偏差の２倍であることが必要である。
【００５５】
また、図２の判断ステップ（Ｓ１１、Ｓ１８）において、値を計算する度に判断する代わりに、計算結果を記録し、全ての計算が完了した後に、計算結果の全てについて、条件を満たすか否かを判断してもよい。
【００５６】
また、測定された検査値が、正規分布または対数正規分布する場合に限定されず、これらの分布に対応しない場合には、変換後の検査値が正規分布となるように変換すればよい。
【００５７】
また、上記では、未処理検査値（医師に報告されていない検査値）の検査過誤を検出する場合を説明したが、これに限定されない。将来の検査値（将来において測定された場合の検査値）の検査過誤を検出することも可能である。また、処理済検査値（医師に報告されて診療に利用された検査値）についても、同様に検査過誤を検出することが可能である。
【実施例】
【００５８】
以下に実施例を示し、本発明の特徴とするところをより一層明確にする。
【００５９】
２０１０年１月１日〜２０１０年６月３０日の間に、大分大学医学部附属病院で検査された検体の検査値に、本発明を適用した。検体の総数は２０８８２１であり、上記で示した生化学検査、血液学検査、血清学検査の合計１１８項目の検査値を用いた。
【００６０】
表２に、図２のステップＳ１〜Ｓ１３に対応する処理を行った結果、得られた相関の高い２つの検査項目（｜ｒ｜＞０．５）の組の一部を示す。全部でなく一部を示すのは相関の高い組合せ数が膨大であるためである。これらの検査項目の組の検査値を用いて、検査過誤の可能性を判断することが望ましい。なお、｜ｒ｜＞０．６の検査項目の組には、右端列に“１”が表示されている。これらの検査項目の組は、別の医療機関においても、検査過誤の判断に使用することができる。
【００６１】
【表２】
【符号の説明】
【００６２】
１ 管理装置
１１ 演算処理部
１２ 一時記憶部
１３ 記録部
１４ 操作部
１５ 表示部
１６ 入出力インタフェース部（入出力ＩＦ部）
１７ 通信インタフェース部（通信ＩＦ部）
１８ 内部バス
２ 測定装置
３ データベース
４ ネットワーク
【特許請求の範囲】
【請求項１】
臨床検査値の検査過誤を検出可能な装置であって、
記録部と演算処理部とを備え、
前記記録部に、複数の処理済検査値及び未処理検査値が被験者、検査項目、及び採取年月日と対応させて記録されており、
前記演算処理部は、
第１又は第２検査項目、且つ第１又は第２採取年月日に対応する前記複数の処理済検査値を、前記記録部から読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目各々における前記第１及び第２採取年月日の前記パラメータ作成用検査値の差を前記被験者ごとに計算し、前記パラメータ作成用検査値の差の相関係数を計算し、前記相関係数の絶対値が第１のしきい値よりも大きい場合、前記第１及び第２検査項目を検査過誤の判断に使用可能な組として決定し、
前記組として決定された第１及び第２検査項目の第２主成分を、前記パラメータ作成用検査値の差を用いて計算し、前記第２主成分の標準偏差を計算し、
前記組として決定された第１及び第２検査項目に該当する前記処理済検査値及び未処理検査値を前記記録部から読み出し、前記第１及び第２検査項目各々における第３採取年月日の前記処理済検査値及び第４採取年月日の前記未処理検査値の差を前記被験者ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差の第２主成分を計算し、前記処理済検査値及び未処理検査値の差の第２主成分を前記標準偏差で除した値の絶対値が、第２のしきい値よりも大きい場合、前記処理済検査値及び未処理検査値の差の第２主成分に対応する前記処理済検査値及び未処理検査値のいずれかに検査過誤があると判断し、
前記第１のしきい値が０．５以上であり、
前記第２のしきい値が２以上である
ことを特徴とする臨床検査値の管理装置。
【請求項２】
前記処理済検査値を読み出してパラメータ作成用検査値群とすることが、
前記演算処理部が、読み出された前記複数の処理済検査値を前記検査項目ごとにトランケーションして、前記各処理済検査値をパラメータ作成用検査値とするか否かを判断すること、及び、
歪度を基に前記検査項目ごとの前記パラメータ作成用検査値群の分布形を判断し、該判断結果に応じて前記各パラメータ作成用検査値を変換する、請求項１に記載の臨床検査値の管理装置。
【請求項３】
前記第１採取年月日と前記第２採取年月日との差、及び前記第３採取年月日と前記第４採取年月日との差がともに許容日差の範囲内にある、請求項１又は２に記載の管理装置。
【請求項４】
臨床検査値の検査過誤を検出可能な方法であって、
複数の処理済検査値及び未処理検査値が、被験者、検査項目、及び採取年月日と対応させて記録された記録部から、第１又は第２検査項目、且つ第１又は第２採取年月日に対応する前記複数の処理済検査値を読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目各々における前記第１及び第２採取年月日の前記パラメータ作成用検査値の差を前記被験者ごとに計算し、前記パラメータ作成用検査値の差の相関係数を計算する第１ステップと、
前記相関係数の絶対値が第１のしきい値よりも大きいか否かを判断し、大きい場合、前記第１及び第２検査項目を検査過誤の判断に使用可能な組として決定する第２ステップと、
前記組として決定された第１及び第２検査項目の第２主成分を、前記パラメータ作成用検査値の差を用いて計算する第３ステップと、
前記第２主成分の標準偏差を計算する第４ステップと、
前記組として決定された第１及び第２検査項目に該当する前記処理済検査値及び未処理検査値を前記記録部から読み出し、前記第１及び第２検査項目各々における第３採取年月日の前記処理済検査値及び第４採取年月日の前記処理済検査値及び未処理検査値の差を前記被験者ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差の第２主成分を計算する第５ステップと、
前記処理済検査値及び未処理検査値の差の第２主成分を前記標準偏差で除した値の絶対値が、第２のしきい値よりも大きい場合、前記処理済検査値及び未処理検査値の差の第２主成分に対応する処理済検査値及び前記未処理検査値のいずれかに検査過誤があると判断する第６ステップと、を含み、
前記第１のしきい値が０．５以上であり、
前記第２のしきい値が２以上である
ことを特徴とする臨床検査値の管理方法。
【請求項５】
臨床検査値の検査過誤を検出可能なコンピュータプログラムであって、
コンピュータに、
複数の処理済検査値及び未処理検査値が、被験者、検査項目、及び採取年月日と対応させて記録された記録部から、第１又は第２検査項目、且つ第１又は第２採取年月日に対応する前記複数の処理済検査値を読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目各々における前記第１及び第２採取年月日の前記パラメータ作成用検査値の差を前記被験者ごとに計算し、前記パラメータ作成用検査値の差の相関係数を計算する第１の機能と、
前記相関係数の絶対値が第１のしきい値よりも大きいか否かを判断し、大きい場合、前記第１及び第２検査項目を検査過誤の判断に使用可能な組として決定する第２の機能と、
前記組として決定された第１及び第２検査項目の第２主成分を、前記処理対象検査値の差を用いて計算する第３機能と、
前記第２主成分の標準偏差を計算する第４機能と、
前記組として決定された第１及び第２検査項目に該当する前記処理済検査値及び未処理検査値を前記記録部から読み出し、前記第１及び第２検査項目各々における第３採取年月日の前記処理済検査値及び第４採取年月日の前記処理済検査値及び未処理検査値の差を前記被験者ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差の第２主成分を計算する第５の機能と、
前記処理済検査値及び未処理検査値の差の第２主成分を前記標準偏差で除した値の絶対値が、第２のしきい値よりも大きい場合、前記処理済検査値及び未処理検査値の差の第２主成分に対応する処理済検査値及び前記未処理検査値のいずれかに検査過誤があると判断する第６の機能と、を実現させ、
前記第１のしきい値が０．５以上であり、
前記第２のしきい値が２以上である
ことを特徴とする臨床検査値の管理プログラム。
【請求項１】
臨床検査値の検査過誤を検出可能な装置であって、
記録部と演算処理部とを備え、
前記記録部に、複数の処理済検査値及び未処理検査値が被験者、検査項目、及び採取年月日と対応させて記録されており、
前記演算処理部は、
第１又は第２検査項目、且つ第１又は第２採取年月日に対応する前記複数の処理済検査値を、前記記録部から読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目各々における前記第１及び第２採取年月日の前記パラメータ作成用検査値の差を前記被験者ごとに計算し、前記パラメータ作成用検査値の差の相関係数を計算し、前記相関係数の絶対値が第１のしきい値よりも大きい場合、前記第１及び第２検査項目を検査過誤の判断に使用可能な組として決定し、
前記組として決定された第１及び第２検査項目の第２主成分を、前記パラメータ作成用検査値の差を用いて計算し、前記第２主成分の標準偏差を計算し、
前記組として決定された第１及び第２検査項目に該当する前記処理済検査値及び未処理検査値を前記記録部から読み出し、前記第１及び第２検査項目各々における第３採取年月日の前記処理済検査値及び第４採取年月日の前記未処理検査値の差を前記被験者ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差の第２主成分を計算し、前記処理済検査値及び未処理検査値の差の第２主成分を前記標準偏差で除した値の絶対値が、第２のしきい値よりも大きい場合、前記処理済検査値及び未処理検査値の差の第２主成分に対応する前記処理済検査値及び未処理検査値のいずれかに検査過誤があると判断し、
前記第１のしきい値が０．５以上であり、
前記第２のしきい値が２以上である
ことを特徴とする臨床検査値の管理装置。
【請求項２】
前記処理済検査値を読み出してパラメータ作成用検査値群とすることが、
前記演算処理部が、読み出された前記複数の処理済検査値を前記検査項目ごとにトランケーションして、前記各処理済検査値をパラメータ作成用検査値とするか否かを判断すること、及び、
歪度を基に前記検査項目ごとの前記パラメータ作成用検査値群の分布形を判断し、該判断結果に応じて前記各パラメータ作成用検査値を変換する、請求項１に記載の臨床検査値の管理装置。
【請求項３】
前記第１採取年月日と前記第２採取年月日との差、及び前記第３採取年月日と前記第４採取年月日との差がともに許容日差の範囲内にある、請求項１又は２に記載の管理装置。
【請求項４】
臨床検査値の検査過誤を検出可能な方法であって、
複数の処理済検査値及び未処理検査値が、被験者、検査項目、及び採取年月日と対応させて記録された記録部から、第１又は第２検査項目、且つ第１又は第２採取年月日に対応する前記複数の処理済検査値を読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目各々における前記第１及び第２採取年月日の前記パラメータ作成用検査値の差を前記被験者ごとに計算し、前記パラメータ作成用検査値の差の相関係数を計算する第１ステップと、
前記相関係数の絶対値が第１のしきい値よりも大きいか否かを判断し、大きい場合、前記第１及び第２検査項目を検査過誤の判断に使用可能な組として決定する第２ステップと、
前記組として決定された第１及び第２検査項目の第２主成分を、前記パラメータ作成用検査値の差を用いて計算する第３ステップと、
前記第２主成分の標準偏差を計算する第４ステップと、
前記組として決定された第１及び第２検査項目に該当する前記処理済検査値及び未処理検査値を前記記録部から読み出し、前記第１及び第２検査項目各々における第３採取年月日の前記処理済検査値及び第４採取年月日の前記処理済検査値及び未処理検査値の差を前記被験者ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差の第２主成分を計算する第５ステップと、
前記処理済検査値及び未処理検査値の差の第２主成分を前記標準偏差で除した値の絶対値が、第２のしきい値よりも大きい場合、前記処理済検査値及び未処理検査値の差の第２主成分に対応する処理済検査値及び前記未処理検査値のいずれかに検査過誤があると判断する第６ステップと、を含み、
前記第１のしきい値が０．５以上であり、
前記第２のしきい値が２以上である
ことを特徴とする臨床検査値の管理方法。
【請求項５】
臨床検査値の検査過誤を検出可能なコンピュータプログラムであって、
コンピュータに、
複数の処理済検査値及び未処理検査値が、被験者、検査項目、及び採取年月日と対応させて記録された記録部から、第１又は第２検査項目、且つ第１又は第２採取年月日に対応する前記複数の処理済検査値を読み出してパラメータ作成用検査値群とし、前記第１及び第２検査項目各々における前記第１及び第２採取年月日の前記パラメータ作成用検査値の差を前記被験者ごとに計算し、前記パラメータ作成用検査値の差の相関係数を計算する第１の機能と、
前記相関係数の絶対値が第１のしきい値よりも大きいか否かを判断し、大きい場合、前記第１及び第２検査項目を検査過誤の判断に使用可能な組として決定する第２の機能と、
前記組として決定された第１及び第２検査項目の第２主成分を、前記処理対象検査値の差を用いて計算する第３機能と、
前記第２主成分の標準偏差を計算する第４機能と、
前記組として決定された第１及び第２検査項目に該当する前記処理済検査値及び未処理検査値を前記記録部から読み出し、前記第１及び第２検査項目各々における第３採取年月日の前記処理済検査値及び第４採取年月日の前記処理済検査値及び未処理検査値の差を前記被験者ごとに計算し、主成分分析によって前記処理済検査値及び未処理検査値の差の第２主成分を計算する第５の機能と、
前記処理済検査値及び未処理検査値の差の第２主成分を前記標準偏差で除した値の絶対値が、第２のしきい値よりも大きい場合、前記処理済検査値及び未処理検査値の差の第２主成分に対応する処理済検査値及び前記未処理検査値のいずれかに検査過誤があると判断する第６の機能と、を実現させ、
前記第１のしきい値が０．５以上であり、
前記第２のしきい値が２以上である
ことを特徴とする臨床検査値の管理プログラム。
【図１】
【図２】
【図３】
【図２】
【図３】
【公開番号】特開２０１２−６４２０３（Ｐ２０１２−６４２０３Ａ）
【公開日】平成２４年３月２９日（２０１２．３．２９）
【国際特許分類】
【出願番号】特願２０１１−１７５０７７（Ｐ２０１１−１７５０７７）
【出願日】平成２３年８月１０日（２０１１．８．１０）
【出願人】（３０４０２８７２６）国立大学法人 大分大学 (181)
【Ｆターム（参考）】
【公開日】平成２４年３月２９日（２０１２．３．２９）
【国際特許分類】
【出願日】平成２３年８月１０日（２０１１．８．１０）
【出願人】（３０４０２８７２６）国立大学法人 大分大学 (181)
【Ｆターム（参考）】
[ Back to top ]