In the present article we study the production of grape molasses. Data drawn from a specified biolaboratory, are properly analyzed in order to detect factors that affect significantly the Brix value and the volatile acidity of the final product. The ground that is used for planting and a variety of grapes have been taken into account. Off-line statistical quality control techniques have been employed and the outcomes are displayed and discussed in detail.
Statistical quality control is a useful tool for ensuring that the quality characteristics of a product are at the nominal or required levels and consists of a set of statistical methods for analyzing data. The overall group of the abovementioned techniques is divided into three main subsets, namely Design of Experiments (off-line techni- ques), Statistical Process Control (on-line monitoring) and Acceptance Sampling. For more details for the on- line monitoring of a productive process, the interested reader is referred to the early work [
The present study focuses on the quality control of the production process of a biological product called grape molasses (or petimezi). This product is made by 100% biological grapes, without using any pesticides. Grape molasses has attracted some research interest in the last decade; for example one may refer to publications [
The present study focuses on the quality improvement of the final bioproduct by examining the effect of some qualitative characteristics to volatile acidity and Brix value of grape molasses. In Section 2, the statistical methodology that is applied throughout this study, is presented in some detail. More specifically, statistical quality control techniques and analysis of variance models are applied in order to investigate the quality of petimezi produced in the aforementioned laboratory. Factors such as the variety of grapes used or the kind of ground, where the cultivation has taken place, are considered and properly analyzed. The results of the analysis are presented comprehensively in Section 3. Finally, Section 4 offers a summary of the conclusions extracted by statistical analysis and formulizes several suggestions that may lead to the quality improvement of the final product.
The study used data drawn during 2014 from the biolaboratory of the company “KARAGGELIS OE”. The response variables of the study are determined to be the Brix value and the volatile acidity of must and grape molasses. Each of the aforementioned quantitative characteristics has been examined whether is strongly affected or not, by a set of qualitative factors such as the variety of grapes and the type of ground that have been used. Analysis of Variance models have been appropriately applied in order to shed light on the above subject. Furthermore, multiple comparisons among the available levels of each factor have been presented by using the Tukey method. In addition, off-line quality control tools, such as Taguchi measures (see, e.g. [
In this section, we shall employ appropriate statistical methodologies in order to analyze data collected from the ad hoc laboratory, mentioned in the previous section, where the bioproduct of grape molasses is brought forth. The main objective of the present data analysis is to detect qualitative factors that have important effect on the volatile acidity and Brix value of the bioproduct of petimezi. Moreover, methods of Statistical Quality Control are also applied in order to determine the appropriate measures of noise and target, while based on the proposed measures, we shall appraise those factors that contribute significantly to the variability and central tendency of the volatile acidity and Brix value of the final bioproduct.
The laboratory that is under our study, is located in continental Greece and more specifically in Peloponnese . The grape molasses that is turned out, is affected basically by two factors; the first one refers to the variety of grape that is planted (factor A, hereafter). The second factor is associated with the ground, where the cultivation is taken place (factor B, hereafter). In the line of grape molasses production, the tilled ground is a three-level factor (lowland, semi-mountain, mountain), while the variety of grape is characterized by five levels (Merlot, Agiorgitiko, Camborne, Corinthian, Phocean).
Data have been appropriately collected so that at least two repetitions for each of the fifteen possible therapies have been accomplished. The number of available repetitions for each therapy, assures that the variance within each therapy could be estimated (see, e.g. [
Let us next denote by Y1, Y2 the Brix value of must and grape molasses respectively, while Y3, Y4 express the volatile acidity of each of them. The main results section is organized as follows. At the beginning, the Brix value of must (Y1-value) is under investigation. Analysis of variance (ANOVA) is carried out and the corresponding results are discussed in detail. In addition, using the appropriate measures, we detect significant factors that control the variability and the central tendency of Y1-value in the bioproductive process. Similar outcomes are also reported for the rest characteristics, e.g. the Brix value of grape molasses and the volatile acidity of must and petimezi.
We next apply ANOVA model in order to divide the total variability of Brix characteristic of must, into two parts; the first one is related to measurable and manageable variables (see, e.g. factor A or B), while the second one comes unmanageable or from random causes.
The factorial experiment that has been taken place, involves two factors with fixed effects (factors A and B) and at least two repetitions for each therapy. The full-factorial model that has been initially applied is given as follows
where μ is the mean of Y1-values, ai (βj) represents the main effect of the i-th (j-th) level of factor A (B),
One may easily conclude that the main effects of factors A, B and the respective interaction term are all significant at significance level 5% (p < 0.001). In words, the Brix value of must that accompanies the bioproduct, is strongly affected by the specific variety of grapes that is planted and the kind of ground that is used. In addition, the effect of variety of grapes on Brix value of must, is strongly differentiated by the choice of the ground and reversely. In order to shed light on the pairwise comparisons between the mean values of different levels of each factor, the following Tukey simultaneous 95% confidence intervals are constructed.
Since the larger is the Brix value, the better is the must, we may readily deduce at significance level 5% (see
・ Camborne variety gains the better of Phocean, Merlot and Agiorgitiko variety.
・ Corinthian variety gains the better of Phocean, Merlot, Camborne and Agiorgitiko variety.
In a similar way, we may conclude that at significance level 5% (see
・ mountain ground gains the better of lowland and semi-mountain ground.
・ semi-mountain ground gains the better of lowland ground.
In the sequel, we shall investigate whether factors A and B can be characterized as control factors of noise or control factors of target of Y1-values. According to Taguchi guidelines (see, e.g. [
Applying the appropriate ANOVA models for the aforementioned measures, the results are recapitulated in
Based on the above tables, we deduce that both factors A and B are characterized as control factors of noise and control factors of target. In order to bring out this confounding, an alternative measure for noise production should be used. More specifically, the proposed measure is given as follows (see, e.g. [
where
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 93.824 | 23.456 | 294.04 | <0.001 |
B | 2 | 88.542 | 44.271 | 554.98 | <0.001 |
AB | 8 | 20.384 | 2.548 | 31.94 | <0.001 |
Error | 29 | 2.313 | 0.080 | ||
Total | 43 | 204.068 |
Difference between levels of factor A | Difference of means | Std. error of difference | Simultaneous 95% confidence intervals | T-value | p-value |
---|---|---|---|---|---|
Merlot-Phocean | −0.164 | 0.136 | (−0.559, 0.231) | −1.21 | 0.748 |
Agiorgitiko-Phocean | 0.022 | 0.130 | (−0.357, 0.401) | 0.17 | 1 |
Corinthian-Phocean | 4.453 | 0.161 | (3.986, 4.920) | 27.69 | <0.001 |
Camborne-Phocean | 1.797 | 0.139 | (1.394, 2.200) | 12.97 | <0.001 |
Agiorgitiko-Merlot | 0.186 | 0.127 | (−0.184, 0.557) | 1.46 | 0.595 |
Corinthian-Merlot | 4.617 | 0.158 | (4.156, 5.077) | 29.13 | <0.001 |
Camborne-Merlot | 1.961 | 0.136 | (1.566, 2.356) | 14.43 | <0.001 |
Corinthian-Agiorgitiko | 4.431 | 0.154 | (3.984, 4.877) | 28.82 | <0.001 |
Camborne-Agiorgitiko | 1.775 | 0.130 | (1.396, 2.154) | 13.62 | <0.001 |
Camborne-Corinthian | −2.656 | 0.161 | (−3.123, −2.188) | −16.52 | <0.001 |
Difference between levels of factor B | Difference of means | Std. Error of Difference | Simultaneous 95% Confidence Intervals | T-value | p-value |
---|---|---|---|---|---|
Semi-mountain-Lowland | 1.891 | 0.104 | (1.634, 2.149) | 18.12 | <0.001 |
Mountain-Lowland | 3.825 | 0.115 | (3.540, 4.110) | 33.17 | <0.001 |
Mountain-Semi-mountain | 1.933 | 0.114 | (1.652, 2.215) | 16.94 | <0.001 |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 5.0522 | 1.2631 | 11.29 | 0.002 |
B | 2 | 4.3907 | 2.1953 | 19.62 | 0.001 |
Error | 8 | 0.8950 | 0.1119 | ||
Total | 14 |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 40.375 | 10.094 | 12.08 | 0.002 |
B | 2 | 32.694 | 16.347 | 19.57 | 0.001 |
Error | 8 | 6.682 | 0.835 | ||
Total | 14 |
the above formula should be estimated via an appropriate regression analysis. Otherwise, the parameter b shall be replaced by zero. If we denote by
Among others, we conclude that there exists a statistically significant linear relation between the logarithms of the aforementioned quantities (p = 0.021) while the estimation of parameter b is equal to 4.18. We next apply Equation (3) in order to construct the alternative measure of noise by replacing properly the parameter b. The analysis of this measure revealed that neither variety of grapes nor the ground could be characterized as control factors of noise anymore. Therefore, our final deduction is that we should adjust properly the factors A and B in order to control the target value of Brix measurements of must.
Employing analogous arguments, one may draw similar conclusions for the remaining quantitative responses of the experiment, such as the Brix value of grape molasses (Y2-value) or the volatile acidity of must (Y3-value) and petimezi (Y4-value). As far it concerns with the Brix value of grape molasses, the ANOVA disclosed that both factors A and B affect significantly the response of each trial (p < 0.001, see
However, there seems to be a negligible interaction between them (p = 0.093). Taking a look at the pairwise comparisons between levels of factor A (see
・ Camborne variety gains the better of Phocean, Agiorgitiko and Merlot variety.
・ Corinthian variety gains the better of Phocean, Merlot and Agiorgitiko variety.
In a similar way, we may conclude that at significance level 5% mountain ground gains the better of lowland and semi-mountain ground (see
According to Taguchi guidelines, the appropriate target measure is the average Brix value of petimezi, while for the noise of production is given in Equation (3).
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
Logarithm of mean | 1 | 0.45518 | 0.45518 | 6.92 | 0.021 |
Error | 13 | 0.85485 | 0.06576 | ||
Regression equation |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 2.036 | 0.5090 | 13.22 | <0.001 |
B | 2 | 1.034 | 0.5170 | 13.43 | <0.001 |
AB | 8 | 0.5946 | 0.0743 | 1.93 | 0.093 |
Error | 29 | 1.1167 | 0.0385 | ||
Total | 43 | 5.0480 |
Difference between levels of factor A | Difference of means | Std. error of difference | Simultaneous 95% confidence intervals | T-value | p-value |
---|---|---|---|---|---|
Merlot-Phocean | 0.0556 | 0.0944 | (−0.2188, 0.3299) | 0.59 | 0.976 |
Agiorgitiko-Phocean | 0.15 | 0.0906 | (−0.1132, 0.4132) | 1.66 | 0.475 |
Corinthian-Phocean | 0.667 | 0.112 | (0.342, 0.991) | 5.97 | <0.001 |
Camborne-Phocean | 0.4278 | 0.0963 | (0.148, 0.7076) | 4.44 | 0.001 |
Agiorgitiko-Merlot | 0.0944 | 0.0886 | (−0.1629, 0.3518) | 1.07 | 0.822 |
Corinthian-Merlot | 0.611 | 0.11 | (0.291, 0.931) | 5.55 | <0.001 |
Camborne-Merlot | 0.3722 | 0.0944 | (0.0978, 0.6466) | 3.94 | 0.004 |
Corinthian-Agiorgitiko | 0.517 | 0.107 | (0.206, 0.827) | 4.84 | <0.001 |
Camborne-Agiorgitiko | 0.2778 | 0.0906 | (0.0146, 0.5410) | 3.07 | 0.035 |
Camborne-Corinthian | −0.239 | 0.112 | (−0.564, 0.086) | −2.14 | 0.232 |
Difference between levels of factor A | Difference of means | Std. error of difference | Simultaneous 95% confidence intervals | T-value | p-value |
---|---|---|---|---|---|
Semi-mountain-Lowland | 0.12 | 0.0725 | (−0.059, 0.299) | 1.65 | 0.24 |
Mountain-Lowland | 0.41 | 0.0801 | (0.2123, 0.6077) | 5.12 | <0.001 |
Mountain-Semi-mountain | 0.29 | 0.0793 | (0.0943, 0.4857) | 3.66 | 0.003 |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 11.172 | 2.793 | 2.23 | 0.156 |
B | 2 | 0.424 | 0.212 | 0.17 | 0.847 |
Error | 8 | 10.038 | 1.255 | ||
Total | 14 | 21.633 |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 0.8184 | 0.2046 | 8.79 | 0.005 |
B | 2 | 0.3841 | 0.19206 | 8.25 | 0.011 |
Error | 8 | 0.1863 | 0.02328 | ||
Total | 14 | 1.3888 |
On the other hand, the volatile acidity of must (Y3-values) is not affected by the ground that has been used for planting, since the corresponding observed significance is equal to 0.979 (see
Grapes variety seems to affect significantly the response Y3 (p = 0.002) of the experiment. Based on multiple comparisons between levels of factor A displayed in
・ Phocean grapes variety gains the better of Agiorgitiko and Camborne grapes variety.
・ Corinthian grapes variety gains the better of Camborne and Agiorgitiko grapes variety.
According to Taguchi guidelines, the appropriate target measure is the average volatile acidity of must, while for the noise of measurements, based on the rule “the smaller the better”, is given as
The analysis of the above measures, leads to the conclusion that grapes variety seems to control both target and noise (see
Following similar methodology as before, we construct an alternative measure for noise (see Equation (3)), where the parameter b has been proved non-significant (p = 0.28,
Therefore the proposed measure m is given as
Based on the above measure, factors A and B do not seem to control noise, while factor A remains as the one that controls central tendency of the Y3 responses (see
Finally, the volatile acidity of petimezi is not affected by the ground that has been used for planting, since the corresponding observed significance is equal to 0.958 (see
On the other hand, grapes variety seems to affect significantly the Y4 response (p = 0.001) of the experiment. Based on multiple comparisons between levels of factor A displayed in
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 0.127454 | 0.031864 | 5.54 | 0.002 |
B | 2 | 0.000248 | 0.000124 | 0.02 | 0.979 |
AB | 8 | 0.076828 | 0.009604 | 1.67 | 0.149 |
Error | 29 | 0.166758 | 0.00575 | ||
Total | 43 | 0.374891 |
Difference between levels of factor A | Difference of means | Std. error of difference | Simultaneous 95% confidence intervals | T-value | p-value |
---|---|---|---|---|---|
Merlot-Phocean | 0.0436 | 0.0365 | (−0.0624, 0.1496) | 1.2 | 0.754 |
Agiorgitiko-Phocean | 0.1094 | 0.035 | (0.0077, 0.2111) | 3.13 | 0.030 |
Corinthian-Phocean | −0.0381 | 0.0432 | (−0.1635, 0.0874) | −0.88 | 0.901 |
Camborne-Phocean | 0.1117 | 0.0372 | (0.0035, 0.2198) | 3.00 | 0.040 |
Agiorgitiko-Merlot | 0.0658 | 0.0342 | (−0.0336, 0.1653) | 1.92 | 0.328 |
Corinthian-Merlot | −0.0817 | 0.0425 | (−0.2053, 0.042) | −1.92 | 0.330 |
Camborne-Merlot | 0.0681 | 0.0365 | (−0.0380, 0.1741) | 1.87 | 0.358 |
Corinthian-Agiorgitiko | −0.1475 | 0.0413 | (−0.2675, −0.0275) | −3.57 | 0.010 |
Camborne-Agiorgitiko | 0.0022 | 0.0350 | (−0.0995, 0.1039) | 0.06 | 1 |
Camborne-Corinthian | 0.1497 | 0.0432 | (0.0243, 0.2752) | 3.47 | 0.013 |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 50.462 | 12.6155 | 4.16 | 0.041 |
B | 2 | 1.922 | 0.9610 | 0.32 | 0.737 |
Error | 8 | 24.246 | 3.0308 | ||
Total | 14 | 76.630 |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 0.04584 | 0.01146 | 3.83 | 0.050 |
B | 2 | 0.00036 | 0.00018 | 0.06 | 0.942 |
Error | 8 | 0.02395 | 0.00299 | ||
Total | 14 | 0.07015 |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
Logarithm of mean | 1 | 0.07986 | 0.07986 | 1.27 | 0.280 |
Error | 13 | 0.81519 | 0.06271 | ||
Regression equation |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 90.12 | 22.53 | 1.38 | 0.323 |
B | 2 | 13.7.19 | 68.60 | 4.20 | 0.057 |
Error | 8 | 130.70 | 16.34 | ||
Total | 14 | 358.02 |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 0.108403 | 0.027101 | 6.10 | 0.001 |
B | 2 | 0.000378 | 0.000189 | 0.04 | 0.958 |
AB | 8 | 0.063741 | 0.007968 | 1.79 | 0.119 |
Error | 29 | 0.128858 | 0.004443 | ||
Total | 43 | 0.304173 |
Difference between levels of factor A | Difference of means | Std. error of difference | Simultaneous 95% confidence intervals | T-value | p-value |
---|---|---|---|---|---|
Merlot-Phocean | 0.0353 | 0.0321 | (−0.0579, 0.1285) | 1.10 | 0.805 |
Agiorgitiko-Phocean | 0.0992 | 0.0308 | (0.0098, 0.1886) | 3.22 | 0.024 |
Corinthian-Phocean | −0.0353 | 0.0379 | (−0.1456, 0.075) | −0.93 | 0.883 |
Camborne-Phocean | 0.1036 | 0.0327 | (0.0086, 0.1987) | 3.17 | 0.027 |
Agiorgitiko-Merlot | 0.0639 | 0.0301 | (−0.0235, 0.1513) | 2.12 | 0.238 |
Corinthian-Merlot | −0.0706 | 0.0374 | (−0.1793, 0.0381) | −1.89 | 0.347 |
Camborne-Merlot | 0.0683 | 0.0321 | (−0.0249, 0.1615) | 2.13 | 0.235 |
Corinthian-Agiorgitiko | −0.1344 | 0.0363 | (−0.2399, −0.0290) | −3.71 | 0.007 |
Camborne-Agiorgitiko | 0.0044 | 0.0308 | (−0.0850, 0.0938) | 0.14 | 1 |
Camborne-Corinthian | 0.1389 | 0.0379 | (0.0286, 0.2492) | 3.66 | 0.008 |
・ Phocean grapes variety gains the better of Agiorgitiko and Camborne grapes variety.
・ Corinthian grapes variety gains the better of Camborne and Agiorgitiko grapes variety.
The analysis of measures proposed by Taguchi, leads to the conclusion that grapes variety seems to control both target and noise. Following similar methodology as before, we construct an alternative measure for noise (see Equation (3)), where the parameter b is non-significant (p = 0.354,
Therefore the proposed measure takes on the form presented in Equation (5). Based on the aforementioned measure, there is no factor that seems to control noise, while factor A remains as the one that controls central tendency of the response Y4 (see
In the present article, we provide an off-line quality control analysis for the bioproduction of grapes molasses. Data are drawn from an ad hoc biolaboratory located in Greece. As response variables have been defined two quantitative characteristics of grapes molasses, called Brix value and volatile acidity of bioproduct. A variety
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
Logarithm of mean | 1 | 0.06647 | 0.06647 | 0.92 | 0.354 |
Error | 13 | 0.93693 | 0.07207 | ||
Regression equation |
Source | Degrees of freedom | Sum of squares | Mean squares | F-value | p-value |
---|---|---|---|---|---|
A | 4 | 72.44 | 18.11 | 0.71 | 0.609 |
B | 2 | 124.13 | 62.06 | 2.42 | 0.150 |
Error | 8 | 204.79 | 25.60 | ||
Total | 14 | 401.36 |
of grapes that are planted and the ground where the cultivation is taken place, have been proved significant factors for the quality of the final product. Corinthian and Camborne varieties of grapes seem to lead to the optimum result as far it concerns the Brix value of must and grape molasses, while Phocean and Corinthian varieties of grapes are the best choices in order to abate their volatile acidity. On the other hand, mountain ground has been proved as the one that is capable to optimize Brix value of must and grape molasses. It is of some interest for future research, to depict and analyze additional characteristics of grapes molasses in order to enhance the quality of the bioproductive process of it.
The author wish to thank the anonymous referees for their useful comments and suggestions which improved the presentation of the paper.
Ioannis S. Triantafyllou, (2016) Off-Line Quality Control in Bioproducts Processing: A Case Study in Greece. Open Journal of Statistics,06,229-238. doi: 10.4236/ojs.2016.62020