2018年3月5日月曜日

<文献検索から>(28) 予測値の不確実性推定を伴う近赤外分光法と部分最小二乗法による水質モニタリング


題名:Monitoring Process Water Quality Using Near Infrared Spectroscopy and Partial Least Squares Regression with Prediction Uncertainty Estimation

著者P. B. Skou, T. A. Berg, S. D. Aunsbjerg, D. Thaysen, M. A. Rasmussen, and F. van den Berg

出典:Applied Spectroscopy 2017, Vol. 71(3) 410–421
膜濾過システムは、乳製品製造のいくつかの工程で使用され、製品および成分の様々なサイズの分子を分離する。膜ろ過技術を使用した乳製品の加工方法は、今後数年で大きく成長すると予測されている。チーズホエイ(Whey)タンパク質は、異なる機能特性を有する広範囲の製品に使用することが出来る。
  <Wheyとは? 全乳または脱脂乳からチーズを作る際に得られる液体。製菓・製パンなどの原料や 家畜飼料として利用される。乳清(にゆうせい)。>
  酪農産業の水使用量は、生産工程または浄化プロセスの下流のプロセス水の効率的な再利用によって大幅に削減できる。 ホエイ及びミルクは、多量の水分を含み、プロセス水の供給源として利用できる可能性がある。しかし、微生物学的な安全性と衛生性は食品や酪農業界にとって非常に重要である。乳糖、尿素、塩類などの有機化合物は逆浸透膜を通過することが報告されている。 プロセス水工程にプロセス分析技術(PAT: Process Analytical Technology)の考え方を適用することにより、リスクを継続的かつリアルタイムに監視することができる。尿素は、逆浸透膜透過液中の主要な化合物としてGC-MS分析によって同定されているが、時間がかかり、測定に労力が必要である。尿素は、近赤外分光法(NIRS)によっても測定可能で、迅速で、労力も少なく、または自動化測定も可能である。
  尿素および乳糖は、水の吸収ピークのはざまのスペクトル範囲20802325nmに特徴的な吸収パターンがある。チーズ製造工程で予想される尿素濃度は100600ppmであり、低濃度と水の吸収ピークの干渉は、検出限界を上げて、予測不確実性を増加させる。 部分最小二乗回帰(PLS)などの多変量モデリング手法は、単変量解析法に比べて予測誤差と検出限界を低くすることができる。PLSに代表される多変量解析法の予測不確実性の信頼できる推定値がなければ、検出限界を確立することはできず、リアルタイムのリスクアセスメントは実現できない。
  リアルタイムモニタリング用の検量線開発をする場合、代表的な下記の2つの方法がある。(1)測定対象となる既知濃度の試薬と測定に潜在的に悪影響を及ぼす溶液を混合し近赤外分析計で測定する。(2)プロセスから収集した実サンプルについてラボで測定し、同時に測定した近赤外分析計のスペクトルと関連付ける。ここでは、酪農プロセス水における尿素濃度予測のための2つの検量線作成法を比較し、予測不確実性を推定する方法を提示する。最後に、NIRSPLSを組み合わせたプロセス水のリアルタイムリスクアセスメントの可能性を示す。
検量線モデルの予測値分散は下記のようにあらわすことが出来る。
 
各サンプル予測値の分散=検量線誤差+モデル誤差+機器誤差+基準値誤差 ……(1)
注)検量線誤差は、モデル、基準値、機器等を含んだ検量線計算結果のバラツキ
  モデル誤差は、検量線の平均二乗誤差から機器誤差と基準値誤差を引いたものと定義。

  プロセス水サンプルは、本格的な生産システムで膜透過液から収集され、可視/紫外分光光度計を用いて340nmでの吸収から濃度測定を行った。ラボサンプルは尿素(9溶液)およびラクトース(3溶液)および脱塩水を様々な比率で組み合わせることによって作成した。ラボサンプルとプロセスサンプルから尿素濃度予測用の2種類の検量線をPLSにより作成した。スペクトルのバラツキは、5回の繰り返し測定スペクトルから平均値に近い3スペクトル及び3日間に測定した同じサンプルの測定スペクトル(ラボサンプルのみ)より計算した。
  ラボサンプルを使用すると何度も繰り返し測定が可能で長期間の再現性を得ることが出来る。一方繰り返し性、再現性が実際のプロセスのものと異なるという危険性がある。プロセスデータでは同じサンプルでも採取日が異なると変成する可能性がある。2つのサンプルから得られた検量線は異なる予測値と不確かさを持つ。検量線作成と検定結果は、検量線の予測誤差はラボサンプルを使用した検量線の方が良い(RMSEC:9.7/10.2ppmLab/Process))が、実プロセスでの検定結果では、プロセスサンプルの方が良い(RMSEP: 13.0/12.1ppmLab/Process))ことを示している。
