**Open Journal of Statistics**

Vol.06 No.02(2016), Article ID:65713,10 pages

10.4236/ojs.2016.62020

Off-Line Quality Control in Bioproducts Processing: A Case Study in Greece

Ioannis S. Triantafyllou

Department of Statistics and Insurance Science, University of Piraeus, Piraeus, Greece

Copyright © 2016 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 22 January 2016; accepted 18 April 2016; published 21 April 2016

ABSTRACT

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.

**Keywords:**

Grape Molasses, Bioproduct, Statistical Quality Control, Analysis of Variance

1. Introduction

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 [1] or to the more recent contributions of [2] - [5] or the books [6] and [7] .

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 [8] - [10] . Throughout this paper, a laboratory located in Peloponnese of Greece, which produces grape mo- lasses by using 100% biological raw material, is under investigation. The aforementioned laboratory, named “KARAGGELIS OE” is involved in foods production since 2011. The company incorporates a production unit of biological grape molasses and carries out products delivery in Greece and other countries. Generally speaking, petimezi is undiluted must, produced by boiled-down grape juice and it is rich in energy, iron and calcium. It is used self-same in cooking and pastry and constitutes an ersatz of sugar or honey. Since bioproducts are not much prevalent in

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.

2. Methods

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. [11] ), or linear regression analysis have been applied in order to detect effectively those factors that control noise or target of the bioproducts processing.

3. Main Results

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

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. [12] ).

Let us next denote by Y_{1}, Y_{2} the Brix value of must and grape molasses respectively, while Y_{3}, Y_{4} express the volatile acidity of each of them. The main results section is organized as follows. At the beginning, the Brix value of must (Y_{1}-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 Y_{1}-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

(1)

where μ is the mean of Y_{1}-values, a_{i} (β_{j}) represents the main effect of the i-th (j-th) level of factor A (B), indicates the interaction term between the i-th level of factor A and the j-th level of factor B, while e_{ijk} is the corresponding estimated residual. The results are summarized in Table 1.

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 Table 2) that

・ 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 Table 3),

・ 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 Y_{1}-values. According to Taguchi guidelines (see, e.g. [11] or [13] ), the appropriate target measure is the average Brix value of must, while for the noise of the production, since the Brix value follows the rule “the larger the better”, the proposed measure, named theta, is given as follows

. (2)

Applying the appropriate ANOVA models for the aforementioned measures, the results are recapitulated in Table 4 and Table 5.

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. [13] )

, (3)

where, s is the sample mean and standard deviation of responses respectively, while the parameter b is related to the possible existence of a functional connection between the population mean (μ) and standard deviation (σ). In other words, if a relation of the form can be statistically affirmed, then the parameter b of

Table 1. ANOVA model: Brix value of must versus ground and variety.

Table 2. Tukey pairwise comparisons: Brix value of must in terms of grapes variety.

Table 3. Tukey pairwise comparisons: Brix value of must in terms of ground.

Table 4. ANOVA model: Theta measure of Brix value of must versus ground and variety.

Table 5. ANOVA model: Target measure of Brix value of must versus ground and variety.

the above formula should be estimated via an appropriate regression analysis. Otherwise, the parameter b shall be replaced by zero. If we denote by, s_{i} the mean and standard deviation of responses in the i-th trial, the outcomes of the regression analysis between them, are summarized in Table 6.

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 (Y_{2}-value) or the volatile acidity of must (Y_{3}-value) and petimezi (Y_{4}-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 Table 7).

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 Table 8), we observe the following

・ 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 Table 9).

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). Table 10 and Table 11 reveal that factors A and B can be characterized only as controllers of target value while noise of the production seems not to be controlled by them.

Table 6. Regression model for Brix value of must: Logarithm of standard deviation versus logarithm of mean.

Table 7. ANOVA model: Brix value of grape molasses versus ground and variety.

Table 8. Tukey pairwise comparisons: Brix value of grape molasses in terms of grapes variety.

Table 9. Tukey pairwise comparisons: Brix value of grape molasses in terms of ground.

Table 10. ANOVA model: Noise measure of Brix value of petimezi versus ground and variety.

Table 11. ANOVA model: Target measure of Brix value of petimezi versus ground and variety.

On the other hand, the volatile acidity of must (Y_{3}-values) is not affected by the ground that has been used for planting, since the corresponding observed significance is equal to 0.979 (see Table 12).

Grapes variety seems to affect significantly the response Y_{3} (p = 0.002) of the experiment. Based on multiple comparisons between levels of factor A displayed in Table 13, one may readily deduce the following

・ 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

(4)

The analysis of the above measures, leads to the conclusion that grapes variety seems to control both target and noise (see Table 14 and Table 15).

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, Table 16).

