JP2017062138A - 中性脂肪濃度検量装置、電子機器、及び、検量線作成方法 - Google Patents

Abstract

【課題】中性脂肪に関する検量を精度良く行う技術を提供する。【解決手段】被検体についての中性脂肪濃度を求める中性脂肪濃度検量装置は、前記被検体についての観測分光データを取得する被検体観測データ取得部と、中性脂肪に対応する独立成分又は主成分と、検量用の単回帰式と、含む検量用データを取得する検量用データ取得部と、前記被検体についての観測分光データと前記検量用データとに基づいて、前記被検体についての中性脂肪に対する混合係数を求める混合係数算出部と、前記中性脂肪に対応する混合係数と中性脂肪濃度との関係を示す前記単回帰式と、前記混合係数算出部によって求められた混合係数とに基づいて、前記中性脂肪濃度を算出する目的成分量算出部と、を含み、前記混合係数算出部は、前記中性脂肪に対する独立成分又は主成分として、予め選択された波長にピークを有する成分を使用する。【選択図】図２

Description

本発明は、中性脂肪を含む被検体の観測データから被検体の中性脂肪濃度を求める技術に関する。

特許文献１には、近赤外光を用いて生体の吸光スペクトルを測定し、この吸光スペクトに応じて中性脂肪濃度の定量を行う技術が記載されている。この従来技術では、１７２８±１０ｎｍより選択した１波長と、１６６１±１０ｎｍあるいは１６６８±１０ｎｍあるいは１６７９±１０ｎｍあるいは１６８５±１０ｎｍあるいは１７２１±１０ｎｍあるいは１７２８±１０ｎｍあるいは１７４２±１０ｎｍあるいは１７５４±１０ｎｍより選択した１波長の少なくとも２波長あるいは２波長帯を選択して、第２の波長域あるいは２波長帯を選択して、第２の波長域あるいは第３の波長域と組み合わせて、中性脂肪濃度を定量している。

特開平１１−６４２１８号公報

しかし、上記従来技術では、中性脂肪の情報を多く含む波長範囲を選択しているものの、他の生体成分の情報が重なりあって含まれていることには変わりはなく、中性脂肪の情報だけを分離抽出することは原理的に難しい。そのため、十分な検量精度が得られない場合があるという問題があった。

本発明は、上述の課題の少なくとも一部を解決するためになされたものであり、以下の形態または適用例として実現することが可能である。

（１）本発明の第１の形態によれば、目的成分としての中性脂肪を含む被検体についての中性脂肪濃度を求める中性脂肪濃度検量装置が提供される。この中性脂肪濃度検量装置は、前記被検体についての観測分光データを取得する被検体観測データ取得部と、中性脂肪に対応する独立成分又は主成分と、検量用の単回帰式と、含む検量用データを取得する検量用データ取得部と、前記被検体についての観測分光データと前記検量用データとに基づいて、前記被検体についての中性脂肪に対する混合係数を求める混合係数算出部と、前記中性脂肪に対応する混合係数と中性脂肪濃度との関係を示す前記単回帰式と、前記混合係数算出部によって求められた混合係数とに基づいて、前記中性脂肪濃度を算出する目的成分量算出部と、を含み、前記混合係数算出部は、前記中性脂肪に対する独立成分又は主成分として、予め選択された波長にピークを有する成分を使用する。
この中性脂肪濃度検量装置によれば、中性脂肪に対応する波長として予め選択された波長にピークを有する独立成分又は主成分を用いて中性脂肪濃度の検量が行われるので、中性脂肪濃度に関する検量精度を高めることが可能である。

（２）上記中性脂肪濃度検量装置において、前記予め選択された波長は、９３０±３０ｎｍ又は９４５±３０ｎｍの一方と、１０６０±３０ｎｍとの少なくとも１つの波長としてもよい。
この構成によれば、中性脂肪に対応する波長として、９３０±３０ｎｍ又は９４５±３０ｎｍの一方と、１０６０±３０ｎｍとの少なくとも１つの波長が選択されるので、中性脂肪濃度に関する検量精度をより高めることが可能である。

（３）上記中性脂肪濃度検量装置において、前記予め選択された波長は、１１５０±３０ｎｍ又は１１８２±３０ｎｍの波長としてもよい。
この構成によれば、中性脂肪に対応する波長として、１１５０±３０ｎｍ又は１１８２±３０ｎｍの波長が選択されるので、中性脂肪濃度に関する検量精度をより高めることが可能である。

（４）本発明の第２の形態によれば、中性脂肪を含む被検体の観測分光データから、前記被検体の中性脂肪濃度を導くことに用いる検量線を作成する検量線作成方法が提供される。この検量線作成方法は、前記被検体の複数のサンプルについての前記観測分光データを取得するサンプル観測データ取得工程と、前記各サンプルについての中性脂肪濃度を取得する中性脂肪濃度取得工程と、前記サンプル毎の観測分光データを複数の独立成分又は主成分に分離したときの複数の独立成分又は主成分を推定する推定工程と、前記複数のサンプルの中性脂肪濃度と、前記サンプル毎の前記観測分光データにおける前記独立成分又は主成分の混合係数とに基づいて、前記検量線の回帰式を求める回帰式算出工程と、を含み、前記推定工程は、前記独立成分又は主成分として、予め選択された波長にピークを有する成分波形信号を選択する工程である。
この検量線作成方法によれば、中性脂肪に対応する波長として予め選択された波長にピークを有する独立成分又は主成分を用いて中性脂肪濃度の検量が行われるので、中性脂肪濃度に関する検量精度を高めることができる。

本発明は、上述した形態以外の種々の形態でも実現可能であり、例えば、上述した装置を含む電子機器や、上述した装置の各部の機能を実現するコンピュータープログラム、及び、コンピュータープログラムを格納する一時的でない記録媒体（non-transitory storage medium）等の形態等で実現することが可能である。

独立成分分析を利用した検量線作成処理の概要を示す説明図。 目的成分の検量処理の概要を示す説明図。 中性脂肪濃度用の検量線作成処理のフローチャート。 中性脂肪濃度検量用の独立成分の例を示す説明図。 中性脂肪濃度検量用の独立成分の例を示す説明図。 中性脂肪濃度検量用の独立成分の例を示す説明図。 検量線作成処理のフローチャート。 検量線作成処理で使用されるコンピューターを示す説明図。 検量線作成処理で使用される装置の機能ブロック図。 独立成分行列算出部の内部構成の一例を示す機能ブロック図。 測定データセットＤＳ１を模式的に示す説明図。 混合係数推定処理のフローチャート。 推定混合行列∧Ａを説明する説明図。 回帰式の算出処理のフローチャート。 目的成分の検量処理に使用される装置の機能ブロック図。 目的成分の検量処理のフローチャート。 本願の発明者による実験結果。

以下、本発明の実施の形態を以下の順序で説明する。
Ａ．検量線作成処理及び検量処理の概要：
Ｂ．中性脂肪濃度用の検量線作成方法：
Ｃ．検量線作成方法：
Ｄ．目的成分の検量方法：
Ｅ．各種のアルゴリズムとその影響：
Ｆ．変形例：

本実施形態では、以下の略語を使用する。
・ＩＣＡ：独立成分分析(Independent Component Analysis)
・ＳＮＶ：標準正規変量変換(Standard Normal Variate transformation)
・ＰＮＳ：零空間射影法(Project on Null Space)
・ＰＣＡ：主成分分析(Principal Components Analysis)
・ＦＡ：因子分析(Factor Analysis)

Ａ．検量線作成処理及び検量処理の概要：
図１は、独立成分分析を利用した検量線作成処理及の概要を示す説明図である。図１（Ａ）は、複数のサンプルについての観測データ（「測定データ」とも呼ぶ）の一例が示されている。この観測データは分光吸光度であり、例えば、中性脂肪などの複数の化学成分を含むサンプルの分光測定によって得ることができる。検量線作成処理で使用する複数のサンプルとしては、目的成分（例えば中性脂肪）の含有量が既知であるサンプルを使用する。この代わりに、複数のサンプルに含まれる目的成分の含有量を、分析装置で測定してもよい。

検量線の作成の際には、まず、観測データに前処理を行うことによって、観測データに含まれる変動や雑音を低減する（図１（Ｂ））。前処理としては、例えば、観測データの正規化を含む第１前処理と、白色化を含む第２前処理とが行われる。第１前処理では、観測データのさまざまな変動要因（サンプルの状態や測定環境の変化など）による影響を減少するために、零空間射影法を実行することが好ましい。次に、前処理後の観測データに対して独立成分分析処理を実行することによって、複数の独立成分ＩＣ１，ＩＣ２…（図１（Ｃ））が得られる。これらの独立成分ＩＣ１，ＩＣ２…は、各サンプルに含まれる個々の物質成分に対応するデータであり、互いに統計的に独立な成分である。各サンプルの観測データは、これらの独立成分ＩＣ１，ＩＣ２…の線形結合として再現可能である。図１（Ｃ）では、２つの独立成分ＩＣ１，ＩＣ２のみが例示されているが、独立成分の数は２以上の任意の数に適宜設定される。なお、実施形態の説明において、「目的成分」という用語は、サンプルに含まれる物質又は化学成分を意味しており、一方、「独立成分」という用語は、サンプルの観測データと同じデータ長を有するデータを意味する。