予測値誤差分散は、(1)式に各誤差を入れて求められた。平均予測誤差(RMSEP)と平均予測不確実性(SPE)を比較すると、プロセスモデルよりラボモデルの方が、差が大きくなった。下表参照

 


RMSEP(ppm)


SPE(ppm)


ラボ


13


15.1


プロセス


12.1


12.4

    原因としては、プロセスサンプルを用いた検証時に予測値が基準値に対し少し偏りを持つためと推測する。SPEは、RMSEPと非常に近いことが分かる。これは、PLSモデルが合理的に作成されていれば、RMSEPSPEを正確に推定できることを示唆している。作成した検量線を限外濾過膜、逆浸透膜、逆浸透膜ポリッシャーの3種の膜を通した後のプロセス水のスペクトルに適用し尿素を測定した。限外濾過膜使用サンプルでは、計算から得られた予測誤差分散を外れることがあったが、その他の濾過膜では予測誤差分散内であった。プロセスモニタリングのための予測モデル構築には、ラボとプロセスサンプルの両方を組み合わせて最良の予測性能が得られることが多い。予測性能向上のためには、ラボサンプルを作り直しプロセスサンプルを長期にわたり収集し追加することで予測性能が向上する。




2018年1月7日日曜日

<文献検索から>(27)近赤外分光法とケモメトリックスモデルを用いたシプロフロキサシン錠の分類


題名:Classification of Ciprofloxacin Tablets Using Near-Infrared Spectroscopy and Chemometric Modeling
著者:N. Fuenffinger, S. Arzhantsev, and C. Gryniewicz-Ruzicka
出典:Applied Spectroscopy 2017, Vol. 71(8) 1927–1937

 
出所不明な医薬品、ラベル間違い、偽造医薬品等は、使用者に重大なリスクをもたらす世界的な問題である。米国の個人が摂取した薬物の約40%は、他国で製造されていて、誤ラベル等による使用者への影響が懸念されている。米国食品医薬品局(FDA)は、医薬品の品質と認証を評価するための迅速で信頼できるスクリーニング方法を開発している。医薬品指紋認証は、製品を特定し、その製品を製造元に追跡し、製品品質の変化を検出するための方法の開発を必要とする。2008年、FDAは、イオン移動度分光法、X線蛍光分光法、ラマン分光法、近赤外(NIR)分光法などの技術を使用して、米国への医薬品輸入の迅速なスクリーニング方法の開発を目的としたプロジェクトを開始した。これまでの研究では、ニューラルネットワーク、SIMCA soft independent modeling)法、K-最近傍法などのパターン認識分類手法を高速液体クロマトグラフィー測定に適用することによって医薬指紋認証をすることが可能であることが示されている。

NIR分光法は、医薬品を安全で迅速かつ効果的にスクリーニング可能なコストパーフォーマンスの高いツールである。この研究の目的は、いくつかの製造業者からのシプロフロキサシン塩酸塩錠剤の分類に関するNIRおよびケモメトリクスの能力を決定することである。シプロフロキサシンはヒト毒性の低い合成フルオロキノロン抗生物質であり、ブドウ球菌、肺炎連鎖球菌、ヘモフィルスおよび緑膿菌などの広範囲の生物に有効であることが証明されている。
米国およびインドの6つの製造業者から53ロットにわたる1025のシプロフロキサシン錠剤のNIR吸収スペクトルを収集し評価した。NIRスペクトルは、フーリエ変換NIRアナライザーを用いて 1製品あたり合計24の錠剤を分析した。各錠剤からのスペクトルは、8.0cm -1の分解能で400010,000cm -1の領域にわたって32回のスキャンを平均することによって計算した。前処理は、Savitzky-Golayアルゴリズムを使用して一次微分しSNVStandard Normal Validate)変換を行った。

