Prev: Laguerre diagram
Next: SGE matlabpool problem
From: alkan alkaya on 19 Nov 2009 10:30
I want to use PCA (Principle component analysis) Q and T statistics test for fault detection control process. How can I use and how can I plot the Q or T2 statistics of the measured data?
From: Peter Perkins on 22 Nov 2009 10:49
alkan alkaya wrote:
> I want to use PCA (Principle component analysis) Q and T statistics test for fault detection control process. How can I use and how can I plot the Q or T2 statistics of the measured data?

The PRINCOMP function in the Statistics Toolbox returns the T^2 statistics; the PCARES function returns residuals that are used to compute the Q statistics. Hope this helps.
From: alkan alkaya on 23 Nov 2009 10:30
Peter Perkins <Peter.Perkins(a)MathRemoveThisWorks.com> wrote in message <hebmj7$m4f$1(a)fred.mathworks.com>...
> alkan alkaya wrote:
> > I want to use PCA (Principle component analysis) Q and T statistics test for fault detection control process. How can I use and how can I plot the Q or T2 statistics of the measured data?
>
> The PRINCOMP function in the Statistics Toolbox returns the T^2 statistics; the PCARES function returns residuals that are used to compute the Q statistics. Hope this helps.


Thank you peter perkins,

I could see the T2 plot by using PRINCOMP function but I couldn't plot the Q statistic. I have 1024*4 matrix 4= variables,1024= samples. When I use PCARES function I have RESIDUAL=R with same dimension (1024*4). for Q statistic as I know Q=R'*R after this calculation Q matrix occurs with 4*4 dimension but I think It should be 1*1024 dimension in order to plot Q test like T2 testing. I think I miss something. Thanks for your help.
From: Peter Perkins on 23 Nov 2009 16:02
alkan alkaya wrote:

> I could see the T2 plot by using PRINCOMP function but I couldn't plot the Q statistic. I have 1024*4 matrix 4= variables,1024= samples. When I use PCARES function I have RESIDUAL=R with same dimension (1024*4). for Q statistic as I know Q=R'*R after this calculation Q matrix occurs with 4*4 dimension but I think It should be 1*1024 dimension in order to plot Q test like T2 testing. I think I miss something. Thanks for your help.

I suspect what you want is diag(R*R'), or equivalently, sum(R.^2,2).
 | 
Pages: 1
Prev: Laguerre diagram
Next: SGE matlabpool problem