# On the Charting Procedures:T2 Chart and DD-Diagram

Authors: Mekki Hajlaoui
Journal: International Journal of Quality, Statistics, and Reliability
 (2011)

---

## Abstract

Multivariate analysis is increasingly used to include all dimensions of quality concept, in
light of rapid development of customer requirements. With the recent advances in information
technology and in recording, large amounts of multivariate data are now needed to be analyzed.
Many charting procedures are based on Mahalanobis distance, but their applicability relies heavily
on the requirement of normality and their performance is related to the choice of a type I
error rate. An alternative charting scheme based on data depth is pursued and its performance
is assessed through a real example. This performance and that of aT2 chart for individual
observations are discussed. Using the centre-outward ranking, this new method named DD-diagram
is used to detect any multivariate quality datum that one of its components exceeds
its limiting variation interval. For a given error-free sample, the DD-diagram can be used to
signal out any point of another observed sample taken from a multivariate quality process. This
new scheme based on data depth uses a properly chosen limiting variation line or Lvalue in
order to evaluate the outlyingness of every point in the observed sample in all directions of the
considered P-variates of quality process.


---

## Body

## 1. Introduction

Control charts are standard tools that are used to monitor quality process to identify instability within the manufacturing process. In practice, the quality of a product is determined by the interaction of multiple characteristics that are correlated, it is a multivariate phenomenon by nature. So adequate techniques need to be used to monitor the multivariate quality process.Multivariate Shewhart control chart was first introduced in 1947. It was based on theT2 test statistic and known as Hotelling’s T2 chart. Then, a number of multivariate control charts were designed to suit different situations such as multivariate CUSUM and multivariate EWMA charts. These classical monitoring charts have been developed under a number of assumptions quoted by [1].The performance of the multivariate control charts relies heavily on the hypothesis that the underlying distribution of the quality process is multivariate normal. It is well known that in practice this hypothesis rarely holds. Alternative procedures are needed to overcome this limit. Based on a data depth notion, [2] refined a visual procedure named DD-diagram which uses data depth plot to monitor any multivariate quality data and does not require any assumptions about the underlying distribution of the process. This graphical method provides a visualization of a change in position and/or in scale between an empirical sample with respect to a reference one. Furthermore, this diagram permits for an easier interpretation and a rapid adjusting action in the multivariate quality process through the use of newly suggested control limits detecting any out-of-control signal.In this paper, an application of bothT2 control chart and DD-diagram is investigated using individual observations taken off a real case of quality process from the Tunisian industry. The data of the samples are collected during two different times of the production process. The reference sample measures are drawn from a production process during which the process is considered in control. However, the empirical sample measures are drawn later in the frame work of a quality control routine. The T2 chart and the DD-diagram are given in Sections 2 and 3, respectively. In Section 4  these monitoring techniques are applied. The results and concluding remarks are provided in Section 5.


## 2. TheT2 Control Chart

Letω be the quality of an item X at time period i, so if this level of quality is characterized by p quality characteristics X1,X2,…,Xp then Xi(ω)=[Xi1(ω),Xi2(ω),…,Xip(ω)]′ or simply(1)Xi=(Xi1,Xi2,…,Xip)′
is a vector-valued output at time i=1,2,…. The components of the sample X1,X2,… are assumed to be independent and identically distributed multivariate normal random variables with mean vector μ=[μ1,μ2,…,μp]′ and covariance matrix Σ.It is clear that theT2 control chart is the suitable monitoring technique for the above vector-valued output series. So, the above assumptions make it possible to say that(2)Qi=(Xi-μ)′Σ-1(Xi-μ)
has a χp2 distribution. The statistic Qi is the Mahalanobis distance of the vector Xi from the mean μ.When the mean vectorμ and the covariance matrix Σ are known, the control chart for the Xi’s series is constructed using the statistic in (2). As stated by [3], its lower control limit (LCL) and upper control limit (UCL) are specified by(3)LCL=χ2(1-α2;p),UCL=χ2(α2;p),
where χ2(α;p) is the 1-α percentile of the chi-distribution with p degrees of freedom. In reality, however, the parameters μ and Σ are rarely known so they have to be estimated from a base period of n observations when the process is in control. If the base period sample referred to as a reference sample of size n is denoted by X1,X2,…,Xn, then the mean vector μ is estimated by μ̂=X¯p×1=(x¯1,x¯2,…,x¯p)′ such that(4)x¯j=1n∑j=1nxj,j=1,2,…,p,
and the covariance matrix Σ is estimated by(5)Σ̂=Sp×p=1n-1∑i=1n(Xi-X¯)(Xi-X¯)′.
Replacing μ and Σ by their respective estimators in (2), the empirical version of the Qi statistic is obtained by(6)Qi=(Xi-X¯)′S-1(Xi-X¯).
When the ith observation Xi is independent of both X¯ and S, then according to [3] the statistic of (6) has an exact F distribution and the control limits of the T2 control chart are given by(7)LCL=p(n-1)(n+1)n(n-p)F(1-α/2;p,n-p),UCL=p(n-1)(n+1)n(n-p)F(α/2;p,n-p),
where F(α;p,n-p) is the 1-α percentile of the F distribution with p and n-p degrees of freedom. The process is declared out of control whenever an observation has its Qi value not lying in the band defined by (7).


## 3. The DD-Diagram