製薬指紋認証の究極の目標は、(1)製造業者を特定することができるモデルの開発、 (2)医薬品の投与量を特定 (3)工程または原材料の変化に起因するロット間の変動、を評価することである。ここでは、最初の2つの目標のいずれかがシプロフロキサチン錠剤のNIRスペクトルから達成可能であるかどうかを特定するための客観的方法として、教師なし分類モデルの階層クラスター分析(Hierarchical Cluster Analysis : HCA)による分析を行った。HCAと二乗ユークリッド距離測定から得られた結果は、製造業者間でNIRスペクトルの最大の違いが生じることを示している。これらの知見に基づいて、シプロフロキサシン錠剤の成分の違いを予測する目的で、主成分分析の主成分を用いてデータ量を縮小した後に二次判別分析(QDA)モデルを構築した。 QDAを使用して、正確に96%以上の精度で907錠の錠剤を正確に分類することができた。

2017年11月30日木曜日

<文献検索から>(26) 医薬品農業サンプルの定量分析へのMCR-ALSの応用


題名:Application of multivariate curve resolution alternating least squares (MCR-ALS) to the quantitative analysis of pharmaceutical and agricultural samples

著者:T. Azzouz, R. Tauler
出典:Talanta 74 (2008) 1201–1210

この研究では、2つのサンプルの定量分析が検討されている。医薬品中の混合成分のUV-vis分光による分析と近赤外分光(NIR)による飼料サンプルの分析であるUV-vis吸収測定は、迅速で安価で医薬品製剤の分析に適しているが、異なる成分の吸収帯が大きく重なっているために一般的な適用に限界がある。医薬製剤は、通常、医薬原薬成分と、同じスペクトル領域で吸収する種々の賦形剤との混合物である。重なり合ったスペクトルの波形分解を基礎にした多変量解析法の開発は、迅速なスペクトル分解と定量化を可能にする。一般的に、多成分試料中の目的の分析対象物すべてについて完全に選択的な分析信号を得ることは非常に困難である。そのため、クロマトグラフィーまたは他の分離分析法による物理的分離、或いは多変量解析を用いた数学的分解法は、特に複雑な天然試料の分析における定量分析に有用である。部分最小二乗法(PLSRのような多変量解析法とは、別に種々のケモメトリック法が存在する。この研究で提案されているMCR-ALSは、分光学的手段を用いて複雑な混合物の定量分析が可能である。

多変量波形分解 交互最小二乗法(MCR-ALS法)
第1Step:予測した成分数の初期濃度ベクトルを作成する。
    (初期条件として、構成成分の濃度情報又は純粋スペクトル情報を必要とする)
2 Step:測定データのマトリックスと、上記濃度の疑似逆行列から予測スペクトル
         マトリックスを得る。
3 Step:予測スペクトルマトリックスと測定データマトリックスより
         疑似成分予測ベクトルを得る。
    あらかじめ設定した拘束条件を満たすまで、順次表記ステップを繰り返す。

拘束条件
   ①非負濃度拘束……成分濃度がマイナスになることはない。
   ②非負スペクトル拘束……スペクトルの吸光度がマイナスになることはない。
   ③相関拘束……ある特定した成分の予測濃度が、初期の濃度と一定の相関がある
                という仮定(この論文で提案された拘束)

部分最小二乗回帰(PLSR
  この回帰方法は、分光データ(データ行列)および濃度ータ(濃度ベクトル)の両方
  に含まれる情報を効率的に使用して、サンプル中の分析対象物濃度の予測モデルを構築
  する。PLSRアルゴリズムは、データ行列と濃度ベクトルとの共分散を最大にする連続す
  る直行因子を選択する。

【実験結果】
<医薬品中の混合成分>
  市販薬MicrogynonNeogynonaのエチニルエストラジオール卵胞ホルモン、前立腺
  癌薬)とレボノルゲストレルステロイド(避妊薬)UVスペクトルをMCR-ALSで解
 析し濃度値をPLSRの結果と比較した。使用測定波長は250-300nmが最適であった。得ら
 れた誤差は、 PLSRが少し良い値を示したが、MCR-ALSPLSRでほぼ同じオーダーで
 あった。同じ成分が入っている褐色の市販薬Triagynonについても定量分析を行い同様の
 結果を得た。MCR-ALSによって得られたスペクトルは、エチニルエストラジオールおよ
 びレボノルゲストレルの純粋な標準のスペクトルと一致する。対照的に、PLS回帰は、
 混合物中に存在する成分の濃度を予測するが混合物の成分の純粋なスペクトルの推定は
 できない。MCR-ALSの明らかな利点は、成分の純粋なスペクトルを含めて定性的情報を
 得られるという点にある。
NIRで測定したライグラス(Ray-Grass:細麦)の湿度とタンパク質含量の定量
  ライグラス中の水分とタンパク質の定性/定量的情報を得るためにMCR-ALSPLSRを適
  用し結果を比較した。水分測定に関しては、PLSRとほぼ同等(PLSRの方が若干良い)
  であったが、タンパク質については、PLSRは明らかにMCR-ALSより優れていた。MCR-
    ALSモデル比較してPLSR多数のコンポーネントFactor使うことが出来るため
  と考えられる。 MCR-ALSによって得られたスペクトルは、飼料中の湿度とタンパク質に
 関する文献で以前に報告されたNIRスペクトルによく似ていた。

【結論】
  未知の混合物および天然試料中の干渉のあるサンプルに新しい相関制約を適用する
  MCR-ALSの能力は、一般に、PLSRを用いて得られた結果と同等であった。
  MCR-ALS使用する主な利点は、測定対象成分と未知物質の差異についての
  スペクトル情報が得られることである。

2017年9月10日日曜日


<文献検索から>(25) 量子化学計算に裏付けされたローズマリー葉に含まれるロズマリン酸の近赤外分光法と赤外全反射分光法による定量測定の評価

題名:Critical Evaluation of NIR and ATRIR Spectroscopic Quantifications of Rosmarinic Acid in Rosmarini folium Supported by Quantum Chemical Calculations
著者:C. G. Kirchler, C. K. Pezzei, K. B. Beć, R. Henn, M. Ishigaki, Y. Ozaki, C. W. Huck1
出典:  Planta Medica April 2017

 用語:ロズマリン酸
ローズマリー、シソなどのシソ科の植物に含まれているポリフェノール。ロズマリン酸は低分子ながらタンニン活性を示す物質としても知られている。ロズマリン酸の抗酸化作用によりアレルギー症状を抑える効果がある。タンニン:タンパク質などと結合する水溶性ポリフェノール
 
【抄録】
本研究では、ローズマリー葉に含まれるロズマリン酸含量の測定における近赤外(NIR)および全反射赤外分光(ATR)の性能を評価した。さらに、量子化学計算を用いて、ロズマリン酸の理論的な赤外スペクトルを得て、測定スペクトルと比較した。
「漢方薬」又は「伝統医学」での生薬の重要性は過去数十年に著しく増加している。2004年以来、欧州医薬品庁(EMA)の一部であるHMPC European Medicines Agency's Committee on Herbal Medicinal Productsは、ハーブの物質、調製物、および組み合わせのためのモノグラフを公開している。 本研究の対象であるローズマリー葉について、HMPCは、胃腸および消化器障害の軽度の痙攣性疾患の症状緩和のための治療に適応するとしている。ローズマリーの葉は約1.12.5%のロズマリン酸(RA)を含む。RAの酸化防止効果は、そのフェノール酸またはヒドロキシケイ皮構造に基づいており、ラジカル消去活性はフェノール性ヒドロキシル部分によると推論されている。 RAの主な効果は、抗ウイルス、抗菌、抗酸化、抗炎症等である。
赤外分光法は、光と物質との相互作用に基づいており、分子振動の励起中の双極子モーメントの変化を必要とする。 MIR領域(4000400cm-1)は、基本振動による吸収によって主に支配される。特にC = OまたはC-Oのような極性基の基本的な振動遷移が最も顕著である。 一方、NIR領域(128214000cm-1)では、吸収は一般にC-HO-H、およびN-H部分の倍音振動および合成振動に由来する。MIRおよびNIR分光法の単純化モデルは、両方の技術の大きな違いを表している。 したがって、本研究では、ATR-IRNIR分光法のローズマリー葉のロズマリン酸(RA)含量を定量する能力の比較と評価を扱う。
粉砕されたローズマリー葉42サンプルは、HPLC-DAD (Diode array detector) 分析で測定されたロズマリン酸標準値(1.1382.199%)と近赤外(NIR)と赤外全反射(ATR-IR)スペクトルとの相関を、部分最小二乗法(PLS)を使用して調べた。HPLC-DAD分析の結果は、RA濃度範囲:1.138-2.199 %で相関係数R2; 0.9997LOD (Limit of Detection)0.067 ; LOQ (Limit of Qualification)0.120 %であった。
クロスバリデーション(CV)とtest-setを使用した評価は、NIRATR-IRで(NIR: test-setによる相関係数0.9CVによる標準誤差:0.06ATR-IRtest-setによる相関係数0.91CVによる標準誤差:0.063)とほぼ同等の結果となった。
測定されたスペクトルをより深く理解するために、量子力学的計算によってRAの理論的IRスペクトルを計算した。RA分子の構造およびその後の振動分析は、ハイブリッド密度汎関数により、信頼できるスペクトルの波数と強度を計算できることが過去の文献に示されているのでDFT密度汎関数)理論で行った。低波数域で純粋なロズマリン酸のスペクトルと良好に一致したが、波数が高い領域では一致度が劣った。




 

 

