JP6272065B2 - NASH diagnostic program, NASH diagnostic device operating method, and NASH diagnostic device - Google Patents

NASH diagnostic program, NASH diagnostic device operating method, and NASH diagnostic device Download PDF

Info

Publication number
JP6272065B2
JP6272065B2 JP2014022570A JP2014022570A JP6272065B2 JP 6272065 B2 JP6272065 B2 JP 6272065B2 JP 2014022570 A JP2014022570 A JP 2014022570A JP 2014022570 A JP2014022570 A JP 2014022570A JP 6272065 B2 JP6272065 B2 JP 6272065B2
Authority
JP
Japan
Prior art keywords
nash
bmipp
patient
curve
probability value
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.)
Active
Application number
JP2014022570A
Other languages
Japanese (ja)
Other versions
JP2015148555A (en
Inventor
和雄 窪田
和雄 窪田
亮吾 南本
亮吾 南本
雄一 野崎
雄一 野崎
明希男 林
明希男 林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nihon Medi Physics Co Ltd
National Center for Global Health and Medicine
Original Assignee
Nihon Medi Physics Co Ltd
National Center for Global Health and Medicine
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nihon Medi Physics Co Ltd, National Center for Global Health and Medicine filed Critical Nihon Medi Physics Co Ltd
Priority to JP2014022570A priority Critical patent/JP6272065B2/en
Publication of JP2015148555A publication Critical patent/JP2015148555A/en
Application granted granted Critical
Publication of JP6272065B2 publication Critical patent/JP6272065B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Nuclear Medicine (AREA)

Description

本発明は、非アルコール性脂肪性肝疾患(NAFLD)のうち、非アルコール性脂肪肝炎(NASH)を診断する方法に関する。   The present invention relates to a method for diagnosing nonalcoholic steatohepatitis (NASH) among nonalcoholic fatty liver disease (NAFLD).

近年、NAFLD患者が急増している。NAFLDの病態の一つとして、NASHがある。NASHは、肝細胞への中性脂肪の沈着に加えて、炎症性細胞浸潤、肝細胞の変性、壊死、繊維化を伴うものである。NASHは、NAFLDの中でも、肝細胞への脂肪沈着のみで他の病変がない単純性脂肪肝(SS)とは区別される。NASH患者のうちの所定割合の患者が、肝硬変、肝細胞癌に移行することが知られているので、NASHとSSとを鑑別し、NASHを早期に発見することが望ましい。   In recent years, the number of NAFLD patients has increased rapidly. One of the pathologies of NAFLD is NASH. NASH is accompanied by inflammatory cell infiltration, hepatocyte degeneration, necrosis, and fibrosis in addition to the deposition of neutral fat in hepatocytes. NASH is distinguished from simple fatty liver (SS) among NAFLDs, in which only fat is deposited on hepatocytes and there are no other lesions. Since it is known that a predetermined percentage of NASH patients shift to cirrhosis and hepatocellular carcinoma, it is desirable to distinguish NASH from SS and detect NASH at an early stage.

従来、NASHの確定診断には肝生検が実施されてきたが、生検は侵襲的であり患者に身体的な負担を強いることから、非侵襲的な方法による診断が望まれてきた。この目的で、これまで、CT、超音波、MRI、トランジェント・エラストグラフィー等の様々な画像診断法や、in vitroバイオマーカー等が検討されてきたが、十分な診断能を有するNASHの診断方法は、これまでに実用化されていなかった。   Conventionally, a liver biopsy has been performed for a definitive diagnosis of NASH. However, since the biopsy is invasive and imposes a physical burden on the patient, a diagnosis by a non-invasive method has been desired. For this purpose, various diagnostic imaging methods such as CT, ultrasound, MRI, transient elastography, in vitro biomarkers, etc. have been studied so far. Until now, it has not been put into practical use.

本発明の発明者らは、心疾患の脂肪酸代謝シンチグラフィ用に市販されている放射性医薬品である123I−BMIPPを用いて、NAFLD患者の進行度評価に関する可能性について検討し、進行したNASH患者では、肝への脂肪酸の集積時間が延長すること、および集積度が低下することを見出した(非特許文献1)。 The inventors of the present invention examined the possibility of evaluating the progression of NAFLD patients using 123 I-BMIPP, a radiopharmaceutical marketed for fatty acid metabolism scintigraphy of heart disease, and advanced NASH patients Then, it has been found that the accumulation time of fatty acids in the liver is prolonged and the accumulation degree is lowered (Non-patent Document 1).

特表2008−538893号公報Special table 2008-538893 gazette

南本亮吾ら「脂肪酸代謝シンチグラフィによるNAFLD患者の進行度診断に関する検討」核医学、第49巻第3号Ryogo Minamimoto et al. “Study on diagnosis of progression of NAFLD patients by fatty acid metabolism scintigraphy” Nuclear Medicine, Vol. 49, No. 3

特許文献1には、123I−BMIPPを用いて体内の脂肪酸β酸化能の評価を行えることが記載されている。しかしながら、特許文献1にはNASHとSSとを鑑別するための具体的な手法は開示されていない。 Patent Document 1 describes that evaluation of fatty acid β-oxidation ability in the body can be performed using 123 I-BMIPP. However, Patent Document 1 does not disclose a specific method for distinguishing NASH from SS.

上記非特許文献1に開示された本発明者らの研究は、NAFLD患者の進行度評価に関する可能性を示すに止まっていた。臨床の現場においては、現在、一千万人以上いると言われている脂肪肝患者の中から、肝生検を要する進行性のNASH患者を効果的に抽出することが急務である。特に、初期のNASHとSSと鑑別するのは困難であり、これを鑑別する方法が求められている。   The present inventors' research disclosed in the above-mentioned Non-Patent Document 1 has only shown the possibility of evaluating the degree of progression of NAFLD patients. In clinical practice, there is an urgent need to effectively extract progressive NASH patients who require liver biopsy from among the fatty liver patients who are said to have more than 10 million people. In particular, it is difficult to discriminate between initial NASH and SS, and a method for discriminating this is required.

そこで、本発明は、NASH、特に早期のNASHとSSとを非侵襲的に診断することができる方法を提供することを目的とする。   Therefore, an object of the present invention is to provide a method capable of noninvasively diagnosing NASH, particularly early NASH and SS.

本発明のプログラムは、検査対象の患者に対する123I−BMIPPを用いたシンチグラフィによって得られた時間放射能曲線のデータに基づいて、検査対象の患者がNASHであるか否かを診断するためのプログラムであって、コンピュータに、NASHであることが既知の少なくとも一の被検者及び非NASHであることが既知の少なくとも一の被検者に対して、それぞれ123I−BMIPPを用いたシンチグラフィを行って得られた時間放射能曲線の各データについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、前記被検者のフィッティングカーブをロジスティック回帰分析することにより、前記説明変数に基づいて患者がNASHである確率値を求める回帰式を生成するステップと、前記検査対象の患者の時間放射能曲線のデータについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、前記フィッティングカーブの説明変数を前記回帰式に代入して、前記検査対象の患者がNASHである確率値を算出するステップと、前記確率値を出力するステップとを実行させる構成を有する。 The program of the present invention is for diagnosing whether or not a patient to be examined is NASH based on data of a time activity curve obtained by scintigraphy using 123 I-BMIPP for a patient to be examined. Scintigraphy using 123 I-BMIPP for at least one subject known to be NASH and at least one subject known to be non-NASH to a computer for each data of the obtained time activity curves performed to a height of 123 I-BMIPP incorporation into liver, the inflow rate of 123 I-BMIPP, the explanatory variables and exit velocity of 123 I-BMIPP The step of obtaining a fitting curve, and logistic regression analysis of the fitting curve of the subject, A step of generating a regression equation for obtaining a probability value that the patient is NASH based on the number, and the time activity curve data of the patient under examination, the height of 123 I-BMIPP uptake into the liver, and 123 I A step of obtaining a fitting curve having the inflow velocity of BMIPP and the outflow velocity of 123 I-BMIPP as explanatory variables, and substituting the explanatory variables of the fitting curve into the regression equation, so that the patient to be examined is NASH It has the structure which performs the step which calculates a certain probability value, and the step which outputs the said probability value.

本発明者らは、上記研究成果により得られた知見に基づいてさらなる実験、研究を重ね、肝シンチグラフィの結果得られる時間放射能曲線が、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブによって適切に表すことができることを見出した。このフィッティングカーブにフィッティングを行うことによって求めた説明変数と、NASH/非NASHのデータに基づいてロジスティック回帰分析を行うことにより、説明変数に基づいてNASHである確率値を求める回帰式を生成することにより、この回帰式を使って検査対象の患者がNASHである確率値を求めることができる本発明を完成させた。フィッティングカーブfは、具体的には、123I−BMIPPの取り込みの高さをA、123I−BMIPPの流出速度をα、123I−BMIPPの流入速度をβとして、f=Ae-βt−Ae-αtが好適である。 The present inventors have repeated further experiments and studies based on the findings obtained from the above research results, and the time activity curve obtained as a result of liver scintigraphy shows that the uptake of 123 I-BMIPP into the liver is high. It has been found that the inflow rate of 123 I-BMIPP and the outflow rate of 123 I-BMIPP can be appropriately represented by a fitting curve having explanatory variables. By generating logistic regression analysis based on the explanatory variables obtained by fitting the fitting curve and the NASH / non-NASH data, a regression equation for obtaining a probability value of NASH based on the explanatory variables is generated. Thus, the present invention has been completed in which the regression value can be used to obtain the probability value that the patient to be examined is NASH. Fitting curve f is specifically the exit velocity of the height of the 123 I-BMIPP uptake A, 123 I-BMIPP α, as the inlet velocity of 123 I-BMIPP β, f = Ae -βt -Ae -αt is preferred.

また、本発明のプログラムは、前記コンピュータに、前記回帰式を用いて前記被検者がNASHである確率値を求め、NASHと非NASHとを鑑別するための閾値を変数として生成されるROC曲線に基づいてNASHと非NASHとを鑑別するための閾値を決定するステップと、前記検査対象の患者がNASHである確率値と前記閾値とに基づいて、検査対象の患者がNASHか非NASHかを鑑別するステップとを実行させてもよい。この構成により、適切な閾値によって、検査対象の患者がNASHか非NASHを鑑別することができる。   In addition, the program of the present invention obtains a probability value that the subject is NASH using the regression equation, and generates an ROC curve generated using the threshold value for distinguishing NASH from non-NASH as a variable. Determining a threshold value for distinguishing NASH from non-NASH based on the step, and determining whether the patient to be tested is NASH or non-NASH based on the probability value that the patient to be examined is NASH and the threshold value. The step of discrimination may be executed. With this configuration, the patient to be examined can be identified as NASH or non-NASH with an appropriate threshold.

本発明の別の態様のプログラムは、上記した方法によって求めた、患者がNASHである確率値を求める回帰式をあらかじめ記憶した記憶部から、前記回帰式を読み出すステップと、前記フィッティングカーブの説明変数を前記回帰式に代入して、前記検査対象の患者がNASHである確率値を算出するステップと、前記確率値を出力するステップとを実行させる。この構成によっても、上記と同様に、検査対象の患者がNASHである確率値を求めることができる。   The program according to another aspect of the present invention includes a step of reading the regression equation from a storage unit that stores in advance a regression equation for obtaining a probability value that the patient is NASH, which is obtained by the above method, and an explanatory variable of the fitting curve. Is substituted into the regression equation, and a step of calculating a probability value that the patient to be examined is NASH and a step of outputting the probability value are executed. Also with this configuration, the probability value that the patient to be examined is NASH can be obtained as described above.

本発明のNASH診断方法は、NASHであることが既知の少なくとも一の被検者及び非NASHであることが既知の少なくとも一の患者に対して、123I−BMIPPを用いたシンチグラフィを行って、時間放射能曲線のデータをそれぞれ取得するステップと、前記被検者の前記時間放射能曲線の各データについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、前記被検者のフィッティングカーブをロジスティック回帰分析することにより、前記説明変数に基づいて患者がNASHである確率値を求める回帰式を生成するステップと、検査対象の患者に対して、123I−BMIPPを用いたシンチグラフィを行って時間放射能曲線のデータを取得するステップと、前記検査対象の患者の時間放射能曲線のデータについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、前記フィッティングカーブの説明変数を前記回帰式に代入して、前記検査対象の患者がNASHである確率値を算出するステップと、前記確率値を出力するステップとを備える。 The NASH diagnostic method of the present invention comprises performing scintigraphy using 123 I-BMIPP on at least one subject known to be NASH and at least one patient known to be non-NASH. Each of the time activity curve data, and for each time activity curve data of the subject, the height of 123 I-BMIPP uptake into the liver, the inflow rate of 123 I-BMIPP, 123, obtaining a fitting curve having an outflow rate of I-BMIPP as an explanatory variable, and logistic regression analysis of the fitting curve of the subject, thereby calculating a probability value that the patient is NASH based on the explanatory variable. and generating a regression formula for determining, for a patient to be examined, the scintigraphy using 123 I-BMIPP Acquiring data in the time activity curves I, the data of the time-activity curve of patients of the test object, the height of the 123 I-BMIPP incorporation into liver, the inflow rate of 123 I-BMIPP, 123 A step of obtaining a fitting curve having an outflow velocity of I-BMIPP as an explanatory variable, and a step of calculating a probability value that the patient to be examined is NASH by substituting the explanatory variable of the fitting curve into the regression equation. And outputting the probability value.

本発明のNASH診断装置は、NASHであることが既知の少なくとも一の被検者及びまたは非NASHであることが既知の少なくとも一の被検者に対して、それぞれ123I−BMIPPを用いたシンチグラフィを行って得られた時間放射能曲線の各データについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求め、前記被検者のフィッティングカーブをロジスティック回帰分析することによって生成された、前記説明変数に基づいて患者がNASHである確率値を求める回帰式を記憶した記憶部と、検査対象の患者に対して、123I−BMIPPを用いたシンチグラフィを行って得られた時間放射能曲線のデータを取得するデータ取得部と、前記検査対象の患者の時間放射能曲線のデータについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求め、前記フィッティングカーブの説明変数を前記回帰式に代入して、前記検査対象の患者がNASHである確率値を算出する演算制御部と、前記確率値を出力する出力部とを備える。 The NASH diagnostic device of the present invention is a scintigraphy using 123 I-BMIPP for at least one subject known to be NASH and / or at least one subject known to be non-NASH. for each data of the obtained time activity curves performing photography, and height of 123 I-BMIPP incorporation into liver, 123 and inflow rate of the I-BMIPP, and the explanatory variables 123 I-BMIPP exit velocity of A storage unit storing a regression equation for obtaining a probability value that the patient is NASH based on the explanatory variable, which is generated by logistic regression analysis of the fitting curve of the subject, and a test object de the relative of the patient, to acquire data of the obtained time activity curves performed scintigraphy using 123 I-BMIPP A data acquisition unit, the data of the time-activity curve of patients of the test object, the height of the 123 I-BMIPP incorporation into liver, the inflow rate of 123 I-BMIPP, the exit velocity of 123 I-BMIPP A calculation curve that obtains a fitting curve as an explanatory variable, substitutes the explanatory variable of the fitting curve into the regression equation, calculates a probability value that the patient to be examined is NASH, and an output that outputs the probability value A part.

本発明の放射性医薬は、肝シンチグラフィにより早期NASHと単純脂肪肝との鑑別に用いられる、123I−BMIPPを含有する。従来、123I−BMIPPを肝シンチグラフィにより早期NASHと単純脂肪肝との鑑別に用いるという用途は知られていなかった。本発明は、非侵襲的な早期NASHの鑑別という当業者が想到し得なかった123I−BMIPPの用途を提供するものである。なお、ここでいう「早期NASH」とは、Brunt病理分類でStage1−2のものをいう。 The radiopharmaceutical of the present invention contains 123 I-BMIPP, which is used for differentiation between early NASH and simple fatty liver by liver scintigraphy. Conventionally, the use of 123 I-BMIPP for differentiation between early NASH and simple fatty liver by liver scintigraphy has not been known. The present invention provides a use of 123 I-BMIPP that could not be conceived by those skilled in the art of non-invasive early NASH discrimination. The term “early NASH” used here refers to a stage 1-2 in the Brunt pathological classification.

本発明によれば、123I−BMIPPを用いた肝シンチグラフィにより、検査対象の患者がNASHであるか否かを適切に鑑別することができる。 According to the present invention, it is possible to appropriately distinguish whether or not a patient to be examined is NASH by liver scintigraphy using 123 I-BMIPP.

本実施の形態のNASH診断装置の構成を示す図である。It is a figure which shows the structure of the NASH diagnostic apparatus of this Embodiment. NASH診断プログラムの構成を示す図である。It is a figure which shows the structure of a NASH diagnostic program. (a)NASH診断装置による回帰式の生成、及びカットオフ値の決定の動作を示すフローチャートである。(b)NASH診断装置によりNASHか否かの鑑別を行う動作を示すフローチャートである。(A) It is a flowchart which shows the operation | movement of the production | generation of the regression type by a NASH diagnostic apparatus, and the determination of a cutoff value. (B) It is a flowchart which shows the operation | movement which discriminate | determines whether it is NASH by a NASH diagnostic apparatus. (a)単純性脂肪肝(SS)の時間放射能曲線を示すグラフである。(b)早期(Brunt病理分類Stage1−2)NASH(eNASH)の時間放射能曲線を示すグラフである。(c)進行(Stage3−4)NASH(aNASH)の時間放射能曲線を示すグラフである。(A) It is a graph which shows the time activity curve of simple fatty liver (SS). (B) It is a graph which shows the time radioactivity curve of early (Brunt pathological classification Stage1-2) NASH (eNASH). (C) It is a graph which shows the time radioactivity curve of progress (Stage3-4) NASH (aNASH).

以下、本発明の実施の形態のNASH診断装置、NASH診断方法およびNASH診断プログラムについて、図面を参照しながら説明する。   Hereinafter, a NASH diagnostic apparatus, a NASH diagnostic method, and a NASH diagnostic program according to embodiments of the present invention will be described with reference to the drawings.

図1は、実施の形態のNASH診断装置1の構成を示す図である。NASH診断装置1は、CPU、RAM、ROM、ハードディスク、モニタ、キーボード、通信インターフェース等を備えたコンピュータによって構成される。NASH診断装置1は、機能的には、データ取得部2と、演算制御部3と、記憶部4と、出力部5とを有している。   FIG. 1 is a diagram illustrating a configuration of a NASH diagnostic apparatus 1 according to the embodiment. The NASH diagnostic apparatus 1 is configured by a computer including a CPU, RAM, ROM, hard disk, monitor, keyboard, communication interface, and the like. The NASH diagnostic device 1 functionally includes a data acquisition unit 2, a calculation control unit 3, a storage unit 4, and an output unit 5.

データ取得部2は、123I−BMIPPを用いたシンチグラフィのデータを取得する機能を有する。本実施の形態においては、データ取得部2は、SPECT装置6と接続されており、SPECT装置6にて撮像したSPECTデータを取得する。なお、データ取得部2は、SPECT装置6と直接接続されている必要はない。他の病院や研究所などで取得したSPECTデータを、ネットワークを介して受信する等する構成とすることも可能である。 The data acquisition unit 2 has a function of acquiring scintigraphic data using 123 I-BMIPP. In the present embodiment, the data acquisition unit 2 is connected to the SPECT device 6 and acquires SPECT data captured by the SPECT device 6. The data acquisition unit 2 does not need to be directly connected to the SPECT device 6. It is also possible to adopt a configuration in which SPECT data acquired at another hospital or laboratory is received via a network.

演算制御部3は、データ取得部2にて取得したSPECTデータを用いて、検査対象の患者がNASHである確率値を求めたり、NASHであるか否かを鑑別する機能を有する。演算制御部3が行う詳しい処理については後述する。   The arithmetic control unit 3 has a function of obtaining a probability value that the patient to be inspected is NASH using the SPECT data acquired by the data acquisition unit 2 and discriminating whether or not the patient is NASH. Detailed processing performed by the arithmetic control unit 3 will be described later.

記憶部4は、NASHか非NASHかが既知の患者のSPECTデータに基づいて、演算制御部3にて生成されるNASHの確率値を求めるための回帰式やNASHを鑑別するための閾値(以下「カットオフ値」という)を記憶する。また、記憶部4には、コンピュータをNASH診断装置1として機能させるためのプログラムが記憶されている。   Based on SPECT data of a patient whose NASH or non-NASH is known, the storage unit 4 is a regression equation for determining the probability value of NASH generated by the arithmetic control unit 3 or a threshold value for distinguishing NASH (hereinafter referred to as NASH). (Referred to as “cut-off value”). The storage unit 4 stores a program for causing the computer to function as the NASH diagnostic device 1.

図2は、NASH診断プログラム10の構成を示す図である。NASH診断プログラムは、フィッティングカーブ生成モジュール11と、回帰式生成モジュール12と、NASH確率値算出モジュール13と、カットオフ値決定モジュール14と、NASH鑑別モジュール15と、出力モジュール16とを有している。CPUがNASH診断プログラム10を読み出して実行することにより、上述した演算制御部3および出力部5の機能が実現される。   FIG. 2 is a diagram showing the configuration of the NASH diagnostic program 10. The NASH diagnostic program includes a fitting curve generation module 11, a regression equation generation module 12, a NASH probability value calculation module 13, a cutoff value determination module 14, an NASH discrimination module 15, and an output module 16. . When the CPU reads and executes the NASH diagnostic program 10, the functions of the arithmetic control unit 3 and the output unit 5 described above are realized.

次に、演算制御部3の機能について説明する。演算制御部の機能は、主に(1)NASHである確率値を求めるための回帰式の生成やカットオフ値の決定を行う機能と、(2)検査対象の患者の123I−BMIPPを用いたシンチグラフィのデータを回帰式に適用して、NASHか否かを鑑別する機能を有している。以下、順に説明する。 Next, functions of the arithmetic control unit 3 will be described. The functions of the arithmetic control unit mainly use (1) a function for generating a regression equation for determining a probability value of NASH and determining a cutoff value, and (2) 123 I-BMIPP of the patient to be examined. It has a function of discriminating whether or not it is NASH by applying the scintigraphy data to the regression equation. Hereinafter, it demonstrates in order.

[回帰式等を生成する機能]
演算制御部3の処理に先立って、データ取得部2が、NASHであることが既知の少なくとも一の被検者、及び、非NASHであることが既知の少なくとも一の被検者の123I−BMIPPを用いたSPECTデータをそれぞれ取得する。このSPECTデータは、早期NASHであることが既知の少なくとも一の被検者、及び、進行NASHであることが既知の少なくとも一の被検者、並びに、単純脂肪肝であることが既知の少なくとも一の被検者から収集されたものであることが好ましい。「早期NASH」とは、Brunt病理分類でStage1−2のものをいい、「進行NASH」とは、Brunt病理分類でStage3−4のものをいう。これらデータは、回帰式等を生成するための学習データである。取得した学習データは、記憶部4に記憶しておいてもよい。
[Function to generate regression equations]
Prior to the processing of the arithmetic control unit 3, the data acquisition unit 2 has 123 I− of at least one subject known to be NASH and at least one subject known to be non-NASH. SPECT data using BMIPP is acquired. The SPECT data includes at least one subject known to be early NASH, at least one subject known to be advanced NASH, and at least one subject known to be simple fatty liver. It is preferably collected from the subject. “Early NASH” refers to that of Stage 1-2 in the Brunt pathology classification, and “advanced NASH” refers to that of Stage 3-4 in the Brunt pathology classification. These data are learning data for generating regression equations and the like. The acquired learning data may be stored in the storage unit 4.

演算制御部3は、データ取得部2にて取得した複数の患者の123I−BMIPPを用いたSPECTデータの時間放射能曲線のそれぞれについて、フィッティングカーブを生成する。本実施の形態では、演算制御部3は、フィッティングカーブのモデル式として、f=Ae-βt−Ae-αtを用いる。この式において、A、α、βが説明変数であり、フィッティングカーブの生成は、これらの説明変数の値を求めることと同義である。 The arithmetic control unit 3 generates a fitting curve for each of the time activity curves of SPECT data using 123 I-BMIPP of a plurality of patients acquired by the data acquisition unit 2. In the present embodiment, the calculation control unit 3 uses f = Ae −βt −Ae −αt as a model equation for the fitting curve. In this equation, A, α, and β are explanatory variables, and generation of a fitting curve is synonymous with obtaining values of these explanatory variables.

演算制御部3は、非線形最小二乗法を用いて、患者の時間放射能曲線のフィッティングカーブを生成する。本実施の形態では、非線形最小二乗法を用いているが、フィッティングカーブを生成する手法は、非線形最小二乗法に限らず、他の方法を用いてもよい。この機能は、フィッティングカーブ生成モジュール11によって実現される。   The arithmetic control unit 3 generates a fitting curve of the time activity curve of the patient using a non-linear least square method. In this embodiment, the nonlinear least square method is used. However, the method for generating the fitting curve is not limited to the nonlinear least square method, and other methods may be used. This function is realized by the fitting curve generation module 11.

演算制御部3は、学習データに係る患者の説明変数A、α、βと、NASH/非NASHのデータのロジスティック回帰分析を行い、説明変数A、α、βに基づいて、患者がNASHである確率値を求める回帰式を生成する。演算制御部3は、生成した回帰式を記憶部4に記憶する。ロジスティック回帰分析には、公知の方法を採用することができる。この機能は、回帰式生成モジュール12によって実現される。   The arithmetic control unit 3 performs logistic regression analysis of the patient explanatory variables A, α, and β related to the learning data and NASH / non-NASH data, and the patient is NASH based on the explanatory variables A, α, and β. A regression equation for obtaining a probability value is generated. The arithmetic control unit 3 stores the generated regression equation in the storage unit 4. A well-known method can be employ | adopted for logistic regression analysis. This function is realized by the regression equation generation module 12.

演算制御部3は、学習データに係る患者のSPECTデータから求めた説明変数A、α、βを回帰式に代入し、それぞれの患者がNASHである確率値を計算する。この機能は、NASH確率値算出モジュール13によって実現される。   The arithmetic control unit 3 substitutes the explanatory variables A, α, and β obtained from the SPECT data of the patient related to the learning data into the regression equation, and calculates the probability value that each patient is NASH. This function is realized by the NASH probability value calculation module 13.

そして、各患者がNASHである確率値と、各患者がNASHか非NASHかのデータに基づいて、ROC曲線を生成し、NASHと非NASHとを鑑別する適切なカットオフ値を決定する。適切なカットオフ値は、例えば、感度が所定値以上で、感度/(1−特異度)が最大となる値というように基準を決めておけば、ROC曲線に基づいてカットオフ値を決定することができる。この機能は、カットオフ値決定モジュール14によって実現される。   Then, based on the probability value that each patient is NASH and the data that each patient is NASH or non-NASH, an ROC curve is generated, and an appropriate cutoff value for distinguishing NASH from non-NASH is determined. The appropriate cutoff value is determined based on the ROC curve, for example, if the standard is determined such that the sensitivity is a predetermined value or more and the sensitivity / (1-specificity) is the maximum. be able to. This function is realized by the cut-off value determination module 14.

ROC曲線とは、縦軸に感度、横軸に(1−特異度)をとったグラフである。感度とは、真陽性率、すなわち、NASH患者をNASH患者として判定する確率である。(1−特異度)は、偽陽性率、すなわち、非NASH患者をNASH患者と判定する確率である。NASHの確率値の高い値をカットオフ値として設定すれば真陽性率も低いが偽陽性率も低く、カットオフ値を低くするに従って真陽性率及び偽陽性率が高くなっていく。演算制御部3は、真陽性率と偽陽性率とに基づいて、NASHと非NASHを鑑別するカットオフ値を決定する。   The ROC curve is a graph with sensitivity on the vertical axis and (1-specificity) on the horizontal axis. Sensitivity is the true positive rate, that is, the probability of determining a NASH patient as a NASH patient. (1-specificity) is the false positive rate, that is, the probability of determining a non-NASH patient as a NASH patient. If a high NASH probability value is set as the cutoff value, the true positive rate is low but the false positive rate is low, and the true positive rate and false positive rate increase as the cutoff value is lowered. The arithmetic control unit 3 determines a cutoff value for discriminating between NASH and non-NASH based on the true positive rate and the false positive rate.

[NASHを鑑別する機能]
演算制御部3の処理に先立って、データ取得部2は、検査対象の患者の123I−BMIPPを用いたSPECTデータを取得する。取得したデータは、記憶部4に記憶してもよい。
[Function for distinguishing NASH]
Prior to the processing of the arithmetic control unit 3, the data acquisition unit 2 acquires SPECT data using 123 I-BMIPP of the patient to be examined. The acquired data may be stored in the storage unit 4.

演算制御部3は、データ取得部2にて取得した患者の123I−BMIPPを用いたSPECTデータのフィッティングカーブを生成する。すなわち、f=Ae-βt−Ae-αtの説明変数A、α、βを求める。この機能は、フィッティングカーブ生成モジュール11によって実現される。 The arithmetic control unit 3 generates a SPECT data fitting curve using the 123 I-BMIPP of the patient acquired by the data acquisition unit 2. That is, the explanatory variables A, α and β of f = Ae −βt −Ae −αt are obtained. This function is realized by the fitting curve generation module 11.

次に、演算制御部3は、記憶部4から回帰式を読み出し、読み出した回帰式に対して、検査対象の患者について求めた説明変数A、α、βの値を代入し、患者が、NASHである確率値を計算する。この機能は、NASH確率値算出モジュール13によって実現される。続いて、演算制御部3は、NASHの確率値とカットオフ値とを比較して、患者がNASHであるか否かを鑑別する。この機能はNASH鑑別モジュール15によって実現される。   Next, the arithmetic control unit 3 reads the regression equation from the storage unit 4 and substitutes the values of the explanatory variables A, α, and β obtained for the patient to be examined into the read regression equation. Compute a probability value that is This function is realized by the NASH probability value calculation module 13. Subsequently, the arithmetic control unit 3 compares the NASH probability value with the cut-off value, and discriminates whether or not the patient is NASH. This function is realized by the NASH discrimination module 15.

演算制御部3は、NASHの鑑別結果を出力部5に送り、出力部5がNASHの鑑別結果を外部に出力する。出力の方法としては、例えば、モニタに表示することができる。この機能は、出力モジュール16によって実現される。   The arithmetic control unit 3 sends the NASH discrimination result to the output unit 5, and the output unit 5 outputs the NASH discrimination result to the outside. As an output method, for example, it can be displayed on a monitor. This function is realized by the output module 16.

図3(a)及び図3(b)は、本実施の形態のNASH診断装置1の動作について説明する。図3(a)は、NASH診断装置1が、回帰式を生成する動作を示すフローチャートであり、図3(b)は、患者がNASHであるか否かを鑑別するフローチャートである。   3A and 3B explain the operation of the NASH diagnostic apparatus 1 according to the present embodiment. FIG. 3A is a flowchart showing the operation of the NASH diagnostic apparatus 1 for generating a regression equation, and FIG. 3B is a flowchart for discriminating whether or not the patient has NASH.

NASH診断装置1は、まず、SPECT装置6から学習データを取得する(S10)。学習データは、上述したとおり、NASHであることが既知の少なくとも一の被検者、及び、非NASHであることが既知の少なくとも一の被検者の123I−BMIPPによるSPECTデータである。NASH診断装置1は、取得したデータを記憶部4に記憶してもよい。 The NASH diagnostic device 1 first acquires learning data from the SPECT device 6 (S10). As described above, the learning data is SPECT data by 123 I-BMIPP of at least one subject known to be NASH and at least one subject known to be non-NASH. The NASH diagnostic device 1 may store the acquired data in the storage unit 4.

NASH診断装置1は、取得した学習データに係るSPECTデータのフィッティングカーブを生成し(S11)、求めた説明変数A、α、βとNASH/非NASHのデータに基づいてロジスティック回帰分析を行い、回帰式を生成する。NASH診断装置1は、生成した回帰式を記憶部4に記憶する(S12)。   The NASH diagnostic apparatus 1 generates a SPECT data fitting curve related to the acquired learning data (S11), performs logistic regression analysis based on the obtained explanatory variables A, α, β and NASH / non-NASH data, and performs regression. Generate an expression. The NASH diagnostic device 1 stores the generated regression equation in the storage unit 4 (S12).

次に、NASH診断装置1は、学習データに係る患者のSPECTデータから求めた説明変数を回帰式に代入し、NASHの確率値を計算する。そして、NASH診断装置1は、上述したとおり、カットオフ値を変数とするROC曲線に基づいて、NASH/非NASHを鑑別するためのカットオフ値を決定し、記憶部4に記憶する(S13)。なお、本実施の形態では、回帰式を生成するために使った学習データを用いているが、NASH/非NASHが既知の患者のデータであれば、回帰式を生成するために使った患者のデータとは異なるデータを用いることも可能である。   Next, the NASH diagnostic device 1 substitutes the explanatory variable obtained from the SPECT data of the patient related to the learning data into the regression equation, and calculates the probability value of NASH. Then, as described above, the NASH diagnostic device 1 determines a cutoff value for distinguishing NASH / non-NASH based on the ROC curve having the cutoff value as a variable, and stores it in the storage unit 4 (S13). . In this embodiment, the learning data used to generate the regression equation is used. However, if the NASH / non-NASH data is known, the data of the patient used to generate the regression equation is used. It is also possible to use data different from the data.

次に、図3(b)を参照して、NASH診断装置1によって、患者がNASHか否かを鑑別する動作について説明する。NASH装置は、まず、SPECT装置6から患者の123I−BMIPPによるSPECTデータを取得する(S20)。NASH診断装置1は、取得したSPECTデータを記憶部4に記憶してもよい。 Next, with reference to FIG. 3B, an operation for identifying whether or not the patient is NASH by the NASH diagnostic apparatus 1 will be described. First, the NASH device acquires SPECT data based on 123 I-BMIPP of the patient from the SPECT device 6 (S20). The NASH diagnostic device 1 may store the acquired SPECT data in the storage unit 4.

NASH診断装置1は、取得したSPECTデータのフィッティングカーブを生成し(S21)、説明変数A、α、βを求める。続いて、NASH診断装置1は、記憶部4から回帰式を読み出し、説明変数A、α、βを回帰式に代入して、患者がNASHである確率値を算出する(S22)。   The NASH diagnostic apparatus 1 generates a fitting curve of the acquired SPECT data (S21), and obtains explanatory variables A, α, and β. Subsequently, the NASH diagnostic apparatus 1 reads the regression equation from the storage unit 4 and substitutes the explanatory variables A, α, and β into the regression equation to calculate a probability value that the patient is NASH (S22).

次に、NASH診断装置1は、NASHの確率値がカットオフ値以上であるか否かによって、患者がNASHであるか、非NASHであるかを鑑別する(S23)。NASH診断装置1は、鑑別した結果を、例えばモニタに表示する等して、外部に出力する(S24)。以上、実施の形態のNASH診断装置1の構成および動作について説明した。   Next, the NASH diagnostic device 1 discriminates whether the patient is NASH or non-NASH depending on whether or not the NASH probability value is equal to or greater than the cutoff value (S23). The NASH diagnostic apparatus 1 outputs the identified result to the outside, for example, by displaying it on a monitor (S24). The configuration and operation of the NASH diagnostic apparatus 1 according to the embodiment have been described above.

図4(a)〜図4(c)は、肝シンチグラフィの結果を示す図である。図4(a)は単純性脂肪肝(SS)の時間放射能曲線、図4(b)は早期(Brunt病理分類Stage1−2)NASH(eNASH)の時間放射能曲線、図4(c)進行(Stage3−4)NASH(aNASH)の時間放射能曲線をそれぞれ示すグラフである。図4(a)〜図4(c)に見られるように、SSとeNASHとaNASHでは、肝への123I−BMIPPの取り込みの高さや、その速度等が異なる。 4 (a) to 4 (c) are diagrams showing the results of liver scintigraphy. FIG. 4 (a) shows the time activity curve of simple fatty liver (SS), FIG. 4 (b) shows the time activity curve of early (Brunt pathological classification Stage1-2) NASH (eNASH), and FIG. 4 (c) progression. (Stage3-4) It is a graph which respectively shows the time activity curve of NASH (aNASH). 4A to 4C, SS, eNASH, and aNASH have different levels of 123 I-BMIPP uptake into the liver, their speed, and the like.

本実施の形態のNASH診断装置1は、学習データに係る患者の123I−BMIPPによるSPECTデータを、f=Ae-βt−Ae-αtをモデル式としてフィッティングカーブを求め、求めたフィッティングカーブを用いて生成した回帰式に基づいてNASHと非NASHとを鑑別する構成により、非侵襲的に適切なNASHの診断を行うことができる。これは、f=Ae-βt−Ae-αtのうち、説明変数Aが肝への123I−BMIPPの取り込みの高さを表し、説明変数αが123I−BMIPPの流出速度を表し、説明変数βが123I−BMIPPの流入速度を表し、NASH、非NASHのカーブの特性の違いが説明変数に適切に現れるためと考えられる。 The NASH diagnostic apparatus 1 according to the present embodiment obtains a fitting curve of SPECT data by 123 I-BMIPP of a patient related to learning data using f = Ae −βt −Ae −αt as a model formula, and uses the obtained fitting curve. With the configuration for distinguishing NASH from non-NASH based on the regression equation generated in this way, it is possible to diagnose NASH appropriately noninvasively. This means that among f = Ae −βt −Ae −αt , the explanatory variable A represents the height of 123 I-BMIPP uptake into the liver, the explanatory variable α represents the outflow rate of 123 I-BMIPP, and the explanatory variable β represents the inflow speed of 123 I-BMIPP, and it is considered that the difference in the characteristics of the NASH and non-NASH curves appears appropriately in the explanatory variables.

また、本実施の形態のNASH診断装置1は、ROC曲線に基づいて、NASHと非NASHとを鑑別するカットオフ値を決定し、このカットオフ値を使って、NASH/非NASHの鑑別を行うので、容易に、NASHの診断を行うことができる。   Further, the NASH diagnostic apparatus 1 according to the present embodiment determines a cutoff value for distinguishing NASH from non-NASH based on the ROC curve, and performs NASH / non-NASD discrimination using this cutoff value. Therefore, NASH can be diagnosed easily.

なお、上述した実施の形態においては、NASH診断装置1、NASH診断方法及びプログラムは、学習データに係る複数の患者の123I−BMIPPによるSPECTデータを取得し、NASHの確率を求める回帰式を生成する段階を含んでいるが、本発明のNASH診断装置1及びNASH診断方法は、あらかじめ学習データを使って求めた回帰式を記憶部4に記憶しておき、当該回帰式を使ってNASH診断を行う装置及び方法も含んでいる。すなわち、図3(b)の動作のみを行う装置、方法およびプログラムも、本発明の範囲に含まれる。 In the above-described embodiment, the NASH diagnostic apparatus 1, the NASH diagnostic method, and the program acquire SPECT data by 123 I-BMIPP of a plurality of patients related to learning data, and generate a regression equation for determining the probability of NASH. However, the NASH diagnosis apparatus 1 and NASH diagnosis method of the present invention store the regression equation obtained beforehand using the learning data in the storage unit 4 and perform NASH diagnosis using the regression equation. Also included is an apparatus and method for performing. That is, an apparatus, a method, and a program that perform only the operation of FIG. 3B are also included in the scope of the present invention.

また、上述した実施の形態においては、NASH診断装置1、NASH診断方法及びプログラムは、NASHの確率値を求めた後、カットオフ値に基づいて、NASHか非NASHかの鑑別を行い、その結果を出力する例について説明したが、本発明のNASH診断装置1、NASH診断方法及びプログラムは、NASHの確率値を出力することとしてもよい。NASHの確率値を提供することにより、NASHの確率値の結果とその他の所見とに基づいて、NASHか否かを医師が総合的に判断することができる。   Further, in the above-described embodiment, the NASH diagnostic device 1, the NASH diagnostic method, and the program, after obtaining the NASH probability value, discriminate between NASH and non-NASH based on the cut-off value. However, the NASH diagnosis apparatus 1, NASH diagnosis method, and program of the present invention may output NASH probability values. By providing the NASH probability value, the doctor can comprehensively determine whether or not the NASH is based on the result of the NASSH probability value and other findings.

次に、本発明の実施例について説明する。
[学習データ]
肝生検にて病理学診断に至ったNAFLD患者14名(男/女=11/3(人)、平均年齢49(歳))、および健常人1名を対象とした。NAFLD患者14名の内訳は、単純性脂肪肝(SS)4名、早期NASH(eNASH)6名、進行NASH(aNASH)4名であった。
Next, examples of the present invention will be described.
[Learning data]
The subjects were 14 NAFLD patients (male / female = 11/3 (person), average age 49 (years)) and one healthy person who had pathological diagnosis by liver biopsy. The breakdown of 14 NAFLD patients was 4 simple fatty liver (SS), 6 early NASH (eNASH), and 4 advanced NASH (aNASH).

[方法]
患者に123I−BMIPP111MBqを急速静注後、2ヘッドガンマカメラを用いて、60分後までの肝シンチグラフィを行った。
[Method]
After rapid intravenous injection of 123 I-BMIPP111MBq to the patient, liver scintigraphy up to 60 minutes later was performed using a 2-head gamma camera.

[結果]
肝シンチグラフィによって得られたSPECTデータのフィッティングカーブを求めた。フィッティングカーブの説明変数A、α、βは、以下のような結果となった。
[result]
A fitting curve of SPECT data obtained by liver scintigraphy was obtained. The explanatory variables A, α, and β of the fitting curve gave the following results.

続いて、上記結果をロジスティック回帰分析を行うと、オッズ比Xbは、
Xb=−14.11−0.58×A+16383×β+2.321×α
と求められた。オッズ比Xbを用いて、NASHの確率値Prを求める式は、
Pr=1÷(e-Xb+1)
で表される。
Subsequently, when the above result is subjected to logistic regression analysis, the odds ratio Xb is
Xb = −14.11−0.58 × A / 16363 × β + 2.321 × α
I was asked. Using the odds ratio Xb, the formula for determining the NAS probability value Pr is:
Pr = 1 / (e −Xb +1)
It is represented by

この回帰式を用いて、15名の患者のNASHの確率値を求めたところ、以下のような結果となった。
Using this regression equation, the probability values of NASH of 15 patients were obtained, and the following results were obtained.

上記の結果に基づいて記載したROC曲線は以下のとおりである。
The ROC curve described based on the above results is as follows.

ROC曲線において、AUC(Area Under Curve)の面積が大きいほど、良い判定方法である。本実施例においては、AUCは、0.980となり、非常に良い鑑別方法であることが示された。   In the ROC curve, the larger the area of AUC (Area Under Curve), the better the determination method. In this example, AUC was 0.980, indicating that this is a very good discrimination method.

本発明は、123I−BMIPPを用いた肝シンチグラフィにより、検査対象の患者がNASHであるか否かを非侵襲的に鑑別することができ、NASHか否かの診断に有用である。 INDUSTRIAL APPLICABILITY According to the present invention, liver scintigraphy using 123 I-BMIPP can noninvasively identify whether or not a patient to be examined is NASH, and is useful for diagnosing whether or not it is NASH.

1 NASH診断装置
2 データ取得部
3 演算制御部
4 記憶部
5 出力部
6 SPECT装置
10 NASH診断プログラム
11 フィッティングカーブ生成モジュール
12 回帰式生成モジュール
13 NASH確率値算出モジュール
14 カットオフ値決定モジュール
15 NASH鑑別モジュール
16 出力モジュール
DESCRIPTION OF SYMBOLS 1 NASH diagnostic apparatus 2 Data acquisition part 3 Operation control part 4 Memory | storage part 5 Output part 6 SPECT apparatus 10 NASH diagnostic program 11 Fitting curve generation module 12 Regression expression generation module 13 NASH probability value calculation module 14 Cutoff value determination module 15 NASH discrimination Module 16 Output module

Claims (10)

123I−BMIPPを用いたシンチグラフィによって得られた時間放射能曲線のデータに基づいて、検査対象の患者がNASHであるか否かを診断するためのプログラムであって、コンピュータに、
前記検査対象の患者の時間放射能曲線のデータについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、
NASHであることが既知の少なくとも一の被検者及び非NASHであることが既知の少なくとも一の被検者に対して、それぞれ123I−BMIPPを用いたシンチグラフィを行って得られた時間放射能曲線の各データについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求め、前記被検者のフィッティングカーブをロジスティック回帰分析することによって生成された、前記説明変数に基づいて患者がNASHである確率値を求める回帰式を記憶した記憶部から、前記回帰式を読み出すステップと、
前記フィッティングカーブの説明変数を前記回帰式に代入して、前記検査対象の患者が NASHである確率値を算出するステップと、
前記確率値を出力するステップと、
を実行させるプログラム。
A program for diagnosing whether or not a patient to be examined is NASH based on data of a time activity curve obtained by scintigraphy using 123 I-BMIPP, comprising:
Fitting of 123 I-BMIPP to the liver, 123 I-BMIPP inflow rate, and 123 I-BMIPP outflow rate as explanatory variables for the time activity curve data of the patient to be examined A step of determining a curve;
Time radiation obtained by performing scintigraphy using 123 I-BMIPP on at least one subject known to be NASH and at least one subject known to be non-NASH, respectively. for each data capacity curve, determined the height of the 123 I-BMIPP incorporation into liver, the inflow rate of 123 I-BMIPP, the fitting curve of an explanatory variable and a flux rate of 123 I-BMIPP, said subject Reading the regression equation from a storage unit that stores a regression equation that is generated by performing logistic regression analysis of the person's fitting curve and obtaining a probability value that the patient is NASH based on the explanatory variable;
Substituting the explanatory variable of the fitting curve into the regression equation to calculate a probability value that the patient to be examined is NASH;
Outputting the probability value;
A program that executes
前記フィッティングカーブfは、123I−BMIPPの取り込みの高さをA、123I−BMIPPの流出速度をα、123I−BMIPPの流入速度をβとして、
f=Ae-βt−Ae-αt
である請求項1に記載のプログラム。
Said fitting curve f is the flow rate of 123 I-BMIPP of the height of the uptake A, 123 I-BMIPP α, as the inlet velocity of 123 I-BMIPP beta,
f = Ae -βt -Ae -αt
The program according to claim 1.
前記コンピュータに、
前記回帰式を用いて前記被検者がNASHである確率値を求め、NASHと非NASHとを鑑別するための閾値を変数として生成されるROC曲線に基づいて決定されたNASHと非NASHとを鑑別するための閾値を記憶した記憶部から、前記閾値を読み出すステップと、
前記検査対象の患者がNASHである確率値と前記閾値とに基づいて、検査対象の患者がNASHか非NASHかを鑑別するステップと、
を実行させる請求項1又は2に記載のプログラム。
In the computer,
A probability value that the subject is NASH is obtained using the regression equation, and NASH and non-NASH determined based on a ROC curve generated using a threshold value for distinguishing NASH from non-NASH as a variable. Reading the threshold value from the storage unit storing the threshold value for discrimination;
Differentiating whether the patient to be examined is NASH or non-NASH based on the probability value that the patient to be examined is NASH and the threshold value;
The program according to claim 1 or 2, wherein the program is executed.
123I−BMIPPを用いたシンチグラフィによって得られた時間放射能曲線のデータに基づいて、検査対象の患者がNASHであるか否かを診断するためのプログラムであって、コンピュータに、
NASHであることが既知の少なくとも一の被検者及び非NASHであることが既知の少なくとも一の被検者に対して、それぞれ123I−BMIPPを用いたシンチグラフィを行って得られた時間放射能曲線の各データについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、
前記被検者のフィッティングカーブをロジスティック回帰分析することにより、前記説明変数に基づいて患者がNASHである確率値を求める回帰式を生成するステップと、
前記検査対象の患者の時間放射能曲線のデータについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、
前記フィッティングカーブの説明変数を前記回帰式に代入して、前記検査対象の患者が NASHである確率値を算出するステップと、
前記確率値を出力するステップと、
を実行させるプログラム。
A program for diagnosing whether or not a patient to be examined is NASH based on data of a time activity curve obtained by scintigraphy using 123 I-BMIPP, comprising:
Time radiation obtained by performing scintigraphy using 123 I-BMIPP on at least one subject known to be NASH and at least one subject known to be non-NASH, respectively. for each data capacity curve, and determining the height of the 123 I-BMIPP incorporation into liver, the inflow rate of 123 I-BMIPP, the fitting curve of an explanatory variable and a flux rate of 123 I-BMIPP,
Generating a regression equation for determining a probability value that the patient is NASH based on the explanatory variable by logistic regression analysis of the fitting curve of the subject;
Fitting of 123 I-BMIPP to the liver, 123 I-BMIPP inflow rate, and 123 I-BMIPP outflow rate as explanatory variables for the time activity curve data of the patient to be examined A step of determining a curve;
Substituting the explanatory variable of the fitting curve into the regression equation to calculate a probability value that the patient to be examined is NASH;
Outputting the probability value;
A program that executes
前記フィッティングカーブfは、123I−BMIPPの取り込みの高さをA、123I−BMIPPの流出速度をα、123I−BMIPPの流入速度をβとして、
f=Ae-βt−Ae-αt
である請求項4に記載のプログラム。
Said fitting curve f is the flow rate of 123 I-BMIPP of the height of the uptake A, 123 I-BMIPP α, as the inlet velocity of 123 I-BMIPP beta,
f = Ae -βt -Ae -αt
The program according to claim 4.
前記コンピュータに、
前記回帰式を用いて前記被検者がNASHである確率値を求め、NASHと非NASHとを鑑別するための閾値を変数として生成されるROC曲線に基づいてNASHと非NASHとを鑑別するための閾値を決定するステップと、
前記検査対象の患者がNASHである確率値と前記閾値とに基づいて、検査対象の患者 がNASHか非NASHかを鑑別するステップと、
を実行させる請求項4又は5に記載のプログラム。
In the computer,
To obtain a probability value that the subject is NASH using the regression equation, and to distinguish NASH from non-NASH based on a ROC curve generated using a threshold value for distinguishing NASH from non-NASH as a variable Determining a threshold of
Differentiating whether the patient to be examined is NASH or non-NASH based on the probability value that the patient to be examined is NASH and the threshold; and
The program according to claim 4 or 5, wherein the program is executed.
NASH診断装置が、NASHであることが既知の少なくとも一の被検者及び非NASHであることが既知の少なくとも一の被検者に対して、123I−BMIPPを用いたシンチグラフィを行って得られた、時間放射能曲線のデータをそれぞれ取得するステップと、
前記NASH診断装置が、前記被検者の前記時間放射能曲線のデータについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、
前記NASH診断装置が、前記被検者のフィッティングカーブをロジスティック回帰分析することにより、前記説明変数に基づいて患者がNASHである確率値を求める回帰式を生成するステップと、
前記NASH診断装置が、検査対象の患者に対して、123I−BMIPPを用いたシンチグラフィを行って時間放射能曲線のデータを取得するステップと、
前記NASH診断装置が、前記検査対象の患者の時間放射能曲線のデータについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求めるステップと、
前記NASH診断装置が、前記フィッティングカーブの説明変数を前記回帰式に代入して、前記検査対象の患者がNASHである確率値を算出するステップと、
前記NASH診断装置が、前記確率値を出力するステップと、
を備えるNASH診断装置の作動方法。
The NASH diagnostic apparatus is obtained by performing scintigraphy using 123 I-BMIPP for at least one subject known to be NASH and at least one subject known to be non-NASH. It was, acquiring each data time activity curve,
The NASH diagnosis apparatus, the data of the time-activity curve of the subject, and the height of the 123 I-BMIPP incorporation into liver, the inflow rate of 123 I-BMIPP, the outflow speed of 123 I-BMIPP Obtaining a fitting curve with as an explanatory variable,
The NASH diagnostic device performs a logistic regression analysis of the subject's fitting curve to generate a regression equation for obtaining a probability value that the patient is NASH based on the explanatory variable;
The NASH diagnostic device performs scintigraphy using 123 I-BMIPP for a patient to be examined to obtain time-activity curve data;
The NASH diagnosis apparatus, the data of the time-activity curve of patients of the test object, the height of the 123 I-BMIPP incorporation into liver, the inflow rate of 123 I-BMIPP, the outflow speed of 123 I-BMIPP Obtaining a fitting curve with as an explanatory variable,
The NASH diagnostic device substituting explanatory variables of the fitting curve into the regression equation to calculate a probability value that the patient to be examined is NASH;
The NASH diagnostic device outputting the probability value;
A method of operating a NASH diagnostic device comprising :
前記フィッティングカーブfは、123I−BMIPPの取り込みの高さをA、123I−BMIPPの流出速度をα、123I−BMIPPの流入速度をβとして、
f=Ae-βt−Ae-αt
である請求項7に記載のNASH診断装置の作動方法。
Said fitting curve f is the flow rate of 123 I-BMIPP of the height of the uptake A, 123 I-BMIPP α, as the inlet velocity of 123 I-BMIPP beta,
f = Ae -βt -Ae -αt
The operation method of the NASH diagnostic apparatus according to claim 7.
前記NASH診断装置が、前記回帰式を用いて前記被検者がNASHである確率値を求め、NASHと非NASHとを鑑別するための閾値を変数として生成されるROC曲線に基づいてNASHと非NASHとを鑑別するための閾値を決定するステップと、
前記NASH診断装置が、前記検査対象の患者がNASHである確率値と前記閾値とに基づいて、検査対象の患者がNASHか非NASHかを鑑別するステップと、
を有する請求項7又は8に記載のNASH診断装置の作動方法。
The NASH diagnosis device obtains a probability value that the subject is NASH using the regression equation, and determines whether the subject is NASH or not based on a ROC curve generated using a threshold for distinguishing NASH from non-NASH as a variable. Determining a threshold for differentiating from NASH;
The NASH diagnostic device discriminates whether the patient to be inspected is NASH or non-NASH based on the probability value that the patient to be inspected is NASH and the threshold; and
The operation | movement method of the NASH diagnostic apparatus of Claim 7 or 8 which has these.
NASHであることが既知の少なくとも一の被検者及び非NASHであることが既知の少なくとも一の被検者に対して、123I−BMIPPを用いたシンチグラフィを行って得られた時間放射能曲線の各データについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求め、前記被検者のフィッティングカーブをロジスティック回帰分析することによって生成された、前記説明変数に基づいて患者がNASHである確率値を求める回帰式を記憶した記憶部と、
検査対象の患者に対して、123I−BMIPPを用いたシンチグラフィを行って得られた時間放射能曲線のデータを取得するデータ取得部と、
前記検査対象の患者の時間放射能曲線のデータについて、肝への123I−BMIPPの取り込みの高さと、123I−BMIPPの流入速度と、123I−BMIPPの流出速度とを説明変数とするフィッティングカーブを求め、前記フィッティングカーブの説明変数を前記回帰式に代入して、前記検査対象の患者がNASHである確率値を算出する演算制御部と、
前記確率値を出力する出力部と、
を備えるNASH診断装置。
Time radioactivity obtained by scintigraphy using 123 I-BMIPP for at least one subject known to be NASH and at least one subject known to be non-NASH for each data curve, the height of the 123 I-BMIPP incorporation into liver, determined the inflow rate of 123 I-BMIPP, the fitting curve of an explanatory variable and a flux rate of 123 I-BMIPP, the subject A storage unit that stores a regression equation that is generated by performing logistic regression analysis of the fitting curve of and calculates a probability value that the patient is NASH based on the explanatory variable;
A data acquisition unit for acquiring data of a time activity curve obtained by performing scintigraphy using 123 I-BMIPP for a patient to be examined;
Fitting of 123 I-BMIPP to the liver, 123 I-BMIPP inflow rate, and 123 I-BMIPP outflow rate as explanatory variables for the time activity curve data of the patient to be examined Calculating a curve, substituting explanatory variables of the fitting curve into the regression equation, and calculating a probability value that the patient to be examined is NASH;
An output unit for outputting the probability value;
A NASH diagnostic apparatus comprising:
JP2014022570A 2014-02-07 2014-02-07 NASH diagnostic program, NASH diagnostic device operating method, and NASH diagnostic device Active JP6272065B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014022570A JP6272065B2 (en) 2014-02-07 2014-02-07 NASH diagnostic program, NASH diagnostic device operating method, and NASH diagnostic device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014022570A JP6272065B2 (en) 2014-02-07 2014-02-07 NASH diagnostic program, NASH diagnostic device operating method, and NASH diagnostic device

Publications (2)

Publication Number Publication Date
JP2015148555A JP2015148555A (en) 2015-08-20
JP6272065B2 true JP6272065B2 (en) 2018-01-31

Family

ID=53892005

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014022570A Active JP6272065B2 (en) 2014-02-07 2014-02-07 NASH diagnostic program, NASH diagnostic device operating method, and NASH diagnostic device

Country Status (1)

Country Link
JP (1) JP6272065B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109313246A (en) * 2016-06-13 2019-02-05 国立大学法人九州大学 The acquisition methods of free radical depletion rate information and the determination method of NASH
CA3134278C (en) * 2018-03-07 2023-08-29 Siemens Medical Solutions Usa, Inc Calibration bias reduction in a pressurized gas ion chanber-based dose calibrator

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008538893A (en) * 2005-04-28 2008-11-13 国立大学法人高知大学 Method for detecting lipid metabolism insufficiency and test agent used therefor

Also Published As

Publication number Publication date
JP2015148555A (en) 2015-08-20

Similar Documents

Publication Publication Date Title
JP7090546B2 (en) Perfusion Digital Subtraction Angiography
JP6873981B2 (en) Mobile FFR simulation
JP6272065B2 (en) NASH diagnostic program, NASH diagnostic device operating method, and NASH diagnostic device
JP2011043489A (en) Image processing apparatus, program and method for evaluating local wall motion state of cardiac muscle
JP2008132320A5 (en)
JP3502513B2 (en) Ultrasonic image processing method and ultrasonic image processing apparatus
JP4964450B2 (en) Image processing method, image processing program, and image processing apparatus
US20150359500A1 (en) Synchronous physiological measurements for cardiac acquisitions
Buisson et al. What is the role of diffusion-weighted imaging in ileocolonic Crohn's disease?
Min et al. New frontiers in CT angiography: physiologic assessment of coronary artery disease by multidetector CT
JPWO2018159708A1 (en) Apparatus and method for blood flow analysis and program
Marcus et al. Myocardial perfusion imaging by computed tomography: today and tomorrow
JP5422543B2 (en) Blood flow estimation apparatus, method, and computer program
JP5710166B2 (en) Image analysis apparatus and method
RU178086U1 (en) Automated device for the diagnosis of heart disease
JP6501432B2 (en) Hemodynamic imaging diagnostic apparatus and diagnostic method
JP6568375B2 (en) Ophthalmic information processing system, image processing apparatus, and image processing method
JP2018077164A (en) Brain image analysis method, brain image analyzer, and program
JP5713868B2 (en) Cerebral blood flow quantification device, cerebral blood flow quantification program, and cerebral blood flow quantification method
Kuzan et al. A new parameter in the diagnosis of acute appendicitis: Ileocolic artery and vein diameter measurements
JP6898969B2 (en) Ophthalmic information processing system and ophthalmic information processing method
JP5489602B2 (en) Cerebral blood flow calculation program, recording medium, and cerebral blood flow calculation method
JP5700712B2 (en) Patrac analysis of 3D dynamic myocardial nuclear medicine image data
JP5449874B2 (en) Medical image processing apparatus and X-ray CT apparatus
JP5695230B1 (en) Cerebral blood flow quantitative analysis device, cerebral blood flow quantitative analysis system and program

Legal Events

Date Code Title Description
RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20160726

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160928

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170828

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170905

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20171020

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20171212

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171228

R150 Certificate of patent or registration of utility model

Ref document number: 6272065

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250