Application of Artificial Neural Network and Response Surface Methodology in Optimization of Cold Fiber Headspace Solid-Phase Microextraction (CF-HS-SPME) for Analysis of 2, 5-Hexandion in Urine Samples

authors:

avatar Yasaman Pourbakhshi 1 , avatar Maryam Farhadian 2 , avatar AbdulRahman Bahrami 1 , * , avatar Leila Tajik 3 , avatar Farshid Ghorbani Shanha 1 , avatar Mohammad Javad Assari ORCID 1

Center of Excellence for Occupational Health, Research Center of Health Sciences, School of Public Health, Hamadan University of Medical Sciences, Hamedan, Iran
Department of Biostatistics, School of Public Health, Hamadan University of Medical Sciences, Hamedan, Iran
Department of Occupational Health, School of Health and Nutrition, Lorestan University of Medical Sciences, Khoramabad, Iran

how to cite: Pourbakhshi Y, Farhadian M, Bahrami A, Tajik L, Ghorbani Shanha F, et al. Application of Artificial Neural Network and Response Surface Methodology in Optimization of Cold Fiber Headspace Solid-Phase Microextraction (CF-HS-SPME) for Analysis of 2, 5-Hexandion in Urine Samples. Jundishapur J Health Sci. 2020;12(3):e100925. https://doi.org/10.5812/jjhs.100925.

Abstract

Background:

Cold-Fiber Solid-Phase Microextraction (CF-SPME) is often used for the extraction of volatile and semi-volatile compounds from complex matrices. Multivariate statistical optimization techniques can save time and chemicals and thus decrease the analytical cost in comparison with the Single Variable Approach (SVA). Over the past few decades, different beneficial mathematical tools have been developed for the optimization of separation processes.

Objectives:

In this study, the Artificial Neural Network (ANN) and Response Surface Methodology (RSM) applications were compared for CF-SPME optimization to determine 2, 5-hexandion in urine samples. 2, 5-Hexanedione is a colorless liquid and the main metabolite of n-hexane as a result of occupational exposure.

Methods:

N-hexane is widely used in the rubber industry, food processing, solvents, and medicinal drugs and has adverse effects such as neurotoxicity on humans. Thus, biological monitoring, analytic methods, and mathematical and statistical techniques concerning 2.5-HD are very important. The RSM and ANN are mathematical and statistical techniques applied for the optimization and process modeling. Designing an experiment based on the Historical Data Design (HDD) of RSM was adopted to evaluate the relationship between independent parameters such as extraction temperature, extraction time, sample volume, and extraction efficiency of 2, 5-hexandion.

Results:

The models were compared for their predictive ability by the coefficient of determination (R2) and Root Mean Square Error (RSME) based on the train and test dataset.

Conclusions:

The results were highly significant (P < 0.05) for the optimization of variables. The ANN model had more generalizability than the RSM model. Also, the ANN had higher predictive accuracy.

1. Background

Solid-phase Microextraction (SPME) is a solvent-free sampling technique that includes fiber and an extracting phase (1). Many published studies have described the role of SPME in the field of occupational and environmental health (2, 3). The principal advantage of SPME is the combination of two extraction and pre-concentration steps in one step. It is also a simple, cost-effective, efficient, and environmentally friendly method (4). However, SPME has certain limitations such as poor thermal, mechanical, and chemical stability of fibers along with low selectivity and sample carry-over (1, 5).

To overcome some of these limitations, novel approaches have been developed based on SPME, including ultrasound-assisted SPME, electrochemically assisted SPME, in-tube SPME, and cold fiber SPME (CF-SPME) (5-7). Increasing temperature has a leading role in facilitating the transfer of the analyte from the matrix to the gas phase during conducting the SPME technique. However, this would contribute to the decrease in distribution coefficients of the analytes between the extraction phase (fiber) and the matrix. According to the procedure, the determination sensitivity will decrease. The CF-SPME method has overcome this limitation by heating the matrix and cooling the fiber simultaneously (8-10). According to this method, the partition coefficient does not decrease. The analyte mass transfer from the gas phase to the extraction phase (fiber) will improve; hence, it increases the sorption capability and detection limits for most volatile and semi-volatile analytes (11-13). On the other hand, the optimization of an analytical procedure is an important step in any method development. Accordingly, attempts have been made for the optimization of this method, and several techniques have been developed, as well (14).