次に、図１（Ｄ）〜（Ｆ）に示すように、前処理後の観測データと、独立成分（例えばＩＣ１）との内積を計算する。図１（Ｄ）の観測データは、図１（Ｂ）と同じものである。１つの観測データと、１つの独立成分ＩＣ１との内積を取ると、その観測データに関して１つの内積値が得られる。従って、複数の観測データに対して同じ独立成分ＩＣ１との内積を計算すると、複数のサンプルに関して同じ独立成分ＩＣ１に対する複数の内積値が得られる。図１（Ｆ）は、複数のサンプルに関する内積値Ｐを横軸にとり、複数のサンプルに含まれている目的成分の既知の含有量Ｃを縦軸に取ってプロットした図である。仮に、内積で使用した独立成分ＩＣ１が目的成分に対応する独立成分である場合には、図１（Ｆ）に示すように、内積値Ｐと、各サンプルの目的成分の含有量Ｃとが強い相関を有する。そこで、図１（Ｃ）で得られた複数の独立成分ＩＣ１，ＩＣ２…のうちで最も強い相関を示す独立成分を、目的成分に対応する独立成分として選択することができる。図１の例では、独立成分ＩＣ１が、検量の目的成分（例えば中性脂肪）に対応する独立成分である。検量線は、図１（Ｆ）のプロットの単回帰式Ｃ＝ｕＰ＋ｖで与えられる直線として表される。なお、内積値Ｐは、各サンプルにおける独立成分ＩＣ１の含有量に比例する値なので、「混合係数」とも呼ぶ。

図２は、検量線を用いた目的成分の検量処理の概要を示す説明図である。検量処理では、図１に示した検量線作成処理で得られた目的成分の独立成分ＩＣ１（図１（Ｅ））と、検量線（図１（Ｆ））とを利用して行われる。検量処理では、まず、目的成分の含有量が未知であるサンプルの観測データを取得する（図２（Ａ））。次に、この観測データに前処理を行う（図２（Ｂ））。この前処理は、検量線の作成時に使用した前処理と同じ処理とすることが好ましい。そして、この前処理後の観測データと、独立成分ＩＣ１（図２（Ｂ））との内積を取ることによって、観測データに関する内積値Ｐを算出する（図２（Ｃ）。この内積値Ｐを、検量線（図２（Ｄ））に適用することによって、目的成分の含有量Ｃを決定することが可能である。なお、図１の検量線作成処理や図２の検量処理の詳細については、後に詳述する。

Ｂ．中性脂肪濃度用の検量線作成方法：
図３は、中性脂肪濃度用の検量線作成処理のフローチャートである。この例では、被検体は人体（例えば手）であり、目的成分は中性脂肪である。なお、図３のステップＴ１１０〜Ｔ１６０の右側には、後述する図５における工程１〜工程５との対応関係が示されている。「工程２」を、「サンプル観測データ取得工程」とも呼び、「工程３」を、「中性脂肪濃度取得工程」とも呼ぶ。また、「工程４」を「推定工程」とも呼び、「工程５」を「回帰式算出工程」とも呼ぶ。

ステップＴ１１０では、人体を分光測定することによってスペクトルデータ（「観測分光データ」又は「測定データ」とも呼ぶ）を取得する。ステップＴ１２０では、ステップＴ１１０と同じ人体から採血を行い、生化学分析を行うことによって、血液中の中性脂肪濃度を正確に求める。これらのステップＴ１１０，Ｔ１２０を複数回繰り返すことによって、複数のサンプルに関するスペクトルデータと中性脂肪濃度とを取得する。

ステップＴ１３０では、図１（Ａ）〜（Ｃ）に示した手順に従って、複数のサンプルのスペクトルデータから複数の独立成分を推定する。ステップＴ１４０では、こうして得られた複数の独立成分の中から、予め選択された特定の波長にピークを有する独立成分を選択する。この特定の波長は、中性脂肪に対応する好ましい波長として予め選択されたものである。特定の波長としては、以下の何れかを使用することが可能である。

図４Ａ及び図４Ｂは、９００〜１１００ｎｍの波長帯のスペクトルデータを検量に用いる場合における好ましい独立成分の例を示している。図４Ａ（１）に示す第１の独立成分ＩＣ１ａは、９３０±α[ｎｍ]にピークを有し、１０６０±α[ｎｍ]にはピークを有さない独立成分である。図４Ａ（２）に示す第２の独立成分ＩＣ１ｂは、９３０±α［ｎｍ］にはピークを有さず、１０６０±α［ｎｍ］にピークを有する独立成分である。図４Ａ（３）に示す第３の独立成分ＩＣ２ｃは、９３０±α［ｎｍ］と１０６０±α［ｎｍ］の両方にピークを有する独立成分である。図４Ａ（４）に示す第４の独立成分ＩＣ２ｃは、９３０±α［ｎｍ］と１０６０±α［ｎｍ］の両方にピークを有し、かつ水に由来する１０００±α［ｎｍ］のピークを有する独立成分である。図４Ｂ（１）に示す第５の独立成分ＩＣ１ｅは、９４５±α[ｎｍ]にピークを有し、１０６０±α[ｎｍ]にはピークを有さない独立成分である。図４Ｂ（２）に示す第６の独立成分ＩＣ１ｆは、９４５±α［ｎｍ］と１０６０±α［ｎｍ］の両方にピークを有する独立成分である。図４Ｂ（３）に示す第７の独立成分ＩＣ２ｇは、９４５±α［ｎｍ］と１０６０±α［ｎｍ］の両方にピークを有し、かつ水に由来する１０００±α［ｎｍ］のピークを有する独立成分である。αは、中性脂肪に特有なピークの範囲を限定するための値である。αの値としては、３０が好ましく、２０が更に好ましい。なお、第３の独立成分ＩＣ１ｃは、第１の独立成分ＩＣ１ａと第２の独立成分ＩＣ１ｂとの和を取ったものに相当する。

９００〜１１００ｎｍの波長帯のスペクトルデータを検量に用いる場合には、図１に示した独立成分分析で得られた複数の独立成分の中から、９３０±３０ｎｍ又は９４５±３０ｎｍの一方と、１０６０±３０ｎｍとの少なくとも１つにピークを有する独立成分が選択される。このような独立成分は、図４Ａ及び図４Ｂに示した７つの独立成分ＩＣ１ａ〜ＩＣ１ｇに近い波形を有する成分波形信号となる。独立成分（成分波形信号）の選択は、例えば、以下のように実行される。

＜ケース１＞
独立成分分析で得られた複数の独立成分の中に、９３０±α[ｎｍ]にピークを有する独立成分は存在するが、１０６０±α[ｎｍ]にピークを有する独立成分が存在しない場合には、９３０±α[ｎｍ]にピークを有する独立成分（図４Ａ（１）の第１の独立成分ＩＣ１ａに近いもの）がステップＴ１４０で選択される。

＜ケース２＞
独立成分分析で得られた複数の独立成分の中に、９３０±α[ｎｍ]にピークを有する独立成分は存在せず、１０６０±α[ｎｍ]にピークを有する独立成分が存在する場合には、１０６０±α[ｎｍ]にピークを有する独立成分（図４Ａ（２）の第２の独立成分ＩＣ１ｂに近いもの）がステップＴ１４０で選択される。

＜ケース３＞
独立成分分析で得られた複数の独立成分の中に、９３０±α[ｎｍ]と１０６０±α[ｎｍ]の両方にピークを有する独立成分が存在する場合には、その独立成分（図４Ａ（３）の第３の独立成分ＩＣ１ｃに近いもの）がステップＴ１４０で選択される。

＜ケース４＞
独立成分分析で得られた複数の独立成分の中に、９３０±α[ｎｍ]と１０００±α[ｎｍ]と１０６０±α[ｎｍ]とにピークを有する独立成分が存在する場合には、９３０±α[ｎｍ] と１０００±α[ｎｍ]と１０６０±α[ｎｍ]とにピークを有する独立成分（図４Ａ（４）の第４の独立成分ＩＣ１ｄに近いもの）がステップＴ１４０で選択される。

＜ケース５＞
独立成分分析で得られた複数の独立成分の中に、９４５±α[ｎｍ]にピークを有する独立成分は存在するが、１０６０±α[ｎｍ]にピークを有する独立成分が存在しない場合には、９４５±α[ｎｍ]にピークを有する独立成分（図４Ｂ（１）の第５の独立成分ＩＣ１eに近いもの）がステップＴ１４０で選択される。

＜ケース６＞
独立成分分析で得られた複数の独立成分の中に、９４５±α[ｎｍ]と１０６０±α[ｎｍ]にピークを有する独立成分が存在する場合には、９４５±α[ｎｍ]と１０６０±α[ｎｍ]にピークを有する独立成分（図４Ｂ（２）の第６の独立成分ＩＣ１ｆに近いもの）がステップＴ１４０で選択される。

＜ケース７＞
独立成分分析で得られた複数の独立成分の中に、９４５±α[ｎｍ]と１０００±α[ｎｍ]と１０６０±α[ｎｍ]とにピークを有する独立成分が存在する場合には、９４５±α[ｎｍ] と１０００±α[ｎｍ]と１０６０±α[ｎｍ]とにピークを有する独立成分（図４Ｂ（３）の第７の独立成分ＩＣ１ｇに近いもの）がステップＴ１４０で選択される。

なお、上記ケース３の場合に、図４Ａに示した３つの独立成分ＩＣ１ａ〜ＩＣ１ｃのそぞれに関して図１（Ｆ）で示したグラフにおける内積値Ｐと成分含有量Ｃ（すなわち中性脂肪濃度）の相関を求め、相関が最も高い１つの独立成分を選択するようにしてもよい。こうすれば、検量精度をより高めることが可能である。上記ケース５の場合も同様である。つまり、図４Ａ及び図４Ｂに示した３つの独立成分ＩＣ１ｂ、ＩＣ１ｅ、ＩＣ１ｆのそぞれに関して図１（Ｆ）で示したグラフにおける内積値Ｐと成分含有量Ｃ（すなわち中性脂肪濃度）の相関を求め、相関が最も高い１つの独立成分を選択するようにしてもよい。