2017年4月23日日曜日

<文献検索から>(24) インライン発酵モニターのための赤外光ファイバーと鏡面反射光学系を使用した減衰全反射(ATR)フーリエ変換赤外(FT-IR)の比較


題名:Comparison of Fiber Optic and Conduit Attenuated Total Reflection (ATR) Fourier Transform Infrared (FT-IR) Setup for In-Line Fermentation Monitoring
著者:Cosima Koch, Andreas E. Posch, Christoph Herwig, and Bernhard Lendl
出典:Applied Spectroscopy 2016, Vol. 70 (12) 1965–1973

【抄録】
    In-situ用の全反射(ATR)プローブの開発に伴い、バイオプロセス分析およびモニタリングのための中赤外領域(400-4000cm-1)フーリエ変換赤外(FT-IR)分光法は、バイオプロセスのインラインモニタリングおよび制御のための重要なプロセス分析ツール)となっている。 ATRプローブは光がサンプルに浸み込む深さが数マイクロメートルと限られているため、通常は発酵培地の主成分である水の吸収の影響で使用が制限されるが、部分最小二乗回帰(PLSR)等の多変量解析モデルを使用できる信頼性の高い分析システムが利用可能であれば、可溶性サンプルの定量が可能である。
    分光器とATRプローブを接続する2つの方法(中空光導管(光導管)と光ファイバー)がある。光導管は、ミラーを調整することによって光のスループットを最大にする必要がある。信頼性の高いスペクトルを得るためには調整位置に変化を加える振動等を避ける必要がある。インライン発酵モニタリングのために光導管プローブを使用する場合、攪拌器からの振動がFT-IR測定を損なう可能性があるため注意が必要である。光ファイバー内では、光は全反射によって透過され、透過範囲は材料特性に依存する。高波数限界は電子吸収(バンドギャップエッジ)によって決定され、低波数限界はマルチフォノンエッジに依存する。多結晶ハロゲン化銀(AgX)は、現在、FT-IR分光プローブのために最も一般的に使用される光ファイバー材料で指紋領域(1500-900cm -1)の測定が可能である。光ファイバーATRプローブを用いたバイオプロセスモニタリングに関する研究はほとんど発表されていない。
    ここでは、インラインバイオプロセスモニタリング用のプロセス分光器に接続された光ファイバー付きATRプローブ(以下RIR15と光導管を用いた乾燥空気パージ付ATRプローブ(RIR45)のシステム性能を比較する。 両方のシステムを同じバイオリアクターに連結し、4つのPenicillium chrysogenumペニシリウムクリソゲナム:青かび)発酵を行い、生成物(ペニシリンV、以下PenV)およびその前駆体(フェノキシ酢酸、以下POX)をPLSRと組み合わせたFT-IR分光法と高速液体クロマトグラフィー(HPLC)によってモニターした。光導管と光ファイバーを使用した ATR FT-IRシステムのそれぞれの利点と潜在的な短所について議論する。
   青カビ菌株4種のフェドバッチ発酵が、プロセス制御と管理システムが統合されたシステムを用いて15Lのオートクレーブを持つ攪拌型バイオリアクターによって実行された。 RIR15RIR45の両方ともバイオリアクターに挿入された。
   PenVPOX濃度は、バイオイアクターに組み込まれたアイソクラティックHPLCによって測定された値を検量線の基準値として使用した。 サンプリングレートはPenVPOX濃度の変化に合わせて調整され、緩やかな変化が予想される場合、3時間ごとに行われたが、POXの添加中にはサンプリング間隔は16分(クロマトグラム実行時間)に短縮された。各々の構成について、オンラインHPLCサンプリング時間に最も近い3つのスペクトルを平均し、検量線作成に使用した。前処理はMean Centering1次微分で検量線はPLSRを使用した。
   RIR45の全体的な性能はRIR15の性能よりも良く、1400-1200cm-1のスペクトル領域において、RIR15RIR45よりもRMSおよびPPノイズがそれぞれ約2倍大きく、予測値もわずかに低いRMSEPが得られた。RIR15で推定された濃度は、PenVPOXの両方のRIR45よりもスペクトル毎のバラツキが大きい。全体的に、二つのシステムは、HPLCの値を基準にして約0.2g L_1以内のPOXおよびPenV濃度推定が可能であり、RIR45が相対的に良い結果となった。 バイオプロセスモニタリングのためのシステムを選択する際には、設置およびメンテナンスの容易さ(乾燥空気パージの有無等)と測定精度とコストとの間のトレードオフがなされなければならない