The traditional optimization approach is usually based on the classical method called one-variable-at-a-time (15, 16). However, this approach has disadvantages, such as not including the interactions between independent variables and increasing the number of experiments, experiment duration, and amount of chemicals (17, 18). Employing multivariate statistical techniques for the optimization of experiments can save time and chemicals and reduce the cost of analytical experiments (19, 20).

In recent decades, different mathematical tools have been developed for the optimization of separation processes, such as Response Surface Methodology (RSM) and Artificial Neural Network (ANN). These methods are powerful mathematical methods suitable for specifying the experimental conditions that produce the best possible analytical performance (5, 21).

Response surface methodology is a significant tool developed by Box and collaborators (17, 22). It is a statistical method based on empirical models for experimentally analyzing the quantitative data, and it is mostly related to the experimental design (23, 24). The RSM evaluates the effects of variables and their interactions on response variables. This standard approach has already been used in many experiments involving volatile and semi-volatile compounds (25, 26). The RSM principal advantages include the reduction of the number of experimental trials and evaluation of parameters and their possible interactions (27-29). Besides, the Artificial Neural Network (ANN) has been used in combination with the experimental design for representing relationships between variables and predicting the optimum conditions based on the results from a few numbers of experiments (14, 20). The ANN has also been reported as a powerful tool for modeling chemical processes (30, 31).

In this study, 2, 5-hexanedione (2, 5-HD) determination by the CF-SPME method in urine was optimized by RSM and ANN methods. 2, 5-HD is a colorless liquid and the biological index of hexane exposure. Design of Experiment (DOE) was adopted based on the Historical Data Design (HDD) of RSM to evaluate the relationship between independent parameters such as extraction temperature, extraction time, and sample volume, on the one hand, and the extraction efficiency of 2, 5-HD, on the other hand.

2. Objectives

Therefore, this paper aimed to compare the ANN and RSM performance for optimizing this process. The models were compared by the coefficient of determination (R2) and Root Mean Square Error (RSME) for their predictive ability based on the train and test dataset. Analysis of Variance (ANOVA) was used to assess any significant lack-of-fit with experimental results in RSM.

3. Methods

3.1. Instrumentation

SPME fiber assembly Carboxen/polydimethylsiloxane (CAR/PDMS - df 75µm) and manual holder were purchased from Supelco (Bellefonte, PA). Clear glass vials (20 ml volume), sealed with silicone septum and aluminum cap, were used for sample extraction.

An analytical method using capillary gas chromatography (Shimadzu GC-2010, Kyoto, Japan) with Flame Ionization Detector (FID) was developed for the analysis of 2, 5-HD.

For heating and stirring the sample, a hot plate-stirrer (Alfa hs-810 - Tehran, Iran) was employed. Temperature controller (Busan, South Korea), electronic power (220VAC/50/60HZ), thermoelectric cooler (TEC), and related heat sink and fan (Busan, South Korea) were used for cooling the CF-SPME system and control the temperature during the procedure.

3.2. Reagents and Materials

All reagents were of analytical grade. 2, 5-HD, and cyclohexanone (internal standard, IS) were purchased from Merck (Darmstadt, Germany). Stock standard solutions (1000 mg/L of 2, 5-HD, and 5 mg/L of IS) and working standard solutions (20 mg/L of 2, 5-HD for extraction process) were prepared with distilled water. The calibration curves for analysis (0.06–20 μg/ml) were established using working standard solutions prepared from the stock standard solutions with distilled water. 2, 5-HD was added to distilled water.

3.3. SPME Method

Samples were put into 20 mL glass vials containing a certain amount of Na2SO4 (20 w/v %) and a small magnetic stirring bar for HS-CF-SPME optimization. Then, the vials were tightly capped and sealed with an aluminum cap and silicone septum. To control the effect of temperature on extraction, a hot plate and a water bath were used.