図４Ｃは、１１００〜１２５０ｎｍの波長帯のスペクトルデータを検量に用いる場合における好ましい独立成分の例を示している。図４Ｃ（１）に示す第１の独立成分ＩＣ２ａは、１１５０±α［ｎｍ］にピークを有し、１１８２±α［ｎｍ］にはピークを有さない独立成分である。図４Ｂ（２）に示す第２の独立成分ＩＣ２ｂは、１１５０±α［ｎｍ］にはピークを有さず、１１８２±α［ｎｍ］にピークを有する独立成分である。図４Ａ（３）に示す第３の独立成分ＩＣ２ｃは、１１５０±α［ｎｍ］と１１８２±α［ｎｍ］の両方にピークを有する独立成分である。図４Ａ（４）に示す第４の独立成分ＩＣ２ｄは、１１５０±α［ｎｍ］と１１８２±α［ｎｍ］の両方にピークが重なり、一つのブロードなピークを有する独立成分である。図４Ａ及び図４Ｂの場合と同様に、αの値は３０が好ましく、２０が更に好ましい。また、第３の独立成分ＩＣ２ｃは、第１の独立成分ＩＣ２ａと第２の独立成分ＩＣ２ｂとの和を取ったものに相当する。

１１００〜１２５０ｎｍの波長帯のスペクトルデータを検量に用いる場合には、図１に示した独立成分分析で得られた複数の独立成分の中から、１１５０±α［ｎｍ］と１１８２±α［ｎｍ］の一方又は両方にピークを有する独立成分が選択される。このような独立成分は、図４Ｂに示した４つの独立成分ＩＣ２ａ〜ＩＣ２ｄに近い波形を有するものとなる。独立成分の選択は、９００〜１１００ｎｍの波長帯を用いた場合と同様に行われる。

なお、ステップＴ１４０における独立成分の選択は、検量線作成装置を操作する操作者が手動で行ってもよく、或いは、検量線作成装置内の独立成分選択部が自動的に実行するようにしてもよい。

ところで、上述した波長のピークは、中性脂肪が持つＣＨ基とＣＨ_２基による光吸収に由来するものと考えられる。一般に、９００〜１１００ｎｍの帯域ではＣＨ基とＣＨ_２基の３倍音の光吸収が現れ、１１００〜１２５０ｎｍの帯域では２倍音の光吸収が現れる。また、例えば人体や血液に対する９００〜１１００ｎｍ帯域の独立成分に現れる２つのピークの波長は、生体分子との相互作用によってＣＨ基やＣＨ_２基の単体の吸収ピークよりもやや長波長側にシフトしたものとなる。２つのピークを個別に扱うのは、水の吸収帯の影響を受けることなどの原因によって、２つのピークがそれぞれ別々の独立成分に現れることがあるからである。本願の発明者の実験によれば、人体を観測したスペクトルデータを用いた独立成分では２つのピークが現れることが多い。従って、人体を被検体とする場合には、２つのピークを有する独立成分ＩＣ１ｃやＩＣ１ｆ（又はＩＣ２ｃ）を用いることが好ましい。但し、独立成分の具体的な波形は、人体からスペクトルデータを測定する際の被測定部位（例えば手や足）の位置にも依存する可能性がある。一方、中性脂肪水溶液を被検体とした場合には、ＣＨ_２基由来のピークはより短波長側にシフトしてしまい、ピークが１つになる可能性がある。

こうして独立成分の選択が終了すると、図３のステップＴ１５０において、選択した独立成分を用いて、複数のサンプルに関する内積値Ｐを算出する（図１（Ｄ）〜（Ｅ）参照）。なお、独立成分の選択の際に内積値Ｐを算出している場合には、ステップＴ１５０は省略可能である。ステップＴ１６０では、複数のサンプルに関する内積値Ｐと、ステップＴ１２０で測定した中性脂肪濃度Ｃとの関係から、検量線の回帰式を算出する。

検量線の単回帰式は、例えば以下の式で表される。

ここで、Ｃは中性脂肪濃度、Ｐは内積値（混合係数）、ｕ，ｖは定数である。

なお、上記ステップＴ１４０で選択された独立成分は、図２で説明した検量処理においてもそのまま使用されて検量が実行される。

また、本実施形態では、中性脂肪に対応する独立成分の特徴が明確になっているので、被検体毎にサンプルデータを一から集めなくても検量を行うことができる。例えば、被験者Ａの測定データを用いて決定した独立成分と検量線とを用いて、別の被験者Ｂの中性脂肪濃度の検量を行うことも可能である。

なお、図３のステップＴ１３０で推定された複数の独立成分の中に上述した特有のピークを有するものがない場合には、中性脂肪情報を十分に含んだ測定データが取得できていないものと判断できる。したがって、測定条件や、独立成分分析におけるアルゴリズムパラメータが適切であるか否かを判定するための指標としても有効である。すなわち、推定された複数の独立成分の中に上述した特有のピークを有するものがない場合には、測定条件や、独立成分分析におけるアルゴリズムパラメータの一部を変更して、図３の処理を再度行うことが好ましい。

Ｃ．検量線作成方法：
図５は、本発明の一実施形態としての検量線作成方法を示すフローチャートである。この検量線作成方法は、工程１から工程５までの５つの工程によって構成される。各工程１〜５はこの順に実行される。各工程１〜５について、順に説明する。なお、このＣ節以降では、特に注記しない限り、図３，図４Ａ，図４Ｂで説明した中性脂肪濃度以外の目的成分の検量においても利用できる方法について説明を行うものとする。

［工程１］
工程１は、準備工程であり、作業者により行なわれるものである。作業者は、目的成分の含有量が相違する同一種類の複数のサンプル（例えば、中性脂肪水溶液や人体）を用意（準備）する。本実施例ではｎ個（ｎは２以上の整数）のサンプルを使用する。人体をサンプルとする場合には、被検体としての人体は１つでも良く、異なる日時における同じ人体を複数のサンプルとして使用することが可能である。

［工程２］
工程２は、スペクトルの測定工程であり、作業者により分光計測器を用いて行なわれるものである。作業者は、工程１で用意した複数のサンプルのそれぞれを分光計測器で撮影することにより、各サンプルについての分光反射率のスペクトルを測定する。分光計測器は、被計測体からの光を分光器に通し、分光器から出力されるスペクトルを撮像素子の撮像面で受けることにより、前記スペクトルを測定する周知の機器である。分光反射率のスペクトルと吸光度のスペクトルとの間には、次式で表される関係が成り立つ。

測定された分光反射率のスペクトルは、式（４）を用いて吸光度スペクトルに変換される。吸光度に変換するのは、後述する独立成分分析において分析される混合信号には線形結合が成立する必要があり、ランベルト・ベールの法則から、吸光度について線形結合が成立するためである。したがって、工程２においては、分光反射率スペクトルの代わりに吸光度スペクトルを測定してもよい。測定結果としては、被計測体の波長に対する特性を示す吸光度分布のデータが出力される。この吸光度分布のデータは、スペクトルデータとも呼ぶ。

なお、分光反射率スペクトルや吸光度スペクトルを分光器で測定する代わりに、これらのスペクトルを他の測定値から推定するようにしてもよい。例えば、サンプルをマルチバンドカメラで測定し、得られたマルチバンド画像から分光反射率や吸光度スペクトルを推定するようにしてもよい。このような推定方法としては、例えば、特開２００１−９９７１０号公報に記載された方法などを利用することができる。

［工程３］
工程３は、目的成分の含有量の測定工程であり、作業者により行なわれるものである。作業者は、工程１で用意した複数のサンプルのそれぞれを化学分析して、各サンプルについての目的成分の含有量（例えば中性脂肪量）を測定する。工程１で準備したサンプルにおける目的成分の含有量が既知の場合には、この工程３は省略可能である。

［工程４］
工程４は、混合係数の推定工程であり、典型的にはコンピューターを用いて行なわれるものである。図６は、工程４および後述する工程５で用いられるコンピューター１００とその周辺装置を示す説明図である。コンピューター１００は、分光計測器２００に電気的に接続されている。

コンピューター１００は、コンピュータープログラムを実行することにより種々の処理や制御を行うＣＰＵ１０と、データの退避場所であるメモリー２０（記憶部）と、コンピュータープログラムやデータを保存するハードディスクドライブ３０と、入力インターフェイス５０と、出力インターフェイス６０とを備えた周知な装置である。

図７は、工程４及び工程５で使用する装置の機能ブロック図である。この装置４００は、サンプル観測データ取得部４１０と、サンプル目的成分量取得部４２０と、混合係数推定部４３０と、回帰式算出部４４０とを有する。混合係数推定部４３２は、独立成分行列算出部４３２、推定混合行列算出部４３４、混合係数選択部４３６、及び、独立成分選択部４３８を含んでいる。なお、サンプル観測データ取得部４１０およびサンプル目的成分量取得部４２０は、例えば図６のＣＰＵ１０が入力Ｉ／Ｆ５０とメモリー２０と協働して実現される。混合係数推定部４３０、独立成分行列算出部４３２、推定混合行列算出部４３４、混合係数選択部４３６、及び、独立成分選択部４３８は、例えば図６のＣＰＵ１０がメモリー２０と協働して実現される。また、回帰式算出部４４０は、例えば図６のＣＰＵ１０がメモリー２０と協働して実現される。なお、これらの各部は、図６に示したコンピューター以外の他の具体的な装置やハードウェア回路によっても実現可能である。

