Introduction n many manufacturing firms where there exists mass production, measurement made on each product is subject to error due to variation from one item to the other. Since there must be variations, it becomes important to study and determine when any observed variation is significant or not. This is the reason why the Federal Government of Nigeria came up with legislations to protect the buyers from buying inferior goods. Increase in consumer buying behaviou towards some selected drinks will directly affect the production of such drinks in our breweries industry. Quality control relies partly upon patronage and some other reliable factors, in beer production process, the measurement of attributes such as fill height and level of co 2 is of paramount important and that is the reason why quality control is evolving in developing systems to ensure standard products or services as well as meeting or exceeding customer's requirements. Walter Shewhart introduced the concept of statistical quality control thereby controlling quality of mass produced goods. Shewhart believed that variation always exists in manufactured products and that the variation can be studied, monitored and controlled using Statistics. Walter Shew hart explained the theories about using statistical quality control charts to improve quality and productivity in which case he developed fourte en points agenda for companies to improve quality and productivity, reduce costs and compete effectively in the world market. # II. # Literature Review Reeves and Bednar (1994) define quality as excellence, value, conformance to specifications, and meeting or exceeding customers' expectation. The term "fitness for use" defined by Juran (1974) is also included in the quality definition presented by Reeves and Bednar (1994). Thus, the customer perspective with respect to quality is the master key that should be understood while determining any term for quality or definition of quality. Deming, W.E (1986). worked on Quality and Productivity Improvement using acceptance sampling method, and he was able to obtain increase in quality and simultaneous reduction in the cost of reducing waste, re write staff attrition and litigation while increasing customer's loyalty. Farhat, B. A. and Al-Darrab, I. (1998). Total quality management is now established and widely used management process. One of its associated features is the application of statistical quality control techniques. A quality product or service is one that meets the customer's needs and provides the value that they want and expect. They are also of the opinion that quality management is a formal approach to management in which the overriding priority of the organization is to deliver a quality product or service and to work towards excellence and continuous improvement in everything it does. Quality can be viewed from the perspectives of design and product in which case; design quality is the different grades or levels of performance, reliability, I serviceability and function that are the results of deliberate engineering and management decision. On the other hand, product quality is the conformance of the product with specifications or expectations of the user in terms of fitness for use and cost. They are also of the opinion that control charts are closely related with statistical test of hypothesis. The control chart is a test of hypothesis that the process is in a state of statistical control. Shres tha and Chalidabhongse (2006) explained over their survey on 300 employees working in 60 Thai companies to what extent job satisfaction is affected by the existing performance appraisal system used by these companies. They concluded that since the performance appraisal system is part of the company's running processes, employees would show lower performance level if the appraisal system is not satisfactorily controlled. Cooper (2008) have emphasized on the impact of TQM practices on job satisfactions. The main aim of their research was to examine the relationship between people-related TQM practices and job satisfaction of service employees. The study triggers the question whether a TQM has an effect on employees' satisfaction. Pitterman (2000)'s findings on Telecordia technologies showed that customer satisfaction figures had gone up from 60% in 1992 to 95% at the time of implementing ISO 9001 quality system. Also, there was a 63% reduction noted in test cost efficiency since 1993 that 98% of major software released by Telcordia between 1995 and 1998 were delivered in time, even though the number of releases had tripled during the four year time. Takala et al. (2006) have gone even further to seeking customer satisfaction by improving and ensuring that customer satisfaction survey is supposed to be well designed and validated in order to be an effective measurement tool for its intended purpose. In their research paper, the purpose was to verify the reliability of customer satisfaction survey in context to three aspects of service; quality, delivery and responsiveness. They concluded that there was a need to work on the flexibility of the customer satisfaction survey to ensure the reliability in the qualitative analysis of the supply chain. Sitko-Lutek et al. (2010) examined the customer complaint handling process with respect to the information quality, thereby suggesting possible areas of improvements in the process. Their research method involved reviewing documents, complaint handling procedures and interviews through a social network analysis (SNA) model. The software used for SNA was UCI net and the results suggested that process engineering leadership played a vital and responsive role in disseminating quality assurance information in identifying potential areas of process improvements, thereby enhance and improve the company's profit and customers satisfaction. # III. # Methodology A control chart is a graphical representation that shows whether a sample data falls within a normal range of variation. It used to know if a process is in statistical quality control or not. It is also a graphical representation of mathematical model used to monitor a process in order to detect changes in parameter of that process. It displays the quality characteristics that has been measured or computed from a sample against the sample number or time. They are simple to construct and to interpret as they employ a center line (denoted as CNL) and two major control limits; an upper control limit (denoted as UCL) and a lower control limit (denoted as LCL). The center line represents the average performance of the process when it is in a state of statistical control-that is, when only common cause variation exists. The upper and lower control limits are horizontal lines situated above and below the center line. These control limits are established so that when the process is in control, almost all plots will be between the upper and lower limits. In practice, -If all observed plot points are between the LCL and UCL and if no unusual pattern of points exists, we have no evidence that assignable causes exist and we assume that the process is in statistical control. In this case, only common causes of the process variation exist, and no action to remove assignable causes is taken on the process. If we were to take such action, we would be unnecessarily tempering with the process. -If we observe one or more plot points outside the control limits, then we have evidence that the process is out of control due to one or more assignable causes. Here we must take action on the process to remove those assignable causes. (1956,1959,1985), Crosier (1988), Hawkins (1991Hawkins ( , 1993b)), Lowry et al. (1992), Lowry and Montgomery (1995), Pignatiello and Runger (1990), Tracy, Young, and Mason (1992), Montgomery andWadsworth (1972), andAlt (1985). This subject is particularly important today, as automatic inspection procedures make it relatively easy to measure many parameters on each unit of product manufactured. control of two or more related quality characteristics is necessary. The process is considered to be in control only if the sample means x 1 and x 2 fall within their respective control limits. Monitoring these two quality characteristics independently can be very misleading. So it is best we use the HotellingT 2 control chart. # b) The Multivariate Normal Distribution In univariate statistical quality control, we generally use the Normal distribution to describe the behaviour of a continuous quality characteristic. The Univariate Normal probability density function is f(x)=1/(?2??^2 ) e^(-1/2 ?((x-?)/?)?^2 ) -? 100, most practitioners use an approximate control limit, either ?????? = ??(?? ? 1) ?? ? ?? ?? ?,??,????? (21) ?????? = ?? ?,?? 2(22) However, for m> 100, equation ( 21) is a reasonable approximation. For phase i, the limits are based on a beta distribution, ?????? = (?? ? 1) 2 ?? ?? ?, ?? 2 , ?? ????1 2 ?????? ?????? =0 (23) Where ?? ?,??/2,(?? ????1)/2 is the upper ? percentage point of a beta distribution with parameters p/2 and (m-p-1)/2. Approximations to the phase I limit based on the F and chi-square distributions are likely to be inaccurate. Basically, the focus will be on the Sub grouped data because it suits the type of data that was collected. # Control Chart for Monitoring Variability Monitoring multivariate process are in two levels, which are to monitor the process mean vector m and to monitor process variability. Process variability is summarized by the p x p covariance matrix ?. The main diagonal elements of this matrix are the variances of the individual process variables, and the off-diagonal elements are the covariances. We can use the approach based on the sample generalized variance, |S|. This statistic, which is the determinant of the sample covariance matrix, is a widely used measure of multivariate dispersion. Another method would be to use the mean and variance of |S|, that is, E(|S|) and V(|S|), and the property that most of the probability distribution of |S| is contained in the interval # ??(|S|) ± 3?(??(|S|)). It can be shown that ??(|??|) = ?? 1 |?| ?????? ??(|??|) = ?? 2 |?| 2 (24) where ?? 1 = 1 (?? ? 1) ?? ?(?? ? ??) ?? ??=1 (25) and ?? 2 = 1 (?? ? 1) 2?? ?(?? ? ??) ?? ??=1 ??(?? ? ?? + 2) ?? ?? = 1 ? ?(?? ? ??) ?? ?? =1 ? (26) Therefore, the parameters of the control charts for |S| would be Investigating the Quality Performance of Production of Some Selected Drinks using Hotelling T-square and Control Chart ) g) ?????? = |?| ??? 1 + 3?? 2 1 2 ? ???? = ?? 1 |?| ?????? = |?|(?? 1 + 3?? 2 1/2 )(27)13 The lower control limit in equation ( 27) is replaced with zero if the calculated value is less than zero. In practice, ? usually will be estimated by a sample covariance matrix S, based on the analysis of preliminary samples. If this is the case, we should replace |?| in equation ( 27 In this study, two measurement quality characteristics are being analyzed using Multivariate statistical quality control. Fill height: It measures the level of liquid in a bottle of drink. The products under study are STAR, MALTINA and GOLDBERG from Nigerian Breweries plc. The standard is always at 60cl. Co 2 level: It measures the level of co 2 in each bottle. The target for corking a bottle of STAR is between (0.52-0.54%wt/wt), that of MALTINA is (0.59-0.61%wt/wt) and GOLDBERG is (0.62-0.64%wt/wt ) where %wt/wt means weight per weight. Data Presentation: The data used for this analysis is shown in the appendix 'A to appendix F. # IV. # Data Analysis and Results In this chapter, the Hotelling T 2 control chart is used for the analysis of fill height and level of co 2 measurements of Star, Maltina and Goldberg using R. # Analysis on the Fill Height Measurement and co 2 Level of Star The fill height of STAR refers to the height of the liquid content in a bottle of a STAR. And the co 2 level refers to the level of co 2 in each bottle of STAR. There can be cases of low fill, high fill and normal fill. The normal or standard fill height of STAR of the company is 60cl. And the standard co 2 level of STAR is between (0.52-0.54%wt/wt). The tables below display analysis carried out using R on various readings on fill height and co 2 level that was observed at different times. Table1.0 shows the Variances and Covariances of the fill height and level of CO 2 of STAR and also the Hotelling T 2 and Variability of each of the 20 samples. The Grand mean, Variance-Covariance Matrix (s) for the control limit used in the Variability plot, and the control Limits for the Hotelling T 2 and Variability plot are represented in the table(s) below. # Grand Mean # h) Interpretation of Star Chart From the Variability plot above, most of the sample variances are on or close to the lower control limit (LCL) while they are very far from the upper control limit, which means that the variability (the variances of the observation from the mean) is in control. Thus, the Hotelling T 2 can be plotted to see if the process is actually in control. From the Hotelling T 2 plotted above also, it can be seen that all the plot point fall within the UCL and LCL, which means that it can be concluded that the fill height and level of Co 2 of STAR is under control. The R code was used for the analysis of STAR. # Analysis on the Fill Height Measurement and co 2 Level of Maltina. The fill height of MALTINA refers to the height of the liquid content in a bottle of a MALTINA. And the co 2 level refers to the level of co 2 in each bottle of MALTINA. There can be cases of low fill, high fill and normal fill. The normal or standard fill height of MALTINA of the company is 60cl. And the standard co 2 level of MALTINA is between (0.59-0.61%wt/wt). The tables below display analysis carried out using R on various readings on fill height and co 2 level that was observed at different times. Table2 shows the Variances and Covariances of the fill height and level of Co 2 of MALTINA and also the Hotelling T 2 and Variability of each of the 20 samples. The Grand mean, Variance-Covariance Matrix (s) for the control limit used in the Variability plot, and the control Limits for the Hotelling T 2 and Variability plot are represented in the table(s) below. From the Hotelling T 2 plotted above, it can be seen that all the plot point fall within the UCL and LCL, which means that the fill height and level of co 2 of MALTINA is under control. # Grand Mean # Variance-Covariance Matrix (S) for the control limit used in the variability plot # Analysis on the Fill Height Measurement and Co 2 Level oif Goldberg The fill height of GOLDBERG refers to the height of the liquid content in a bottle of a GOLDBERG. And the co 2 level refers to the level of co 2 in each bottle of GOLDBERG. There can be cases of low fill, high fill and normal fill. The normal or standard fill height of GOLDBERG of the company is 60cl. And the standard level of co 2 GOLDBERG is between (0.62-0.64%wt/wt). The tables below display analysis carried out using R on various readings on fill height and co 2 level that was observed at different times. From the Hotelling T2 plotted above also, it shows that all the plot point fall within the UCL and LCL, which means that the fill height and level of Co of GOLDBERG is under control. Based on the results obtained from the analysis so far for all the drinks considered, none of the characteristics examined and analyzed fall within the control which invariably means we do not have sufficient evidence to reject the null hypothesis hence we Accept the null hypothesis for both the fill height and the level of Co2. V. # Conclusion The results obtained from the method used show that the components for the production of the beer under consideration(fill height and level of Co2) shows that the variability of the three products are in control, and this information helped in proceeding to check if the two quality characteristics are in control, also, using the Hotelling T2 control chart of Sub grouped data, the values were all within the lower and upper control limit for the three products, which helps to affirm the fact that the quality characteristics of STAR, MALTINA AND GOLDBERG are in control. This shows that the Quality Control Unit of the Company should not relent in carrying out their test on the products, all these will help the company to maintain the required standard and survive competition with other likely products from other company. ![Investigating the Quality Performance of Production of Some Selected Drinks using Hotelling T-square and Control Chart](image-2.png "") ![??,???? ??? ???+1 ?????? Global Journal of Management and Business Research Volume XX Issue X Version I Investigating the Quality Performance of Production of Some Selected Drinks using Hotelling T-square and Control Chart ?????? = 0 (18)](image-3.png "") ![) by |S|/b 1](image-4.png "") 1![Figure 1: HotellingT-Square plot of Star](image-5.png "ControlFigure 1 :") 1MeansVariance and CovariancesControl Chart StatisticsSampleFill Height (?? ? 1k )Level of CO 2S 2 1kS 2 2kS 12kHotelling T 2 k|S k |Number k(?? ? 2k )160.00.5181.5 0.000370.02003.8339370 0.00015500259.60.5301.3 0.000350.01000.1962245 0.00035500359.80.5200.7 0.000050.00002.6772350 0.00003500460.00.5301.0 0.000100.00250.1181133 0.00009375559.80.5462.2 0.000130.00157.7073474 0.00028375660.00.5261.0 0.00013-0.00250.4357597 0.00012375759.40.5161.3 0.00008-0.00056.4532623 0.00010375860.00.5221.5 0.000070.00751.6743676 0.00004875960.00.5482.5 0.000120.000010.0856025 0.000300001059.40.5280.3 0.00017-0.00650.8187689 0.000008751160.00.5141.0 0.00013-0.00256.9144678 0.000123751259.60.5340.8 0.00003-0.00300.6801556 0.000015001360.00.5481.5 0.000170.000010.0856025 0.000255001459.80.5362.7 0.00053-0.02851.1678815 0.000618751560.00.5161.0 0.00023-0.01255.2590822 0.000073751659.60.5500.3 0.00005-0.002511.8254949 0.000008751760.00.5242.5 0.000030.00250.9399435 0.000068751859.80.5341.7 0.00008-0.00150.5507094 0.000133751959.80.5240.7 0.000130.00350.9184544 0.000078752060.00.5281.5 0.00057-0.02250.1618163 0.00034875Averages59.830.52961.5 0.00057-0.0225 2MeansVariances and CovariancesControl Chart StatisticsSample Number k Fill Height (?? ? 1k ) Level of CO 2 (?? ? 2k ) S 2 1kS 2 2kS 12kHotelling T 2 k|S k |160.00.5961.0 0.00003 0.00251.76067480.00002375259.60.5901.3 0.00010 0.01008.02084860.00003000359.80.5981.2 0.00007 0.00700.32965290.00003500460.00.5980.5 0.00017 0.00750.61618940.00002875559.60.5981.3 0.00002 -0.00100.42355470.00002500660.00.5940.5 0.00003 0.00253.58753380.00000875759.40.6021.3 0.00002 -0.00351.46221190.00001375 3MeansVariances and CovariancesControl Chart StatisticsSample Number kFill Height (?? ? 1k )Level of CO 2 (?? ? 2k ) S 2 1kS 2 2kS 12kHotelling k T 2|S k |160.00.6421.0 0.00007 0.0075000 2.28910604 0.00001375259.60.6361.3 0.00008 0.0030000 0.28786007 0.00009500359.80.6300.7 0.00025 -0.0025000 1.31283293 0.00016875460.00.6340.5 0.00008 0.0025000 0.25264057 0.00003375559.60.6260.3 0.00013 -0.0045000 4.15963167 0.00001875660.00.6301.0 0.00010 -0.0075000 1.58110868 0.00004375759.40.6380.8 0.00007 0.0060000 1.37548756 0.00002000860.00.6320.5 0.00037 0.0010000 0.72131622 0.00008500960.00.6463.5 0.00008 0.0000000 5.65403962 0.000280001059.40.6341.3 0.00003 -0.0045000 0.88317403 0.000018751160.00.6441.5 0.00008 -0.0025000 3.77601443 0.000113751259.60.6320.3 0.00007 0.0010000 0.66321829 0.000020001360.00.6320.5 0.00002 0.0000000 0.72131622 0.000010001459.80.6462.7 0.00008 0.0015000 5.72751642 0.000213751560.00.6461.0 0.00008 0.0050000 5.65403962 0.000055001659.60.6281.3 0.00007 0.0040000 2.60304373 0.000075001760.00.6340.5 0.00008 0.0000000 0.25264057 0.000040001859.80.6300.7 0.00010 -1.1 x 10 -22 1.31283293 0.000070001959.80.6340.7 0.00018 0.0010000 0.06980296 0.000125002060.00.6300.5 0.00005 0.00250000 1.58110868 0.00001875 3Investigating the Quality Performance of Production of Some Selected Drinks using Hotelling T-squareand Control ChartYear 202017Volume XX Issue X Version I( ) BGrand Mean Variance-Covariance Matrix (S) for the control limit used in the variability plot Fill Height Level of CO 2 Fill Height 1.0300000 0.0011250 Fill Height 59.820 Level of CO 2 0.6352Global Journal of Management and Business ResearchLevel of CO 20.00112500.0001035Control Limits for the Hotelling T 2 and Variability plotLCLCLUCLHotelling T 20-14.52384Variability00.0001053 0.000492© 2020 Global Journals B39:00am5960606159410:00am6061615959511:00am6059605862612 noon615960596171:00pm605959586182:00pm626059605993:00pm5862605961104:00pm6059596059115:00pm6159616059126:00pm6059596159137:00pm5860606161148:00pm5760616160159:00pm59615960611610:00pm60596060591711:00pm58606259611812:00am5960596259191:00am5960596061202:00am6060595962Sample NoTimeABCDE17:00am0.520.510.500.550.5128:00am0.550.500.540.530.5339:00am0.530.520.520.520.51410:00am0.540.540.520.520.53511:00am0.530.560.550.540.55612noon0.510.520.540.530.5371:00pm0.520.510.530.510.5182:00pm0.530.520.520.530.5193:00pm0.540.540.540.560.56104:00pm0.520.540.530.510.54115:00pm0.510.530.510.520.50126:00pm0.530.540.540.530.53137:00pm0.550.560.530.540.56148:00pm0.560.550.530.500.54159:00pm0.540.510.520.510.501610:00pm0.550.560.550.540.551711:00pm0.520.520.530.530.521812:00am0.540.540.520.530.54191:00am0.520.510.520.540.53202:00am0.510.540.530.560.50Note: A, B, C, D, and E are the numbers of observations for each samples respectively. DSample No.TimeABCDE17:00am0.600.600.590.600.5928:00am0.580.580.590.600.6039:00am0.590.600.600.610.59410:00am0.590.610.610.600.58511:00am0.600.600.600.600.59612noon0.590.590.590.600.6071:00pm0.600.600.610.600.6082:00pm0.600.610.600.590.5993:00p m0.590.610.600.590.60104:00pm0.600.600.600.600.61115:00pm0.620.610.610.600.60126:00pm0.610.610.610.600.61137:00pm0.610.620.600.610.60148:00pm0.620.600.600.610.62159:00pm0.600.610.600.610.601610:00pm0.610.590.600.610.601711:00pm0.600.600.600.600.611812:00am0.600.610.590.590.59191:00am0.590.600.580.600.61 FYear 202022Volume XX Issue X Version I( ) BGlobal Journal of Management and Business ResearchSample No 1 2 3 4 5 6 7 8 9 10 11 12 13 14Time 7:00am 8:00am 9:00am 10:00am 11:00am 12noon 1:00pm 2:00pm 3:00pm 4:00pm 5:00pm 6:00pm 7:00pm 8:00pmA 0.65 0.63 0.61 0.64 0.63 0.62 0.63 0.64 0.65 0.64 0.64 0.62 0.63 0.64B 0.64 0.63 0.63 0.64 0.62 0.64 0.65 0.66 0.66 0.64 0.63 0.64 0.63 0.66C 0.64 0.64 0.65 0.64 0.61 0.64 0.64 0.62 0.64 0.63 0.65 0.63 0.63 0.64D 0.65 0.65 0.62 0.63 0.64 0.63 0.63 0.63 0.64 0.63 0.65 0.64 0.63 0.64E 0.63 0.63 0.64 0.62 0.63 0.62 0.64 0.61 0.64 0.63 0.65 0.63 0.64 0.65159:00pm0.640.660.640.650.641610:00pm0.630.620.620.630.641711:00pm0.640.640.640.620.631812:00am0.620.640.640.620.63191:00am0.620.640.650.620.64202:00am0.620.630.630.630.64© 2020 Global Journals © 2020 Global Journals ## Appendix * The importance of intrinsic and extrinsic cues to expected and experienced quality: An empirical application for beef LBAcebrón DCDopico Food Quality and Preference 11 2000 * Understanding attitudes and predicting social behavior IAjzen MFishbein 1980 Prentice-Hall Englewood Cliffs, NJ 1st ed * The evolution of credence goods: A transaction approach to product specification and quality control ESAndersen 1994 MAPP working paper no 21 * Food-related life styles in Singapore: Preliminary testing of a Western European research instrument in Southeast Asia SAskegaard K&brunsø Journal of Euro marketing 7 1999 * Undersøgelseafmulighederne for eksportafø kologiskemael keproduktertil Nordtyskland ABaadsgaard KGGrunert SCGrunert HSkytte 1994 MAPP project paper * Integrating the theory of planned behaviour with means-end chain theory -A study of possible improvements in predictive ability TBech-Larsen KGGrunert 1998. May Stockholm 27th EMAC Conference * Production and Operations Management A life cycle approach 5th RBChaser N JAquilano 1989 * Irwin Series in Quantitative Analysis for Business Edn 164 * Integrated Materials Management; Therm. and E. 2ndEdn. Handbook Series RJCarter PMPrice 1995 McGraw-Hill Inc 648 Dervirsiotis KN. Operations Management; International Student Edn * CJames Benneyan & David Birn Methods in Infection Control and Hospital Epidemiology. Hospital Epidemiology 19 3 1998 * Statistical Quality Control of High Resolution Winds of Different Radios ode Types for Chmatology Analysis KarimHouchi Journal of Atmospheric and Oceanic Technology 10 32 2015. 2015 * Production Control for Supervisors; 1st Edn. Collins, London and Glasgow King-ScottPeter 1971 95 * Process Variability Reduction Through Statistica Process Control for Quality Improvement BPMahosh MSPrabhuswany Int. J. Quality Rese 4 3 2010 * Statistical quality control. A case study research MariaEmiliaCamergo WalterPaesnitzFilho AngelaIsabel Dos SantosBulius Jose Maurcio CarreMaciel 4th IEEE International Conference Mayer Raymond R1975 McGraw-Hill book Company Inc 2008. 2008. 2008 347 Mittag HJ and Rinne H1993 * An essential ingredient for improving service and manufacturing quality BMason AJJiju 2007 Statistical process control * Analysis of the Effect of TQM on competitive Advantage EShenawy TBaker International Journal of Quality and Reliability management 2007 * Statistics and quantitative methods for construction and business managers JIEze The Nigerian institute of Building 2005 * An economic model for X and R chart with time-varying parameters HOhta AKimura MARahim Quality and Reliability Engineering-International 18 2002 * Joint economic design of X and R charts under Weibull shock model MARahim AF BCosta International Journal of Production . Research 28 2000