Jmlimit
Contents |
Purpose
Confidence limits for Q residuals via Jackson-Mudholkar.
Synopsis
- rescl = jmlimit(pc,s,cl)
- cl = jmlimit(pc,s,Q,2);
Description
JMLIMIT estimates confidence limits for Q residuals based on the Jackson-Mudholkar method. See Jackson, J.E., "A User's Guide to Principal Components", John Wiley & Sons, New York, NY (1991), and the discussion in the Chemometrics Tutorial on PCA.
Inputs are the number of PCs used pc, the vector of eigenvalues s, and the confidence limit cl expressed as a fraction (e.g. 0.95). Note that for a PCA model structure, model, that the eigenvalues can be found in model.detail.ssq(:,2).
The output rescl is the confidence limit based on the method of Jackson and Mudholkar. See chilimit for an alternate method of residual limit calculation based on chi squared.
Do inverse calculation with flag=2:cl = jmlimit(pc,s,Q,2);where Q is the sum of squares residuals from a model and "2" is a flag indicating the inverse calculation
Examples
rescl = jmlimit(2,ssq(:,2),0.95);
For a PCA model contained in the structure model:
rescl = jmlimit(4,model.detail.ssq(:,2),0.99);