Many charting procedures are based on Mahalanobis distance, but their applicability relies heavily on the requirement of normality and their performance is related to the choice of a type I error rate. An alternative charting scheme based on data depth is pursued. For a given error free sample, this new method named DD-diagram can be used to signal out any point of another observed sample taken from a multivariate quality process. This new scheme based on data depth uses a properly chosen limiting variation line orLvalue in order to evaluate the outlyingness of every point in the observed sample in all directions of the considered P-variates of quality process.LetF be a probability distribution in ℜp,  p≥1. Throughout the following, unless stated otherwise, we assume that F is absolutely continuous and also that the reference sample (X1,X2,…,Xn) is derived from F. Therefore, if the quality of the  ith observed unit is denoted by ω, then Xi(ω)=(Xi1,Xi2,…,Xip)′.According to [4], a data depth is a way of measuring how deep or central a given point x(ω)∈ℜp is with respect to F or to a given data cloud (X1,X2,…,Xn). Then using (2), the Mahalanobis depth at x(ω) with respect to F is defined to be(8)D(F,x(ω))=11+(x(ω)-μ)′Σ-1(x(ω)-μ).
The sample version of (8) is obtained by replacing the mean vector μ and the covariance matrix Σ of F as in [4], with their respective sample estimates in (4) and (5), then(9)Dn(F,x(ω))=11+(x(ω)-X¯)′S-1(x(ω)-X¯).Henceforth,D or Dn will be used to indicate the datum depth notion and a larger value of D•(F,x(ω)) always implies a deeper (or more central) x(ω) with respect to F.Given a notion of data depth, one can compute the depths of all quality measures(X1,X2,…,Xn) and order them according to decreasing depth values. This gives a ranking of the sample point associated with the ith highest depth value. Consequently, X[1],X[2],…,X[n] is the order statistics, with X[1] being the deepest or the most central point or simply the centre, and X[n] the most outlying point. The implication is that a larger rank is associated with a more outlying position with respect to the data cloud. These order statistics induced by a data depth are different from the usual order statistics on the real line, since the latter are ordered from the smallest sample point to the largest, while the former start from the middle sample point and move outwards in all directions, see [4].Given definitions (8) or (9), the sample becomes X[1],X[2],…; there is a natural choice of location parameter for the observed distribution. Specifically, the centre ν(ω) is the most central point so(10)ν(ω)=X[1].When the depth-equivalence class contains more than one point measure,X[11],X[12],…,X[1K]  according to [4] the median or the centre is the average of the deepest points, so in this case(11)ν(ω)=1K∑k=1KX[1k].On this basis and using data depth, (10) and (11) fix out a centre or a multivariate median. Moreover, [4] stated that if Mahalanobis depth is used, the central point defined in (10) and (11) turns out to be the mean of the observed data. This suggests concepts of location which are intermediate between the mean and the median.A data depth plot is a graphical comparison between two multivariate distributions based on data depth. So in addition to the reference sample(X1,X2,…,Xn), let (Y1,Y2,…,Yn) be another sample referred to as an empirical sample of distribution G characterizing the observed quality process. The reference sample X(ω) is taken when the production process is in control and the empirical sample Y(ω) is drawn when inspecting the process during another period of time later. The DD-diagram is obtained when plotting data depth values of a sample versus the other. Precisely, the DD-diagram is defined by(12)DD(F,G)={(D(F,x(ω)),D(F,y(ω))),foranyx(ω),y(ω)∈Rp(Dn(F,x(ω)),Dn(F,y(ω)))}.SinceDD(·,·) is a subset of ℜ2, the resulting graph is one-dimensional curve in the plane. If the two distributions are identical, the DD-diagram in (12) turns out to be the diagonal line from the point (0,0) to (1,1). Different patterns of deviations from the diagonal line in the DD-diagram are indications of differences in specific characteristics of F and G.In general, the distributions are rarely known so instead we use an empirical version of the DD-diagram. IfF and G are unknown distributions for the samples X(ω) and Y(ω), then the DD-diagram is obtained when plotting(13)DD(Fn,Gn)={(Dn(F,x(ω)),Dn(F,y(ω))),foranyx(ω),y(ω)∈Rp(Dn(F,x(ω)),Dn(F,y(ω)))}.
if (9) is used to compute the data depth.Ifp≥2 and if F and G are both absolutely continuous, then DD-diagram corresponds to a region with nonzero area. The area of this region can serve as a measure of the discrepancy between F and G, see [4]. If the two distributions are identical, the data cloud of the DD-plot should be concentrated along the diagonal line. Other patterns are indications of differences in specific characteristics of F and G, that is, in position, in scale, in skewness, and so forth.In most cases, the departure from the diagonal line usually takes the form of pulling down from the point(0.5,0.5) to the origin, leaving the upper right corner empty and spreading out the points as a scatter plot diagram pointing at (0,0). In order to bring out scale differences, the centre of the samples should be equalized first by subtracting the obtained centre of (10) or (11) from the data points. Suppose that G is more spread out than the reference sample F, then the points in DD-diagram tend to arch toward the F sample around the origin.In analogous manner to the multivariate Shewhart chart, [2] has suggested the control limits Lmax=D•(F,x(ω))[2-D•(F,x(ω))] and  Lmin=1-1-D•2(F,x(ω)), in order to detect visually the shifts in location and/or in dispersion. The region under control is located between the limits Lmax and Lmin. This marked region turns out to be so large. In order to make the DD-diagram more sensitive, it is best to include the centre defined in (10) or in (11) in detecting any out of control shift so we suggest to use instead the following control limit:(14)Lvalue=1(p-1)×[D•(F,ν(ω))+Log(n+p-1)-1],
as a least critical value for a data depth according to which the corresponding point will be considered to have components not satisfying the variation intervals.This control limit,Lvalue, is conversely proportional to the sum of the depth of the reference sample centre, D•(F,ν(ω)), and the Log function of modified degrees of freedom of the empirical sample, Log(n+p-1), minus the depth of the closest point to the centre, 1, all multiplied by (p-1) with p being the number of the considered variables affecting quality. So, Lvalue is the defined function of the p-variates that characterize the quality of a product, the sample length, n, and the data depth of the deepest point (or the centre) of the reference observed distribution F. It appears to be sensitive to any change in the values of these three parameters. The role of this limiting value line is to decide graphically if the production process is in control or not. Therefore, if a point computed using (12) or (13) is lower than this Lvalue, the observed process is out of control and an investigation must be made in order to point out which characteristic is responsible among the considered p-variates.