独立成分選択部４３８は、図３のステップＴ１４０で説明した独立成分の選択処理を実行する機能を有する。独立成分選択部４３８は、検量線作成装置４００を操作する操作者が入力部（図示省略）を用いて選択指示を入力すると、この入力に応じて独立成分の選択を実行するように構成されていてもよい。或いは、操作者の指示を必要とせずに独立成分選択部４３８が自動的に独立成分の選択を実行するように構成されていてもよい。なお、図３の手順で独立成分の選択を行う場合には、図７に示した推定混合号列算出部４３４と混合係数選択部４３６は省略可能である。以下では、混合係数推定部４３０を単に「推定部４３０」とも呼び、また、サンプル目的成分量取得部４２０を「中性脂肪濃度取得部」とも呼ぶ。

図８は、独立成分行列算出部４３２の内部構成の一例を示す機能ブロック図である。独立成分行列算出部４３２は、第１前処理部４５０と、第２前処理部４６０と、独立成分分析処理部４７０とを有している。これらの３つの処理部４５０，４６０，４７０は、この順番に処理対象データ（本実施形態では吸光度スペクトル）を処理することによって、独立成分行列（後述）を求める。これらの各部の処理内容については後述する。

図６に示した分光計測器２００は、工程２で使用されたものである。コンピューター１００は、工程２で分光計測器２００により測定された分光分布から得られた吸光度スペクトルを、スペクトルデータとして入力Ｉ／Ｆ５０を介して取得する（図７のサンプル観測データ取得部４１０に対応）。また、コンピューター１００は、工程３で測定された目的成分の含有量を、作業者によるキーボードの操作等によって取得する（図７のサンプル目的成分量取得部４２０に対応）。

上記のスペクトルデータと目的成分含有量の取得の結果、コンピューター１００のハードディスクドライブ３０には、スペクトルデータと目的成分含有量とを含むデータセット（以下、「測定データセット」と呼ぶ）ＤＳ１が保存される。

図９は、ハードディスクドライブ３０に保存された測定データセットＤＳ１を模式的に示す説明図である。図示するように、測定データセットＤＳ１は、工程１で用意した複数のサンプルを識別するためのサンプル番号Ｂ1，Ｂ2，…，Ｂnと、各サンプルについての目的成分含有量Ｃ1，Ｃ2，…，Ｃnと、各サンプルについてスペクトルデータＸ1，Ｘ2，…，Ｘnとを含むデータ構造である。測定データセットＤＳ１において、目的成分含有量Ｃ1，Ｃ2，…，ＣnおよびスペクトルデータＸ1，Ｘ2，…，Ｘnは、いずれのサンプルについてのものであるかが判るように、サンプル番号Ｂ1，Ｂ2，…，Ｂnと対応づけがなされている。

ＣＰＵ１０は、ハードディスクドライブ３０に格納された所定のプログラムをメモリー２０にロードし、そのプログラムを実行することで、工程４の作業である、混合係数を推定する処理を行う。ここで、前記所定のプログラムは、外部からインターネット等のネットワークを用いてダウンロードする構成とすることもできる。工程４において、ＣＰＵ１０は図７の混合係数推定部４３０として機能する。

図１０は、ＣＰＵ１０で実行される混合係数推定処理を示すフローチャートである。処理が開始されると、ＣＰＵ１０は、まず、独立成分分析を行う（ステップＳ１１０）。

独立成分分析(Independent Component Analysis, ICA) は、多次元信号解析法の１つであり、独立な信号が重なり合った混合信号をいくつかの異なる条件で観測し、それを基に独立な原信号を分離する技術である。独立成分分析を用いれば、工程２により得られたスペクトルデータを、目的成分を始めとするｍ個の独立成分（未知）が混合されたものと捕らえることで、工程２により得られたスペクトルデータ（観測データ）から独立成分のスペクトルを推定することが可能となる。

本実施形態において、独立成分分析は、図８に示した３つの処理部４５０，４６０，４７０がこの順に処理を行うことによって実行される。第１前処理部４５０は、標準正規変量変換（ＳＮＶ）４５２と、零空間射影法（ＰＮＳ）４５４のうちの一方又は両方を用いた前処理を実行することが可能である。ＳＮＶ４５２は、処理対象データの平均値を減算し、その標準偏差で除算することによって、平均値が０で標準偏差が１である正規化されたデータを得る処理である。ＰＮＳ４５４は、処理対象データに含まれるベースライン変動を取り除くための処理である。スペクトルの測定では、さまざまな要因により、測定データごとにデータの平均値が上下するなどのベースライン変動と呼ばれるデータ間のばらつきが発生する。そのため、独立成分分析処理を行う前に、この変動要因を取り除くことが好ましい。ＰＮＳは、任意のベースライン変動を取り除くことの出来る前処理として使用することが可能である。なお、ＰＮＳについては、例えば、Zeng-Ping Chen, Julian Morris, and Elaine Martin, “Extracting Chemical Information from Spectral Data with Multiplicative Light Scattering Effects by Optical Path-Length Estimation and Correction”, 2006に説明されている。

なお、図５の工程２で得られたスペクトルデータに対してＳＮＶ４５２を行う場合にはＰＮＳ４５４による処理を行う必要は無い。一方、ＰＮＳ４５４による処理を行う場合には、その後に何らかの正規化処理（例えばＳＮＶ４５２）を行うことが好ましい。

なお、第１前処理として、ＳＮＶやＰＮＳ以外の処理を行うようにしても良い。第１前処理では何らかの正規化処理を行うことが好ましいが、正規化処理を省略してもよい。以下では、第１前処理部４５０を「正規化処理部」とも呼ぶ。これらの２つの処理４５２，４５４の内容については、更に後述する。なお、独立成分行列算出部４３２に与えられる処理対象データが正規化済みのデータである場合には、第１前処理を省略することも可能である。

第２前処理部４６０は、主成分分析（ＰＣＡ）４６２と、因子分析（ＦＡ）４６４のうちのいずれか一方を用いた前処理を実行することが可能である。また、第２前処理として、ＰＣＡやＦＡ以外の処理を用いるようにしても良い。以下では、第２前処理部４６０を「白色化処理部」とも呼ぶ。一般的なＩＣＡの手法では、第２前処理として、処理対象データの次元圧縮と、無相関化を行う。この第２前処理によって、ＩＣＡで求めるべき変換行列が直交変換行列に限定されるため、ＩＣＡの計算量を削減できる。このような第２前処理を「白色化」とよび、多くの場合にＰＣＡが用いられる。しかし、ＰＣＡは、処理対象データにランダムノイズが含まれる場合、その影響を受けて結果に誤差が生じる場合がある。そこで、ランダムノイズの影響を低減するために、ＰＣＡの代わりに、ノイズに対するロバスト性を持つＦＡを用いて白色化を行うことが好ましい。図８の第２前処理部４６０は、ＰＣＡとＦＡのいずれか一方を選択して白色化を実行することが可能である。これらの２つの処理４６２，４６４の内容については、更に後述する。なお、白色化処理は省略しても良い。

独立成分分析処理部（ＩＣＡ処理部）４７０は、第１前処理と第２前処理が行われたスペクトルデータに対して、ＩＣＡを実行することによって、独立成分のスペクトルを推定する。ＩＣＡ処理部４７０は、独立性指標として尖度を用いる第１処理４７２と、独立性指標としてβダイバージェンスを用いる第２処理４７４のうちのいずれか一方を用いた分析を実行することが可能である。ＩＣＡは、一般に、独立成分の分離のための指標として、分離したデータ同士の独立性を表す高次統計量を独立性指標として用いる。尖度は、典型的な独立性指標である。しかし、処理対象データにスパイクノイズのような外れ値が入っている場合には、その外れ値も含めた統計量が独立性指標として計算されてしまう。このため、処理対象データについての本来の統計量と、算出された統計量との間に誤差が生じ、分離精度の低下を引き起こす場合がある。そこで、処理対象データ中の外れ値からの影響を低減するために、その影響を受けにくい独立性指標を使用することが好ましい。この様な特性を持つ独立性指標として、βダイバージェンスを使用可能である。尖度とβダイバージェンスの内容については、更に後述する。なお、ＩＣＡの独立性指標としては、尖度やβダイバージェンス以外の指標を利用するようにしても良い。

独立成分分析の典型的な処理内容について、次に詳しく説明する。ｍ個の未知成分（ソース）のスペクトルＳ（以下、このスペクトルを単に「未知成分」と呼ぶこともある）が下記の式（５）のベクトルで与えられたとし、工程２により得られたｎ個のスペクトルデータＸが下記の式（６）のベクトルで与えられたとする。なお、式（５）に含まれる各要素（Ｓ1、Ｓ2、・・・Ｓm）は、それぞれがベクトル（スペクトル）となっている。すなわち、例えば要素Ｓ1は式（７）のように表される。式（６）に含まれる要素（Ｘ1、Ｘ2、・・・、Ｘn）もベクトルであり、例えば要素Ｘｊは式（８）のように表される。要素Ｘｊの添え字のｊはスペクトルを測定した波長帯の数である。なお、未知成分のスペクトルＳの要素数ｍは、１以上の整数であり、サンプルの種類に応じて予め経験的又は実験的に決められている。

各未知成分は統計的に独立であるとする。これら未知成分ＳとこれらスペクトルデータＸとの間に、次式の関係が成立する。

式（９）におけるＡは混合行列であり、次式（１０）で示すこともできる。なお、ここでも“Ａ”の文字は、式（１０）に示すように太字で示す必要があるが、明細書の使用文字の制限からここでは普通文字で表す。以下、行列を表す他の太字について、同様に普通文字で表すものとする。

