"一定の波長間隔で反射率を算出して反射スペクトルを求める工程と、前記反射スペクトルを高次関数に近似する工程と、前記高次関数から最小の反射率を含む波長区間を選択する工程と、前記波長区間で前記高次関数をそれより低次の2次関数に近似する工程と、前記2次関数を波長で微分してその値が0となる解を求める工程と、を備えることを特徴とする分子間相互作用の検出方法。一定の波長間隔で反射率を算出して反射スペクトルを求める工程と、前記反射スペクトルをフィルタリングする工程と、フィルタリング後の前記反射スペクトルから最小の反射率を含みかつ変曲点を境界とする第1の波長区間を選択する工程と、前記第1の波長区間でフィルタリング後の前記反射スペクトルを高次関数に近似する工程と、前記高次関数から最小の反射率を含む第2の波長区間を選択する工程と、前記第2の波長区間で前記高次関数をそれより低次の2次関数に近似する工程と、前記2次関数を波長で微分してその値が0となる解を求める工程と、を備えることを特徴とする分子間相互作用の検出方法。リガンドを有する検出器と、前記検出器に対し白色光を照射する白色光源と、前記検出器からの反射光を検出する分光器と、前記白色光源と前記分光器とを制御する制御装置であって、一定の波長間隔で反射率を算出して反射スペクトルを求め、前記反射スペクトルを高次関数に近似し、前記高次関数から最小の反射率を含む波長区間を選択し、前記波長区間で前記高次関数をそれより低次の2次関数に近似し、前記2次関数を波長で微分してその値が0となる解を求める前記制御装置と、を備えることを特徴とする分子間相互作用の検出装置。リガンドを有する検出器と、前記検出器に対し白色光を照射する白色光源と、前記検出器からの反射光を検出する分光器と、前記白色光源と前記分光器とを制御する制御装置であって、一定の波長間隔で反射率を算出して反射スペクトルを求め、前記反射スペクトルをフィルタリングし、フィルタリング後の前記反射スペクトルから最小の反射率を含みかつ変曲点を境界とする第1の波長区間を選択し、前記第1の波長区間でフィルタリング後の前記反射スペクトルを高次関数に近似し、前記高次関数から最小の反射率を含む第2の波長区間を選択し、前記第2の波長区間で前記高次関数をそれより低次の2次関数に近似し、前記2次関数を波長で微分してその値が0となる解を求める前記制御装置と、を備えることを特徴とする分子間相互作用の検出装置。ところで、ボトムピーク波長の変化の推移を観測する上記RIfS方式においては、実際に検出した反射光の波長間に真の反射率ボトムピーク波長が存在することが多いため、実際に検出した反射光の波長間隔よりも小さい狭い波長間隔でボトムピーク波長の推移を追跡する必要がある。これに対しては、通常、反射率の波長分布を高次多項式により近似し、その高次多項式からその解最小値を求めてこれをボトムピーク波長とすることが行われている。しかしながら、高次多項式による近似はフィッティングという観点では有用であるものの、高次ゆえにその多項式から解を求めるには時間がかかり、ボトムピーク波長を容易かつ短時間で算出特定することが困難である。また、ボトムピーク波長の変化の推移を観測する上記RIfS方式においては、複数種類の分子を取り扱うことも可能となっており、複数種類の分子を積み重ねたり、サンプル溶液中に複数種類のアナライト130を混入させたりしても、分子間相互作用を検出することができる。例えば、図12に示すとおり、複雑な抗原抗体反応といったタンパク質の分子間相互作用を検出することができる。図12の例によれば、光学薄膜104の表面をアミノ化しa参照、NHSPEG4を介してビオチン150を結合させb参照、アビジン152を結合させc参照、BSA154でブロッキングして抗体156を結合させd参照、抗原158を検出するe参照。ボトムピーク波長の変化の推移を経時的に観測すると、図13に示すとおり、区間160でビオチン150,アビジン152の結合を、区間162においてBSA154のブロッキングを、区間164において抗体156の結合を、区間166で濃度に応じた抗原158の結合を、それぞれ検出することができる。実際に複数種類の分子の相互作用を検出しようとする場合において、ボトムピーク波長を算出特定するときには、例えば、図14a,bに示すような反射スペクトル170を得てこれを高次関数で近似し、ボトムピーク波長を算出特定する。しかしながら、この場合に、複数種類の分子を扱うがゆえに、反射スペクトル170には極小点172が存在し、反射スペクトル170を近似すると、極小点172の影響を受けた反射スペクトル174が得られ、それに伴いボトムピーク波長がずれる。そのため、反射スペクトル170における正確なボトムピーク波長を算出特定できないことがある。したがって、本発明の主な目的は、反射スペクトルの極小点の有無に関らず、ボトムピーク波長を容易かつ短時間で算出特定することができる分子間相互作用の検出方法およびその検出装置を提供することにある。本発明の一態様によれば、一定の波長間隔で反射率を算出して反射スペクトルを求める工程と、前記反射スペクトルを高次関数に近似する工程と、前記高次関数から最小の反射率を含む波長区間を選択する工程と、前記波長区間で前記高次関数をそれより低次の2次関数に近似する工程と、前記2次関数を波長で微分してその値が0となる解を求める工程と、を備えることを特徴とする分子間相互作用の検出方法が提供される。本発明の他の態様によれば、一定の波長間隔で反射率を算出して反射スペクトルを求める工程と、前記反射スペクトルをフィルタリングする工程と、フィルタリング後の前記反射スペクトルから最小の反射率を含みかつ変曲点を境界とする第1の波長区間を選択する工程と、前記第1の波長区間でフィルタリング後の前記反射スペクトルを高次関数に近似する工程と、前記高次関数から最小の反射率を含む第2の波長区間を選択する工程と、前記第2の波長区間で前記高次関数をそれより低次の2次関数に近似する工程と、前記2次関数を波長で微分してその値が0となる解を求める工程と、を備えることを特徴とする分子間相互作用の検出方法が提供される。本発明の他の態様によれば、リガンドを有する検出器と、前記検出器に対し白色光を照射する白色光源と、前記検出器からの反射光を検出する分光器と、前記白色光源と前記分光器とを制御する制御装置であって、一定の波長間隔で反射率を算出して反射スペクトルを求め、前記反射スペクトルを高次関数に近似し、前記高次関数から最小の反射率を含む波長区間を選択し、前記波長区間で前記高次関数をそれより低次の2次関数に近似し、前記2次関数を波長で微分してその値が0となる解を求める前記制御装置と、を備えることを特徴とする分子間相互作用の検出装置が提供される。本発明の他の態様によれば、リガンドを有する検出器と、前記検出器に対し白色光を照射する白色光源と、前記検出器からの反射光を検出する分光器と、前記白色光源と前記分光器とを制御する制御装置であって、一定の波長間隔で反射率を算出して反射スペクトルを求め、前記反射スペクトルをフィルタリングし、フィルタリング後の前記反射スペクトルから最小の反射率を含みかつ変曲点を境界とする第1の波長区間を選択し、前記第1の波長区間でフィルタリング後の前記反射スペクトルを高次関数に近似し、前記高次関数から最小の反射率を含む第2の波長区間を選択し、前記第2の波長区間で前記高次関数をそれより低次の2次関数に近似し、前記2次関数を波長で微分してその値が0となる解を求める前記制御装置と、を備えることを特徴とする分子間相互作用の検出装置が提供される。本発明によれば、反射率の波長分布反射スペクトルを高次関数で近似した後に、波長区間を選択してその波長区間で反射スペクトルを2次関数で再度近似するから、その2次関数を微分してその値が0となる解を求めるといった簡単な処理でボトムピーク波長が求められ、反射スペクトルの極小点の有無に関らず、ボトムピーク波長を容易かつ短時間で算出特定することができる。第1の実施形態以下、図面を参照しながら本発明の好ましい第1の実施形態について説明する。図1に示すとおり、検出装置1は主に、検出器10,白色光源20,分光器30,光ファイバ40,制御装置50などから構成されている。検出器10は基本的にはセンサーチップ12,フローセル14から構成されている。図2に示すとおり、センサーチップ12は矩形状を呈したシリコン基板12aを有している。シリコン基板12a上にはSiN膜12b窒化シリコンが蒸着されている。SiN膜12bは光学薄膜の一例である。フローセル14はシリコーンゴム製の透明な部材である。フローセル14には溝14aが形成されている。フローセル14をセンサーチップ12に密着させると、密閉流路14bが形成される図1参照。溝14aの両端部はフローセル14の表面から露出しており、一方の端部がサンプル溶液の流入口14cとして、他方の端部がその流出口14dとしてそれぞれ機能するようになっている。溝14の底部にはリガンド16が結合されている図1参照。検出器10では、センサーチップ12に対しフローセル14を貼り替え可能となっており、フローセル14はディスポーザブル使い捨て使用が可能となっている。センサーチップ12の表面には、シランカップリング剤などにより、表面修飾をおこなってもよく、この場合フローセル14の貼り替えが容易となる。図1に示すとおり、フローセル14の密閉流路14bの上方には光ファイバ40が設置されている。光ファイバ40の一方の端部には白色光源20が接続されている。白色光源20としては例えばハロゲン光源が使用される。光ファイバ40の他方の端部には分光器30が接続されている。白色光源20が点灯すると、その光が光ファイバ40を介して密閉流路14bに照射され、その反射光が再度光ファイバ40を介して分光器30で検出される。白色光源20や分光器30は制御装置50に接続され、制御装置50はこれら部材の動作を制御するようになっている。特に検出装置1では、制御装置50により、分光器30による検出結果からその分光特性や反射特性などを解析可能となっている。続いて、検出装置1を用いた分子間相互作用の検出方法について説明する。図1に示すとおり、アナライト62を含むサンプル溶液60を、流入口14cから密閉流路14bを経て流出口14dに流通させる。アナライト62とは、リガンド16と特異的に結合する物質であり、検出しようとする目的の分子である。アナライト62としては、例えばタンパク質,核酸,脂質,糖などの生体分子や、薬剤物質,内分泌錯乱化学物質などの生体分子と結合する外来物質などが使用される。サンプル溶液60が密閉流路14bを流通している間に、白色光源20を点灯させる。白色光はフローセル14を透過してセンサーチップ12に照射され、その反射光が分光器30で検出される。分光器30により検出された反射光の検出強度は制御装置50に送信される。この場合に、図3に示すとおり、サンプル溶液60中のアナライト62がリガンド16と結合すると、光学的厚さ70が変化し、干渉色分光器30による検出強度が最も小さくなる波長が変化する。制御装置50は、分光器30による検出結果を受けてアナライト62とリガンド16との結合前後での反射スペクトルの各ボトムピーク波長λbottomを算出特定する。ボトムピーク波長λbottomの算出特定方法は、主には、図4のS1~S4の工程から構成され、S1~S4の工程の処理は、制御装置50において実行される。S1反射スペクトルを求め高次関数に近似する。S2高次関数から任意の波長区間を選択する。S3選択した波長区間で高次関数をそれより低次の2次関数に近似する。S42次関数からボトムピーク波長λbottomを求める。なお、ボトムピーク波長λbottomの算出特定方法はアナライト62とリガンド16との結合前後で同様である。S1の工程では、分光器30の検出結果から、波長に対する反射率を一定の波長間隔例えば0.35nm間隔で算出し、図5aに示すとおり、測定波長範囲例えば400~800nmにおける波長と反射率との関係反射スペクトル72を求める。この段階では、反射スペクトル72の波形は微小な凹凸が繰り返されるような不規則な形状を呈しており、ボトムピーク波長λbottomを算出特定するのが困難な状態となっている。その後、反射スペクトル72を20次程度の高次関数で近似し、図5bに示すとおり、反射スペクトル72の波形を滑らかにする。当該近似は公知のいずれの手法によるものであってもよく、例えば線形の関数と擬似voigt関数の和で実現することができる。S2の工程では、図5bに示すとおり、近似後の反射スペクトル72高次関数において、反射率の最小値をとる波長を基準波長λ0として任意に波長区間74を選択決定する。すなわち、S1の工程では、高次関数においても測定波長に対応する反射率が算出されているから、それら反射率のなかから、最小値をとる波長を特定してこれを基準波長λ0とする図5c参照。その後、基準波長λ0を中心としてその前後で数点の波長を選択し、その選択した区間を波長区間74とする。波長区間74の選択は任意であり、少なくとも基準波長λ0を含んだ3点を選択すればよい。S3の工程では、波長区間74において近似後の反射スペクトル72高次関数を2次関数で近似する。当該近似は公知のいずれの手法によるものであってもよく、例えば最小二乗法で実現することができる。この場合、波長区間74における波長の反射率は、波長をx,反射率をyとして、式1で表される。y=ax2bxc...1式1中、a,b,cは任意の定数である。S4の工程では、式1の2次関数から、ボトムピーク波長λbottomを求める。式1の2次関数は、下に凸の放物線を描くものであり、ある波長における傾きが0ゼロであれば、その波長が最小の反射率をとりうる波長であってボトムピーク波長λbottomに相当する。したがって、S4の工程では、式1の2次関数を波長xで微分し、その微分後の値を0とする方程式を導き出し、その解を求める。微分後の値を0とする方程式は式2で表され、その解は式3で表される。2axb=0...2x=-b/2a...3このような処理から、ボトムピーク波長λbottomを-b/2aと算出特定することができる。以上の第1の実施形態によれば、S1の工程において分光器30の検出結果から反射スペクトル72を高次関数で近似した後に、S2,S3の工程において波長区間74を決定し反射スペクトル72を2次関数で再度近似するから、S4の工程においてはその2次関数を微分してその値が0となる解を求めるといった簡単な処理を実行することにより、ボトムピーク波長λbottomが算出特定され、ボトムピーク波長λbottomを容易かつ短時間で算出特定することができる。第2の実施形態以下、図面を参照しながら本発明の好ましい第2の実施形態について説明する。第2の実施形態における検出装置の構成は第1の実施形態の検出装置1と同じである。第2の実施形態が第1の実施形態と異なるのは、ボトムピーク波長λbottomの算出特定方法である。以下、第2の実施形態におけるボトムピーク波長λbottomの算出特定方法について説明する。ボトムピーク波長λbottomの算出特定方法は、主には、図7のS11~S16の工程から構成され、S11~S16の工程の処理は、制御装置50において実行される。S11反射スペクトルを求めフィルタリングする。S12フィルタリング後の反射スペクトルから任意の波長区間を選択する。S13選択した波長区間でフィルタリング後の反射スペクトルを高次関数に近似する。S14高次関数から任意の波長区間を選択する。S15選択した波長区間で高次関数をそれより低次の2次関数に近似する。S162次関数からボトムピーク波長λbottomを求める。なお、ボトムピーク波長λbottomの算出特定方法はアナライト62とリガンド16との結合前後で同様である。S11の工程では、分光器30の検出結果から、波長に対する反射率を一定の波長間隔例えば0.35nm間隔で算出し、図8aに示すとおり、測定波長範囲例えば400~800nmにおける波長と反射率との関係反射スペクトル72を求める。この段階では、反射スペクトル72の波形は微小な凹凸が繰り返されるような不規則な形状を呈しており、ボトムピーク波長λbottomを算出特定するのが困難な状態となっている。その後、反射スペクトル72をフィルタリング平滑化処理し、図8bに示すとおり、反射スペクトル72の波形を滑らかにする。当該フィルタリングでは、一定の波長間隔で各波長における反射率を平均化する。S12の工程では、図8bに示すとおり、フィルタリング後の反射スペクトル72において、反射率の最小値をとる波長を基準波長λ0として任意に波長区間76を選択決定する。すなわち、S11の工程では、一定の波長間隔で平均化された反射率が算出されているから、それら反射率のなかから、最小値をとる波長を特定してこれを基準波長λ0とする。その後、基準波長λ0を中心としてその前後で波長を選択し、その選択した区間を波長区間76とする。基準波長λ0の短波長側や長波長側に変曲点78が存在する場合には、変曲点78を波長区間76の境界とする。変曲点78より数点2~3点分の波長だけ基準波長λ0側に寄せた地点を、波長区間76の境界としてもよい。S13の工程では、波長区間76においてフィルタリング後の反射スペクトル72を20次程度の高次関数で近似し、図8cに示すとおり、反射スペクトル72の波形をさらに滑らかにする。当該近似は公知のいずれの手法によるものであってもよく、例えば線形の関数と擬似voigt関数の和で実現することができる。S14の工程では、図8c,図8dに示すとおり、近似後の反射スペクトル72高次関数において、反射率の最小値をとる波長を基準波長λ1として任意に波長区間74を選択決定する。すなわち、S13の工程では、高次関数においても測定波長に対応する反射率が算出されているから、それら反射率のなかから、最小値をとる波長を特定してこれを基準波長λ1とする図8d参照。その後、基準波長λ1を中心としてその前後で数点の波長を選択し、その選択した区間を波長区間74とする。波長区間74の選択は任意であり、少なくとも基準波長λ1を含んだ3点を選択すればよい。S15の工程では、波長区間74において近似後の反射スペクトル72高次関数を2次関数で近似する。当該近似は公知のいずれの手法によるものであってもよく、例えば最小二乗法で実現することができる。この場合、波長区間74における波長の反射率は、波長をx,反射率をyとして、前述の式1で表される。S16の工程では、式1の2次関数から、ボトムピーク波長λbottomを求める。式1の2次関数は、下に凸の放物線を描くものであり、ある波長における傾きが0ゼロであれば、その波長が最小の反射率をとりうる波長であってボトムピーク波長λbottomに相当する。したがって、S16の工程では、式1の2次関数を波長xで微分し、その微分後の値を0とする方程式を導き出し、その解を求める。微分後の値を0とする方程式は前述の式2で表され、その解は前述の式3で表される。このような処理から、ボトムピーク波長λbottomを-b/2aと算出特定することができる。以上の第2の実施形態によれば、S11の工程において分光器30の検出結果から反射スペクトル72を求めてフィルタリングした後に、S12の工程において変曲点78を境界とする波長区間76を選択決定するから、反射スペクトル72に極小点が複数存在する場合であっても、その極小点を回避した範囲除く範囲でボトムピーク波長λbottomを求められ、ボトムピーク波長λbottomを正確に算出特定することができる。さらに、第2の実施形態によれば、S13の工程においてフィルタリング後の反射スペクトル72を高次関数で近似した後に、S14,S15の工程において波長区間74を決定し反射スペクトル72を2次関数で再度近似するから、S16の工程においてはその2次関数を微分してその値が0となる解を求めるといった簡単な処理を実行することにより、ボトムピーク波長λbottomが算出特定され、ボトムピーク波長λbottomを容易かつ短時間で算出特定することができる。"