利用主成分回歸分析建立微量成分與菜籽油煙點的數學方程

方 冰 欒 霞 李曉寧 王瑩輝

(國家糧食局科學研究院油脂化學與增值加工技術研究組，北京 100037)

利用主成分回歸分析建立微量成分與菜籽油煙點的數學方程

方 冰 欒 霞 李曉寧 王瑩輝

(國家糧食局科學研究院油脂化學與增值加工技術研究組，北京 100037)

隨著食用植物油適度精煉理念的提出和被廣泛認可，植物油國標中對酸值、色澤等指標開始放寬，而這種放寬會導致植物油中其他的非甘油三酯組分的含量，尤其是對煙點有影響的磷脂、溶劑、過氧化物、水分及揮發物等組分發生變化，進而影響植物油的煙點。盡管針對游離脂肪酸這一將放寬的質量指標，其與煙點間數學關系已有研究，但酸值的放寬引起的其他組分的變化對煙點的影響并未能綜合考量。因此，通過實驗室模擬的方法，采集不同游離脂肪酸、磷脂、過氧化物、溶劑殘留、水分及揮發物的菜籽油，進行復配，獲得組分不同的植物油樣品，測定煙點及各指標含量，利用主成分回歸的方法得到各因素與煙點間的回歸方程：ln煙點=5.006-0.029×ln酸值-0.041×ln過氧化值-0.003×ln溶殘-0.039×ln含磷量-0.044×ln水分及揮發物，R2=0.926，結果表明，各因素的影響程度大小由大到小分別為水分及揮發物、過氧化值、含磷量、酸值、溶殘。

游離脂肪酸 過氧化物 溶劑殘留 磷 水分及揮發物 多重共線性

在我國，植物油的食用方式主要為烹炒和煎炸，油脂在煎炸和烹飪過程中會產生油煙，危害健康。因而，植物油除滿足包括酸值、過氧化值、含磷量等質量指標外，還應具備一定的熱穩定性，來滿足消費者對于植物油使用過程中有較少油煙的要求。煙點這一指標便可用于衡量油脂的熱穩定性。

煙點是指油脂在受熱時肉眼能看見樣品的熱分解物或雜質連續揮發的最低溫度[1-2]，反映了油脂中的甘油三酯及非甘油三酯組分因沸點不同在加熱過程中先后揮發的特性。未精煉油脂的煙點遠低于精煉植物油的煙點，這是由于未經精煉處理的植物油中，諸如游離脂肪酸、磷脂、氧化產物、水分、溶劑等沸點低于甘油三酯的組分含量較高，因此，隨著精煉過程中對這些非甘油三酯組分的去除，植物油的煙點逐漸升高[2-4]。

然而，隨著適度精煉概念的提出，植物油國標中對酸值、色澤等指標開始放寬，這種放寬會導致植物油中非甘油三酯組分的含量，尤其是對煙點有影響的磷脂、溶劑、過氧化物、水分及揮發物等組分發生變化，進而影響植物油的煙點。王月娥等[5]研究發現，煙點與色澤成反比關系，即煙點高的油脂色澤淺，煙點低的則色澤深。王興國等[4]研究發現，油脂酸值若由0.2升至0.5 mgKOH/g油，其煙點將下降16 ℃；含磷量每增加30 mg/kg，煙點下降 4 ℃；氧化產物、單甘酯及抗氧化劑如BHA、BHT及TBHQ對煙點也都有不同程度的影響。盡管針對游離脂肪酸這一變化較大的指標已開展其與煙點間數學關系的研究，不管是唐萍華等[6]及吳政宙等[7]基于現行國標中0.2 mgKOH/g油的酸值范圍內得到的煙點與游離脂肪酸百分含量間的回歸方程還是本課題組前期研究得到的針對4.0 mgKOH/g油的酸值范圍內得到游離脂肪酸百分含量與煙點降低百分比的回歸方程，均不能很好地反映適度精煉條件下，各組分含量變化對煙點的影響。因此，本研究通過實驗室模擬的方法，通過采集不同游離脂肪酸、磷脂、過氧化物、溶劑殘留、水分及揮發物的菜籽油，進行復配，獲得組分不同的植物油樣品，測定煙點及各指標含量后，利用主成分回歸的方法建立各因素與煙點間的回歸方程，從而為適度精煉條件下的工藝參數及組分含量的建立提供理論基礎。

