Multivariate Control Charts: The Hotelling T2 Control Chart

October 2019

(Note: all the previous publications in the variable control charts category are listed on the right-hand side. Select “Publications” to go to the SPC Knowledge Base homepage. Selectthis linkfor information on the SPC for Excel software.)

A control chart normally monitors one variable over time. Perhaps this variable is machine uptime, a product characteristic, or on-time delivery. There are times, however, when the simultaneous monitoring of two or more related variables is important. The group of control charts that do this are called multivariate control charts. The most familiar one of these is the Hotelling T² control chart or just the T² control chart. This control chart is introduced in this publication.

In this issue:

• Introduction to the T² Control Chart
• Constructing a T² Control Chart
• Out of Control Points
• Summary
• Quick Links

You may download a pdf copy of this publication at this link. Please feel free to leave a comment at the end of the publication.

Introduction the T² Control Chart

In 1947, Harold Hotelling introduced a statistic which allowed multivariate observations to be plotted on a single chart. This statistic is now called Hotelling’s T² statistic. The statistic combines information from the mean as well as the dispersion of more than one variable. The calculations, which include some matrix algebra, are more difficult than those of “normal” control charts. This was a barrier to using multivariate control charts until software that could perform the calculations came along.

The T² control chart is used to detect shifts in the mean of more than one interrelated variable. The data can be in subgroups (like the X-R control chart) or the data can be individual observations (like in the X-mR control charts).

A few words of caution. The T² control chart, like other multivariate control charts, plots a value on the chart that you really can’t explain too well. Suppose you have two variables that are important in an adhesive process. The adhesive is characterized by two variables, pH and viscosity, which need to be controlled. Data for 20 batches are shown below for pH and viscosity. The data in this example is adapted from “Advanced Topics in Statistical Process Control” by Donald Wheeler (www.spcpress.com).

Table 1: pH and Viscosity Data

One way of monitoring the pH and the viscosity is by using a control chart for each variable. An individuals chart (X-mR) could be used for both pH and the viscosity. Figure 1 is the X chart for pH and Figure 2 is the X chart for viscosity. The moving range charts are not shown here.

Figure 1: X Chart for pH

Figure 2: X Chart for Viscosity

Figures 1 and 2 are both in statistical control. Look at the figures again. Notice anything? The two figures have very similar patterns. This implies that the two variables are correlated. A scatter diagram for the two variables is shown in Figure 3.

Figure 3: Scatter Diagram of pH vs Viscosity

The scatter diagram shows that the two variables are related. As pH increases, the viscosity tends to increase. Note the point at the bottom of the scatter diagram for a pH of 8.2. This data pair looks like it might be an outlier, despite the X charts being in control. This is where the T² control chart can be used. It is designed to see the impact of multiple variables at the same time.

Figure 4 is the T² control chart for the data in Table 1.

Figure 4: T² Control Chart for pH and Viscosity

The T² control chart shows an out of control point for batch 13. That out of control point does correspond to the low point on the scatter diagram above.

Look at the T² control chart in Figure 4. One of the first things you may notice about this control chart is that the value of T² has no resemblance to the original pH and viscosity data. So, looking at the value of T² really tells you nothing about the original data.

The only test for out of control points on this type of chart is points beyond the upper control limit (UCL). There is no lower control limit (LCL). There is one point beyond the control limit – but you don’t know which variable (pH or viscosity) caused the out of control point by looking at the control chart. Information on how to determine which variable(s) is responsible is given below.

Constructing a T² Control Chart

As stated before, the T² control chart can be used with data in subgroups or data that are individual observations. We will use individual observations to show how to construct a T² control chart. Note that the calculations are different for data in subgroups.

The first step in creating a T² control chart is to calculate the values of T². This is where matrix algebra comes in. The value of T² is given by:

T² = (x–x)^‘S^-1(x –x)

where xis the sample mean vector and S is the sample covariance vector. The bolded characters represent vectors. This is clearly now a little more complicated than the calculations for the basic control charts, but we will give you the general idea of how it is done for the case of two variables.

We will start with the(x–x)term in the T² equation. To create this term, you subtract the average value for the variable from each individual value for the variable. The average for pH and viscosity are given below.

xfor pH = 8.005

xfor viscosity = 4.9475

From Table 1, the first point for pH and viscosity are 7.75 and 5.48 respectively. Then, for the first point,

This will be eventually calculated for each point. Note that it is a matrix with one row for each point and a column for each variable.

Now,(x–x)’is the transpose of (x–x). The transpose of (x–x)is:

The transposed matrix has a column for each point and a row for each variable.

The S matrix is a little more difficult to find. It is found from the vector of moving differences for each variable. For each variable, v_i is found where

v_i = x_i+1 – x_i

So, for pH, v₁ is the difference between the second pH value and the first pH value:

v₁ = 8.50 – 7.75 = 0.75

This is done for both variables and the vector V contains the results for both variables:

where m = number of samples. The components of V are given below for the two variables.

So V has a column for each variable and m – 1 rows. S can now be found using the following:

where V’ is the transpose of V. You can easily multiply two matrices in Excel using the function MMULT. If you multiply V’ by V, you get the following:

S is then found by dividing each term in V’V by 2(m – 1) = 2(19) = 38.

Remember that T² is given by:

T²= (x–x)^‘S^-1(x–x)

S^-1 is the inverse of S. You can use the Excel function MINVERSE to find the inverse. Using that function gives:

Now we have all three terms:(x–x)^‘,S^-1, and(x–x)

So, for the first point:

T²= (x–x)^‘S^-1(x–x)

You can use the Excel MMULT function to multiply these three matrices together. The result is given below.

T² = 3.006896

This is the value of T² for the first data point. The above process can be used to generate the T² values for the rest of the data. The values of T² are shown in the table below.

Table 2: Values of T² for pH and Viscosity

These are the values that are plotted in the T² control chart in Figure 4. The other calculation that is needed is for the UCL. And, of course, this is a little more complicated here as well. The upper control limit based on individual observations is given by the following:

where m = number of samples, p = number of variables,? = Beta distribution,? = the confidence level and

You can use the BETAINV function in Excel to determine the value of the beta distribution. Note that the control limits do not depend on the values of T². Using? = 0.0027 gives an UCL of 12.6 as shown in Figure 4. (Note: if you use the BETAINV function,? is 1 – 0.0027.) There is no lower control limit.

Out of Control Points

When there is an out of control point on a T² control chart, more work is needed to determine which variable could be the cause of it. If there are points out of control, the results are decomposed to find out which variables could be responsible for the out of control situation. Each variable is removed from the calculations and the T² calculations repeated. The difference between the value of T² with the variable present and without the variable present is calculated. The larger the difference, the more likely the variable is to be the cause of the out of control point.

For the out of control point in Figure 4, the difference due to pH is 9.34 and the difference due to viscosity is 13.51. Based on this analysis, it appears that viscosity is the variable responsible for the out of control point.

Summary

This publication has introduced the T² control chart. This control chart is used to monitor multiple variables on one chart. This publication demonstrated how to do the calculations for the case of two variables. This included the calculation of T² as well as the upper control limit. It was shown how to determine which variable might be responsible for an out of control point. Computer software, like SPC for Excel, can easily handle the calculations presented in this publication.