T-Squared Q residuals and Contributions: Difference between revisions

From Eigenvector Research Documentation Wiki
Jump to navigation Jump to search
imported>Jeremy
No edit summary
 
(8 intermediate revisions by 2 users not shown)
Line 1: Line 1:
Hotelling's T<sup>2</sup> (T-Squared), Q residuals, and their corresponding contributions, are summary statistics which help describe how well a model is describing the data in a given sample, and why that sample has its observed scores for a given model. These statistics are generally only available for factor-based models like [[Pca|Principal Components Analysis]], [[Pls| Partial Least Squares]], [[Pcr|Principal Components Regression]], or [[Mcr| Multivariate Curve Resolution]] models. But some non-factor-based models (e.g. [[Mlr|Multiple Linear Regression]]) also provide similar statistics.
__TOC__
==Introduction==
 
Hotelling's T<sup>2</sup> (T-Squared), Q residuals, and their corresponding contributions, are summary statistics which help explain how well a model is describing a given sample, and why that sample has its observed scores in a given model.
 
The statistics are described for factor-based models like [[Pca|Principal Components Analysis]] and [[Pcr|Principal Components Regression]]. It should be noted that similar statistics can be calculated for [[Pls| Partial Least Squares]], nonfactor-based models (e.g. [[Mlr|Multiple Linear Regression]]) and non-orthogonal models such as classical least squares but these are not shown here.


Given a model which describes data generically using loadings (<b>P</b>) and scores (<b>T</b>), the following equation describes the relationship between the data and the model:
Given a model which describes data generically using loadings (<b>P</b>) and scores (<b>T</b>), the following equation describes the relationship between the data and the model:
Line 6: Line 11:


Where X is a matrix of ''m'' rows (samples) and ''n'' columns (variables). Based on this equation, the Q residuals, Hotelling's T<sup>2</sup>, and contributions for each can be calculated. This page describes the calculation, theory, and software how-to for calculating and using these statistics.
Where X is a matrix of ''m'' rows (samples) and ''n'' columns (variables). Based on this equation, the Q residuals, Hotelling's T<sup>2</sup>, and contributions for each can be calculated. This page describes the calculation, theory, and software how-to for calculating and using these statistics.
The following discussion is generally written with respect to a Principal Components Analysis (PCA) model. When the statistics are calculated for other models, the calculations are somewhat different, but their interpretation is the same.


==Q Residuals==
==Q Residuals==
Line 14: Line 21:


where <math>e_i</math> is the i<sup>th</sup> row of <math>E</math>, <math>P_k</math> is the matrix of the k loadings vectors retained in the model (where each vector is a column of <math>P_k</math>) and I is the identity matrix of appropriate size (n by n). The Q statistic indicates how well each sample conforms to the model. It is a measure of the difference, or residual, between a sample and its projection into the k factors retained in the model.
where <math>e_i</math> is the i<sup>th</sup> row of <math>E</math>, <math>P_k</math> is the matrix of the k loadings vectors retained in the model (where each vector is a column of <math>P_k</math>) and I is the identity matrix of appropriate size (n by n). The Q statistic indicates how well each sample conforms to the model. It is a measure of the difference, or residual, between a sample and its projection into the k factors retained in the model.
''Note:'' For models in which the <math>P</math> loadings are not orthogonal, the second part of the above equation cannot be used. Instead, an orthogonalization must be used.


===Q Contributions===
===Q Contributions===
Line 25: Line 34:
:* Right-clicking the figure will make the figure "interactive" allowing you to view different subsets of samples when more than one has been selected.
:* Right-clicking the figure will make the figure "interactive" allowing you to view different subsets of samples when more than one has been selected.


====Relative Q Contributions====
===Relative Q Contributions===


