Regression Models for Description of Roasted Ground Coffee Powder Color Change during Secondary Shelf-Life as Related to Storage Conditions and Packaging Material

Besides sensory attributes, color is a parameter affecting consumers’ perception of the powdered coffee product or brew. The aim of this study was to develop and compare non-linear and linear regression models for the description of experimentally determined color changes during 6 months of storage in two different packaging materials. Model parameters were estimated using two software packages: Eureqa Formulize (Nutonian, Inc., Boston, MA, USA) and Statistica 10.0 (StatSoft, Palo Alto, CA, USA) and compared based on their R2 goodness of fit. Both non-linear and linear models used in this study pointed to a significant influence of intrinsic (sample moisture content) and external (relative humidity (RH) and temperature) factors on ground roasted coffee color change. Non-linear model was the most suitable for description of color changes during storage. Based on lower moisture sorption of the sample packed in triplex bag, triplex packaging is proposed as the more suitable one.


Introduction
Ground roasted coffee, due to its exquisite taste, aroma and psychoactive properties, represents the most widely used traditional beverage.It is consumed on regular basis as part of a daily routine by many people, mostly due to the stimulatory effect of caffeine.However, it also contains carbohydrates, lipids, vitamins, minerals, and phenolic compounds with numerous health effects.When discussing the health effect of coffee, numerous conflicting research findings have arisen [1].For example, high levels of caffeine cause a rise in blood pressure and increased metabolic rate, which can have a negative impact on the overall physiological state.On the other hand, polyphenolic compounds present in coffee, mostly chlorgenic acid, is known for its beneficial health effects such as reduction of blood pressure and anti-inflammatory effects [2].
The quality of coffee is strictly related to the chemical composition of the coffee beans, which is affected by factors such as growth climate, and harvesting and post-harvesting conditions (drying, roasting, grinding, and storage) [3].After initial harvesting and drying, coffee beans are subjected to roasting and, subsequently, grinding.One of the methods used for the assessment of coffee quality and the progress of the roasting procedure is color analysis [4].Besides roasting, color analysis can also be used as a quality predictor during secondary shelf life, where, due to package opening and susceptibility to ambient light, air and temperature, coffee degradation reactions can occur [5].Prediction of changes occurring in coffee products during storage is of great importance since these changes have the greatest influence on the consumers' perception on the product quality.Furthermore, consumers' perception and sensory evaluation of coffee is not related solely to aroma and odor, but also to the color and appearance of the coffee powder as well as the prepared brew [6,7].
Color degradation in different food products during processing and storage has been widely studied.Most of these studies refer to color degradation in relation to processing conditions [8][9][10][11], while some of them are oriented towards creating kinetic models for color degradation which occur during storage [5,12,13].The importance of modelling the secondary shelf life of food products in general, lies in the understanding of their behavior after package opening and, thus, optimizing their properties, quality, and packaging material in order to maximize its shelf life, where the role of prediction models is especially emphasized [5].
The aim of this study was to develop and compare regression models for experimentally determined color changes during 6 months of storage of ground and roasted coffee packed in two types of packaging: tin cans and triplex (PE/Al/ PET) bags.Linear and non-linear regression models were proposed for changes in L*, a*, b*, ∆E, Chroma, and Hue.Model parameters were estimated using two software packages: Eureqa Formulize (Nutonian, Inc., Boston, MA, USA) and Statistica 10.0 (StatSoft, Palo Alto, CA, USA).The obtained models were compared based on their R 2 goodness of fit.

Ground Coffee Samples
Ground roasted coffee samples were obtained from Frank d.o.o.(Zagreb, Croatia).Samples were originally packed in two types of packaging: tin cans (L1KMK) and triplex (PE/Al/PET) bags (V1KMK).After the first package opening, samples were stored at ambient conditions in a dark cupboard, in order to simulate home storage.Packaging was opened once a month, for 6 subsequent months of storage, during which moisture content of the samples, relative air humidity, air temperature, and color parameters were determined.

Moisture Content
Moisture content was determined according to the modified method of Haugaard-Sørensen et al. [14].Samples were dried at 105 • C for 3 h in an oven dryer (INKO, Zagreb, Croatia) and weighed on an analytical balance (Sartorius, Göttingen, Germany) prior to and after drying.The difference in weight before and after drying was recorded as the mass of the water contained in the sample.Measurements were done in triplicate.