2017年2月23日木曜日

<文献検索から>(23)2D-COSおよびPCMW-2D解析を用いたベンチトップFTIRとミニチュアハンドヘルド近赤外分光器を用いたポリアミド11の温度変化の比較研究


題名:Comparative Variable Temperature Studies of Polyamide II with a Benchtop Fourier Transform and a Miniature Handheld Near-Infrared Spectrometer Using 2D-COS and PCMW-2D Analysis

著者:Miriam Unger, Frank Pfeifer, and Heinz W. Siesler
出典:Applied Spectroscopy 2016, Vol. 70(7) 1202–1208
【抄録】
  この論文の主な目的は、小型ハンドヘルド近赤外(NIR)分光器とフーリエ変換近赤外線(FT-NIR)分光器の性能を比較することである。ポリウレタンおよび脂肪族ポリアミドの場合、温度上昇の関数としての水素結合の変化は、それらの熱的および機械的特性に大きな影響を及ぼす。この論文で調査されたポリアミド11PA11)の可変温度NIRスペクトルの分析は、カルボニルとNH官能基との間の水素結合の温度変化に関する重要な情報を提供する。外部摂動によって誘発されるスペクトル変化は、2次元相関分光法(2DCOS)および摂動相関移動窓二次元(PCMW2D)分析法によって分析された。同期2DCOSは、外部摂動(ここでは温度の上昇)の変化の間の2つの別個のスペクトル間の全体的な類似性または同時の変化を表す。非同期2DCOSは、同期して相関されていないスペクトル強度の変動を表す。PCMW2D相関法は、スペクトル変数(ここでは波数)軸と摂動(ここでは温度)軸との2次元マップにプロットされた同期スペクトルと非同期スペクトルのペアを提供する。両方の技術の重要な利点は、外部摂動中の重なり合ったバンドの複雑なスペクトル強度変動を検出することができ、これらのスペクトル強度変動を引き起こす現象の連続的な順序を決定できることである。
  FT-NIR分光器と、小型NIR分光器の性能を比較するために、半結晶PA11シートの可変温度スペクトルを同じ実験条件下で測定し、2D-COSおよびPCMW2Dを使用して分析、評価した。60507500cm-1の波数域がFT-NIRおよび小型NIR分光器でそれぞれ測定された。この波数範囲は、2D-COSおよびPCMW2Dを使用して温度区間30-190℃における脂肪族セグメントの再結晶化および水素結合構造の変化に関する情報を抽出するために選択された。この波数域は、両方の分析器で利用可能なため、それぞれの2D-COSPCMW2Dの結果を直接比較することができる。
