fig2 shows a hyperspatial cluster formed by a training set comprising 1000 spectra of three hypothetical randomly mixed compounds a , b , and c recorded at two wavelengths . in fig2 bivariate - normal noise with the variance of the horizontal axis set four times greater than that of the vertical axis has been added to the noiseless ( wavelength 1 , wavelength 2 ) reflectances of ( 5 , 5 ): ( 5 , 10 ); and ( 10 , 10 ). the distribution of points about the center of the cluster is clearly not symmetric . a cross - section of points through the mean or center of the cluster shows a definite right ( or positive ) skew , indicating a leverage effect by points with higher values . ( a frequency distribution polygon of such a cross - section appears in fig3 .) real near - ir reflectance data also show asymmetry ( see fig7 ). fig4 demonstrates the effect of such a skew on the determination of confidence limits . because of the balancing property of the mean , a small number of points located some distance away from the others exert a considerable leverage on the value of the mean . when these distant points are lower values than the rest as in fig4 the distribution is said to be left ( or negatively ) skewed , and the mean is shifted downward . in fig4 the outlying points should not be discarded because the solid line represents the real underlying distribution of the total population from which the samples were drawn . for this situation , the usual symmetric statement of confidence limits , e . g . : confidence limits for anything = mean +/- t a , ψ s anything does not provide an adequate description . see d . g . peters et al ., &# 34 ; chemical separations and measurements ,&# 34 ; p . 25 , w . b . saunders , philadelphia 1974 . for example , suppose the contour level of a distribution is given by a = z a σ ( see fig4 ) where z a = φ - 1 ( a ) and φ ( a )=( 2π ) - 1 / 2 - ∞ ∫ a e - t . spsp . 2 . sbsp ./ 2 dt . clearly , the probability of a point appearing at mean + a is different from the probability of a point appearing at mean - a . the beast ( bootstrap error - adjusted single - sample technique ) is a nonparametric method ( i . e ., a distribution - free method whose properties hold under very few assumptions about the population from which the data are obtained ) of assigning an asymmetric direction - dependent standard deviation to a multidimensional point estimate . the beast constructs a multi - dimensional cluster in hyperspace using the reflectances of each training - set sample at a predetermined number of wavelengths . new samples are projected into this space and a nonparametric confidence test is performed to determine whether the new sample is part of the training - set cluster . since it is possible that the distribution underlying a set of samples is skewed , an asymmetric nonparametric method of setting confidence limits is employed in accordance with the present invention . this method does not assume that the best point estimate of the mean lies at the center of the interval between the confidence limits . the asymmetric estimation capability is derived in part from the bootstrap , which can be summarized as follows . given an unknown distribution function f and some parameter of interest ( such as the mean or median ) that is a function of a number of independent identically distributed samples from f , the bootstrap determines the standard error in the parameter of interest from the samples in a training set which is a representative sampling of the unknown distribution . the real standard error of the parameter is a function of the unknown distribution , the size of the training sample set , and the form of the parameter . nevertheless , knowing the number of samples and the form of the parameter , the bootstrap estimate of the standard error is a function of only the unknown distribution f . in other words , although the actual distribution f is unknown , it can be estimated using the empirical probability distribution represented by the training set . the training set , composed of a number n of equally weighted samples from the unknown distribution , is created so that the empirical probability distribution adequately describes the variation in the unknown distribution f , and remains unchanged for the remainder of the estimation procedure . in one aspect of the invention , the training set is carefully constructed from unadulterated pharmaceutical samples in a fashion that captures all possible sample variations . for example , each sample member of the training set is illuminated with near - infrared light , and its reflectance is measured at a number d of wavelengths with a suitable infrared spectrometer . the spectra are recorded and transformed or mapped to produce a set or cluster of points in hyperspace . ( see fig5 step a ). bootstrap observations are randomly and independently drawn from the training set , with replacement from the training set so that the training set is never depleted or changed , to form a bootstrap set or replicate training set with the same total number as the training set : a number m of replicate training sets are created . ( see fig5 steps b and c ). the sampling distribution of the parameter of interest ( e . g ., the center of a group or cluster of spectral points in hyperspace ) for the unknown distribution f is approximated by the multidimensional distribution of that parameter in the m bootstrap or replicate training sets ( hereafter this distribution is termed simply the &# 34 ; bootstrap distribution &# 34 ;). in the present invention , the center of each replicate training set is determined by averaging the points obtained ( i . e ., by finding the group mean ), and the centers of the replicate training sets form a bootstrap distribution b ( see fig5 step d ). after the bootstrap distribution b is developed from the training set , the center cb of the bootstrap distribution is determined by averaging the centers . ( see fig5 step e ) the calculation of the bootstrap distribution is preferably arrived at by a monte carlo approximation through which a large number ( m ) of bootstrap sample sets is generated by randomly drawing observations from the training set ( again , with replacement from the training set ) to form bootstrap sets of the same size as the training set . the empirical distribution of the corresponding values of the parameter of interest calculated using the bootstrap sets is taken as an approximation to the bootstrap distribution , and its quantiles are used to estimate the value and precision of the parameter for the original population . the monte carlo method is by far the most common , but not the only , method for obtaining the bootstrap distribution . another method is direct theoretical calculation in special cases . taylor - series expansion methods can also be used to estimate the mean and variance of the bootstrap distribution of the parameter of interest . a training set of sample spectral values ( e . g ., reflectance , absorbance , etc .) recorded at d wavelengths from n training samples of a population can be represented by the n - by - d matrix t . ( generally , another n - by - d matrix v of spectra from validation samples is also assembled from the same source as the training set . the validation sample set v serves as an indicator of how well the training set describes the overall population variation .) the training set is carefully constructed from known samples ( samples that have been analyzed or identified by a suitable reference procedure ) that adequately describe all possible sample variations . this step is common to many spectral analysis procedures . a predetermined number m of bootstrap replications are calculated from t , and the result is the bootstrap distribution , the m - by - d array b . the calculation can be initiated by filling a matrix p with sample numbers to be used in bootstrap sample sets b . sub . ( s ) : where r is a random number between zero and unity . the values in p are scaled to the training - set size by : a bootstrap sample set b . sub . ( s ) is then created for each row i of the m - by - d bootstrap distribution b by : where k are the elements of the i - th rows of p and t are the elements of t . the q - th row of b is filled by the center of the q - th bootstrap sample set by : ## equ2 ## and the center of the bootstrap distribution is ## equ3 ## the foregoing operations can then be repeated using v . thus , a randomly selected set of samples ( containing the same number of elements as the training set t ) is drawn from the training set , with replacement from the training set , to form a bootstrap distribution by monte carlo approximation . it will be appreciated that a bootstrap distribution is used to estimate the true population distribution for a sample set . bootstrap distribution quantiles are readily converted into confidence intervals and are therefore useful in defining the boundaries of a training set in the hyperspace of spectral points , as described in co - pending u . s . patent application ser . no . 07 / 359 , 084 filed may 30 , 1989 , incorporated herein by reference . using the method described above in which the parameter of interest is the bootstrap - set center , selecting any two bootstrap - distribution percentiles gives the corresponding confidence limits for the center - parameter ( e . g ., selecting the 16th and 84th percentiles of the bootstrap distribution produces the central 68 % confidence limits ). this method fails to capture the asymmetry of a distribution when the median of the training set is noticeably different from the mean of the bootstrap distribution . in accordance with the present invention , a bias correction is provided in the calculation of the confidence limits that recognizes that the mean value of a distribution drifts in the direction of the skew ( measured with respect to the median of the distribution ). large differences between the median and mean suggest a skew that casts doubt upon the validity of a simple symmetric confidence interval . when no difference between the mean and the median exists , the value of the correction is zero . otherwise , new confidence limits are obtained with a improved representation of any skew present . as described above , each monitored wavelength is considered a dimension in a hyperspace , and each point in the hyperspace represents an entire spectrum , translated from the origin by amounts that correspond to the magnitude of the reflectance observed at each wavelength . the distribution of reflectances on each wavelength axis gives projections of the clusters of points . valid samples are defined as those that fall inside a cluster of training - set points , while false samples are those that fall outside the cluster . confidence limits , set along any linear combination of wavelengths ( dimensions ), define the surface of the cluster at a specified confidence level . these confidence limits are obtained by using the bootstrap to arrive at an estimate of the real - sample population distribution based upon the training - set distribution . the center of the real - sample distribution is estimated by the center of the bootstrap - set centers ( the bootstrap distribution ). when a new sample is tested , a vector or hyperline is formed in hyperspace between the new spectral point and the estimated center point of the real - sample distribution . a hypercylinder formed about this line contains a number of replicate - training - set centers . it will be appreciated that the hypercylinder is an approximation to the hyperline since hardly any replicate centers would be expected to fall on the hyperline . when the coordinates of these points are transformed into distances from the estimated center cb of the bootstrap distribution , a univariate distribution is constructed that is used in the confidence test of the present method . the reliance on nonparametric techniques produces a false - sample test that operates without assumptions about the shape , size , symmetry , or orientation of the spectral - point cluster in hyperspace . the spectrum x of a new sample ( such as an inhomogeneous or adulterated sample ) is recorded and projected as a point into the hyperspace . ( see fig6 step f ). a vector or hyperline is formed in hyperspace between the new hyperspatial point and the computed center point cb of the bootstrap distribution . ( see fig6 step g ). a hypercylinder is formed around this hyperline ( the radius r of the hypercylinder is typically 2 to 3 orders of magnitude smaller than its length ) containing a number of replicate training set centers from the bootstrap distribution b ( a minimum number of 50 points preferably should be used ) that are projected onto the hyperline by triangulation and transformed into distances from the estimated center cb to form a univariate distribution q . ( see fig6 step h ). the univariate distribution q is used to construct the confidence limits . the euclidean distances between the center cb of the bootstrap distribution and the points in the hypercylinder are sorted , and the upper u and lower l confidence limits are scaled in accordance with a measure of the skew of q . ( see fig6 step i ). the euclidean distance from the center cb to the new sample spectrum x is scaled by a skew - adjusted standard deviation ( sd ) calculated from the sorted distances . ( see fig6 step j ). it will be understood that the reproducibility of the training set points may be determined by obtaining the spectra of a second or validation set of samples having the same composition as the training set and transforming the spectra obtained into points in hyperspace . if the validation set hyperspatial points fall within the cluster of the training set , the training set cluster has been validated . the reliance on nonparametric techniques ( techniques that assume no particular underlying distribution ) produces an analytical method that functions without assumptions about the shape , size , symmetry or orientation of spectral - point clusters in hyperspace . this freedom is important because cluster characteristics have been shown to be unpredictable . thus , the present method can be used both for quality control applications ( i . e .. to identify empty capsules , or rust , metal shavings or other contaminants in capsules ) and for identification of adulterants placed in capsules ( i . e .. where a capsule has been the subject of deliberate tampering ). under such widely varying conditions , the unpredictability of the cluster characteristics makes the use of nonparametric techniques even more important . the calculation of the distance between cb and the new sample involves finding the hyperline connecting the center cb ( group mean ) of the bootstrap distribution b to the new or test sample spectrum x and determining the probability of x belonging to the training set cluster t based upon the number of points ( rows of the bootstrap distribution b ) within a predetermined distance r h of the hyperline ( in effect , taking the points or rows of b that fall within a hypercylinder of radius r h ). a preferable implementation proceeds toward the point - density of the hypercylinder by forming a plurality of planes connecting x , cb , and the rows of b . the use of a plurality of planes allows a complex structure like a hypercylinder in d - dimensional hyperspace to be represented in a simple manner regardless of the number of dimensions ( in fact , the number of planes is completely independent of the spatial dimension ). the 3 points that specify a plane in 3 - dimensional space also specify a triangle whose sides are readily determined from : ## equ4 ## where s . sub . ( 02 ) is the euclidean distance from cb to x , s . sub . ( or ) i is the euclidean distance from the i - th point or row of b to cb , and s . sub . ( 2r ) i is the euclidean distance from the i - th point of b to x . ( in general , a vector or array will be represented by an uppercase letter , and the vector &# 39 ; s elements will be represented by the lowercase letter thus , it will be understood that s . sub . ( or ) and s . sub . ( 2r ) are m - element arrays .) once the triangles have been formed , finding the rows of b , i . e ., the points of the bootstrap distribution , that fall inside the hypercylinder is now carried out by a procedure beginning : ## equ5 ## the elements of s . sub . ( p ) are euclidean distances , from the center cb to each point in the bootstrap distribution b , projected on the hyperline connecting cb to the new sample spectrum x . constructing a hypothetical hyperplane through cb such that the hyperline from x to cb is normal to the hyperplane allows the s . sub . ( p ) distances to be given a direction along the hyperline . points in the bootstrap distribution that are on the same side of the hyperplane as x are assigned positive distances in s . sub . ( p ). the remainder of the elements of s . sub . ( p ) are assigned negative values . this directional assignment is preferably accomplished by multiplying the elements s . sub . ( p ) i for which { s . sub . ( 02 ) 2 + s . sub . ( or ) i 2 & lt ; s . sub . ( 2r ) i 2 } by - 1 . the values of s . sub . ( p ) representing points in b that are outside the hypercylinder , i . e ., those points for which s . sub . ( r ) i & gt ; r h , are discarded for the remainder of the calculations , and the other values of s . sub . ( p ) are rank - ordered , i . e . : with n h the number of elements of s . sub . ( q ), i . e ., an array representation of the univariate distribution q . for a symmetric normal distribution , a 1 - sd contour is expressed as upper u and lower l confidence limits l =[ 0 . 16n h ] and u =[ 0 . 84n h ], where the square brackets indicate the greatest integer . thus , the symmetric confidence interval along the hyperline connecting x and cb is : { s . sub . ( q ) l & lt ; cb & lt ; s . sub . ( q ) u }. as described above , if n h is less than about 50 , the interval will not be very precise . the uncorrected standard deviation σ of the univariate distribution q can be found by either : ## equ6 ## or calculating the standard deviation from all of the points of s . sub . ( q ) and multiplying by n h 1 / 2 . once σ is known , the distance to the test spectrum in uncorrected sds ( suitable for unskewed training sets ) is simply : ## equ7 ## since many training sets are skewed , l and u are adjusted in accordance with the present invention to compensate for the skew before finding s . sub . ( q ) l and s . sub . ( q ) u . the number m of points in the bootstrap distribution b influences the selection of an adequate contour level . for m & lt ; 1000 the contour should probably be one sd , so a = φ (- 1 ); for m & gt ;& gt ; 1000 , e . g ., 10 5 or 10 6 , contours of two or three sds could be used . it will be understood that as m increases s (. sub . ( q ) &# 34 ; plateaus ,&# 34 ; or does not change , at increased contour size . setting z a = φ - 1 ( a ) and locating the median of the training set by c . sub . ( t ) j = m j ( t ij ), where t ij denotes the i - th sample of the training set , initiates the adjustment of the confidence limits to compensate for skew . c . sub . ( t ) tends to lie in hyperspace in the direction opposite to the direction of the skew ( with respect to cb ) because of the leverage effect of skewed points on the mean . this fact is the basis of the confidence - limit adjustment in accordance with the present invention , and the calculation of the adjustment continues with a determination of the distance and direction of the difference s . sub . ( cp ) between cb and c . sub . ( t ) projected on the hyperline connecting cb to x : ## equ8 ## the directional sign given to s . sub . ( cp ) is opposite that given to s . sub . ( p ) : if { s . sub . ( 02 ) 2 + s . sub . ( cor ) 2 & gt ; s . sub . ( c2r ) } then s . sub . ( cp ) is multiplied by - 1 . it will be understood that if the mean and median of s . sub . ( q ) are substantially different , s . sub . ( q ) may be skewed . since the central limit theorem applies to s . sub . ( q ), a skew probably indicates that n h points are not enough to create a stable confidence - limit adjustment . if skew is present , two options are available ( 1 ) use a larger number m of replicate training sets for the bootstrap distribution b ; or ( 2 ) increase r h and recalculate the univariate distribution ( note that this option may cause a loss of directional selectivity that can bias the quantiles of s . sub . ( q )). finally , it will be understood that s . sub . ( q ) is ordered so the use of some common ways of calculating m ( s . sub . ( q )) will result in very poor running times for the algorithms . to efficiently find m ( s . sub . ( q )), simply select the ( n h / 2 + 1 / 2 ) th element of s . sub . ( q ) when n h is odd , and the mean of the ( n h / 2 ) th and ( n h / 2 + 1 ) th elements when n h is even . in accordance with the invention , s . sub . ( cp ) is replaced by s . sub . ( cp ) δ + m ( s . sub . ( q )) to perform well as an adjustment in a computational environment where almost any axis scale or skew is possible . the addition of m ( s . sub . ( q )), the median of the univariate distribution q , helps to assure that the correction s . sub . ( cp ) and its array - analog z o ( defined below ) have the same sign ( direction ), and the factor δ provides a skew sensitivity adjustment . typically , δ has a value between zero and unity that is set empirically for each combination of t and b to keep the absolute magnitude of the adjustment inside the values of s . sub . ( q ). the calculation of the z o adjustment from s . sub . ( cp ) proceeds as follows : where [ r ( g ( x ))] denotes the greatest integer root of a function g ( x ) found by , e . g ., trapezoidal interpolation f i are the discrete samples of the distribution function f , i . e ., the n h elements of s . sub . ( q ) corrected using c t ; and i . sub . ( n ) is a vector of n h independent variables paired with f to locate the roots of f . if | 2z o |& gt ;| z a | then δ should be decreased and the calculation redone . otherwise , skew - corrected l and u values for s . sub . ( q ) are calculated from : as in the case of the uncorrected sd σ , the confidence interval along the hyperline connecting x and cb is { s . sub . ( q ) l . sbsb . c & lt ; cb & lt ; s . sub . ( q ) u . sbsb . c }. the upper confidence limit is preferably the one closest to the test spectrum x , and the skew - adjusted standard deviation σ c is simply s . sub . ( q ) u . sbsb . c . the distance in skew - adjusted sds from the center of the training set t to the test spectrum x is : ## equ9 ## the algorithms previously described were implemented in programs using vax - 11 basic language ( version 2 . 4 , digital equipment corporation ) and speakeasy iv delta ( vms version , speakeasy computing corporation , chicago ). these programs were run on vax - 11 / 780 and vax - 11 / 785 computers . spectral data at 18 discrete wavelengths were collected using a technicon instruments corp . infraalyzer 400 filter spectrophotometer . this spectrophotometer was directly connected to a vax - 11 / 780 computer using custom interface and graphics programs . three types of demonstrations were conducted , beginning with the common compound - identification problem and proceeding to the more difficult problem of detecting contaminated mixtures . the three tests were : 1 . the present method was used qualitatively to distinguish among four pure benzoic acid derivatives and aluminum oxide . each of the compound names in table i represents three spectra of that compound obtained by rotating the closed - sample cup in the spectrophotometer drawer and scanning through the 18 filters . 2 . in the second test , a training set was composed of 40 random mixtures of salicylic , benzoic , isophthalic , and p - aminobenzoic acids . a total of 80 samples were prepared for the training set in this test : of these , 20 were retained for use in test 3 below . the proportion of each of the four benzoic acid derivative components was varied from 0 to 25 % ( by weight ) in each sample , and the remainder was made up of aluminum - oxide as a diluent . a conventional random - mixing algorithm was used to generate the amounts to be used for each component in each sample , and each sample was examined at three wavelengths . the purpose of test 2 was to develop a large training - set spectral cluster that contained variations from a number of sources and determine whether pure samples could still be correctly identified as being different ( or false ) using only the training - set spectral data . the results for a few compounds from the laboratory shelf appear in table ii . 3 . the 40 - sample training set from test 2 was used again in conjunction with the 20 mixture samples that were held out . ten of those 20 samples were &# 34 ; contaminated &# 34 ; with a false - sample compound to determine whether those ten samples would be flagged . the other ten were unaltered and served as validation samples . fig8 shows a frequency distribution polygon for a cross - section through the center of the training set used for tables ii and iii , projected on a wavelength axis . the logarithm of the inverse reflectance value at the wavelength is on the abscissa , and the distribution is right - skewed . the benzoic acid derivatives used in test 1 and in the training set for tests 2 and 3 were analytical reagent - grade salicylic acid , p - aminobenzoic acid , isophthalic acid , and benzoic acid . reagent - grade aluminum oxide was used as a diluent in mixtures where a range of component concentrations was desired ( tests 2 and 3 ). before each sample ( or sample mixture ) was read in the infraalyzer 400 spectrophotometer , the sample was ground and mixed in a spex mixer / mill . the powder was then sifted through a 100 - mesh sieve and packed into the closed sample cup provided with the spectrophotometer . as described above , three readings were taken on each sample , each successive reading after a 120 - degree rotation of the closed cup . the ability of the present method to differentiate among four similar benzoic acid derivatives ( test 1 ) is demonstrated by the data in table i . the pure - compound spectra form four clusters in hyperspace whose intercluster distances ( expressed in skew - adjusted standard deviations ) appear in the table . the distances tabulated have units of asymmetric nonparametric central 68 % confidence intervals ( equivalent to one standard deviation if the underlying distribution were gaussian ) determined in accordance with the present invention . the intervals were calculated from the distance ( euclidean metric ) between the 0 . 16 and 0 . 84 quantiles of the bootstrap distribution , following the projection of this distribution onto the hyperline connecting the center of the bootstrap distribution and the new - sample spectral point . the compound heading each column was designated as the training set , so that the distance from a column - heading compound to a row - heading compound is given in terms of nonparametric standard deviations ( sds ) of the column - heading compound . for example , the distance between the clusters representing benzoic acid and salicylic acid can be expressed in two ways : in terms of the standard deviation of benzoic acid in the direction of salicylic acid ( 39 sds ), or in terms of the standard deviation of salicylic acid in the direction of benzoic acid ( 137 sds ). the difference in sds , of course , reflects the difference in the variances ( sizes ) of the benzoic and salicylic acid clusters in the direction of each other . benzoic acid exhibited the largest overall variance of any of the compounds tested because it was difficult to grind . incompletely ground crystals could slip through a 100 - mesh sieve and thereby introduce an additional source of variation into the benzoic acid data cluster . the variance in the other compounds can be attributed primarily to orientation effects caused by packing peculiarities , because this was the only factor that was varied between replicate measurements of these substances . aluminum oxide ( included in this test because it was used as the diluent in tests 2 and 3 ) had the smallest overall variance probably because it was the only reagent that was available in a powder fine enough to pass through the 100 - mesh sieve . other distances in table i ar also in accordance with what might be predicted aluminum oxide and isophthalic acid are 4208 sds apart because chemically they have little in common . in fact , aluminum oxide is the most distinct from the other compounds benzoic acid and salicylic acid form the closest spectral clusters since they differ by only a single oxygen atom . for the pure components , only a small number of bootstrap replications are required ( m = 50 ) because the variance of the training set spectra is relatively small , and the distance between adjacent clusters ( measured in 18 - dimensional space ) is relatively large . this clear disparity allows the samples to be distinguished computationally in a fraction of a second . considering the distances ( in terms of sds ) involved , there would be little danger of misidentifying the pure compounds . a large number of sample components , varying over a relatively broad range of concentrations , can have the effect of &# 34 ; filling &# 34 ; the spectral hyperspace provided by only a few analytical wavelengths . this phenomenon is demonstrated by the 3 - wavelength training set used in obtaining table ii . the euclidean distances ( not shown ) for the data in tables i and ii were within an order of magnitude or so of each other ; however , the distances in terms of standard deviation units are far different . the reason that the distances in sds shrink so much from table i to table ii is that the cluster size in table ii increased by about a factor of 100 . this spread causes some apparently unrelated compounds , like sucrose and whole wheat flour , to appear to be similar to the training set of benzoic acid derivatives . the increase in training - set cluster size also highlights the need to be aware of the exact shape of the training - set cluster when small distances are to be used for identifying valid or false samples . finally , mixing a number of different components in each sample dilutes the contribution of each component to the total sample spectrum , in effect compressing the available analytical space and increasing the probability of mixture - clusters overlapping . ( this behavior can be seen in fig2 .) the third test of the beast investigated a worst - case mixture - overlap situation by employing the same training set described in test 2 to train the beast algorithm . the algorithm was then presented with the remaining 20 sample mixtures . ten of these 20 mixtures were &# 34 ; contaminated &# 34 ; with acetylsalicylic acid ( randomly varying over a range from 1 to 20 %). the results appear in table iii . even in this worst - case example it is apparent that the contaminated samples were likely to be detected as false . three out of ten of the contaminated samples failed a 3 - sd test , and none of the uncontaminated samples were incorrectly identified . one might wonder why the distance ( in terms of sd ) for similar contaminant concentrations varies so widely among contaminants in table iii and why contaminant concentration does not appear to correlate well to distance . when components similar to acetylsalicylic acid rise in concentration , a 1 % contribution to the sample spectrum from acetylsalicylic acid becomes relatively smaller , and the sample spectral point appears to move closer to the training set . table iv shows the distance response of the beast for two groups of similar acetylsalicylic acid contaminant concentrations . the concentration of the diluent ( aluminum oxide ) is inversely related to the concentrations of the benzoic acid derivatives , so lower diluent concentrations indicate an increase in the concentrations of the other ( non - contaminant ) compounds . the relative increase in the concentrations of the non - contaminant compounds , in turn , correlates to a decrease in the distance of the sample from the training set . table i______________________________________using the beast for qualitative analysisdifferentiating among four benzoic - acidderivatives and a diluent ( distances inskew - adjusted sds . sup . a ) isoh2p . sup . b benzo . sup . c al . sub . 2 o . sub . 3 . sup . d salcyl . sup . e paba . sup . f______________________________________isoh2p 0 188 4208 638 1072benzo 1254 0 2889 137 483al . sub . 2 o . sub . 3 3930 405 0 1488 3041salcyl 1197 39 2991 0 449paba 1052 72 3198 234 0______________________________________ . sup . a the distance between any two compoundcluster combinations is given in terms of the columnheading compound in the direction of the rowheading compound . . sup . b isophthalic acid . sup . c benzoic acid . sup . d aluminum oxide ( the diluent in later tables ) . sup . e salicylic acid . sup . f paminobenzoic acid table ii______________________________________distance ( in skew - adjusted sds ) of fourpure ground ( 100 mesh ) compounds from atraining set . sup . a of mixtures of compoundscompound distance______________________________________acetylsalicylic acid 8 . 20dextrose 46 . 33whole wheat flour 12 . 71sucrose 3 . 88______________________________________ . sup . a training set composed of mixtures of benzoic acid , isophthalic acid , salicylic acid , paminobenzoic acid , and aluminum oxide diluent . table iii______________________________________distance ( in skew - adjusted sds ) of realand false ( contaminated ) mixtures fromthe mixture training set used in table iisample real samples . sup . a false samples . sup . b contaminantno . ( validation set ) ( contaminated set ) (%) ______________________________________1 1 . 29 6 . 59 1 . 12 . 61 3 . 25 8 . 03 . 76 5 . 05 3 . 84 . 47 4 . 55 14 . 85 1 . 97 2 . 66 4 . 26 1 . 54 1 . 77 19 . 97 1 . 65 6 . 53 1 . 78 1 . 18 3 . 34 1 . 09 . 79 3 . 88 10 . 710 1 . 10 2 . 32 4 . 5______________________________________ ( note : the cluster surface is customarily defined as being three standard deviations from the center of the cluster .) . sup . a &# 34 ; real &# 34 ; mixtures contain benzoic acid , isophthalic acid , salicylic acid , paminobenzoic acid , and aluminum oxide ; the real samples form a validation set because they contain the same compounds as the training se samples . . sup . b &# 34 ; false &# 34 ; mixtures contain the same components as real mixtures and acetylsalicylic acid . table iv______________________________________distance response of the beast ( in skew - adjustedstandard deviations ) for two groups of similar contaminantconcentrations in the false samples for table iiisample acetylsalicylicno . distance ( sds ) diluent (%) acid (%) ______________________________________1 . 6 . 59 67 1 . 17 . 6 . 53 64 1 . 78 . 3 . 34 50 1 . 03 . 5 . 05 57 3 . 85 . 2 . 66 41 4 . 210 . 2 . 32 44 4 . 5______________________________________ the invention has been described illustratively , not restrictively , and those of ordinary skill in the art will recognize modifications that are intended to be included within the spirit and scope of the invention that is to be delimited solely by the appended claims .