混合行列Ａに含まれる混合係数ａi jは、未知成分Ｓj（j＝１〜ｍ）が、観測データであるスペクトルデータＸi（i＝１〜ｎ）へ寄与する度合いを示す。

混合行列Ａが既知の場合、未知成分Ｓの最小二乗解は、Ａの擬似逆行列Ａ+を用いてＡ+・Ｘと簡単に求めることができるが、本実施形態の場合、混合行列Ａも未知なので、観測データＸのみから、未知成分Ｓと混合行列Ａを推定しなければならない。すなわち、次式（１１）に示すように、観測データＸのみから、ｍ×ｎの分離行列Ｗを用いて、独立成分のスペクトルを示す行列（以下、「独立成分行列」と呼ぶ）Ｙを算出する。次式（１１）における分離行列Ｗを求めるアルゴリズムとしては、Infomax、FastICA (Fast Independent Component Analysis)、JADE (Joint Approximate Diagonalization of Eigenmatrices) など種々のものを採用することが可能である。

独立成分行列Ｙは未知成分Ｓの推定値に該当する。よって、下記の式（１２）を得ることができ、式（１２）を変形して下記の式（１３）を得ることができる。

式（１３）で得られる推定混合行列∧Ａ（明細書の使用文字の制限からこのように表記しただけであり、実際は式（１３）の左辺の符号付き文字を意味する。他の文字も同じ）は、次式で示すこともできる。

図１０のステップＳ１１０では、ＣＰＵ１０は、上述した分離行列Ｗを求める処理までを行う。詳しくは、工程２で得られハードディスクドライブ３０に予め保存したサンプル毎のスペクトルデータＸを入力として、この入力に基づいて、前述したInfomax、FastICA、JADE などのいずれかのアルゴリズムを用いて、分離行列Ｗを求める。なお、前述した図８に示すように、独立成分分析の前処理として、第１前処理部４５０による正規化処理と、第２前処理部４６０による白色化処理とを行うことが好ましい。

ステップＳ１１０の実行後、ＣＰＵ１０は、その分離行列Ｗと、工程２で得られハードディスクドライブ３０に予め保存したサンプル毎のスペクトルデータＸとに基づいて、独立成分行列Ｙを算出する処理を行う（ステップＳ１２０）。この算出処理は、前述した式（１１）に従う演算を行うものである。ステップＳ１１０、Ｓ１２０の処理において、ＣＰＵ１０は図７の独立成分行列算出部４３２として機能する。

続いて、ＣＰＵ１０は、前記ハードディスクドライブ３０に予め保存したサンプル毎のスペクトルデータＸと、ステップＳ１２０で算出した独立成分行列Ｙとに基づいて、推定混合行列∧Ａを算出する処理を行う（ステップＳ１３０）。この算出処理は、前述した式（１３）に従う演算を行うものである。

図１１は、推定混合行列∧Ａを説明するための説明図である。この表ＴＢは、縦方向に各サンプル番号Ｂ1，Ｂ2，…，Ｂnをとり、横方向に独立成分行列Ｙの各要素（以下、「独立成分要素」と呼ぶ）Ｙ1，Ｙ2，…，Ｙmをとったものである。サンプル番号Ｂｉ（ｉ＝１〜ｎ）と独立成分要素Ｙj（j＝１〜ｍ）から定まる表ＴＢ中の要素は、推定混合行列∧Ａの要素∧ａi j（式（１４）参照）と同一である。この表ＴＢからも、推定混合行列∧Ａの要素∧ａi jは、サンプルのそれぞれにおける独立成分要素Ｙ1，Ｙ2，…，Ｙmの比率を表したものであることがわかる。図１１に例示した目的成分順位ｋについては後述する。ステップＳ１３０の処理において、ＣＰＵ１０は図７の推定混合行列算出部４３４として機能する。

ステップＳ１３０までの処理によって、推定混合行列∧Ａが得られる。すなわち、推定混合行列∧Ａの要素（推定混合係数）∧ａi jが得られる。なお、推定混合係数∧ａi jは、図１の例では、図１（Ｄ）〜（Ｆ）で算出される内積値Ｐに相当する。その後、ステップＳ１６０又はステップＳ１４０に進む。

ステップＳ１６０では、図３のステップＴ１４０で説明した方法に従って、独立成分が選択される。この選択は、図７の独立成分選択部４３８によって行われる。ステップＴ１７０では、選択した独立成分を用いて、複数のサンプルに関する内積値Ｐを算出し、これらの内積値Ｐを要素とする混合係数ベクトル∧αを求める。この混合係数ベクトル∧αは、図１１に示す１列分の混合係数に相当する。こうして、図１０の処理も終了する。なお、ステップＳ１６０における独立成分の選択が実行されない場合には、以下で説明するステップＳ１４０，Ｓ１５０によって独立成分が選択される。ステップＳ１６０〜Ｔ１７０の処理とステップＳ１４０〜Ｓ１５０の処理のいずれが実行されるかについては、操作者が任意に指定できるように、検量線作成装置４００（図７）が構成されていることが好ましい。

ステップＳ１４０では、ＣＰＵ１０は、工程３で測定された目的成分含有量Ｃ1，Ｃ2，…，Ｃnと、ステップＳ１３０で算出された推定混合行列∧Ａに含まれる各列の成分（以下、混合係数ベクトル∧αと呼ぶ）との間の相関（類似性の度合い）を求める。詳しくは、目的成分含有有量Ｃ（Ｃ1，Ｃ2，…，Ｃn）と第１列目の混合係数ベクトル∧α1（∧ａ11，∧ａ21，…，∧ａn1）との相関を求め、次いで、目的成分含有量Ｃ（Ｃ1，Ｃ2，…，Ｃn）と第２列目の混合係数ベクトル∧α2（∧ａ12，∧ａ22，…，∧ａn2）との相関を求め、こうして順次各列について目的成分含有量Ｃに対する相関を求める。

上記相関の高さを示す指標としては、次式に従う相関係数Ｒを使用できる。この相関係数Ｒは、ピアソンの積率相関係数と呼ばれるものである。

図１０のステップＳ１４０の結果、独立成分（独立成分スペクトル）Ｙｊ毎の相関係数Ｒj（ｊ＝１，２，……，ｍ）が得られる。その後、ＣＰＵ１０は、ステップＳＳ１４０で得られた相関係数Ｒjの中から最も相関が高いもの、すなわち値が１に近いものを特定する。そして、最も高い相関係数Ｒが得られた列ベクトル∧αを、推定混合行列∧Ａの中から選択する（ステップＳ１５０）。

ステップＳ１５０における選択は、図１１の表ＴＢでいえば、複数の列の中から一の列を選択することである。この選択された列の要素が、目的成分に対応する独立成分の混合係数である。前記選択の結果、混合係数ベクトル∧αk（∧ａ1k，∧ａ2k，…，∧ａnk）が得られる。ここで、ｋは１〜ｍのいずれかの整数をとる。なお、このｋの値を、何番目の独立成分が目的成分に当たるかを示す目的成分順位としてメモリー２０に一時的に保存してもよい。この混合係数ベクトル∧αkに含まれる要素∧ａ1k，∧ａ2k，…，∧ａnkが「目的成分に対応する混合係数」に相当する。なお、図１１の例では、目的成分順位ｋ＝２が、独立成分Ｙ２に対応する混合係数ベクトル∧α2＝（∧ａ12，∧ａ22，…，∧ａn2）を示している。なお、本明細書において、「順位」という語句は、「行列内の位置を示す値」という意味で使用されている。ステップＳ１４０、Ｓ１５０の処理において、ＣＰＵ１０は図７の推定係数選択部４３６として機能する。ステップＳ１５０の実行後、ＣＰＵは、この混合係数の算出処理を終える。この結果、工程４が完了し、その後、工程５に進む。

［工程５］
工程５は、回帰式の算出工程であり、工程４を実行した時と同じくコンピューター１００を用いて行なわれるものである。工程５では、コンピューター１００は、検量線の回帰式を算出する処理を実行する。なお、工程５は、工程４までのデータを別のコンピューターや装置に移して実行してもよい。

図１２は、コンピューター１００のＣＰＵ１０で実行される回帰式の算出処理を示すフローチャートである。処理が開始されると、ＣＰＵ１０は、まず、工程３で測定された目的成分含有量Ｃ（Ｃ1，Ｃ2，…，Ｃn)と、ステップＳ１５０（又はステップＳ１７０）で得られた混合係数ベクトル∧αk（∧ａ1k，∧ａ2k，…，∧ａnk）とに基づいて、回帰式を算出する（ステップＳ２１０）。この回帰式は、次式（１６）にて表すことができる。ステップＳ２１０では、式（１６）における定数ｕ，ｖが求められることになる。

ここで、Ｃは目的成分含有量、Ｐは測定データと独立成分との内積値、ｕ，ｖは定数である。

ステップＳ２１０の実行後、ＣＰＵ１０は、ステップＳ２１０で求められた回帰式の定数ｕ，ｖと、ステップＳ１５０で決定された目的成分順位ｋ（図１１）に対応する独立成分Ｙｋ（又はステップＳ１６０で選択された独立成分）とを、検量用データセットＤＳ２としてハードディスクドライブ３０に保存する（ステップＳ２２０）。その後、ＣＰＵ１０は、「リターン」に抜けて、この回帰式の算出処理を一旦終了する。この結果、検量線の回帰線を求めることができ、図５に示した検量線作成方法も終了する。ステップＳ２１０、Ｓ２２０の処理において、ＣＰＵ１０は図７の回帰式算出部４４０として機能する。