Before starting each extraction, the solutions were under a certain extraction temperature and stirring rate of 1000 rpm for 10 min to reach equilibrium. To initiate the extraction process, the SPME fiber was exposed to the sample headspace. Following the extraction process, the fiber injected into the GC injection port directly that was withdrawn into the needle. The CF-SPME system was performed based on a thermoelectric device that created a cooling source for the extraction procedure. The device was running during the procedure.

A copper plate was joined to the heat sink of the Thermoelectric Cooling Device (TEC), and a groove of 0.7 mm depth was made on the surface of the copper plate. The needle from the SPME syringe was placed into the groove, and the fiber was exposed to the headspace of the solution for 17.5 min. A K-type thermocouple was fixed on this copper plate. The temperature of the fiber was around 5 °C. A small fan inside the chamber was used to facilitate air convection (32).

3.4. Experimental Design

Cold fiber headspace solid-phase microextraction (CF-HS-SPME) based on a thermoelectric cooling method for the analysis of 2, 5 hexandion in urine samples was optimized using Response Surface Methodology (RSM) and Artificial Neural Network (ANN). The RSM and ANN are mathematical and statistical techniques applied for optimization and process modeling.

3.5. Historical Data Design (HDD) of RSM

Based on the HDD of RSM, to assess the effects of the main independent variables on the extraction efficiency of 2, 5-HD, the CF-HS-SPME optimization was carried out using RSM. The levels of independent parameters investigated in this research including extraction time, extraction temperature, and sample volume are given in Table 1. The levels of three independent variables were selected based on preliminary work and previous data. According to the quadratic model, the response prediction and coefficient estimations were performed by the least-squares regression.

Y=β0+i=1kβiXi+i=1kβiiXi2+i=1k-1j=2kβijXiXj+e
Table 1.

Experimental Variables and Levels

VariableSymbolUnitMin.Max.
Extraction TemperatureACelsius5080
Extraction TimeBminute530
Sample VolumeCmilliliter515

Where Y is the predicted extraction efficiency of 2, 5-HD, βo is the model intercept coefficient, βi, βii, and βij are the linear, quadratic, and interaction coefficients, respectively, Xi and Xj are the independent variables, and e is the error. The Analysis of Variance (ANOVA) was used to determine the statistical significance of each regression coefficient.

3.6. Artificial Neural Network

A neural network trained by the Back Propagation Algorithm (BPA) was applied. The layers of the network include the input layer, hidden layer, and output layer. According to the research method, a feed-forward back propagation neural network with three layers was used.

In the feed-forward neural network, information flows from input to output. This progress was without feedback. Learning nonlinear and linear relationships between input and output vectors to network happens with multiple layers of neurons with nonlinear transfer functions. According to research, multilayer ANN models with only one hidden layer have universal applications.

Figure 1 illustrates ANN (8: 4: 1) for modeling of 2, 5-HD analysis in urine samples. To determining the number of neurons in the hidden layer, different neural networks with different topologies were fitted to the data, and finally, the best neural network architecture was identified. It is noted that linear transfer function considered for the output layer, and Hyperbolic tangent is considered as a transfer function of the hidden layer. (Neural networks toolbox of Matlab 7.12.0).

Neural Network Architecture (3-4-1) Trained for the Extraction Efficiency of 2, 5-hexandion
Neural Network Architecture (3-4-1) Trained for the Extraction Efficiency of 2, 5-hexandion

4. Results and Discussion

2,5-HD extraction condition was optimized by the RSM. Therefore, three factors were optimized in this design. Table 2 shows the entire 48 runs. According to the quadratic model, the process response and three independent factors of empirical relationships are as follows:

Chromatographyarea=589.90+94.54A+86.56B-244.84C+34.63AB-2.44AC+10.63BC-115.85A2-63.16B2-52.87C2

The order of effectiveness of all model terms is as follows: C > A2 > A > B > B2 > C2 > AB > BC > AC (irrespective of the sign of model coefficients).

