Improvement in calibration

Abstract: An improved calibration method for calibrating the measured value of an analyte in a liquid sample by referring to a calibration curve drawn by plotting the interrelation between the calibration values of n ( n is an integer of not less than 2) calibrators containing different contents of the analyte and the measured values of the analyte contained in respective calibrators. In the improved method provided by the invention, the function of the calibration curve is determined by a first step of estimating an imaginal point, the imaginal point either corresponding to a calibration value larger than the uppermost value of the n calibrators and vicinal to the upper limit of the determination range of the analyte or corresponding to a calibration value smaller than the lowermost value of the n calibrators and vicinal to the lower limit of the determination range of the analyte or a calibration value vicinal to the zero value, and a second step of determining the function of the calibration curve by plotting the imaginal point with two to n measured values.

Description: The invention relates to a method for correcting measured values of an analyte in a liquid sample according to the preamble of claim 1. More generally, the present invention relates to an improvement in calibration of measured values obtained by analyzing liquid samples. More particularly, it provides a calibration method wherein the measured values distributed in a wide range can be calibrated or corrected even when standard samples or calibrators containing an analyte in the concentration of limited range can be prepared. The preamble of claim 1 is based on the conventional quantitative analysis of an analyte contained in a liquid sample, wherein the content of analyte in a test sample is determined by measuring an optical density (absorbancy or optical density of the transmitting or reflected light) of the liquid sample after it is subjected to a proper chemical or enzymatic reaction. In such a method, it is a common practice to determine the content of analyte by using a standard curve (commonly referred to as "calibration curve") which has been preliminarily drawn by plotting the interrelation between the known contents of analyte in the standard samples and the optical densities of the standard samples. However, a number of standard samples must be prepared and examined to draw a standard curve, and thus cumbersome and time-consuming operations are needed therefor. Particularly when a liquid sample obtained from a living body, such as whole blood, plasma or serum, is used as in the case of a clinical assay, it is difficult to prepare always many standard samples (body fluids) each containing a known content of an analyte. To obviate such a difficulty, it has been proposed and practised to prepare a tentative standard curve which is corrected or calibrated by finding the difference or error from the standard curve. In detail, only a few standard samples or calibrators are examined to find measured values which are used to correct the standard curve to prepare a calibration curve which is used for the determination of real measured values of respective samples. For example, calibrators containing standard values (calibration values; The term "calibration value" used throughout the specification means the given value for which calibration is to be made, and may be referred to as "real value" or "correct value") of H, M and L are subjected to the same assaying procedure to find optical densities or other proper parameters from which real measured values h , m and l are found by referring to the standard curve. The interrelation between the standard values H, M and L and the real measured values h , m and l are plotted to draw a calibration curve as shown in Fig. 2, and the measured value of each examined sample is corrected by using the thus drawn calibration curve. In an automated analyzer, the fucntion of the calibration curve (generally in the form of a quadratic equation) is determined from the data (l, L), (m, M) and (h, H) by the least squares method to obtain the following correction equation of: Y = α + βX + γX² Then, the measured value is corrected by using the thus obtained correction equation. When the differences in content between the standard values (calibration values) L, M and H are sufficiently large as is the case illustrated in Fig. 2, the effect of curvature of the calibration curve ( i.e. the contribution of the term of second order) is small even if the measured values l , m and h are varied independently, so that the measured values distributed in a wide range can be appropriately corrected by the use of the calibration curve. However, if the differences in content of the standard correct values L, M and H are not large enough, the correct value above the uppermost standard value H and the correct value below the lowermost standard value L tend to contain large errors. Particularly when the measured values l , m and h are varied in the manner as shown in Fig. 3, the calibration curve has an extreme value within the determination range to make it impossible to correct the measured values. Such a problem arises, for example, in a dry analysis method wherein a dry analysis element is used to analyze an analyte in a liquid sample taken from a living body fluid, such as whole blood, plasma or serum. The dry analysis element is generally composed of plural layers including, for example, a reagent layer, a porous spreading layer, etc., laminated sequentially on a transparent support. An aqueous liquid sample is spotted on the spreading layer to migrate into the reagent layer which is colored by the action of the analyte contained in the sample. The density of the colored reagent layer is determined by measuring the optical density of the reflected light, and the quantity or content of the analyte contained in the sample is determined by means of colorimetric method. Specific examples of dry analysis element are disclosed in US-A-2 846 808, US-A-3 036 893, US-A-3 368 872 and US-A-3 992 928, Unexamined Japanese Patent Publication Nos. 53888/1974 (corresponding to US-A-3 992 158), 164356/1980 (corresponding to US-A-4 292 272), 222769/1985 (corresponding to EP-A-0 162 302), 4959/1986 (corresponding to EP-A-0 166 365) and 90859/1980 (corresponding to US-A-4 258 001), Clinical Chemistry, 24 , 1335 - 1350 (1978), Analytical Chemistry, 55 (4), 498 - 514 (1983) and Clinical Chemistry, 27 , 1287 - 1290 (1981). The liquid sample spotted on the spreading layer spreads on the spreading layer to cover a generally circular zone having an area substantially in proportion to the volume of the spotted sample, and then liquid ingredients migrate into the reagent layer while the solid ingredients are filtered off. As a result, a substantially constant volume of aqueous liquid sample is fed to each unit area of the reagent layer. This function is known as a spreading function or metering function. However, the spreading function is significantly affected by the properties ( e.g. viscosity, specific gravity, pH, etc.) of the spotted aquaeous liquid sample. Accordingly, when an analyte contained in a body fluid, such as whole blood, plasma, serum or urine, is to be analyzed, it is preferable to use a calibrator which has properties resembling the body fluid taken from a living body. For example, when albumin or a total protein in blood is quantitatively analyzed, used calibrators are prepared by dissolving lyophilized human blood serum to have proper standard contents. However, the properties of the calibrators are significantly differentiated from those of the body fluids of natural origin as the content of albumin or total protein is varied in a wide range. For this reason, the difference between the lowermost content of analyte and the uppermost content of analyte in the calibrators cannot be set large enough to cover a wide determination range. As a result, when a calibration curve is drawn by using the least squares method, adequate calibration cannot be made within a wide range since an extreme value is found in the thus drawn calibration curve at the worst case. Furthermore, Patent Abstracts of Japan, Vol. 012390, Gr P772 of 18/10/1988 and JP-A-63-132 166 discloses a method for setting a standard curve by using a decreased number of measurements of standard samples as compared with previous methods, however, this document does not relate to a method for correcting measured values. Although the word "calibration curve" is mentioned in the English language abstract of this document, such "calibration curve" does not correspond to that of the present invention, but rather to a standard curve as mentioned at the beginning of the above discussion of the conventional quantitative analysis of an analyte contained in a liquid sample. According to the method of JP-A-63-132 166, the standard curve must be prepared in every assay process. Furthermore, in the method of JP-A-63-132 166 it is additionally necessary to determine the linear function of the locus of the maxima of the standard curve for every lot of the standard samples, since two of the essential constants used in the method of this document fluctuate in each lot of the standard samples or reagents used. An object of this invention is to provide an improved method for calibrating a measured value of an analyte contained in a liquid sample accurately within a wide range. A more specific object of this invention is to provide an improved method for calibrating an measured value of an analyte contained in a liquid sample by using a calibration curve which is drawn by plotting the interrelation between the calibration values of a few calibrators containing different contents of the analyte and the measured values of the analyte contained in respective calibrators, the calibration curve being modified by extrapolating an imaginal point. The object of the invention is achieved by the features of claim 1. The present invention thus provides an improved method in the calibration method for calibrating the measured value of an analyte in a liquid sample by referring to a calibration curve drawn by plotting the interrelation between the calibration values of n ( n is an integer of not less than 2) calibrators containing different contents of said analyte and the measured values of said analyte contained in respective calibrators, the improved method wherein the function of said calibration curve is determined by a first step of estimating an imaginal point, said imaginal point either corresponding to a calibration value larger than the uppermost value of said n calibrators and vicinal to the upper limit of the determination range of said analyte or corresponding to a calibration value smaller than the lowermost value of said n calibrators and vicinal to the lower limit of the determination range of said analyte or a calibration value vicinal to the zero value, and a second step of determining the function of said calibration curve by plotting said imaginal point with two to n measured values. In a preferred embodiment, said calibration value of said imaginal point is the upper limit of the determination range. In another preferred embodiment, said calibration value of said imaginal point is zero value. BRIEF DESCRIPTION OF THE DRAWINGS : The principle and a presently preferred embodiment of this invention will be described with reference to appended drawings in which: Fig. 1 is a graph illustrating the principle of this invention; Fig. 2 is a graph showing a calibration curve drawn by the conventional technology when calibrators covering a wide range can be prepared; Fig. 3 is a graph showing a calibration curve drawn by the conventional technology when calibrators covering a wide range cannot be prepared; Fig. 4 is a graph showing a calibration curve drawn by the method of this invention; and Fig. 5 is a graph showing a calibration curve drawn by a Comparative Example according to the conventional technology. DETAILED DESCRIPTION OF THE INVENTION : The principle of this invention will be initially described with reference to Fig. 1 for better understanding of the invention. Three calibrators containing known standard calibration values H, M and L of an analyte are assayed in the same process for the quantitative analysis of the analyte to find the measured values l , m and h , and the co-ordinates (H, h), (M, m) and (L, l) are plotted as shown in Fig. 1. The co-ordinates (imaginal calibration value, imaginal measured value) of an imaginal point is then taken. When it is desired to take an imaginal point above the uppermost point of the calibrators, the imaginal calibration value of the imaginal point is preferably close to the upper limit of the determination range, and more particularly it is preferred that the imaginal calibration value of the imaginal point is the upper limit of the determination range. When it is desired to take an imaginal point below the lowermost point of the calibrators, it is preferable that the imaginal calibration value of the imaginal point is close to zero value of the content of analyte or close to the lower limit of the determination range. In general, it is particularly preferred that the calibration value of the imaginal point is zero value when the subsequent mathematical processing is taken into account. The imaginal measured value corresponding to the imaginal calibration value will be referrd to as "estimated measured value" which will be estimated by the method described hereinafter. The imaginal measured value corresponding to the calibration value H max close to the upper limit of the determination range will be represented by h max , and the imaginal measured value at the zero value (correct value = 0) will be represented by C₀. The upper limit of the determination range means the upper limit value of the analyte within the determination range, and the zero value means that the content of analyte is zero. The imaginal measured values h max and C₀ are measured. When the calibration curve to be drawn is represented by a function of m order, (m+1) values including (l,L), ---, (m,M), ---, (h,H), are measured (or plural measured values may be obtained by repeating the analyzing operations for plural times for the three points (l,L), (m,M) and (h,H) depending on the order of the calibration curve) to determine the function of a curve of (m-1) order by the least squares method. From the thus determined function, values h max and C₀ are found. In a preferred embodiment of this invention, the calibration curve is represented by a function of second order. Generally in case where the calibration curve is represented by a function of second order, an equation of the first order, i.e. a function representing a linear line, is determined from three points (l,L), (m,M) and (h,H), and then the co-ordinates of the point (h max , H max ) may be fixed by finding the point at which the linear line intersects with the line represented by Y = H max and the co-ordinates of the point (C₀, 0) may be fixed by finding the point at which the linear line intersects with the X axis. Alternatively, the point (C₀, 0) is fixed by extending the line connecting the two points (l,L) and (m,M) and then finding the point at which the extension intersects with the X axis; and the point (h max , H max ) is fixed by extending the line connecting the two points (m,M) and (h,H) and then finding the point at which the extension intersects with the line represented by Y = H max . By connecting the five points (C₀, 0), (l,L), (m,M), (h,H) and (h max , H max ), a calibration curve (shown by the real line in Fig. 1) may be drawn. Measured values lying out of the range of from l to h can be corrected by using the thus drawn calibration curve. Two alternative methods for taking an imaginal point have been described. However, it has been found that the calibration values of the points beyond the lowermost and uppermost values of the calibrators are deviated to have errors when all of the n calibration values of the calibrators are used to draw the calibration curve due to the influence of the lowermost and uppermost calibration values of the calibrators. In order to eliminate such an error, when an imaginal point is taken in a range above the uppermost point of the calibrators, it is recommended to exclude the lowermost pair of the calibration value and measured value; and similarly the uppermost pair of the calibration value and measured value should be excluded when an imaginal point is taken in a range below the lowermost point of the calibrators. In other words, when n calibrators are used, the measured value of an imaginal point taken in a range above the uppermost measured value of the used n calibrators should be estimated by drawing a calibration curve by using 2 to (n-1) pairs of the calibration and measured values of the calibrators excluding the lowermost pair of the measured and calibration values; whereas the measured value of an imaginal point taken in a range below the lowermost value of the used n calibrators should be estimated by drawing a calibration curve using 2 to (n-1) pairs of the calibration and measured values of the calibrators excluding the uppermost pair of the measured and calibration values. When the imaginal measured value is estimated only by using the higher two points of the n calibrators, the imaginal measured value X max at the upper limit of the determination range may be calculated from the following equation of: wherein X max is the imaginal measured value at the upper limit of the determination range, Y max is the upper limit of the determination range, Y n is the calibration value of the calibrator containing maximum content of said analyte, Y n-1 is the correct value of the calibrator containing next to the maximum content of said analyte, X n is the measured value of the calibrator containing maximum content of said analyte, and X n-1 is the measured value of the calibrator containing next to the maximum content of said analyte. The measured value X₀ of the imaginal point having a calibration value of zero may be estimated by the following equation of: wherein X₀ is the imaginal measured value corresponding to the zero value of the calibration value, Y₁ is the calibration value of the calibrator containing the minimum content of said analyte, Y₂ is the calibration value of the calibrator containing next to the minimum content of said analyte, X₁ is the measured value of the calibrator containing the minimum content of said analyte, and X₂ is the measured value of the calibrator containing next to the minimum content of said analyte. EXAMPLES : A multi-layered analysis element for the quantitative analysis of albumin was corrected as follows while using Analyzer FDC 5000 which was an automated analyzer for analyzing ingredients contained in blood and produced and sold by Fuji Photo Film Co., Ltd. The used multi-layered analysis element was the one which is described in Example 1 of Unexamined Japanese Patent Publication No. 49962/1989 (corresponding to USSN 07/234,465 and EP-A-0304052). The multi-layered analysis element comprises a transparent support layer, a water-absorbing layer laminated on the support layer, and a spreading layer made of a knitted fabric and laminated on the water-absorbing layer, the spreading layer being impregnated with a reagent containing Bromocresol Green as an indicator to develop color by the presence of albumin. Calibrators used for the calibration were prepared as follows. The level L calibrator (Calibration Value: 2.4 g/dl) was prepared by re-dissolving a lyophilized product which had been produced by adding distilled water to the pooled human serum to prepare a solution having the desired content of albumin followed by lyophilization. The level M calibrator (Calibration Value: 3.4 g/dl) was prepared by re-dissolving a product which had been produced from the pooled human serum without any treatment followed by lyophilization. The level H calibrator (Calibration Value: 4.8 g/dl) was prepared by re-dissolving a product which had been produced by adding an ingredient to be analyzed to the pooled human blood serum followed by lyophilization. A standard curve drawn by plotting co-ordinates determined by using standard samples having known contents was stored in the Analyzer FDC 5000, and a content of the analyte ( i.e. measured value) in the examined sample was determined by measuring the optical density of the reflected light found in the multi-layered analysis element. Each of the level L, M and H calibrators was spotted on the analysis element and incubated at 37°C for 6 minutes, and then the optical density of the reflected light showing the color change in the analysis element was measured at a wavelength of 640 nm. The test was repeated for three times for each standard sample to obtain the data set forth below. Calibration Value Measured Value Mean Value Level L Calibrator 2.4 g/dl 2.3 g/dl 2.20 g/dl 2.2 g/dl 2.1 g/dl Level M Calibrator 3.4 g/dl 3.6 g/dl 3.73 g/dl 3.8 g/dl 3.8 g/dl Level H Calibrator 4.8 g/dl 4.6 g/dl 4.63 g/dl 4.8 g/dl 4.5 g/dl Then, the values C₀ and h max were calculated by using the following equation. Meantime, the upper limit H max of the analysis element was 6.0 g/dl. Imaginal Measured Value (C₀) at the zero point: wherein Y L is the calibration value of the level L calibrator; Y M is the calibration value of the level M calibrator; X ℓ is the measured value of the level L calibrator; and X m is the measured value of the level M calibrator. Imaginal Measured Value (h max ) at the upper limit of the determination range: wherein H max is the upper limit of the determination range; Y H is the calibration value of the level H calibrator; Y M is the calibration value of the level M calibrator; X h is the measured value of the level H calibrator; and X m is the measured value of the level M calibrator. The results of caluculation are as follows: C₀ = - 1.48 g/dl h max = 5.40 g/dl The function of the calibration curve was computed while using the thus calculated date to draw a calibration curve shown in Fig. 4. By plotting the measured values along the abscissa (X axis) and plotting the calibration values along the ordinate (Y axis), the function of the calibration curve was represented by the following equation. Y = 0.07608X² + 0.52366X + 0.68746 Correlation Coefficient γ = 0.98917 Deviation S yx = 0.25505 For the comparison purpose, the function of the calibration curve was computed while using only three points indicating the levels L, M and H to find that the function of the calibration curve was represented by the following equation. Y = 0.31172X² - 1.16337X + 3.4459 Correlation Coefficient γ = 0.98897 Deviation S yx = 0.16536 The calibration curve drawn only by using the three points indicating the level L, M and H has an extreme value as shown in Fig. 5. The measured values of various samples having different contents of the analyte were corrected by using the thus computed functions. The Bromocresol Green colorimetric method was used as the standard method for the determination of content of albumin. As will be seen from the Table set forth above, according to the Comparative Example, although effective correction can be made within the medium content range, errors induced by correction are too large within the low and high content ranges to indicate that pertinent correction cannot be made within a wide range. On the contrary, according to the Example of this invention, effective correction can be made within a wide range covering the range lower than the lowermost standard sample and higher than the uppermost standard sample. As will be appreciated from the foregoing, measured values distributed in a range out of the content range, within which standard calibrators can be prepared, can be accurately corrected according to this invention. Although the present invention has been described by referring to an embodiment in which a multi-layered analysis element is used, the invention is not limited thereto but may be applied for correcting the measured values within a wide range over which calibrators cannot be prepared due to problems caused by properties of the sample liquids. Accordingly, it is intended to include all such applications within the scope of the invention. A calibration curve is drawn, according to the present invention, by plotting points including an imaginal point defined by an imaginal calibration value and an estimated measured value corresponding to the imaginal calibration value, the imaginal calibration value being larger than the uppermost calibration value of the used calibrators or smaller than the lowermost calibration value of the used calibrators. The measured values can be calibrated accurately even if the calibration values corresponding to them are out of the range covered by the calibration values of the used calibrators.

