US20080249717A1 - Method and apparatus for chromatographic data processing - Google Patents

Method and apparatus for chromatographic data processing Download PDF

Info

Publication number
US20080249717A1
US20080249717A1 US12/155,812 US15581208A US2008249717A1 US 20080249717 A1 US20080249717 A1 US 20080249717A1 US 15581208 A US15581208 A US 15581208A US 2008249717 A1 US2008249717 A1 US 2008249717A1
Authority
US
United States
Prior art keywords
unit
peaks
chromatogram
peak
fitting
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US12/155,812
Inventor
Masahito Ito
Kisaburo Deguchi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US09/959,575 external-priority patent/US6748333B1/en
Application filed by Individual filed Critical Individual
Priority to US12/155,812 priority Critical patent/US20080249717A1/en
Publication of US20080249717A1 publication Critical patent/US20080249717A1/en
Priority to US12/611,181 priority patent/US8078411B2/en
Priority to US12/983,714 priority patent/US8204696B2/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/86Signal analysis
    • G01N30/8624Detection of slopes or peaks; baseline correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/86Signal analysis
    • G01N30/8624Detection of slopes or peaks; baseline correction
    • G01N30/8631Peaks
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/86Signal analysis
    • G01N30/8624Detection of slopes or peaks; baseline correction
    • G01N30/8631Peaks
    • G01N30/8637Peak shape
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/86Signal analysis
    • G01N30/8624Detection of slopes or peaks; baseline correction
    • G01N30/8641Baseline