## 4. Application ofT2 Chart and DD-Diagram

The basic idea is to collect individual observations data from a production process during which the process is considered in control. These observations are used to identify the parameters of the distributionF of the observed process. Then, another series of individual observations are drawn from the same process when its distribution has drifted to G. Both series of observations are used to construct and to argue the performance of the monitoring schemes: the T2 chart and the DD-diagram.


### 4.1. Quality Processing of the Firm

Kairouan Tobacco Manufacture (KTM) and National Tobacco and Matches Corporation represent the monopoly of the tobacco in Tunisia. This study focuses on the quality of cigarettes manufactured by KTM. Aiming to improve quality, to satisfy consumers, and to comply with the legislation standards, both companies have equipped their laboratories with equipments that enable a rigorous quality control during all phases of manufacturing process. KTM produces several types of cigarettes, but in this application only the quality data of “Cristal light” cigarettes will be seen. The quality control activities cover all steps of the production process. It begins with the supervision of manufacturing sheets of tobacco, until conditioning the cigarettes in packages. In the tobacco industry, quality is defined by a plenty of requirements imposed by state. These requirements are considered when setting up the cigarettes specifications. Hence, tolerance limits are defined for each measure of quality characteristic of the cigarettes. The firm KTM uses five (5) characteristics that ensure the quality of cigarettes. They are weight, module, humidity rate of tobacco, pulling resistance, and folding capacity. These characteristics are hermetically measured to cope with the State norms.(i)
The weight of a cigarette is made up of the tobacco, the filter, and the cigarette paper weights. It varies between 0.965 and 1 gram.(ii)
The module of a cigarette corresponds to its diameter, varying from 6.75 to 8.0 mm.(iii)
The humidity rate of tobacco is the proportion of water contained in a cigarette. It is considered acceptable if it varies between 11.5% and 13.5%.(iv)
Pulling resistance of a cigarette is defined by the difference in pressure between the two extremities of a cigarette when a quantity of air is passed through it. The pulling resistance is considered acceptable when it varies from 100 to 115 CE (colonne d’eau).(v)
The folding capacity or density corresponds to the volume occupied by the mass of the tobacco inside a cigarette. It is tolerable to belong to450±20 cm3.


### 4.2. Collecting and Processing the Data