Ｄ．目的成分の検量方法：
目的成分の検量方法について、次に説明する。被検体は、検量線を作成したときに用いたサンプルと同じ成分で構成されるものとする。具体的には、目的成分の検量方法は、コンピューターを用いて行なわれるものである。なお、ここでのコンピューターは、検量線を作成する際に用いたコンピューター１００であってもよいし、他のコンピューターであってもよい。

図１３は、目的成分の検量を行う際に使用する装置の機能ブロック図である。この中性脂肪濃度検量装置５００は、被検体観測データ取得部５１０と、検量用データ取得部５２０と、混合係数算出部５３０と、目的成分量算出部５４０と、不揮発性記憶装置５５０とを有する。混合係数算出部５３０は、前処理部５３２を含んでいる。この前処理部５３２は、図８の第１前処理部４５０と第２前処理部４６０の両者の機能を有している。混合計算算出部５３０は、図２（Ａ）〜（Ｃ）で説明した内積演算を行う機能を有しているので、「内積演算部」と呼ぶことも可能である。被検体観測データ取得部５１０は、例えば図６のＣＰＵ１０が入力Ｉ／Ｆ５０とメモリー２０と協働して実現される。検量用データ取得部５２０は、例えば図６のＣＰＵ１０がメモリー２０とハードディスクドライブ３０と協働して実現される。混合係数算出部５３０、および目的成分量算出部５４０は、例えば図６のＣＰＵ１０がメモリー２０と協働して実現される。不揮発性記憶装置５５０には、検量用データセットＤＳ２（独立成分及び回帰式の定数ｕ，ｖ）が格納されている。なお、図１３の装置は、図６のコンピュータとは異なる別の装置又は電子機器として実装してもよい。この場合に、図１３の装置又はこれを備える電子機器は、分光計測器を有することが好ましい。

図１４は、コンピューター１００のＣＰＵ１０で実行する目的成分検量処理を示すフローチャートである。この目的成分検量処理は、ＣＰＵ１０がハードディスクドライブ３０に格納された所定のプログラムをメモリー２０にロードし、そのプログラムを実行することで、実現される。まず、ＣＰＵ１０は、被検体を分光計測器で撮影する処理を行う（ステップＳ３１０）。ステップＳ３１０による撮影は工程２と同様にして行うことができ、この結果、被検体の吸光度スペクトルＸpが得られる。検量処理で使用する分光計測器は、誤差を抑制するため検量線の作成に使用した分光計測器と同一機種であることが好ましい。誤差をさらに抑制するためには、同一の機体であることがより好ましい。なお、図５の工程２と同様に、分光反射率スペクトルや吸光度スペクトルを分光器で測定する代わりに、これらのスペクトルを他の測定値から推定するようにしてもよい。一の被検体を一回撮像した場合に得られる被検体の吸光度のスペクトルＸpは、次式のようにベクトルで表される。

ステップＳ３１０の処理において、ＣＰＵ１０は図１３の被検体観測データ取得部５１０として機能する。次いで、ＣＰＵ１０は、ハードディスクドライブ３０（図１３の不揮発性記憶装置５５０）から検量用データセットＤＳ２を取得し、メモリー２０に格納する（ステップＳ３２０）。ステップＳ３２０の処理において、ＣＰＵ１０は図１３の検量用データ取得部５２０として機能する。

ステップＳ３２０の実行後、ステップＳ３１０で得られた被検体の観測データ（吸光度スペクトルＸp）に対して前処理を実行する（ステップＳ３３０）。この前処理としては、検量線の作成時に図５の工程４（より具体的には図１０のステップＳ１１０）において使用された前処理（すなわち第１前処理部４５０による正規化処理及び第２前処理部４６０による白色化処理）と同じ処理を実行することが好ましい。

その後、ＣＰＵ１０は、検量用データセットＤＳ２に含まれる独立成分とステップＳ３３０で得られた前処理済みのスペクトル（前処理済みの観測データ）との内積値Ｐを求める（ステップＳ３４０）。ステップＳ３４０の処理は、前述した図２（Ｂ），（Ｃ）の処理に相当する。なお、この内積値Ｐは、検量線の作成時において、図１０のステップＳ１３０で算出した混合係数に相当する。従って、内積値Ｐを「混合係数」とも呼ぶ。

ステップＳ３３０，Ｓ３４０の処理において、ＣＰＵ１０は図１３の混合係数算出部５３０として機能することになる。

続いて、ＣＰＵ１０は、ハードディスクドライブ３０（図１３の不揮発性記憶装置５５０）から検量用データセットＤＳ２に含まれる回帰式の定数ｕ，ｖを読み出し、この定数ｕ，ｖとステップＳ３４０で得られた内積値Ｐとを前述した式（１６）の右辺に代入することで、目的成分の含有量Ｃを求める（ステップＳ３５０）。この際、必要に応じて定数ｕ，ｖを調整してもよい。含有量Ｃは、例えば、被検体の単位容積又は単位質量当たり（例えば、１ｄＬ当たり、又は、１００グラム当たり）の目的成分の質量として求められる。ステップＳ３５０の処理において、ＣＰＵ１０は図１３の目的成分量算出部５４０として機能することになる。その後、「リターン」に抜けて、この目的成分検量処理を終了する。

なお、本実施形態では、ステップＳ３５０によって求められた含有量Ｃを被検体の目的成分の含有量としていたが、これに替えて、ステップＳ３５０によって求められた含有量Ｃを、ステップＳ３３０による正規化で用いた正規化係数で補正して、この補正後の値を求めるべき含有量としてもよい。具体的には、含有量Ｃに標準偏差を乗算することによって、含有量の絶対値（グラム）を求めるようにしてもよい。この構成によれば、目的成分の種類によっては、含有量Ｃをより精度の高いものとすることが可能となる。

以上の検量方法によれば、被検体の実測値である一のスペクトルから目的成分の含有量を高精度に求めることができる。

［実験結果］
図１５は、本願の発明者による実験結果である。図１５（Ａ）は、独立成分回帰分析を用いた結果を示し、図１５（Ｂ）は、ＰＬＳ回帰分析を用いた結果を示す。図１５（Ａ）及び図１５（Ｂ）において、横軸を真の中性脂肪濃度（ｍｇ／ｄｌ）とし、縦軸を各分析により得られた検量線濃度（ｍｇ／ｄｌ）とする。本実験では、人体を被検体とした観測分光データとそのデータを取得した人体の血液中における中性脂肪濃度とのセットを３０セット用意し、検量回帰式作成に１５セット用い、残りの１５セットを検量精度（ＳＥＰ）算出に用いた。検量線作成処理を行った結果、図４Ａ（１）に示した第１の独立成分ＩＣ１ａを用いて、図１５（Ａ）に示すように、検量精度ＳＥＰ（予測標準偏差）として１８．６５ｍｇ／ｄＬが得られた。一方、図１５（Ｂ）に示すように、ＰＬＳ回帰分析の中性脂肪濃度の検量では、検量精度ＳＥＰとして２１．５６ｍｇ／ｄＬが得られた。この結果から、独立成分回帰分析を用いた場合、ＰＬＳ回帰分析を用いた場合に比べて高い検量精度を達成することが可能であることがわかった。

Ｅ．各種のアルゴリズムとその影響：
以下では、図８に示した第１前処理部４５０と、第２前処理部４６０と、独立成分分析処理部４７０とにおいて利用される各種のアルゴリズムについて順次説明する。

Ｅ−１．第１前処理（ＳＮＶ／ＰＮＳを利用した正規化処理）：
第１前処理部４５０が行う第１前処理としては、ＳＮＶ（標準正規変量変換）とＰＮＳ（零空間射影法）を利用可能である。

ＳＮＶは、以下の式で与えられる。

ここで、ｚは処理後のデータ、ｘは処理対象データ（本実施例では吸光度スペクトル）、ｘaveは処理対象データｘの平均値、σは処理対象データｘの標準偏差である。標準正規変量変換の結果、平均値が０で標準偏差が１である正規化されたデータｚが得られる。

ＰＮＳを行うと、処理対象データに含まれるベースライン変動を低減することが可能である。処理対象データ（本実施形態では吸光度スペクトル）の測定では、さまざまな要因により、測定データごとにデータの平均値が上下するなどのベースライン変動と呼ばれるデータ間のばらつきが発生する。そのため、ＩＣＡ（独立成分分析）を行う前に、この変動要因を取り除くことが好ましい。ＰＮＳは、処理対象データのベースライン変動を低減することのできる前処理として使用することが可能である。特に、赤外領域を含む吸収光スペクトル又は反射光スペクトルの測定データについては、このようなベースライン変動が多いので、ＰＮＳを適用する利点が大きい。以下では、測定で得られるデータ（単に「測定データｘ」とも呼ぶ）に含まれるベースライン変動を、ＰＮＳによって取り除く原理について説明する。また、典型的な例として、測定データが赤外領域を含む吸収光スペクトル又は反射光スペクトルである場合について説明を行う。但し、他の種類の測定データ（例えば音声データなど）についても同様にＰＮＳを適用可能である。

一般に、理想系において、測定データｘ（処理対象データｘ）は、ｍ個（ｍは２以上の整数）の独立成分ｓｉ（ｉ＝１〜ｍ）とそれぞれの混合比ｃｉを用いて、以下の式で表わされる。

ここで、Ａは混合比ｃｉで形成される行列（混合行列）である。

ＩＣＡ（独立成分分析）においても、このモデルを前提として処理が実行される。しかし、実際の測定データには、さまざまな変動要因（試料の状態や測定環境の変化など）が存在している。そこで、それらを考慮したモデルとして、以下の式によって測定データｘを表現するモデルが考えられる。