スペクトルで観察される最も顕著な強度変化は、水素結合および遊離NH官能基(64006800cm -1)の倍音(NH)領域で起こる。スペクトル分解能の大きな違い (データポイント間隔:4cm-143cm-1 ) にもかかわらず、68006520cm -1間および65206450cm -1間の両方のPCMW2Dは、水素結合種を犠牲にしてNHフリー官能基が生成することを示している。
  PCMW2Dデータとは対照的に、FT-NIRのより高いスペクトル分解能は、小型NIR分光器と比較してはるかに複雑な2DCOSクロスピークを示す。 FT-NIRでは同期および非同期の2D相関スペクトルにおいて、それぞれ3つおよび8つのクロスピークが観測されるが、小型分光計では、2つおよび3つのクロスピークのみが見える。
  FT-NIRおよび小型NIR分光器で測定されたPA11の可変温度NIRスペクトルの2D相関およびPCMW2D分析結果から、温度間隔30-190℃での脂肪族鎖構造における同期および変化の順序ならびにこのポリマーの水素結合に関する情報が得られた。40℃~110℃間の脂肪族鎖セグメントの再結晶の後、190℃での融点までに配座の乱れが起こる。スペクトル分解能の大きな差にもかかわらず、FT-NIRおよび小型NIR分光器で測定された可変温度スペクトルから調査されたポリマーの構造変化に関して同じような結論を引き出すことができることが判明した。