As introduced earlier, the collection of data is made up at two different periods of time. In fact, 60 measures are drawn as the reference sample on February 1, 2004, when the process is in control, whereas the empirical sample, of size 60 measures, is drawn on September 30, 2009. These collected samples are given in Table2.Like any multivariate quality process, the observation involves the simultaneous measurement of five variables as indicated: (1) the weight, (2) the module, (3) the humidity rate, (4) the pulling resistance, and (5) the density.Processing the observed data begins with the start-up stage that consists of estimating the parameters ofF, constructing the control limits of T2 chart, and determining the centre of F as reference sample.According to (4), the vector of sample mean is(15)X¯5×1=160∑j=160Xj=(0.9859,7.6501,0.1206,107.1833,451.5267)′
and to (5), the sample covariance matrix is(16)S5×5=159∑j=160(Xj-X¯)(Xj-X¯)′=[0.0001.1⋅10-30.000-8.2⋅10-36⋅10-3*8.89⋅10-2-1.4⋅10-3-0.71210.4635**0.0001.2⋅10-2-7.5⋅10-3***10.7963-3.8999****2.9820].To construct a multivariate control chart, control limits are determined according to (7) with α=0.05 as follows:(17)LCL=5(60-1)(60+1)60(60-5)0.1631=0.88926,UCL=5(60-1)(60+1)60(60-5)2.8073=15.3082.Figure1 shows the corresponding T2 chart. The observations were examined individually to determine a possible assignable cause and no observation is detected lying outside the in-control region specified by the above computed control limits.Figure 1
Multivariate control chart for the reference sample.To determine the centre, data depths of all observations of the reference sample are calculated using (9). As recorded in Table 2, the highest value of the statistic data depth is D60(F,x(ω))=0.38293 for a sample of size n=60. It corresponds to observation X45=(0.982,7.349,0.126,111,450)′ in the sample F so the most central point is X45=ν(ω) defining a depth-equivalence class of order one containing a single cigarette with rank 45.In order to detect graphically any point that is not satisfying the limiting variation intervals, theLvalue for the minimum data depths is calculated according to(14)(18)Lvalue=1(5-1)×(0.38293+Log(60+5-1)-1),=0.0705
as the least acceptable data depth value and below which the corresponding point is considered out of control, that is, at least one of the p-characteristics exceeds its limiting variation interval.When constructing DD-diagram in Figure2, there is no restriction made to the data distribution generated by the production process. So to obtain both subplots of Figure 2, (13) is used and instead of data depth of G distribution, the data depth of F distribution is employed in the x-axis and y-axis, respectively. This work is done twice before and after centring measures with respect to the identified centre of F distribution.DD-diagram for the sampleF versus itself and the limiting value.
(a)(b)If the observed measures are centred with respect to the centre ofF, this most central point will be characterized by the highest depth (or by the smallest Mahalanobis distance). So, it is clear that in the reference sample, the point of order 45 is characterized by the maximum data depth in either case before and after centering with respect to the computed vector-valued centre ν(ω). This fact explains why Lvalue is affine invariant and is used simultaneously in DD-diagram and DD-diagram of centred measures, see [4]. This deepest point in the reference sample is marked by a circle in both subplots of Figure 2.The second stage consists of using both control schemes to evaluate the stability of the observed production process when an empirical sample is drawn. At this phase, the parameters of the reference sampleF obtained in the start-up stage are used to monitor any taken empirical sample in the future. Specifically, after drawing the empirical sample G given in Table 2 the vector-valued mean X¯ and the covariance matrix S are used to assess the charting statistic Qi and to construct the T2 chart.Figure3 shows the corresponding multivariate control chart. Observations 19, 33, 50, 52, and 54 of the empirical sample lie outside the in-control region. These points are examined thoroughly to determine which characteristics are causing this drift in quality. It indicates that cigarettes 19 and 54 are considered out of control because the “module" (X2) exceeds its specified measure. The other observations 33, 50, and 52 are considered out of control because the humidity rate (X3) is lower than its minimum value.Figure 3
Multivariate control chart for the empirical sample.Henceforth, the vector-valued centreν(ω) and the limiting variation value Lvalue are used to evaluate the data depths of all measures making up the G sample and also to identify any change in location and/or in scale of the process.Figure4 shows the corresponding DD-diagram for both cases before and after centering measures. Both subplots (left and right) of Figure 4 sketch out any change in location and/or a scale increase when moving from the distribution F to the distribution G in the multivariate quality process.DD-diagram for the reference sampleF versus empirical sample G.Before centering measures, the DD-diagram gives a larger set of out-of-control observations than theT2 chart (with type I error rate α=0.05), not only the points 19, 33, 50, 52, and 54 but also 2, 3, 27, 29, 39, and 43 of the empirical sample G. These points have data depths lower than the computed limiting variation value because of a change in location and/or in scale, respectively. They are indicated with red stars in the left subplot of Figure 4. In order to determine the characteristics responsible for this drift, refer to Table 1.Table 1
Out-of-control observations detected by DD-diagram.
Obs.X1X2X3X4X5#[0.965, 1][6.75, 8][0.115, 0.135][100, 115][430, 470]20.9917.9560.113105454.6030.9827.9000.114101454.30190.9968.0000.117101456.13270.9877.8500.114102453.00290.9997.9500.114107452.10330.9917.9510.113109454.21390.9937.9500.115111454.00430.9967.9400.114109453.88500.9827.9590.111106454.12520.9877.9400.113107455.00540.9878.0100.115106450.10Table 2
SamplesF and G and the charting statistics for T2 and DD-diagram.
Obs.Reference sampleFEmpirical sampleGCharting statistics#X1X2X3X4X5X1X2X3X4X5QiD60(F,X)D60(F,Y)D60(F,Xc)D60(F,Yc)10.9857.3690.126107450.50.9897.9730.115105453.51.70790.19100.36930.06110.014420.9777.3590.126107449.30.9917.9560.113105454.613.29640.14740.06990.03410.009030.9827.9020.115104453.60.9827.9000.114101454.314.82990.12270.06320.38290.382940.9867.3530.125111450.20.9917.9560.115105453.21.86770.27440.34870.03860.009050.9787.3970.124110449.60.9917.9730.114107453.15.97690.22710.14330.05260.009060.9887.3850.126111450.10.9917.9730.114105453.64.12950.20800.19500.01910.009070.9927.9920.116103452.60.9927.9900.117103454.710.14620.19490.08970.00740.007480.9927.3600.124109449.90.9837.9800.115106454.16.84540.10010.12750.00740.212290.9897.9580.115107453.30.9907.9580.115108453.54.90730.23180.16930.01440.0113100.9917.3660.125109449.70.9927.9450.115105453.83.14940.14280.24100.00900.0074110.9807.4010.125108451.10.9917.9570.115104454.65.23960.23450.16030.17090.0090120.9927.9750.115107454.30.9928.0000.114107454.67.79740.15800.11370.00740.0074130.9857.3670.125110448.90.9947.9400.116108453.95.48390.20000.15420.06110.0052140.9897.9570.115105452.90.9917.9570.114105453.14.83440.35010.17140.01440.0090150.9817.3640.124109449.60.9918.0000.114105453.42.91410.24700.25550.33430.0090160.9827.3780.124108449.40.9917.9560.115105453.11.87810.24030.34740.38290.0090170.9907.9730.115106451.90.9917.9730.115106454.34.39630.16550.18530.01130.0090180.9827.3770.126109450.30.9907.9500.115107454.15.37300.30640.15690.38290.0113190.9967.9750.116101454.00.9968.0000.117101456.125.03700.11330.03840.00390.0039200.9827.3770.125106450.80.9897.9550.115107453.33.39950.20410.22730.38290.0144210.9817.3740.126110450.00.9917.9800.114103454.55.66200.33270.15010.33430.0090220.9907.9720.115102452.80.9917.9670.115106453.22.16880.24830.31560.01130.0090230.9877.9940.116105453.10.9897.9900.116103453.03.06210.25770.24620.02640.0144240.9797.4000.125112450.10.9927.9650.115103453.12.26590.23180.30620.08930.0074250.9787.3630.127115449.80.9897.9560.114105453.24.70440.09050.17530.05260.0144260.9857.9960.116106452.20.9858.0000.117106453.17.64410.14600.11570.06110.0611270.9877.3420.124107450.60.9877.8500.114102453.013.54890.13560.06870.02640.0264280.9897.9560.115102454.70.9907.9560.115107454.04.69220.13160.17570.01440.0113290.9837.3480.125107449.30.9997.9500.114107452.115.19340.21890.06180.21220.0027300.9907.9570.115105452.60.9917.9530.115107453.63.84670.30190.20630.01130.0090310.9827.9770.117103452.30.9807.9770.117105454.012.62230.08750.07340.38290.1709320.9887.4150.126111450.40.9917.9570.115105453.92.67830.17310.27190.01910.0090330.9817.3280.124108450.20.9917.9510.113109454.219.57490.17250.04860.33430.0090340.9937.9840.116110453.10.9947.9800.116110454.19.56290.11780.09470.00620.0052350.9867.3710.126110448.90.9917.9600.114105453.04.74490.17030.17410.03860.0090360.9877.9430.116103452.10.9867.9400.116105453.11.73260.21110.36600.02640.0386370.9817.4470.127109450.70.9917.9560.114106453.96.53030.07360.13280.33430.0090380.9907.9470.115104453.70.9917.9400.115104454.34.72040.29510.17480.01130.0090390.9837.3600.124108450.20.9937.9500.115111454.013.82000.25750.06750.21220.0062400.9817.3440.125115449.90.9897.9700.11497453.412.98970.10530.07150.33430.0144410.9907.9730.115103453.90.9907.9690.115106453.52.24210.29060.30840.01130.0113420.9857.3880.125108449.20.9917.9520.115105453.02.00240.22150.33310.06110.0090430.9997.9490.117109453.50.9967.9400.114109453.915.04540.09170.06230.00270.0039440.9757.3470.124110450.80.9937.9740.115105453.12.12980.10850.31950.01740.0062450.9827.3490.126111450.00.9907.9700.115105453.11.56880.38290.38930.38290.0113460.9917.9570.115105452.10.9917.9570.115105454.13.44050.20110.22520.00900.0090470.9907.9740.115106452.40.9907.9740.115106453.02.06850.24070.32590.01130.0113480.9967.9420.116102454.60.9937.9400.116103454.14.14160.10360.19450.00390.0062490.9827.3860.125111450.60.9897.9580.115105453.11.73230.34440.36600.38290.0144500.9827.9780.115101452.10.9827.9590.111106454.133.05670.09610.02940.38290.3829510.9907.9720.115106452.60.9917.9730.115106453.62.35670.27370.29790.01130.0090520.9867.9480.116107454.60.9877.9400.113107455.021.01870.11200.04540.03860.0264530.9927.9580.116106452.40.9927.9580.116106454.13.49980.27730.22220.00740.0074540.9877.9760.115104454.60.9878.0100.115106450.121.07300.15660.04530.02640.0264550.9807.3850.125107451.20.9927.9560.115105453.12.10560.18730.32200.17090.0074560.9827.3950.125108451.11.0007.9700.115107454.110.03720.26880.09060.38290.0024570.9977.9880.117106453.70.9977.9800.115106454.77.93130.11780.11200.00340.0034580.9777.4030.124107450.00.9917.9700.115106453.52.29700.21370.30330.03410.0090590.9807.3770.125113450.20.9917.9560.115106453.02.41040.18940.29320.17090.0090600.9937.9560.115102453.30.9907.9560.115106453.12.28810.22070.30410.00620.0113For centered measures of both samples with respect to the centre ofF, the DD-diagram in the right subplot of Figure 4 shows the observations under the effect of a scale change only. It can be deduced that observations 3 and 50 are out of control because of a change in location only. When measures are centered, these two observations disappear from the out of control region in the right subplot of Figure 4. But, observations 2, 19, 27, 29, 33, 39, 43, 52, and 54 are out-of-control under the effect of a change in dispersion and of a location shift, respectively. These points are indicated with red stars in the DD-diagram in the right subplot of Figure 4.