Definitions

  • the present invention relates to a chromatography technology such as a liquid chromatography technology, and particularly to a data processing method.
  • chromatographs such as a liquid chromatographic analyzer, a gas chromatographic analyzer or the like
  • a sample to be measured is let to pass through a column to be separated into components, and each of the separated components is detected as an output value at each elapsing time using a photometer such as a chromatographic detector.
  • Signals output from the chromatographic detector are recorded as time sequential data having a time interval of several hundreds ms. This is what is called a chromatogram having signal intensity in the ordinate and retention time in abscissa. In general, the signal intensity is converted to a digital value I j every an arbitrary time interval (time t j ) to execute data processing.
  • FIG. 3 shows an example of a chromatogram obtained by executing a body fluid amino acid analysis.
  • peaks of 11 components from Gly (glycine) to Tyr (tyrosine) densely exist in the range of retention time from 23 to 34 (min).
  • area-quantitative calculation is conventionally executed using a vertically dividing method in which a vertical line is drawn from each minimum point between peaks, that is, what is called “a root”.
  • this method produces an error as large as several tens % to cause an incorrect result when the peaks are strongly overlapped with each other. Therefore, when the chromatogram of such a kind needs to be quantitatively analyzed in a high accuracy, it has been general that the analyzing time is lengthened to improve the separation degree.
  • At least three independent parameters (A: area, T R : retention time, ⁇ : standard deviation) are used as variables in order to execute fitting of a peak for one component. Therefore, in order to execute fitting of peaks for a plurality of components, it is necessary to calculate three parameters of A i , T Ri , ⁇ i for each of the components (i).
  • overlapping peaks on a chromatogram are curve-fit using a waveform function such as the Gaussian function or an EMG function (exponentially modified Gaussian function) which can express an asymmetric waveform of a peak.
  • a waveform function such as the Gaussian function or an EMG function (exponentially modified Gaussian function) which can express an asymmetric waveform of a peak.
  • the overlapping peaks can be separated into individual peak waveforms, and the quantitative calculation can be performed by obtaining peak sizes such as a peak area and so on corresponding to a component of each of the peaks.
  • the separation of the peak waveforms is applied to two or three overlapping peaks, but not applied to the more number of overlapping peaks.
  • An object of the present invention is to provide a chromatographic analyzer capable of easily executing curve fitting using the non-linear least-square method to a chromatogram having a plurality of overlapping peaks.
  • the present invention to attain the above object is characterized by a chromatographic data processor for executing data processing of a chromatogram obtained by separating a sample to be measured using a column and detecting the separated sample, wherein fitting processing is executed to each peak in an arbitrary time region having the plurality of peaks of the chromatogram starting from the front side of the time region or from the back side of the time region, and the processed peaks are subtracted from the time region of the chromatogram so that the plurality of peaks in the chromatogram can be separated from one another.
  • FIG. 1 is a flowchart showing fitting processing in accordance with the present invention.
  • FIG. 2 is a diagram showing the outline of a chromatographic analyzer.
  • FIG. 3 is a chart showing a chromatogram for 70 minutes of a body fluid amino acid analysis method.
  • FIG. 4 is a chart showing an example of a display of a chromatogram in which a time window is set.
  • FIG. 5 is a view showing an example of a display of a dialogue box for setting a waveform function.
  • FIG. 6 is diagrams showing examples of weighting functions.
  • FIG. 7 is a view showing an example of a display of a dialogue box for selecting a weighting function and a direction of fitting.
  • FIG. 8 is a view showing a dialogue box for instructing execution of the fitting processing.
  • FIG. 9 is an illustrating chart showing a chromatogram in which a first fitting region is set.
  • FIG. 10 is an illustrating chart showing the chromatogram in which a first component is cut out.
  • FIG. 11 is a view showing an example of a display of a fitting result.
  • FIG. 12 is a chart showing an example of graphically setting the weighting function.
  • FIG. 2 is a diagram showing the outline of a liquid chromatographic analyzer to which the present invention is applied.
  • An eluting solution 1 is initially pumped to a column 3 using a fluid pump 2 by an instruction from a controller 5 .
  • a sample supply portion 8 is arranged between the fluid pump 2 and the column 3 , and a sample is supplied from a sampler 7 containing the sample to the eluting solution by an instruction of the controller 5 .
  • the sample is separated by the column 3 to be detected using a detector 4 such as a UV detector.
  • a chromatogram of the detected data is transmitted to the controller 5 to be data-analyzed, and the result is displayed on a display 6 or printed by a printer 9 .
  • the data processing of the chromatogram is executed mainly according to the following procedure.
  • Step 1 specification of a time interval to execute fitting thereto.
  • Step 2 selection of a weighting pattern.
  • Step 3 selection of a fitting direction.
  • Step 4 clicking of a fitting execution button.
  • Step 5 displaying and outputting of the result.
  • FIG. 1 is a flowchart showing the detailed process of the data processing of the chromatogram described above. In the present embodiment, description will be made on a case where the fitting processing is executed two components by two components.
  • a chromatogram to be executed fitting is selected from detected results obtained to start the processing ( 100 ).
  • an interval of retention time to execute fitting processing thereto (a time window) is set ( 101 ).
  • the peak group of 11 components from Gly to Tyr in the chromatogram of FIG. 4 is set to the time window.
  • the setting of time window is performed by dragging the time abscissa of the chromatogram displayed on the display 9 using a cursor. Otherwise, a peak starting and a peak ending the fitting may be selected by picking the peaks Gly and Tyr of the chromatogram displayed. Further, it is possible to employ a method of inputting a starting time and an ending time.
  • a waveform function used for the fitting calculation is set.
  • the waveform function for example, “Gaussian” or “EMG” may be selected. This selection can be performed by selecting a waveform function used using a dialogue box on the display 9 shown in FIG. 5 .
  • the weighting function is a function for weighting each of the two peaks.
  • various kinds of weighting functions may be employed. Some of the examples are described below.
  • FIG. 6 ( b ) a weighting function having a curvilinear gradient.
  • FIG. 6 ( c ) a weighting function having an S-shaped curvilinear gradient.
  • FIG. 6 ( d ) a weighting function having a linear gradient from a peak start point of a first peak to a peak end point of a second peak.
  • FIG. 6 ( e ) a weighting function having a weight of 0.5 at an end point of a second peak.
  • FIG. 6 ( f ) a weighting function having a linear gradient from a maximum point of a second peak to an end point of the second peak.
  • FIG. 6 ( g ) a weighting function having a linear gradient from a maximum point of a first peak to an end point of a second peak by particularly attaching importance to the left-hand side of the first peak.
  • the weight of 0.5 is set to an end point of the first peak (a start point of the second peak.
  • FIG. 6 ( h ) a weight function having a flat weight of 0.5 from an end point of a first peak (a start point of a second peak) to an end point of the second peak.
  • FIG. 6 ( i ) a weighting function which is equivalent to a case without any weight.
  • the weighting function is set by displaying a dialogue box shown in FIG. 7 on the display 9 and selecting a weighting pattern.
  • a weighting pattern Forin, by selecting an “ARBITRARY SET” and pushing an “OPTION” button, the above-described various kinds of pattern function type graphs shown in FIG. 6 can be set.
  • the various kinds of patterns of FIG. 6 are pre-stored in a memory such as a hard disk in the controller 5 , and can be easily read out from the memory to be used for the processing by being specified using a dialogue box of FIG. 7 .
  • a dialogue box shown in FIG. 8 is displayed after selecting an “OK” button in the dialogue box of FIG. 7 , and the fitting calculation processing is started by pushing an “EXECUTION” button.
  • inflection points inside the set time window are detected to obtain maximum points and minimum points, and an interval for each of the peaks is defined by determining the maximum point as the apex and the minimum points as the end point and the starting point of the peak.
  • the roots (the minimum points) of the peak and a retention time having the points are stored in relation to each other.
  • a width of the weighting function set in Process 103 is adjusted so as to meet with a width between the roots detected in Process 105 ( 108 ).
  • FIG. 9 is an illustrating chart showing a chromatogram in which a region of first fitting processing is set. The width of the weighting function is adjusted so as to meet with the width of the fitting region as shown in FIG. 9 .
  • the fitting processing is executed to the two components of the front side in the fitting direction ( 109 ).
  • the Gaussian or EMG function is used as the waveform function f(t) for the fitting.
  • the waveform function is shown in Equation 1.
  • variable i is incremented, and the processing is returned to Process 107 .
  • the fitting processing for the peaks Ala and Cit is executed to determine the second component waveform of the second peak Ala according to the same procedure as that of the fitting processing for the first peak.
  • the processing proceeds to Process 113 to determine the peak waveforms of the final two components by executing the peak fitting for the two components without weighting ( 114 ). Then, the peak waveforms for the two components are subtracted from the original chromatogram ( 115 ) to complete the processing.
  • the uncertainties (deviations) of the parameters for each of the components can be obtained by the above-described processing, the uncertainty of each of the measured quantitative values can be calculated based on the above-described uncertainties.
  • the point that the uncertainties of the measured quantitative values can be calculated is also one of the advantages of the present embodiment.
  • the calculation of obtaining the uncertainties of the measured quantitative values is executed according to the error propagation equation.
  • the quality of the fitting can be judged from the magnitude of the uncertainty of the measured quantitative value (it can be said that the quality is better when the uncertainty is smaller).
  • the uncertainty can be used for judging quality of the calculation result of the overlapping portion when the fitting processing executed from both of the front side end and the back side end which is to be described later.
  • the above-described processing executing separating the peaks while the peaks are successively being cutting out is called as “sequencing”, and the chromatographic data processor for executing the sequencing is called as “a chromatographic peak sequencer”.
  • the fitting calculation described above is executed from the front side end
  • the fitting calculation may be executed from the back side end by setting the fitting direction of the dialogue box of FIG. 7 .
  • the function g 2 (t) expressing the second peak is subtracted from the chromatogram, and the shape of the weighting function is reversed left to right.
  • the peaks are cut out from the peaks having the larger retention time, that is, in order of the peaks Tyr, Leu, Ile (isoleucine) in the case of FIG. 4 .
  • the fitting processing may be executed from either of the front side and the back side. However, actually, it is preferable that the fitting processing is executed from the both sides toward a peak having the strongest overlapping intensity.
  • the button “AUTOMATIC JUDGMENT” in the dialogue box of FIG. 7 is selected.
  • the forward sequencing processing is executed from the peak Gly to the peak Cysthi (Cystathionine) or Ile and the backward sequencing processing is executed to the peak Lue or Ile because overlapping of the peaks Lue and Ile is strong.
  • the peaks Lue and Ile which fitting result should be employed and determined can be judged by comparing the appropriateness of the fitting results between the forward processing and the backward processing.
  • the judgment can be executed by using the uncertainties (the deviations) of the fitting parameters A i , t Ri and ⁇ i themselves as the indexes for judging the appropriateness or by using the statistical index ⁇ 2 or number of repetitive calculating times for watching the convergence rate of ⁇ 2 . Further, more comprehensive judgment can be performed by combining these parameters. Any way, the peak area value A i needs to be accurately calculated in to execute quantitative value calculation, and accordingly importance should be placed on the uncertainty of the peak area value A i .
  • the weighting function set in the above-described process 102 has been set on the dialogue box, the weighting function may be set through a program setting method in which a table method shown in Table 1 is used or through a graphical setting method in which setting is performed by displaying and superposing the weighting function on a chromatogram as shown in FIG. 12 .
  • Weighting Characteristic point value Shape 0 starting point of first 1 drawing peak a vertical line horizontal line apex of first peak 1 gradient (straight line) root between peaks 0.5 gradient (convex) end point of second peak 0
  • the setting program of Table 1 shows the example of setting the weighting function of FIG. 6 ( g ).
  • the weighting function is defined by that a measurer inputs “weighting values” and “shapes”. Table 1 is displayed on the display.
  • FIG. 12 shows the method of setting the weighting function by superposing on the chromatogram on the display. For example, a cursor is moved to a point near a node 31 using a pointing device, and a characteristic point of the first peak starting point is displayed by clicking on the point. By doing so, the node 31 is determined. When the node 31 needs to be moved, the node is picked and then released. That is, the node may be dragged. Next, a node 32 is determined by clicking on a point near the node 32 and vertically dragging the node 31 to a point having a weight of 1. Therein, a vertical line is automatically drawn from the node 31 to the node 32 and determined. Similarly, nodes 33 and 34 are determined to set the weighting function.
  • changing of setting can be performed by clicking an arbitrary node to highlight the node. Further, a node can be deleted. Furthermore, a new node can be added by clicking an arbitrary point.
  • a line section between nodes is changed to a curvilinear line, the line section to be changed is clicked to be highlighted (the displaying color is changed), and the set attribution is changed by light-hand side clicking.
  • concave, convex or the like can be specified.
  • the measurement is once performed to obtain a chromatogram, and after that the data processing is performed.
  • the simultaneous sequencing processing with measurement is progressed as follows.
  • Step 1 setting of a time interval (a time window) to execute fitting thereto.
  • Step 2 setting of a weighting pattern.
  • Step 3 setting of a fitting direction.
  • Step 4 measuring of a sample.
  • Step 5 executing of fitting processing.
  • Step 6 displaying and outputting of the result.
  • a time program is set before starting measurement as shown in Table 2.
  • a time window is set by specifying a time period from the fitting start time to the fitting end time.
  • the time window may be also set by selecting a fitting start peak and a fitting end peak as by specifying peak names “Gly” to “Tyr”.
  • fitting of a chromatogram having a plurality of overlapping peaks can be automatically executed only by setting some number of conditions.
  • a highly accurate base line can be calculate at the same time.
  • the chromatographic peak sequencer in accordance with the present invention is considered to be an excellent data processor.
  • the fitting is executed by assuming the shape of peak waveform as a shape of Gaussian or EMG function, an attention should be paid on this point when a shape of an actual peak waveform is significantly different from the shapes of these analytical functions. In such a case, it is necessary to pre-store specific waveform functions by obtaining an isolated peak waveform for each component from a standard sample.
  • the peaks in a chromatogram having a plurality of overlapping peaks, particularly, having three or more overlapping peaks, the peaks can be easily separated only by setting some number of setting conditions. By doing so, the accuracy of quantitative analysis and qualitative analysis can be improved. Further, the capability of determining the base line can be also improved. Furthermore, measurer's labor for the data processing can be substantially reduced.

Abstract

A chromatographic analyzer is provided for facilitating curve fitting by means of the linear least-square method for a chromatogram that contains a plurality of overlapping peaks. The present invention is characterized by a chromatographic data processor for executing data processing of a chromatogram obtained by separating a sample to be measured using a column and detecting the separated sample, wherein fitting processing is executed to each peak in an arbitrary time region having the plurality of peaks of the chromatogram starting from the front side of the time region or from the back side of the time region, and the processed peaks are subtracted from the chromatogram in the time region so that the plurality of peaks in the chromatogram can be separated from one another. Thus, the plurality of overlapping peaks, particularly three or more overlapping peaks in the chromatogram can be easily separated from one another only by defining some setting conditions.

Description

    TECHNICAL FIELD
  • The present invention relates to a chromatography technology such as a liquid chromatography technology, and particularly to a data processing method.
  • BACKGROUND ART
  • In chromatographs such as a liquid chromatographic analyzer, a gas chromatographic analyzer or the like, a sample to be measured is let to pass through a column to be separated into components, and each of the separated components is detected as an output value at each elapsing time using a photometer such as a chromatographic detector.
  • Signals output from the chromatographic detector are recorded as time sequential data having a time interval of several hundreds ms. This is what is called a chromatogram having signal intensity in the ordinate and retention time in abscissa. In general, the signal intensity is converted to a digital value Ij every an arbitrary time interval (time tj) to execute data processing.
  • FIG. 3 shows an example of a chromatogram obtained by executing a body fluid amino acid analysis.
  • As shown in FIG. 3, peaks of 11 components from Gly (glycine) to Tyr (tyrosine) densely exist in the range of retention time from 23 to 34 (min). In such a case, area-quantitative calculation is conventionally executed using a vertically dividing method in which a vertical line is drawn from each minimum point between peaks, that is, what is called “a root”. However, this method produces an error as large as several tens % to cause an incorrect result when the peaks are strongly overlapped with each other. Therefore, when the chromatogram of such a kind needs to be quantitatively analyzed in a high accuracy, it has been general that the analyzing time is lengthened to improve the separation degree.
  • On the other hand, in order to perform quantitative calculation without lengthening the analyzing time even if peaks are overlapped so strongly with each other, quantitative calculation is tried to be performed using numerical analysis in a manner like data processing. This method is, for example, a non-linear least-square method.
  • In the case of using the non-linear least-square method, at least three independent parameters (A: area, TR: retention time, σ: standard deviation) are used as variables in order to execute fitting of a peak for one component. Therefore, in order to execute fitting of peaks for a plurality of components, it is necessary to calculate three parameters of Ai, TRi, σi for each of the components (i).
  • The conventional examples of using the non-linear least-square method are disclosed in Japanese Patent Application Laid-Open No.6-324029 and Japanese Patent Application Laid-Open No.63-151851.
  • These examples disclose that overlapping peaks on a chromatogram are curve-fit using a waveform function such as the Gaussian function or an EMG function (exponentially modified Gaussian function) which can express an asymmetric waveform of a peak. As shown in these examples, the overlapping peaks can be separated into individual peak waveforms, and the quantitative calculation can be performed by obtaining peak sizes such as a peak area and so on corresponding to a component of each of the peaks.
  • DISCLOSURE OF INVENTION
  • However, in the conventional examples using the non-linear least-square method, the separation of the peak waveforms is applied to two or three overlapping peaks, but not applied to the more number of overlapping peaks.
  • The reason is that in the case of separating the overlapping peak waveforms using curve fitting through the non-linear least-squire method, as the number of peak components is increased to 3, 4, 5, . . . , there occurs a phenomenon that the fitting processing is difficult to be converged or that the separation of peaks can not correctly performed (the error is increased).
  • For example, in the case of the chromatogram shown in FIG. 3, when fitting is tried to the 11 components from Gly to Tyr at a time, the 33 parameters of 11×3 must be determined at a time. This is very difficult calculation processing to the non-linear least-square method, and accordingly various kinds of techniques are necessary in order to solve this problem. Therefore, when the curve-fitting is executed using the non-linear least-square method in the case of existing many overlapping peaks on a chromatogram, a measuring operator must specify calculation regions (time windows) for 2 or 3 peaks seeming to be converged one by one. This process expenses much time and much effort, and in addition, there is a problem in the reliability of the calculation result because the regions are artificially determined.
  • An object of the present invention is to provide a chromatographic analyzer capable of easily executing curve fitting using the non-linear least-square method to a chromatogram having a plurality of overlapping peaks.
  • The present invention to attain the above object is characterized by a chromatographic data processor for executing data processing of a chromatogram obtained by separating a sample to be measured using a column and detecting the separated sample, wherein fitting processing is executed to each peak in an arbitrary time region having the plurality of peaks of the chromatogram starting from the front side of the time region or from the back side of the time region, and the processed peaks are subtracted from the time region of the chromatogram so that the plurality of peaks in the chromatogram can be separated from one another.
  • The object, the operation and the effect of the present invention will be described in detail in the section of DESCRIPTION OF THE PREFERRED EMBODIMENTS to be described later.
  • BRIEF DESCRIPTION OF DRAWINGS
  • FIG. 1 is a flowchart showing fitting processing in accordance with the present invention.
  • FIG. 2 is a diagram showing the outline of a chromatographic analyzer.
  • FIG. 3 is a chart showing a chromatogram for 70 minutes of a body fluid amino acid analysis method.
  • FIG. 4 is a chart showing an example of a display of a chromatogram in which a time window is set.
  • FIG. 5 is a view showing an example of a display of a dialogue box for setting a waveform function.
  • FIG. 6 is diagrams showing examples of weighting functions.
  • FIG. 7 is a view showing an example of a display of a dialogue box for selecting a weighting function and a direction of fitting.
  • FIG. 8 is a view showing a dialogue box for instructing execution of the fitting processing.
  • FIG. 9 is an illustrating chart showing a chromatogram in which a first fitting region is set.
  • FIG. 10 is an illustrating chart showing the chromatogram in which a first component is cut out.
  • FIG. 11 is a view showing an example of a display of a fitting result.
  • FIG. 12 is a chart showing an example of graphically setting the weighting function.
  • BEST MODE IN WHICH THE INVENTION IS CARRIED OUT
  • An embodiment of the present invention will be described below.
  • FIG. 2 is a diagram showing the outline of a liquid chromatographic analyzer to which the present invention is applied. An eluting solution 1 is initially pumped to a column 3 using a fluid pump 2 by an instruction from a controller 5. A sample supply portion 8 is arranged between the fluid pump 2 and the column 3, and a sample is supplied from a sampler 7 containing the sample to the eluting solution by an instruction of the controller 5. The sample is separated by the column 3 to be detected using a detector 4 such as a UV detector. A chromatogram of the detected data is transmitted to the controller 5 to be data-analyzed, and the result is displayed on a display 6 or printed by a printer 9.
  • Data processing in the controller 5 will be described below when the chromatogram of FIG. 3 is obtained as the detected result.
  • In the present embodiment, the data processing of the chromatogram is executed mainly according to the following procedure.
  • Step 1: specification of a time interval to execute fitting thereto.
  • Step 2: selection of a weighting pattern.
  • Step 3: selection of a fitting direction.
  • Step 4: clicking of a fitting execution button.
  • Step 5: displaying and outputting of the result.
  • FIG. 1 is a flowchart showing the detailed process of the data processing of the chromatogram described above. In the present embodiment, description will be made on a case where the fitting processing is executed two components by two components.
  • Initially, a chromatogram to be executed fitting is selected from detected results obtained to start the processing (100). Here, the variable i is set to 1 (one), that is i=1.
  • Next, an interval of retention time to execute fitting processing thereto (a time window) is set (101). Here, the peak group of 11 components from Gly to Tyr in the chromatogram of FIG. 4 is set to the time window.
  • The setting of time window is performed by dragging the time abscissa of the chromatogram displayed on the display 9 using a cursor. Otherwise, a peak starting and a peak ending the fitting may be selected by picking the peaks Gly and Tyr of the chromatogram displayed. Further, it is possible to employ a method of inputting a starting time and an ending time.
  • Next, a waveform function used for the fitting calculation is set. As the waveform function, for example, “Gaussian” or “EMG” may be selected. This selection can be performed by selecting a waveform function used using a dialogue box on the display 9 shown in FIG. 5.
  • Although the fitting in the present embodiment is executed two components by two components, it is necessary to make the influence of signals of a second and later waveforms in order to accurately determine a waveform of a peak in the front side of the fitting direction, that is, a first peak. Therefore, a weighting function shown in FIG. 6 (a) is set (103).
  • The weighting function is a function for weighting each of the two peaks. The example of FIG. 6 (a) is a function that the weight w is set to 1 (one), that is, w=1 from a peak start time 21 to a peak end (root) time 23 of the first peak, and then the weight w is linearly slanted to be reduced from 1 to 0 from a peak start time 23 to a peak end (root) time 25 of the second peak. As the weighting function, various kinds of weighting functions may be employed. Some of the examples are described below.
  • FIG. 6 (b): a weighting function having a curvilinear gradient. FIG. 6 (c): a weighting function having an S-shaped curvilinear gradient. FIG. 6 (d): a weighting function having a linear gradient from a peak start point of a first peak to a peak end point of a second peak. FIG. 6 (e): a weighting function having a weight of 0.5 at an end point of a second peak. FIG. 6 (f): a weighting function having a linear gradient from a maximum point of a second peak to an end point of the second peak. FIG. 6 (g): a weighting function having a linear gradient from a maximum point of a first peak to an end point of a second peak by particularly attaching importance to the left-hand side of the first peak. The weight of 0.5 is set to an end point of the first peak (a start point of the second peak. FIG. 6 (h): a weight function having a flat weight of 0.5 from an end point of a first peak (a start point of a second peak) to an end point of the second peak. FIG. 6 (i): a weighting function which is equivalent to a case without any weight.
  • In the concrete, the weighting function is set by displaying a dialogue box shown in FIG. 7 on the display 9 and selecting a weighting pattern. Therein, by selecting an “ARBITRARY SET” and pushing an “OPTION” button, the above-described various kinds of pattern function type graphs shown in FIG. 6 can be set. The various kinds of patterns of FIG. 6 are pre-stored in a memory such as a hard disk in the controller 5, and can be easily read out from the memory to be used for the processing by being specified using a dialogue box of FIG. 7.
  • Next, a direction of executing fitting calculation is determined (104).
  • That is, describing the example of FIG. 3, it is decided whether the processing is executed from the front side (the Gly side) or the back side (the Tyr side). In the concrete, the setting of the direction executing the fitting calculation is set using the dialogue box shown in FIG. 7.
  • An example of executing the processing from the Gly side will be described below.
  • By the above, setting of the conditions necessary for the fitting processing has been completed.
  • Then, a dialogue box shown in FIG. 8 is displayed after selecting an “OK” button in the dialogue box of FIG. 7, and the fitting calculation processing is started by pushing an “EXECUTION” button.
  • As the calculation processing is started, number of peaks within the set time window is detected to be set to Imax (106).
  • In the concrete, inflection points inside the set time window are detected to obtain maximum points and minimum points, and an interval for each of the peaks is defined by determining the maximum point as the apex and the minimum points as the end point and the starting point of the peak. Therein, the roots (the minimum points) of the peak and a retention time having the points are stored in relation to each other.
  • Next, a difference between Imax and i is calculated (107). Therein, if the difference is larger than 3 (three components), the processing proceeds to Process 108.
  • When it is judged in Process 107 that the difference is larger than 3 (three components), a width of the weighting function set in Process 103 is adjusted so as to meet with a width between the roots detected in Process 105 (108).
  • FIG. 9 is an illustrating chart showing a chromatogram in which a region of first fitting processing is set. The width of the weighting function is adjusted so as to meet with the width of the fitting region as shown in FIG. 9.
  • Next, the fitting processing is executed to the two components of the front side in the fitting direction (109).
  • The fitting processing will be described below in detail.
  • As described above, the Gaussian or EMG function is used as the waveform function f(t) for the fitting. Firstly, in a case of employing the Gaussian function, the waveform function is shown in Equation 1.
  • f ( t ) = g 1 ( t ) + g 2 ( t ) + at + b = A 1 2 π σ 1 exp { - ( t - t R 1 ) 2 2 σ 1 2 } + A 2 2 π σ 2 exp { - ( t - t R 2 ) 2 2 σ 2 2 } + at + b ( g i ( t ) : the suffix i denotes a component i . ) ( Equation 1 )
  • There, the term (at+b) denotes a base line. Though the waveform function becomes Equation 1 shown above because the fitting processing is executed two components by two components in the present embodiment, the fitting processing may be executed three components by three components if the capability of non-linear least-square method can be secured. Of course, extension to the EMG function may be possible. In a case of employing the EMG function, the term gi(t) of the waveform function in Equation 1 is replaced by Equation 2.
  • g 1 EMG ( t ) = A i 2 π σ i τ i 0 t exp { - ( t - t Ri - t ) 2 2 σ i 2 - t τ i } t ( Equation 2 )
  • In the least-square method, each of the fitting parameters Ai, tRI, σi and τi is determined so as to minimize the following S1 of Equation 3 or the following S2 of Equation 4. There, Ij is a signal intensity of a measured chromatogram, j is a suffix expressing time tj, and N is number of data points in a time interval. Either of S1 and S2 may be used.
  • S 1 = j = 1 N w j ( f ( t j ) - I j ) 2 ( Equation 3 ) S 2 = j = 1 N w j 2 ( f ( t j ) - I j ) 2 ( Equation 4 )
  • The function gi(t) having a smaller value of tR (that is, the earlier retention time) in the f(t) of Equation 1 obtained here corresponds to the first peak Gly, and accordingly the waveform of the first component is determined (110).
  • Next, the waveform of the first component gi(t) is subtracted from the original measured chromatogram to form a chromatogram cut out the first peak Gly (111). FIG. 10 is an illustrating chart showing the chromatogram in which the first peak Gly has been cut out. It is preferable not to subtract the base line at+b from the original chromatogram because it does not influence the peak waveform.
  • Next, the variable i is incremented, and the processing is returned to Process 107.
  • Then, the fitting processing for the peaks Ala and Cit is executed to determine the second component waveform of the second peak Ala according to the same procedure as that of the fitting processing for the first peak.
  • After that, the peaks Ala, Cit . . . are successively subtracted until the peaks for two components Leu (leucine) and Tyr in the back side end. When the remaining becomes two components (i=2), the processing proceeds to Process 113 to determine the peak waveforms of the final two components by executing the peak fitting for the two components without weighting (114). Then, the peak waveforms for the two components are subtracted from the original chromatogram (115) to complete the processing.
  • By the series of the processing described above, all the waveforms in the time window have been separated from one another. Further, according to the above-described processing, since the remaining chromatogram after cutting out all the peak waveforms becomes a base line, the base line can be obtained at a time.
  • After completion of the above-described processing, quantitative calculation (for example, concentration) is executed using parameters for each of the peaks obtained by the above-described processing. Then, a display for each of the peaks (each of the components) shown in FIG. 11 is displayed on the display 9 as the calculation result.
  • There, because the uncertainties (deviations) of the parameters for each of the components can be obtained by the above-described processing, the uncertainty of each of the measured quantitative values can be calculated based on the above-described uncertainties.
  • The point that the uncertainties of the measured quantitative values can be calculated is also one of the advantages of the present embodiment. the calculation of obtaining the uncertainties of the measured quantitative values is executed according to the error propagation equation.
  • Further, the quality of the fitting can be judged from the magnitude of the uncertainty of the measured quantitative value (it can be said that the quality is better when the uncertainty is smaller). For example, the uncertainty can be used for judging quality of the calculation result of the overlapping portion when the fitting processing executed from both of the front side end and the back side end which is to be described later.
  • In the present invention, the above-described processing executing separating the peaks while the peaks are successively being cutting out is called as “sequencing”, and the chromatographic data processor for executing the sequencing is called as “a chromatographic peak sequencer”.
  • Although the fitting calculation described above is executed from the front side end, the fitting calculation may be executed from the back side end by setting the fitting direction of the dialogue box of FIG. 7. In this case, the function g2(t) expressing the second peak is subtracted from the chromatogram, and the shape of the weighting function is reversed left to right. In the processing direction described above, the peaks are cut out from the peaks having the larger retention time, that is, in order of the peaks Tyr, Leu, Ile (isoleucine) in the case of FIG. 4.
  • Ideally, the fitting processing may be executed from either of the front side and the back side. However, actually, it is preferable that the fitting processing is executed from the both sides toward a peak having the strongest overlapping intensity. When such fitting processing is performed, the button “AUTOMATIC JUDGMENT” in the dialogue box of FIG. 7 is selected.
  • For example, in the case of the chromatogram of FIG. 3, the forward sequencing processing is executed from the peak Gly to the peak Cysthi (Cystathionine) or Ile and the backward sequencing processing is executed to the peak Lue or Ile because overlapping of the peaks Lue and Ile is strong. By doing so, in regard to the peaks Lue and Ile, which fitting result should be employed and determined can be judged by comparing the appropriateness of the fitting results between the forward processing and the backward processing.
  • Therein, the judgment can be executed by using the uncertainties (the deviations) of the fitting parameters Ai, tRi and σi themselves as the indexes for judging the appropriateness or by using the statistical index χ2 or number of repetitive calculating times for watching the convergence rate of χ2. Further, more comprehensive judgment can be performed by combining these parameters. Any way, the peak area value Ai needs to be accurately calculated in to execute quantitative value calculation, and accordingly importance should be placed on the uncertainty of the peak area value Ai.
  • Another Embodiment of Setting the Weighting Function
  • Although the weighting function set in the above-described process 102 has been set on the dialogue box, the weighting function may be set through a program setting method in which a table method shown in Table 1 is used or through a graphical setting method in which setting is performed by displaying and superposing the weighting function on a chromatogram as shown in FIG. 12.
  • TABLE 1
    Weighting program (expressing FIG. 6 (g))
    Weighting
    Characteristic point value Shape
    0
    starting point of first 1 drawing
    peak a vertical line
    horizontal line
    apex of first peak 1 gradient
    (straight line)
    root between peaks 0.5 gradient (convex)
    end point of second peak 0
  • The setting program of Table 1 shows the example of setting the weighting function of FIG. 6 (g). The weighting function is defined by that a measurer inputs “weighting values” and “shapes”. Table 1 is displayed on the display.
  • FIG. 12 shows the method of setting the weighting function by superposing on the chromatogram on the display. For example, a cursor is moved to a point near a node 31 using a pointing device, and a characteristic point of the first peak starting point is displayed by clicking on the point. By doing so, the node 31 is determined. When the node 31 needs to be moved, the node is picked and then released. That is, the node may be dragged. Next, a node 32 is determined by clicking on a point near the node 32 and vertically dragging the node 31 to a point having a weight of 1. Therein, a vertical line is automatically drawn from the node 31 to the node 32 and determined. Similarly, nodes 33 and 34 are determined to set the weighting function. When the weighting function needs to be changed, changing of setting can be performed by clicking an arbitrary node to highlight the node. Further, a node can be deleted. Furthermore, a new node can be added by clicking an arbitrary point. When a line section between nodes is changed to a curvilinear line, the line section to be changed is clicked to be highlighted (the displaying color is changed), and the set attribution is changed by light-hand side clicking. Therein, concave, convex or the like can be specified.
  • Another Embodiment of Sequencing
  • In the embodiment of sequencing described above, the measurement is once performed to obtain a chromatogram, and after that the data processing is performed. However, it is possible to simultaneously execute the sequencing processing with measuring of the sample. An embodiment in such a case will be described below.
  • The simultaneous sequencing processing with measurement is progressed as follows.
  • Step 1: setting of a time interval (a time window) to execute fitting thereto.
  • Step 2: setting of a weighting pattern.
  • Step 3: setting of a fitting direction.
  • Step 4: measuring of a sample.
  • Step 5: executing of fitting processing.
  • Step 6: displaying and outputting of the result.
  • In steps 1 to 3, a time program is set before starting measurement as shown in Table 2. A time window is set by specifying a time period from the fitting start time to the fitting end time. In addition, the time window may be also set by selecting a fitting start peak and a fitting end peak as by specifying peak names “Gly” to “Tyr”.
  • TABLE 2
    Time Program
    Time Command Status
    0.0 Waveform function Gaussian
    0.0 Weighting function Straight line gradient
    0.0 Fitting direction Automatic judgment
    23.0 Fitting start ON
    34.0 Fitting end ON
  • According to the chromatographic peak sequencer in accordance with the present invention, fitting of a chromatogram having a plurality of overlapping peaks can be automatically executed only by setting some number of conditions.
  • Further, a highly accurate base line can be calculate at the same time. From the viewpoint of existing the highly accurate base line determining function, the chromatographic peak sequencer in accordance with the present invention is considered to be an excellent data processor. However, because the fitting is executed by assuming the shape of peak waveform as a shape of Gaussian or EMG function, an attention should be paid on this point when a shape of an actual peak waveform is significantly different from the shapes of these analytical functions. In such a case, it is necessary to pre-store specific waveform functions by obtaining an isolated peak waveform for each component from a standard sample.
  • According to the present invention, in a chromatogram having a plurality of overlapping peaks, particularly, having three or more overlapping peaks, the peaks can be easily separated only by setting some number of setting conditions. By doing so, the accuracy of quantitative analysis and qualitative analysis can be improved. Further, the capability of determining the base line can be also improved. Furthermore, measurer's labor for the data processing can be substantially reduced.

Claims (2)

1-12. (canceled)
13. A liquid chromatographic analyzer comprising:
a sample separating unit which separates a sample to be measured into components, the sample separating unit comprising: a fluid pump which feeds an eluting solution at a constant rate; an eluting solution storing unit which stores the eluting solution; a sample storing unit which stores the sample; and a separation column to which the sample and the eluting solution are fed to perform an elution which separates the sample into the components;
a detecting unit which detects the separated sample, the detecting unit comprising a photometer, the photometer being capable of detecting each of the components as an output value at each elapsing time;
a data processing unit which executes data processing by obtaining detected result from the detecting unit, the data processing unit comprising: a unit which specifies an arbitrary retention time region having a plurality of peaks of a chromatogram obtained as the detected result of the detecting unit; a unit which specifies whether fitting processing is executed from a front side end or a back side end of the specified retention time region; and a unit which specifies a weighting function and a waveform function used in the fitting processing; and
a display unit which displays results obtained from the data processing unit,
wherein a new retention time region, which has peaks other than those peaks at which the fitting processing is completed to define their peak shapes, is displayed on the display unit, and the specifying of the weighting function is performed with conceptual graphics of candidate weighting functions superposed over the chromatogram displayed on the display unit.
US12/155,812 1999-09-27 2008-06-10 Method and apparatus for chromatographic data processing Abandoned US20080249717A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US12/155,812 US20080249717A1 (en) 2001-10-30 2008-06-10 Method and apparatus for chromatographic data processing
US12/611,181 US8078411B2 (en) 1999-09-27 2009-11-03 Method and apparatus for chromatographic data processing
US12/983,714 US8204696B2 (en) 1999-09-27 2011-01-03 Method and apparatus for chromatographic data processing

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US09/959,575 US6748333B1 (en) 1999-09-27 1999-09-27 Method and apparatus for chromatographic data processing, and chromatograph
US10/831,189 US6907355B2 (en) 1999-09-27 2004-04-26 Method and apparatus for chromatographic data processing
US11/061,610 US7200494B2 (en) 2001-10-30 2005-02-22 Method and apparatus for chromatographic data processing
US11/654,577 US7403859B2 (en) 1999-09-27 2007-01-18 Method and apparatus for chromatographic data processing
US12/155,812 US20080249717A1 (en) 2001-10-30 2008-06-10 Method and apparatus for chromatographic data processing

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US11/654,577 Continuation US7403859B2 (en) 1999-09-27 2007-01-18 Method and apparatus for chromatographic data processing

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US12/611,181 Continuation US8078411B2 (en) 1999-09-27 2009-11-03 Method and apparatus for chromatographic data processing

Publications (1)

Publication Number Publication Date
US20080249717A1 true US20080249717A1 (en) 2008-10-09

Family

ID=34704459

Family Applications (5)

Application Number Title Priority Date Filing Date
US11/061,610 Expired - Fee Related US7200494B2 (en) 1999-09-27 2005-02-22 Method and apparatus for chromatographic data processing
US11/654,577 Expired - Fee Related US7403859B2 (en) 1999-09-27 2007-01-18 Method and apparatus for chromatographic data processing
US12/155,812 Abandoned US20080249717A1 (en) 1999-09-27 2008-06-10 Method and apparatus for chromatographic data processing
US12/611,181 Expired - Fee Related US8078411B2 (en) 1999-09-27 2009-11-03 Method and apparatus for chromatographic data processing
US12/983,714 Expired - Fee Related US8204696B2 (en) 1999-09-27 2011-01-03 Method and apparatus for chromatographic data processing

Family Applications Before (2)

Application Number Title Priority Date Filing Date
US11/061,610 Expired - Fee Related US7200494B2 (en) 1999-09-27 2005-02-22 Method and apparatus for chromatographic data processing
US11/654,577 Expired - Fee Related US7403859B2 (en) 1999-09-27 2007-01-18 Method and apparatus for chromatographic data processing

Family Applications After (2)

Application Number Title Priority Date Filing Date
US12/611,181 Expired - Fee Related US8078411B2 (en) 1999-09-27 2009-11-03 Method and apparatus for chromatographic data processing
US12/983,714 Expired - Fee Related US8204696B2 (en) 1999-09-27 2011-01-03 Method and apparatus for chromatographic data processing

Country Status (1)

Country Link
US (5) US7200494B2 (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7200494B2 (en) * 2001-10-30 2007-04-03 Hitachi, Ltd. Method and apparatus for chromatographic data processing
US20060255258A1 (en) * 2005-04-11 2006-11-16 Yongdong Wang Chromatographic and mass spectral date analysis
FR2898283B1 (en) * 2006-03-08 2011-07-15 Novasep METHOD AND DEVICE FOR SEPARATING FRACTIONS FROM A MIXTURE
DE102009035587B3 (en) * 2009-07-31 2011-03-24 Siemens Aktiengesellschaft Method for filtering a chromatogram
WO2011021198A2 (en) 2009-08-20 2011-02-24 Spectrosense Ltd. Gas chromatographic analysis method and system
JP5947567B2 (en) * 2012-03-02 2016-07-06 株式会社日立ハイテクノロジーズ Mass spectrometry system
CN102735783B (en) * 2012-06-25 2015-02-18 上海申瑞继保电气有限公司 Method for calculating chromatogram area of gas in oil
WO2015052721A1 (en) * 2013-10-09 2015-04-16 Spectrosense Ltd. Modified data representation in gas chromatographic analysis
US11244818B2 (en) * 2018-02-19 2022-02-08 Agilent Technologies, Inc. Method for finding species peaks in mass spectrometry
CN113742652A (en) * 2021-09-06 2021-12-03 杭州凯曼健康科技有限公司 Immunofluorescence chromatography curve analysis method and device and electronic equipment

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4807148A (en) * 1987-05-29 1989-02-21 Hewlett-Packard Company Deconvolving chromatographic peaks
US4941101A (en) * 1987-09-01 1990-07-10 Hewlett-Packard Company Method for analyzing chromatograms
US5277871A (en) * 1989-10-20 1994-01-11 Hitachi, Ltd. Liquid chromatographic analyzer, sample feeder and prelabeling reaction treating method
US5291426A (en) * 1991-02-27 1994-03-01 The Perkin-Elmer Corporation Method of correcting spectral data for background
US6029115A (en) * 1996-10-03 2000-02-22 Perkin Elmer Llc Analyzing spectrometric data
US6393368B1 (en) * 1994-03-28 2002-05-21 Hitachi, Ltd. Method and apparatus for analyzing multi-channel chromatogram
US6529836B2 (en) * 1995-02-27 2003-03-04 Hitachi, Ltd. Data processing apparatus for chromatograph
US6549861B1 (en) * 2000-08-10 2003-04-15 Euro-Celtique, S.A. Automated system and method for spectroscopic analysis
US6815668B2 (en) * 1999-07-21 2004-11-09 The Charles Stark Draper Laboratory, Inc. Method and apparatus for chromatography-high field asymmetric waveform ion mobility spectrometry
US7200494B2 (en) * 2001-10-30 2007-04-03 Hitachi, Ltd. Method and apparatus for chromatographic data processing

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS4745556B1 (en) 1970-09-07 1972-11-16
JPS5321313B2 (en) 1973-02-26 1978-07-01
JPS6022659A (en) 1983-07-18 1985-02-05 Asahi Chem Ind Co Ltd Calculator for peak value of chromatogram
JPS6285850A (en) 1985-10-11 1987-04-20 Shimadzu Corp Separation of component peak of energy spectra
JPS63151851A (en) 1986-12-17 1988-06-24 Hitachi Ltd Data processor for chromatograph
JP3071085B2 (en) 1993-03-15 2000-07-31 株式会社日立製作所 Chromatogram analysis display method and apparatus
SE503289C2 (en) 1993-09-28 1996-05-13 Anders Hanning Ways to determine the concentration of substances in a liquid
JPH0954071A (en) 1995-08-17 1997-02-25 Kao Corp Gas chromatograph
WO2001023880A1 (en) * 1999-09-27 2001-04-05 Hitachi, Ltd. Method and apparatus for chromatographic data processing, and chromatograph
JP4745556B2 (en) 2001-08-20 2011-08-10 キヤノン株式会社 Positioning apparatus, exposure apparatus, and device manufacturing method

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4807148A (en) * 1987-05-29 1989-02-21 Hewlett-Packard Company Deconvolving chromatographic peaks
US4941101A (en) * 1987-09-01 1990-07-10 Hewlett-Packard Company Method for analyzing chromatograms
US5277871A (en) * 1989-10-20 1994-01-11 Hitachi, Ltd. Liquid chromatographic analyzer, sample feeder and prelabeling reaction treating method
US5291426A (en) * 1991-02-27 1994-03-01 The Perkin-Elmer Corporation Method of correcting spectral data for background
US6393368B1 (en) * 1994-03-28 2002-05-21 Hitachi, Ltd. Method and apparatus for analyzing multi-channel chromatogram
US6529836B2 (en) * 1995-02-27 2003-03-04 Hitachi, Ltd. Data processing apparatus for chromatograph
US6029115A (en) * 1996-10-03 2000-02-22 Perkin Elmer Llc Analyzing spectrometric data
US6815668B2 (en) * 1999-07-21 2004-11-09 The Charles Stark Draper Laboratory, Inc. Method and apparatus for chromatography-high field asymmetric waveform ion mobility spectrometry
US7403859B2 (en) * 1999-09-27 2008-07-22 Hitachi, Ltd. Method and apparatus for chromatographic data processing
US6549861B1 (en) * 2000-08-10 2003-04-15 Euro-Celtique, S.A. Automated system and method for spectroscopic analysis
US7200494B2 (en) * 2001-10-30 2007-04-03 Hitachi, Ltd. Method and apparatus for chromatographic data processing

Also Published As

Publication number Publication date
US20100057379A1 (en) 2010-03-04
US20110098940A1 (en) 2011-04-28
US8078411B2 (en) 2011-12-13
US7403859B2 (en) 2008-07-22
US7200494B2 (en) 2007-04-03
US20070118298A1 (en) 2007-05-24
US20050143948A1 (en) 2005-06-30
US8204696B2 (en) 2012-06-19

Similar Documents

Publication Publication Date Title
US8204696B2 (en) Method and apparatus for chromatographic data processing
US6907355B2 (en) Method and apparatus for chromatographic data processing
CN105492903B (en) Chromatograph mass spectrum analysis device
US10386346B2 (en) System and method for processing chromatogram data
US20160025691A1 (en) Chromatography/mass spectrometry data processing device
JP3888389B2 (en) Chromatographic data processing apparatus, chromatographic data processing method, and chromatographic analyzer
US8831892B2 (en) Chromatograph analyzing device
JP5211753B2 (en) Chromatographic data processor
US10481138B2 (en) Chromatogram data processing device and processing method
JP6060793B2 (en) Peak detector
US20210293763A1 (en) Peak analyzing method and waveform processing device
US10684263B2 (en) Chromatogram data processing device and program
JP2019211301A (en) Chromatograph mass analysis system and method for displaying measurement condition
JP7088306B2 (en) Chromatograph device
JP7334788B2 (en) WAVEFORM ANALYSIS METHOD AND WAVEFORM ANALYSIS DEVICE
EP3137891B1 (en) Multi-trace quantitation
CN109073614A (en) Data processing equipment and data processing method
JPH09269319A (en) Data processor for quantitative analysis
JPH11326304A (en) Data processing device for chromatograph
CN115398226A (en) Waveform information estimation method and device, and peak waveform processing method and device

Legal Events

Date Code Title Description
STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION