(57)【要約】【課題】    ニューラルネットワークを利用して被検試料中の複数の化学物質の濃度を測定する公知の方法よりもさらに正確に該化学物質の濃度を測定することができる、被検試料中の複数の化学物質の濃度の測定方法を提供すること。【解決手段】  測定された実データの系の確率分布を推定し、推定分布に従って生成した多数の擬似データ生成し、これをニューラルネットワークに学習させることにより、限られた実測データ数でも高精度に被検試料中の複数の化学物質の濃度を測定することが可能である。【選択図】  なし【特許請求の範囲】
【請求項１】
各既知濃度の複数の化学物質を含む標準試料に、該化学物質の濃度を測定する１又は複数のセンサを適用し、センサ応答を測定することにより、複数の該化学物質の濃度とそれに対応するセンサ応答の組みからなる化学データを得る第１工程と、
得られた化学データを滑らかに補間する曲面を作成する第２工程と、
第２工程における化学データの補間誤差を確率変数として捉え、補間誤差の確率密度関数を推定する第３工程と、
前記第３工程で得られる補間誤差の推定確率密度関数に従う擬似乱数ベクトルを生成し、前記第２工程において得られた補間曲面上のベクトルからランダムに選び出したベクトル擬似補間ベクトルに、擬似補間誤差ベクトルを足し合わせて新たなデータベクトルを生成することを多数繰り返すことで、補間曲面・補間誤差の特徴を反映した多数のデータベクトルである擬似データを生成する第４工程と、
前記第４工程で生成した前記擬似データをニューラルネットワークに学習させる第５工程と、
未知の被検試料について前記センサを適用し、センサ応答を測定する第６工程と、
前記第５工程で学習済みとなったニューラルネットワークに、第６工程で測定されたセンサ応答を入力し、前記ニューラルネットワークの出力より複数の該化学物質の未知濃度を推定する第７工程とを含む、被検試料中の複数の化学物質濃度の測定方法。
【請求項２】
前記複数の化学物質は、目的とする１つの標的化学物質と、その測定にとっての妨害物質であり、前記センサは、該標的化学物質を測定するセンサである請求項１記載の方法。
【請求項３】
前記センサ応答の測定は、前記標的物質を基質とする酵素を含む酵素電極を作用電極として用いるボルタンメトリーにおいて複数の電位における電流を測定することにより行われる請求項２記載の方法。
【請求項４】
前記酵素がオキシダーゼである請求項３記載の方法。
【請求項５】
前記オキシダーゼがグルコースオキシダーゼである請求項４記載の方法。
【請求項６】
前記妨害物質が、アスコルビン酸、アセトアミノフェン及び尿酸から成る群より選ばれる少なくとも１種である請求項５記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被検試料中に含まれる複数の化学物質の測定方法に関する。
【背景技術】
【０００２】
被検試料中の化学物質の濃度を測定するためには、目的とする化学物質の濃度を測定できるセンサが種々用いられている。しかし、一般に、センサの特異性は必ずしもそれほど高くはなく、被検試料中に複数の化学物質が含まれる場合には、それらが互いに影響を及ぼすため、単にセンサ応答を測定するだけでは正確な化学物質濃度を測定することができない。
【０００３】
本願発明者らは、先にニューラルネットワークを利用して、複数の化学物質の濃度の測定方法を発明し、特許出願した（特許文献１）。また、本願発明者らは、先に妨害物質が共存する被検試料中の被検物質の濃度を簡便、かつ正確に測定することができる、ニューラルネットワークを用いた被検物質の定量方法を提供した（非特許文献１）。
【０００４】
【非特許文献１】日本分析化学会第５３年会講演要旨集
【非特許文献２】文献「ニューロコンピューティング入門」(森北出版 1997年、坂和正敏・田中雅博 共著)
【非特許文献３】Using Additive Noise in Backpropagation Training   L.Holmstrom, P.Koistinen,  IEEE Trans. Neural Networks, Vol.3, No.1, pp.24-38, Jan (1992))
【非特許文献４】「乱数とモンテカルロ法」(森北出版、1978年、宮武修、脇本和昌  共著)
【特許文献１】国際公開公報WO 03/044498 A1
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上記のとおり、ニューラルネットワークを利用して被検試料中の複数の化学物質の濃度を測定する方法は、既に知られている。しかし、これらの公知の方法よりもさらに正確に化学物質濃度を測定することができれば望ましいことは言うまでもない。センサ応答に基づき、ニューラルネットワークを駆使して化学物質濃度を推定する場合、データ数が多いほど正確な推測が可能になるが、化学データはいちいち計測する必要があるため、データの数には限りがある。限られたデータ数でもより正確に化学物質濃度を測定することが可能な方法が望まれる。
【０００６】
したがって、本発明の目的は、ニューラルネットワークを利用して被検試料中の複数の化学物質の濃度を測定する公知の方法よりもさらに正確に該化学物質の濃度を測定することができる、被検試料中の複数の化学物質の濃度の測定方法を提供することである。
【課題を解決するための手段】
【０００７】
本願発明者らは、鋭意研究の結果、測定された実データの系の確率分布を推定し、推定分布に従って生成した多数の擬似データ生成し、これをニューラルネットワークに学習させることにより、限られた実測データ数でも高精度に被検試料中の複数の化学物質の濃度を測定することが可能であることを見出し、本発明を完成した。
【０００８】
すなわち、本発明は、各既知濃度の複数の化学物質を含む標準試料に、該化学物質の濃度を測定する１又は複数のセンサを適用し、センサ応答を測定することにより、複数の該化学物質の濃度とそれに対応するセンサ応答の組みからなる化学データを得る第１工程と、
得られた化学データを滑らかに補間する曲面を作成する第２工程と、
第２工程における化学データの補間誤差を確率変数として捉え、補間誤差の確率密度関数を推定する第３工程と、
前記第３工程で得られる補間誤差の推定確率密度関数に従う擬似乱数ベクトルを生成し、前記第２工程において得られた補間曲面上のベクトルからランダムに選び出したベクトル擬似補間ベクトルに、擬似補間誤差ベクトルを足し合わせて新たなデータベクトルを生成することを多数繰り返すことで、補間曲面・補間誤差の特徴を反映した多数のデータベクトルである擬似データを生成する第４工程と、
前記第４工程で生成した前記擬似データをニューラルネットワークに学習させる第５工程と、
未知の被検試料について前記センサを適用し、センサ応答を測定する第６工程と、
前記第５工程で学習済みとなったニューラルネットワークに、第６工程で測定されたセンサ応答を入力し、前記ニューラルネットワークの出力より複数の該化学物質の未知濃度を推定する第７工程とを含む、被検試料中の複数の化学物質濃度の測定方法を提供する。
【発明の効果】
【０００９】
本発明により、限られた実測データ数でも高精度に被検試料中の複数の化学物質の濃度を測定することが可能な複数の化学物質濃度の新規な測定方法が提供された。本発明の方法によれば、測定された実データの系の確率分布を推定し、推定分布に従って生成した多数の擬似データ生成し、これをニューラルネットワークに学習させるので、限られたデータ数でもより正確に化学物質濃度を測定することが可能である。
【発明を実施するための最良の形態】
【００１０】
工程１
工程１では、各既知濃度の複数の化学物質を含む標準試料に、該化学物質の濃度を測定する１又は複数のセンサを適用し、センサ応答を測定することにより、複数の該化学物質の濃度とそれに対応するセンサ応答の組みからなる化学データを得る。この工程自体は公知の方法により行うことができる。
【００１１】
センサとしては、目的とする化学物質濃度を測定することができる、公知のいずれのセンサをも利用することができる。下記実施例では、その濃度の測定を目的とする１つの標的物質を基質とする酵素を含む酵素電極を作用電極として用いるボルタンメトリーにおいて複数の電位における電流を測定することにより化学データを得たが、必ずしもこの方法に限定されるものではなく、測定原理を異にする他の種々の公知のセンサを用いることができる。また、下記実施例では、目的とする標的物質の濃度を測定できる１つのセンサを用い、その測定を妨害する妨害物質の存在下で標的物質の濃度及び妨害物質の濃度を測定したが、本発明の方法は必ずしもこのような場合に限定されるものではなく、求めたい複数の化学物質の濃度をそれぞれ異なるセンサで測定することもできる。
【００１２】
工程２
本工程では、「工程1」で得られたデータの補間を行う。以下、データの補間に先立ち行われる前処理について、「工程2-1」、「工程2-2」に記す。また前処理後のデータの補間について「工程2-3」に記す。
【００１３】
[工程2-1 センサ応答に対する前処理]
2006343124.tif 000002
【数１】
2006343124.tif 000003
2006343124.tif 000004
下記実施例では120成分の応答を扱うので、となる。
各
【数２】
2006343124.tif 000005
に対して
【数３】
2006343124.tif 000006
のように変換を施す。
【数４】
2006343124.tif 000007
2006343124.tif 000008
【００１４】
以上の変換を前処理として行い、[0,1]の範囲に正規化されたセンサ応答ベクトル
【数５】
2006343124.tif 000009
を導出する。
【００１５】
[工程2-2 化学物質の濃度に対する前処理]
2006343124.tif 000010
【数６】
2006343124.tif 000011
【００１６】
下記実施例では、4種類の化学物質の濃度推定を行うので、J=4となる。前処理として、
【数７】
2006343124.tif 000012
に対して
【数８】
2006343124.tif 000013
として変換を行う。
【００１７】
この変換により、[0.1,0.9]の範囲に正規化されたベクトル
【数９】
2006343124.tif 000014
を導出する。
【００１８】
下記実施例では、データの構成は下記表１に示すとおりである。
【００１９】
【表１】
2006343124.tif 000015
【００２０】
[工程2-3 データの補間]
2006343124.tif 000016
【００２１】
【数１０】
2006343124.tif 000017
【００２２】
各データセットにおいて、データセットの各要素間を滑らかに内挿する。内挿法としは階層型ニューラルネットワークが望ましい。階層型ニューラルネットワークについては、非特許文献２に記載がある。階層型ニューラルネットワークは結合加重wを調整することで、非線形関数を近似する能力を持つ。
【００２３】
本工程では、各データセット
【数１１】
2006343124.tif 000018
ごとに階層型ニューラルネットワークを用意する。データセットiを学習し、データセットiの各要素間を補間する階層型ニューラルネットワークをANNiInterpolationとする。
【００２４】
ANNiInterpolationによる、データセット
【数１２】
2006343124.tif 000019
の学習では、各要素
【数１３】
2006343124.tif 000020
中の
【数１４】
2006343124.tif 000021
をANNiInterpolationの入力とし、
【数１５】
2006343124.tif 000022
2006343124.tif 000023
【数１６】
2006343124.tif 000024
からランダムに選ばれたデータ
【数１７】
2006343124.tif 000025
の学習誤差
【数１８】
2006343124.tif 000026
2006343124.tif 000027
【数１９】
2006343124.tif 000028
ηは学習係数と呼ばれるパラメータである。woldは前回の学習で得た結合加重を、wnewは今回の学習で更新した結合荷重を表す。学習を行った後、ANNiInterpolationの入出力関係によってデータセット
【数２０】
2006343124.tif 000029
の各要素間は滑らかに補間される。階層型ニューラルネットワークの構造および学習に関するパラメータは下記表２のように設定した。
【数２１】
2006343124.tif 000030
はANNiInterpolationに
【数２２】
2006343124.tif 000031
を入力した時の出力を表す。wiはANNiInterpolationの結合加重を表す。
【００２５】
【表２】
2006343124.tif 000032
【００２６】
以上の学習を各データセット
【数２３】
2006343124.tif 000033
について行い、各データセットの要素を補間する
【数２４】
2006343124.tif 000034
の結合加重
【数２５】
2006343124.tif 000035
を導出する。以後、学習により導出した結合加重
【数２６】
2006343124.tif 000036
を用いて、学習済みのANNiInterpolationの出力を
【数２７】
2006343124.tif 000037
として表す。
【００２７】
工程3
本工程では、「工程2」のデータ補間における補間誤差を確率変数として捉え、補間誤差の確率密度関数を推定する。以下、「工程3-1」において補間誤差の算出法について、「工程3-2」において補間誤差の前処理について、「工程3-3」では補間誤差の分布の推定について記す。
【００２８】
[工程 3-1 補間誤差の算出]
データセット
【数２８】
2006343124.tif 000038
に対するANNiInterpolationの補間誤差について、d番目のデータ
【数２９】
2006343124.tif 000039
に対する補間誤差をeidとして下記のように導出する。
【数３０】
2006343124.tif 000040
【００２９】
[工程3-2 補間誤差の前処理]
各補間誤差
【数３１】
2006343124.tif 000041
に対し、下記のようなスケール変換を前処理として行う。
【数３２】
2006343124.tif 000042
【００３０】
[工程3-3 補間誤差の分布推定]
各データセット
【数３３】
2006343124.tif 000043
2006343124.tif 000044
【数３４】
2006343124.tif 000045
を定義する。
【００３１】
補間誤差ベクトルを確率変数ベクトル
【数３５】
2006343124.tif 000046
として考える。
【数３６】
2006343124.tif 000047
2006343124.tif 000048
【００３２】
【数３７】
2006343124.tif 000049
ここでsは補間誤差ベクトルの標本共分散行列を表している。
ただしsが正則、つまりS-1が存在するために下の式のように
対角成分に小さな正数Δを加える。
【数３８】
2006343124.tif 000050
ここでEはIxIの単位行列を表す。
Δは0.01に設定した。
パラメータhは以下のように設定する。
【００３３】
2006343124.tif 000051
【数３９】
2006343124.tif 000052
に対して、
【数４０】
2006343124.tif 000053
を算出し、
【数４１】
2006343124.tif 000054
を最大化するhmaxをパラメータhとして採用する。
【００３４】
この設定法はCross-Validated-Methodと呼ばれ、非特許文献３に記載がある。hmaxの探索は以下のアルゴリズムで行う。
【００３５】
[工程3-3-2 hmaxの探索法]
次のようなアルゴリズムでhmaxを探索する。
【００３６】
(i)初期化
探索範囲の下限α及び上限β(α＜β)、探索幅
【数４２】
2006343124.tif 000055
を決定。
【数４３】
2006343124.tif 000056
とする。
【００３７】
(ii)繰り返し
(ii-ア) hに
【数４４】
2006343124.tif 000057
を加える。h>βなら終了。そうでなければ(ii-イ)へ進む。
【００３８】
(ii-イ)
【数４５】
2006343124.tif 000058
なら
【数４６】
2006343124.tif 000059
とする。(ii-ア)へ戻り、
さらに繰り返す。
【００３９】
パラメータは下記の表３のように設定した。
【００４０】
【表３】
2006343124.tif 000060
【００４１】
工程4
本工程では、「工程2」で得られた補間および「工程3」で推定した補間誤差の分布についての特徴を反映した擬似データを生成する。以下、「工程4-1」において、補間誤差ベクトルの推定確率密度関数P(ε)に従う擬似乱数である「擬似補間誤差ベクトル」を生成する手法を記す。「工程4-2」においては、「工程2」で補間された曲面上に存在するベクトルから、一様な擬似乱数に基づいて「擬似補間ベクトル」を生成する方法を、「工程4-3」では、擬似補間誤差ベクトルおよび擬似補間ベクトルを用いて擬似データを生成する手法を記す。
【００４２】
[工程4-1 擬似補間誤差ベクトルの生成]
補間誤差ベクトルの推定確率密度関数P(ε)に従う擬似的な乱数により擬似補間誤差ベクトルを生成する。P(ε)に従う擬似補間誤差ベクトル
【数４７】
2006343124.tif 000061
を生成する手法は以下の通りである。
【００４３】
まず、独立な(0,1)の一様分布から2I個の独立した乱数
【数４８】
2006343124.tif 000062
を生成し、
【数４９】
2006343124.tif 000063
について
【数５０】
2006343124.tif 000064
を算出する。そして、0からNまでの整数の中から等確率で1つの整数を選び出す。その整数Zにより
【数５１】
2006343124.tif 000065
として、擬似補間誤差ベクトル
【数５２】
2006343124.tif 000066
を算出する。
ここで
【数５３】
2006343124.tif 000067
DはS=DTDとなるような三角行列である。Dはcholesky分解により求める。cholesky分解については非特許文献４に記載がある。
さらにアフィン変換
【数５４】
2006343124.tif 000068
を前処理として行い、前処理済みの擬似補間誤差ベクトル
【数５５】
2006343124.tif 000069
を算出する。
【００４４】
以上の作業をある上限a1回繰り返すことで、P(ε)に従うa1個の擬似補間誤差ベクトルを前処理したベクトルの集合
【数５６】
2006343124.tif 000070
を得る。
【００４５】
[工程4-2 擬似補間ベクトルの生成]
データを補間する階層型ニューラルネットワーク
【数５７】
2006343124.tif 000071
を用いて、データ間を補間する曲面上において、ランダムにベクトル：擬似補間ベクトルを生成する。
【００４６】
擬似補間ベクトルは以下のように生成される。
まず、
【数５８】
2006343124.tif 000072
2006343124.tif 000073
【数５９】
2006343124.tif 000074
を生成し、ランダムな入力ベクトル
【数６０】
2006343124.tif 000075
とする。
【００４７】
次に、各
【数６１】
2006343124.tif 000076
に入力し、各ANNiInterpolationから出力
【数６２】
2006343124.tif 000077
を導出し、ベクトル
【数６３】
2006343124.tif 000078
を生成する。このようにして、データを補間する曲面上のベクトル
【数６４】
2006343124.tif 000079
を擬似補間ベクトルとして生成する。
【００４８】
以上の操作を上限回数a2まで行うことにより、擬似補間ベクトルの集合
【数６５】
2006343124.tif 000080
を導出する。
【００４９】
[工程4-3 擬似データの生成]
「工程4-1」、「工程4-2」で生成した擬似補間誤差ベクトルと擬似補間ベクトルにより擬似データを生成する。
【００５０】
擬似データの生成アルゴリズムは以下の通りである。
(i)初期化
b=1とする。
【００５１】
(ii)繰り返し
(ii-1) 1からa1までの整数をランダムに選択(すべての整数が選ばれる確率は等確率)。
選ばれた整数をA1とする。
(ii-2) 1からa2までの整数をランダムに選択(すべての整数が選ばれる確率は等確率)。
選ばれた整数をA2とする。
(ii-3) 
【数６６】
2006343124.tif 000081
とする。
【００５２】
(ii-2) bが上限a3を上回れば終了。そうでなければbに1を加え、(ii-1)に戻り、         さらに繰り返す。
【００５３】
以上のアルゴリズムにより擬似データの集合
【数６７】
2006343124.tif 000082
【００５４】
パラメータa1、a2、a3は下記表４のように設定した。
【００５５】
【表４】
2006343124.tif 000083
【００５６】
工程5
本工程では、センサ応答を入力とし、センサ応答に対応する複数の該化学物質の未知濃度を出力する階層型ニューラルネットワークANNquantificationの学習を行う。ANNquantificationの学習には、擬似データの集合
【数６８】
2006343124.tif 000084
を用いる。
【００５７】
[工程5-1 擬似データの前処理]
ANNquantificationに擬似データの集合
【数６９】
2006343124.tif 000085
を学習させる前に、前処理としてスケール変換を行い、
【数７０】
2006343124.tif 000086
を算出する。
【数７１】
2006343124.tif 000087
【００５８】
[工程5-2 階層型ニューラルネットワークによる擬似データの学習]
前処理済みの擬似データの集合
【数７２】
2006343124.tif 000088
をANNquantificationに学習させる。入力には
【数７３】
2006343124.tif 000089
を、出力への教師信号として
【数７４】
2006343124.tif 000090
を与える。学習はオンラインモードのバックプロパゲーションで行う。オンラインモードのバックプロパゲーションでは、一回の学習につきランダムに選ばれたi番目のデータに対する学習誤差
【数７５】
2006343124.tif 000091
を小さくするように、階層型ニューラルネットワークANNquantificationの結合加重wを以下の式のように調整していく。
【数７６】
2006343124.tif 000092
構造及び学習に関するパラメータは下の表５のように設定した。
【００５９】
【表５】
2006343124.tif 000093
【００６０】
ここで、
【数７７】
2006343124.tif 000094
は結合加重wを持つANNquantificationに
【数７８】
2006343124.tif 000095
を入力したときのANNquantificationの出力ベクトルを表す。
【００６１】
以上の学習によりANNquantificationの結合加重がw*に調整される。学習済みのANNquantification
の出力を
【数７９】
2006343124.tif 000096
表す。
【００６２】
工程6
本工程では、未知濃度の複数の化学物質を含む未知試料に対して、工程１と同様にセンサを適用し、センサ応答を測定する。
【００６３】
工程７
本工程では、「工程4」で擬似データを学習した階層型ニューラルネットワークANNquantificationを利用して未知資料のセンサ応答から未知試料に含まれる複数の化学物質の濃度を推定する。
【００６４】
[工程7-1階層型ニューラルネットワークによる未知濃度推定]
濃度が未知の試料に対するセンサ応答
【数８０】
2006343124.tif 000097
から、擬似データを学習した階層型ニューラルネットワークを利用して試料の未知濃度を推定する。まずセンサ応答を正規化するための前処理を下記の式のように行う。
【数８１】
2006343124.tif 000098
正規化されたセンサ応答
【数８２】
2006343124.tif 000099
を学習済みのANNquantificationに入力し、ANNquantificationの出力
【数８３】
2006343124.tif 000100
により、センサ応答に対する試料の正規化された未知濃度を推定する。gjはANNquantificationのj番目の出力ユニットからの出力である。下記の式のように、正規化された形から元のスケールになるよう変換を行う。
【数８４】
2006343124.tif 000101
そして、
【数８５】
2006343124.tif 000102
をi番目の化学物質の未知濃度として、複数の化学物質の未知濃度を推定する。
【００６５】
以下、本発明を実施例に基づきより具体的に説明する。もっとも、本発明は下記実施例に限定されるものではない。
【実施例】
【００６６】
工程１
(1) 試薬
本研究に用いた試薬は以下の通りである。
GOD: Aspergillus niger 由来Type II-S (EC 1.1.3.4) Sigma-Aldrich社製
ポリ (エチレングリコール)ジグリシジルエーテル (分子量400): Aldrich社製
アスコルビン酸(AA):和光純薬工業製
アセトアミノフェン(AAP): 和光純薬工業製
尿酸(          UA): 和光純薬工業製
リン酸緩衝食塩水 (PBS), pH7.4: Gibco BRL製
【００６７】
(2)  オスミウムメディエーターの作製
ポリ (1-ビニルイミダゾール)n-オスミウム（PVIn-Os）はAdam Hellerらの手法（Anal. Chem. 1993, 65, 3512-3517）に従って合成した。具体的には次のようにして合成した。すなわち、6 mlの1-ビニルイミダゾールと0.5 gの2,2'-アゾビス（イソブチロニトリル）をアルゴン気流下で２時間、70℃で加熱した。反応溶液を放冷後、生じた沈殿をメタノールに溶解し、その溶液をアセトンに攪拌しながら滴下した。溶液を濾過し、黄色い沈殿（ポリビニルイミダゾール）を濾取した。ポリイミダゾールとビス（2,2'-ビピリジン）ジクロロオスミウムをエタノール中で３日間還流し、オスミウムメディエーター（PVIn-Os）を合成した。
【００６８】
(3)  装置
電気化学測定には以下のポテンショスタットおよび電極を用いた。
ALS CH Instruments Electrochemical Analyzer Model 1000
φ3 mm PFCEカーボン電極（BAS社製）
Ptカウンター電極（BAS社製）
飽和KCl銀塩化銀参照電極（BAS社製）
【００６９】
(4)  グルコース電極の作製
メディエーターPVIn-Osおよび架橋剤ポリ (エチレングリコール)ジグリシジルエーテル（peg 400）をそれぞれ5 mg/ml，2.5 mg/mlになるよう、純水(Milli-Q(商品名）水）にて溶解した。また、GODは4 mg/ml になるよう、PBSにて溶解した。
【００７０】
研磨したカーボン電極表面にPVIn-Os溶液2μl，GOD溶液2μl，peg 400溶液1.2μlをマイクロシリンジにて滴下し、シリンジの先を用いて溶液をよく混合した。室温・遮光条件で48時間乾燥させ、電極表面上に酵素固定化膜を作製した。
【００７１】
(5)  測定
グルコースおよび共存妨害物質の電気化学測定は酵素を固定した作用電極、Pt対極、Ag|AgCl参照極を用いた３電極系によるCVにて行った。装置の設定値は以下のとおりであった。
最低電位: 0 V
最高電位: 0.6 V
スキャン速度: 0.5 V/s
試料間隔(sample interval): 0.001 V
感度: 1x10-5 A/V
【００７２】
CVは、具体的には次のようにして行なった。すなわち、電気化学測定用セルに測定溶液を入れ、ポテンショスタットにつないだ作用電極、対極、参照極を差し入れた。電極に溶液をなじませて３分間静止した後、ポテンショスタットにて作用電極に0-600 mVの間で1 mV間隔に観測される電流を検出した。
【００７３】
測定では、10 mlのPBSにグルコースと共存妨害物質の溶液を添加した後３分間静置して各濃度のサンプルに対する電流を測った。
【００７４】
下記の表６のような、センサ応答、複数の化学物質の組計30個をデータとしてサンプリングした。
【００７５】
【表６】
2006343124.tif 000103
【００７６】
より詳細には、前記複数の３０組のデータは、下記の表７に示す通りであった。
【００７７】
【表７】
2006343124.tif 000104
【００７８】
工程２〜７
工程２〜７は、先に詳述した方法により行なった。
【００７９】
結果
30個のデータを用い、leave-one-outクロスバリデーションを行った結果、以下の表８のような濃度推定結果が得られた。
【００８０】
【表８】
2006343124.tif 000105
【００８１】
表８の推定値と実験値の相違を平均相対誤差で評価した結果を表９に示す。
【００８２】
【表９】
2006343124.tif 000106
【００８３】
以上のように、本手法では少ないデータを利用した場合において、複数の化学物質の濃度を精度良く推定することができる。
【００８４】
工程2についてのソースコードを図１〜図１６に示す。工程3,工程4についてのソースコードを図１７〜図２６に示す。 (ただし工程3-3における標本共分散行列の逆行列の算出は図２７〜図３１、工程4-1のcholesky分解についてのソースコードは図３２〜図３６である。)工程5,工程6についてのソースコードを図３７〜図５１に示す。
【図面の簡単な説明】
【００８５】
【図１】工程２に用いたニューラルネットワークのソースコードを示す図である。
【図２】図１の続きを示す図である。
【図３】図２の続きを示す図である。
【図４】図３の続きを示す図である。
【図５】図４の続きを示す図である。
【図６】図５の続きを示す図である。
【図７】図６の続きを示す図である。
【図８】図７の続きを示す図である。
【図９】図８の続きを示す図である。
【図１０】図９の続きを示す図である。
【図１１】図１０の続きを示す図である。
【図１２】図１１の続きを示す図である。
【図１３】図１２の続きを示す図である。
【図１４】図１３の続きを示す図である。
【図１５】図１４の続きを示す図である。
【図１６】図１５の続きを示す図である。
【図１７】工程３及び４に用いたニューラルネットワークのソースコードを示す図である。
【図１８】図１７の続きを示す図である。
【図１９】図１８の続きを示す図である。
【図２０】図１９の続きを示す図である。
【図２１】図２０の続きを示す図である。
【図２２】図２１の続きを示す図である。
【図２３】図２２の続きを示す図である。
【図２４】図２３の続きを示す図である。
【図２５】図２４の続きを示す図である。
【図２６】図２５の続きを示す図である。
【図２７】工程3-3に用いたニューラルネットワークのソースコードを示す図である。
【図２８】図２７の続きを示す図である。
【図２９】図２８の続きを示す図である。
【図３０】図２９の続きを示す図である。
【図３１】図３０の続きを示す図である。
【図３２】工程4-1に用いたニューラルネットワークのソースコードを示す図である。
【図３３】図３２の続きを示す図である。
【図３４】図３３の続きを示す図である。
【図３５】図３４の続きを示す図である。
【図３６】図３５の続きを示す図である。
【図３７】工程３及び４に用いたニューラルネットワークのソースコードを示す図である。
【図３８】工程5及び工程6に用いたニューラルネットワークのソースコードを示す図である。
【図３９】図３８の続きを示す図である。
【図４０】図３９の続きを示す図である。
【図４１】図４０の続きを示す図である。
【図４２】図４１の続きを示す図である。
【図４３】図４２の続きを示す図である。
【図４４】図４３の続きを示す図である。
【図４５】図４４の続きを示す図である。
【図４６】図４５の続きを示す図である。
【図４７】工程3-3に用いたニューラルネットワークのソースコードを示す図である。
【図４８】図４７の続きを示す図である。
【図４９】図４８の続きを示す図である。
【図５０】図４９の続きを示す図である。
【図５１】図５０の続きを示す図である。
【図１】
2006343124.tif 000107
【図２】
2006343124.tif 000108
【図３】
2006343124.tif 000109
【図４】
2006343124.tif 000110
【図５】
2006343124.tif 000111
【図６】
2006343124.tif 000112
【図７】
2006343124.tif 000113
【図８】
2006343124.tif 000114
【図９】
2006343124.tif 000115
【図１０】
2006343124.tif 000116
【図１１】
2006343124.tif 000117
【図１２】
2006343124.tif 000118
【図１３】
2006343124.tif 000119
【図１４】
2006343124.tif 000120
【図１５】
2006343124.tif 000121
【図１６】
2006343124.tif 000122
【図１７】
2006343124.tif 000123
【図１８】
2006343124.tif 000124
【図１９】
2006343124.tif 000125
【図２０】
2006343124.tif 000126
【図２１】
2006343124.tif 000127
【図２２】
2006343124.tif 000128
【図２３】
2006343124.tif 000129
【図２４】
2006343124.tif 000130
【図２５】
2006343124.tif 000131
【図２６】
2006343124.tif 000132
【図２７】
2006343124.tif 000133
【図２８】
2006343124.tif 000134
【図２９】
2006343124.tif 000135
【図３０】
2006343124.tif 000136
【図３１】
2006343124.tif 000137
【図３２】
2006343124.tif 000138
【図３３】
2006343124.tif 000139
【図３４】
2006343124.tif 000140
【図３５】
2006343124.tif 000141
【図３６】
2006343124.tif 000142
【図３７】
2006343124.tif 000143
【図３８】
2006343124.tif 000144
【図３９】
2006343124.tif 000145
【図４０】
2006343124.tif 000146
【図４１】
2006343124.tif 000147
【図４２】
2006343124.tif 000148
【図４３】
2006343124.tif 000149
【図４４】
2006343124.tif 000150
【図４５】
2006343124.tif 000151
【図４６】
2006343124.tif 000152
【図４７】
2006343124.tif 000153
【図４８】
2006343124.tif 000154
【図４９】
2006343124.tif 000155
【図５０】
2006343124.tif 000156
【図５１】
2006343124.tif 000157