Ambient Conditions
During each color measurement, ambient conditions (relative air humidity-RH and air temperature-T) were recorded using a digital multimeter (Hyundai corp., Seoul, Korea).Measurements were done in triplicate.

Color Measurements
Color measurements were performed using the CM-700d spectrophotometer (Konica Minolta Sensing INC., Osaka, Japan).Prior to measurement white calibration was conducted using a white calibration plate supplied by the manufacturer.Color was recorded based on the CIE-lab L*, a*, b* scale throughout 6 months of storage.According to Bicho et al. [4] and McGuire [15], the coordinated L* represents lightness (contribution of black or white varying between 0 and 100), a* represents the contribution of green and red (positive or negative), and b* represents the contribution of blue or yellow (negative or positive).Besides L*, a* and b* values, total color difference, Chroma and Hue angle were calculated according to Equation (1) (∆E), Equation (2) (Chroma) and Equation (3) (Hue angle): where subscript "0" refers to the initial color reading of the sample after the first package opening [16][17][18].All measurements were done in triplicate.

Nonlinear Regression Models
To describe the coffee color components (L*, a*, b*, ∆E, Chroma, Hue angle) change during storage, non-linear regression models were proposed.Color components change was modelled as a function of sample moisture content SMC (%), relative air humidity, RH (%), and air temperature, T ( • C) (Equations ( 4) and ( 5)).Regression models for two packaging types were proposed.
Parameters of nonlinear regression models were estimated using evolution algorithm implemented in Eureqa Formulize 0.96 Beta (Nutonian INC, Boston, MA, USA) and Levenberg-Marquardt algorithm implemented in Statistica 10.0 (StatSoft, Palo Alto, CA, USA).Evolutionary algorithms are heuristics inspired by Darwin's theory of evolution.They use the principles of selection and reproduction of the best adapted individuals.The aim of an evolutionary algorithm is to optimize a given objective function F over a given search space.F is called the fitness function.Algorithm creates a population of "individuals", where each individual is a set of parameters in the search space [19].The Levenberg-Marquardt algorithm ensures numerical solutions in function parameter space using the least squares method.Calculation in Statistica 10.0 was performed with 50 iterations with convergence of 10 −6 and confidence interval of 95%.

Linear Regression Models
Natural logarithm form of Equation ( 5) was used to obtain linear regression models for description of color components change (Equation ( 6)).