1 材料與方法

1.1 試驗材料

試驗用一級、四級菜籽油取自國內植物油生產廠家，N-N-二甲基乙酰胺購于Sigma公司，六號溶劑為國內市售標準物質，乙醚、乙醇、異辛烷、冰醋酸、鹽酸、濃硫酸、氧化鋅、氫氧化鉀、鉬酸鈉、硫酸聯氨、磷酸二氫鉀等均為市售分析純樣品(北化)。

1.2 不同組分油樣的配制

根據采集的不同批次的菜籽油的游離脂肪酸、磷脂、過氧化物、溶劑殘留、水分及揮發物的含量，按照不同質量比進行植物油的復配，獲得游離脂肪酸、磷脂、過氧化物、溶劑殘留、水分及揮發物組分差異較大的植物油樣品，密封，-20 ℃儲藏待用。

1.3 基本理化指標的測定

試驗中所用菜籽油樣品的酸值、過氧化值、含磷量、水分及揮發物、溶劑殘留分別按照GB/T 5530—2005、GB/T 5538—2005、GB/T 5537—2008、GB/T 5528—2008、GB/T 5009.37—2003中方法測定，每個樣品重復3次。

1.4 煙點的測定

植物油樣品的煙點按照GB/T 20795—2006中方法測定，油脂煙點儀(HLY-Ⅲ)為杭州錢江儀器設備有限公司，每個樣品重復3次。

2 結果與討論

2.1 不同組分的菜籽油樣品的組成及煙點

從國內廠家采集不同等級的菜籽油，測定其初始酸值、過氧化值、含磷量、溶劑殘留、水分及揮發物，按照不同質量比進行配制，獲得組分差異較大的菜籽油樣品，所得樣品中各組分含量及對應煙點如表1所示。由表1中可以看出，通過復配得到了酸值范圍在0.120～2.760 mgKOH/g 油，過氧化值范圍在0.500～3.900 mmol/kg，溶劑殘留范圍在0.000～3.700 mg/kg，含磷量在0.247～2.280 mg/100 g油，水分及揮發物范圍在0.010%～0.080%內的88個菜籽油樣品，其煙點范圍為151.0～218.5 ℃。

表1 菜籽油樣品的煙點與各組分含量

表1(續)

表1(續)

注：數據為3次測定的平均值。

2.2 煙點各影響因素間多重共線性的檢驗

在進行回歸分析時，為排除不同變量間數值上的差異，先將所有數據取對數，一方面可以降低樣本異方差程度，另一方面可以減少變量的波動，與其他變量的波動水平相適應[8]，再進行后續分析。因各精煉工序段均會影響游離脂肪酸、氧化產物、磷脂、溶劑殘留、水分及揮發物等組分的含量，表1中各指標間很可能存在交互作用，即影響煙點的各因素間具有多重共線性。多重共線性是影響回歸模型穩定性與結論可靠性的重要因素，故對原數據首先進行共線性的檢驗，結果如表2所示。

自變量的容忍度越小，該自變量與其他自變量的共線程度就越大(接近1較好)；方差膨脹因子不小于10時，該自變量與其他自變量間存在嚴重的多重共線性[9]。由表2可知，溶殘這一變量被排除，其他4個自變量酸值、過氧化值、含磷量、水分及揮發物的容忍度均較小，方差膨脹因子也遠遠超過10，表明影響煙點的因素中，酸值、過氧化值、含磷量、水分及揮發物這些自變量間存在極大的共線性。克服多重共線性有擴大樣本容量、剔除存在多重共線性的變量、嶺回歸、主成分回歸、偏最小二乘法等多種方法[10]，其中，主成分是標準化原始變量的線性或非線性組合，具有可逆性，可還原為原始變量的回歸方程，故本課題則選用主成分回歸的方法進行分析。

表2 回歸系數表

注：被解釋變量為煙點。

2.3 煙點各影響因素的主成分回歸分析