## 4.1. Quality Processing of the Firm

Kairouan Tobacco Manufacture (KTM) and National Tobacco and Matches Corporation represent the monopoly of the tobacco in Tunisia. This study focuses on the quality of cigarettes manufactured by KTM. Aiming to improve quality, to satisfy consumers, and to comply with the legislation standards, both companies have equipped their laboratories with equipments that enable a rigorous quality control during all phases of manufacturing process. KTM produces several types of cigarettes, but in this application only the quality data of “Cristal light” cigarettes will be seen. The quality control activities cover all steps of the production process. It begins with the supervision of manufacturing sheets of tobacco, until conditioning the cigarettes in packages. In the tobacco industry, quality is defined by a plenty of requirements imposed by state. These requirements are considered when setting up the cigarettes specifications. Hence, tolerance limits are defined for each measure of quality characteristic of the cigarettes. The firm KTM uses five (5) characteristics that ensure the quality of cigarettes. They are weight, module, humidity rate of tobacco, pulling resistance, and folding capacity. These characteristics are hermetically measured to cope with the State norms.(i)
The weight of a cigarette is made up of the tobacco, the filter, and the cigarette paper weights. It varies between 0.965 and 1 gram.(ii)
The module of a cigarette corresponds to its diameter, varying from 6.75 to 8.0 mm.(iii)
The humidity rate of tobacco is the proportion of water contained in a cigarette. It is considered acceptable if it varies between 11.5% and 13.5%.(iv)
Pulling resistance of a cigarette is defined by the difference in pressure between the two extremities of a cigarette when a quantity of air is passed through it. The pulling resistance is considered acceptable when it varies from 100 to 115 CE (colonne d’eau).(v)
The folding capacity or density corresponds to the volume occupied by the mass of the tobacco inside a cigarette. It is tolerable to belong to450±20 cm3.