There was a coefficient of determination (R2) of 0.897. The R2 from the above quadratic equation was explained by independent variables of extraction temperature (A), extraction time (B), and sample volume (C).

Table 2.

The Coded Historical Design of Independent Factors and Their Corresponding Experimental Values, the RSM Model Predicted and ANN Model Predicted Values of Chromatography Area

RunAaBbCcY (Observedd)RSM PredictedANN PredictedRunABCY (Observedd)RSM PredictedANN Predicted
180515182.45116.002137.812560510433421.507436.04
2801015130.65192.078123.1526601010483.48478.926486.86
3802015297.7295.217291.7427602010652.17543.752659.71
4803015291333.006298.0228603010640.03545.228626.96
570515104.95183.687122.992950510165293.277187.41
6701015153.33249.171147.6930501010202.7340.14194.85
7702015324.25331.127329.0331502010217.23384.747212.87
8703015315347.733357.9632503010210.67364.039209.55
96051569156.729151.46338055600591.332584.9
10601015117.17211.621131.423480105638.3672.461639.3
11602015286.13272.393280.943580205787.23785.706803.04
12603015297267.817293.193680305756.68833.602754.81
135051565.2335.127152.18377055612.7645.761602.95
1450101595.279.428117.233870105654.3716.298659.61
15502015107.25119.017113.753970205818.5858.36822.83
16503015147.193.25673.324070305797.2835.073763.12
1780510449.7394.036415.26416055578.8605.536582.35
18801010496.7472.639491.694260105618.2665.492620.67
19802010662.5580.831666.454360205778.1736.371770.28
20803010649.7623.674644.754460305751.3741.9720.44
2170510502.03455.094439.01455055565.3470.688547.75
22701010545.6523.105520.984650105586.9520.042586.58
23702010694.8610.114678.54750205611.73569.737601.16
24703010658.53631.773660.284850305597.2554.083602.64

The results of the ANOVA are shown in Table 3. All the independent variables were significant (P < 0.05). The F values and P values represent the high significance of the model and no significant lack-of-fit. The variables were not significantly correlated with each other (P > 0.05). The predicted values matched well with the experimental values (R2 = 0.9325 and Adj. R2 = 0.9165).

Table 3.

Analysis of Variance (ANOVA) for Process Responsea,b

ParametersStatistics
Sum of squaresDegree of freedomMean squareF ValueP Value
Model2.576×10692.862×10558.30< 0.0001
A2.346×10512.346×10547.79< 0.0001
B2.121×10512.121×10543.20< 0.0001
C1.886×10611.886×106384.27< 0.0001
AB18872.33118872.333.840.0573
AC106.041106.040.0220.8839
BC2134.4612134.460.430.5136
A21.273×10511.273×10525.92< 0.0001
B233066.55133066.556.740.0134
C229815.63129815.636.070.0184
Residuals1.865×10538
Cor Total2.762×10647

As presented in Figure 2A, the residuals are normally distributed and show a random scatter. Also, there is a linear correlation between predicted and actual response values that indicates a sufficient agreement between real data and model data (Figure 2B). These plots show that the model fits well to optimize the independent variables for the chromatography area prediction.

The model describes the entire experimental range of the study. It shows the comparison of the graphical representation of actual versus predicted values (Figure 2). The results of ANOVA of the linear model for the analysis of 2, 5-HD in urine samples are shown in Table 3.

The model came up with an F value of 69.90 and a P value of < 0.0001. P values of less than 0.05 were considered significant. According to the results, all factors were found to be highly significant (P < 0.05). According to the non-significant lack-of-fit, the model shows good predictability.

Comparison of graphical representation of actual versus predicted values
Comparison of graphical representation of actual versus predicted values

The predicted values were found to match well with the experimental values (R2 = 0.9325). The R2 value was in a reasonable agreement with Adj. R2 (Adj. R2=0.9165), which indicates fair predictability of the model.

The lowest training and verification errors in architecture 3-4-1 of an MLP neural network were permissible. According to the results, it was a suitable network for the prediction of optimization for the analysis of 2, 5-HD in urine samples.