ln(colour components change
Parameters of linear regression models were also estimated using the evolutionary algorithm implemented into the Eureqa Formulize 0.96 Beta (Nutonian, INC) and using the least squares algorithm implemented into the Statistica 10.0 (StatSoft, Palo Alto, CA, USA).

Experimental Data
One of the main characteristics of coffee is a significant change in color during the roasting process [20].After roasting, ground roasted coffee samples can be classified according to their color as light (L* > 35), medium (25 < L* < 35), or dark (L* < 25) roast [21].During secondary shelf-life, changes occur in lipid composition, which affect sensory properties and color of the coffee.Furthermore, while during roasting temperature and exposure time play a crucial role in sensory and color changes, during storage moisture content, relative humidity and temperature seem to be the most important degradation factors [22].Experimental data determined during 6 months of storage of ground roasted coffee samples is shown in Table 1.Based on the experimental data shown in Table 1, L* values of all ground roasted coffee samples ranged from medium to dark roast.a* values ranged from 11.00 to 12.62 for the coffee sample packed in triplex bag; and between 11.25 and 12.46 for coffee sample packed in tin can, which was slightly higher in comparison to the triplex bag sample.a* and b* values for both packaging materials were positive and concentrated near the center of the planar plane, exhibiting the dominant red and yellow components of the tested samples.The sample packed in triplex bag exhibited a decrease in a* values, while the coffee sample packed in the tin can exhibited the opposite effect.Both samples exhibited a decrease in b* values.According to Rendon et al. [23], an increase in L*, a* and b* values during storage of roasted coffee beans was observed, which means that the brown color component fades during storage.In this case, different behavior of the a* color parameter can be attributed to different packaging material.According to Kallio et al. [24], types of packaging and permeability of packaging are of primary importance when coffee storage stability is addressed.When stored in a permeable package, coffee is especially subjected to rapid changes, mostly in the contents of volatile compounds responsible for its aroma.
During total color change calculation, L*, a* and b* parameters determined immediately after the first package opening were used as reference values.As shown in Table 1, total color difference was more pronounced for the coffee sample packed in the triplex bag.
The Chroma value indicates the degree of saturation of color and is proportional to the strength of the color [18].Chroma values decreased for both packaging materials used in this study, with a more pronounced decrease for triplex bag packaging.
The Hue angle is another parameter frequently used to characterize color in food products.An angle of 0 • or 360 • represents red Hue, while angles of 90 • , 180 • and 270 • represent yellow, green and blue Hue, respectively [16].Hue angle of the coffee samples determined in this study, expressed in degrees, ranged from approximately 55 • to 60 • , which corresponded to the red Hue.
Moisture content values exhibited a consecutive rise throughout storage.The coffee sample packed in the tin can absorbed higher amounts of moisture due to larger headspace air volume present in the tin can packaging.The presence of air in the headspace is also known to contribute to a loss of coffee flavor [25].Ambient conditions (relative humidity and temperature values) are also listed in Table 1, and were used for regression modelling, together with the moisture content measured directly in the samples by an oven drying method.

Non-Linear Regression Models
To analyze how the color components (L*, a*, b*, ∆E, Chroma, and Hue angle) of the roasted ground coffee change during secondary shelf-life, regression models were proposed.According to Equation (4) change in the color components depends on the sample moisture content, relative air humidity and air temperature.The influence of the packaging was also taken into account and regression models were developed for each packaging material separately.Parameters of both non-linear and linear regression models were estimated using two numerical methods, evolutionary algorithm implemented into the Eureqa Formulize software, and the least square algorithm implemented into the Statistica 10.0.software.The estimated values of model parameters are presented in Tables 2-5 and results of the dispersion of color components in comparison to model predictions are presented in Figures 1-4.
Eureqa Formulize estimated the model parameter values using an evolutionary algorithm.Evolutionary algorithms use mechanisms inspired by biological evolution (e.g., reproduction, mutation, recombination, etc.) and depend on the performance of the individual structures.Each individual in the population receives a measure of its fitness in the environment and selection focuses attention on high fitness individuals exploiting the available fitness information.On the other hand, the least square method is based on the minimization of the sum of squared residuals, with residuals being the differences between an observed data and model predicted results.In the case of non-linear regression, the Levenberg-Marquardt (LM) algorithm was applied.The Levenberg-Marquardt algorithm is an iterative technique that locates the minimum of the multivariate function that is expressed as a sum of squares of non-linear real-valued functions [26][27][28].The non-linear regression model coefficients describing the color components change during storage of the L1KMK coffee sample packed in tin can are presented in Tables 2 and 3. Results presented in Table 2 were obtained using evolutionary algorithm implemented into Eureqa Formulize.By comparison of the obtained R 2 values, it was noticed that the best agreement between experimental data and model was obtained for the a* color component with R 2 = 0.842.The smallest value of R 2 was obtained for the Hue angle R 2 = 0.488.Eureqa Formulize works on a concept of development of the best possible fit to describe the defined set of the experimental data.In most cases obtained models are over parameterized and applicable only on the tested data set.Due to that, in this work predefined expressions (non-linear and linear) were used to describe the color components change during storage.When using the Eureqa Formulize software it is not possible to get information about the error of the estimation of the parameters and their significance.To overcome the mentioned limitations Stataistica 10.0 was used.regression models were developed for each packaging material separately.Parameters of both nonlinear and linear regression models were estimated using two numerical methods, evolutionary algorithm implemented into the Eureqa Formulize software, and the least square algorithm implemented into the Statistica 10.0.software.The estimated values of model parameters are presented in Tables 2-5 and results of the dispersion of color components in comparison to model predictions are presented in Figures 1-4.
Eureqa Formulize estimated the model parameter values using an evolutionary algorithm.Evolutionary algorithms use mechanisms inspired by biological evolution (e.g., reproduction, mutation, recombination, etc.) and depend on the performance of the individual structures.Each individual in the population receives a measure of its fitness in the environment and selection focuses attention on high fitness individuals exploiting the available fitness information.On the other hand, the least square method is based on the minimization of the sum of squared residuals, with residuals being the differences between an observed data and model predicted results.In the case of non-linear regression, the Levenberg-Marquardt (LM) algorithm was applied.The Levenberg-Marquardt algorithm is an iterative technique that locates the minimum of the multivariate function that is expressed as a sum of squares of non-linear real-valued functions [26][27][28].The non-linear regression model coefficients describing the color components change during storage of the L1KMK coffee sample packed in tin can are presented in Tables 2 and 3. Results presented in Table 2 were obtained using evolutionary algorithm implemented into Eureqa Formulize.By comparison of the obtained R 2 values, it was noticed that the best agreement between experimental data and model was obtained for the a* color component with R 2 = 0.842.The smallest value of R 2 was obtained for the Hue angle R 2 = 0.488.Eureqa Formulize works on a concept of development of the best possible fit to describe the defined set of the experimental data.In most cases obtained models are over parameterized and applicable only on the tested data set.Due to that, in this work predefined expressions (non-linear and linear) were used to describe the color components change during storage.When using the Eureqa Formulize software it is not possible to get information about the error of the estimation of the parameters and their significance.To overcome the mentioned limitations Stataistica 10.0 was used.Parameters of the non-linear regression model for the color change during storage of L1KMK coffee sample packed in tin can were estimated using Statistica 10.0, are given in Table 3.By analyzing the obtained p-values, it can be seen that all estimated parameters were significant for models describing L*, b*, Chroma and Hue angle.In case of model describing a*, p-value indicated that the b 2 parameter value was not significant.This parameter was connected to the relative air humidity in the model.Interestingly, in the case of ∆E there was only one significant parameter, b 1 , connected to sample moisture content.By comparing the estimated parameters values obtained using two algorithms presented in Tables 2 and 3, it can be seen that the largest differences were obtained for the estimation of the b 0 parameter for all six color characteristics.
Results presented in Figures 1 and 2, showed that the dispersion of color components comparing to model predictions was quite uniform with only few outlayers, especially evident for the Hue angle description.Using Levenberg-Marquardt algorithm for model parameters estimation very good agreement between model and experimental data was obtained.The estimated R 2 values were R 2 (L*) = 0.918, R 2 (a*) = 0.936, R 2 (b*) = 0.783, R 2 (∆E) = 0.819, R 2 (Chroma) = 0.794, R 2 (Hue angle) = 0.843, what represented a significant improvement in comparison to those obtained using Eureqa Formulize (Table 2 and Figure 1).
The effect of storage conditions on the coffee colour components change was also analysed for the V1KMK coffee sample packed in the triplex bag.Parameters of the non-linear model were first estimated using the Eureqa Formulize.Results are presented in Table 4.The evolutionary algorithm used in the Eureqa Formulize was not able to estimate the optimal values of the parameters describing the change of the L* colour component, probably due to numerical instability.Regarding the R 2 value, the best agreement between experimental data and model was obtained for description of the change of the a* colour component with R 2 = 0.977, while the smallest value of R 2 was obtained for the Hue angle R 2 = 0.642 which was in analogy with the results obtained for the L1KMK sample (Table 2).Comparing the R 2 values for both packaging, it was noticed that the proposed non-linear models described colour components change in triplex packaging (V1KMK) much better.As for the tin can packaging (L1KMK), parameters of the non-linear models for description of the colour components change during storage in V1KMK were also estimated using Statistica 10.0.Results are given in Table 5.The number or significant parameters based on p-value is reduced compared to the models developed for the coffee sample packed in the tin can.In the case of the L* colour component there was no significant parameter, while in case of ∆E, only b 1 was shown to be significant.Only for description of the Hue angle change all parameters were significant.The proposed model described the experimental data very well (Figure 4), R 2 (a*) = 0.990, R 2 (b*) = 0.974, R 2 (∆E) = 0.918, R 2 (Chroma) = 0.978, R 2 (Hue angle) = 0.918.The smallest value of R 2 was obtained for the L* colour component with R 2 = 0.471.As mentioned before, the evolutionary algorithm was not able to estimate the parameters of the proposed model, so the results obtained by the Levenberg-Marquardt algorithm presented a significant improvement.The dispersion of colour components compared to model predictions (Figures 3 and 4) was quite uniform, except for the L* component visible in Figure 4.The number or significant parameters based on p-value is reduced compared to the models developed for the coffee sample packed in the tin can.In the case of the L* colour component there was no significant parameter, while in case of ΔE, only b1 was shown to be significant.Only for description of the Hue angle change all parameters were significant.The proposed model described the experimental data very well (Figure 4), R 2 (a*) = 0.990, R 2 (b*) = 0.974, R 2 (ΔE) = 0.918, R 2 (Chroma) = 0.978, R 2 (Hue angle) = 0.918.The smallest value of R 2 was obtained for the L* colour component with R 2 = 0.471.As mentioned before, the evolutionary algorithm was not able to estimate the parameters of the proposed model, so the results obtained by the Levenberg-Marquardt algorithm presented a significant improvement.The dispersion of colour components compared to model predictions (Figures 3 and 4) was quite uniform, except for the L* component visible in Figure 4.

Linear Regression Models
To test the relationship between storage conditions (sample moisture content, relative air humidity, and air temperature) and coffee colour components change during storage linear regression models were proposed and the logarithmic linearization was used (Equation ( 6)).Parameters of the models were estimated using the evolutionary algorithm and variance minimization (least square method) (Tables 6-9).Linear regression models were developed for both packaging materials separately.
Coefficients of linear regression models for description of colour parameters change during the storage of L1KMK estimated using Eureqa Formulize software are given in Table 6.The best agreement between model and experimental data was obtained for the a* component with R 2 = 0.856.The smallest value of R 2 was obtained for the b* component with R 2 = 0.400.In the case of non-linear regression models (Table 2), the best agreement between model and experiment was also obtained for a* colour component and the smallest value of R 2 was obtained for the Hue angle.As mentioned before, parameters of the linear regression model were also estimated using variance minimization method implemented into Statistica 10.0.Results are given in Table 7.
It can be noticed that for the change of the L*, b* and Chroma, all four parameters are significant (p < 0.05).In case of the a* component, all parameters except β 2 (combined with relative air humidity) are significant, while in the case of Hue angle, all parameters except β 3 (combined with air temperature) are significant.Interestingly, p-value analysis revealed that in case of ∆E only parameter β 3 proved to be significant, which meant that temperature was the key factor affecting total colour difference.Cardelli and Labudza [29] also showed that storage temperature affects shelf-life of roasted and ground coffee.The best agreement between model and experimental data regarding R 2 value was obtained for the a* colour component with R 2 = 0.937, while the smallest value of R 2 was obtained for the ∆E component, with R 2 = 0.457.As in the case of non-linear model for coffee sample L1KMK packed in tin can, models developed using Statistica 10.0 described experimental data more precisely than models developed using Eureqa Formulize.
precisely than models developed using Eureqa Formulize.
By analysing the results presented in Figures 5-6, it can be seen that the dispersion of observed colour component values compared to the model predictions followed a linear trend with few outliers.Outliers are mostly evident for regression models describing the change of a* and Chroma colour components estimated using the Eureqa Formulize, and for regression models describing ΔE and Chroma change estimated using the Statistica 10.0.
As in the case of non-linear regression models (Table 4), the Eureqa Formulize was not able to estimate the parameters of the linear regression model for description of L* colour component change during storage for coffee sample V1KMK in triplex packaging (Table 8).By analysing the results presented in Figures 5 and 6, it can be seen that the dispersion of observed colour component values compared to the model predictions followed a linear trend with few outliers.Outliers are mostly evident for regression models describing the change of a* and Chroma colour components estimated using the Eureqa Formulize, and for regression models describing ∆E and Chroma change estimated using the Statistica 10.0.
As in the case of non-linear regression models (Table 4), the Eureqa Formulize was not able to estimate the parameters of the linear regression model for description of L* colour component change during storage for coffee sample V1KMK in triplex packaging (Table 8).The best agreement between model and experimental data was obtained for the a* colour component with R 2 = 0.979 and the smallest R 2 value was obtained for Hue angle (R 2 = 0.786).Comparing the R 2 values of both non-linear and linear regression models for colour components change during storage of V1KMK sample, it can be noticed that non-linear models could be selected as more suitable for the colour change description.
Linear regression model parameters for V1KMK were also estimated using Statistica 10.0 (Table 9).All parameters describing the change of b*, ∆E and Chroma were significant.In case of the L* colour component all parameters except β 2 (connected to relative air humidity) were significant while in the case of the L* colour component parameters β 0 and β 2 were significant.Interestingly, p-value analysis revealed that in the case of the Hue angle there were no significant parameters, while in the case on non-linear regression model for description of the Hue angle change all parameters were significant.The phenomena could be explained by the fact that model linearization was performed using logarithm values of the experimental data.By introducing logarithmic values, disruption between orders of magnitude of the measured data can occur.The best agreement between model and experimental data was obtained for the a* component (R 2 = 0.984) and the smallest value of R 2 was obtained for the L* component (R 2 = 0.504).Dispersion of colour components comparing to model predictions followed a linear trend (Figures 7 and 8).The largest data dispersion was evident for models describing the change of the Hue angle estimated using the Eureqa Formulize and for models describing the L* component estimated using the Statistica 10.0.
Both non-linear and linear models used in this study, pointed to a significant influence of intrinsic (sample moisture content) and external (RH and temperature) factors mostly on L*, b*, Chroma and Hue values.Previous research has also shown a significant influence of storage conditions (e.g., temperature, headspace volume, water activity) on coffee quality [25,29].Furthermore, it is important to emphasize that, in home storage conditions, models should be created for all influences combined together, rather than separately.The nature of the relationship between variables for description of the colour change during storage is non-linear, so non-linear model were the most suitable for process description.Further, non-linear model used in this study uses raw experimental data for parameter prediction, while the linear one introduces ln value to calculation, which causes disturbance between orders of magnitudes of measured data.On the other hand, result of the linear form of the model can usually be much easier to interpret.Based on the obtained result it can be concluded that due to its numerical stability Statistica 10.0 was a better choice for estimation of the parameter of both non-linear and linear models.The difference between packaging materials could only be estimated based on sample moisture content-the coffee sample packed in tin can (L1KMK) absorbed more moisture due to a larger amount of air entrapped in the tin can.Since moisture is well known to cause a deterioration in food powder properties [30], triplex bag is proposed as the more suitable packaging type for roasted ground coffee.

Conclusions
Both non-linear and linear models used in this study, pointed to a significant influence of intrinsic (sample moisture content) and external (RH and temperature) factors on ground roasted coffee color change.Based on lower moisture sorption of the sample packed in triplex bag, triplex packaging is proposed as more suitable.The nature of the relationship between variables for description of the color change during storage is non-linear, so non-linear model were the most suitable for process description.Due to numerical stability, Statistica 10.0 was a better choice for estimation of parameter of both non-linear and linear models.

Conclusions
Both non-linear and linear models used in this study, pointed to a significant influence of intrinsic (sample moisture content) and external (RH and temperature) factors on ground roasted coffee color change.Based on lower moisture sorption of the sample packed in triplex bag, triplex packaging is proposed as more suitable.The nature of the relationship between variables for description of the color change during storage is non-linear, so non-linear model were the most suitable for process description.Due to numerical stability, Statistica 10.0 was a better choice for estimation of parameter of both non-linear and linear models.

Conclusions
Both non-linear and linear models used in this study, pointed to a significant influence of intrinsic (sample moisture content) and external (RH and temperature) factors on ground roasted coffee color change.Based on lower moisture sorption of the sample packed in triplex bag, triplex packaging is proposed as more suitable.The nature of the relationship between variables for description of the color change during storage is non-linear, so non-linear model were the most suitable for process description.Due to numerical stability, Statistica 10.0 was a better choice for estimation of parameter of both non-linear and linear models.

Table 1 .
Color components of roasted ground coffee samples during 6 months of storage.Results are shown as average ± SD.

Table 2 .
Non-linear regression models coefficients for description of color component change during storage for L1KMK (α = 95%) estimated using Eureqa Formulize software.

Table 2 .
Non-linear regression models coefficients for description of color component change during storage for L1KMK (α = 95%) estimated using Eureqa Formulize software.

Table 3 .
Non-linear regression models coefficients for description of color component change during storage for tin can packaging (L1KMK) (α = 95%) estimated using Statistica 10.0.† marked values are significant at p < 0.05.

Table 5 .
Non-linear regression models coefficients for description of colour component change during storage for V1KMK (α = 95%) estimated using Statistica 10.0.† marked values are significant at p < 0.05.

Table 5 .
Non-linear regression models coefficients for description of colour component change during storage for V1KMK (α = 95%) estimated using Statistica 10.0.† marked values are significant at p < 0.05.

Table 6 .
Linear regression models coefficients for description of colour components change during storage of L1KMK (α = 95%) estimated using the Eureqa Formulize software.

Table 8 .
Linear regression models coefficients for description of colour component change during storage of V1KMK (α = 95%) estimated using the Eureqa Formulize.

Table 7 .
Linear regression models coefficients for description of colour component change during the storage of L1KMK (α = 95%) estimated using Statistica 10.0.† significant at p < 0.05.

Table 8 .
Linear regression models coefficients for description of colour component change during storage of V1KMK (α = 95%) estimated using the Eureqa Formulize.

Table 9 .
Linear regression models coefficients for description of colour components change during storage of V1KMK (α = 95%) estimated using Statistica 10.0.† marked values are significant at p < 0.05.