2017年1月11日水曜日

<スペクトルあれこれ>(12)マハラノビス距離とユークリッド距離


スペクトルを解析する際にマハラノビス距離をしばしば使用することがある。使用している割には、基本的な理解が欠けているかなと思いと、あらためてマハラノビス距離について見直してみた。2点間の距離を測定する際、一般的に使用するのはユークリッド距離だ。一方マハラノビス距離は、統計的な観点から距離をより一般化したものといってよい。
一次元の簡単な例を挙げる。仮に異常/正常を検査する簡易検査法があり、正常の標準値は2.5で異常の標準値は1であるとする。検査結果が1.8の場合は、ユークリッド距離を用いて、正常/異常の標準値の中間を境界(1.75)とすると、正常とみなされる。50回の検査結果について、他の方法で正常/異常を正確に調べたところ図Ar12-1となったとする。この結果からは正常/異常の境界は指標2あたりが妥当であり、1.8の検査結果は異常に分類する方が妥当であることがわかる。


これをマハラノビス距離で表したのが図Ar12-220回の正常な検査結果(指標2.46)の平均値からのマハラノビス距離を表している。異常と正常がきれいに別れ、正常/異常が明確に判別可能である。指標2はマハラノビス距離で7、1.8はマハラノビス距離で10となり異常に分類するのが妥当であることが確認できる。この場合のマハラノビス距離(MD)の計算式は12AR-1式となる。 



次にスペクトルのマハラノビス距離を考える。図Ar12-3と図Ar12-4に疑似吸光スペクトル例をしめす。図Ar12-320の疑似スペクトルの重ね書きである。図Ar12-4には図Ar12-3で示した20スペクトルの平均スペクトルと波長シフトと吸光度変化をした5スペクトルをしめした。

波長の10,15μmの吸光度を2次元のグラフにして図Ar12-5に示した。図Ar12-5では、重なりを防ぐためS2, S4の矢印を省略したが、ユークリッド距離は平均値(重心、黒丸●)からの距離を示す。一方マハラノビス距離は少し複雑な計算になる。12AR-1式は、個々のデータが独立の場合に使用できるが、スペクトルデータは一般に波長間のデータが相関を持つために12AR-1式は使用できない。多変量データに使用できる一般的なマハラノビス距離(マハラノビス汎距離)は12AR-2式で表せる。



Ar12-4のデータを使用して12AR-2式から得たマハラノビス距離とユークリッド距離を図AR12-6に示す。S1S5のスペクトルのマハラノビス距離が大きく拡大されていることが判る。この性質によりマハラノビス距離はスペクトル多変量解析においてアウトライヤー(外れ値)検知に使用される。多変量解析においてはスペクトルから直接マハラノビス距離を計算することは少なく、主成分分析などを行いその結果(スコア値)を使ってマハラノビス距離を計算しアウトライヤー検知を行うことが多い。



 


































……12AR-1