Therefore the proposed measure m is given as

(5)

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 Y_{3} responses (see Table 17).

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 Table 18).

On the other hand, grapes variety seems to affect significantly the Y_{4} response (p = 0.001) of the experiment. Based on multiple comparisons between levels of factor A displayed in Table 19, one may readily deduce the following

Table 12. ANOVA model: Volatile acidity of must versus ground and variety.

Table 13. Tukey pairwise comparisons: Volatile acidity of must in terms of grapes variety.

Table 14. ANOVA model: Noise measure h of volatile acidity of must versus ground and variety.

Table 15. ANOVA model: Target measure of volatile acidity of must versus ground and variety.

Table 16. Regression model for volatile acidity of must: Logarithm of standard deviation versus logarithm of mean.

Table 17. ANOVA model: Noise measure m of volatile acidity of must versus ground and variety.

Table 18. ANOVA model: Volatile acidity of petimezi versus ground and variety.

Table 19. Tukey pairwise comparisons: Volatile acidity of petimezi in terms of grapes variety.

・ 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, Table 20).

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 Y_{4} (see Table 21).

4. Conclusion

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

Table 20. Regression model for volatile acidity of petimezi: Logarithm of standard deviation versus logarithm of mean.

Table 21. ANOVA model: Noise measure of volatile acidity of petimezi versus ground and variety.

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.

Acknowledgements

The author wish to thank the anonymous referees for their useful comments and suggestions which improved the presentation of the paper.

Cite this 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

References

- 1. Shewhart, W.A. (1926) Quality Control Charts. Bell System Technical Journal, 2, 593-603.

http://dx.doi.org/10.1002/j.1538-7305.1926.tb00125.x - 2. Balakrishnan, N., Triantafyllou, I.S. and Koutras, M.V. (2009) Nonparametric Control Charts Based on Runs and Wilcoxon-Type Rank-Sum Statistics. Journal of Statistical Planning and Inference, 139, 3177-3192.

http://dx.doi.org/10.1016/j.jspi.2009.02.013 - 3. Balakrishnan, N., Triantafyllou, I.S. and Koutras, M.V. (2010) A Distribution-Free Control Chart Based on Order Statistics. Communication in Statistics—Theory and Methods, 39, 3652-3677.

http://dx.doi.org/10.1080/03610920903324858 - 4. Chakraborti, S. (2011) Nonparametric (Distribution-Free) Quality Control Charts. Encyclopedia of Statistical Sciences, 1-27.

http://dx.doi.org/10.1002/0471667196.ess7150 - 5. Chowdhury, S., Mukherjee, A. and Chakraborti, S. (2015) Distribution-Free Phase II CUSUM Control Chart for Joint Monitoring of Location and Scale. Quality and Reliability Engineering International, 31, 135-151.

http://dx.doi.org/10.1002/qre.1677 - 6. Gibbons, J.D. and Chakraborti, S. (2010) Nonparametric Statistical Inference. 5th Edition, Marcel Dekker, New York.
- 7. Balakrishnan, N. and Ng, H.K.T. (2006) Precedence-Type Tests and Applications. John Wiley & Sons, Hoboken.

http://dx.doi.org/10.1002/0470037849 - 8. Toker, O.S. and Dogan, M. (2013) Effect of Temperature and Starch Concentration on the Creep/Recovery Behavior of the Grape Molasses: Modelling with ANN, ANFIS and Response Surface Methodology. European Food Research and Technology, 236, 1049-1061.

http://dx.doi.org/10.1007/s00217-013-1959-0 - 9. Goksel, M., Dogan, M., Toker, O.S., Ozgen, S., Sarioglu, K. and Oral, R.A. (2013) The Effect of Starch Concentration and Temperature on Grape Molasses: Rheological and Textural Properties. Food and Bioprocess Technology, 6, 259-271.

http://dx.doi.org/10.1007/s11947-011-0705-5 - 10. Erdogrul, O. (2008) Pesticide Residues in Liquid Pekmez (Grape Molasses). Environmental Monitoring and Assessment, 144, 323-328.

http://dx.doi.org/10.1007/s10661-007-9995-5 - 11. Taguchi, G. (1986) Introduction to Quality Engineering: Designing Quality into Products and Processes. Asian Productivity Organization, Tokyo.
- 12. Montgomery, D.C. (2012) Design and Analysis of Experiments. 8th Edition, John Wiley & Sons, New York.
- 13. Logothetis, N. and Wynn, H.P. (1989) Quality through Design: Experimental Design, Off-Line Quality Control and Taguchi’s Contributions. Clarendon Press, Oxford.