主成分回歸，是指通過適當的線性或非線性變換，將若干相關性較強的解釋變量綜合為彼此獨立的新變量，再根據新變量的方差大小(代表有效信息量)，選取若干方差最大的新變量(主成分)為解釋變量對被解釋變量進行回歸[11]。主成分回歸建立在主成分分析之上，將主成分分析得到的互不相關的若干主成分與因變量進行回歸，得到主成分回歸模型，再通過數據轉換，將標準化的y及主成分，變換回原水平的y和xi，得到一般線形回歸模型。

在進行回歸分析前，應先對表1中數據進行標準化處理(表3)，對各標準化后的數據利用SPSS軟件進行主成分回歸分析，得到總方差解釋表(表4)和主成分法提取因子載荷陣初始解(表5)。選擇特征值大于1的成分作為主成分，由表4可知，成分1的特征值為3.961，貢獻率達到79.220%，故選擇成分1作為變量，利用SPSS的“Compute Variable”程序，按照表5輸入各成分的公式，可得到：C1=0.966×Z(ln酸值)+0.984×Z(ln過氧化值)+0.959×Z(ln溶殘)+0.995×Z(ln含磷量)+0.387×Z(ln水分及揮發物)式(1)，與Z(ln 煙點)建立回歸方程，結果如表6所示。由表6可以得出回歸方程：Z(ln煙點)=-0.247×C1(R2=0.926)，將式(1)帶入可得，Z(ln煙點)=-0.238×Z(ln酸值)-0.239×Z(ln過氧化值)-0.085×Z(ln溶殘)-0.237×Z(ln含磷量)-0.245×Z(ln水分及揮發物)式(2)

表3 描述性統計量表

表4 總方差解釋表

表5 主成分法提取因子載荷陣初始解

表6 線性回歸方程系數表及模型評估

ln煙點=5.006-0.029×ln酸值-0.041×ln過氧化值-0.003×ln溶殘-0.039×ln含磷量-0.044×ln水分及揮發物。

由方程可以看出，各自變量按照對煙點的影響程度從大到小的順序分別為水分及揮發物、過氧化值、含磷量、酸值、溶殘。

表7 各變量的相關系數表

注: ** 相關性系數在0.01的水平下具有顯著性差異。

3 結論

本研究考慮到精煉過程中，植物油中非甘油三酯組分的變化對煙點的綜合影響，建立了煙點與酸值、過氧化值、溶殘、含磷量、水分及揮發物間的回歸方程，并明確影響因素的影響程度的高低，最終發現各因素的對數值存在以下的數學關系：ln煙點=5.006-0.044×ln水分及揮發物-0.041×ln過氧化值-0.039×ln含磷量-0.029×ln酸值-0.003×ln溶殘(R2=0.926)。本研究為適度精煉工藝條件下，植物油中組分變化與其熱穩定性評價指標煙點的界定的一致性提供了理論依據。

[1]陸健，黃衛寧，冀聰偉，等. 煎炸油熱穩定性及氧化穩定性影響因素分析[J]. 農業機械，2012，3： 47-49 Lu J, Huang W N, Ji C W, et al. Analysis of the factors that influence the thermal stability and oxidation stability of frying oil[J].Farm Machinery, 2012,3:47-49

[2] 韓景生，葉彩文，韓磊. 淺談油脂的煙點與品質[J]. 中國油脂，1997， 22(2)：14-15 Han J S, Ye C W, Han L. The relationship between the smoke point and quality of oil[J].China Oils and Fats, 1997,22(2):14-15

[3] 韓山山，孟宗，李進偉，等. 磷脂與游離脂肪酸對油脂煙點和氧化穩定性的影響[J]. 中國油脂，2014，39(4)：23-26 Han S S, Meng Z, Li J W, et al. Effects of phospholipids and free fatty acids on smoke point and oxidative stability of oils[J].China Oils and Fats, 2014,39(4): 23-26

[4] 王興國，周勝利，徐霞，等. 非甘油三酯成分對油品煙點的影響[J]. 中國油脂，1999，24(1)：38-39 Wang X G, Zhou S L, Xu X, et al. Effect of non-triglyceride compounds on the smoke point of oil[J].China Oils and Fats, 1999, 24(1): 38-39