Another form of contributions are "relative" Q contributions. Normally, the Q contributions are calculated as the residuals relative to the model. In some cases, comparing the Q contributions of two or more samples relative to each other can be useful in determining whether the same measurement phenomenon is responsible for both samples' Q residuals. In this case, relative contributions can be used in which one or more samples are selected as a "reference" set, then Q contributions
Another form of contributions are "relative" Q contributions. Normally, the Q contributions are calculated as the residuals relative to the model. In some cases, comparing the Q contributions of two or more samples relative to each other can be useful in determining whether the same measurement phenomenon is responsible for both samples' Q residuals. In this case, relative contributions can be used in which one or more samples are selected as a "reference" set, then Q contributions are calculated for a sample of interest and the two are subtracted, indicating what is different in their residuals.


Relative Q contributions can be viewed when viewing a Scores plot. On the Plot Controls, select the "Q con Ref." button and, when prompted, select the sample(s) you want to use as the reference. When complete, a new class set will be created showing you which sample(s) are being used for the relative contributions. Next, use the "Q con" button to select the samples for which you want to view contributions. The result will be a contributions plot with relative contributions calculated and noted on the Y-axis.
Relative Q contributions can be plotted when viewing a Scores plot. On the Plot Controls, select the "Q con Ref." button and, when prompted, select the sample(s) you want to use as the reference. When complete, a new class set will be created showing you which sample(s) are being used for the relative contributions. Next, use the "Q con" button to select the samples for which you want to view contributions. The result will be a contributions plot with relative contributions calculated and noted on the Y-axis.


==Hotelling's T-Squared==
==Hotelling's T-Squared==
Line 40: Line 49:
::<math>T_i^2 = t_i\lambda^{-1}t_i^T = x_iP_k\lambda^{-1}P_k^Tx_i^T</math>
::<math>T_i^2 = t_i\lambda^{-1}t_i^T = x_iP_k\lambda^{-1}P_k^Tx_i^T</math>


where <math>t_i</math> refers to the i<sup>th</sup> row of <math>T_k</math>, the ''m'' by ''k'' matrix of scores from the model, and <math>\lambda</math> is a diagonal matrix containing the eigenvalues (<math>\lambda_1</math> through <math>\lambda_k</math>) corresponding to ''k'' principal components, factors, or latent variables retained in the model.
where <math>t_i</math> refers to the i<sup>th</sup> row of <math>T_k</math>, the ''m'' by ''k'' matrix of scores from the model, and <math>\lambda</math> is a diagonal matrix containing the eigenvalues (<math>\lambda_1</math> through <math>\lambda_k</math>) corresponding to ''k'' principal components, factors, or latent variables retained in the model where <math>\lambda = T_k^TT_k/(m-1)</math> accounts for scaling by <math>(m-1)</math>.
 
''Note:'' The equations are specific to orthogonal loadings <math>P</math> and diagonal <math>\lambda</math>.


===T-Squared Contributions===
===T-Squared Contributions===
T<sup>2</sup> contributions describe how individual variables contribute to the Hotelling's T<sup>2</sup> value for a given sample. The contributions to <math>T_i^2</math> for the i<sup>th</sup> sample, <math>t_{con,i}</math>, is a 1 by ''k'' vector calculated from:
T<sup>2</sup> contributions (a.k.a. T Contributions) describe how individual variables contribute to the Hotelling's T<sup>2</sup> value for a given sample. The contributions to <math>T_i^2</math> for the i<sup>th</sup> sample, <math>t_{con,i}</math>, is a 1 by ''k'' vector calculated from:


::<math>t_{con,i} = t_i\lambda^{-1/2}P_k^T = x_iP_k\lambda^{-1/2}P_k^T</math>
::<math>t_{con,i} = t_i\lambda^{-1/2}P_k^T = x_iP_k\lambda^{-1/2}P_k^T</math>
Line 50: Line 61:
data are scaled to equalize the variance captured by each factor. Note that this formulation for T<sup>2</sup> contributions has the property that the sum-squared contributions give <math>T_i^2</math> for the given sample. That is:
data are scaled to equalize the variance captured by each factor. Note that this formulation for T<sup>2</sup> contributions has the property that the sum-squared contributions give <math>T_i^2</math> for the given sample. That is:


::<math>T_i^2 = t_{con,i}t_{con,i}^T
::<math>T_i^2 = t_{con,i}t_{con,i}^T</math>
 
T<sup>2</sup> contributions are plotted while viewing a [[Model_Building:_Plotting_Scores|scores plot]]. Use the "T con" button on the Plot Controls to request T<sup>2</sup> contributions. You will be prompted to select one or more samples for which you wish to see contributions. Select the samples and a new figure will be created showing the contributions.
 
:'''Hints:'''
:* Left-click on a variable in the figure to view the raw data for that variable (column of X).
:* Right-clicking the figure will make the figure "interactive" allowing you to view different subsets of samples when more than one has been selected.
 
===Relative T-Squared Contributions===
 
Normally, T<sup>2</sup> contributions are calculated relative to the center of the model which is often the average sample. They describe how the original variables combined to give the scores that were observed for a given sample. In some cases, it is of more interest to calculate what makes a sample different from the other samples ''within the model''. That is the utility of relative T<sup>2</sup> contributions. Relative contributions show how two or more samples differ from each other. Put another way, relative T<sup>2</sup> contributions show what changes in the underlying X variables to cause a sample to move arbitrarily from one point to another the "score space" of the model.
 
Numerically, relative contributions are simply the T<sup>2</sup> contributions for a reference sample (or the mean contributions from a group of reference samples) subtracted from the T<sup>2</sup> contributions for a sample of interest.
 
Relative T<sup>2</sup> contributions can be plotted when viewing a Scores plot. On the Plot Controls, select the "T con Ref." button and, when prompted, select the sample(s) you want to use as the reference. When complete, a new class set will be created showing you which sample(s) are being used for the relative contributions. Next, use the "T con" button to select the samples for which you want to view contributions. The result will be a contributions plot with relative contributions calculated and noted on the Y-axis.

Latest revision as of 08:51, 4 August 2022

Introduction

Hotelling's T2 (T-Squared), Q residuals, and their corresponding contributions, are summary statistics which help explain how well a model is describing a given sample, and why that sample has its observed scores in a given model.

The statistics are described for factor-based models like Principal Components Analysis and Principal Components Regression. It should be noted that similar statistics can be calculated for Partial Least Squares, nonfactor-based models (e.g. Multiple Linear Regression) and non-orthogonal models such as classical least squares but these are not shown here.

Given a model which describes data generically using loadings (P) and scores (T), the following equation describes the relationship between the data and the model:

Where X is a matrix of m rows (samples) and n columns (variables). Based on this equation, the Q residuals, Hotelling's T2, and contributions for each can be calculated. This page describes the calculation, theory, and software how-to for calculating and using these statistics.

The following discussion is generally written with respect to a Principal Components Analysis (PCA) model. When the statistics are calculated for other models, the calculations are somewhat different, but their interpretation is the same.

Q Residuals

Q Residuals are a lack-of-fit statistic calculated as the sum of squares of each row (sample) of E; i.e., for the ith sample in , :

where is the ith row of , is the matrix of the k loadings vectors retained in the model (where each vector is a column of ) and I is the identity matrix of appropriate size (n by n). The Q statistic indicates how well each sample conforms to the model. It is a measure of the difference, or residual, between a sample and its projection into the k factors retained in the model.

Note: For models in which the loadings are not orthogonal, the second part of the above equation cannot be used. Instead, an orthogonalization must be used.

Q Contributions

A single row of the m by n E matrix, , represents the Q contributions for a given sample. The Q contributions show how much each variable contributes to the overall Q statistic for the sample while retaining the sign of the deviation in the given variable. Such contributions can be useful in identifying the variables which contribute most to a given sample's sum-squared residual error. Such residuals can often be used to determine whether the lack-of-fit is due to systematic deviation of that sample or simple random variation.