Claims: A method for correcting measured values of an analyte in a liquid sample comprising the steps of: a) analyzing the sample to determine a measured value of the analyte; b) subjecting n (n is an integer and is equal to or greater than 2) calibrators containig varying known amounts of the analyte to assay in the same process for the quantitative analysis of the analyte so as to obtain the measured values of said n calibrators by using a pre-determined standard curve, said varying known amounts of analyte being referred as the calibration values of said calibrators; c) determining a function of a calibration curve by plotting the coordinates (calibration value, measured value) of said calibrators; d) comparing the measured value of the analyte to the calibration curve; and e) obtaining a corrected measured value of the analyte; characterized in that: said function of said calibration curve is determined by plotting at least one of a first and second imaginal points with the coordinates of said calibrators; wherein said first imaginal point has a coordinate of the first imaginal calibration value and the first imaginal measured value; said first imaginal calibration value being selected to be greater than the uppermost value of said calibration values of said n calibrators but not above the upper limit of a pre-determined determination range for said analyte; said first imaginal measured value being corresponding to said first imaginal calibration value and being estimated by referring a function of a first tentative calibration curve, which is determined by using two to (n-1) pairs of the calibration value and measured value of said n calibrators excluding a pair having the minimum of the calibration and measured value; and wherein said second imaginal point has a coordinate of the second imaginal calibration value and the second imaginal measured value; said second imaginal calibration value being selected to be i) less than the lowermost value of said calibration values of said n calibrators but not less than the lower limit of the pre-detemined determination range of said analyte; or ii) zero value; said second imaginal measured value being corresponding to said second imaginal calibration value and being estimated by referring a function of a second tentative calibration curve, which is determined by using two to (n-1) pairs of the calibration value and measured value of said n calibrators excluding a pair having the maximum of the calibration and measured value. The method of claim 1, wherein said function of said calibration curve is determined by plotting of said first imaginal point with the coordinates of said calibrators, and wherein said first imaginal calibration value is the upper limit correct value within the pre-determined determination range. The method of claim 2, wherein said first tentative calibration curve is linear, and wherein said first imaginal measured value is estimated by the following equation: wherein; X max is the first imaginal measured value at the upper limit of the determination range; Y max is the first imaginal calibration value which is the upper limit of the determination range; Y n is the calibration value of the calibrator containing the maximum content of said analyte; Y n-1 is the calibration value of the calibrator containing next to the maximum content of said analyte; X n is the measured value of the calibrator containing maximum content of said analyte; and X n-1 is the measured value of the calibrator containing next to the maximum content of said analyte. The method of claim 2 or 3, wherein said calibration curve is determined through the least squares method. The method of claim 2, wherein said liquid sample is analyzed by a dry analysis process. The method of claim 1, wherein said function of said calibration curve is determined by plotting of said second imaginal point with the coordinates of said calibrators. The method of claim 6, wherein said second imaginal calibration value is the zero value. The method of claim 7, wherein said second tentative calibration curve is linear, and wherein said second imaginal measured value is estimated by the following equation: wherein; X₀ is the second imaginal measured value corresponding to the zero value of the second tentative calibration curve; Y₁ is the calibration value of the calibrator containing the minimum content of said analyte; Y₂ is the calibration value of the calibrator containing next to the minimum content of said analyte; X₁ is the measured value of the calibrator containing the minimum content of said analyte; and X₂ is the measured value of the calibrator containing next to the minimum content of said analyte. The method of any of claims 6 to 8, wherein said calibration curve is determined through the least squares method. The method of claim 6, wherein said liquid sample is analyzed by a dry analysis process. The method of claim 1, wherein said function of said calibration curve is determined by plotting both of said first and second imaginal points with the coordinates of said calibrators. The method of claim 11, wherein said first imaginal calibration value is the upper limit correct value within the pre-determined determination range. The method of claim 12, wherein said first tentative calibration curve is linear, and wherein said first imaginal measured value is estimated by the following equation: wherein; X max is the first imaginal measured value at the upper limit of the determination range; Y max is the first imaginal calibration value which is the upper limit of the determination range; Y n is the calibration value of the calibrator containing the maximum content of said analyte; Y n-1 is the calibration value of the calibrator containing next to the maximum content of said analyte; X n is the measured value of the calibrator containing maximum content of said analyte; and X n-1 is the measured value of the calibrator containing next to the maximum content of said analyte. The method of claim 11, wherein said second imaginal calibration value is the zero value. The method of claim 14, wherein said second tentative calibration curve is linear, and wherein said second imaginal measured value is estimated by the following equation: wherein; X₀ is the second imaginal measured value corresponding to the zero value of the second tentative calibration curve; Y₁ is the calibration value of the calibrator containing the minimum content of said analyte; Y₂ is the calibration value of the calibrator containing next to the minimum content of said analyte; X₁ is the measured value of the calibrator containing the minimum content of said analyte; and X₂ is the measured value of the calibrator containing next to the minimum content of said analyte. The method of any of claims 11 to 15, wherein said calibration curve is determined through the least squares method. The method of claim 11, wherein said liquid sample is analyzed by a dry analysis process.

Inventors: MAKIUCHI HAJIME, TERASHIMA MASAAKI, MAKIUCHI, HAJIME, TERASHIMA, MASAAKI
Applicants: FUJI PHOTO FILM CO LTD, FUJI PHOTO FILM CO., LTD.
Application year: 1990
Main classification: G01N21
IPC codes: G01N21