(57)【要約】【課題】  化合物の品質劣化予測を、非線形解析法（薄板スプライン法）を用いて精度良く行う。【解決手段】  実験値データファイル３２に記憶されたデータを読込（Ｓ０２）。薄板スプライン補間処理を行い、応答曲面モデルを生成（Ｓ０４）。予測値を得るために必要な条件データを入力（Ｓ０６）。条件データを、スプライン補間（薄板スプライン法）によって生成した各応答曲面モデルの式にＸ１〜Ｘ５の値だけを与え、予測値Ｙ１’、Ｙ２’を得る（Ｓ０８）。【選択図】  図１
2011021903.tif 000002
【特許請求の範囲】
【請求項１】
化合物の品質劣化を予測するための品質劣化予測方法であって、
予測対象である酸化反応及び加水分解反応による劣化を起こす化合物について、少なくとも時間を１つの入力指標として含む特定指標の値を実測して得た実験値データを蓄積する実験値蓄積工程、
前記実験値蓄積工程により得た実験値データに基づき薄板スプライン法を用いて非線形解析を行うことにより応答曲面モデルを生成し、所定条件で保存した後の品質劣化出力指標の値を当該応答曲面モデルに基づいて算出する出力指標算出工程、
を備えたこと、を特徴とする品質劣化予測方法。
【請求項２】
請求項１の品質劣化予測方法において、
劣化予測の対象となる前記化合物が、パルミチン酸、ステアリン酸、オレイン酸、リノール酸、リノレン酸のエチルエステル化合物、またはヨード化されたこれらの化合物類およびその混合物であること、を特徴とする品質劣化予測方法。
【請求項３】
請求項１の品質劣化予測方法において、
劣化予測の対象となる前記化合物が、ヨード化ケシ油脂肪酸エチルエステルであること、を特徴とする品質劣化予測方法。
【請求項４】
請求項１〜３の何れかの品質劣化予測方法において、
前記特定指標として、酸価、水分を含む、
ことを特徴とする品質劣化予測方法。
【請求項５】
請求項１〜４の何れかの品質劣化予測方法において、
前記特定指標として、酸素濃度を含み、
当該酸素濃度の値が所定範囲内に含まれる実験値データのみを用いて非線形解析を行う、
ことを特徴とする品質劣化予測方法。
【請求項６】
化合物の品質劣化を予測するための品質劣化データ算出プログラムであって、
コンピュータに、
予測対象である酸化反応及び加水分解反応による劣化を起こす化合物について、少なくとも時間を１つの入力指標として含む特定指標の値を実測して得た実験値データを蓄積する実験値蓄積処理、
前記実験値蓄積処理により得た実験値データに基づき薄板スプライン法を用いて非線形解析を行うことにより応答曲面モデルを生成し、所定条件で保存した後の品質劣化出力指標の値を当該応答曲面モデルに基づいて算出する出力指標算出処理、
を実行させることを特徴とする品質劣化データ算出プログラム。
【請求項７】
化合物の品質劣化を予測するための品質劣化予測装置であって、
予測対象の化合物である酸化反応及び加水分解反応による劣化を起こす化合物について、少なくとも時間を１つの入力指標として含む特定指標の値を実測して得た実験値データを蓄積する実験値蓄積手段、
前記実験値蓄積手段により得た実験値データに基づき薄板スプライン法を用いて非線形解析を行うことにより応答曲面モデルを生成し、所定期間保存した後の品質劣化出力指標の値を当該応答曲面モデルに基づいて算出する出力指標算出手段、
を備えたこと、を特徴とする品質劣化予測装置。
【請求項８】
クライアント端末からの要求に応じて、化合物の品質劣化を品質劣化予測サーバが予測するための品質劣化予測システムであって、
前記品質劣化予測サーバが、
予測対象の化合物である酸化反応及び加水分解反応による劣化を起こす化合物について、少なくとも時間を１つの入力指標として含む特定指標の値を実測して得た実験値データを蓄積する実験値蓄積データベース、
前記実験値蓄積データベースにより得た実験値データに基づき薄板スプライン法を用いて非線形解析を行うことにより応答曲面モデルを生成し、所定期間保存した後の品質劣化出力指標の値を当該応答曲面モデルに基づいて算出して、クライアント端末に出力する出力手段、
を備えたこと、を特徴とする品質劣化予測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、化合物（特に、酸化反応及び加水分解反応により劣化を起こす脂肪酸エステル化合物類）の品質劣化を予測するための技術に関する。
【背景技術】
【０００２】
食品、医薬品などに含まれる化合物は、経時的に品質劣化を起こすことが一般に知られている。しかし、品質劣化は化合物ごとに異なる種々の要因によって生じるため、将来における品質劣化の予測は容易ではなかった。
このため、油脂系化合物の経時的な品質劣化を予測するためには、保存温度や経過時間などを所定条件とした下で、製造した検体について品質劣化を示す特定の指標を実際に測定する必要があった。なお、対象物の品質劣化を予測する方法は従来から存在する（特許文献１、２）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００１−１３２４２５号公報
【特許文献２】特開２００５−１７１３２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、これら従来の劣化予測方法（特許文献１、２）は、品質劣化の要因である１つの指標に基づいて近似式を生成して経時的に劣化を判断するものに過ぎず、多要因によって劣化を生じる対象物には適用できなかった。また、多要因に基づく予測が可能なニューラルネットワーク法を用いた場合でも、油脂系化合物の品質劣化予測では良好な実験結果は得られなかった。
【０００５】
この発明は、多要因による複雑なメカニズムによって経時的に品質劣化を起こす化合物について、少なくとも時間を１つの入力指標として含む特定指標の値を測定した実測値データを薄板スプライン補間処理することにより応答曲面モデルを生成し、これに基づいて将来における化合物の品質劣化を精度良く予測することを目的とする。
【課題を解決するための手段】
【０００６】
（１）この発明の品質劣化予測方法は、
化合物の品質劣化を予測するための品質劣化予測方法であって、
予測対象である酸化反応及び加水分解反応による劣化を起こす化合物について、少なくとも時間を１つの入力指標として含む特定指標の値を実測して得た実験値データを蓄積する実験値蓄積工程、
前記実験値蓄積工程により得た実験値データに基づき薄板スプライン法を用いて非線形解析を行うことにより応答曲面モデルを生成し、所定条件で保存した後の品質劣化出力指標の値を当該応答曲面モデルに基づいて算出する出力指標算出工程、
を備えたこと、を特徴とする。
【０００７】
これにより、多要因による複雑なメカニズムによって経時的に品質劣化を起こす化合物について、将来の劣化予測を正確に行うことができる。また、多くのサンプルデータが得られないような場合でも、将来の劣化予測を行うことができる。
【０００８】
（２）この発明の品質劣化予測方法は、
劣化予測の対象となる前記化合物が、パルミチン酸、ステアリン酸、オレイン酸、リノール酸、リノレン酸のエチルエステル化合物、またはヨード化されたこれらの化合物類およびその混合物であること、を特徴とする。
【０００９】
これにより、多要因による複雑なメカニズムによって経時的に品質劣化を起こすエチルエステル化合物について、将来の劣化予測を正確に行うことができる。
【００１０】
（３）この発明の品質劣化予測方法は、
劣化予測の対象となる前記化合物が、ヨード化ケシ油脂肪酸エチルエステルであること、を特徴とする。
【００１１】
これにより、多要因による複雑なメカニズムによって経時的に品質劣化を起こすヨード化ケシ油脂肪酸エチルエステルについて、将来の劣化予測を正確に行うことができる。
【００１２】
（４）この発明の品質劣化予測方法は、
前記特定指標として、酸価、水分を含む、ことを特徴とする。
【００１３】
これにより、ヨード化ケシ油脂肪酸エチルエステルのような多要因による複雑なメカニズムによって経時的に品質劣化を起こす化合物についても、将来の劣化予測を正確に行うことができる。
【００１４】
（５）この発明の品質劣化予測方法は、
前記特定指標として、酸素濃度を含み、
当該酸素濃度の値が所定範囲内に含まれる実験値データのみを用いて非線形解析を行う、ことを特徴とする。
【００１５】
これにより、ヨード化ケシ油脂肪酸エチルエステルのような多要因による複雑なメカニズムによって経時的に品質劣化を起こす化合物についても、将来の劣化予測をより正確に行うことができる。
【図面の簡単な説明】
【００１６】
【図１】品質劣化のメカニズムを示す図である。
【図２】品質劣化予測装置のハードウェア構成を示す図である。
【図３】ＣＳＶファイルに記憶されたデータの具体例を示す図である。
【図４】品質劣化予測プログラムが行う処理のフローチャートである。
【図５】酸価および水分量の予測値、予測値を算出するために入力する条件データなどを示す図である。
【図６】酸価および水分量の予測値と実験値とを比較したグラフである。
【図７】クロスバリデーションによるモデル評価によって得られたデータの具体例を示す図である。
【図８】model0（全データ）により得られた予測値と、検体を実際に測定して得られた実験値との相関を示すグラフである。
【図９】model1（２５℃のデータのみ抽出）により得られた予測値と、検体を実際に測定して得られた実験値との相関を示すグラフである。
【図１０】model2（酸素濃度０．４％以下のデータのみ抽出）により得られた予測値と、検体を実際に測定して得られた実験値との相関を示すグラフである。
【図１１】model3（保存期間３６５日までのデータのみ抽出）により得られた予測値と、検体を実際に測定して得られた実験値との相関を示すグラフである。
【発明を実施するための形態】
【００１７】
[エステル化合物の品質劣化メカニズム]
まず、油脂系化合物の品質劣化予測の困難性について、ヨード化ケシ油脂肪酸エチルエステル（商品名：リピオドールウルトラフルイド（商標）、購入元はフランスのゲルベ（Guerbet）（商標）社）を例に、図１を用いて説明する。
【００１８】
なお、ヨード化ケシ油脂肪酸エチルエステルの経時的な品質劣化を判断する指標としては、酸価、水分、ヨウ素が挙げられるが、特に、酸価は品質劣化の基準が法律（医薬品や食品に関する）などで定められており重要である。また、水分についても、後述するように、酸価の生成プロセスに密接に関係しており、また、それ自体が造影剤の粘性を特徴付ける指標でもあるため品質劣化の指標として重要であると考えられる。
【００１９】
図１は、ヨード化ケシ油脂肪酸エチルエステルの品質劣化のメカニズムを示す図である。
【００２０】
ヨード化ケシ油脂肪酸エチルエステルには、製造時の不純物として、元々酸素Ｍ１や水分Ｍ２（図１に斜線で示す）が少量含まれている。
【００２１】
この酸素Ｍ１が、図１に示すように、ヨード化ケシ油脂肪酸エチルエステルＭ３と反応することでまず酸化物Ｍ４が生成される（図１に示す一次酸化Ｐ１）。一次酸化Ｐ１により生成された酸化物Ｍ４が、さらに酸化して過酸化物Ｍ５と水Ｍ６が生成される（図１に示す二次酸化Ｐ２）。また、二次酸化Ｐ２により生成された水Ｍ６が、さらにヨード化ケシ油脂肪酸エチルエステルＭ３と反応することによりカルボン酸（遊離脂肪酸）Ｍ７が生成される（図１に示す加水分解Ｐ３）。元々含まれる水Ｍ２によっても、同様に、ヨード化ケシ油脂肪酸エチルエステルＭ３と反応することでカルボン酸（遊離脂肪酸）Ｍ８が生成される（図１に示す加水分解Ｐ４）。
【００２２】
このようにカルボン酸（遊離脂肪酸）Ｍ７，Ｍ８が経時的に生成されることで、酸価が増加する。酸価とは、エステル（例えば、ヨード化ケシ油脂肪酸エチルエステル）に水が加わって加水分解することで生成されるカルボン酸（遊離脂肪酸）の量であり、カルボン酸が生成されると異臭を生じ、成分が分解される等の問題が生じる。つまり、カルボン酸の量が多くなると酸価の値が大きくなり品質劣化が進行していることになる。
【００２３】
また、ヨード化ケシ油脂肪酸エチルエステルに含まれる水分が増え、物性が変化（粘度が変わる等）してしまうことで、医薬品として投薬した場合に所望の効果が得られないことが有る。具体的には、リピオドールがＸ線の造影剤（酸価などで所定の基準が定められている）として他の薬と併せて用いられる場合に、懸濁液の物性が変わることで、例えばターゲットとしている局所に滞留させることができないという問題が生じる。
【００２４】
以上のように、ヨード化ケシ油脂肪酸エチルエステルの品質劣化は、多要因による複雑なメカニズムによって生じており、酸価、水分（水分量）、酸素濃度（酸素量）などの要素は、ヨード化ケシ油脂肪酸エチルエステルの劣化を判断する上で重要な指標となることが分かる。
【００２５】
[品質劣化予測装置のハードウェア構成]
図２は、この発明の品質劣化予測方法を実行する品質劣化予測装置１００のハードウェア構成を示す図である。図２に示すように、品質劣化予測装置１００は、ＣＰＵ１０、ＲＡＭ１２、ディスプレイ１４、ハードディスク１６、キーボード／マウス１８、記録媒体リーダー２０などを備える。
【００２６】
図１に示すように、ハードディスク１６には品質劣化予測プログラム（品質劣化データ算出プログラム）３０および実験値データファイル３２が記憶されている。なお、品質劣化予測プログラム３０のハードディスク１６へのインストールは、ＣＤ−ＲＯＭなどの記録媒体から記録媒体リーダー２０を介してデータを読み込ませて行われる。
【００２７】
品質劣化予測プログラム３０としては、少なくとも薄板スプライン法に基づく補間処理を与えられた実験データについて、応答曲面モデルを生成することが可能な市販のデータ解析用ソフトウェア（例えば、株式会社山武が販売するdataNESIA／visualNESIA（商標））を用いることができる。
【００２８】
実験値データファイル３２は、化合物について実際に測定した酸価（初期値）、保存時間などの実験値データを記憶したものである。図３に、実験値データファイル３２に記憶されるデータの例を示す。なお、実験値蓄積工程（処理、手段）とは、測定して得た上記実験値データを実験値データファイル３２に記憶する工程（処理、手段）をいう。  
【００２９】
図３に示すように、実験値データファイル３２には、各試料について実際に測定した、保存開始時の水分Ｘ１（mg/mL）、保存開始時の酸価Ｘ２（KOHmg/g）、保存開始時の酸素濃度Ｘ３（％）の値と、保存条件である保存温度Ｘ４（℃）、保存期間Ｘ５（日）の値（入力パラメータ）が記憶されており、さらにこれらの条件（保存温度Ｘ４（℃）、保存期間Ｘ５（日））で保存した後で実際に測定した水分Ｙ１（mg/mL）および酸価Ｙ２（KOHmg/g）の値（出力パラメータ）が記憶されている。このように、実験値データファイル３２に記憶されるデータには、特定の入力指標（入力パラメータ）として少なくとも時間（経過期間）が含まれている。なお、特定指標とは、応答曲面モデルを生成するために必要な特定の指標であり、例えば、上記Ｘ1〜Ｘ5の入力パラメータおよびＹ1、Ｙ2の出力パラメータが相当する。
【００３０】
なお、図３に示すデータ取得のために測定を行ったヨード化ケシ油脂肪酸エチルエステルは、天然のケシ油成分（パルミチン酸9%、ステアリン酸2%、オレイン酸19%、リノール酸68%、リノレン酸1%の混合物）のエチルエステル混合物につき、ヨード化（ヨウ素付加）した化合物である。
【００３１】
図３に示す酸価Ｘ２（カルボン酸の量）のデータを得るための測定は、ビュレットまたは市販の酸価測定装置を用いて、例えば、中和滴定法により行うことができる。水分X１のデータを得るための測定は、市販の水分測定装置を用いて、例えば、カールフィッシャー法により行うことができる。酸素濃度Ｘ３（酸素量）のデータを得るための測定は、市販の酸素濃度装置（例えば、ニードル式酸素濃度計）またはガスクロマトグラフィーを用いて行うことができる。
【００３２】
[品質劣化予測プログラム３０が行う処理]
図４を用いて、品質劣化予測プログラム３０が行う処理について説明する。
【００３３】
まず、ユーザーが、キーボード／マウス１８を介して、ディスプレイ１４上に表示された品質劣化予測プログラム３０を起動し、実験値データファイル３２のロード処理を実行するよう操作する。これにより、ＣＰＵ１０は、実験値データファイル３２に記憶されたデータを読み込む（ステップＳ０２）。
【００３４】
例えば、各試薬について測定された図３に示すＸ１〜Ｘ５およびＹ１、Ｙ２の値が全て読み込まれる。なお、実験値データファイル３２は、予めキーボード／マウス１８を介してデータを入力し、ＣＳＶ形式などに保存して生成しておく。
【００３５】
さらに、品質劣化予測プログラム３０は、ステップＳ０２で読み込まれたデータに基づいて薄板スプライン補間処理を行い、Ｙ軸に設定された変数（応答）Ｙ１、Ｙ２のそれぞれについて自動的に応答曲面モデルを生成する（ステップＳ０４）。つまり、図３に示す入力パラメータＸ１〜Ｘ５と出力パラメータＹ１について水分に関する１つの応答曲面モデルが生成され、図３に示す入力パラメータＸ１〜Ｘ５と出力パラメータＹ２について酸価に関する１つの応答曲面モデルが生成される。
【００３６】
薄板スプライン補間処理について、以下に説明する。なお、重調和スプライン補間は入力変数が増えて多入力になると、曲面モデル同定演算が発散するという欠点があるが、薄板スプライン（Thin Plate Spline：ＴＰＳ）補間を適用した場合には誤差を含む実験データであっても自然な補間が出来るため、これを回避することができる。
【００３７】
スプラインによる曲面生成の問題は、梁又は板のような弾性体の変形問題に置き換えることができる。すなわち、収集された実験データ（サンプル点）を弾性体のある部位に与えられた変位とみなせば、最も滑らかなスプライン曲面を求めることはある弾性体にいくつかの変位が与えられた場合に発生する下記式（１）に示す内部歪エネルギーＥを最小化する問題とみなすことができる。
【００３８】
【数１】
2011021903.tif 000003
上記式（１）を解く方法のうち、多変数データ、ランダムデータへの適用が容易なものとしてグリーン関数を用いた解法があり、これにて導かれたものが下記式（２）に示す重調和スプラインである。
【００３９】
【数２】
2011021903.tif 000004
ここで、ｎは与えられたデータ点数、ｄｉは、与えられたデータｉのＸ座標と任意のＸ座標とのユークリッド距離、ｇ（ｄｉ）は、ｄｉを変数とするグリーン関数であり、関数定義は、入力変数Ｘの次元数により異なる。例えば、２次元Ｘ＝（ｘ１，ｘ２）の場合、上記式（２）におけるグリーン関数は、下記式（３）となる。また、αｉは、係数であり、与えられたデータから線形マトリクス演算で計算することができる。
【００４０】
【数３】
2011021903.tif 000005
また、スムージング機能を有する薄板スプラインの場合は、最も滑らかなスプライン曲面を求めることは、下記式（４）のエネルギーＥ’式を最小化する問題とみなすことができる。
【００４１】
【数４】
2011021903.tif 000006
上記式（４）から薄板スプライン法による下記式（５）が導かれる。ここで、λはスムージングパラメータ、ｃｊは係数、ｐは入力変数の次元数を示す。なお、ｃ０は薄板スプライン補間における薄板部分を表す線形多項式の定数項である。
【００４２】
【数５】
2011021903.tif 000007
つまり、入力変数がＸ１〜Ｘ５で５次元の場合、ｐ＝５とし、上記式（５）の係数ｃｊ、定数項ｃ０などの値をコンピュータ演算により算出すれば、応答Ｙ１、Ｙ２それぞれについて応答曲面モデルが得られる。なお、スムージングパラメータλの値は、滑らかな曲面モデルが生成できるように調節される。
【００４３】
つぎに、ユーザーは、キーボード／マウス１８を介して予測値を得るために必要な条件データを入力する。これにより、ＣＰＵ１０は条件データを読み込む（ステップＳ０６）。例えば、図５に示すような予測値を算出したい条件（Ｘ１〜Ｘ５の値）だけを入力したデータ（条件データ）をＣＳＶ形式のファイルとして作成しておき、これを読み込むことでデータを入力することができる。
【００４４】
さらに、ＣＰＵ１０は、読み込んだ条件データを、スプライン補間（薄板スプライン法）によって生成した各応答曲面モデルの式に条件データ（Ｘ１〜Ｘ５の値だけ）を与える。これにより、各出力指標の予測値Ｙ１’、Ｙ２’が得られる（ステップＳ０８）。なお、品質劣化出力指標とは、化合物の品質劣化を表す出力指標のことであり、品質劣化出力指標の値とは、例えば、上記Ｘ１〜Ｘ５の条件データを各応答曲面モデルの式に与えることにより得られる上記水分の予測値Ｙ１’や酸価の予測値Ｙ２’が相当する。
【００４５】
ヨード化ケシ油脂肪酸エチルエステルの複数のロット（ロット１、２、３）について劣化予測を行った結果を、図５に示す。なお、図５Ａは、全データ（図３に示す実験値データファイル３２に記憶された全てのデータ）を用いて生成した応答曲面モデル（model0）を用いて劣化予測を行ったときの結果を示す図である。なお、ロット１、２、３の酸素濃度X3[0.4]、保存温度[25]、保存期間[365]は全て同じであるが、水分X1および酸価X2のパラメータがそれぞれ異なっている。
【００４６】
図５Ａに示すように、ロット１については条件データX1[0.52], X[0.4], X3[0.4], X4[25], X5[365]の値を応答曲面モデルに入力することで、図５に示す予測値Y1'[0.47], Y2'[0.52]が得られた。ロット１についての予測値Y1'[0.47], Y2'[0.52]は、実験値Y1[0.49], Y2[0.54]と比較して差が0.2しかないため比較的良好な結果と考えられる。
【００４７】
その他、温度（図３のＸ４）が２５℃のデータのみ全データから抽出して応答曲面モデルを生成した場合（model1）、保存期間（図３のＸ５）が３６５日までのデータのみ全データから抽出して応答曲面モデルを生成した場合（model2）、保存開始時の酸素濃度（図３のＸ３）が０．４％以下のデータのみ全データから抽出して応答曲面モデルを生成した場合（model3）についても同様の実験を行った。図５Ｂ、Ｃ、Ｄにそれらの実験結果を示す。
【００４８】
なお、図５Ａ〜Ｄでは、保存期間１２ヶ月での予測結果を示しているが、保存期間だけを変更した予測値が得たい場合、保存期間X5だけを所望の値に置き換えればよい。例えば、３６ヶ月後の劣化予測を行いたい場合には、図５に示す条件データのうちX5の値「365」だけを「1095」に置き換えればよい。
【００４９】
図６は、図５Ａ〜Ｄに示す酸価及び水分量の予測値と実験値の相関を示すグラフである。図６に示すように、何れのモデル（model0-3）についても予測値、実験値についてプロットした点が正解を示す線に沿っており、比較的良好な劣化予測が行えたことが分かる。
【００５０】
[クロスバリデーションによる評価]
つぎに、クロスバリデーションにより各モデルの評価を行った結果を、以下に図７〜図１１を用いて説明する。
【００５１】
クロスバリデーションとは、モデルの汎化性を調べるための方法の1つであり、n個のサンプル試料から１個のサンプル試料を除いたn-1個のセットを用いて，応答曲面モデルの生成を行い，除いた1個のサンプル試料を用いて得られた予測値によりモデルの評価を行うことである。
【００５２】
例えば、図３の実験値データファイル３２に記憶されたデータのうち、試料No.1のデータ以外のデータ（サンプル試料の総数n=125より、No.2〜No125までの１２４個）を用いて応答曲面モデルをまず生成し、これに試料No.1の条件データX1[0.87], X[0.42], X3[0], X4[40], X5[30]を入力することで、図７に示す予測値Y1'[0.80], Y2'[0.53]が得られる。No.2〜No125の試料についても順に同様の処理を行っていくことで、図７に示すようにすべての各試料について予測値Y1', Y2'が得られる。なお、図７に示すＬＯＯは「Leave-One-Out」の略語である。
【００５３】
図８〜図１１は、図５Ａ〜Ｄに示す各モデル（model0-3）についてクロスバリデーションを行った結果、各試料について得られた予測値と実験値の相関性を示すグラフである。
【００５４】
全データを用いて応答曲面モデルを生成した場合（model0）、図８に示すような結果となった。クロスバリデーションにより得られる決定係数Ｒ２は、水分については０．９１４９であり、酸価については０．６５５７である。これより、水分の予測精度が良いことが分かった。
【００５５】
温度（図３のＸ４）が２５℃のデータのみを全データから抽出して応答曲面モデルを生成した場合（model1）、図９に示すような結果となった。クロスバリデーションにより得られる決定係数Ｒ２は、水分については０．８６８３であり、酸価については０．８３３７である。これより、抽出するデータの温度範囲を２５℃に限定した場合には、水分と酸価の何れについても予測精度が比較的良好であることが分かる。
【００５６】
保存期間（図３のＸ５）が３６５日までのデータのみを全データから抽出して応答曲面モデルを生成した場合（model2）、図１０に示すような結果となった。クロスバリデーションにより得られる決定係数Ｒ２は、水分については０．９４２１であり、酸価については０．７５８３である。これより、抽出するデータの保存期間の範囲を３６５日までに限定した場合には、特に水分の予測精度が良好であることが分かる。
【００５７】
保存開始時の酸素濃度（図３のＸ３）が０．４％以下のデータのみを全データから抽出して応答曲面モデルを生成した場合（model3）、図１１に示すような結果となった。クロスバリデーションにより得られる決定係数Ｒ２は、水分については０．９６６５であり、酸価については０．８３０３である。これより、抽出するデータの酸素濃度の範囲を０．４％以下に限定した場合には、特に水分の予測精度が良好であり、また、酸価についても予測精度が比較的良好であり、全てのモデルと比較しても総合的に最も良い結果が得られることが分かる。
【００５８】
このように、応答曲面モデルを生成するために抽出するデータを所定の数値範囲に限定することで、全データを用いて応答曲面モデルを生成した場合（model0）よりも、特定指標についての予測精度の改善がみられた。特に、酸素量０．４％以下のデータのみを用いた場合（model3）に、酸価の予測精度が大きく改善した。
【００５９】
[その他の実施形態]
なお、上記実施形態では、ヨード化ケシ油脂肪酸エチルエステルの劣化予測を行うこととしたが、これに限定されるものではなく、酸化反応及び加水分解反応による劣化を起こす脂肪酸エステル化合物類（混合物を含む）であれば同様の手法により劣化予測を行うことが可能である。
【００６０】
具体的には、脂肪酸（パルミチン酸、ステアリン酸、オレイン酸、リノール酸、リノレン酸など）のエチルエステル化合物が挙げられ、これらの化合物類にはヨード化されたもの（及びその混合物）も含まれる。
【００６１】
なお、上記実施形態では、応答曲面モデルに条件データを与えて酸価および水分量などの予測値を得ることにより品質劣化評価を行うこととしたが、最適値探索プログラム（例えば、株式会社山武が販売するdataNESIA／designNESIA（商標））を用いて所望の品質劣化条件（酸価および水分量など）を満たすような保存条件（温度、期間など）の値を得るようにしてもよい。
【００６２】
なお、上記実施形態では、品質劣化の予測をスタンドアローンのパーソナルコンピューターにより行ったが、これに限定されるものではなく、インターネットなどのネットワークを介して品質劣化の予測を行うサーバから端末に対して行うようにしてもよい。例えば、医薬品を購入した顧客に対して、製造番号の入力を受けた医薬品販売会社のサーバから品質劣化の予測結果をクライアントの端末に出力するよう構成してもよい。
【符号の説明】
【００６３】
１００・・・・品質劣化予測装置
３０・・・・品質劣化予測プログラム
３２・・・・実験値データファイル
【図１】
2011021903.tif 000008
【図２】
2011021903.tif 000009
【図３】
2011021903.tif 000010
【図４】
2011021903.tif 000011
【図５】
2011021903.tif 000012
【図６】
2011021903.tif 000013
【図７】
2011021903.tif 000014
【図８】
2011021903.tif 000015
【図９】
2011021903.tif 000016
【図１０】
2011021903.tif 000017
【図１１】
2011021903.tif 000018