Q contributions are plotted while viewing a scores plot. Use the "Q con" button on the Plot Controls to request Q contributions. You will be prompted to select one or more samples for which you wish to see contributions. Select the samples and a new figure will be created showing the contributions.

Hints:
  • Left-click on a variable in the figure to view the raw data for that variable (column of X).
  • Right-clicking the figure will make the figure "interactive" allowing you to view different subsets of samples when more than one has been selected.

Relative Q Contributions

Another form of contributions are "relative" Q contributions. Normally, the Q contributions are calculated as the residuals relative to the model. In some cases, comparing the Q contributions of two or more samples relative to each other can be useful in determining whether the same measurement phenomenon is responsible for both samples' Q residuals. In this case, relative contributions can be used in which one or more samples are selected as a "reference" set, then Q contributions are calculated for a sample of interest and the two are subtracted, indicating what is different in their residuals.

Relative Q contributions can be plotted when viewing a Scores plot. On the Plot Controls, select the "Q con Ref." button and, when prompted, select the sample(s) you want to use as the reference. When complete, a new class set will be created showing you which sample(s) are being used for the relative contributions. Next, use the "Q con" button to select the samples for which you want to view contributions. The result will be a contributions plot with relative contributions calculated and noted on the Y-axis.

Hotelling's T-Squared

Whereas Q residuals represent the magnitude of the variation remaining in each sample after projection through the model, the Hotelling's T2 values represent a measure of the variation in each sample within the model. It indicates how far each sample is from the center (scores = 0) of the model.

Hotelling’s T2 is the sum of the normalized squared scores and is defined as:

where refers to the ith row of , the m by k matrix of scores from the model, and is a diagonal matrix containing the eigenvalues ( through ) corresponding to k principal components, factors, or latent variables retained in the model where accounts for scaling by .

Note: The equations are specific to orthogonal loadings and diagonal .

T-Squared Contributions

T2 contributions (a.k.a. T Contributions) describe how individual variables contribute to the Hotelling's T2 value for a given sample. The contributions to for the ith sample, , is a 1 by k vector calculated from:

The term can be considered a scaled version of the data within the model. The data are scaled to equalize the variance captured by each factor. Note that this formulation for T2 contributions has the property that the sum-squared contributions give for the given sample. That is:

T2 contributions are plotted while viewing a scores plot. Use the "T con" button on the Plot Controls to request T2 contributions. You will be prompted to select one or more samples for which you wish to see contributions. Select the samples and a new figure will be created showing the contributions.

Hints:
  • Left-click on a variable in the figure to view the raw data for that variable (column of X).
  • Right-clicking the figure will make the figure "interactive" allowing you to view different subsets of samples when more than one has been selected.

Relative T-Squared Contributions

Normally, T2 contributions are calculated relative to the center of the model which is often the average sample. They describe how the original variables combined to give the scores that were observed for a given sample. In some cases, it is of more interest to calculate what makes a sample different from the other samples within the model. That is the utility of relative T2 contributions. Relative contributions show how two or more samples differ from each other. Put another way, relative T2 contributions show what changes in the underlying X variables to cause a sample to move arbitrarily from one point to another the "score space" of the model.

Numerically, relative contributions are simply the T2 contributions for a reference sample (or the mean contributions from a group of reference samples) subtracted from the T2 contributions for a sample of interest.

Relative T2 contributions can be plotted when viewing a Scores plot. On the Plot Controls, select the "T con Ref." button and, when prompted, select the sample(s) you want to use as the reference. When complete, a new class set will be created showing you which sample(s) are being used for the relative contributions. Next, use the "T con" button to select the samples for which you want to view contributions. The result will be a contributions plot with relative contributions calculated and noted on the Y-axis.