ここで、ｂはスペクトルの振幅方向の変動量を表すパラメーター、ａは定数ベースライン変動Ｅ（「平均値変動」とも呼ぶ）の量を表すパラメーター、ｂ１…ｂｇは波長に依存するｇ個（ｇは１以上の整数）の変動ｆ１（λ）〜ｆｇ（λ）の量を表すパラメーターであり、εはそれ以外の変動成分である。また、定数ベースライン変動Ｅは、Ｅ＝｛１，１，１，…１｝T（右肩のＴは転置を示す）で与えられ、そのデータ長が測定データｘのデータ長Ｎ（波長域の区分数）と等しい定数ベクトルである。波長を表す変数λとしては、１からＮまでのＮ個の整数が使用される。すなわち、この変数λは、測定データｘのデータ長Ｎ（Ｎは２以上の整数）の序数に相当する。このとき、波長に依存する変動ｆ１（λ）…ｆｇ（λ）は、ｆ１（λ）＝｛ｆ１（１）, ｆ１（２）, …ｆ１（Ｎ）｝T，…，ｆｇ（λ）＝｛ｆｇ（１）, ｆｇ（２）, …ｆｇ（Ｎ）｝Tで与えられる。これらの変動は、ＩＣＡや検量における誤差要因となるため、事前に取り除くことが望ましい。

関数ｆ（λ）としては、λの値が１からＮまでの範囲においてλの増加に応じて関数ｆ（λ）の値が単調に増加する１変数関数を使用することが好ましい。零空間射影法において、指数αが整数であるλのべき乗関数λα以外の他の関数を使用することによって、測定データに含まれる変動をより低減することが可能である。

好ましい関数ｆ（λ）の関数形とその個数ｇと決定する方法としては、実験的なトライアンドエラーを採用してもよく、或いは、既存のパラメーター推定アルゴリズム（例えばＥＭ（期待値最大化法）アルゴリズム）を用いてもよい。

ＰＮＳでは、上述したそれぞれのベースライン変動成分Ｅ，ｆ１（λ）〜ｆｇ（λ）からなる空間を考え、測定データｘを、それら変動成分を含まない空間（零空間）に射影することで、ベースライン変動成分Ｅ，ｆ１（λ）〜ｆｇ（λ）を低減したデータを得ることができる。具体的な演算として、ＰＮＳによる処理後のデータｚは、以下の式で算出される。

ここで、Ｐ＋はＰの擬似逆行列である。ｋｉは、式（２０）の構成成分ｓｉを、変動成分を含まない零空間に射影したものである。また、ε*は、式（２０）の変動成分εを零空間に射影したものである。

なお、ＰＮＳの処理後に、正規化（例えばＳＮＶ）を行えば、式（２０）におけるスペクトルの振幅方向の変動量ｂの影響も取り除くことができる。

このようなＰＮＳによる前処理を行ったデータに対してＩＣＡを行うと、得られる独立成分は、式（２１）の成分ｋｉの推定値となり、真の構成成分ｓｉとは異なるものとなる。しかし、混合比ｃｉは、元の式（２０）における値から変化しないため、混合比ｃｉを使用した検量処理（図２，図１４）には影響しない。このように、ＩＣＡの前処理としてＰＮＳを実行すると、ＩＣＡによって真の構成成分ｓｉを得ることはできないので、通常は、ＩＣＡの前処理にＰＮＳを適用するという発想は生じ得ない。一方、本実施形態では、ＰＮＳをＩＣＡの前処理として行っても検量処理には影響しないので、ＰＮＳを前処理として実行すれば、より精度良く検量を行うことが可能である。

なお、ＰＮＳの詳細については、例えば、Zeng-Ping Chen, Julian Morris, and Elaine Martin, “Extracting Chemical Information from Spectral Data with Multiplicative Light Scattering Effects by Optical Path-Length Estimation and Correction”, 2006に説明されている。

Ｅ−２．第２前処理（ＰＣＡ／ＦＡを利用した白色化処理）：
第２前処理部４６０が行う第２前処理としては、ＰＣＡ（主成分分析）とＦＡ（因子分析）を利用可能である。

一般的なＩＣＡの手法では、前処理として、処理対象データの次元圧縮と、無相関化を行う。この前処理によって、ＩＣＡで求めるべき変換行列が直交変換行列に限定されるため、ＩＣＡの計算量を削減できる。このような前処理を「白色化」とよび、多くの場合にＰＣＡが用いられる。ＰＣＡを用いた白色化については、例えばAapo Hyvarinen, Juha Karhumen, Erkki Oja, "Independent Comonent Analysis", 2001, John Wiley &Sons, Inc.（「独立成分分析」、２００５年２月、東京電気大学出版部発行）の第６章に詳述されている。

しかしながら、ＰＣＡでは、処理対象データにランダムノイズが含まれる場合、そのランダムノイズの影響を受けて、処理結果に誤差が生じる場合がある。そこで、ランダムノイズの影響を低減するために、ＰＣＡの代わりに、ノイズに対するロバスト性を持つＦＡ（因子分析）を用いて白色化を行うことが好ましい。以下では、ＦＡによる白色化の原理を説明する。

前述したように、一般にＩＣＡでは、処理対象データｘを、構成成分ｓｉの線形和として表す線形混合モデル（上記式（１９））を仮定し、混合比ｃｉと構成成分ｓｉとを求める。しかし、実際のデータには、構成成分ｓｉ以外のランダムノイズが付加されていることが多い。そこで、ランダムノイズを考慮したモデルとして、以下の式によって測定データｘを表現するモデルが考えられる。

ここで、ρはランダムノイズである。

そして、このノイズ混合モデルを考慮した白色化を行い、その後、ＩＣＡを行って混合行列Ａと独立成分ｓｉの推定を得ることができる。

本実施形態のＦＡでは、独立成分ｓｉとランダムノイズρとがそれぞれ正規分布Ｎ（０，Ｉｍ），Ｎ（０，Σ）に従うと仮定する。なお、一般に知られているように、正規分布Ｎ（ｘ１，ｘ２）の第１のパラメーターｘ１は期待値を示し、第２のパラメーターｘ２は標準偏差を示す。このとき、処理対象データｘは、正規分布に従う変数の線形和となるので、処理対象データｘもやはり正規分布に従う。ここで、処理対象データｘの共分散行列をＶ［ｘ］とすると、処理対象データｘが従う正規分布はＮ（０，Ｖ［ｘ］）と表現できる。このとき、処理対象データｘの共分散行列Ｖ［ｘ］に関する尤度関数を、下記の手順で計算できる。

まず、独立成分ｓｉが互いに直交すると仮定すると、処理対象データｘの共分散行列Ｖ［ｘ］は以下の式で計算される。

ここで、Σはノイズρの共分散行列である。

このように、共分散行列Ｖ［ｘ］は、混合行列Ａとノイズの共分散行列Σとで表すことができる。このとき、対数尤度関数Ｌ（Ａ，Σ）は以下の式で与えられる。

ここで、ｎはデータｘのデータ個数、ｍは独立成分の個数であり、演算子ｔｒは行列のトレース（対角成分の和）を示し、演算子ｄｅｔは行列式を示す。また、Ｃはデータｘから標本計算で求められる標本共分散行列であり、以下の式で計算される。

上記式（２４）の対数尤度関数Ｌ（Ａ，Σ）を用いた最尤法により、混合行列Ａとノイズの共分散行列Σとを求めることができる。この混合行列Ａとしては、上記式（２２）のランダムノイズρの影響がほとんど無いものが得られる。これがＦＡの基本的な原理である。なお、ＦＡのアルゴリズムとしては、最尤法以外のアルゴリズムを利用した種々のアルゴリズムが存在する。本実施形態においても、このような各種のＦＡを利用可能である。

ところで、ＦＡで得られる推定値は、あくまでＡＡＴの値であり、この値に適合する混合行列Ａを決めた場合に、ランダムノイズの影響を低減しつつデータを無相関化することが可能だが、回転の自由度が残るために複数の構成成分ｓｉを一意に決めることができない。一方、ＩＣＡは、複数の構成成分ｓｉが互いに直交するように複数の構成成分ｓｉの回転の自由度を減らす処理である。そこで、本実施形態では、ＦＡで求める混合行列Ａの値を白色化行列（白色化済の行列）として用い、残された回転に対する任意性をＩＣＡにより特定する。これにより、ランダムノイズにロバストな白色化処理を行った後に、ＩＣＡを実行することにより、互いに直交する独立な構成成分ｓｉを決定することが可能である。また、このような処理の結果、ランダムノイズの影響を低減して、構成成分ｓｉに関する検量精度を向上させることができる。

Ｅ−３．ＩＣＡ（独立性指標としての尖度及びβダイバージェンス）：
ＩＣＡ（独立成分分析）では、一般に、独立成分の分離のための指標として、分離したデータ同士の独立性を表す高次統計量を独立性指標として用いられる。尖度は、典型的な独立性指標である。独立性指標として尖度を用いたＩＣＡについては、例えばAapo Hyvarinen, Juha Karhumen, Erkki Oja, "Independent Comonent Analysis", 2001, John Wiley &Sons, Inc.（「独立成分分析」、２００５年２月、東京電気大学出版部発行）の第８章に詳述されている。

