Next Page Previous Page Home Tools & Aids Search Handbook
6. Process or Product Monitoring and Control
6.5. Tutorials
6.5.4. Elements of Multivariate Analysis
6.5.4.3. Hotelling's T squared

6.5.4.3.1.

T2 Chart for Subgroup Averages -- Phase I

Estimate mu with xbarbar Since mu is generally unknown, it is necessary to estimate mu analogous to the way that mu is estimated when an Xbar chart is used. Specifically, when there are rational subgroups, mu is estimated by Xdoublebar, with
    xbarbar = [xbarbar(1); xbarbar(2); ... ; xbarbar(p)]
Obtaining the xbarbari Each xbarbari, i = 1, 2, ..., p, is obtained the same way as with an Xbar chart, namely, by taking k subgroups of size n and computing
    xbarbar(i) = (1/k)*SUM[l=1 to k][xbar(il)].
Here xbarbar(il) is used to denote the average for the lth subgroup of the ith variable. That is,
    xbar(il) = SUM[r=1 to n][x(ilr)]
with xilr denoting the rth observation (out of n) for the ith variable in the lth subgroup.
Estimating the variances and covariances The variances and covariances are similarly averaged over the subgroups. Specifically, the sij elements of the variance-covariance matrix S are obtained as
    s(ij) = (1/k)*SUM[l=1 to k][s(ijl]
with sijl for i <> j denoting the sample covariance between variables Xi and Xj for the lth subgroup, and sij for i = j denotes the sample variance of Xi. The variances s(il)^2 (= siil) for subgroup l and for variables i = 1, 2, ..., p are computed as
    (1/(n-1))*SUM[r=1 to n][(x(ilr) - xbar(il))^2] .
Similarly, the covariances sijl between variables Xi and Xj for subgroup l are computed as
    (1/(n-1))*SUM[r=1 to n][(x(ilr) - xbar(il))*(x(jlr) - xbar(jl))] .
Compare T2 against control values As with an chart (or any other chart), the k subgroups would be tested for control by computing k values of T2 and comparing each against the UCL. If any value falls above the UCL (there is no lower control limit), the corresponding subgroup would be investigated.
Formula for plotted T2 values Thus, one would plot
    T(j)^2 = n*(xbar(j) - xbarbar)'*S(p)^(-1)*(xbar(j) - xbarbar)
for the jth subgroup (j = 1, 2, ..., k), with xbar denoting a vector with p elements that contains the subgroup averages for each of the p characteristics for the jth subgroup. (S(p)^(-1) is the inverse matrix of the "pooled" variance-covariance matrix, S(p), which is obtained by averaging the subgroup variance-covariance matrices over the k subgroups.)
Formula for the upper control limit Each of the k values of T(j)^2 given in the equation above would be compared with
    UCL = [(k*n*p - k*p - n*p + p)/(k*n - k - p + 1)]*F(alpha,p,k*n-k-p+1)
Lower control limits A lower control limit is generally not used in multivariate control chart applications, although some control chart methods do utilize a LCL. Although a small value for T(j)^2 might seem desirable, a value that is very small would likely indicate a problem of some type as we would not expect every element of xbar(j) to be virtually equal to every element in xbarbar.
Delete out-of-control points once cause discovered and corrected As with any Phase I control chart procedure, if there are any points that plot above the UCL and can be identified as corresponding to out-of-control conditions that have been corrected, the point(s) should be deleted and the UCL recomputed. The remaining points would then be compared with the new UCL and the process continued as long as necessary, remembering that points should be deleted only if their correspondence with out-of-control conditions can be identified and the cause(s) of the condition(s) were removed.
Home Tools & Aids Search Handbook Previous Page Next Page