[5]王月娥，王立霞，王月華. 關于大豆油煙點的高低與大豆油質量的關系[J]. 黑龍江糧油科技，1997，3： 60-61

Wang Y E, Wang L X, Wang Y H. The relationship betweensmoke pointand quality of soybean oil[J].Journal of Heilongjiang Cereals and Oils Science and Technology, 1997,3:60-61

[6]唐萍華. 精制大豆油煙點與游離脂肪酸含量的關系模型[J]. 糧食與油脂，2003，增刊: 7-9 Tang P H. Relationshipmodel between the smoking point and the content of free fatty acid in the refined soybean oil[J]. Cereals and Oils, 2003, S: 7-9

[7]吳政宙，張軍華. 影響大豆油煙點的因素及控制[J]. 糧油食品科技，2006，14(2)：20-21 Wu Z Z, Zhang J H. Influence and control the smoking point in the refined soybean oil[J].Science and Technology of Cereals,Oils and Foods, 2006, 14(2): 20-21

[8]施久玉，柴艷有. 灰色成分數據模型在中國產業結構分析預測中的應用[J]. 統計與信息論壇，2006，22(1)：32-35 Shi J Y, Chai Y Y. The grey compositional data model and its application in the analysis and forecast of China′s industrial structure[J].Statistics and Information Tribune, 2006,22(1): 32-35

[9]曾小彬，劉凌娟. 城鄉居民收入差距影響因素及其作用的再分析——基于“一連串事件”邏輯闡述的實證分析[J]. 財經研究，2008，34(12)：118-128 Zeng X B, Liu L J. Re-analysis on influencing factors of the incomegap between urban and rural residents in chinaand their effects: an empirical analysis based on thelogic explanation of''a series of events''[J].Journal of Finance and Economics, 2008,34(12): 118-128

[10]高惠旋. 處理多元線性回歸中自變量共線性的幾種方法[J]. 數理統計與管理，2000，20(5)：49-55 Gao H X. Some method on treating the collinearity of independent variables in multiple linear regression[J].Application of Statistics and Management, 2000,20(5): 49-55

[11]劉潤幸，蕭燦培，宮齊，等. 利用SPSS進行主成分回歸分析[J]. 數理醫藥學雜志，2001，14(2)：103-105 Liu R X, Xiao C P, Gong Q, et al. Principal component regression analysis using SPSS[J].Journal of Mathematical Medicine, 2001,14(2): 103-105.

Establishment of Mathematical Equations of Microconstituents and Smoke Point in Rapeseed Oil Based on Principal Component Regression Analysis

Fang Bing Luan Xia Li Xiaoning Wang Yinghui
(Academy of State Administration of Grain, Lipid Chemistry and Value Added Processing Technology Research Group, Beijing 100037)

As the raising and wide acceptance of concept of moderate refining of edible vegetable oil, indexes for acid value, color, etc. in the vegetable oil national standard starts to be loosened, leading to the content change of other components except for triglyceride in the vegetable, especially phospholipid, solvent, peroxide, moisture, volatile matter, etc, which have an effect on smoke point. Although the mathematical relationship between smoke point and free fatty acid has been studied aiming at free fatty acid belonging to quality index to be loosened, influence of change of other components caused by loosening of acid value on smoke points haven't gained comprehensive consideration. Therefore, in this study, rapeseed oils with different contents of free fatty acid, phospholipid, peroxide, residual solvent, moisture and volatile matter were collected and compounded to generate vegetable oil samples with different components. After measuring smoke point and contents of various indexes, use principal component regression method to gain regression equation between various element and smoke point: Ln (smoke point) =5.006-0.029×Ln (acid value)-0.041×Ln (peroxide value) -0.003×Ln (residual solvent content)-0.039×Ln (phosphorus content)-0.044×Ln (moisture and volatile matter) (R2=0.926). The impact of various factors was the moisture and volatile matter, peroxide value, phosphorus content, acid value and solvent residue in an order from big to small.

free fatty acid, peroxide, residual solvent, phosphorus, water and volatile matter, multicollinearity

糧食公益性行業科研專項(201313011-7)

2015-08-17

方冰，女，1987年出生，助理研究員，油脂加工與營養

TS221

A

1003-0174(2017)04-0127-06