The predictive performance of ANN and RSM models are presented in Table 4. The results showed that the ANN model in both tests and training sets had better performance in terms of R2 and RMSE compared to the RSM model. Smaller values of RMSE for ANN indicate less prediction error for this model. Considering the perfect prediction line, the predictions of the ANN model were closer than those of RSM models. Therefore, in terms of generalization capacity, the ANN model had more generalizability than the RSM model. Also, The ANN had a higher predictive accuracy, which might be related to its general ability to estimate the system nonlinearity. The RSM limitation is in a second-order polynomial.

Table 4.

Comparison Between RSM and ANN

MeasureNeural NetworksRSM
TrainTestTrainTest
RMSE9.0038.7971.8865.74
R20.9990.9920.9530.961

A large number of iterative calculations are required in the ANN model. In contrast, there is a single- step calculation for the response surface model

The most important limitation of the RSM method is that the experimental data are fitted to a polynomial model at the second order, while in practice using a second-order polynomial model is not compatible with all systems. In comparison, although the neural network model is a black box in nature, it can take into account the complexities of systems well, even with the limited number of experiments.

5. Conclusions

In this study, the main factors, extraction temperature (A), extraction time (B), and sample volume (C), were found to be highly significant (P < 0.05). The non-significant lack-of-fit showed good predictability. The ANN and RSM models were the same in predicted values. The results showed that the ANN method is preferable for recording nonlinear behavior and prediction. It had more generalizability, cost, and calculation time than the RSM model. Because of the ANN general ability to approximate the system nonlinearity, it had higher predictive accuracy. Unlike the response surface model, a large number of iterative calculations is required in the ANN model.

Acknowledgements