## 4.2. Collecting and Processing the Data

As introduced earlier, the collection of data is made up at two different periods of time. In fact, 60 measures are drawn as the reference sample on February 1, 2004, when the process is in control, whereas the empirical sample, of size 60 measures, is drawn on September 30, 2009. These collected samples are given in Table2.Like any multivariate quality process, the observation involves the simultaneous measurement of five variables as indicated: (1) the weight, (2) the module, (3) the humidity rate, (4) the pulling resistance, and (5) the density.Processing the observed data begins with the start-up stage that consists of estimating the parameters ofF, constructing the control limits of T2 chart, and determining the centre of F as reference sample.According to (4), the vector of sample mean is(15)X¯5×1=160∑j=160Xj=(0.9859,7.6501,0.1206,107.1833,451.5267)′
and to (5), the sample covariance matrix is(16)S5×5=159∑j=160(Xj-X¯)(Xj-X¯)′=[0.0001.1⋅10-30.000-8.2⋅10-36⋅10-3*8.89⋅10-2-1.4⋅10-3-0.71210.4635**0.0001.2⋅10-2-7.5⋅10-3***10.7963-3.8999****2.9820].To construct a multivariate control chart, control limits are determined according to (7) with α=0.05 as follows:(17)LCL=5(60-1)(60+1)60(60-5)0.1631=0.88926,UCL=5(60-1)(60+1)60(60-5)2.8073=15.3082.Figure1 shows the corresponding T2 chart. The observations were examined individually to determine a possible assignable cause and no observation is detected lying outside the in-control region specified by the above computed control limits.Figure 1
Multivariate control chart for the reference sample.To determine the centre, data depths of all observations of the reference sample are calculated using (9). As recorded in Table 2, the highest value of the statistic data depth is D60(F,x(ω))=0.38293 for a sample of size n=60. It corresponds to observation X45=(0.982,7.349,0.126,111,450)′ in the sample F so the most central point is X45=ν(ω) defining a depth-equivalence class of order one containing a single cigarette with rank 45.In order to detect graphically any point that is not satisfying the limiting variation intervals, theLvalue for the minimum data depths is calculated according to(14)(18)Lvalue=1(5-1)×(0.38293+Log(60+5-1)-1),=0.0705
as the least acceptable data depth value and below which the corresponding point is considered out of control, that is, at least one of the p-characteristics exceeds its limiting variation interval.When constructing DD-diagram in Figure2, there is no restriction made to the data distribution generated by the production process. So to obtain both subplots of Figure 2, (13) is used and instead of data depth of G distribution, the data depth of F distribution is employed in the x-axis and y-axis, respectively. This work is done twice before and after centring measures with respect to the identified centre of F distribution.DD-diagram for the sampleF versus itself and the limiting value.
(a)(b)If the observed measures are centred with respect to the centre ofF, this most central point will be characterized by the highest depth (or by the smallest Mahalanobis distance). So, it is clear that in the reference sample, the point of order 45 is characterized by the maximum data depth in either case before and after centering with respect to the computed vector-valued centre ν(ω). This fact explains why Lvalue is affine invariant and is used simultaneously in DD-diagram and DD-diagram of centred measures, see [4]. This deepest point in the reference sample is marked by a circle in both subplots of Figure 2.The second stage consists of using both control schemes to evaluate the stability of the observed production process when an empirical sample is drawn. At this phase, the parameters of the reference sampleF obtained in the start-up stage are used to monitor any taken empirical sample in the future. Specifically, after drawing the empirical sample G given in Table 2 the vector-valued mean X¯ and the covariance matrix S are used to assess the charting statistic Qi and to construct the T2 chart.Figure3 shows the corresponding multivariate control chart. Observations 19, 33, 50, 52, and 54 of the empirical sample lie outside the in-control region. These points are examined thoroughly to determine which characteristics are causing this drift in quality. It indicates that cigarettes 19 and 54 are considered out of control because the “module" (X2) exceeds its specified measure. The other observations 33, 50, and 52 are considered out of control because the humidity rate (X3) is lower than its minimum value.Figure 3
Multivariate control chart for the empirical sample.Henceforth, the vector-valued centreν(ω) and the limiting variation value Lvalue are used to evaluate the data depths of all measures making up the G sample and also to identify any change in location and/or in scale of the process.Figure4 shows the corresponding DD-diagram for both cases before and after centering measures. Both subplots (left and right) of Figure 4 sketch out any change in location and/or a scale increase when moving from the distribution F to the distribution G in the multivariate quality process.DD-diagram for the reference sampleF versus empirical sample G.Before centering measures, the DD-diagram gives a larger set of out-of-control observations than theT2 chart (with type I error rate α=0.05), not only the points 19, 33, 50, 52, and 54 but also 2, 3, 27, 29, 39, and 43 of the empirical sample G. These points have data depths lower than the computed limiting variation value because of a change in location and/or in scale, respectively. They are indicated with red stars in the left subplot of Figure 4. In order to determine the characteristics responsible for this drift, refer to Table 1.Table 1
Out-of-control observations detected by DD-diagram.
Obs.X1X2X3X4X5#[0.965, 1][6.75, 8][0.115, 0.135][100, 115][430, 470]20.9917.9560.113105454.6030.9827.9000.114101454.30190.9968.0000.117101456.13270.9877.8500.114102453.00290.9997.9500.114107452.10330.9917.9510.113109454.21390.9937.9500.115111454.00430.9967.9400.114109453.88500.9827.9590.111106454.12520.9877.9400.113107455.00540.9878.0100.115106450.10Table 2
SamplesF and G and the charting statistics for T2 and DD-diagram.
Obs.Reference sampleFEmpirical sampleGCharting statistics#X1X2X3X4X5X1X2X3X4X5QiD60(F,X)D60(F,Y)D60(F,Xc)D60(F,Yc)10.9857.3690.126107450.50.9897.9730.115105453.51.70790.19100.36930.06110.014420.9777.3590.126107449.30.9917.9560.113105454.613.29640.14740.06990.03410.009030.9827.9020.115104453.60.9827.9000.114101454.314.82990.12270.06320.38290.382940.9867.3530.125111450.20.9917.9560.115105453.21.86770.27440.34870.03860.009050.9787.3970.124110449.60.9917.9730.114107453.15.97690.22710.14330.05260.009060.9887.3850.126111450.10.9917.9730.114105453.64.12950.20800.19500.01910.009070.9927.9920.116103452.60.9927.9900.117103454.710.14620.19490.08970.00740.007480.9927.3600.124109449.90.9837.9800.115106454.16.84540.10010.12750.00740.212290.9897.9580.115107453.30.9907.9580.115108453.54.90730.23180.16930.01440.0113100.9917.3660.125109449.70.9927.9450.115105453.83.14940.14280.24100.00900.0074110.9807.4010.125108451.10.9917.9570.115104454.65.23960.23450.16030.17090.0090120.9927.9750.115107454.30.9928.0000.114107454.67.79740.15800.11370.00740.0074130.9857.3670.125110448.90.9947.9400.116108453.95.48390.20000.15420.06110.0052140.9897.9570.115105452.90.9917.9570.114105453.14.83440.35010.17140.01440.0090150.9817.3640.124109449.60.9918.0000.114105453.42.91410.24700.25550.33430.0090160.9827.3780.124108449.40.9917.9560.115105453.11.87810.24030.34740.38290.0090170.9907.9730.115106451.90.9917.9730.115106454.34.39630.16550.18530.01130.0090180.9827.3770.126109450.30.9907.9500.115107454.15.37300.30640.15690.38290.0113190.9967.9750.116101454.00.9968.0000.117101456.125.03700.11330.03840.00390.0039200.9827.3770.125106450.80.9897.9550.115107453.33.39950.20410.22730.38290.0144210.9817.3740.126110450.00.9917.9800.114103454.55.66200.33270.15010.33430.0090220.9907.9720.115102452.80.9917.9670.115106453.22.16880.24830.31560.01130.0090230.9877.9940.116105453.10.9897.9900.116103453.03.06210.25770.24620.02640.0144240.9797.4000.125112450.10.9927.9650.115103453.12.26590.23180.30620.08930.0074250.9787.3630.127115449.80.9897.9560.114105453.24.70440.09050.17530.05260.0144260.9857.9960.116106452.20.9858.0000.117106453.17.64410.14600.11570.06110.0611270.9877.3420.124107450.60.9877.8500.114102453.013.54890.13560.06870.02640.0264280.9897.9560.115102454.70.9907.9560.115107454.04.69220.13160.17570.01440.0113290.9837.3480.125107449.30.9997.9500.114107452.115.19340.21890.06180.21220.0027300.9907.9570.115105452.60.9917.9530.115107453.63.84670.30190.20630.01130.0090310.9827.9770.117103452.30.9807.9770.117105454.012.62230.08750.07340.38290.1709320.9887.4150.126111450.40.9917.9570.115105453.92.67830.17310.27190.01910.0090330.9817.3280.124108450.20.9917.9510.113109454.219.57490.17250.04860.33430.0090340.9937.9840.116110453.10.9947.9800.116110454.19.56290.11780.09470.00620.0052350.9867.3710.126110448.90.9917.9600.114105453.04.74490.17030.17410.03860.0090360.9877.9430.116103452.10.9867.9400.116105453.11.73260.21110.36600.02640.0386370.9817.4470.127109450.70.9917.9560.114106453.96.53030.07360.13280.33430.0090380.9907.9470.115104453.70.9917.9400.115104454.34.72040.29510.17480.01130.0090390.9837.3600.124108450.20.9937.9500.115111454.013.82000.25750.06750.21220.0062400.9817.3440.125115449.90.9897.9700.11497453.412.98970.10530.07150.33430.0144410.9907.9730.115103453.90.9907.9690.115106453.52.24210.29060.30840.01130.0113420.9857.3880.125108449.20.9917.9520.115105453.02.00240.22150.33310.06110.0090430.9997.9490.117109453.50.9967.9400.114109453.915.04540.09170.06230.00270.0039440.9757.3470.124110450.80.9937.9740.115105453.12.12980.10850.31950.01740.0062450.9827.3490.126111450.00.9907.9700.115105453.11.56880.38290.38930.38290.0113460.9917.9570.115105452.10.9917.9570.115105454.13.44050.20110.22520.00900.0090470.9907.9740.115106452.40.9907.9740.115106453.02.06850.24070.32590.01130.0113480.9967.9420.116102454.60.9937.9400.116103454.14.14160.10360.19450.00390.0062490.9827.3860.125111450.60.9897.9580.115105453.11.73230.34440.36600.38290.0144500.9827.9780.115101452.10.9827.9590.111106454.133.05670.09610.02940.38290.3829510.9907.9720.115106452.60.9917.9730.115106453.62.35670.27370.29790.01130.0090520.9867.9480.116107454.60.9877.9400.113107455.021.01870.11200.04540.03860.0264530.9927.9580.116106452.40.9927.9580.116106454.13.49980.27730.22220.00740.0074540.9877.9760.115104454.60.9878.0100.115106450.121.07300.15660.04530.02640.0264550.9807.3850.125107451.20.9927.9560.115105453.12.10560.18730.32200.17090.0074560.9827.3950.125108451.11.0007.9700.115107454.110.03720.26880.09060.38290.0024570.9977.9880.117106453.70.9977.9800.115106454.77.93130.11780.11200.00340.0034580.9777.4030.124107450.00.9917.9700.115106453.52.29700.21370.30330.03410.0090590.9807.3770.125113450.20.9917.9560.115106453.02.41040.18940.29320.17090.0090600.9937.9560.115102453.30.9907.9560.115106453.12.28810.22070.30410.00620.0113For centered measures of both samples with respect to the centre ofF, the DD-diagram in the right subplot of Figure 4 shows the observations under the effect of a scale change only. It can be deduced that observations 3 and 50 are out of control because of a change in location only. When measures are centered, these two observations disappear from the out of control region in the right subplot of Figure 4. But, observations 2, 19, 27, 29, 33, 39, 43, 52, and 54 are out-of-control under the effect of a change in dispersion and of a location shift, respectively. These points are indicated with red stars in the DD-diagram in the right subplot of Figure 4.