しかし、処理対象データにスパイクノイズのような外れ値が入っている場合には、その外れ値も含めた統計量が独立性指標として計算されてしまう。このため、処理対象データについての本来の統計量と、算出された統計量との間に誤差が生じ、分離精度の低下を引き起こす場合がある。そこで、処理対象データ中の外れ値からの影響を受けにくい独立性指標を使用することが好ましい。この様な特性を持つ独立性指標として、βダイバージェンスを使用可能である。以下では、ＩＣＡにおける独立性指標としてのβダイバージェンスの原理を説明する。

前述したように、一般にＩＣＡでは、処理対象データｘを、構成成分ｓｉの線形和として表す線形混合モデル（上記式（１９））を仮定し、混合比ｃｉと構成成分ｓｉとを求める。ＩＣＡにより求められる構成成分ｓの推定値ｙは、分離行列Ｗを用いてｙ＝Ｗ・ｙと表わされる。このとき分離行列Ｗは混合行列Ａの逆行列であることが望まれる。

ここで、分離行列Ｗの推定値＾Ｗの対数尤度関数Ｌ（＾Ｗ）は、以下の式で表すことができる。

ここで、積算記号Σの要素は、各データ点ｘ（ｔ）における対数尤度である。この対数尤度関数Ｌ（＾Ｗ）を、ＩＣＡにおける独立性指標として用いることができる。βダイバージェンスの手法は、この対数尤度関数Ｌ（＾Ｗ）に適当な関数を作用させることにより、データ中のスパイクノイズのような外れ値に対して、その影響を抑えるように対数尤度関数Ｌ（＾Ｗ）を変換しようとする方法である。

βダイバージェンスを独立性指標として利用する場合には、まず、予め選択された関数Φβを用いて対数尤度関数Ｌ（＾Ｗ）を以下の式で変換する。

そして、この関数ＬΦ（＾Ｗ）を新たな尤度関数と考える。

スパイクノイズのような外れ値の影響を小さくするための関数Φβとしては、対数尤度の値（関数Φβの括弧内の値）が小さくなるほど、関数Φβの値が指数関数的に減衰するような関数が考えられる。このような関数Φβとしては、例えば以下を使用することができる。

この関数では、βの値が大きいほど、各データ点ｚ（上述の式（２７）では対数尤度）に対する関数値が小さくなる。このβの値は経験的に決定することができ、例えば約０．１に設定することが可能である。なお、この関数Φβとしては、式（２８）のものに限らず、βの値が大きいほど各データ点ｚに対する関数値が小さくなるような他の関数を利用することも可能である。

このようなβダイバージェンスを独立性指標として使用すると、スパイクノイズのような外れ値の影響を適切に抑えることができる。上記式（２７）のような尤度関数ＬΦ（＾Ｗ）を考える場合に、この尤度の最大化に対応して最小化される確率分布間の擬似距離がβダイバージェンスである。このようなβダイバージェンスを独立性指標として使用したＩＣＡを実行すれば、スパイクノイズのような外れ値の影響を低減して、構成成分ｓｉに関する検量精度を向上させることができる。

なお、βダイバージェンスを用いたＩＣＡについては、例えば、Minami Mihoko, Shinto Eguchi, “Robust Blind Source Separation by β-Divergence”, 2002に説明されている。

Ｆ．変形例：
この発明は前記実施例やその変形例に限られるものではなく、その要旨を逸脱しない範囲において種々の態様において実施することが可能であり、例えば次のような変形も可能である。

・変形例１：
前記実施形態では、未知成分のスペクトルＳの要素数ｍは予め経験的又は実験的に決められるとしたが、未知成分のスペクトルＳの要素数ｍはＭＤＬ（Minimum Description Length）やＡＩＣ（Akaike Information Criteria）として知られる情報量基準などによって決定してもよい。ＭＤＬなどを用いる場合、未知成分のスペクトルＳの要素数ｍはサンプルの観測データから演算によって自動的に決めることができる。なお、ＭＤＬについては、例えば”Independent component analysis for noisy data ? MEG data analysis, 2000”に説明されている。

・変形例２：
前記実施形態では、検量処理の対象となる被検体は、検量線を作成したときに用いたサンプルと同じ成分で構成されるとしたが、被検体に検量線を作成したときに用いたサンプルと同じ成分以外の未知成分が含まれる、としてもよい。独立成分同士の内積は0と仮定するため、未知成分に対応する独立成分とも内積0と考えられるためであり、内積で混合係数を求める場合には未知成分の影響は無視できる。

・変形例３：
前記実施形態において用いたコンピューターは、専用の装置として構成してもよい。例えば、ハードウェア回路のみで図７や図１３に示した装置を実現してもよい。或いは、図７や図１３に示した装置の機能の一部をハードウェア回路で実現し、他の一部をソフトウェアで実現してもよい。

・変形例４：
前記実施形態では、サンプルや被検体についての分光反射率のスペクトルの入力を、分光計測器により測定されたスペクトルを入力することで行っていたが、本発明はこれに限られない。例えば、波長帯域の相違する複数のバンド画像から分光スペクトルを推定し、この分光スペクトルを入力する構成としてもよい。前記バンド画像は、例えば、透過波長帯域を変更可能なフィルターを備えるマルチバンドカメラによってサンプルや被検体を撮影することで得られる。

・変形例５：
上記実施形態では、独立成分分析を用いて独立成分を選択する際に中性脂肪に適した波長のピークを有するものを選択していたが、この代わりに、主成分分析やＰＬＳ回帰分析等の手法を用いて主成分を選択する際に、中性脂肪に適した波長のピークを有するものを選択するようにしてもよい。この場合にも、独立成分の代わりに主成分を用いて中性脂肪濃度の検量線作成処理や検量処理を実行することによって検量精度を高めることが可能である。

なお、前述した各実施例および各変形例における構成要素の中の、独立請求項で記載された要素以外の要素は、付加的な要素であり、適宜省略可能である。

１０…ＣＰＵ、２０…メモリー、３０…ハードディスクドライブ、５０…入力インターフェイス、６０…出力インターフェイス、１００…パーソナルコンピューター、２００…分光計測器、４００…検量線作成装置、４１０…サンプル観測データ取得部、４２０…サンプル目的成分量取得部（中性脂肪濃度取得部）、４３０…混合係数推定部（推定部）、４３２…独立成分行列算出部、４３４…推定混合行列算出部、４３６…混合係数選択部、４３８…独立成分選択部、４４０…回帰式算出部、４５０…第１前処理部（正規化処理部）、４５２…標準正規変量変換（ＳＮＶ）、４５４…零空間射影法（ＰＮＳ）、４６０…第２前処理部（白色化処理部）、４６２…主成分分析（ＰＣＮ）、４６４…因子分析（ＦＡ）、４７０…独立成分分析処理部、４７２…第１処理、４７４…第２処理、５００…検量装置、５１０…被検体観測データ取得部、５２０…検量用データ取得部、５３０…混合係数算出部（内積演算部）、５３２…前処理部、５４０…目的成分量算出部、５５０…不揮発性記憶装置

Claims (7)

1. 目的成分としての中性脂肪を含む被検体についての中性脂肪濃度を求める中性脂肪濃度検量装置であって、
   前記被検体についての観測分光データを取得する被検体観測データ取得部と、
   中性脂肪に対応する独立成分又は主成分と、検量用の単回帰式と、含む検量用データを取得する検量用データ取得部と、
   前記被検体についての観測分光データと前記検量用データとに基づいて、前記被検体についての中性脂肪に対する混合係数を求める混合係数算出部と、
   前記中性脂肪に対応する混合係数と中性脂肪濃度との関係を示す前記単回帰式と、前記混合係数算出部によって求められた混合係数とに基づいて、前記中性脂肪濃度を算出する目的成分量算出部と、
   を含み、
   前記混合係数算出部は、前記中性脂肪に対する独立成分又は主成分として、予め選択された波長にピークを有する成分を使用する、中性脂肪濃度検量装置。
2. 請求項１に記載の中性脂肪濃度検量装置であって、
   前記予め選択された波長は、
   ９３０±３０ｎｍ又は９４５±３０ｎｍの一方と、１０６０±３０ｎｍとの少なくとも１つの波長である、中性脂肪濃度検量装置。
3. 請求項１に記載の中性脂肪濃度検量装置であって、
   前記予め選択された波長は、
   １１５０±３０ｎｍ又は１１８２±３０ｎｍの波長である、中性脂肪濃度検量装置。
4. 請求項１から請求項３のいずれか１項に記載の中性脂肪濃度検量装置を備える電子機器。
5. 中性脂肪を含む被検体の観測分光データから、前記被検体の中性脂肪濃度を導くことに用いる検量線を作成する検量線作成方法であって、
   前記被検体の複数のサンプルについての前記観測分光データを取得するサンプル観測データ取得工程と、
   前記各サンプルについての中性脂肪濃度を取得する中性脂肪濃度取得工程と、
   前記サンプル毎の観測分光データを複数の独立成分又は主成分に分離したときの複数の独立成分又は主成分を推定する推定工程と、
   前記複数のサンプルの中性脂肪濃度と、前記サンプル毎の前記観測分光データにおける前記独立成分又は主成分の混合係数とに基づいて、前記検量線の回帰式を求める回帰式算出工程と、
   を含み、
   前記推定工程は、前記独立成分又は主成分として、予め選択された波長にピークを有する成分波形信号を選択する工程である、検量線作成方法。
6. 請求項５に記載の検量線作成方法であって、
   前記予め選択された波長は、
   ９３０±３０ｎｍ又は９４５±３０ｎｍの一方と、１０６０±３０ｎｍとの少なくとも１つの波長である、検量線作成方法。
7. 請求項５に記載の検量線作成方法であって、
   前記予め選択された波長は、
   １１５０±３０ｎｍ又は１１８２±３０ｎｍの波長である、検量線作成方法。

Images