References

  • 1.

    Shahvandi SK, Banitaba MH, Ahmar H. Development of a new pH assisted homogeneous liquid-liquid microextraction by a solvent with switchable hydrophilicity: Application for GC-MS determination of methamphetamine. Talanta. 2018;184:103-8. [PubMed ID: 29674019]. https://doi.org/10.1016/j.talanta.2018.02.115.

  • 2.

    Llompart M, Celeiro M, García-Jares C, Dagnac T. Environmental applications of solid-phase microextraction. TrAC Trends in Analytical Chemistry. 2019;112:1-12. https://doi.org/10.1016/j.trac.2018.12.020.

  • 3.

    Zeverdegani SK, Bahrami A, Rismanchian M, Shahna FG. Analysis of xylene in aqueous media using needle-trap microextraction with a carbon nanotube sorbent. J Sep Sci. 2014;37(14):1850-5. [PubMed ID: 24802025]. https://doi.org/10.1002/jssc.201400262.

  • 4.

    Bakaikina NV, Kenessov B, Ul'yanovskii NV, Kosyakov DS. Quantification of transformation products of rocket fuel unsymmetrical dimethylhydrazine in soils using SPME and GC-MS. Talanta. 2018;184:332-7. [PubMed ID: 29674050]. https://doi.org/10.1016/j.talanta.2018.02.047.

  • 5.

    Memarian E, Hosseiny Davarani SS, Nojavan S, Movahed SK. Direct synthesis of nitrogen-doped graphene on platinum wire as a new fiber coating method for the solid-phase microextraction of BXes in water samples: Comparison of headspace and cold-fiber headspace modes. Anal Chim Acta. 2016;935:151-60. [PubMed ID: 27543024]. https://doi.org/10.1016/j.aca.2016.06.048.

  • 6.

    Kumar A, Malik AK, Tewary DK, Singh B; Gaurav. A review on development of solid phase microextraction fibers by sol-gel methods and their applications. Anal Chim Acta. 2008;610(1):1-14. [PubMed ID: 18267134]. https://doi.org/10.1016/j.aca.2008.01.028.

  • 7.

    Piri-Moghadam H, Ahmadi F, Pawliszyn J. A critical review of solid phase microextraction for analysis of water samples. TrAC Trends in Analytical Chemistry. 2016;85:133-43. https://doi.org/10.1016/j.trac.2016.05.029.

  • 8.

    Asfaw AA, Aspromonte J, Wolfs K, Van Schepdael A, Adams E. Overview of sample introduction techniques prior to GC for the analysis of volatiles in solid materials. J Sep Sci. 2019;42(1):214-25. [PubMed ID: 30290055]. https://doi.org/10.1002/jssc.201800711.

  • 9.

    Fedotov PS, Malofeeva GI, Savonina E, Spivakov B. Solid-Phase Extraction of Organic Substances: Unconventional Methods and Approaches. Journal of Analytical Chemistry. 2019;74(3):205-12. https://doi.org/10.1134/s1061934819030043.

  • 10.

    Xu S, Shuai Q, Pawliszyn J. Determination of Polycyclic Aromatic Hydrocarbons in Sediment by Pressure-Balanced Cold Fiber Solid Phase Microextraction. Anal Chem. 2016;88(18):8936-41. [PubMed ID: 27562230]. https://doi.org/10.1021/acs.analchem.6b01944.

  • 11.

    Kenessov B, Koziel JA, Bakaikina NV, Orazbayeva D. Perspectives and challenges of on-site quantification of organic pollutants in soils using solid-phase microextraction. TrAC Trends in Analytical Chemistry. 2016;85:111-22. https://doi.org/10.1016/j.trac.2016.04.007.

  • 12.

    Abd El-Kader EM, Serag A, Aref MS, Ewais EE, Farag MA. Metabolomics reveals ionones upregulation in MeJA elicited Cinnamomum camphora (camphor tree) cell culture. Plant Cell, Tissue and Organ Culture (PCTOC). 2019;137(2):309-18. https://doi.org/10.1007/s11240-019-01572-z.

  • 13.

    Dias CM, Menezes HC, Cardeal ZL. Environmental and biological determination of acrolein using new cold fiber solid phase microextraction with gas chromatography mass spectrometry. Anal Bioanal Chem. 2017;409(11):2821-8. [PubMed ID: 28204886]. https://doi.org/10.1007/s00216-017-0226-9.

  • 14.

    Hernández-Borges J, Rodríguez-Delgado MA, García-Montelongo FJ. Optimization of the Microwave-Assisted Saponification and Extraction of Organic Pollutants from Marine Biota Using Experimental Design and Artificial Neural Networks. Chromatographia. 2006;63(3-4):155-60. https://doi.org/10.1365/s10337-006-0722-0.

  • 15.

    Hu YY, Zheng P, He YZ, Sheng GP. Response surface optimization for determination of pesticide multiresidues by matrix solid-phase dispersion and gas chromatography. J Chromatogr A. 2005;1098(1-2):188-93. [PubMed ID: 16289085]. https://doi.org/10.1016/j.chroma.2005.09.093.

  • 16.

    Haji A, Qavamnia SS. Response surface methodology optimized dyeing of wool with cumin seeds extract improved with plasma treatment. Fibers and Polymers. 2015;16(1):46-53. https://doi.org/10.1007/s12221-015-0046-5.

  • 17.

    Bezerra MA, Santelli RE, Oliveira EP, Villar LS, Escaleira LA. Response surface methodology (RSM) as a tool for optimization in analytical chemistry. Talanta. 2008;76(5):965-77. [PubMed ID: 18761143]. https://doi.org/10.1016/j.talanta.2008.05.019.

  • 18.

    Wu Y, Cui SW, Tang J, Gu X. Optimization of extraction process of crude polysaccharides from boat-fruited sterculia seeds by response surface methodology. Food Chemistry. 2007;105(4):1599-605. https://doi.org/10.1016/j.foodchem.2007.03.066.

  • 19.

    Xi J, Xue Y, Xu Y, Shen Y. Artificial neural network modeling and optimization of ultrahigh pressure extraction of green tea polyphenols. Food Chem. 2013;141(1):320-6. [PubMed ID: 23768364]. https://doi.org/10.1016/j.foodchem.2013.02.084.

  • 20.

    Bingol D, Hercan M, Elevli S, Kilic E. Comparison of the results of response surface methodology and artificial neural network for the biosorption of lead using black cumin. Bioresour Technol. 2012;112:111-5. [PubMed ID: 22425399]. https://doi.org/10.1016/j.biortech.2012.02.084.

  • 21.

    Haddadi SH, Niri VH, Pawliszyn J. Study of desorption kinetics of polycyclic aromatic hydrocarbons (PAHs) from solid matrices using internally cooled coated fiber. Anal Chim Acta. 2009;652(1-2):224-30. [PubMed ID: 19786184]. https://doi.org/10.1016/j.aca.2009.05.026.

  • 22.

    Wang S, Chen F, Wu J, Wang Z, Liao X, Hu X. Optimization of pectin extraction assisted by microwave from apple pomace using response surface methodology. Journal of Food Engineering. 2007;78(2):693-700. https://doi.org/10.1016/j.jfoodeng.2005.11.008.

  • 23.

    Li H, Zhang Z, Hou T, Li X, Chen T. Optimization of ultrasound-assisted hexane extraction of perilla oil using response surface methodology. Industrial Crops and Products. 2015;76:18-24. https://doi.org/10.1016/j.indcrop.2015.06.021.

  • 24.

    Li F, Yang L, Zhao T, Zhao J, Zou Y, Zou Y, et al. Optimization of enzymatic pretreatment for n-hexane extraction of oil from Silybum marianum seeds using response surface methodology. Food and Bioproducts Processing. 2012;90(2):87-94. https://doi.org/10.1016/j.fbp.2011.02.010.

  • 25.

    Akalın MK, Tekin K, Akyüz M, Karagöz S. Sage oil extraction and optimization by response surface methodology. Industrial Crops and Products. 2015;76:829-35. https://doi.org/10.1016/j.indcrop.2015.08.005.

  • 26.

    Ma QL, Hamid N, Bekhit AED, Robertson J, Law TF. Optimization of headspace solid phase microextraction (HS-SPME) for gas chromatography mass spectrometry (GC–MS) analysis of aroma compounds in cooked beef using response surface methodology. Microchemical Journal. 2013;111:16-24. https://doi.org/10.1016/j.microc.2012.10.007.

  • 27.

    Wei ZJ, Liao AM, Zhang HX, Liu J, Jiang ST. Optimization of supercritical carbon dioxide extraction of silkworm pupal oil applying the response surface methodology. Bioresour Technol. 2009;100(18):4214-9. [PubMed ID: 19414250]. https://doi.org/10.1016/j.biortech.2009.04.010.

  • 28.

    Rezzoug SA, Boutekedjiret C, Allaf K. Optimization of operating conditions of rosemary essential oil extraction by a fast controlled pressure drop process using response surface methodology. Journal of Food Engineering. 2005;71(1):9-17. https://doi.org/10.1016/j.jfoodeng.2004.10.044.

  • 29.

    Guo X, Zou X, Sun M. Optimization of extraction process by response surface methodology and preliminary characterization of polysaccharides from Phellinus igniarius. Carbohydrate Polymers. 2010;80(2):344-9. https://doi.org/10.1016/j.carbpol.2009.11.028.

  • 30.

    Khajeh M, Moghaddam MG, Shakeri M. Application of artificial neural network in predicting the extraction yield of essential oils of Diplotaenia cachrydifolia by supercritical fluid extraction. The Journal of Supercritical Fluids. 2012;69:91-6. https://doi.org/10.1016/j.supflu.2012.05.006.

  • 31.

    Torrecilla JS, Mena ML, Yáñez-Sedeño P, García J. Application of artificial neural network to the determination of phenolic compounds in olive oil mill wastewater. Journal of Food Engineering. 2007;81(3):544-52. https://doi.org/10.1016/j.jfoodeng.2006.12.003.

  • 32.

    Tajik L, Bahrami A, Ghiasvand A, Shahna FG. Determination of benzene, toluene, ethylbenzene and xylene in field and laboratory by means of cold fiber SPME equipped with thermoelectric cooler and GC/FID method. Polish Journal of Chemical Technology. 2017;19(3):9-15. https://doi.org/10.1515/pjct-2017-0041.