IDENTIFICATION OF PARAMETERS FOR HILL CURVE USED IN GENERAL ANESTHESIA
|
|
- Cory Russell
- 5 years ago
- Views:
Transcription
1 IDENTIFICATION OF PARAMETERS FOR HILL CURVE USED IN GENERAL ANESTHESIA Diego F. Sendoya-Losada Department of Electronic Engineering, Faculty of Engineering, Surcolombiana University, Neiva, Huila, Colombia ABSTRACT In this contribution the data collected from clinical trials in open-loop are used to identify the parameters of the nonlinear and linear Hill Curve models. Initially, the parameters are found when only the Propofol signal is considered as an input and the signal is an output - SISO case. Later, Propofol and Remifentanil signals are considered inputs and the signal is an output - MISO case. The algorithms used to find the unknown parameters of the model as well as the results are presented. Keywords: anesthesia control, hill curve identification, intensive care unit, pharmacodynamic model, pharmacokinetic model. 1. INTRODUCTION The modeling from one system s behavior and/or from physical laws aims to establish the structure of the equations which rule this system s behavior and to establish a priori value of its parameters (lengths, masses, inertia, capacities, resistance, friction and so on). But it s rarely possible to obtain a priori an accurate knowledge of all the parameters of the model. In order to sharpen and complete this knowledge, it is therefore necessary to make an identification of the system: from its reactions to given and well-known promptings, if the system can be observed, we can identify still unknown parameters. The identification, although it represents a big chapter in the automatics, can not anymore be considered as only an integral part of this discipline, with only as aim leading or diagnostic strategies. Its use in the engineering science, sometimes with a different terminology, shows it is a full discipline from the experimental sciences linked to the modeling. In automatics, two approaches for the systems identification can be distinguished: Offline and Online. In an offline mode (or batch mode), an experiment is carried out and afterwards all the data are processed simultaneously. The methods employed for offline system identification are thus based on information from the plant which has been obtained previously. In an online mode (or recursive mode), the data are used as soon as they are available. The parameters can thus be continuously estimated during the experiment. It is employed not only with control algorithms, but also for many filtering problems and signal processing. The offline mode is described in this work. The data collected from clinical trials in open-loop will be used to identify the Single Input-Single Output (SISO) and the Multiple Inputs-Single Output (MISO) Hill curve models. The algorithms used to identify the unknown parameters of the model as well as the results will be presented. 2. MATERIALS AND METHODS The relationship between the Propofol infusion rate and its effect can be described by pharmacokinetic (PK) and pharmacodynamic (PD) models. PK model describes the distribution of Propofol in the patient body and PD model describes the relationship between Propofol blood concentration and its clinical effect. PK PD V2 MUSCLE x2 EFECT SITE COMPARTMENT xe ke0 Cl2 k12 k21 DRUG INFUSION (input) INTRAVASCULAR BLOOD V1 x1 x3 Nonlinear Drug/Effect Relation Figure-1. Multi-compartmental model. In the past, different patient models have been developed for Propofol, however one popular model is the Propofol three-compartmental Schnider model (Schnider et al., 1998). This PK model can be described by a threecompartmental model (see upper part Figure-1) and it is represented in space states by the following equations: x t k + k + k k k x t [ x t] = [ k k ] [ x t] + [ ] ut x t k k x t = = [ ] [ ] k1e where denotes the amount of drug in the central compartment (blood) and its units are milligrams (mg). The peripheral compartments model the drug exchange between the blood and the others body tissues. The amount of drug in these compartments is denoted by (muscle tissue) and (fat mass), respectively. The constants, for, denote the drug transfer rate from the th tothe th compartment. The constant is the constant rate for the process that irreversibly removes Propofol from the central and peripheral compartments; its k10 Clearance: Cl1 Cl3 k13 k31 V3 FAT Real Time Monitoring (output) 27
2 units are min -1. Finally, is the Propofol infusion rate in the central compartment (blood) and represents the input of the system, its units are milligrams per second (mg/s) and it is limited to 3.33 mg/s by the pump. The output of the PK model is the plasma concentration of Propofol ( ), which is calculated dividing by, the central compartment volume. Similarly, the concentrations in the other compartments are obtained dividing and by and, respectively (Ionescu et al., 2008). In this case, the units for the compartment volume are liters (l) and these values are given by: =. =..ge = The constants are obtained as follows: = / = / = / = / = / where (clearance) is the rate at which a substance is removed from the body by the kidneys and its units are liters per minute (l/min). These values are obtained as follows: C =. +.weight.lbm +.height =. +.ge =. with (lean body mass): A =.egh egh hegh A =.egh egh hegh As can be observed in the above equations, some of the values depend on the mass (in kilogram), height (in centimeter), age (in years), and gender of the patients. Regarding the PD model (see lower part Figure- 1), an additional hypothetical effect compartment was proposed in order to represent the lag between the plasma concentration of Propofol and the patient response to the drug. The effect site compartment receives Propofol from the central compartment by a first-order model: = + where and are constants and is the amount of Propofol in the effect compartment. If the effect compartment is supposed very small compared to the other compartments, then will be a very small fraction of (Shafer et al., 1998). The apparent concentration of Propofol in the effect site compartment can be calculated once is known, because this will characterize the temporal effects of equilibration between the plasma concentration of Propofol and its corresponding effect. Thus, the equation is often used as: = [ ] where is the effect site compartment concentration of Propofol and the value of is 0.6 min -1. Therefore, the effect site compartment model can be represented as follows: C () p t Plasma concentration 1 1 s 1 k e0 Figure-2. Effect-site compartment model. In this model, / represents the time constant of the system (τ) and its value is τ =. seconds for all of the patients. The measured can be related to the effect-site concentration by the empirical static but time-varying nonlinear relationship (Bailey and Haddad, 2005), called also the Sigmoid Hill Curve: γ IS = E E ax γ γ 5 e () C t Effect concentration E is the value when the patient is awake, by convention, a value of 100 is typically assigned; E ax is the maximum effect that can be achieved by the infusion of Propofol; 5 is the Propofol concentration at half maximum effect and represents the patient sensitivity to the drug; and determines the degree of nonlinearity (steepness of the curve) (Haddad et al., 2003). Thus, the combination of the PK and PD models allows a complete modelling of the patient. The parameters involved in the PK model are known because they depend of the biometric values for each patient.regarding the PD model, the value of is also known, but not the parameters of the Sigmoid Hill Curve. The group of interest corresponds to 10 patients to whom the Propofol infusion rate was administrated by the nurse. The detailed biometric values of the patients in this group are shown in Table
3 Table-1. Biometric values of the patients. Patient Age (years) Length (cm) Weight (kg) Gender M M F M M M F M M M The average age of the group is.10 years with a standard deviation of ± 6.94 years. The average length of the patients is cm with a standard deviation of 9.88 cm. The average weight of the patients is Kg with a standard deviation of Kg. 3. RESULTS AND DISCUSSIONS 3.1Hill curve model for SISOCase In the previous section, it was established that a patient could be modeled through a multi compartmental system. This system includes pharmacodynamic and pharmacokinetic models as well as a nonlinear interaction model which describes the relationship between and the Bispectral Index. It is also known that the monitor introduces some time delay on the signal. In order to represent this situation the Sigmoid Hill Curve can be written as follows: γ τ IS = E E ax γ γ τ 5 Data of patient 10 were used in this case. In order to appreciate in a better way and to compare the real signal, it was filtered using a 3 rd order low pass filter and a cutoff frequency of 2.5 mhz to eliminate the noise. from Propofol and filtered signals are presented below. Normalized amplitude x 10 4 Figure-3. Normalized Filtered and Pr signals. If one looks at the filtered signal, several peaks can be observed, which correspond to leg movement and coughing. Since these peaks are the effect of disturbances, and not the effect of Propofol. Taking into account this, the relationship between and the Bispectral Index can be re-written as follows. γ τ IS = E E ax γ γ τ 5 + q + q Where the last part represents the disturbance model and it is a filtered positive unit step; this step occurs at time Nonlinear Hill curve In order to estimate the different parameters of the previous equation, a small part of the signals is taken. It can be observed that the difference in time between and Pr signals is around 100 seconds, i.e. the time delay in this case is 10 samples (τ = samples). 2757
4 Normalized amplitude Figure-4. Small part of filtered and C epr signals taken for identification The parameters are estimated using the MATLAB functionlsqnonlin. The results are given below Figure-5. Filtered and estimated for nonlinear SISO case. Because the number of parameters is higher than the number of equations used for estimation, the lsqnonlin function finds a local solution for this problem. The parameter values in this case are: E =. E ax =. 5 =. [μg/ml] =. X: 7790 Y: X: 7890 Y: 0.77 The algorithm needs around 3.8 seconds to estimate all the parameters. t = E E ax + γe ax γe ax C C t τ 5 + b q + a q ut t The same small part of the signals is taken and the result of the estimation is shown in Figure-6. Figure-6. Filtered and estimated for linear SISO case. In this case the number of parameters is also higher than the number of equations used for estimation. In this way, the lsqnonlinfunction finds a local solution for this problem. The parameter values found in this case are: E = E ax =. 5 =. [μg/ml] =. =. In this case, the algorithm calculates all of the parameters in 1.9 seconds Simplified linear Hill curve The linearized Hill curve model can be simplified in order to reduce the number of parameter to be estimated. IS = τ + q + q The same small part of the signals is taken and the result of the estimation is shown below Linear Hill curve Because in the Intensive Care Unit the Bispectral Index is generally between and, the previous equation is linearized around 5Pr. 2758
5 Figure-7. Filtered and estimated for simplified linear SISO case. In this case the number of parameters to be estimated is decreased but it is still higher than the number of equations used for estimation. In this way, the lsqnonlin function finds also a local solution for this problem. The parameter values found in this case are: =. =. =. Because the number of parameters is lower, the algorithm needs 1.1 seconds to find them. For all of the three models the result is good, but it can be improved if the Hill curve model takes into account the Remifentanil signal like a second input. 3.2 Hill curve model for MISOCase Pharmacodynamic drug interaction is usually described by mathematical models. The combination of Propofol and Remifentanil can be represented by a point on the graph. If the graph is straight, then the interaction is additive. If the graph bows toward the origin, the interaction is synergistic (or supra-additive), which means that the effect of the two drugs taken together is greater than the sum of their separate effect at the same doses. If the graph bows away from the origin, the interaction is infra-additive (greater amounts of both drugs are needed to produce the drug effect when administered together). Figure-8. Intravenous anesthetic synergy. Propofol and Remifentanil have a supra-additive interaction. To obtain the interaction model ( depending on the Pr and C er ), first the concentrations were normalized to their respective potencies 5Pr (Propofol effect concentration at half of the maximum effect) and 5R (Remifentanil effect concentration at half of the maximum effect). Pr = Pr [μg/ml] 5Pr [μg/ml] R = R [ng/ml] 5R [ng/ml] Remifentanil effect concentration at half of the maximum effect is calculated with the following formula (Ionescu et al., 2008): 5R =..ge [ng/ml] The ratio of the interacting drugs (regardless of the interaction type) can be expressed by: θ = Pr Pr + R θ is the concentration ratio of the new drug and ranges from 0 (Remifentanil only) to 1 (Propofol only). The new combined drug has its own sigmoid concentrationresponse relation, as depicted in the Figure-9 (general case). The concentration-response relation of the two drugs can be described as: t [U Pr t + U R t] γθ = E E ax θ [U Pr t + U R t] γθ + [U 5 θ] γθ where θ is the concentration ratio of the two drugs; Pr + R is the new drug concentration; θ 2759
6 is the steepness of the concentration-response relation at ratio θ; 5 θ is the number of units () associated with % of maximum effect at ratio θ; E ax θ is the maximum possible drug effect at ratio θ. According to (Price et al., 19), 5 θ can be expressed by the following quadratic polynomial: 5 θ = θ + θ where is a parameter which ranges between 0 and 1. Table-2 presents the possible values of for the three drug interaction types. When is equal to zero, the interaction is additive, whereas, if is positive or negative, supraand infra-additive respectively, interaction is obtained. In this case, the difference in time between and Pr signals is also around 100 seconds, i.e. the time delay is 10 samples (τ = samples). The equation due to the two drugs is rewritten as: t = E E ax [U Pr t + U R t] γ [U Pr t + U R t] γ + [ βθ + βθ ] γ + b q + a q ut t And the parameters are estimated using the MATLAB lsqnonlin function. The results are given below. Figure-9. Sigmoid concentration-response relation for Propofol and Remifentanil. Table-2. Interaction depending of. Interaction 5 θ Additive 0 1 Supra-additive >0 <1 Infra-additive <0 > Nonlinear Hill curve In order to estimate the different parameters of the concentration-response relation of the two drugs, a small part of the signals is taken. Normalized amplitude (c) Figure-10. Small part of Filtered, C epr and (c) C erem signals taken for identification. Figure-11. Filtered and estimated for nonlinear MISO case. Because the number of parameters is higher than the number of equations used for estimation, the lsqnonlin function finds a local solution for this problem. The parameter values in this case are: E =. E ax =. 5Pr =. [μg/ml] 5R =. [ng/ml] =. =. =. The algorithm needs around 5.6 seconds to estimate all of the parameters Linear Hill curve Because in the Intensive Care Unit the Bispectral Index is generally between and, the previous equation is linearized around 5Pr and 5R. 27
7 γ E ax γ γe ax t = E γ + γ + β γ + γ β γ γe ax [ γ + γ] β γ γe ax + [ γ + γ] ( β C 5Pr γ γe ax + [ γ + γ] ( β C 5R γ γe ax [ γ + β γ ] C 5Pr ) γ γe ax [ γ + β γ ] C 5R ) C Pr t τ C R t τ + b q + a q ut t The same small part of the signals is taken and the result of the estimation is shown below. t = b + m C Pr t τ + m C R t τ + b q + a q ut t The same small part of the signals is taken and the result of the estimation is shown below. Figure-12. Filtered and estimated for linear MISO case. In this case the number of parameters is also higher than the number of equations used for estimation. In this way, the lsqnonlin function finds a local solution for this problem. The parameter values found in this case are: E =. E ax =. 5Pr =. [μg/ml] 5R =. [ng/ml] =. =. =. In this case, the algorithm calculates all of the parameters in 0.4 seconds Simplified linear Hill curve The linearized Hill curve model can be simplified in order to reduce the number of parameter to be estimated. Figure-13. Filtered and estimated for simplified linear MISO case. In this case the number of parameters to be estimated is decreased but it is still higher than the number of equations used for estimation. In this way, the lsqnonlin function finds also a local solution for this problem. The parameter values found in this case are: =. =. =. =. Because the number of parameters is lower, the algorithm needs 0.4 seconds to find them. For all of the three models the result is better than when only Propofol is used, however the identification was made using a small part of the signal and only for one patient. Thus, in the future, a reliable MISO model should be identified in order to be used in the controller. The identification procedure implemented in this contribution 2761
8 is offline, but for future works an online procedure should be implemented. Before applying the Hill curve model identified using MISO case in real-life, several investigations must be carried out. Indeed, to handle the inter- and intra-patient variation of the hypnotic state, the model-based controller with adaptation on prediction model is a suitable solution. In order to apply adaptive control algorithms, it is advisable that the number of model parameters to be adapted is as low as possible. Since the initial patient model is complex, the first step is to obtain a reduced model, such as to have a minimal number of parameters. Further on, the reduced model could be used effectively in a closed-loop predictive-adaptive control strategy for Propofol and Remifentanil delivery during anesthesia. It will be thus possible to exploit the benefits of multi-drug anesthesia in an automatic control. ACKNOWLEDGEMENT This work has been funded by IWT TBM project number 0776 supported by the Flemish Institute for Innovation through Science and Technology. REFERENCES Bailey J. M., Haddad W. M Drug dosing control in clinical pharmacology. IEEE Control Syst. Mag. 25, Chilcoat R A review of the control of depth of anesthesia. Transactions of the Institute of Measurement and Control. 2(1): 38-. Cullen D., Eger E., Stevens W., Smith N., Cromwell T., Cullen B., Gregory G., Bahlman S., Dolan W., Stoelting R., Fourcade H Clinical signs of anesthesia, Anesthesiology. 36(1): Goldmann L Information-processing under general anesthesia: a review, Journal of the Royal Society of Medicine. 81: Ionescu C., De Keyser R., Torrico B. C., De Smet T., Struys M., Normey-Rico J. E Robust Predictive Control Strategy Applied for Propofol Dosing using as a Controlled Variable during Anaesthesia. IEEE Transactions on Biomedical Engineering. (3): O Hara D. A., Bogen D. K., Noordergral A The use of computers for controlling the delivery of anesthesia, Anesthesiology. 77(3): Price H. L., Kovnat P. J., Safer J. N., Conner E. H., Price M. L. 19. The uptake of thiopental by body tissues and its relation to the duration of narcosis. Clin. Pharmacol. Ther. 1: Prys-Roberts C Anaesthesia: a practical or impossible construct?, British Journal of Anaesthesia. 59: Schnider T. W., Minto C. F., Gambus P. L., Andresen C., Goodale D. B., Youngs E. J The influence of method of administration and covariates on the pharmacokinetics of propofol in adult volunteers, Anesthesiology. 88(5): Schwilden H., Stoeckel H. (eds). 1995, Control and Automation in Anaesthesia, Springer-Verlag, Berlin. Sendoya-Losada D. F., Robayo Betancourt F., Salgado Patrón J Resampling methods for signals recorded at variable sample-time intervals during general anesthesia, ARPN Journal of Engineering and Applied Sciences. 11(17): Sendoya-Losada D. F., Robayo Betancourt F., Salgado Patrón J Validation of the mathematical model for patients used in general anesthesia, ARPN Journal of Engineering and Applied Sciences. 11(19): Shafer S. L., Longnecker D. E., Tinker J. H., Morgan G. E Principles of pharmacokinetics and pharmacodynamics, In: Principles and Practice of Anesthesiology, McGraw-Hill (2nd ed.), New York,USA. Viby-Mogensen J., Ostergaard D., Donati F. D. F., Hunter J., Kampmann J., Kopman A., Proost J., Rasmussen S., Skovgaard L. F. V., Wright P Pharmacokinetic studies of neuromuscular blocking agents: good clinical research practice (GCRP), Acta Anaesthesiologica Scandmavica. 10(44): Linkens D. A., Hacisalihzade S. S Computer control systems and pharmacological drug administration: a survey, Journal of Medical Engineering & Technology. 14(2): Minto C. F., Schnider T. W., Short T. G., Gregg K. M., Gentilini A., Shafer S. L Response surface model for anesthetic drug interactions, Anesthesiology. 92(6):
A simple PK/PD model identification procedure for controller design in anesthesia
A simple PK/PD model identification procedure for controller design in anesthesia Filipa N. Nogueira, Teresa Mendonça and Paula Rocha Abstract In this paper a new estimation procedure for the parameters
More informationOn the use of fractional order PK-PD models
Journal of Physics: Conference Series PAPER OPEN ACCESS On the use of fractional order PK-PD models To cite this article: Clara Ionescu and Dana Copot 2017 J. Phys.: Conf. Ser. 783 012050 View the article
More informationMinto Model. Remifentanil
Minto Model Remifentanil This work was supported by FEDER founds through COMPETE-Operational Programme Factors of Competitiveness ( Programa Operacional Factores de Competitividade") and by Portuguese
More informationAdvanced Model Based Control Studies for the Induction and Maintenance of Intravenous Anaesthesia
> REPLAE THIS LINE WITH YOUR PAPER IDENTIFIATION NUMBER (DOUBLE-LIK HERE TO EDIT) < Advanced Model Based ontrol Studies for the Induction and Maintenance of Intravenous Anaesthesia Ioana Nașcu, Alexandra
More informationMarsh Model. Propofol
Marsh Model Propofol This work was supported by FEDER founds through COMPETE-Operational Programme Factors of Competitiveness ( Programa Operacional Factores de Competitividade") and by Portuguese founds
More informationMachine learning techniques for decision support in anesthesia
Machine learning techniques for decision support in anesthesia Olivier Caelen 1, Gianluca Bontempi 1, and Luc Barvais 2 1 Machine Learning Group, Département d Informatique, Université Libre de Bruxelles,
More informationAs mentioned in the introduction of the manuscript, isoboles are commonly used to analyze
Appendix 1: Review of Common Drug Interaction Models As mentioned in the introduction of the manuscript, isoboles are commonly used to analyze anesthetic drug interactions. Isoboles show dose combinations
More informationlaw with Probabilistic Certification for Propofol-Based Control of BIS During Anesthesia.
Feedback law with Probabilistic Certification for Propofol-Based Control of BIS During Anesthesia Mazen Alamir, Mirko Fiacchini, Isabelle Queinnec, Sophie Tarbouriech, Michel Mazerolles To cite this version:
More informationDirect Adaptive Control of Nonnegative and Compartmental Dynamical Systems with Time Delay
1 Direct Adaptive Control of Nonnegative and Compartmental Dynamical Systems with Time Delay V. Chellaboina, W. M. Haddad, J. Ramakrishnan and T. Hayakawa Mechanical and Aerospace Engineering, University
More informationThe general concept of pharmacokinetics
The general concept of pharmacokinetics Hartmut Derendorf, PhD University of Florida Pharmacokinetics the time course of drug and metabolite concentrations in the body Pharmacokinetics helps to optimize
More informationOnline Nonlinear Identification of the Effect of Drugs in Anaesthesia Using a Minimal Parameterization and BIS Measurements
2010 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 30-Jul2, 2010 ThC13.5 Online Nonlinear Identification of the Effect of Drugs in Anaesthesia Using a Minimal Parameterization
More informationEngineering Pharmacology: Pharmacokinetic Models Using Recursive Finite Difference Equations
Engineering Pharmacology: Pharmacokinetic Models Using Recursive Finite Difference Equations GLEN ATLAS Dept of Anesthesiology University of Medicine Dentistry of New Jersey Newark, NJ USA atlasgm@umdnj.edu
More informationBenchmark Problem: A PK/PD Model and Safety Constraints for Anesthesia Delivery
Benchmark Problem: A PK/PD Model and Safety Constraints for Anesthesia Delivery Victor Gan, Guy A. Dumont and Ian M. Mitchell March 31, 01 Abstract Hypnosis, or depth of unconsciousness, is one of the
More informationSystem identification and control for general anesthesia based on parsimonious Wiener models
IT Licentiate theses 2012-007 System identification and control for general anesthesia based on parsimonious Wiener models MARGARIDA MARTINS DA SILVA UPPSALA UNIVERSITY Department of Information Technology
More informationIEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 16, NO. 2, MARCH
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 16, NO. 2, MARCH 2005 387 Passivity-Based Neural Network Adaptive Output Feedback Control for Nonlinear Nonnegative Dynamical Systems Tomohisa Hayakawa, Member,
More informationPHA 4123 First Exam Fall On my honor, I have neither given nor received unauthorized aid in doing this assignment.
PHA 4123 First Exam Fall 1998 On my honor, I have neither given nor received unauthorized aid in doing this assignment. TYPED KEY Name Question 1. / pts 2. /2 pts 3. / pts 4. / pts. / pts 6. /10 pts 7.
More informationMulticompartment Pharmacokinetic Models. Objectives. Multicompartment Models. 26 July Chapter 30 1
Multicompartment Pharmacokinetic Models Objectives To draw schemes and write differential equations for multicompartment models To recognize and use integrated equations to calculate dosage regimens To
More informationAn optimal control approach to reference level tracking in general anesthesia
An optimal control approach to reference level tracking in general anesthesia J. Almeida 1 almeidajfc@gmail.com L.T. Paiva 1 ltpaiva@fe.up.pt T. Mendonça tmendo@fc.up.pt P. Rocha 1 mprocha@fe.up.pt Abstract
More informationCurrent approaches to covariate modeling
Current approaches to covariate modeling Douglas J. Eleveld, Ph.D., Assistant Professor University of Groningen, University Medical Center Groningen, Department of Anesthesiology, Groningen, The Netherlands
More informationMachine learning techniques to enable closed-loop control in anesthesia
Machine learning techniques to enable closed-loop control in anesthesia Olivier Caelen 1,2, Gianluca Bontempi 2, Eddy Coussaert 3, Luc Barvais 4 Université Libre de Bruxelles Bruxelles, Belgium {ocaelen,gbonte,coussae,lbarvais}@ulb.ac.be
More informationA Novel Screening Method Using Score Test for Efficient Covariate Selection in Population Pharmacokinetic Analysis
A Novel Screening Method Using Score Test for Efficient Covariate Selection in Population Pharmacokinetic Analysis Yixuan Zou 1, Chee M. Ng 1 1 College of Pharmacy, University of Kentucky, Lexington, KY
More informationNonlinear pharmacokinetics
5 Nonlinear pharmacokinetics 5 Introduction 33 5 Capacity-limited metabolism 35 53 Estimation of Michaelis Menten parameters(v max andk m ) 37 55 Time to reach a given fraction of steady state 56 Example:
More informationMultivariable Fuzzy Logic/Self-organizing for Anesthesia Control
Journal of Medical and Biological Engineering, 30(5): 297-306 297 Multivariable Fuzzy Logic/Self-organizing for Anesthesia Control Yuan-Chao Chou 1 Maysam F. Abbod 2,* Jiann-Shing Shieh 1 Chia-Yu Hsu 3
More informationPHARMACOKINETIC DERIVATION OF RATES AND ORDERS OF REACTIONS IN MULTI- COMPARTMENT MODEL USING MATLAB
IJPSR (2016), Vol. 7, Issue 11 (Research Article) Received on 29 May, 2016; received in revised form, 07 July, 2016; accepted, 27 July, 2016; published 01 November, 2016 PHARMACOKINETIC DERIVATION OF RATES
More informationTotal Mass TCI driven by Parametric Estimation
Total Mass TCI driven by Parametric Estimation Margarida M. Silva, Cláudia Sousa, Raquel Sebastião, João Gama, Teresa Mendonça, Paula Rocha and Simão Esteves DMA - Faculdade de Ciências da Universidade
More informationThe Modi ed Chain Method for Systems of Delay Di erential Equations with Applications
The Modi ed Chain Method for Systems of Delay Di erential Equations with Applications Theses of a PhD Dissertation Beáta Krasznai Supervisor: Professor Mihály Pituk University of Pannonia Faculty of Information
More informationHomework 1 (PHA 5127)
Homework 1 (PHA 5127) 1. The elimination rate constant of a drug is 1 hr -1 : k e =1 hr -1 A. Half-life: t 1/2 = ln(2)/k e = 0.693/1 hr -1 = 0.693 hr B. C 1 =5ng/ml First-order elimination: k e = (ln(c
More informationOn drug transport after intravenous administration
On drug transport after intravenous administration S.Piekarski (1), M.Rewekant (2) Institute of Fundamental Technological Research Polish Academy of Sciences (1), Medical University of Warsaw, Poland Abstract
More informationNoncompartmental vs. Compartmental Approaches to Pharmacokinetic Data Analysis Paolo Vicini, Ph.D. Pfizer Global Research and Development David M.
Noncompartmental vs. Compartmental Approaches to Pharmacokinetic Data Analysis Paolo Vicini, Ph.D. Pfizer Global Research and Development David M. Foster., Ph.D. University of Washington October 28, 2010
More informationNoncompartmental vs. Compartmental Approaches to Pharmacokinetic Data Analysis Paolo Vicini, Ph.D. Pfizer Global Research and Development David M.
Noncompartmental vs. Compartmental Approaches to Pharmacokinetic Data Analysis Paolo Vicini, Ph.D. Pfizer Global Research and Development David M. Foster., Ph.D. University of Washington October 18, 2012
More informationOptimal drug dosing control for intensive care unit sedation by using a hybrid deterministic stochastic pharmacokinetic and pharmacodynamic model
OPTIMAL CONTROL APPLICATIONS AND METHODS Optim. Control Appl. Meth. 2013; 34:547 561 Published online 28 June 2012 in Wiley Online Library (wileyonlinelibrary.com)..2038 Optimal drug dosing control for
More informationBeka 2 Cpt: Two Compartment Model - Loading Dose And Maintenance Infusion
MEDSCI 719 Pharmacometrics Beka 2 Cpt: Two Compartment Model - Loading Dose And Maintenance Infusion Objective 1. To observe the time course of drug concentration in the central and peripheral compartments
More informationOn approximate solutions in pharmacokinetics
On separation of time scales in pharmacokinetics Piekarski S, Rewekant M. IPPT PAN, WUM Abstract A lot of criticism against the standard formulation of pharmacokinetics has been raised by several authors.
More informationLove and differential equations: An introduction to modeling
Nonlinear Models 1 Love and differential equations: An introduction to modeling November 4, 1999 David M. Allen Professor Emeritus Department of Statistics University of Kentucky Nonlinear Models Introduction
More informationNonlinear controller for bispectral index tracking: Robustness and on-line retuning
CONROL ENGINEERING PRACICE Vol. 58 p. 343-353, 207 DOI: 0.06/j.conengprac.206.03.02 Nonlinear controller for bispectral index tracking: Robustness and on-line retuning Filipa N. Nogueira a,b, eresa Mendonça
More informationEstimating terminal half life by non-compartmental methods with some data below the limit of quantification
Paper SP08 Estimating terminal half life by non-compartmental methods with some data below the limit of quantification Jochen Müller-Cohrs, CSL Behring, Marburg, Germany ABSTRACT In pharmacokinetic studies
More information3 Results. Part I. 3.1 Base/primary model
3 Results Part I 3.1 Base/primary model For the development of the base/primary population model the development dataset (for data details see Table 5 and sections 2.1 and 2.2), which included 1256 serum
More informationChapter One. Introduction
Chapter One Introduction With the ever-increasing influence of mathematical modeling and engineering on biological, social, and medical sciences, it is not surprising that dynamical system theory has played
More informationMODELING THE ABSORPTION, CIRCULATION, AND METABOLISM OF TIRAPAZAMINE
MODELING THE ABSORPTION, CIRCULATION, AND METABOLISM OF TIRAPAZAMINE Olivia Nguyen and Xiaoyu Shi Problem Statement: Tirapazamine (TPZ; 1,2,4-benzotriazin-3-amine 1,4-dioxide) is a hypoxia-activated prodrug
More informationPractice of SAS Logistic Regression on Binary Pharmacodynamic Data Problems and Solutions. Alan J Xiao, Cognigen Corporation, Buffalo NY
Practice of SAS Logistic Regression on Binary Pharmacodynamic Data Problems and Solutions Alan J Xiao, Cognigen Corporation, Buffalo NY ABSTRACT Logistic regression has been widely applied to population
More informationIntegration of SAS and NONMEM for Automation of Population Pharmacokinetic/Pharmacodynamic Modeling on UNIX systems
Integration of SAS and NONMEM for Automation of Population Pharmacokinetic/Pharmacodynamic Modeling on UNIX systems Alan J Xiao, Cognigen Corporation, Buffalo NY Jill B Fiedler-Kelly, Cognigen Corporation,
More informationRADIOLOGIV TECHNOLOGY 4912 COMPREHENSEIVE REVIEW/MRI WORSHEET #1- PATIENT CARE AND SAFETY/PHYSICAL PRINCIPLES
RADIOLOGIV TECHNOLOGY 4912 COMPREHENSEIVE REVIEW/MRI WORSHEET #1- PATIENT CARE AND SAFETY/PHYSICAL PRINCIPLES 1. What are potential consequences to patients and personnel should there be a release of gaseous
More informationAn Iterative Algorithm for the Subspace Identification of SISO Hammerstein Systems
An Iterative Algorithm for the Subspace Identification of SISO Hammerstein Systems Kian Jalaleddini R. E. Kearney Department of Biomedical Engineering, McGill University, 3775 University, Montréal, Québec
More informationFitting PK Models with SAS NLMIXED Procedure Halimu Haridona, PPD Inc., Beijing
PharmaSUG China 1 st Conference, 2012 Fitting PK Models with SAS NLMIXED Procedure Halimu Haridona, PPD Inc., Beijing ABSTRACT Pharmacokinetic (PK) models are important for new drug development. Statistical
More informationPET Tracer Kinetic Modeling In Drug
PET Tracer Kinetic Modeling In Drug Discovery Research Applications Sandra M. Sanabria-Bohórquez Imaging Merck & Co., Inc. Positron Emission Tomography - PET PET is an advanced d imaging i technique permitting
More informationMIDW 125 Math Review and Equation Sheet
MIDW 125 Math Review and Equation Sheet 1. The Metric System Measures of weight are based on the gram (g): 1 kilogram = 1 kg = 1000 gram = 1000 g = 10 3 g 1 milligram = 1 mg = 10 3 g 1 microgram = 1 g
More informationThe Calculus Behind Generic Drug Equivalence
The Calculus Behind Generic Drug Equivalence Stanley R. Huddy and Michael A. Jones Stanley R. Huddy srh@fdu.edu, MRID183534, ORCID -3-33-231) isanassistantprofessorat Fairleigh Dickinson University in
More informationRELAY CONTROL WITH PARALLEL COMPENSATOR FOR NONMINIMUM PHASE PLANTS. Ryszard Gessing
RELAY CONTROL WITH PARALLEL COMPENSATOR FOR NONMINIMUM PHASE PLANTS Ryszard Gessing Politechnika Śl aska Instytut Automatyki, ul. Akademicka 16, 44-101 Gliwice, Poland, fax: +4832 372127, email: gessing@ia.gliwice.edu.pl
More informationNonlinear identification of a minimal
Nonlinear identification of a minimal 1 NeuroMuscular Blockade model in anaesthesia Margarida Martins da Silva, Torbjörn Wigren, Senior Member, IEEE, and Teresa Mendonça, Member, IEEE Abstract The paper
More informationIntroduction to Pharmacokinetic Modelling Rationale
Introduction to Pharmacokinetic Modelling Rationale Michael Weiss michael.weiss@medizin.uni-halle.de Martin Luther University Halle-Wittenberg Philippus von Hohenheim, known as Paracelsus (1493-1541) Everything
More informationDESIGN OF AN ON-LINE TITRATOR FOR NONLINEAR ph CONTROL
DESIGN OF AN ON-LINE TITRATOR FOR NONLINEAR CONTROL Alex D. Kalafatis Liuping Wang William R. Cluett AspenTech, Toronto, Canada School of Electrical & Computer Engineering, RMIT University, Melbourne,
More informationDescription of UseCase models in MDL
Description of UseCase models in MDL A set of UseCases (UCs) has been prepared to illustrate how different modelling features can be implemented in the Model Description Language or MDL. These UseCases
More informationAnalysis of Allosterism in Functional Assays a
JPET This Fast article Forward. has not been Published copyedited and on formatted. July 26, The 2005 final as version DOI:10.1124/jpet.105.090886 may differ from this version. Analysis of Allosterism
More informationPredicting the Onset of Nonlinear Pharmacokinetics
Citation: CPT Pharmacometrics Syst Pharmacol (208) 7, 670 677; doi:0002/psp4236 ARTICLE Predicting the Onset of Nonlinear Pharmacokinetics Andrew M Stein and Lambertus A Peletier 2 When analyzing the pharmacokinetics
More informationHow Critical is the Duration of the Sampling Scheme for the Determination of Half-Life, Characterization of Exposure and Assessment of Bioequivalence?
How Critical is the Duration of the Sampling Scheme for the Determination of Half-Life, Characterization of Exposure and Assessment of Bioequivalence? Philippe Colucci 1,2 ; Jacques Turgeon 1,3 ; and Murray
More informationPubH 7405: REGRESSION ANALYSIS INTRODUCTION TO LOGISTIC REGRESSION
PubH 745: REGRESSION ANALYSIS INTRODUCTION TO LOGISTIC REGRESSION Let Y be the Dependent Variable Y taking on values and, and: π Pr(Y) Y is said to have the Bernouilli distribution (Binomial with n ).
More informationCorrection of the likelihood function as an alternative for imputing missing covariates. Wojciech Krzyzanski and An Vermeulen PAGE 2017 Budapest
Correction of the likelihood function as an alternative for imputing missing covariates Wojciech Krzyzanski and An Vermeulen PAGE 2017 Budapest 1 Covariates in Population PKPD Analysis Covariates are defined
More informationUsing ODEs to Model Drug Concentrations within the Field of Pharmacokinetics
Augustana College Augustana Digital Commons Mathematics: Student Scholarship & Creative Works Mathematics Spring 2016 Using ODEs to Model Drug Concentrations within the Field of Pharmacokinetics Andrea
More informationNonlinear Mixed Effects Models
Nonlinear Mixed Effects Modeling Department of Mathematics Center for Research in Scientific Computation Center for Quantitative Sciences in Biomedicine North Carolina State University July 31, 216 Introduction
More informationFixed-Complexity Piecewise Ellipsoidal Representation of the Continual Reachability Set Based on Ellipsoidal Techniques
Fixed-Complexity Piecewise Ellipsoidal Representation of the Continual Reachability Set Based on Ellipsoidal Techniques Shahab Kaynama, Meeko Oishi, Ian M. Mitchell, and Guy A. Dumont Abstract In a previous
More informationConfidence and Prediction Intervals for Pharmacometric Models
Citation: CPT Pharmacometrics Syst. Pharmacol. (218) 7, 36 373; VC 218 ASCPT All rights reserved doi:1.12/psp4.12286 TUTORIAL Confidence and Prediction Intervals for Pharmacometric Models Anne K ummel
More informationStudy of different spectral parameters to measure depth of anesthesia.
Biomedical Research 2017; 28 (15): 6533-6540 ISSN 0970-938X www.biomedres.info Study of different spectral parameters to measure depth of anesthesia. Elizabeth George *, Glan Devadas G, Avinashe KK, Sudharsana
More informationJUSTIFICATION OF INPUT AND OUTPUT CONSTRAINTS INCORPORATION INTO PREDICTIVE CONTROL DESIGN
JUSTIFICATION OF INPUT AND OUTPUT CONSTRAINTS INCORPORATION INTO PREDICTIVE CONTROL DESIGN J. Škultéty, E. Miklovičová, M. Mrosko Slovak University of Technology, Faculty of Electrical Engineering and
More informationPerformance assessment of MIMO systems under partial information
Performance assessment of MIMO systems under partial information H Xia P Majecki A Ordys M Grimble Abstract Minimum variance (MV) can characterize the most fundamental performance limitation of a system,
More informationTK Solver Case Study: Pharmacokinetics
TK Solver Case Study: Pharmacokinetics The content for this example was provided by Professor Prasad Dhurjati of the University of Delaware. It is used there in the class Applied Mathematics for Chemical
More informationONE-COMPARTMENT KINETICS
British Journal of Anaesthesia 1992; 69: 387-396 REVIEW ARTICLE ONE-COMPARTMENT KINETICS J. NORMAN Pharmacokinetics may appear daunting, especially if the background mathematical development is presented
More informationNonlinear Modeling and Feedback Control of Drug Delivery in Anesthesia
Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 1187 Nonlinear Modeling and Feedback Control of Drug Delivery in Anesthesia MARGARIDA M. SILVA ACTA UNIVERSITATIS
More informationNonlinear PD Controllers with Gravity Compensation for Robot Manipulators
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 4, No Sofia 04 Print ISSN: 3-970; Online ISSN: 34-408 DOI: 0.478/cait-04-00 Nonlinear PD Controllers with Gravity Compensation
More informationPrinciples of Covariate Modelling
1 Principles of Covariate Modelling Nick Holford Dept of Pharmacology & Clinical Pharmacology University of Auckland New Zealand 2 What is a Covariate? A covariate is any variable that is specific to an
More informationPlant and animal cells (eukaryotic cells) have a cell membrane, cytoplasm and genetic material enclosed in a nucleus.
4.1 Cell biology Cells are the basic unit of all forms of life. In this section we explore how structural differences between types of cells enables them to perform specific functions within the organism.
More informationLecture 5 Classical Control Overview III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 5 Classical Control Overview III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore A Fundamental Problem in Control Systems Poles of open
More informationIndex. INDEX_p /15/02 3:08 PM Page 765
INDEX_p.765-770 11/15/02 3:08 PM Page 765 Index N A Adaptive control, 144 Adiabatic reactors, 465 Algorithm, control, 5 All-pass factorization, 257 All-pass, frequency response, 225 Amplitude, 216 Amplitude
More informationOptimization of Model-Reference Variable-Structure Controller Parameters for Direct-Current Motor
Journal of Computations & Modelling, vol., no.3, 1, 67-88 ISSN: 179-765 (print), 179-885 (online) Scienpress Ltd, 1 Optimization of Model-Reference Variable-Structure Controller Parameters for Direct-Current
More informationMA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems
MA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems Integrating Mathematical and Statistical Models Recap of mathematical models Models and data Statistical models and sources of
More informationBody Size is the most important quantitative determinant of drug dose
1 2 Babies, Children, Adults Are All One Species Integrated Rational Dosing Nick Holford Dept Pharmacology & Clinical Pharmacology University of Auckland Body Size is the most important quantitative determinant
More informationF. Combes (1,2,3) S. Retout (2), N. Frey (2) and F. Mentré (1) PODE 2012
Prediction of shrinkage of individual parameters using the Bayesian information matrix in nonlinear mixed-effect models with application in pharmacokinetics F. Combes (1,2,3) S. Retout (2), N. Frey (2)
More informationANN Control of Non-Linear and Unstable System and its Implementation on Inverted Pendulum
Research Article International Journal of Current Engineering and Technology E-ISSN 2277 4106, P-ISSN 2347-5161 2014 INPRESSCO, All Rights Reserved Available at http://inpressco.com/category/ijcet ANN
More informationModelling a complex input process in a population pharmacokinetic analysis: example of mavoglurant oral absorption in healthy subjects
Modelling a complex input process in a population pharmacokinetic analysis: example of mavoglurant oral absorption in healthy subjects Thierry Wendling Manchester Pharmacy School Novartis Institutes for
More informationReducing System Response Time and Noise of Electrochemical Gas Sensors - Discussed for Propofol Monitoring in Breathing Gas
Reducing System Response Time and Noise of Electrochemical Gas Sensors - Discussed for Propofol Monitoring in Breathing Gas Dammon Ziaian, Philipp Rostalski, Andreas Hengstenberg and Stefan Zimmermann
More informationNonparametric Bayes approach for a semimechanistic pharmacodynamic model
University of New Mexico UNM Digital Repository Mathematics & Statistics ETDs Electronic Theses and Dissertations 6-23-2015 Nonparametric Bayes approach for a semimechanistic pharmacokinetic and pharmacodynamic
More informationSelf-Tuning Control for Synchronous Machine Stabilization
http://dx.doi.org/.5755/j.eee.2.4.2773 ELEKTRONIKA IR ELEKTROTECHNIKA, ISSN 392-25, VOL. 2, NO. 4, 25 Self-Tuning Control for Synchronous Machine Stabilization Jozef Ritonja Faculty of Electrical Engineering
More informationAPPLICATION OF LAPLACE TRANSFORMS TO A PHARMACOKINETIC OPEN TWO-COMPARTMENT BODY MODEL
Acta Poloniae Pharmaceutica ñ Drug Research, Vol. 74 No. 5 pp. 1527ñ1531, 2017 ISSN 0001-6837 Polish Pharmaceutical Society APPLICATION OF LAPLACE TRANSFORMS TO A PHARMACOKINETIC OPEN TWO-COMPARTMENT BODY
More informationCompartmental Systems - definition. Compartmental models: system analysis and control issues. Compartmental Systems - general model
Compartmental Systems - definition Compartmental models: system analysis and control issues Compartmental system Finite number of compartments which exchange material with each other in such a way that
More informationModeling biological systems - The Pharmaco-Kinetic-Dynamic paradigm
Modeling biological systems - The Pharmaco-Kinetic-Dynamic paradigm Main features:. Time profiles and complex systems 2. Disturbed and sparse measurements 3. Mixed variabilities One of the most highly
More informationA STUDY ON THE STATE ESTIMATION OF NONLINEAR ELECTRIC CIRCUITS BY UNSCENTED KALMAN FILTER
A STUDY ON THE STATE ESTIMATION OF NONLINEAR ELECTRIC CIRCUITS BY UNSCENTED KALMAN FILTER Esra SAATCI Aydın AKAN 2 e-mail: esra.saatci@iku.edu.tr e-mail: akan@istanbul.edu.tr Department of Electronic Eng.,
More informationAdaptive Predictive Observer Design for Class of Uncertain Nonlinear Systems with Bounded Disturbance
International Journal of Control Science and Engineering 2018, 8(2): 31-35 DOI: 10.5923/j.control.20180802.01 Adaptive Predictive Observer Design for Class of Saeed Kashefi *, Majid Hajatipor Faculty of
More informationIdentification in closed-loop, MISO identification, practical issues of identification
Identification in closed-loop, MISO identification, practical issues of identification CHEM-E7145 Advanced Process Control Methods Lecture 4 Contents Identification in practice Identification in closed-loop
More informationThM06-2. Coprime Factor Based Closed-Loop Model Validation Applied to a Flexible Structure
Proceedings of the 42nd IEEE Conference on Decision and Control Maui, Hawaii USA, December 2003 ThM06-2 Coprime Factor Based Closed-Loop Model Validation Applied to a Flexible Structure Marianne Crowder
More informationContributions to the Automatic Control of Drug Administration for Anesthesia
Contributions to the Automatic Control of Drug Administration for Anesthesia Juliana Filipa Carvalho de Almeida Supervisor: Maria Paula Macedo Rocha Malonek (FEUP) Co-Supervisor: Teresa Maria de Gouveia
More informationProcess Identification for an SOPDT Model Using Rectangular Pulse Input
Korean J. Chem. Eng., 18(5), 586-592 (2001) SHORT COMMUNICATION Process Identification for an SOPDT Model Using Rectangular Pulse Input Don Jang, Young Han Kim* and Kyu Suk Hwang Dept. of Chem. Eng., Pusan
More informationRegression techniques provide statistical analysis of relationships. Research designs may be classified as experimental or observational; regression
LOGISTIC REGRESSION Regression techniques provide statistical analysis of relationships. Research designs may be classified as eperimental or observational; regression analyses are applicable to both types.
More informationFurther Results on Model Structure Validation for Closed Loop System Identification
Advances in Wireless Communications and etworks 7; 3(5: 57-66 http://www.sciencepublishinggroup.com/j/awcn doi:.648/j.awcn.735. Further esults on Model Structure Validation for Closed Loop System Identification
More informationState Observers and the Kalman filter
Modelling and Control of Dynamic Systems State Observers and the Kalman filter Prof. Oreste S. Bursi University of Trento Page 1 Feedback System State variable feedback system: Control feedback law:u =
More information10-2: Exponential Function Introduction
Math 95 10-2: Exponential Function Introduction Quadratic functions, like y = x B, are made of a sum of powers of the independent variable. In a quadratic function, the variables are in the base of the
More informationUsing the Kalman Filter to Estimate the State of a Maneuvering Aircraft
1 Using the Kalman Filter to Estimate the State of a Maneuvering Aircraft K. Meier and A. Desai Abstract Using sensors that only measure the bearing angle and range of an aircraft, a Kalman filter is implemented
More informationThe Art of Mathematical Modeling
The Art of Mathematical Modeling Lisa Rogers (1) The Art of Mathematical Modeling For these exercises, please also refer to the Algorithm for Constructing a Mathematical Model.
More informationRecent Advances in Positive Systems: The Servomechanism Problem
Recent Advances in Positive Systems: The Servomechanism Problem 47 th IEEE Conference on Decision and Control December 28. Bartek Roszak and Edward J. Davison Systems Control Group, University of Toronto
More informationTrue False. Question The memory capacity of a flash drive is measured in gigabytes so that the capacity can be expressed using simple integers.
1 of 8 TEST BANK > CONTRO PANE > POO MANAGER > POO CANVAS Pool Canvas Add, modify, and remove questions. Select a question type from the Add Question drop-down list and click Go to add questions. Use Creation
More informationON SEPARATION PRINCIPLE FOR THE DISTRIBUTED ESTIMATION AND CONTROL OF FORMATION FLYING SPACECRAFT
ON SEPARATION PRINCIPLE FOR THE DISTRIBUTED ESTIMATION AND CONTROL OF FORMATION FLYING SPACECRAFT Amir Rahmani (), Olivia Ching (2), and Luis A Rodriguez (3) ()(2)(3) University of Miami, Coral Gables,
More informationNumerical Simulations of Stochastic Circulatory Models
Bulletin of Mathematical Biology (1999) 61, 365 377 Article No. bulm.1998.94 Available online at http://www.idealibrary.com on Numerical Simulations of Stochastic Circulatory Models Gejza Wimmer Mathematical
More information