(57)【要約】【課題】多試料の分析に適した発酵産物原料の乳酸緩衝能推定方法を提供すること【解決手段】乳酸緩衝能既知の発酵産物原料に対して蛍光Ｘ線分析を行うことで得られる蛍光Ｘ線スペクトルに基づくデータと前記既知の乳酸緩衝能との相関から検量線を作成し、乳酸緩衝能未知の発酵産物原料に対して蛍光Ｘ線分析を行うことで得られる蛍光Ｘ線スペクトルに基づくデータを検量線に当てはめて乳酸緩衝能を推定する発酵産物原料の乳酸緩衝能推定方法である。  これにより、検量線を用いることで乳酸緩衝能未知の発酵産物原料に関して実際に乳酸緩衝能を測定する労力が不要となる。また、蛍光Ｘ線分析によれば前処理を含む分析時間が大きく短縮され、多試料の分析に適している。【選択図】図１
2016176830.tif 000002
【特許請求の範囲】
【請求項１】
乳酸緩衝能が既知である発酵産物原料から放出される光の分光スペクトル及びイオンの少なくとも一方を測定して得られるデータと前記既知の乳酸緩衝能との相関から検量線を作成する検量線作成工程と、
乳酸緩衝能が未知である発酵産物原料から放出される前記分光スペクトル及びイオンの少なくとも一方を測定して得られるデータを前記検量線に当てはめて乳酸緩衝能を推定する乳酸緩衝能推定工程と、
を有することを特徴とする発酵産物原料の乳酸緩衝能推定方法。
【請求項２】
前記分光スペクトル及びイオンの少なくとも一方が、前記乳酸緩衝能が既知である発酵産物原料に対して比色分析、原子吸光分析、ＩＣＰ発光分析、ＩＣＰ質量分析及び蛍光Ｘ線分析から選択される少なくとも１つの分析を行って得られるものであり、
乳酸緩衝能が既知である発酵産物原料から放出される光の分光スペクトル及びイオンの少なくとも一方を測定して得られる前記データが、該分光スペクトル及びイオンの少なくとも一方により特定される前記発酵産物原料中の元素の割合であることを特徴とする請求項１に記載の発酵産物原料の乳酸緩衝能推定方法。
【請求項３】
前記分光スペクトル及びイオンの少なくとも一方が、前記乳酸緩衝能が既知である発酵産物原料に対して蛍光Ｘ線分析を行って得られる蛍光スペクトルであることを特徴とする請求項２に記載の発酵産物原料の乳酸緩衝能推定方法。
【請求項４】
前記元素の割合が、Ｋ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ及びＳから選択される少なくとも２種以上の元素の割合であることを特徴とする請求項２又は３に記載の発酵産物原料の乳酸緩衝能推定方法。
【請求項５】
前記分光スペクトル及びイオンの少なくとも一方が、前記乳酸緩衝能が既知である発酵産物原料に対して近赤外分析を行って得られる近赤外スペクトルであり、
前記データは、前記近赤外スペクトルに対して２次微分処理を行って得られるものであることを特徴とする請求項１に記載の発酵産物原料の乳酸緩衝能推定方法。
【請求項６】
前記検量線が、重回帰分析を用いて作成されることを特徴とする請求項１〜５の何れか１項に記載の発酵産物原料の乳酸緩衝能推定方法。
【請求項７】
前記検量線が、ＰＬＳ回帰分析を用いて作成されることを特徴とする請求項１〜５の何れか１項に記載の発酵産物原料の乳酸緩衝能推定方法。
【請求項８】
前記発酵産物原料がサイレージ原料であることを特徴とする請求項１〜７の何れか１項に記載の発酵産物原料の乳酸緩衝能推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、発酵産物原料の乳酸緩衝能推定方法に関し、特に、乳酸緩衝能が未知の発酵産物原料に対する乳酸緩衝能の実測を不要とするための発酵産物原料の乳酸緩衝能推定方法に関する。
【背景技術】
【０００２】
発酵とは、微生物の働きで有機物が分解され、特定の物質を生成する現象をいい、これにより、サイレージ、漬物、チーズ、酒、醤油等の発酵産物が得られる。これら発酵産物によれば、素材の保存性が向上し、物により栄養価の向上し、有機物が分解されることでヒトを含む動物にとって消化しやすくなる等、様々な利点が得られる。ところで、良質な発酵産物を製造するためには、発酵産物原料のｐＨを速やかに低下させて雑菌の繁殖を抑制して目的となる発酵を進めることが重要である。例えば、良質なサイレージ（発酵産物）を調整するには、青刈り作物や生の牧草（発酵産物原料）をサイロ内で速やかに乳酸発酵させ、原料となる牧草のｐＨを速やかに低下させることが重要である。
【０００３】
しかし、原料となる牧草には緩衝力が存在し、牧草によってｐＨを低下させるのに必要な酸の量が異なる。すなわち、牧草の緩衝力が小さいと速やかにｐＨを低下させて発酵を進めることができ、この緩衝力は良質なサイレージを調整するために重要な要因となっている。
【０００４】
非特許文献１には、従来の緩衝力を測定する方法が開示されている。具体的には、牧草の乾燥物の抽出液を作成し、それらに０．１Ｎ乳酸液を添加してｐＨ４となる値から、牧草の乾燥物１ｇをｐＨ４とするために必要な乳酸量（以下、乳酸緩衝能と記す）を測定する方法が開示されている。
【０００５】
この方法によれば、牧草（発酵産物原料）の乳酸緩衝能を比較することで、良質なサイレージ、すなわち、発酵産物をその牧草から得ることができるか否かを事前に知ることができる。
【０００６】
なお、牧草の乾燥物に異なる濃度の乳酸溶液を同量ずつ加えて各区分のｐＨを測定して得られる緩衝曲線（縦軸をｐＨ、横軸を添加した０．１Ｎ乳酸液濃度で示す）から乳酸緩衝能を得ることも可能である。
【０００７】
非特許文献１においては、牧草に存在する緩衝力は、有機酸とたんぱく質に起因すると考察されている。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】Ｊ．Ｓｃｉ．Ｆｏｏｄ  Ａｇｒｉｃ．， １９６２，Ｖｏｌ．１３，Ｊｕｌｙ ｐ．３９５−４００
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかし、乳酸緩衝能の測定は、１個の牧草試料について複数区分の抽出液を作成したうえでｐＨ測定を行って緩衝曲線を作成するか、あるいは実際に滴定作業を行う必要があり、非常に労力を伴う。したがって、多試料の分析への対応が困難である。
【００１０】
また、緩衝力に影響を与える有機酸やたんぱく質を測定することにより、牧草の乳酸緩衝能を推定することも考えられる。しかし、有機酸を測定するにはガスクロマトグラフィ（ＧＣ）やＨＰＬＣ等が必要となり１試料あたりの分析時間は必然的に長くなり、多試料への対応は困難である。また、たんぱく質を測定するにはケルダール法等を行わなければならず、試料の前準備に非常に労力を必要とし、且つ測定時間も非常に長くなり、多試料への対応は困難である。尚、有機酸やたんぱく質を測定することにより乳酸緩衝能が推定可能か否かを検証した知見はない。
【００１１】
本発明は、上記課題に鑑みてなされたものであり、その目的は、多試料の分析に適した発酵産物原料の乳酸緩衝能推定方法を提供することにある。
【課題を解決するための手段】
【００１２】
上記目的を達成するための請求項１に記載の発酵産物試料の乳酸緩衝能推定方法は、乳酸緩衝能が既知である発酵産物原料から放出される光の分光スペクトル及びイオンの少なくとも一方を測定して得られるデータと前記既知の乳酸緩衝能との相関から検量線を作成する検量線作成工程と、乳酸緩衝能が未知である発酵産物原料から放出される前記分光スペクトル及びイオンの少なくとも一方を測定して得られるデータを前記検量線に当てはめて乳酸緩衝能を推定する乳酸緩衝能推定工程と、を有することを特徴とする。
【００１３】
この構成によれば、乳酸緩衝能既知の発酵産物原料から放出される光の分光スペクトル等を測定して得られるデータと既知の乳酸緩衝能との相関から検量線が作成されるので、この検量線を用いることで、乳酸緩衝能が未知の発酵産物原料に対して実際に乳酸緩衝能の測定を行う労力が不要となる。
【００１４】
また、光の分光スペクトル等の測定は、発酵産物原料から放出される光を分光し、あるいは放出されるイオンを検出することでなされるので、発酵産物原料に対する前処理を含む分析時間が大きく短縮され、多試料の分析に適している。
【００１５】
本発明に係る好ましい態様は以下の通りである。
【００１６】
（１）分光スペクトル及びイオンの少なくとも一方が、前記乳酸緩衝能が既知である発酵産物原料に対して比色分析、原子吸光分析、ＩＣＰ発光分析、ＩＣＰ質量分析及び蛍光Ｘ線分析から選択される少なくとも１つの分析を行って得られるものであり、
乳酸緩衝能が既知である発酵産物原料から放出される光の分光スペクトル及びイオンの少なくとも一方を測定して得られる前記データが、該分光スペクトル及びイオンの少なくとも一方により特定される前記発酵産物原料中の元素の割合である。これにより、迅速に測定可能な発酵産物原料中の元素の割合に基づき、迅速な検量線の作成及び乳酸緩衝能の推定が可能となる。さらに、元素によっては他の用途でも別途測定され得る元素もあることから、発酵産物原料について既に測定され、データとして蓄積された元素の割合の値を用いることで新たに測定することなく検量線を得ることが可能となる。
【００１７】
（２）分光スペクトル及びイオンの少なくとも一方が、前記乳酸緩衝能が既知である発酵産物原料に対して蛍光Ｘ線分析を行って得られる蛍光スペクトルである。これにより、発酵産物原料の前処理を含む分析がより迅速なものとなり、さらに迅速な検量線の作成及び乳酸緩衝能の推定が可能となる。
【００１８】
（３）元素の割合が、Ｋ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ及びＳから選択される少なくとも２種以上の元素の割合である。これにより、より精度の高い検量線の作成及び乳酸緩衝能の推定が可能となる。
【００１９】
（４）分光スペクトル及びイオンの少なくとも一方が、前記乳酸緩衝能が既知である発酵産物原料に対して近赤外分析を行って得られる近赤外スペクトルであり、
前記データは、前記近赤外スペクトルに対して２次微分処理を行って得られるものである。これにより、迅速に測定可能な近赤外スペクトルに基づき、迅速な検量線の作成及び乳酸緩衝能の推定が可能になるとともに、２次微分処理によりデータの精度が向上し、より精度の高い検量線の作成及び乳酸緩衝能の推定が可能となる。また、発酵産物原料の光学的な吸収情報を網羅的に捉えることで検量線を作成することができ、発酵産物原料中の特定の成分に着目してこれを対比する必要もない。
【００２０】
（５）検量線が、重回帰分析を用いて作成されるものである。
【００２１】
（６）検量線が、ＰＬＳ回帰分析を用いて作成されるものである。
【００２２】
（７）発酵産物原料が、サイレージ原料である。
【発明の効果】
【００２３】
本発明によれば、乳酸緩衝能が未知の発酵産物原料に対して実際に乳酸緩衝能の測定を行う労力が不要となる。また、発酵産物原料に対する前処理を含む分析時間が大きく短縮される。したがって、多試料の分析に適した方法であるとともに、光学的な分析が行われることから使用する試薬や廃液の量も低減するか、あるいは無くなり、環境への負荷も小さいものとなる。
【図面の簡単な説明】
【００２４】
【図１】本実施の形態にかかる発酵産物原料の乳酸緩衝能推定方法の概要を示すフローチャートである。
【図２】発酵産物原料に対する蛍光Ｘ線分析により得られる蛍光Ｘ線スペクトルを示す図である。
【発明を実施するための形態】
【００２５】
次に、本発明の実施の形態について図に基づいて詳細に説明する。本発明の発酵産物原料の乳酸緩衝能推定方法を、図１及び図２を参照して説明する。図１は、本発明の実施の形態に係る発酵産物原料の乳酸緩衝能推定方法の概要を示すフローチャートであり、図２は発酵産物原料に対する蛍光Ｘ線分析により得られる蛍光Ｘ線スペクトルを示す図である。
【００２６】
図１に示すように、本実施の形態に係る発酵産物原料の乳酸緩衝能推定方法は、乳酸緩衝能が既知である発酵産物原料の測定工程（Ｓ１０１）、この測定により得られるデータと上記既知の乳酸緩衝能との相関から検量線を作成する検量線作成工程（Ｓ１０２）、検量線の精度を評価する検量線評価工程（Ｓ１０３）、乳酸緩衝能が未知である発酵産物原料の測定工程（Ｓ１０４）、及び乳酸緩衝能を推定する乳酸緩衝能推定工程（Ｓ１０５）からなる。
【００２７】
本実施の形態においては、発酵産物はサイレージであり、発酵産物原料はサイレージ原料である。サイレージ原料は、乳酸発酵によりサイレージとなる原料であり、例えば、干草、牧草、とうもろこし等が挙げられる。尚、サイレージとは、青刈り作物や生の牧草をサイロ内で乳酸発酵させて貯蔵した飼料をいう（大辞林、第三版）。
【００２８】
サイレージ原料（発酵産物原料）の測定は、サイレージ原料の前処理を含む測定が迅速である測定であればどのようなものであってもよいが、具体的には、サイレージ原料に対して比色分析、原子吸光分析、ＩＣＰ発光分析、蛍光Ｘ線分析、ＩＣＰ質量分析及び近赤外分析が挙げられる。
【００２９】
したがって、サイレージ原料（発酵産物原料）から放出される光の分光スペクトルとは、サイレージ原料に対して比色分析、原子吸光分析、ＩＣＰ発光分析、蛍光Ｘ線分析及び近赤外分析から選択される少なくとも１つの分析を行って得られる分光スペクトルであり、また、サイレージ原料から放出されるイオンとは、サイレージ原料に対してＩＣＰ質量分析を行って得られるイオンである。
【００３０】
本実施の形態においては、乳酸緩衝能が既知であるサイレージ原料に対して蛍光Ｘ線分析を行うことによりサイレージ原料から放出される蛍光Ｘ線スペクトルを得る。
【００３１】
［乳酸緩衝能が既知である発酵産物原料の測定工程（Ｓ１０１）］
まず、ステップＳ１０１では、上述のとおり、乳酸緩衝能が既知であるサイレージ原料に対して蛍光Ｘ線分析を行い、サイレージ原料から放出される蛍光Ｘ線スペクトルを得る。蛍光Ｘ線分析は、市場で入手できる蛍光Ｘ線分析装置を用いて行うことができる。
【００３２】
蛍光Ｘ線分析装置では、実際に測定して得られた蛍光Ｘ線スペクトルに対してファンダメンタルパラメータ法（ＦＰ法）を用いることにより変換したデータを得る。
【００３３】
ＦＰ法とは、試料を構成する元素の種類と組成濃度が分かれば、試料から放出される各元素の蛍光Ｘ線の理論強度を計算することができるということを利用して、未知試料を測定して得られた各元素の蛍光Ｘ線強度に一致するような組成を推定する方法である。
【００３４】
一般に、蛍光Ｘ線分析装置には、予め純粋な元素、又は純粋な元素を準備できない元素ではその元素を含む組成が既知の化合物から蛍光Ｘ線を測定し、理論強度を測定強度に変換する装置係数が元素ごとに求められており、データライブラリとして登録されている。したがって、ＦＰ法によれば、実際に測定して得られた蛍光Ｘ線スペクトルから、標準試料を準備することなくサイレージ原料中の元素全体の中での特定の元素の質量割合（ｗｔ％ＤＭ）を推定することができる。
【００３５】
具体的には、図２に示すように、実際に測定して得られた蛍光Ｘ線スペクトルは、エネルギー順位（ｋｅＶ）を横軸とし、蛍光Ｘ線強度（ＣＰＳ）を縦軸とするグラフとして表される。
【００３６】
したがって、あるサイレージ原料ｉに対して蛍光Ｘ線分析を行った場合に、特定のエネルギー順位の値（ｘ）に対応する蛍光Ｘ線強度の値をＦＰ法により換算した値をサイレージ原料の特定の元素の質量割合Ｅｉ（ｗｔ％ＤＭ）として表現することができ、かかる質量割合はｘを変数として定まる関数Ｅｉ（ｘ）と表現することができる。尚、ｗｔ％ＤＭのＤＭとは、サイレージ原料の乾物（Ｄｒｙ Ｍａｔｔｅｒ）を意味する。
【００３７】
上記サイレージ原料に対する蛍光Ｘ線分析は、検量線作成用に収集されたｐ点の乳酸緩衝能が既知であるサイレージ原料全てに対して行う。
【００３８】
［検量線作成工程（Ｓ１０２）］
次に、検量線作成工程（Ｓ１０２）について説明する。本工程では、ステップＳ１０１において得られた、サイレージ原料中の試料ｉの特定の元素の質量割合Ｅｉ（ｘ）とサイレージ原料の既知の乳酸緩衝能ＬＢとから検量線（乳酸緩衝能の推定式）を作成する。
【００３９】
具体的には、サイレージ原料ｉについて複数（ｎ個）のエネルギー順位ｘｊ（０≦ｘｊ≦３０）を選択する（但し、ｊは１≦ｊ≦ｎを充足する数でｊ番目のエネルギー順位を意味する）。そして、エネルギー順位ｘｊが０〜３０の範囲の元素の質量の総和に対する選択されたエネルギー順位ｘｊにおける特定の元素の質量割合Ｅｉ（ｘｊ）と、係数Ａｊと、を用いて以下の式に示す線形結合
【数１】
2016176830.tif 000003
を得る。すなわち、検量線（推定式）は回帰分析に用いる回帰式である。なお、エネルギー順位の値としては、図２に示す蛍光Ｘ線強度のピークが存在する位置の値を選択する。これにより、Ｅｉ（ｘ）は特定の元素の質量割合を示すこととなる。また、エネルギー順位ｘの値としては、０〜３０の範囲で選択することが可能であるが、例えば、Ｅｉ（ｘ）がＫ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ、Ｓ、Ｎａ、Ａｌ、Ｍｎ、Ｆｅ、Ｚｎ、Ｓｉ、Ｃｕの質量割合を示すようなｘの値が挙げられ、好ましくは、Ｅｉ（ｘ）がＫ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ、Ｓの質量割合を示すようなｘの値が選択される。
【００４０】
さらに、選択するエネルギー順位の個数ｎについて特に制限はないが、後述する回帰式の評価精度の観点からｎは２以上であることが好ましく、また、オーバーフッティングの観点から多すぎない個数をとる必要がある。特に好ましくは、ｎ＝２〜６である。
【００４１】
そして、サイレージ原料ｉの既知の乳酸緩衝能ＬＢをＬＢｉで表せば、以下の式で表される残差
【数２】
2016176830.tif 000004
が最小となるように最小二乗法を実行し、Ａｊを決定する。尚、残差を最小にする方法は、最小二乗法に限られず、他の種々の数学的アルゴリズムを用いることが可能である。これにより、サイレージ原料の乳酸緩衝能ＬＢｉが、数１で表されるサイレージ原料ｉの特定の元素の質量割合Ｅｉ（ｘｊ）の線形結合で近似されるように係数Ａｊが決定される。
【００４２】
そして、決定された係数Ａｊを数１の式に代入し、サイレージ原料の乳酸緩衝能の推定値＜ＬＢｉ＞と表し、以下の式で表される回帰式
【数３】
2016176830.tif 000005
を定める。すなわち、当該回帰式は、サイレージ原料ｉの特定の元素の質量割合Ｅｉ（ｘｊ）が説明変数であり、測定対象のサイレージ原料の乳酸緩衝能推定値＜ＬＢｉ＞が目的変数である。本実施の形態では、この数３の式が、乳酸緩衝能が未知であるサイレージ原料の蛍光Ｘ線強度からサイレージ原料の乳酸緩衝能を推定する検量線（推定式）として用いられる（以上、検量線作成工程Ｓ１０２）。
【００４３】
［検量線の精度評価工程（Ｓ１０３）］
次に、作成した検量線の精度評価を行う。この精度評価には、以下の各数値を算出する。
【００４４】
本工程では、数３の式の検量線（推定式）に乳酸緩衝能が既知のｐ点のサイレージ原料のスペクトルデータ（ＦＰ法により変換後のもの）を代入し、検量線による理論的な乳酸緩衝能推定値＜ＬＢｉ＞を求める。この乳酸緩衝能推定値＜ＬＢｉ＞と上記既知のＬＢｉとの重相関係数の二乗Ｒ２（以下、検量線の相関係数という）、及び検量線の標準誤差ＳＥＣ（ｓｔａｎｄａｒｄ ｅｒｒｏｒ ｏｆ ｃａｌｉｂｒａｔｉｏｎ）を算出する。
【００４５】
検量線の相関係数Ｒ２は、以下の式により算出される。
【００４６】
【数４】
2016176830.tif 000006
【００４７】
ただし、
【数５】
2016176830.tif 000007
は、ＬＢｉの標本平均を示している。
【００４８】
一方、標準偏差ＳＥＣは、以下の式により算出される。
【００４９】
【数６】
2016176830.tif 000008
【００５０】
さらに、公知の方法で乳酸緩衝能がＬＢｍと測定されている検量線評価用のサイレージ原料ｍのｎ個のエネルギー順位に対する蛍光Ｘ線強度を測定してＦＰ法により特定の元素の質量割合に変換し、係数Ａｊが定まった数３の式に代入することで乳酸緩衝能の推定値を算出する。ここで、この算出される乳酸緩衝能の推定値を［ＬＢｍ］とする。なお、検量線評価用のサイレージ原料ｍはｑ点存在する。
【００５１】
算出された乳酸緩衝能の推定値［ＬＢｍ］と、サイレージ原料ｍの乳酸緩衝能ＬＢｍとの間の重相関係数ｒの二乗ｒ２（以下、評価用試料の相関係数という）、及び予測標準誤差ＳＥＰ（ｓｔａｎｄａｒｄ ｅｒｒｏｒ ｏｆ ｐｒｅｄｉｃｔｉｏｎ）を算出する。評価用試料の相関係数ｒ２の算出は、本工程における検量線作成用試料の相関係数Ｒ２の算出とほぼ同様の方法で行うことが可能である。
【００５２】
一方で、ＳＥＰは、以下の式で算出される。
【００５３】
【数７】
2016176830.tif 000009
【００５４】
式中、ｄｍは乳酸緩衝能の推定値［ＬＢｍ］と既知のＬＢｍとの差であり、
【数８】
2016176830.tif 000010
は、ｄｍの標本平均、すなわち、以下の式
【数９】
2016176830.tif 000011
で算出される値である。
【００５５】
また、検量線を評価するため、ＥＩ（Ｅｖａｌｕａｔｉｏｎ  Ｉｎｄｅｘ）を用いる。ＥＩは、作成した検量線の推定精度の指標であり、以下の式
ＥＩ＝１００×２×ＳＤＰ／Ｒａｎｇｅ
により算出される。
【００５６】
式中、ＳＤＰ（ｓｔａｎｄａｒｄ  ｄｅｖｉａｔｉｏｎ  ｏｆ  ｐｒｅｄｉｃｔｉｏｎ）は以下の式
【数１０】
2016176830.tif 000012
で表わされる、ｑ点の検量線評価用のサイレージ原料ｍでの誤差の標準偏差であり、Ｒａｎｇｅはｑ点の検量線評価用のサイレージ原料ｍの実際の乳酸緩衝能の上限から下限までの範囲である。すなわち、Ｒａｎｇｅは、ｑ点の検量線評価用のサイレージ原料ｍの乳酸緩衝能ＬＢｍの最大値と最小値の差である。例えば、最大値が６０（ｍｇ／ｇＤＭ）、最小値が２０（ｍｇ／ｇＤＭ）であれば、Ｒａｎｇｅは４０（ｍｇ／ｇＤＭ）となる。
【００５７】
ＥＩは、その値により、
１２．５未満（Ａ：非常に精度が高い）、
１２．５〜２５．０未満（Ｂ：精度が高い）、
２５．０〜３７．５未満（Ｃ：精度が低い）、
３７．５〜５０．０未満（Ｄ：非常に精度が低い）、及び
５０．０以上（Ｅ：精度の推定不能）
にランク分けされる。
【００５８】
検量線の相関係数Ｒ２及び評価用試料の相関係数ｒ２は、ともに値が高ければ検量線（推定式）の精度が高いと判断され、検量線の標準誤差ＳＥＣ及び予測標準誤差ＳＥＰは、ともに値が小さいほど検量線（推定式）の精度が高いと判断される。また、ＥＩについては、上述のとおりである（以上、検量線の精度評価工程Ｓ１０３）。
【００５９】
［乳酸緩衝能が未知である発酵産物原料の測定工程（Ｓ１０４）］
ステップＳ１０３において作製した検量線（推定式）の精度をステップＳ１０４において確認した後、その検量線を用いて乳酸緩衝能が未知であるサイレージ原料に対して蛍光Ｘ線分析を行い、ＦＰ法によりサイレージ原料の特定の元素の質量割合を求める。求める特定の元素の質量割合は、ステップＳ１０２における検量線の作成に用いた元素と同じ元素のものである。
【００６０】
なお、本工程において確認した検量線（推定式）の精度は、次のステップＳ１０５において推定される乳酸緩衝能の値の信頼性を判断するための一つの指標となる（以上、乳酸緩衝能が未知であるサイレージ原料の測定工程Ｓ１０４）。
【００６１】
［乳酸緩衝能推定工程（Ｓ１０５）］
最後に、ステップＳ１０４において求められた、乳酸緩衝能が未知のサイレージ原料の特定の元素の質量割合をステップＳ１０３で作成した検量線（推定式）に代入し、得られた値を乳酸緩衝能の推定値として得る（以上、乳酸緩衝能推定工程Ｓ１０５）。
【００６２】
したがって、本実施の形態に係る発酵産物原料の乳酸緩衝能の推定方法によれば、乳酸緩衝能既知のサイレージ原料に対して蛍光Ｘ線分析を行って得られた蛍光Ｘ線スペクトルをＦＰ法により変換した特定の元素の割合と既知の乳酸緩衝能との相関から検量線が作成されるので、この検量線を用いることで、乳酸緩衝能が未知のサイレージ原料（発酵産物原料）に対して実際に乳酸緩衝能の測定を行う労力が不要となる。
【００６３】
また、蛍光Ｘ線分析は、サイレージ原料の風乾物試料に対してほぼそのまま一次Ｘ線を照射し、試料から放出される試料の元素特有の蛍光Ｘ線を検出することによりなされるので、サイレージ原料の測定に要する時間が大きく短縮され、多試料の分析に適している。特に、検量線の作成及びその評価には多くの試料を測定することが要求されるところ、本方法によれば、多試料を迅速に測定可能であるから、迅速な検量線の作成とその評価、及び乳酸緩衝能未知の試料について乳酸緩衝能の迅速な推定が可能となる。
【００６４】
さらに、サイレージ原料の所定の元素は、乳酸緩衝能の推定の目的以外にも分析され得る。したがって、例えば、サイレージ原料について既に測定され、データとして蓄積された所定の元素の割合の値を用いることで、新たに測定を行うことなく乳酸緩衝能を測定するための検量線（推定式）を作成することができ、乳酸緩衝能が未知のサイレージ原料についてより迅速な乳酸緩衝能の推定を行うことが可能となる。
【００６５】
なお、上記ステップＳ１０２では、回帰式（検量線）を作成する手法として重回帰分析を用いる方法を説明したが、これに限られず、回帰式を定めることができる分析手法であれば、例えば、ＰＬＳ回帰分析、主成分分析やクラスター分析等の他の種々のものを用いても良い。例えば、ＰＬＳ回帰分析を用いる場合は、サイレージ原料ｉの蛍光Ｘ線強度をＦＰ法により変換した特定の元素の割合に基づいて定められる潜在変数を設定し、これを説明変数として回帰式が設定される。このＰＬＳ回帰分析は、回帰式を定めるために選択される元素数が多い場合（すなわち、選択されるエネルギー順位ｘの個数が多い場合）、或いは選択される元素間で元素の割合のデータに強い相関が生じている場合に、オーバーフッティングや多重共線性を防止するために有効な分析手法である。
【００６６】
また、ステップＳ１０１における測定対象は特定の元素の割合に限られない。例えば、サイレージ原料に対して近赤外分析を行う場合には、単に反射光スペクトルのデータに基づいて（すなわち、サイレージ中の特定の成分の割合を比較することなく）検量線（推定式）を作成し、乳酸緩衝能を推定することが可能である。
【００６７】
以下、サイレージ原料に対して近赤外分析を行う場合を例に、上記実施の形態の変形例を、上記実施の形態とは異なる部分について主に説明する。なお、上記実施の形態に対応するステップについては符号（´）を付けて示す。
【００６８】
［乳酸緩衝能が既知であるサイレージ原料の測定工程（Ｓ１０１´）］
まず、ステップＳ１０１では、乳酸緩衝能が既知であるサイレージ原料に対して近赤外分析を行い、サイレージ原料から反射される反射光スペクトルを得る。反射光スペクトルを採取する波長範囲は、例えば、２００〜２５００ｎｍであり、好ましくは４００ｎｍ〜２５００ｎｍである。近赤外分析は、市場で入手できる近赤外光分光光度計を用いて行うことができる。
【００６９】
近赤外分析を用いた測定は、乳酸緩衝能が既知であるサイレージ原料ｉに対して４００〜２５００ｎｍの波長の光を照射し、その反射光として近赤外光分光光度計の検出器で強度スペクトルデータＩｉ（λ）を得る。
【００７０】
得られた強度スペクトルデータＩｉ（λ）に対して、ピーク位置を検出するため２次微分処理が行われる。
【００７１】
具体的に、２次微分処理は、以下の式
【数１１】
2016176830.tif 000013
を計算することにより行われる。なお、ピーク位置の検出に１次微分処理を用いてもよい（以上、乳酸緩衝能が既知であるサイレージ原料の測定工程Ｓ１０１´）。
【００７２】
［検量線作成工程（Ｓ１０２´）］
次に、検量線作成工程（Ｓ１０２´）について説明する。得られた波長λから複数の波長（ｎ個）を選択し、ｊ番目の波長における２次微分処理した強度スペクトルデータＩｉ（λｊ）と、サイレージ原料の既知の乳酸緩衝能とから重回帰分析によって検量線（乳酸緩衝能の推定式）を求める。この式は上記数１の式からＥｉ（ｘｊ）をＩｉ（λｊ）に置き換えたものである。
【００７３】
選択される波長は、好ましくは１４００ｎｍ〜１５００ｎｍ、１７００ｎｍ〜１８００ｎｍ及び２１００ｎｍ〜２２００ｎｍの区間からそれぞれ選択される（尚、重回帰分析に代えてＰＬＳ回帰分析が行われる場合、選択される波長は好ましくは１１００ｎｍ〜２５００ｎｍの範囲であるか、あるいは１４５０ｎｍ〜１８２０ｎｍ及び１９００ｎｍ〜２２００ｎｍの範囲であることが好ましい）。
【００７４】
また、選択される波長の個数は、選択する波長の個数ｎについて複数であれば特に制限はないが、オーバーフッティングの観点から多すぎない個数をとる必要がある。好ましくは、ｎ＝３〜４の範囲である。
【００７５】
以降、上記実施の形態同様に本工程（Ｓ１０２´）の残りの部分を行い、さらに検量線の精度評価工程（Ｓ１０３´）、乳酸緩衝能が未知であるサイレージ原料の測定工程（Ｓ１０４´）及び乳酸緩衝能推定工程（Ｓ１０５´）を行い、近赤外分析による乳酸緩衝能の推定を迅速に行うことができる。
【００７６】
この変形例によれば、サイレージ原料中の特定の成分に着目してこれを対比することなく、サイレージ原料の光学的な吸収情報を網羅的に捉えることで検量線を作成することができ、したがって、迅速に乳酸緩衝能未知のサイレージ原料の乳酸緩衝能を推定することができる。
【００７７】
なお、本発明は上記実施の形態に限定されることはなく、発明の趣旨を逸脱しない範囲で種々変更可能である。
【００７８】
例えば、上記実施の形態において、検量線の精度評価工程（Ｓ１０３）が存在するが、当該ステップＳ１０３は必須というわけではない。
【００７９】
検量線の精度評価工程が無いことにより検量線の精度に差は生じるとしても、迅速な分析手法を用いることにより乳酸緩衝能を実際に測定することなく推定することは可能であり、したがって、多試料の分析に適したサイレージ原料の乳酸緩衝能推定方法を提供するという本願発明の目的を達成することが可能である。
【００８０】
また、上記実施の形態においては発酵産物原料はサイレージ原料（すなわち、干草、牧草、とうもろこし等）としているが、これに限られるものではない。例えば、漬物の発酵産物原料となる野菜（ナス、大根、キュウリ、キャベツ、ニンジン等）、チーズの発酵産物原料となる乳（牛、羊、山羊等の乳）、酒の発酵産物原料となる穀物や果実、醤油の発酵産物原料となる大豆等、ヒトを含む動物が摂取可能な種々の発酵産物を得るための発酵産物原料が挙げられる。
【実施例】
【００８１】
以下、本発明を実施例により説明する。
【００８２】
（実施例１−１）
北海道内で収穫されたイネ科主体牧草１７９点の風乾粉砕試料を供試試料（サイレージ原料）とした。これら供試試料を１点につき２ｇに小分けした区分を複数区分作成し、これらの区分に対して０．０９％，０．１３％，０．１７％，０．２６％の乳酸溶液４０ｍｌをそれぞれ添加した（一部の供給試料については、さらに０．３４％，０．５１％の乳酸溶液４０ｍｌを添加した区分を設けた）。その後、各区分を３分間振盪した後に各区分のｐＨを測定した。
【００８３】
次に、そのｐＨ測定結果を縦軸とし、使用した乳酸溶液の濃度を横軸として緩衝曲線を作成し、この緩衝曲線から供給試料１ｇのｐＨを４．０とするために必要な乳酸量を測定した。その後、風乾粉砕試料中に含まれる水分を考慮して、供試試料の乾物１ｇのｐＨを４．０とするために必要な乳酸量（乳酸緩衝能）に換算した。
【００８４】
得られた乳酸緩衝能の範囲は、１７９点の供試試料において１２．６〜８３．７ｍｇ／ＤＭ（供試試料の乾物１ｇ）であった。
【００８５】
次に、供試試料３ｇを直径３０ｍｍのタブレット状に圧縮形成し、当該タブレットに対して蛍光Ｘ線分析を行った。具体的には、エネルギー分散型蛍光Ｘ線分析装置（品番ＪＳＸ−３１００ＲＩＩ、日本電子株式会社製）で蛍光Ｘ線スペクトルを測定し、得られた値をＦＰ法により換算することで、測定結果から得られた全元素の質量に対するＫ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ及びＳの元素の質量割合（ｗｔ％ＤＭ）を求めた。
【００８６】
なお、蛍光Ｘ線スペクトルの測定は、真空条件下（３０Ｐａ以下）、管電圧３０ｋＶでＸ線係数率２５，０００ＣＰＳとなるように管電流を設定してタブレット状とした供試試料にＸ線を照射し、０〜３０ｋｅＶのエネルギー範囲の蛍光Ｘ線スペクトルを測定した。
【００８７】
ランダムに検量線作成用の試料１４０点と検量線の精度評価用の試料３９点に供試試料を区分し、検量線作成用試料のＫ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ及びＳの元素の質量割合（ｗｔ％ＤＭ）の測定結果から重回帰分析で検量線（すなわち、乳酸緩衝能の推定式）を決定した。次に、同様に検量線の精度評価用試料３９点のＫ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ及びＳの元素の質量割合（ｗｔ％ＤＭ）を上記決定した検量線（推定式）に当てはめることで検量線の精度評価用試料３９点の乳酸緩衝能を求めた。
【００８８】
検量線（推定式）の精度評価は、検量線の相関係数Ｒ２および標準誤差ＳＥＣを求めることにより、並びに、検量線の精度評価用試料３９点の元素（Ｋ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ及びＳ）の質量割合（ｗｔ％ＤＭ）を上記検量線に当てはめて推定された評価用試料の相関係数ｒ２、標準誤差ＳＥＰ及びＥＩを求めることにより、行った。
【００８９】
決定した乳酸緩衝能（推定値）［ＬＢ］を求めるための検量線（推定式）を以下に示し、この検量線（推定式）の精度評価の結果を表１に示す。
【００９０】
［ＬＢ］＝１３．７６×Ｋ（ｗｔ％ＤＭ）＋（−１９．７０）×Ｃｌ（ｗｔ％ＤＭ）＋４７．８３×Ｍｇ（ｗｔ％ＤＭ）＋４．７２×Ｃａ（ｗｔ％ＤＭ）＋１８．３８×Ｐ（ｗｔ％ＤＭ）＋（−１３．８８）×Ｓ（ｗｔ％ＤＭ）−１．６６
【００９１】
（実施例１−２）
実施例１−１においてランダムに区分された検量線作成用の試料１４０点のＫ、Ｃｌ、Ｍｇ、Ｃａ及びＳの元素の質量割合（ｗｔ％ＤＭ）の測定結果から重回帰分析で検量線（すなわち、乳酸緩衝能の推定式）を決定した。
【００９２】
同じく、実施例１−１においてランダムに区分された検量線の精度評価用の試料３９点のＫ、Ｃｌ、Ｍｇ、Ｃａ及びＳの元素の質量割合（ｗｔ％ＤＭ）を上記検量線（推定式）に当てはめることで検量線の精度評価用試料３９点の乳酸緩衝能を求めた。
【００９３】
検量線（推定式）の精度評価は、実施例１−１と同様に行った。
【００９４】
決定した検量線（推定式）を以下に示し、この検量線（推定式）の精度評価の結果を表１に示す。
【００９５】
［ＬＢ］＝１３．５７×Ｋ（ｗｔ％ＤＭ）＋（−１９．７０）×Ｃｌ（ｗｔ％ＤＭ）＋４６．２２×Ｍｇ（ｗｔ％ＤＭ）＋４．９４×Ｃａ（ｗｔ％ＤＭ）＋１１．０８×Ｓ（ｗｔ％ＤＭ）−１．０４
【００９６】
（実施例１−３）
実施例１−１においてランダムに区分された検量線作成用の試料１４０点のＫ、Ｃｌ、Ｍｇ及びＣａの元素の質量割合（ｗｔ％ＤＭ）の測定結果から重回帰分析で検量線（すなわち、乳酸緩衝能の推定式）を決定した。
【００９７】
同じく、実施例１−１においてランダムに区分された検量線の精度評価用の試料３９点のＫ、Ｃｌ、Ｍｇ及びＣａの元素の質量割合（ｗｔ％ＤＭ）を上記検量線（推定式）に当てはめることで検量線の精度評価用試料３９点の乳酸緩衝能を求めた。
【００９８】
検量線（推定式）の精度評価は、実施例１−１と同様に行った。
【００９９】
決定した検量線（推定式）を以下に示し、この検量線（推定式）の精度評価の結果を表１に示す。
【０１００】
［ＬＢ］＝１４．０２×Ｋ（ｗｔ％ＤＭ）＋（−１９．４８）×Ｃｌ（ｗｔ％ＤＭ）＋５１．８７×Ｍｇ（ｗｔ％ＤＭ）＋４．５９×Ｃａ（ｗｔ％ＤＭ）−０．０２
【０１０１】
（実施例１−４）
実施例１−１においてランダムに区分された検量線作成用の試料１４０点のＫ、Ｃｌ及びＭｇの元素の質量割合（ｗｔ％ＤＭ）の測定結果から重回帰分析で検量線（すなわち、乳酸緩衝能の推定式）を決定した。
【０１０２】
同じく、実施例１−１においてランダムに区分された検量線の精度評価用の試料３９点のＫ、Ｃｌ及びＭｇの元素の質量割合（ｗｔ％ＤＭ）を上記検量線（推定式）に当てはめることで検量線の精度評価用試料３９点の乳酸緩衝能を求めた。
【０１０３】
検量線（推定式）の精度評価は、実施例１−１と同様に行った。
【０１０４】
決定した検量線（推定式）を以下に示し、この検量線（推定式）の精度評価の結果を表１に示す。
【０１０５】
［ＬＢ］＝１４．２７×Ｋ（ｗｔ％ＤＭ）＋（−２０．４８）×Ｃｌ（ｗｔ％ＤＭ）＋５９．４９×Ｍｇ（ｗｔ％ＤＭ）＋０．３４
【０１０６】
（実施例１−５）
実施例１−１においてランダムに区分された検量線作成用の試料１４０点のＫ及びＣｌの元素の質量割合（ｗｔ％ＤＭ）の測定結果から重回帰分析で検量線（すなわち、乳酸緩衝能の推定式）を決定した。
【０１０７】
同じく、実施例１−１においてランダムに区分された検量線の精度評価用の試料３９点のＫ及びＣｌの元素の質量割合（％ＤＭ）を上記検量線（推定式）に当てはめることで検量線の精度評価用試料３９点の乳酸緩衝能を求めた。
【０１０８】
検量線（推定式）の精度評価は、実施例１−１と同様に行った。
【０１０９】
決定した検量線（推定式）を以下に示し、この検量線（推定式）の精度評価の結果を表１に示す。
【０１１０】
［ＬＢ］＝１２．８２×Ｋ（ｗｔ％ＤＭ）＋（−２１．９０）×Ｃｌ（ｗｔ％ＤＭ）＋１６．２８
【０１１１】
（比較例１）
実施例１−１においてランダムに区分された検量線作成用の試料１４０点のＫの元素の質量割合（ｗｔ％ＤＭ）の測定結果から重回帰分析で検量線（すなわち、乳酸緩衝能の推定式）を決定した。
【０１１２】
同じく、実施例１−１においてランダムに区分された検量線の精度評価用の試料３９点のＫの元素の質量割合（ｗｔ％ＤＭ）を上記検量線（推定式）に当てはめることで検量線の精度評価用試料３９点の乳酸緩衝能を求めた。
【０１１３】
検量線（推定式）の精度評価は、実施例１−１と同様に行った。
【０１１４】
決定した検量線（推定式）を以下に示し、この検量線（推定式）の精度評価の結果を表１に示す。
【０１１５】
［ＬＢ］＝１１．６１×Ｋ（ｗｔ％ＤＭ）＋４．６６
【０１１６】
【表１】
2016176830.tif 000014
【０１１７】
［結果］
表１に示すように、選択した６種の元素（Ｋ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ、Ｓ）のうち少なくとも２種以上の元素の質量割合（ｗｔ％ＤＭ）を用いた実施例１−１〜実施例１−５については、相関（Ｒ２）及び相関（ｒ２）が高い値となり、且つ、ＳＥＣ、ＳＥＰ及びＥＩが低い値に維持されているため、検量線（推定式）の精度が高いことが示された。
【０１１８】
一方、１種の元素（Ｋ）のみを選択した比較例１については、実施例１−１〜実施例１−５と比較して相関（Ｒ２）及び相関（ｒ２）が大きく低下し、ＳＥＣ、ＳＥＰ及びＥＩの値も大きく上昇しており、検量線（推定式）の精度が低いことが示された。
【０１１９】
実施例１−１〜実施例１−５間の対比においては、実施例１−５については、実施例１−１〜実施例１−４と比較して、相関（Ｒ２）及び相関（ｒ２）の値が若干小さく、且つ、ＳＥＣ、ＳＥＰ及びＥＩの値が若干大きいという結果が得られた。したがって、元素の割合は、Ｋ、Ｃｌ、Ｍｇ、Ｃａ、Ｐ、Ｓから選択される少なくとも２元素以上の元素の割合であることが好ましく、３元素以上の元素の割合であることが特に好ましいことがわかった。
【０１２０】
（実施例２−１）
実施例１−１において乳酸緩衝能を測定した供試試料１７９点を石英ガラスの窓が付いたリング状容器（標準セル）に約４ｇづつ充填し、４００ｎｍ〜２５００ｎｍの拡散反射光スペクトルを近赤外分析計（品番Ｍｏｄｅｌ-６５００、ＦＯＳＳＮＩＲＳｙｓｔｅｍｓ社製）にて測定した。
【０１２１】
測定したデータは、セグメントサイズ１０ｎｍ、ギャップサイズ０で２次微分処理を行った後、実施例１−１と同様に検量線作成用の試料１４０点と検量線の精度評価用の試料３９点に区分した。
【０１２２】
２次微分処理後のデータから波長１４５４ｎｍ、１７３８ｎｍ、２１５２ｎｍ及び１７５８ｎｍを選択し、検量線作成用の試料１４０点について重回帰分析を行い、検量線（すなわち、乳酸緩衝能の推定式）を決定した。次に、評価用試料３９点の２次微分処理後の拡散反射光スペクトルのデータから波長１４５４ｎｍ、１７３８ｎｍ、２１５２ｎｍ及び１７５８ｎｍのスペクトル強度の値を上記決定された検量線（推定式）に代入し、評価用試料３９点の乳酸緩衝能を求めた。
【０１２３】
検量線（推定式）の精度評価の結果を表２に示す。
【０１２４】
（実施例２−２）
実施例２−１において４００ｎｍ〜２５００ｎｍの拡散反射光スペクトルを測定し、２次微分処理を行った１７９点の供試試料を検量線作成用の試料１４０点と検量線の精度評価用の試料３９点に区分した。
【０１２５】
２次微分処理後のデータから波長１１５６ｎｍ〜２４４２ｎｍを選択し、検量線作成用の試料１４０点についてＰＬＳ回帰分析を行い、検量線（乳酸緩衝能の推定式）を決定した。次に、評価用試料３９点の２次微分処理後のデータから波長１１５６ｎｍ〜２４４２ｎｍのスペクトル強度の値を上記決定された検量線に代入し、評価用試料３９点の乳酸緩衝能を求めた。
【０１２６】
検量線（推定式）の精度評価の結果を表２に示す。
【０１２７】
（実施例２−３）
実施例２−１において４００ｎｍ〜２５００ｎｍの拡散反射光スペクトルを測定し、２次微分処理を行った１７９点の供試試料を検量線作成用の試料１４０点と検量線の精度評価用の試料３９点に区分した。
【０１２８】
２次微分処理後のデータから波長１４５０ｎｍ〜１８２０ｎｍ，１９００ｎｍ〜２２００ｎｍを選択し、検量線作成用の試料１４０点についてＰＬＳ回帰分析を行い、検量線（乳酸緩衝能の推定式）を決定した。次に、評価用試料３９点の２次微分処理後のデータから波長１４５０ｎｍ〜１８２０ｎｍ，１９００ｎｍ〜２２００ｎｍのスペクトル強度の値を上記決定された検量線に代入し、評価用試料３９点の乳酸緩衝能を求めた。
【０１２９】
検量線（推定式）の精度評価の結果を表２に示す。
【０１３０】
（比較例２）
実施例２−１において４００ｎｍ〜２５００ｎｍの拡散反射光スペクトルを測定し、２次微分処理を行った１７９点の供試試料を検量線作成用の試料１４０点と検量線の精度評価用の試料３９点に区分した。
【０１３１】
２次微分処理後のデータから波長８５６ｎｍ〜１０４２ｎｍを選択し、検量線作成用の試料１４０点についてＰＬＳ回帰分析を行い、検量線（乳酸緩衝能の推定式）を決定した。次に、評価用試料３９点の２次微分処理後のデータから波長８５６ｎｍ〜１０４２ｎｍのスペクトル強度の値を上記決定された検量線に代入し、評価用試料３９点の乳酸緩衝能を求めた。
【０１３２】
検量線（推定式）の精度評価の結果を表２に示す。
【０１３３】
【表２】
2016176830.tif 000015
【０１３４】
［結果］
実施例２−１〜実施例２−３の何れにおいても、相関（Ｒ２）及び相関（ｒ２）が高い値となり、且つ、ＳＥＣ、ＳＥＰ及びＥＩが低い値に維持されているため、検量線（推定式）の精度が高いことが示された。
【０１３５】
また、実施例２−１と実施例２−２及び実施例２−３との対比によれば、実施例２−２及び実施例２−３は、実施例２−１よりも相関（Ｒ２）及び相関（ｒ２）が高い値となり、且つ、ＳＥＣ、ＳＥＰ及びＥＩが低い値となることから、広範な波長範囲を選択できるＰＬＳＲによる検量線の作成が重回帰分析よりも有利であることが示された。
【０１３６】
さらに、実施例２−２及び実施例２−３と比較例２との対比によれば、比較例２は、実施例２−２及び実施例２−３よりも相関（Ｒ２）及び相関（ｒ２）が大幅に小さく、且つ、ＳＥＣ、ＳＥＰ及びＥＩが大きく上昇していることから、ＰＬＳＲによっても８５６ｎｍ〜１０４２ｎｍの範囲の波長を選択すると精度の高い検量線が得られないことが確認できた。
【０１３７】
なお、実施例２−３を実施例２−２と対比すると、実施例２−３においては若干ＳＥＰ及びＥＩの値が実施例２−２と比べて大きくなるものの、他の値は実施例２−２とほぼ変わらない値となっており、検量線の精度において実施例２−３は実施例２−２と遜色ないことがわかった。
【０１３８】
したがって、ＰＬＳＲにより検量線を作成する場合、選択波長範囲を１１００ｎｍ〜２５００ｎｍとすること、あるいは、１４５０ｎｍ〜１８２０ｎｍ及び１９００ｎｍ〜２２００ｎｍの範囲とすることが好ましいことがわかった。
【０１３９】
よって、サイレージ（発酵産物）原料に対して近赤外分析を行って得られたデータに対して重回帰分析及びＰＬＳ回帰分析を行うことで精度の高い検量線を得ることができることがわかった。
【０１４０】
なお、本発明は上記実施の形態に限定されることはなく、発明の趣旨を逸脱しない範囲で種々変更可能である。
【図１】
2016176830.tif 000016
【図２】
2016176830.tif 000017