## 5. Results and Concluding Remarks

The DD-diagram is a graphical comparison that exhibits location shifts and/or scale increase when moving from the distributionF of the reference sample to the distribution G of the empirical one. And to use this diagram, we do not need any requirement about the nature of the observed multivariate quality process distribution. Although, this procedure looks like a nonparametric method, DD-diagram does not require large samples. It suffices to have a size of the samples that goes beyond 30 to ensure a reasonable performance. So, whenever this size goes bigger, the DD-diagram improves in performance.The above application allows us to say that DD-diagram performs better thanT2 chart because its out-of-control signal does not depend on an error rate α as for the case of T2 control chart. Prior to the use of DD-diagram to monitor any observed process, the reference sample is tested if it is in control or not by using it twice as a reference sample and an empirical one in both x-axis and y-axis, respectively. In the above application, DD-diagram detects 11 points indicating that their components exceed their specified limits, whereas the T2 control chart gives only 5 points corresponding to an error rate α=5%. In order to detect the same points given by DD-diagram, the T2 control chart ought to use a type I error rate α=10%.When a multivariate quality process changes its distribution fromF to G and if the location shift is eliminated, that is, centering the measures with respect to the centre or the deepest point of F, DD-diagram makes it possible to distinguish between the out-of-control observations that were drifted because of location shifts and scale increase, respectively, and those that were drifted under the effect of variations in dispersion only. This fact is not feasible when using T2 chart.In general consider the test of a null hypothesis asserting stability of a production process versus an alternative one that concerns the existence of shifts in location and/or in scale then the empirical sample has higher dispersion than that of the reference one. This is deduced from the fact that the resulting clouds, of centered measures or not, are located under the limiting variation level line atLvalue. Then DD-diagram enables to present and detect graphically any out-of-control observation that the components exceed their specified limits. In addition, the DD-diagram sends out an out-of-control signal when the outlyingness of a point exceeds a specified value in all directions.



---

*Source: 830764-2012-01-17.xml*