US20010053959A1 - Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays - Google Patents
Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays Download PDFInfo
- Publication number
- US20010053959A1 US20010053959A1 US09/918,214 US91821401A US2001053959A1 US 20010053959 A1 US20010053959 A1 US 20010053959A1 US 91821401 A US91821401 A US 91821401A US 2001053959 A1 US2001053959 A1 US 2001053959A1
- Authority
- US
- United States
- Prior art keywords
- profile
- time
- optical
- dependent measurement
- congenital
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/483—Physical analysis of biological material
- G01N33/487—Physical analysis of biological material of liquid biological material
- G01N33/49—Blood
- G01N33/4905—Determining clotting time of blood
-
- C—CHEMISTRY; METALLURGY
- C07—ORGANIC CHEMISTRY
- C07K—PEPTIDES
- C07K16/00—Immunoglobulins [IGs], e.g. monoclonal or polyclonal antibodies
- C07K16/18—Immunoglobulins [IGs], e.g. monoclonal or polyclonal antibodies against material from animals or humans
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/86—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving blood coagulating time or factors, or their receptors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N35/00—Automatic analysis not limited to methods or materials provided for in any single one of groups G01N1/00 - G01N33/00; Handling materials therefor
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/40—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for data related to laboratory analysis, e.g. patient specimen analysis
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2800/00—Detection or diagnosis of diseases
- G01N2800/22—Haematology
- G01N2800/226—Thrombotic disorders, i.e. thrombo-embolism irrespective of location/organ involved, e.g. renal vein thrombosis, venous thrombosis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2800/00—Detection or diagnosis of diseases
- G01N2800/38—Pediatrics
- G01N2800/385—Congenital anomalies
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
Definitions
- Blood clots are the end product of a complex chain reaction where proteins form an enzyme cascade acting as a biologic amplification system. This system enables relatively few molecules of initiator products to induce sequential activation of a series of inactive proteins, known as factors, culminating in the production of the fibrin clot. Mathematical models of the kinetics of the cascade's pathways have been previously proposed.
- Thrombosis and hemostasis testing is the in vitro study of the ability of blood to form clots and to break clots in vivo.
- Coagulation (hemostasis) assays began as manual methods where clot formation was observed in a test tube either by tilting the tube or removing fibrin strands by a wire loop. The goal was to determine if a patient's blood sample would clot after certain materials were added. It was later determined that the amount of time from initiation of the reaction to the point of clot formation in vitro is related to congenital disorders, acquired disorders, and therapeutic monitoring.
- Two assays are widely used to screen for abnormalities in the coagulation system, although several other screening assays can be used, e.g. protein C, fibrinogen, protein S and/or thrombin time. If screening assays show an abnormal result, one or several additional tests are needed to isolate the exact source of the abnormality.
- the PT and APTT assays rely primarily upon measurement of time required for clot time, although some variations of the PT also use the amplitude of the change in optical signal in estimating fibrinogen concentration.
- Blood coagulation is affected by administration of drugs, in addition to the vast array of internal factors and proteins that normally influence clot formation.
- heparin is a widely-used therapeutic drug that is used to prevent thrombosis following surgery or under other conditions, or is used to combat existing thrombosis.
- the administration of heparin is typically monitored using the APTT assay, which gives a prolonged clot time in the presence of heparin. Clot times for PT assays are affected to a much smaller degree. Since a number of other plasma abnormalities may also cause prolonged APTT results, the ability to discriminate between these effectors from screening assay results may be clinically significant.
- Baumann, et al [4] showed that a ratio of two coefficients was unique for a select group of blood factor deficiencies when fibrinogen was artificially maintained by addition of exogenous fibrinogen to a fixed concentration, and that same ratio also correlates heparin to FII deficiency and FXa deficiencies.
- the present invention makes it possible to predict a congenital or acquired imbalance or therapeutic condition for clinical samples from a time-dependent measurement profile without artificial manipulation of samples.
- the present invention was conceived of and developed for predicting the presence of congenital or acquired imbalances or therapeutic conditions of an unknown sample based on one or more time-dependent measurement profiles, such as optical time-dependent measurement profiles, where a set of predictor variables are provided which define characteristics of profile, and where in turn a model is derived that represents the relationship between a congenital or acquired imbalance or therapeutic condition and the set of predictor variables (so as to, in turn, utilize this model to predict the existence of the congenital or acquired imbalance or therapeutic condition in the unknown sample).
- time-dependent measurement profiles such as optical time-dependent measurement profiles
- the present invention is directed to a method and apparatus for predicting the presence of at least one congenital or acquired imbalance or therapeutic condition from at least one time-dependent measurement profile.
- the method and apparatus include a) performing at least one assay on an unknown sample and measuring a respective property over time so as to derive a time-dependent measurement profile, b) defining a set of predictor variables which sufficiently define the data of the time-dependent profile, c) deriving a model that represents the relationship between a diagnostic output and the set of predictor variables, and d) utilizing the model to predict the existence of a congenital or acquired imbalance or therapeutic condition in the unknown sample relative to the diagnostic output.
- training data is provided by performing a plurality of assays on known samples
- the model is a multilayer perceptron
- the relationship between the diagnostic output and the set of predictor variables is determined by at least one algorithm
- the at least one algorithm is a back propagation learning algorithm.
- the relationship between the diagnostic output and the set of predictor variables is derived by a set of statistical equations.
- time-dependent measurement profiles are derived, which time-dependent measurement profiles can be optical time-dependent measurement profiles such as ones provided by a automated analyzer for thrombosis and hemostasis, where a plurality of optical measurements are taken over time, and where the plurality of optical measurements are normalized.
- the optical profiles can include one or more of a PT profile, a fibrinogen profile, an APTT profile, a TT profile, a protein C profile, a protein S profile and a plurality of other assays associated with congenital or acquired imbalances or therapeutic conditions.
- FIG. 1 is a general neuron diagram relating to the embodiment of the present invention utilizing a neural network
- FIG. 2 is a diagram of a multilayer perceptron for predicting congenital or acquired imbalances or therapeutic conditions, relating to the neural network embodiment of the present invention
- FIG. 3 is an optical profile with first and second derivatives of a normal clotting sample
- FIG. 4 is an illustration of two learning curves
- FIG. 5 is an illustration of an unstable learning curve
- FIG. 6 is a graph showing a comparison of training and cross-validation learning curves
- FIG. 7 is a graph showing a comparison of training error for training tolerances of 0.0 and 0.1;
- FIG. 8 is a ROC illustrating the effect of decision boundary on classification
- FIG. 9 is a Table comparing hidden layer size with prediction error
- FIG. 10 is a receiver operator characteristic plot related to predicting an abnormality in relation to Factor VIII
- FIG. 11 is a graph demonstrating the ability to predict actual Factor VIII activity
- FIG. 12 is a receiver operator characteristic plot related to predicting an abnormality in relation to Factor X.
- FIG. 13 is a chart listing examples of predictor variables for use in the present invention.
- both a method and apparatus are provided for predicting the presence of at least one congenital or acquired imbalance or therapeutic condition.
- one or more time-dependent measurements are performed on an unknown sample.
- time-dependent measurement is referred to herein to include measurements derived from assays (e.g. PT, APTT, fibrinogen, protein C, protein S, TT, ATIII, plasminogen and factor assays).
- assays e.g. PT, APTT, fibrinogen, protein C, protein S, TT, ATIII, plasminogen and factor assays.
- unknown sample and “clinical sample” refer to a sample, such as one from a medical patient, where a congenital or acquired imbalance or therapeutic condition associated with thrombosis/hemostasis is not known (or, if suspected, has not been confirmed).
- a coagulation property is measured over time so as to derive a time-dependent measurement profile.
- the time-dependent measurement is an optical measurement for deriving an optical profile.
- a PT profile, a fibrinogen profile, a TT profile, an APTT profile and/or variations thereof can be provided where, an unknown sample is analyzed for clot formation based on light transmittance over time through the unknown sample.
- two (or more) optical profiles are provided, such as both a PT profile and an APTT profile.
- a set of predictor variables are defined which sufficiently define the data of the time-dependent profile.
- One or more predictor variables comprise the set. And, in one embodiment, three or more, and in a preferred embodiment, four or more predictor variables were found to desirably make up the set.
- the characteristics of the time-dependent measurement profile could best be defined by one or more predictor variables, including the minimum of the first derivative of the optical profile, the time index of this minimum, the minimum of the second derivative of the optical profile, the time index of this minimum, the maximum of the second derivative, the time index of this maximum, the overall change in transmittance during the time-dependent measurement, clotting time, slope of the optical profile prior to clot formation, and slope of the optical profile after clot formation.
- a model is derived which represents the relationship between a congenital or acquired imbalance or therapeutic condition and the set of predictor variables.
- This model can be derived from a neural network in one embodiment of the present invention. In another embodiment, the model is derived via a set of statistical equations.
- Neural networks represent a branch of artificial intelligence that can be used to learn and model complex, unknown systems given some known data from which it can train.
- neural networks that make them an attractive alternative for modeling complex systems are:
- Neural networks are formed from multiple layers of interconnected neurons like that shown in FIG. 1. Each neuron has one output and receives input i 1 . . . i n from multiple other neurons over connecting links, or synapses. Each synapse is associated with a synaptic weight, w j .
- An adder ⁇ or linear combiner sums the products of the input signals and synaptic weights i j *w j .
- the linear combiner output sum, and ⁇ 1 (a threshold which lowers or a bias which raises the output) are the input to the activation function f().
- the synaptic weights are learned by adjusting their values through a learning algorithm.
- the model is utilized to predict the existence of a congenital or acquired imbalance or therapeutic condition in the unknown sample relative to the time-dependent measurement profile(s).
- a congenital or acquired imbalance or therapeutic condition can be predicted.
- Conditions which can be predicted as being abnormal in the present invention can include, among others, a) factor deficiencies, e.g. fibrinogen, Factors II, V, VII, VIII, IX, X, XI and XII, as well as ATIII, plasminogen, protein C, protein S, etc., b) therapeutic conditions, e.g.
- the method is performed on an automated analyzer.
- the time-dependent measurement profile such as an optical data profile
- the time-dependent measurement profile can be provided automatically by the automated analyzer, where the unknown sample is automatically removed by an automated probe from a sample container to a test well, one or more reagents are automatically added to the test well so as to initiate the reaction within the sample.
- a property over time is automatically optically monitored so as to derive the optical profile.
- the predicted congenital or therapeutic condition can be automatically stored in a memory of an automated analyzer and/or displayed on the automated analyzer, such as on a computer monitor, or printed out on paper.
- the predicted congenital or acquired imbalance or therapeutic condition is an abnormal condition
- one or more assays for confirming the existence of the abnormal condition are performed on the automated analyzer.
- the one or more confirming assays are automatically ordered and performed on the analyzer once the predicted condition is determined, with the results of the one or more confirming assays being stored in a memory of the automated analyzer and/or displayed on the analyzer.
- This example shows a set of predictor variables that adequately describe screening assay optical profiles, develops an optimal neural network design, and determines the predictive capabilities of an abnormal condition associated with thrombosis/hemostasis (in this case for the detection of heparin) with a substantial and well-quantified test data set.
- SimplastinTM L, PlatelinTM L, calcium chloride solution (0.025 M), imidazole buffer were obtained from Organon Teknika Corporation, Durham, N.C., 27712, USA. All plasma specimens were collected in 3.2% or 3.8% sodium citrate in the ratio of one part anticoagulant to nine parts whole blood. The tubes were centrifuged at 2000 g for 30 minutes and then decanted into polypropylene tubes and stored at ⁇ 80° C. until evaluated. 757 specimens were prepared from 200 samples.
- Samples represented normal patients, a variety of deficiencies, and therapeutic conditions, of the specimen population 216 were positive for heparin determined by a heparin concentration greater than 0.05 units/ml measured with a chromogenic assay specific for heparin
- PT and APTT screening assays were performed on each specimen utilizing two automated analyzers (MDATM 180s) and multiple reagent and plasma vials (Organon Teknika Corporation, Durham N.C. 27712, USA ) over a period of five days.
- MDATM 180s automated analyzers
- multiple reagent and plasma vials Organon Teknika Corporation, Durham N.C. 27712, USA
- clot-based coagulation assays are performed by an automated optically-based analyzer such as the MDA 180, data are collected over time that represents the normalized level of light transmission through a sample as a clot forms (the optical profile). As the fibrin clot forms, the transmission of light is decreased. The optical profile was stored from each test.
- MLP multilayer perceptron
- a normal optical profile is shown in FIG. 3.
- the set of predictor variables were chosen with the intent of describing optical profiles as completely as possible with a minimum number of variables. They are summarized in Table 1 where t is time from initiation of reaction, T is normalized light transmission through the reaction mixture, and pv jk is the kth predictor variable of assay j.
- the predictor variables were scaled to values between 0 and 1, based on the range of values observed for each variable for assay type k
- i j f ( pv jk , ( pv j ⁇ n,k ) min , ( pv j ⁇ n,k ) max ).
- the input variable set includes i 1 . . . 7 for both a PT assay and APTT assay for each specimen.
- heparin samples with results of greater than 0.05 units/ml were considered positive and assigned a value of 1 while negative samples were assigned a value of 0.
- the error-correction learning rule is an iterative process used to update the synaptic weights by a method of gradient descent in which the network minimizes the error as pattern associations (known input-output pairs) in the training set are presented to the network Each cycle through the training set is known as an epoch. The order or presentation of the pattern associations was the same for all epochs.
- the learning algorithm consists of six steps which make up the forward pass and the backward pass. In the forward pass, the hidden layer neuron activations are first determined
- h is the vector of hidden-layer neurons
- i the vector of input-layer neurons
- W1 the weight matrix between the input and hidden layers
- Fo the activation function
- d i s the desired output. If each element of e o is less than some predefined training error tolerance vector TE tol , than the weights are not updated during that pass and the process continues with the next pattern association. A training error tolerance of 0.1 was used in all experiments unless otherwise specified. Otherwise, the local gradient at the output layer is then computed:
- [0066] is the learning rate
- ⁇ is the momentum factor
- m is the learning iteration.
- the first layer of weights is adjusted in a similar manner
- the forward pass and backward pass are repeated for all of the pattern associations in the training set, referred to as an epoch, 1000 times.
- the trained network is applied to the cross-validation set.
- the learning curve is defined as the plot of E versus epoch
- the percent classification, ⁇ describes the percent of the total test set (training and cross-validation) that is correctly classified based on some defined decision boundary, ⁇ .
- Receiver-Operating Characteristic (ROC) plots have also been utilized to describe trained networks' ability to discriminate between the alternative possible outcome states. In these plots, measures of sensitivity and specificity are shown for a complete range of decision boundaries.
- FIG. 9 shows the effect of the hidden layer size on the training and cross validation error and the percent correct classification for the optimal decision boundary, defined as the decision boundary which yielded the lowest total number of false positives and false negatives from the total test set.
- the hidden layer size As the hidden layer size is increased, the error is decreased. However, the ability to generalize does not increase after a hidden layer size of 6. The most significant benefit in terms of both error and percentage correct classification is between 4 and 6. A hidden layer size of 6 was used for the remainder of the experiments.
- FIG. 4 shows the learning curves for two of the best combinations of parameters.
- FIG. 5 shows an example learning curve when the learning rate is so high it leads to oscillations and convergence to a higher E.
- the network converged to a lower E and as ⁇ 1, the rate of convergence improved.
- value of E converged too increased and oscillations increased.
- ⁇ 1 exacerbated the oscillations.
- FIG. 8 shows the ROC plot for networks trained with the predictor variables from each of the two screening assays with that of them combined.
- the hidden layer size was 3. While using the data from one assay does lead to some success, using the information from both assays makes a significant improvement in the ability of the network to correctly predict the presence of heparin.
- This graph indicates that a 90% true positive proportion can be achieved with a false positive proportion of 15%.
- a 60-70% true positive proportion can be achieved with a false positive proportion of approximately 15%.
- FIG. 10 is a receiver operator characteristic plot related to predicting an abnormality in relation to Factor VIII.
- everything below 30% activity was indicated as positive, and everything above 30% was indicated as negative. Cutoff values other than 30% could also be used.
- the activity percentage has a known accuracy of approximately + or ⁇ 10%.
- the actual percent activity was utilized as the output.
Abstract
A method and apparatus are disclosed for predicting the presence of at least one congenital or acquired imbalance or therapeutic condition associated with thrombosis/hemostasis from at least one time-dependent measurement profile. At least one time-dependent measurement on an unknown sample is performed and a respective property of said sample is measured over time so as to derive a time-dependent measurement profile. A set of a plurality of predictor variables are defined which sufficiently define the data of the time-dependent measurement profile. A model is then derived that represents the relationship between the congenital or acquired imbalance or therapeutic condition, and the set of predictor variables. Subsequently, the model is utilized to predict the existence of the congenital or acquired imbalance or therapeutic condition in the unknown sample.
Description
- This application is a continuation of U.S. patent application Ser. No. 08/477,389 to Givens et al. filed Jun. 7, 1995, the subject matter of which is incorporated herein by reference. This application is also related to the following publications, the subject matter of each also being incorporated herein by reference:
- 1. B. Pohl, C. Beringer, M. Bomhard, F. Keller, The quick machine—a mathematical model for the extrinsic activation of coagulation,Haemostasis, 24, 325-337 (1994).
- 2. J. Brandt, D. Triplett, W. Rock, E. Bovill, C. Arkin, Effect of lupus anticoagulants on the activated partial thromboplastin time,Arch Pathol Lab Med, 115, 109-14 (1991).
- 3. I. Talstad, Which coagulation factors interfere with the one-stage prothrombin time?,Haemostasis, 23, 19-25 (1993).
- 4. P. Baumann, T. Jurgensen, C. Heuck, Computerized analysis of the in vitro activation of the plasmatic clotting system,Haemostasis, 19, 309-321 (1989).
- 5. C. Heuck, P. Baumann, Kinetic analysis of the clotting system in the presence of heparin and depolymerized heparin,Haemostasis, 21, 10-18 (1991).
- 6. M. Astion and P. Wilding, The application of backpropagation neural networks to problems in pathology and laboratory medicine,Arch Pathol Lab Med, 116, 995-1001 (1992).
- 7. M. Astion, M. Wener, R. Thomas, G. Hunder, and D. Bloch, Overtraining in neural networks that interpret clinical data,Clinical Chemistry, 39, 1998-2004 (1993).
- 8. J Furlong, M. Dupuy, and J. Heinsimer, Neural network analysis of serial cardiac enzyme data,A.J.C.P., 96, 134-141 (1991).
- 9. W. Dassen, R. Mulleneers, J. Smeets, K. den Dulk, F. Cruz, P. Brugada, and H. Wellens, Self-learning neural networks in electrocardiography,J. Electrocardiol, 23, 200-202 (1990).
- 10. E. Baum and D. Haussler, What size net gives valid generalization?Advances in Neural Information Processing Systems, Morgan Kauffman Publishers, San Mateo, Calif., 81-90 (1989).
- 11. A. Blum,Neural Networks in C++, John Wiley & Sons, New York, (1992).
- 12. S. Haykin,Neural Networks A Comprehensive Foundation, Macmillan College Publishing Company, New York, (1994).
- 13. J. Swets, Measuring the accuracy of diagnostic systems,Science, 240, 1285-1293 (1988)
- 14. M. Zweig and G. Campbell, Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicine,Clinical Chemistry, 39, 561-577 (1993)
- 15. D. Bluestein, L. Archer, The sensitivity, specificity and predictive value of diagnostic information: a guide for clinicians,Nurse Practitioner, 16, 39-45 (1991).
- 16. C. Schweiger, G. Soeregi, S. Spitzauer, G. Maenner, and A. Pohl, Evaluation of laboratory data by conventional statistics and by three types of neural networks,Clinical Chemistry, 39, 1966-1971 (1993).
- Blood clots are the end product of a complex chain reaction where proteins form an enzyme cascade acting as a biologic amplification system. This system enables relatively few molecules of initiator products to induce sequential activation of a series of inactive proteins, known as factors, culminating in the production of the fibrin clot. Mathematical models of the kinetics of the cascade's pathways have been previously proposed.
- In [1], a dynamic model of the extrinsic coagulation cascade was described where data were collected for 20 samples using quick percent, activated partial thromboplastin time (APTT), thrombin time (TT), fibrinogen, factor(F) II, FV, FVII, FX, anti-thrombin III (ATIII), and factor degradation product (FDP) assays. These data were used as input to the model and the predictive output compared to actual recovered prothrombin time (PT) screening assay results. The model accurately predicted the PT result in only 11 of 20 cases. These coagulation cascade models demonstrate: (1) the complexity of the clot formation process, and (2) the difficulty in associating PT clot times alone with specific conditions.
- Thrombosis and hemostasis testing is the in vitro study of the ability of blood to form clots and to break clots in vivo. Coagulation (hemostasis) assays began as manual methods where clot formation was observed in a test tube either by tilting the tube or removing fibrin strands by a wire loop. The goal was to determine if a patient's blood sample would clot after certain materials were added. It was later determined that the amount of time from initiation of the reaction to the point of clot formation in vitro is related to congenital disorders, acquired disorders, and therapeutic monitoring. In order to remove the inherent variability associated with the subjective endpoint determinations of manual techniques, instrumentation has been developed to measure clot time, based on (1) electromechanical properties, (2) clot elasticity, (3) light scattering, (4) fibrin adhesion, and (5) impedance. For light scattering methods, data is gathered that represents the transmission of light through the specimen as a function of time (an optical time-dependent measurement profile).
- Two assays, the PT and APTT, are widely used to screen for abnormalities in the coagulation system, although several other screening assays can be used, e.g. protein C, fibrinogen, protein S and/or thrombin time. If screening assays show an abnormal result, one or several additional tests are needed to isolate the exact source of the abnormality. The PT and APTT assays rely primarily upon measurement of time required for clot time, although some variations of the PT also use the amplitude of the change in optical signal in estimating fibrinogen concentration.
- Blood coagulation is affected by administration of drugs, in addition to the vast array of internal factors and proteins that normally influence clot formation. For example, heparin is a widely-used therapeutic drug that is used to prevent thrombosis following surgery or under other conditions, or is used to combat existing thrombosis. The administration of heparin is typically monitored using the APTT assay, which gives a prolonged clot time in the presence of heparin. Clot times for PT assays are affected to a much smaller degree. Since a number of other plasma abnormalities may also cause prolonged APTT results, the ability to discriminate between these effectors from screening assay results may be clinically significant.
- Using a sigmoidal curve fit to a profile, Baumann, et al [4] showed that a ratio of two coefficients was unique for a select group of blood factor deficiencies when fibrinogen was artificially maintained by addition of exogenous fibrinogen to a fixed concentration, and that same ratio also correlates heparin to FII deficiency and FXa deficiencies. However, the requirement for artificially fixed fibrinogen makes this approach inappropriate for analysis of clinical specimens. The present invention makes it possible to predict a congenital or acquired imbalance or therapeutic condition for clinical samples from a time-dependent measurement profile without artificial manipulation of samples.
- The present invention was conceived of and developed for predicting the presence of congenital or acquired imbalances or therapeutic conditions of an unknown sample based on one or more time-dependent measurement profiles, such as optical time-dependent measurement profiles, where a set of predictor variables are provided which define characteristics of profile, and where in turn a model is derived that represents the relationship between a congenital or acquired imbalance or therapeutic condition and the set of predictor variables (so as to, in turn, utilize this model to predict the existence of the congenital or acquired imbalance or therapeutic condition in the unknown sample).
- The present invention is directed to a method and apparatus for predicting the presence of at least one congenital or acquired imbalance or therapeutic condition from at least one time-dependent measurement profile. The method and apparatus include a) performing at least one assay on an unknown sample and measuring a respective property over time so as to derive a time-dependent measurement profile, b) defining a set of predictor variables which sufficiently define the data of the time-dependent profile, c) deriving a model that represents the relationship between a diagnostic output and the set of predictor variables, and d) utilizing the model to predict the existence of a congenital or acquired imbalance or therapeutic condition in the unknown sample relative to the diagnostic output. In one embodiment, training data is provided by performing a plurality of assays on known samples, the model is a multilayer perceptron, the relationship between the diagnostic output and the set of predictor variables is determined by at least one algorithm, and the at least one algorithm is a back propagation learning algorithm. In a second embodiment of the present invention, the relationship between the diagnostic output and the set of predictor variables is derived by a set of statistical equations.
- Also in the present invention, a plurality of time-dependent measurement profiles are derived, which time-dependent measurement profiles can be optical time-dependent measurement profiles such as ones provided by a automated analyzer for thrombosis and hemostasis, where a plurality of optical measurements are taken over time, and where the plurality of optical measurements are normalized. The optical profiles can include one or more of a PT profile, a fibrinogen profile, an APTT profile, a TT profile, a protein C profile, a protein S profile and a plurality of other assays associated with congenital or acquired imbalances or therapeutic conditions.
- FIG. 1 is a general neuron diagram relating to the embodiment of the present invention utilizing a neural network;
- FIG. 2 is a diagram of a multilayer perceptron for predicting congenital or acquired imbalances or therapeutic conditions, relating to the neural network embodiment of the present invention;
- FIG. 3 is an optical profile with first and second derivatives of a normal clotting sample;
- FIG. 4 is an illustration of two learning curves;
- FIG. 5 is an illustration of an unstable learning curve;
- FIG. 6 is a graph showing a comparison of training and cross-validation learning curves;
- FIG. 7 is a graph showing a comparison of training error for training tolerances of 0.0 and 0.1;
- FIG. 8 is a ROC illustrating the effect of decision boundary on classification;
- FIG. 9 is a Table comparing hidden layer size with prediction error;
- FIG. 10 is a receiver operator characteristic plot related to predicting an abnormality in relation to Factor VIII;
- FIG. 11 is a graph demonstrating the ability to predict actual Factor VIII activity;
- FIG. 12 is a receiver operator characteristic plot related to predicting an abnormality in relation to Factor X; and
- FIG. 13 is a chart listing examples of predictor variables for use in the present invention.
- In the present invention, both a method and apparatus are provided for predicting the presence of at least one congenital or acquired imbalance or therapeutic condition. As one of the first seeps of the method, one or more time-dependent measurements are performed on an unknown sample. The term “time-dependent measurement” is referred to herein to include measurements derived from assays (e.g. PT, APTT, fibrinogen, protein C, protein S, TT, ATIII, plasminogen and factor assays). The terms “unknown sample” and “clinical sample” refer to a sample, such as one from a medical patient, where a congenital or acquired imbalance or therapeutic condition associated with thrombosis/hemostasis is not known (or, if suspected, has not been confirmed). In the present invention, a coagulation property is measured over time so as to derive a time-dependent measurement profile. In a preferred embodiment, the time-dependent measurement is an optical measurement for deriving an optical profile. For example, a PT profile, a fibrinogen profile, a TT profile, an APTT profile and/or variations thereof can be provided where, an unknown sample is analyzed for clot formation based on light transmittance over time through the unknown sample. In another preferred embodiment, two (or more) optical profiles are provided, such as both a PT profile and an APTT profile.
- After the time-dependent measurement profiles are provided, a set of predictor variables are defined which sufficiently define the data of the time-dependent profile. One or more predictor variables comprise the set. And, in one embodiment, three or more, and in a preferred embodiment, four or more predictor variables were found to desirably make up the set. It was found that the characteristics of the time-dependent measurement profile could best be defined by one or more predictor variables, including the minimum of the first derivative of the optical profile, the time index of this minimum, the minimum of the second derivative of the optical profile, the time index of this minimum, the maximum of the second derivative, the time index of this maximum, the overall change in transmittance during the time-dependent measurement, clotting time, slope of the optical profile prior to clot formation, and slope of the optical profile after clot formation.
- After defining the set of predictor variables, a model is derived which represents the relationship between a congenital or acquired imbalance or therapeutic condition and the set of predictor variables. This model can be derived from a neural network in one embodiment of the present invention. In another embodiment, the model is derived via a set of statistical equations.
- Neural networks represent a branch of artificial intelligence that can be used to learn and model complex, unknown systems given some known data from which it can train. Among the features of neural networks that make them an attractive alternative for modeling complex systems are:
- 1. They can handle noisy data well and recognize patterns even when some of the input data are obscured or missing.
- 2. It is unnecessary to determine what factors are relevant a priori since the network will determine during the training phase what data are relevant, assuming there are at least some meaningful parameters in the set.
- Neural networks are formed from multiple layers of interconnected neurons like that shown in FIG. 1. Each neuron has one output and receives input i1 . . . in from multiple other neurons over connecting links, or synapses. Each synapse is associated with a synaptic weight, wj. An adder Σ or linear combiner sums the products of the input signals and synaptic weights ij*wj. The linear combiner output sum, and θ1 (a threshold which lowers or a bias which raises the output) are the input to the activation function f(). The synaptic weights are learned by adjusting their values through a learning algorithm.
- After deriving the model, whether based on neural networks or statistical equations, the model is utilized to predict the existence of a congenital or acquired imbalance or therapeutic condition in the unknown sample relative to the time-dependent measurement profile(s). As such, a congenital or acquired imbalance or therapeutic condition can be predicted. Conditions which can be predicted as being abnormal in the present invention can include, among others, a) factor deficiencies, e.g. fibrinogen, Factors II, V, VII, VIII, IX, X, XI and XII, as well as ATIII, plasminogen, protein C, protein S, etc., b) therapeutic conditions, e.g. heparin, coumadin, etc., and c) conditions such as lupus anticoagulant. In one embodiment of the present invention, the method is performed on an automated analyzer. The time-dependent measurement profile, such as an optical data profile, can be provided automatically by the automated analyzer, where the unknown sample is automatically removed by an automated probe from a sample container to a test well, one or more reagents are automatically added to the test well so as to initiate the reaction within the sample. A property over time is automatically optically monitored so as to derive the optical profile. The predicted congenital or therapeutic condition can be automatically stored in a memory of an automated analyzer and/or displayed on the automated analyzer, such as on a computer monitor, or printed out on paper. As a further feature of the invention, if the predicted congenital or acquired imbalance or therapeutic condition is an abnormal condition, then one or more assays for confirming the existence of the abnormal condition are performed on the automated analyzer. In fact, in a preferred embodiment, the one or more confirming assays are automatically ordered and performed on the analyzer once the predicted condition is determined, with the results of the one or more confirming assays being stored in a memory of the automated analyzer and/or displayed on the analyzer.
- This example shows a set of predictor variables that adequately describe screening assay optical profiles, develops an optimal neural network design, and determines the predictive capabilities of an abnormal condition associated with thrombosis/hemostasis (in this case for the detection of heparin) with a substantial and well-quantified test data set.
- Simplastin™ L, Platelin™ L, calcium chloride solution (0.025 M), imidazole buffer were obtained from Organon Teknika Corporation, Durham, N.C., 27712, USA. All plasma specimens were collected in 3.2% or 3.8% sodium citrate in the ratio of one part anticoagulant to nine parts whole blood. The tubes were centrifuged at 2000 g for 30 minutes and then decanted into polypropylene tubes and stored at −80° C. until evaluated. 757 specimens were prepared from 200 samples. These specimens were tested by the following specific assays: FII, FV, FVII, FVIII, FIX, FX, FXI, FXII, heparin, fibrinogen, plasminogen, protein C, and AT-III. Samples represented normal patients, a variety of deficiencies, and therapeutic conditions, of the
specimen population 216 were positive for heparin determined by a heparin concentration greater than 0.05 units/ml measured with a chromogenic assay specific for heparin The remaining specimens, classified as heparin-negative, included normal specimens, a variety of single or multiple factor deficiencies, and patients receiving other therapeutic drugs. Positive heparin samples ranged to 0.54 units/ml. - PT and APTT screening assays were performed on each specimen utilizing two automated analyzers (MDA™ 180s) and multiple reagent and plasma vials (Organon Teknika Corporation, Durham N.C. 27712, USA ) over a period of five days. When clot-based coagulation assays are performed by an automated optically-based analyzer such as the
MDA 180, data are collected over time that represents the normalized level of light transmission through a sample as a clot forms (the optical profile). As the fibrin clot forms, the transmission of light is decreased. The optical profile was stored from each test. - The network configuration chosen, a multilayer perceptron (MLP) maps input predictor variables from the PT and APTT screening assays to one output variable (see FIG. 2) which represents a single specified condition. A similar network was also employed for PT-only variables and APTT-only variables. This specific MLP consists of three layers: the input layer, one hidden layer, and the output layer.
- A normal optical profile is shown in FIG. 3. The set of predictor variables were chosen with the intent of describing optical profiles as completely as possible with a minimum number of variables. They are summarized in Table 1 where t is time from initiation of reaction, T is normalized light transmission through the reaction mixture, and pvjk is the kth predictor variable of assay j.
- The predictor variables were scaled to values between 0 and 1, based on the range of values observed for each variable for assay type k
- i j =f(pv jk, (pv j−n,k)min, (pv j−n,k)max).
- The input variable set includes i1 . . . 7 for both a PT assay and APTT assay for each specimen. For known output variable values, heparin samples with results of greater than 0.05 units/ml were considered positive and assigned a value of 1 while negative samples were assigned a value of 0.
- As the ratio of training set sample to the number of weights in a network decreases, the probability of generalizing decreases, reducing the confidence that the network will lead to correct classification of future samples taken from the same distribution as the training set. Thus, small samples sizes, then can lead to artificially high classification rates. This phenomenon is known as overtraining. In order to achieve a true accuracy rate of 80%, a guideline for the number of samples in the training set is approximately five times the number of weights in the network. For most of this work, a 14-6-1 network was used, leading to an upward bound on the sample size of 0(450). To monitor and evaluate the performance of the network and its ability to generalize, a cross-validation set is processed at the end of each training epoch. This cross-validation set is a randomly determined subset of the known test set that is excluded from the training set.
- Once the input predictor variables and output values were determined for all specimen optical profiles, the 757 sets of data were randomly distributed into two groups: 387 were used in the training set and 370 were used in the cross-validation set. These same two randomly determined sets were used throughout all the experiments.
- All synaptic weights and threshold values were initialized at the beginning of each training session to small random numbers.
- The error-correction learning rule is an iterative process used to update the synaptic weights by a method of gradient descent in which the network minimizes the error as pattern associations (known input-output pairs) in the training set are presented to the network Each cycle through the training set is known as an epoch. The order or presentation of the pattern associations was the same for all epochs. The learning algorithm consists of six steps which make up the forward pass and the backward pass. In the forward pass, the hidden layer neuron activations are first determined
- h=F(iW1+θh)
-
- Then the output-layer neurons are computed
- o=F(hW2+θo)
- where o represents the output layer, h the hidden layer and W2 the matrix of synapses connecting the hidden layer and output layers. The backward pass begins with the computation of the output layer error
- e o=(o−d),
- where d i s the desired output. If each element of eo is less than some predefined training error tolerance vector TEtol, than the weights are not updated during that pass and the process continues with the next pattern association. A training error tolerance of 0.1 was used in all experiments unless otherwise specified. Otherwise, the local gradient at the output layer is then computed:
- g o =o(1−o)e o.
- Next, the hidden-layer local gradient is computed:
- g h =h(1−h)W2g o.
- Once the hidden layer error is calculated, the second layer of weights is adjusted
- W2m =W2m−1 +ΔW2
- where
- ΔW2=ηhg o +γΔW2m−1.
- is the learning rate, γ is the momentum factor, and m is the learning iteration. The first layer of weights is adjusted in a similar manner
- W1m =W1m−1 +ΔW1
- where
- ΔW1=ηie+γΔW1m−1.
- The forward pass and backward pass are repeated for all of the pattern associations in the training set, referred to as an epoch, 1000 times. At the end of each epoch, the trained network is applied to the cross-validation set.
-
- The learning curve is defined as the plot of E versus epoch The percent classification, φ, describes the percent of the total test set (training and cross-validation) that is correctly classified based on some defined decision boundary, β. Receiver-Operating Characteristic (ROC) plots have also been utilized to describe trained networks' ability to discriminate between the alternative possible outcome states. In these plots, measures of sensitivity and specificity are shown for a complete range of decision boundaries. The sensitivity, or true-positive fraction is defined as
-
- These ROC plots represent a common tool for evaluating clinical laboratory test performance.
- Using the test set described, experiments were performed to determine if the presence of heparin could be predicted with this method. First, experiments were conducted to determine optimal error-correction backpropagation learning parameters: (1) hidden layer size, (2) learning rate, and (3) momentum. Additional experiments were also conducted to compare the performance of networks based on PT and APTT assays alone with that of one combining the results of both, the effect of the training error tolerance, and the decision boundary selection.
- FIG. 9 shows the effect of the hidden layer size on the training and cross validation error and the percent correct classification for the optimal decision boundary, defined as the decision boundary which yielded the lowest total number of false positives and false negatives from the total test set. As the hidden layer size is increased, the error is decreased. However, the ability to generalize does not increase after a hidden layer size of 6. The most significant benefit in terms of both error and percentage correct classification is between 4 and 6. A hidden layer size of 6 was used for the remainder of the experiments.
- A series of experiments were conducted with η{0.01,0.1,0.5,0.9} and γ={0.0,0.1,0.5,0.9}. FIG. 4 shows the learning curves for two of the best combinations of parameters. FIG. 5 shows an example learning curve when the learning rate is so high it leads to oscillations and convergence to a higher E. In general, as η→0 the network converged to a lower E and as γ→1, the rate of convergence improved. As η→1, value of E converged too increased and oscillations increased. In addition, as η→1, γ→1 exacerbated the oscillations.
- FIG. 6 shows a comparison of the learning curve for the training set and cross-validation set for η=0.5 and γ=0.1. It is a primary concern when developing neural networks, and it has been previously shown that it is important to look not only at the error in the training set for each cycle, but also the cross-validation error.
- FIG. 7 shows the learning curve η=0.5 and γ=0.1 and a learning tolerance of 0.0 and 0.1. These results suggest that a small learning tends to smoothen the convergence of the learning process.
- FIG. 8 shows the ROC plot for networks trained with the predictor variables from each of the two screening assays with that of them combined. In the single assay cases, the hidden layer size was 3. While using the data from one assay does lead to some success, using the information from both assays makes a significant improvement in the ability of the network to correctly predict the presence of heparin. This graph indicates that a 90% true positive proportion can be achieved with a false positive proportion of 15%. Using a single assay, a 60-70% true positive proportion can be achieved with a false positive proportion of approximately 15%.
- Similar tests were run as in Example 1. As can be seen in FIGS. 10 and 11, two training sessions were conducted for predicting a Factor VIII condition in an unknown sample. FIG. 10 is a receiver operator characteristic plot related to predicting an abnormality in relation to Factor VIII. In FIG. 10, everything below 30% activity was indicated as positive, and everything above 30% was indicated as negative. Cutoff values other than 30% could also be used. In this Example, the activity percentage has a known accuracy of approximately + or −10%. In FIG. 11, the actual percent activity was utilized as the output.
- As can be seen in FIG. 12, the method of the present invention was run similar to that as-in Example 2, where here an abnormality in Factor X concentration was predicted from unknown samples. Everything below 30% activity was indicated as positive, and everything above 30% was indicated as negative. Cutoff values other than 30% could also be used.
- The results of the cross-validation sample sets throughout the experiments indicate that the sample size was sufficient for the network to generalize. While the random distribution of the training and cross-validation sets were held constant throughout the experiments presented, other distributions have been used. These distributions, while all yielding different results, still lead to the same general conclusion.
- Many alternatives for or additions to the set of predictor variables were explored. This included coefficients of a curve fitted to the data profile, pattern recognition, and clot time-based parameters. Low order functions tend to lose information due to their poor fit, and high order functions tend to lose information in their multiple close solutions. Clot-based parameters, such as clot time, slope in the section prior to the initiation of clot formation, and afterwards, are often available, but not always (because in some samples, the clot time is not detectable). The successful results observed indicate that the set of predictor variables used are effective for predicting congenital or acquired imbalances or therapeutic conditions.
- The optimization of the network learning algorithm's parameters made significant differences in its performance. In general, performance was best with low learning rates, high momentum rates, some small training error tolerance, and a hidden layer size approximately half of the size of the input layer.
- It is to be understood that the invention described and illustrated herein is to be taken as a preferred example of the same, and that various changes in the method and apparatus of the invention may be resorted to, without departing from the spirit of the invention or scope of the claims.
Claims (40)
1. A method for predicting the presence of at least one congenital or acquired imbalance or therapeutic condition associated with thrombosis/hemostasis from at least one time-dependent measurement profile, comprising:
a) performing at least one time-dependent measurement on an unknown sample and measuring a respective property over time so as to derive a time-dependent measurement profile;
b) defining a set of a plurality of predictor variables which sufficiently define the data of the time-dependent measurement profile;
c) deriving a model that represents the relationship between the congenital or acquired imbalance or therapeutic condition, and the set of predictor variables; and
d) utilizing the model of step c) to predict the existence of the congenital or acquired imbalance or therapeutic condition in the unknown sample.
2. A method according to , wherein said at least one time-dependent measurement profile is at least one optical profile.
claim 1
3. A method according to , wherein said at least one optical profile is provided by an automated analyzer for thrombosis and hemostasis testing.
claim 2
4. A method according to , wherein a plurality of optical measurements at one or more wavelengths are taken over time so as to derive said at least one optical profile, said optical measurements corresponding to changes in light scattering and/or light absorption in the unknown sample.
claim 2
5. A method according to , wherein a plurality of optical measurements are taken over time so as to derive said at least one optical profile, and wherein said plurality of optical measurements are each normalized to a first optical measurement.
claim 2
6. A method according to , wherein in step a) said at least one optical profile is provided automatically by said analyzer, whereby said unknown sample is automatically removed by an automated probe from a sample container to a test well, one or more reagents are automatically added to said test well so as to initiate said property changes within said sample, and the development of said property over time is automatically optically monitored so as to derive said optical data profile.
claim 3
7. A method according to , wherein after step d), a predicted congenital or acquired imbalance or therapeutic condition is automatically stored in a memory of said automated analyzer and/or displayed on said automated analyzer.
claim 6
8. A method according to , wherein in step d), one or more assays for confirming the existence of said congenital or acquired imbalance or therapeutic condition is automatically performed.
claim 6
9. A method according to , wherein said one or more confirming assays are automatically ordered and performed on said analyzer, with results of said one or more assays being stored in a memory of said automated analyzer and/or displayed on said analyzer.
claim 8
10. A method according to , further comprising:
claim 1
before step a), providing a set of data from known samples, which data is used in step c) for deriving said model.
11. A method according to , wherein said data from known samples is provided by performing a plurality of assays on said known samples.
claim 10
12. A method according to , wherein said model of step c) is a neural network.
claim 10
13. A method according to , wherein said relationship in step c) is determined via at least one automated algorithm.
claim 1
14. A method according to , wherein in step a), a plurality of time-dependent measurement profiles are derived for use in step b).
claim 1
15. A method according to , wherein said plurality of time dependent measurement profiles includes at least two profiles from assays initiated with PT reagents, APTT reagents, fibrinogen reagents and TT reagents.
claim 14
16. A method according to , wherein said model is a multilayer perceptron, and wherein said at least one algorithm is a back propagation learning algorithm.
claim 13
17. A method according to , wherein said set of predictor variables includes a plurality of: a minimum of the first derivative of the profile, a time index of the minimum of the first derivative, a minimum of the second derivative of the profile, a time index of the minimum of the second derivative, a maximum of the second derivative of the profile, a time index of the maximum of the second derivative, an overall change in the coagulation parameter during the time-dependent measurement on the unknown sample, a clotting time, a slope of the profile prior to clot formation, and a slope of the profile after clot formation.
claim 1
18. A method according to , wherein three or more of said predictor variables are within said set.
claim 17
19. A method according to , wherein more than three of said predictor variables are within said set.
claim 18
20. A method according to , wherein said unknown sample is a sample from a medical patient, and wherein in step d), both said model and additional patient medical data are utilized for predicting the existence of said congenital or acquired imbalance or therapeutic condition.
claim 1
21. An apparatus for performing at least one time-dependent measurement on an unknown sample to derive at least one time-dependent measurement profile, and predicting the presence of at least one congenital or acquired imbalance or therapeutic condition associated with thrombosis/hemostasis from the at least one time-dependent measurement profile, comprising:
means for performing at least one time-dependent measurement on an unknown sample and measuring a respective property over time so as to derive a time-dependent measurement profile;
means for defining a set of a plurality of predictor variables which sufficiently define the data of the time-dependent measurement profile;
means for deriving a model that represents the relationship between the congenital or acquired imbalance or therapeutic condition, and the set of predictor variables; and
means for utilizing the model of step c) to predict the existence of the congenital or acquired imbalance or therapeutic condition in the unknown sample.
22. An apparatus according to , wherein said means for performing at least one time-dependent measurement comprises an optical system for performing at least one optical measurement over time and so as to derive an at least one optical profile.
claim 21
23. An apparatus according to , wherein said optical system is part of an automated analyzer for thrombosis and hemostasis testing.
claim 22
24. An apparatus according to , wherein said optical means comprises a means for performing a plurality of optical measurements at one or more wavelengths over time so as to derive said at least one optical profile, said optical measurements corresponding to changes in light scattering and/or light absorption in the unknown sample.
claim 22
25. An apparatus according to , wherein in said optical system, a plurality of optical measurements are taken over time so as to derive said at least one optical profile, and wherein said plurality of optical measurements are each normalized to a first optical measurement.
claim 22
26. An apparatus according to which is an automated analyzer for thrombosis and hemostasis testing, and wherein said at least one optical profile is provided automatically by said analyzer, whereby said unknown sample is automatically removed by an automated probe from a sample container to a test well, one or more reagents are automatically added to said test well so as to initiate said property changes within said sample, and the development of said property over time is automatically optically monitored so as to derive said optical data profile.
claim 23
27. An apparatus according to , further comprising at least one of a memory and a display wherein a predicted congenital or acquired imbalance or therapeutic condition is automatically stored in said memory of said automated analyzer and/or displayed on said display of said automated analyzer.
claim 26
28. An apparatus according to , further comprising means for automatically performing one or more assays for confirming the existence of said congenital or acquired imbalance or therapeutic condition.
claim 26
29. An apparatus according to , wherein said means for performing one or more confirming assays is an automatic performing means wherein said confirming assays are automatically ordered and performed on said analyzer, with results of said one or more assays being stored in a memory of said automated analyzer and/or displayed on a display of said analyzer.
claim 28
30. An apparatus according to , further comprising means for providing a set of data from known samples, which data is used in step c) for deriving said model.
claim 21
31. An apparatus according to , wherein said data from known samples is provided by said means for performing a plurality of assays on said known samples.
claim 30
32. An apparatus according to , wherein said means for deriving a model is a means for deriving a model by means of a neural network.
claim 30
33. An apparatus according to , wherein said relationship determined by said deriving means comprises a means for determining said relationship via at least one automated algorithm.
claim 21
34. An apparatus according to , wherein said means for performing at least one time-dependent measurement is capable of performing a plurality of time-dependent measurement profiles.
claim 21
35. An apparatus according to , wherein said means for performing a plurality of time dependent measurement profiles includes a means for performing at least two profiles from assays initiated with PT reagents, APTT reagents, fibrinogen reagents and TT reagents.
claim 34
36. An apparatus according to , wherein said model is a multilayer perceptron, and wherein said at least one algorithm is a back propagation learning algorithm.
claim 33
37. An apparatus according to , wherein said set of predictor variables includes a plurality of: a minimum of the first derivative of the profile, a time index of the minimum of the first derivative, a minimum of the second derivative of the profile, a time index of the minimum of the second derivative, a maximum of the second derivative of the profile, a time index of the maximum of the second derivative, an overall change in the coagulation parameter during the time-dependent measurement on the unknown sample, a clotting time, a slope of the profile prior to clot formation, and a slope of the profile after clot formation.
claim 21
38. An apparatus according to , wherein three or more of said predictor variables are within said set.
claim 37
39. An apparatus according to , wherein more than three of said predictor variables are within said set.
claim 38
40. An apparatus according to , wherein said unknown sample is a sample from a medical patient, and wherein said utilizing means comprising a means for utilizing both said model and additional patient medical data for predicting the existence of said congenital or acquired imbalance or therapeutic condition.
claim 21
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/918,214 US20010053959A1 (en) | 1995-06-07 | 2001-07-30 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US08/477,839 US5708591A (en) | 1995-02-14 | 1995-06-07 | Method and apparatus for predicting the presence of congenital and acquired imbalances and therapeutic conditions |
US08/859,773 US6101449A (en) | 1995-06-07 | 1997-05-21 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
US09/517,496 US6269313B1 (en) | 1995-06-07 | 2000-03-02 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
US09/918,214 US20010053959A1 (en) | 1995-06-07 | 2001-07-30 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US09/517,496 Continuation US6269313B1 (en) | 1995-06-07 | 2000-03-02 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
Publications (1)
Publication Number | Publication Date |
---|---|
US20010053959A1 true US20010053959A1 (en) | 2001-12-20 |
Family
ID=23897564
Family Applications (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US08/477,839 Expired - Lifetime US5708591A (en) | 1995-02-14 | 1995-06-07 | Method and apparatus for predicting the presence of congenital and acquired imbalances and therapeutic conditions |
US08/859,773 Expired - Fee Related US6101449A (en) | 1995-06-07 | 1997-05-21 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
US09/517,496 Expired - Fee Related US6269313B1 (en) | 1995-06-07 | 2000-03-02 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
US09/918,214 Abandoned US20010053959A1 (en) | 1995-06-07 | 2001-07-30 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
Family Applications Before (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US08/477,839 Expired - Lifetime US5708591A (en) | 1995-02-14 | 1995-06-07 | Method and apparatus for predicting the presence of congenital and acquired imbalances and therapeutic conditions |
US08/859,773 Expired - Fee Related US6101449A (en) | 1995-06-07 | 1997-05-21 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
US09/517,496 Expired - Fee Related US6269313B1 (en) | 1995-06-07 | 2000-03-02 | Method for predicting the presence of congenital and therapeutic conditions from coagulation screening assays |
Country Status (12)
Country | Link |
---|---|
US (4) | US5708591A (en) |
EP (1) | EP0834145B1 (en) |
JP (1) | JP3534415B2 (en) |
KR (1) | KR100462688B1 (en) |
AT (1) | ATE206226T1 (en) |
AU (1) | AU711998B2 (en) |
CA (1) | CA2221712C (en) |
DE (1) | DE69615556T2 (en) |
DK (1) | DK0834145T3 (en) |
ES (1) | ES2164888T3 (en) |
PT (1) | PT834145E (en) |
WO (1) | WO1996041291A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110016073A1 (en) * | 2008-02-14 | 2011-01-20 | Orion Diagnostica Oy | Method for predicting a future property |
Families Citing this family (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5708591A (en) * | 1995-02-14 | 1998-01-13 | Akzo Nobel N.V. | Method and apparatus for predicting the presence of congenital and acquired imbalances and therapeutic conditions |
US6429017B1 (en) * | 1999-02-04 | 2002-08-06 | Biomerieux | Method for predicting the presence of haemostatic dysfunction in a patient sample |
US6898532B1 (en) | 1995-06-07 | 2005-05-24 | Biomerieux, Inc. | Method and apparatus for predicting the presence of haemostatic dysfunction in a patient sample |
US6321164B1 (en) | 1995-06-07 | 2001-11-20 | Akzo Nobel N.V. | Method and apparatus for predicting the presence of an abnormal level of one or more proteins in the clotting cascade |
US6108493A (en) * | 1996-10-08 | 2000-08-22 | Regents Of The University Of Minnesota | System, method, and article of manufacture for utilizing implicit ratings in collaborative filters |
US6016475A (en) * | 1996-10-08 | 2000-01-18 | The Regents Of The University Of Minnesota | System, method, and article of manufacture for generating implicit ratings based on receiver operating curves |
US5842199A (en) * | 1996-10-18 | 1998-11-24 | Regents Of The University Of Minnesota | System, method and article of manufacture for using receiver operating curves to evaluate predictive utility |
US6063028A (en) | 1997-03-20 | 2000-05-16 | Luciano; Joanne Sylvia | Automated treatment selection method |
US6502040B2 (en) * | 1997-12-31 | 2002-12-31 | Biomerieux, Inc. | Method for presenting thrombosis and hemostasis assay data |
DE19858278A1 (en) * | 1998-12-17 | 2000-06-21 | Dade Behring Marburg Gmbh | Method for locating defects in the coagulation system |
EP1522860A1 (en) * | 1999-02-04 | 2005-04-13 | bioMérieux, Inc. | A method and apparatus for predicting the presence of haemostatic dysfunction in a patient sample |
JP4486260B2 (en) * | 1999-02-04 | 2010-06-23 | バイオメリュー・インコーポレイテッド | Method and apparatus for predicting the presence of hemostatic dysfunction in patient samples |
AU2004201011B2 (en) * | 1999-06-30 | 2005-12-22 | Biomerieux, Inc | Method and apparatus for presenting thrombosis and hemostasis assay data |
US6484122B1 (en) * | 2000-05-16 | 2002-11-19 | Siemens Aktiengesellschaft | System and method for evaluating characteristics for suitability in classification |
ATE302947T1 (en) * | 2000-06-09 | 2005-09-15 | Bio Merieux Inc | METHOD FOR DETECTING A LIPOPROTEIN-ACUTE PHASE PROTEIN COMPLEX |
US7179612B2 (en) | 2000-06-09 | 2007-02-20 | Biomerieux, Inc. | Method for detecting a lipoprotein-acute phase protein complex and predicting an increased risk of system failure or mortality |
US7003403B1 (en) | 2000-06-15 | 2006-02-21 | The United States Of America As Represented By The Department Of Health And Human Services | Quantifying gene relatedness via nonlinear prediction of gene |
US7107287B2 (en) * | 2000-07-27 | 2006-09-12 | Canberra Industries | Method, system and storage medium for automated independent technical review |
US6743596B1 (en) | 2000-10-27 | 2004-06-01 | Biomerieux, Inc. | Method of determining global coagulability hemostatic potential |
WO2003083490A2 (en) * | 2001-06-29 | 2003-10-09 | Biomerieux, Inc. | Method for predicting antiphospholid syndrome |
JP2005519267A (en) | 2002-02-27 | 2005-06-30 | バイオメリュー・インコーポレイテッド | Methods for diagnosing and monitoring hemostatic dysfunction, severe infections and systemic inflammatory response syndrome |
US20040039556A1 (en) * | 2002-08-22 | 2004-02-26 | Ibex Process Technology, Inc. | Filter models for dynamic control of complex processes |
AU2003267147A1 (en) | 2002-09-10 | 2004-04-30 | Placor, Inc. | Method and device for monitoring platelet function |
CN1738909A (en) | 2002-11-12 | 2006-02-22 | 贝克顿迪金森公司 | Diagnosis of sepsis or SIRS using biomarker profiles |
CA2505921A1 (en) | 2002-11-12 | 2004-05-27 | Becton, Dickinson And Company | Diagnosis of sepsis or sirs using biomarker profiles |
JP2008538007A (en) | 2005-04-15 | 2008-10-02 | ベクトン,ディッキンソン アンド カンパニー | Diagnosis of sepsis |
JP2008539438A (en) * | 2005-04-25 | 2008-11-13 | プラコア・インコーポレーテッド | Method and device for monitoring platelet function |
US20070020764A1 (en) * | 2005-07-20 | 2007-01-25 | Miller Kerry L | Method for processing chemistry and coagulation test samples in a laboratory workcell |
US7815809B2 (en) * | 2005-12-13 | 2010-10-19 | Gambro Lundia Ab | Method for conductivity calculation in a treatment fluid upstream and downstream a filtration unit in apparatuses for the blood treatment |
JP4925703B2 (en) * | 2006-03-30 | 2012-05-09 | シスメックス株式会社 | Blood coagulation analyzer |
US20100099130A1 (en) * | 2006-10-25 | 2010-04-22 | Placor Inc. | Methods and devices for monitoring platelet function |
US8669113B2 (en) | 2008-04-03 | 2014-03-11 | Becton, Dickinson And Company | Advanced detection of sepsis |
US8046175B2 (en) * | 2008-10-13 | 2011-10-25 | Actherm Inc | Analytical strip reading apparatus and the analyical strip used therein |
CN102216758B (en) * | 2008-12-05 | 2013-07-17 | 红电医学科技股份有限公司 | A testing piece reader provided with a removable firmware |
US8645306B2 (en) * | 2010-07-02 | 2014-02-04 | Idexx Laboratories, Inc. | Automated calibration method and system for a diagnostic analyzer |
US8606750B2 (en) * | 2010-11-12 | 2013-12-10 | Life Technologies Corporation | Systems and methods for laboratory assay validation or verification |
JP6211806B2 (en) * | 2013-05-31 | 2017-10-11 | 株式会社日立ハイテクノロジーズ | Automatic analyzer |
WO2018063914A1 (en) | 2016-09-29 | 2018-04-05 | Animantis, Llc | Methods and apparatus for assessing immune system activity and therapeutic efficacy |
JP6873833B2 (en) * | 2017-06-09 | 2021-05-19 | シスメックス株式会社 | Blood sample judgment method, blood sample analyzer and computer program |
GB202004047D0 (en) | 2020-03-20 | 2020-05-06 | Ainger Phill | Point of care sepsis assay device and method |
Family Cites Families (61)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3307392A (en) * | 1964-05-04 | 1967-03-07 | Research Corp | Automatic prothrombin timer apparatus and method |
US3458287A (en) * | 1965-04-29 | 1969-07-29 | Medical Laboratory Automation | Method and means of determining endpoint times in blood clotting tests |
US3658480A (en) * | 1970-04-13 | 1972-04-25 | Bio Data Corp | Coagulation timing apparatus, and method |
US4047890A (en) * | 1973-11-01 | 1977-09-13 | Bio/Data Corporation | Method and apparatus for determining deficiencies in enzymatic reactors particularly clotting factor levels in blood plasmas |
SU590665A1 (en) * | 1976-02-11 | 1978-01-30 | Новосибирский государственный медицинский институт | Method of investigating |
DE2635081C3 (en) * | 1976-08-04 | 1980-08-07 | Bio/Data Corp., Willow Grove, Pa. (V.St.A.) | Method and apparatus for determining blood plasma coagulation factor levels |
FR2364453A1 (en) * | 1976-09-08 | 1978-04-07 | Lacombe Pierre | Measurement of blood plasma coagulation time - using optical density method, suitable for oxalated or citrated plasma |
US4199748A (en) * | 1976-11-01 | 1980-04-22 | Rush-Presbyterian-St. Luke's Medical Center | Automated method and apparatus for classification of cells with application to the diagnosis of anemia |
JPS5451893A (en) * | 1977-09-30 | 1979-04-24 | Sankyo Co | Measuring of blood coagulating time |
JPS5822703B2 (en) * | 1978-08-30 | 1983-05-10 | 三共株式会社 | Blood coagulation measurement method |
US4289498A (en) * | 1979-01-08 | 1981-09-15 | Ortho Diagnostics, Inc. | One-stage prothrombin assay and compositions useful therein |
US4766083A (en) * | 1982-04-04 | 1988-08-23 | Wako Pure Chemical Industries, Ltd. | Method for the photometric determination of biological agglutination |
SU1076086A1 (en) * | 1982-12-28 | 1984-02-29 | Всесоюзный Научно-Исследовательский И Конструкторский Институт Медицинской Лабораторной Техники | Hemocoagulograph |
US4705756A (en) * | 1983-01-26 | 1987-11-10 | University Of Medicine And Dentistry Of New Jersey | Method of determining the existence and/or the monitoring of a pathological condition in a mammal |
EP0115459A3 (en) * | 1983-01-26 | 1986-11-20 | University Of Medicine And Dentistry Of New Jersey | Hematological prediction of sepsis and other disease states |
JPS59203959A (en) * | 1983-05-04 | 1984-11-19 | Toa Medical Electronics Co Ltd | Measuring device for blood clotting time |
JPH0695095B2 (en) * | 1983-11-25 | 1994-11-24 | 株式会社島津製作所 | Fully automatic blood analyzer |
DE3502878A1 (en) * | 1985-01-29 | 1986-07-31 | Boehringer Mannheim Gmbh, 6800 Mannheim | METHOD FOR DETERMINING THE FIBRINOLYTIC STATE OF PLASMA |
JPS61272655A (en) * | 1985-05-28 | 1986-12-02 | Iryo Kogaku Kenkyusho:Kk | Quantitative discrimination method for platelet agglutination power |
US4965725B1 (en) * | 1988-04-08 | 1996-05-07 | Neuromedical Systems Inc | Neural network based automated cytological specimen classification system and method |
US5156974A (en) * | 1988-05-27 | 1992-10-20 | Biodata Corporation | Method for determining the fibrinogen level of a blood sample |
WO1991001383A1 (en) * | 1989-07-14 | 1991-02-07 | Michigan State University | Method for diagnosing blood clotting disorders |
US5051357A (en) * | 1989-07-14 | 1991-09-24 | Board Of Trustees Operating Michigan State University | Method and assay using inactivation of factors Va and VIIIa by activated Protein C to diagnose thrombic disease or assay for Protein C and kit therefor |
AU6067190A (en) * | 1989-07-20 | 1991-02-22 | Analytical Control Systems, Inc. | Improved stable coagulation controls |
SU1691741A1 (en) * | 1989-08-25 | 1991-11-15 | Донецкий государственный медицинский институт | Method for prognostication of hemorrhagic complications secondary to traumatosis |
US4998535A (en) * | 1989-09-05 | 1991-03-12 | Univ. of Washington New England Medical Center Hospitals, Inc. | Thrombolysis predictive instrument |
IE904244A1 (en) * | 1989-12-01 | 1991-06-05 | Akzo Nv | Direct fibrinogen assay |
US5169786A (en) * | 1989-12-19 | 1992-12-08 | Ortho Diagnostic Systems, Inc. | Method of determining levels of extrinsic and intrinsic clotting factors and protein c |
AU663683B2 (en) * | 1990-04-17 | 1995-10-19 | Analytical Control Systems, Inc. | Coagulation assays and reagents |
US5862304A (en) * | 1990-05-21 | 1999-01-19 | Board Of Regents, The University Of Texas System | Method for predicting the future occurrence of clinically occult or non-existent medical conditions |
US5218529A (en) * | 1990-07-30 | 1993-06-08 | University Of Georgia Research Foundation, Inc. | Neural network system and methods for analysis of organic materials and structures using spectral data |
RU2012877C1 (en) * | 1991-04-01 | 1994-05-15 | Белорусский научно-исследовательский институт охраны материнства и детства | Method of prognosticating pulmonary disorder syndrome in premature neonates |
DE69122752T2 (en) * | 1991-08-02 | 1997-04-10 | Grifols Grupo Sa | Method and device for determining the clotting time of blood and plasma |
US5716795A (en) * | 1991-10-04 | 1998-02-10 | Matschiner; John T. | Thrombomodulin-based coagulometric assay of the protein C system |
WO1993007491A1 (en) * | 1991-10-04 | 1993-04-15 | Board Of Regents Of The University Of Nebraska | A soluble thrombomodulin-based one-stage assay for vitamin k-dependent coagulation-inhibiting proteins |
JPH05180835A (en) * | 1991-12-27 | 1993-07-23 | S R L:Kk | Measurement of loops anti-coagulant |
JP2934557B2 (en) * | 1992-07-10 | 1999-08-16 | 国際試薬株式会社 | Blood coagulation time measuring method and apparatus |
JPH0636061A (en) * | 1992-07-21 | 1994-02-10 | Fujitsu Ltd | Learning system of hierarchic neural network |
US5388164A (en) * | 1992-08-19 | 1995-02-07 | Olympus Optical Co., Ltd. | Method for judging particle agglutination patterns using neural networks |
JP3224607B2 (en) * | 1992-09-11 | 2001-11-05 | 株式会社アズウェル | Method for measuring blood coagulation factor XIII activity and reagent kit for the measurement |
US5369484A (en) * | 1992-11-09 | 1994-11-29 | Akzo N.V. | Multiple discrete analyzer test apparatus and method |
RU2070327C1 (en) * | 1992-12-25 | 1996-12-10 | Безруков Александр Васильевич | Device for measurement of coagulation time |
US5344754A (en) * | 1993-01-13 | 1994-09-06 | Avocet Medical, Inc. | Assay timed by electrical resistance change and test strip |
JP2938302B2 (en) * | 1993-02-25 | 1999-08-23 | 国際試薬株式会社 | Blood coagulation time measuring method and apparatus |
JPH06249857A (en) * | 1993-02-26 | 1994-09-09 | Hitachi Ltd | Automatic analyzing device |
RU2061953C1 (en) * | 1993-03-10 | 1996-06-10 | Тюменский государственный медицинский институт МЗ РФ | Method for qualitatively determining general coagulation activity of blood platelets |
JP3476826B2 (en) * | 1993-08-16 | 2003-12-10 | アクゾ・ノベル・エヌ・ベー | Method and apparatus for automatically performing an analysis involving thrombus and hemostasis |
US5472852A (en) * | 1993-09-15 | 1995-12-05 | Oklahoma Medical Research Foundation | Assay for detection of selective Protein C inhibition by patients |
US5473732A (en) * | 1993-11-02 | 1995-12-05 | Chang; Hou-Mei H. | Relational artificial intelligence system |
US5553616A (en) * | 1993-11-30 | 1996-09-10 | Florida Institute Of Technology | Determination of concentrations of biological substances using raman spectroscopy and artificial neural network discriminator |
US5856114A (en) * | 1994-03-21 | 1999-01-05 | The University Of Vermont | Immunologic detection of factor VA fragments in hemorrhagic and thrombotic clinical settings |
WO1995030154A1 (en) * | 1994-04-28 | 1995-11-09 | Dade International Inc. | Calibrator for prothrombin time (pt) assays |
US5504011A (en) * | 1994-10-21 | 1996-04-02 | International Technidyne Corporation | Portable test apparatus and associated method of performing a blood coagulation test |
SE9403833D0 (en) * | 1994-11-08 | 1994-11-08 | Global Hemostasis Inst Mgr Ab | Analysis procedure and kit |
WO1996021740A1 (en) * | 1995-01-10 | 1996-07-18 | Hendrik Coenraad Hemker | Methods of determining endogenous thrombin potential (etp) and thrombin substrates for use in said methods |
US5708591A (en) * | 1995-02-14 | 1998-01-13 | Akzo Nobel N.V. | Method and apparatus for predicting the presence of congenital and acquired imbalances and therapeutic conditions |
US5780255A (en) * | 1995-06-09 | 1998-07-14 | Instrumentation Laboratory, S.P.A. | Protein C pathway screening test |
US5766869A (en) * | 1995-11-30 | 1998-06-16 | Ahs Hospital Corp. | Factor V ratio blood test for susceptibility to thromboembolism |
ATE281238T1 (en) * | 1996-03-22 | 2004-11-15 | Dade Behring Inc | COMBINED REAGENT COLLECTION AND TEST ARRANGEMENT |
JP3876022B2 (en) * | 1996-09-26 | 2007-01-31 | 生化学工業株式会社 | Method for quantifying substances that affect blood clotting time |
EP0841566A1 (en) * | 1996-11-11 | 1998-05-13 | Pentapharm A.G. | Positive control plasma for lupus anticoagulant |
-
1995
- 1995-06-07 US US08/477,839 patent/US5708591A/en not_active Expired - Lifetime
-
1996
- 1996-06-05 AU AU59847/96A patent/AU711998B2/en not_active Ceased
- 1996-06-05 PT PT96917181T patent/PT834145E/en unknown
- 1996-06-05 EP EP96917181A patent/EP0834145B1/en not_active Expired - Lifetime
- 1996-06-05 AT AT96917181T patent/ATE206226T1/en not_active IP Right Cessation
- 1996-06-05 ES ES96917181T patent/ES2164888T3/en not_active Expired - Lifetime
- 1996-06-05 DE DE69615556T patent/DE69615556T2/en not_active Expired - Fee Related
- 1996-06-05 JP JP50136597A patent/JP3534415B2/en not_active Expired - Fee Related
- 1996-06-05 WO PCT/US1996/008905 patent/WO1996041291A1/en active IP Right Grant
- 1996-06-05 KR KR1019970708887A patent/KR100462688B1/en not_active IP Right Cessation
- 1996-06-05 DK DK96917181T patent/DK0834145T3/en active
- 1996-06-05 CA CA002221712A patent/CA2221712C/en not_active Expired - Fee Related
-
1997
- 1997-05-21 US US08/859,773 patent/US6101449A/en not_active Expired - Fee Related
-
2000
- 2000-03-02 US US09/517,496 patent/US6269313B1/en not_active Expired - Fee Related
-
2001
- 2001-07-30 US US09/918,214 patent/US20010053959A1/en not_active Abandoned
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110016073A1 (en) * | 2008-02-14 | 2011-01-20 | Orion Diagnostica Oy | Method for predicting a future property |
US8548935B2 (en) | 2008-02-14 | 2013-10-01 | Orion Diagnostica Oy | Predicting a future property using reagents by measuring properties at points in time |
Also Published As
Publication number | Publication date |
---|---|
EP0834145B1 (en) | 2001-09-26 |
US6269313B1 (en) | 2001-07-31 |
US5708591A (en) | 1998-01-13 |
DK0834145T3 (en) | 2002-01-21 |
CA2221712C (en) | 2005-01-11 |
KR19990022407A (en) | 1999-03-25 |
KR100462688B1 (en) | 2005-06-13 |
JPH11507131A (en) | 1999-06-22 |
EP0834145A4 (en) | 2000-01-26 |
AU5984796A (en) | 1996-12-30 |
DE69615556T2 (en) | 2002-06-06 |
ES2164888T3 (en) | 2002-03-01 |
JP3534415B2 (en) | 2004-06-07 |
CA2221712A1 (en) | 1996-12-19 |
DE69615556D1 (en) | 2001-10-31 |
PT834145E (en) | 2002-03-28 |
EP0834145A1 (en) | 1998-04-08 |
US6101449A (en) | 2000-08-08 |
AU711998B2 (en) | 1999-10-28 |
ATE206226T1 (en) | 2001-10-15 |
WO1996041291A1 (en) | 1996-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0834145B1 (en) | A method and apparatus for predicting the presence of congenital and acquired imbalances and therapeutic conditions | |
US6321164B1 (en) | Method and apparatus for predicting the presence of an abnormal level of one or more proteins in the clotting cascade | |
EP1147423B1 (en) | A method and apparatus for predicting the presence of haemostatic dysfunction in a patient sample | |
US6898532B1 (en) | Method and apparatus for predicting the presence of haemostatic dysfunction in a patient sample | |
Scala et al. | Comparison of ROTEM sigma to standard laboratory tests and development of an algorithm for the management of coagulopathic bleeding in a tertiary center | |
EP1188060B1 (en) | Method and apparatus for presenting thrombosis and hemostasis assay data | |
Givens et al. | Predicting the presence of plasma heparin using neural networks to analyze coagulation screening assay optical profiles | |
EP1522860A1 (en) | A method and apparatus for predicting the presence of haemostatic dysfunction in a patient sample | |
Givens et al. | Classification of factor deficiencies from coagulation assays using neural networks | |
AU2004201011B2 (en) | Method and apparatus for presenting thrombosis and hemostasis assay data | |
RU2449281C1 (en) | Method of diagnosing pathologies of hemostasis system by means of neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: BIOMERIEUX, NORTH CAROLINA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:AKZO NOBEL N.V.;REEL/FRAME:012987/0141 Effective date: 20011221 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |