This article has Open Peer Review reports available.

# Integration of modeling and simulation into hospital-based decision support systems guiding pediatric pharmacotherapy

- Jeffrey S Barrett
^{1}Email author, - John T Mondick
^{1}, - Mahesh Narayan
^{1}, - Kalpana Vijayakumar
^{2}and - Sundararajan Vijayakumar
^{2}

**8**:6

https://doi.org/10.1186/1472-6947-8-6

© Barrett et al; licensee BioMed Central Ltd. 2008

**Received: **30 July 2007

**Accepted: **28 January 2008

**Published: **28 January 2008

## Abstract

### Background

Decision analysis in hospital-based settings is becoming more common place. The application of modeling and simulation approaches has likewise become more prevalent in order to support decision analytics. With respect to clinical decision making at the level of the patient, modeling and simulation approaches have been used to study and forecast treatment options, examine and rate caregiver performance and assign resources (staffing, beds, patient throughput). There us a great need to facilitate pharmacotherapeutic decision making in pediatrics given the often limited data available to guide dosing and manage patient response. We have employed nonlinear mixed effect models and Bayesian forecasting algorithms coupled with data summary and visualization tools to create drug-specific decision support systems that utilize individualized patient data from our electronic medical records systems.

### Methods

Pharmacokinetic and pharmacodynamic nonlinear mixed-effect models of specific drugs are generated based on historical data in relevant pediatric populations or from adults when no pediatric data is available. These models are re-executed with individual patient data allowing for patient-specific guidance via a Bayesian forecasting approach. The models are called and executed in an interactive manner through our web-based dashboard environment which interfaces to the hospital's electronic medical records system.

### Results

The methotrexate dashboard utilizes a two-compartment, population-based, PK mixed-effect model to project patient response to specific dosing events. Projected plasma concentrations are viewable against protocol-specific nomograms to provide dosing guidance for potential rescue therapy with leucovorin. These data are also viewable against common biomarkers used to assess patient safety (e.g., vital signs and plasma creatinine levels). As additional data become available via therapeutic drug monitoring, the model is re-executed and projections are revised.

### Conclusion

The management of pediatric pharmacotherapy can be greatly enhanced via the immediate feedback provided by decision analytics which incorporate the current, best-available knowledge pertaining to dose-exposure and exposure-response relationships, especially for narrow therapeutic agents that are difficult to manage.

## Background

Decision making in a hospital environment occurs at multiple levels of the organization and in a variety of departmental settings. Likewise, numerous stakeholders including hospital administration, staffing planners, facilities management, pharmacy administration, caregivers and physicians as well as healthcare providers are reliant on decision support systems (DSS) to facilitate decision making in their specific areas. Ultimately, the outcome, at least in theory, should be better decisions yielding more efficient provision of services and optimal (most appropriate and cost-effective) patient care. Informing today's decision makers is a cadre of tools and decision analytics. Historically, hospital environments have not been the hallmark of innovation in decision analytics and the discrepancy between the hospital environment and other industries has received a great deal of attention recently. In late 2005, the National Academies of Engineering and Institute of Medicine issued a joint report that cited the urgency and importance of bringing contemporary System Engineering techniques to healthcare. Annual gross waste of a staggering 30–40% of every dollar spent in healthcare and continued medical errors that cause nearly 100,000 patient deaths and serious injuries yearly were only two of the more serious problems covered in the report [1, 2]. This situation is changing in a dramatic manner as information technology, engineering, clinical and informatics scientists collaborate on analytic approaches that address decision requirements of the current inpatient environment.

Diversity in modeling and simulation applications in the hospital setting

M&S Fields | Application | Approach |
---|---|---|

Treatment Outcomes | • Medical folder management system – physician clinical decision making [3] • Cancer pharmacotherapy multi-drug decision support [4] • Predictive model to predict Clostridium difficile infection (diahhrea) outbreaks [5] • Hospital-wide surveillance for nosocomial infection to assess patient risk [6] • Methicillin-resistant Staphylococcus aureus transmission among hospitalized patients – risk factors and prediction [7] | • DSS interfaced to EMRS • KITT model and decision tree • Reversible jump MCMC model • Logistic regression model • Monte Carlo simulation |

Healthcare Costs | • Health care costs of geriatric inpatients [8] • Hospital-acquired infection costs [9] • Costs and outcomes of cardiovascular surgery [10] | • Bayesian Network Theory/Model • Monte Carlo simulation model • Systems dynamic model STELLA |

Patient Flow/Occupancy | • Patient flow in a pediatric emergency department [11] • Critical care planning capacity [12] • Healthcare facility patient flow [13] • Hospital patient flow [13] | • Discrete event simulation • CART analysis • Queuing network system • Clustering |

Hospital Operations | • Hospital operations for emergency response [15] • Length of stay in the ICU [16] • Directly observed therapy in newly diagnosed HIV infection [17] • System-level investigation of emergency department (ED) operations [18] • Healthcare quality improvement via simulation [19] • Optimum operating room staffing needs for trauma centers [20] • ICU duration/length of stay analysis [21] | • Transient modeling regression approach • Linear regression • Probabilistic Markov Model • Discrete event simulation (EDSIM) • Multivariate simulation models • Queuing and simulation methods • Class probability tree |

Pharmacotherapy is generally concerned with the safe and effective management of drug administration. It implies an understanding of drug pharmacokinetics (PK) and pharmacodynamics (PD) so that individual dosing guidance, when necessary, can be provided to optimize patient response within their individual therapeutic window. Pediatric pharmacotherapy can be challenging due to developmental changes that may alter drug kinetics, pathophysiologic differences that may alter pharmacodynamics, disease etiologies that may be different from adults, and other factors that may result in great variation in safety and efficacy outcomes. The situation becomes more convoluted when one considers children and the paucity of well-controlled pediatric clinical trials. This situation, despite the efforts of the Food and Drug Administration and the US Congress, is not likely to improve substantially due to the economic reality of the pediatric market.

We have begun to interface our electronic medical record system with decision support analytics (drug dashboards) that summarize individual patient records and assemble the most relevant clinical data associated with drug therapy. Data visualization tools summarize patient profiles of lab values, vital signs, and associated biomarkers into tables and plots based on user-defined requirements. More importantly, this data populates models that predict future events as drug therapy continues. A prototype methotrexate dashboard is demonstrated. A population-based pharmacokinetic model (nonlinear mixed effect model) is used to simulate individual patient MTX drug concentrations based on that patient's current dosing regimen and compares expected exposures with nomograms that predict toxicity. The compilation of additional dashboards is planned for the construction of a larger pediatric knowledgebase (PKB) currently under construction. Our objective is to define the general approach of pharmacostatistical model building while demonstrating how such models can be interfaced to electronic medical record data front-ended by a web-based decision support system.

## Methods

### Population Pharmacokinetics

Population pharmacokinetics is the study of the sources and correlates of variability in drug concentrations among individuals who represent the target patient population receiving clinically relevant doses of a drug of interest. The goal of population pharmacokinetic analysis to identify pathophysiologic factors that cause changes in the dose-concentration relationship and the extent of these changes so that, if such changes are associated with clinically significant shifts in the therapeutic index, dosage can be appropriately modified. Population pharmacokinetic models adhere into a hierarchical structure. At the initial stage, the relationship between concentration and time (pharmacokinetics) is modeled for an individual patient. At the second stage, pharmacokinetic parameters that define each of the individual patients' drug concentration profiles are assigned some distributional form, after accounting for relevant covariate information. A primary aim of a population analysis is to determine covariates that are important predictors of pharmacokinetic parameters. A Bayesian model then requires a third stage in which prior distributions are specified for the parameters defining the second-stage distributional form and the intra-individual variance parameters. Such a Bayesian model defines and estimates the variability that is observed both in individual concentrations and between different individuals' pharmacokinetic parameters. This model framework makes it is possible to determine appropriate patient-specific dosage regimens that ensure the attainment of desirable drug concentrations.

*likelihood*or

*log*-likelihood function (or to minimize the negative log-likelihood function). In general terms, the likelihood function is defined as:

The probability (now called *L*, the *likelihood*) is predicted in the sample data, given the respective regression model. Provided that all observations are independent of each other, this likelihood is the geometric sum (Π, across *i* = 1 to *n* cases) of probabilities for each individual observation (*i*) to occur, given the respective model and parameters (*θ*'s) for the *x* values. As it is customary to express this function as a natural logarithm, the geometric sum becomes a regular arithmetic sum (Σ, across *i* = 1 to *n* cases). The larger the likelihood of the model, the larger is the probability of the dependent variable values to occur in the sample and the better is the fit of the model to the data. If all assumptions for standard multiple regression are met, then the standard least squares estimation method will yield results identical to the maximum likelihood method. If the assumption of equal error variances across the range of the *x* variable(s) is violated, then the weighted least squares method will yield maximum likelihood estimates.

The typical structural model is chosen from one of several compartmental models which incorporate the route of administration as a fixed input into the model with certain assumptions (i.e., linear or zero order input). Compartmental models are, for the most part, empirical even though they may incorporate some mechanistic assumptions so they appear more realistic. Numerically, they are generally easier to handle as opposed to mechanistic models. Complex mechanistic and/or highly parameterized structural models can be accommodated as well of course. The prediction engine discussed herein is not limited by the nature of the model definition.

The framework for the mixed effect modeling approach to population pharmacokinetic analysis can be defined as follows: for i = 1, ... n individuals in a population of interest, let x_{ij}, j = 1, ... n_{j} represent the design points on which the y_{ij} responses are observed. In the pharmacokinetic (PK) setting, x_{ij} are typically the sampling time points and y_{ij} are the observed concentrations in the biologic matrix of interest (usually plasma or blood). Hence, the PK response can be described by

y_{ij} = *f*(*θ*
_{i}, x_{ij}) + *ε*
_{ij}

where the function *f* denotes the structural model; *θ*
_{i} is the p × 1 parameter vector for the i^{th} individual and *ε*
_{ij} are the independently and identically distributed (*i.i.d*.) error terms assumed to be normal random variables with a zero mean and a variance (*σ*
_{e}
^{2}) which may depend on the mean concentration. The *ε*
_{ij}'s account for the intraindividual variability and may incorporate model misspecification or other unresolved (or incorrect) error partitioning. In most population pharmacokinetic software the structural model is chosen from a library of compartmental models, expressed as a closed form system of equations or defined via differential equations. The probability density function which accounts for the within-individual variability only as

p(y_{ij} | *θ*
_{i}, x_{ij})

The intra-individual variation about the ith individual is defined when the distribution *ε*
_{i} of is specified. The second stage model defines the between-individual variability in the parameters as follows

*θ*
_{i} = *θ* + *η*
_{i}

where *θ* is the mean parameter vector for the population and *η*
_{i} are the individual deviations assumed to be *i.i.d*. and normal with zero mean vector and covariance matrix Ω. The expression of *θ*
_{i} shown (additive) is one of numerous ways that individual *θ*'s can be defined. In addition, the population *θ* can be expressed as a function of covariates (β_{i}). The covariate matrix Ω captures both the variance and covariance among the *η*'s. The density of the second stage model can then be defined as

p(*θ*
_{i} | *θ*, Ω, β_{i})

where β_{i} represents the individual patient covariate data (i.e., age, sex, race, etc.). The third stage of the mixed effect model approach would represent a Bayesian representation in which the model would contain the prior distributions of the population parameters as mentioned previously.

### Prediction Models for Toxicity and Adverse Events

*i*at sampling time

*j*is given by

*p*

_{ ij }. The ratio of the probability of that event occurring vs. the probability that it does not occur is given by the Odds Ratio; the log of the Odds Ratio is known as the logit function (

*λ*

_{ ij }):

Independent variables and covariates (*x*
_{
ij
}) will be incorporated into the model via the logit function, with population typical population (*θ*) and individual random effect (*η*
_{
i
}) parameters to be estimated:

*λ*
_{
ij
}= *f*(*x*
_{
ij
}, *θ*, *η*
_{
i
})

Covariate effects, and random effects can modulate the predicted probability in a positive or negative direction, with the probability constrained between the values of 0 and 1.

*η*

_{ i }) will be fixed at a value of zero and a naïve-pooled data analysis will be conducted. This approach has also been suggested as a check for nonlinear mixed-effects models of dichotomous outcome data[28]. The predicted likelihood (

*l*

_{ ij }for individual

*i*and sampling time

*j*) of the data (

*y*

_{ ij }), given the model and parameters will be described by a binomial probability density function:

The likelihood for the entire population PD data set is simply the product of likelihoods across all individuals and data points. Diagnostic plots and the minimum value of the objective function are used to guide model building and assess goodness-of-fit.

### Methotrexate (MTX) Model

The administration of methotrexate (MTX) to children with cancer was chosen as our initial setting to develop the first drug dashboard prototype. The difficulty in effectively administering high-dose MTX to oncology patients lie in balancing efficacy and safety. Increased MTX exposure has been shown to be predictive of greater efficacy [29–31], while increased MTX concentrations and prolonged exposure time have also been linked to toxicity [32]. Due to the high inter- and intra-patient variability in methotrexate pharmacokinetics, monitoring of methotrexate plasma concentrations in individual patients has become a standard procedure in order to identify patients at risk of toxicity. Typically, patient plasma concentrations are monitored starting at 24 hours post infusion until MTX plasma concentrations fall below 0.1 to 1 μM [33–37], with adjunct rescue therapy implemented as needed.

The occurrence of methotrexate-induced renal toxicity further complicates chemotherapy administration. Although methotrexate-induced nephrotoxicity is a relatively rare occurrence, it is none-the-less a life threatening complication of methotrexate therapy [38]. Since methotrexate is mainly cleared from systemic circulation via glomerular filtration and renal secretion, delayed drug elimination is a product of this nephrotoxicity. This results in prolonged drug exposure and elevated plasma concentrations. As a result of this increased exposure, severe adverse events such as myelosuppression, mucositis, and hepatitis become more prevalent and severe.

Numerous studies have been conducted to examine the feasibility and reliability of applying Bayesian forecasting approaches to predicting MTX pharmacokinetics. The goals of these studies have been to predict MTX concentrations at later times or the time that MTX concentrations fall below a threshold value [39–41], MTX dose adjustment [42, 43], or providing guidance for rescue administration in the case of elevated MTX concentrations for prolonged time periods [44]. The Bayesian prediction models developed thus far have concentrated on those patients with normal renal function, and are not applicable in the case of severe renal dysfunction secondary to high-dose MTX administration.

We have developed a population pharmacokinetic model to implement as a Bayesian predictor of MTX concentrations in patients with normal renal function and MTX-induced renal dysfunction. Plasma concentrations from patients with normal renal function and patients with MTX-induced renal dysfunction were obtained from standard clinical monitoring. The model was constructed from methotrexate dosing histories and monitored drug concentrations in 240 patients. The original dataset contained 2176 observations covering a range of one to 56 observations per patient (an average of 9 observations per patient). The age range was from 1 to 80 years with a weight range of 6.6 to 157 kg. The gender distribution was approximately 48% male (52% female). Hence, our underlying patient diversity allowed us to include and consider relevant size and demographic dependencies. The model was developed using NONMEM version VI[45].

where:

*P*
_{
i
}is the estimated parameter value for individual *i*

$\widehat{P}$ is the typical population value (geometric mean) of the parameter

*η*
^{
Pi
}are individual-specific interindividual random effects for individual *i* and parameter *P* and are assumed to be distributed: *η* ~ *N(0, ω*
^{
2
}
*)* with covariance defined by the inter-individual covariance matrix Ω.

where:

*C*
_{
ij
}is the *j* th measured observation in individual *i*

${\widehat{C}}_{ij}$ is the *j* th model predicted value in individual *i*

*ε*
_{
ij
}is the additive residual random error for individual *i* and measurement *j* and is assumed to be independently and identically distributed

The Bayesian forecasting model utilizes the NONMEM PRIOR subroutine to incorporate population priors into the model. Fixed effects parameters obtained from the final pop PK model were implemented for the initial Bayesian model. Prior distributions of the fixed effects parameters were obtained from the variance-covariance matrix from the final pop PK model as well. Prior distributions for random effects parameters were specified as an inverse Wishart distribution. Clearance was implemented as a mixture model, where a patient is assigned to a population (normal or impaired clearance) based on the probability of that patient belonging to either population given their MTX plasma concentrations. The Bayesian forecasting model was evaluated using MTX plasma concentrations that were not used during model construction. The model reliably predicts future MTX plasma concentrations from two prior concentrations in all patients except a small number who develop renal toxicity at delayed times (> 48 hours). In these patients, the addition of a third concentration after 48 hours increases the precision of the prediction of concentrations at later times.

### Dashboard Design and System Architecture

The middle tier consists of rules and processing logic required to collect and prepare data for user presentation (alerts, filters, aggregations, derived values and predictions). Predictions are conducted in an external computational platform – our modeling and simulation (M&S) workbench. This platform can execute code in a variety of languages provided they can run in a batch mode. Of note, the M&S workbench can currently accommodates many of the standard prediction engines used to forecast PK and PK/PD relationships (NONMEM, SAS, SPLUS and R). Details of analytical run processing using NONMEM with the workbench are described below. While the workbench can perform various data processing functions and analytics including generation of plots and figures, it is important to note that all PKB related analytics are gated in the middle tier through logic to ensure that minimally required data sets are available for each patient or sets of patients for meaningful analysis (e.g., appropriate data density to make predictions, etc).

The user interface is a currently web-based and utilizes a combination of HTML, JavaScript and XML content. We are in the process of migrating to an AJAX paradigm and implementing Web 2.0 standards (the new standard for live HTML content). Upon successful transition to Web 2.0, it is possible to readily support more user interface types such as java applet, java application or stay completely within the browser. The Workbench component can co-reside on the same server or installed on another server within the intranet or even hosted anywhere on the internet. The system is designed to transfer data and messages in an authenticated and completely secure manner. When a physician invokes a given forecasting function within PKB, the data required of that patient for the analysis of interest is extracted and written to a flat file in the required application format. The system then edits a previously prepared control file (with job settings and analytical inputs) with necessary changes and makes a call to the Workbench requesting an execution. The communication and data transfer mechanisms are implemented via message queues, thereby ensuring guaranteed delivery, fault tolerance and scalability. The Workbench execution module listening to the queue receives the request and submits the jobs to NONMEM or whatever application has been specified for that analysis by the experts and calls the process manager module to monitor for completion. Upon completion, the process manager calls an application specific parser module to parse and return the results back to the calling application (middle-tier). The middle-tier receives the results and formats it into an XML object and returns the same to the browser for display. If the output of the analysis calls for graphical display of results, the middle-tier sends the unformatted output to the graphing engine and forwards the generated graphs back to the browser for display. The graphs can be setup to support drill downs, mouseovers or invoke further analysis. This ensures that analysis of any kind can be performed via the M&S Workbench to support drug/disease dash boards of varying complexities.

## Results

Current population pharmacokinetic parameter priors used to forecast methotrexate plasma concentrations in pediatric patients

Parameter | Units | Estimate | BSV |
---|---|---|---|

CL | L/h | 8.13 | 41.2% |

CL | L/h | 2.59 | 82.0% |

V1 | L | 39.6 | 21.0% |

V2 | L | 3.94 | 47.6% |

Q | L/h | 0.113 | 7.56% |

Figure 4B is the view projected after the dosing guidance menu button is selected. The plot shows the observed 24 and 48 hour MTX plasma concentrations along with the model-predicted MTX exposure (solid line) at these time points and at a subsequent, extrapolated 96 hour point. This extrapolated time point is estimated by calling the MTX population PK model and running the model executable with the patient's previous (observed data) incorporated. The refitted model with updated "patient-specific priors" is then used to project (simulate) the exposure at the 96 hour point, a time when blood collection for MTX plasma concentration determination is normally scheduled in compliance with formulary monitoring practice for MTX.

## Discussion

The development of drug-specific dashboards to educate patient caregivers on principals of clinical pharmacology and guide pharmacotherapy in pediatric populations is likely to yield superior clinical outcomes (fewer medication errors, reduced toxicity, reduced length of hospital stay, etc). While this approach has been advocated for some time and pioneering work by Jelliffe and others[25, 26, 46] has long demonstrated the clinical benefit of model-based dosing guidance, this research has not yielded any sustainable impact. The dashboard system proposed herein will rely heavily on the integration of modeling and simulation approaches in order to provide meaningful decision analytics to the end user. Our prototype methotrexate dashboard assembles the most clinically relevant patient data into an interface that allows end users to assess relevant biomarkers against drug exposure and forecast future exposures from a given dosing regimen as opposed to waiting to measure such levels via traditional TDM approach. Hence, an earlier assessment of the potential for nephrotoxicity can be made. The underlying population pharmacokinetic model was defined based on limited pediatric data which will be rechallenged prior to final model qualification/validation and production release of the MTX dashboard.

Future considerations for dashboard concepts will include functionality to predict the likelihood of drug interaction with co-administered drugs. By simulating virtual drug interaction studies, we will have the ability to report potential adverse events based on data mining and correlation analysis. We also envision the necessity of expanding the dashboard design construct to accommodate multiple agents considered as treatment options for targeted indications. In this instance, the choice of agent would be a decision criterion to be evaluated prior to the initiation of drug therapy. Our concept here would be to create workflows that considered patient status and previous pharmacotherapy outcomes along with criteria for ranking agent choice depending on the selection attributes. This situation is quite common with antibiotic therapy, which is already receiving attention with respect to commercially available solutions to tracking and prescribing. TheraDoc [47–49], Cereplex [50] and MedMined [51, 52] all represent commercial solutions in this arena. TheraDoc mines patient data for trends in infections and suggests courses of action for particular patients, while Cereplex searches for unusual infection patterns and identifies patients requiring changes in therapy and MedMined uses data-mining algorithms to tease out unusual patterns and correlations from patient records and lab tests. The development of decision support systems for managing antibiotic therapy spans several decades now[6, 53–55]. Much of the impetus for such systems has been the desire to respond clinically to dynamic changes in local or global bacterial prevalence as well as develop strategies to combat resistance. Likewise, as the landscape of therapeutic options changes with the introduction of new antibiotics, new data on additional indications (e.g., efficacy against new bacterial strains), epidemiologic data on cure rates, global/regional resistance development, and/or the exodus of agents from the market, such systems need some level of continued support beyond information technology. Ownership, governance and preventive maintenance efforts must become formalized for such systems to continue to provide the same level of guidance as when they were first implemented.

The broad array of decision support systems currently employed in hospital settings coupled with those in development highlights the need for robust data integration and flexible decision analytics validated against all possible conditions of use and practice. While the concept of modeling and simulation integration is relatively straightforward, the details of ensuring the performance of these systems, particularly those that impart clinical guidance are complex and require input from IT, clinical pharmacology, pharmacy and clinical practice. The governance of our efforts is overseen by our IRB and therapeutics standards committee but this alone does not ensure the practical issues associated with guiding pharmacotherapy. Given the paucity of information often available to guide pediatric pharmacotherapy, there is a strong desire to "fill-in" such gaps with the best available information available. Likewise, the void in data and knowledge today does not imply that such gaps will remain and decision analytics provided to guide present pharmacotherapy must be revisited as new information becomes available.

## Conclusion

The integration of modeling and simulation algorithms with hospital-based networks to guide the pharmacotherapeutic management of individual patients has great potential to improve outcomes. The benefits of such a system should include improved therapy (efficacy), reduced medication errors, greater appreciation for drug interaction potential, earlier identification of toxicity, and earlier guidance on rescue therapy. Our prototype dashboard concept is part of a broader initiative to develop a pediatric knowledgebase of which dashboards are only one component. At present, dashboards for methotrexate, tacrolimus and vancomycin are at various stages of development. Collaborations with other institutions and investigators should allow the generation of additional dashboards beyond our existing capacity. Such systems, as they are developed, will require a level of support beyond which many hospitals are accustomed as mentioned previously. More importantly, they imply a continued effort from clinical pharmacologists and engineers to implement relevant, new research into these systems to ensure that they continue to perform up to expectations and evolve with advances in drug therapy in pediatrics.

## Declarations

### Acknowledgements

This work was funded in part from an internal grant, the Pediatrics Chair's Initiative, The Children's Hospital of Philadelphia.

## Authors’ Affiliations

## References

- IOM: Building a Better Delivery system: A New Engineering/Health Care Partnership. 2005, National Academy of Engineering and Institute of MedicineGoogle Scholar
- IOM: Report: Patient safety – achieving a new standard for care. Acad Emerg Med. 2005, 12 (10): 1011-2.Google Scholar
- Ferri F: User modeling techniques as support in the clinical decision-making process. Medinfo. 1995, 8 (Pt 2): 926-30.PubMedGoogle Scholar
- Gardner SN: Modeling multi-drug chemotherapy: tailoring treatment to individuals. J Theor Biol. 2002, 214 (2): 181-207. 10.1006/jtbi.2001.2459.View ArticlePubMedGoogle Scholar
- Starr JM, Campbell A: Mathematical modeling of Clostridium difficile infection. Clin Microbiol Infect. 2001, 7 (8): 432-7. 10.1046/j.1198-743x.2001.00291.x.View ArticlePubMedGoogle Scholar
- Broderick A: Nosocomial infections: validation of surveillance and computer modeling to identify patients at risk. Am J Epidemiol. 1990, 131 (4): 734-42.PubMedGoogle Scholar
- Raboud J: Modeling transmission of methicillin-resistant Staphylococcus aureus among patients admitted to a hospital. Infect Control Hosp Epidemiol. 2005, 26 (7): 607-15. 10.1086/502589.View ArticlePubMedGoogle Scholar
- Shaw B, Marshall AH: Modeling the health care costs of geriatric inpatients. IEEE Trans Inf Technol Biomed. 2006, 10 (3): 526-32. 10.1109/TITB.2005.863821.View ArticlePubMedGoogle Scholar
- Graves N, Nicholls TM, Morris AJ: Modeling the costs of hospital-acquired infections in New Zealand. Infect Control Hosp Epidemiol. 2003, 24 (3): 214-23. 10.1086/502192.View ArticlePubMedGoogle Scholar
- Anderson JG: Modeling the costs and outcomes of cardiovascular surgery. Health Care Manag Sci. 2002, 5 (2): 103-11. 10.1023/A:1014472731382.View ArticlePubMedGoogle Scholar
- Hung GR: Computer modeling of patient flow in a pediatric emergency department using discrete event simulation. Pediatr Emerg Care. 2007, 23 (1): 5-10. 10.1097/PEC.0b013e31802c611e.View ArticlePubMedGoogle Scholar
- Costa AX: Mathematical modelling and simulation for planning critical care capacity. Anaesthesia. 2003, 58 (4): 320-7. 10.1046/j.1365-2044.2003.03042.x.View ArticlePubMedGoogle Scholar
- Koizumi N, Kuno E, Smith TE: Modeling patient flows using a queuing network with blocking. Health Care Manag Sci. 2005, 8 (1): 49-60. 10.1007/s10729-005-5216-3.View ArticlePubMedPubMed CentralGoogle Scholar
- Isken MW, Rajagopalan B: Data mining to support simulation modeling of patient flow in hospitals. J Med Syst. 2002, 26 (2): 179-97. 10.1023/A:1014814111524.View ArticlePubMedGoogle Scholar
- Paul JA: Transient modeling in simulation of hospital operations for emergency response. Prehospital Disaster Med. 2006, 21 (4): 223-36.View ArticlePubMedGoogle Scholar
- Van Houdenhoven M: Optimizing intensive care capacity using individual length-of-stay prediction models. Crit Care. 2007, 11 (2): R42-10.1186/cc5730.View ArticlePubMedPubMed CentralGoogle Scholar
- Braithwaite RS, Roberts MS, Justice AC: Incorporating quality of evidence into decision analytic modeling. Ann Intern Med. 2007, 146 (2): 133-41.View ArticlePubMedPubMed CentralGoogle Scholar
- Connelly LG, Bair AE: Discrete event simulation of emergency department activity: a platform for system-level operations research. Acad Emerg Med. 2004, 11 (11): 1177-85.View ArticlePubMedGoogle Scholar
- Slovensky DJ, Morin B: Learning through simulation: the next dimension in quality improvement. Qual Manag Health Care. 1997, 5 (3): 72-9.View ArticlePubMedGoogle Scholar
- Lucas CE: Mathematical modeling to define optimum operating room staffing needs for trauma centers. J Am Coll Surg. 2001, 192 (5): 559-65. 10.1016/S1072-7515(01)00829-8.View ArticlePubMedGoogle Scholar
- Verduijn M: Modeling length of stay as an optimized two-class prediction problem. Methods Inf Med. 2007, 46 (3): 352-9.PubMedGoogle Scholar
- Chien JY: Pharmacokinetics/Pharmacodynamics and the stages of drug development: role of modeling and simulation. Aaps J. 2005, 7 (3): E544-59. 10.1208/aapsj070355.View ArticlePubMedPubMed CentralGoogle Scholar
- Gobburu JV, Sekar VJ: Application of modeling and simulation to integrate clinical pharmacology knowledge across a new drug application. Int J Clin Pharmacol Ther. 2002, 40 (7): 281-8.View ArticlePubMedGoogle Scholar
- Meibohm B: Population pharmacokinetic studies in pediatrics: issues in design and analysis. Aaps J. 2005, 7 (2): E475-87. 10.1208/aapsj070248.View ArticlePubMedPubMed CentralGoogle Scholar
- Bondareva IB: Nonparametric population modeling of valproate pharmacokinetics in epileptic patients using routine serum monitoring data: implications for dosage. J Clin Pharm Ther. 2004, 29 (2): 105-20. 10.1111/j.1365-2710.2003.00538.x.View ArticlePubMedGoogle Scholar
- Jelliffe R: Goal-oriented, model-based drug regimens: setting individualized goals for each patient. Ther Drug Monit. 2000, 22 (3): 325-9. 10.1097/00007691-200006000-00016.View ArticlePubMedGoogle Scholar
- Agresti A: A survey of models for repeated ordered categorical response data. Stat Med. 1989, 8 (10): 1209-24. 10.1002/sim.4780081005.View ArticlePubMedGoogle Scholar
- Yano I, Beal SL, Sheiner LB: The need for mixed-effects modeling with population dichotomous data. J Pharmacokinet Pharmacodyn. 2001, 28 (4): 389-412. 10.1023/A:1011586814601.View ArticlePubMedGoogle Scholar
- Borsi J, Revesz T, Schuler D: [Prognostic significance of systemic clearance of methotrexate in acute lymphoid leukemia in childhood]. Orv Hetil. 1986, 127 (8): 439-42.PubMedGoogle Scholar
- Evans WE: Methotrexate systemic clearance influences probability of relapse in children with standard-risk acute lymphocytic leukaemia. Lancet. 1984, 1 (8373): 359-62. 10.1016/S0140-6736(84)90411-2.View ArticlePubMedGoogle Scholar
- Evans WE: Clinical pharmacodynamics of high-dose methotrexate in acute lymphocytic leukemia. Identification of a relation between concentration and effect. N Engl J Med. 1986, 314 (8): 471-7.View ArticlePubMedGoogle Scholar
- Goldie JH, Price LA, Harrap KR: Methotrexate toxicity: correlation with duration of administration, plasma levels, dose and excretion pattern. Eur J Cancer. 1972, 8 (4): 409-14.View ArticlePubMedGoogle Scholar
- Tattersall MH, Brown B, Frei E: The reversal of methotrexate toxicity by thymidine with maintenance of antitumour effects. Nature. 1975, 253 (5488): 198-200. 10.1038/253198a0.View ArticlePubMedGoogle Scholar
- Isacoff WH: Pharmacokinetics of high-dose methotrexate with citrovorum factor rescue. Cancer Treat Rep. 1977, 61 (9): 1665-74.PubMedGoogle Scholar
- Stoller RG: Use of plasma pharmacokinetics to predict and prevent methotrexate toxicity. N Engl J Med. 1977, 297 (12): 630-4.View ArticlePubMedGoogle Scholar
- Nirenberg A: High-dose methotrexate with citrovorum factor rescue: predictive value of serum methotrexate concentrations and corrective measures to avert toxicity. Cancer Treat Rep. 1977, 61 (5): 779-83.PubMedGoogle Scholar
- Perez C: Significance of the 48-hour plasma level in high-dose methotrexate regimens. Cancer Clin Trials. 1978, 1 (2): 107-11.PubMedGoogle Scholar
- Widemann BC: High-dose methotrexate-induced nephrotoxicity in patients with osteosarcoma. Cancer. 2004, 100 (10): 2222-32. 10.1002/cncr.20255.View ArticlePubMedGoogle Scholar
- Piard C: A limited sampling strategy to estimate individual pharmacokinetic parameters of methotrexate in children with acute lymphoblastic leukemia. Cancer Chemother Pharmacol. 2007, 60 (4): 609-20. 10.1007/s00280-006-0394-3.View ArticlePubMedGoogle Scholar
- Rousseau A: Bayesian estimation of methotrexate pharmacokinetic parameters and area under the curve in children and young adults with localised osteosarcoma. Clin Pharmacokinet. 2002, 41 (13): 1095-104. 10.2165/00003088-200241130-00006.View ArticlePubMedGoogle Scholar
- Odoul F: Prediction of methotrexate elimination after high dose infusion in children with acute lymphoblastic leukaemia using a population pharmacokinetic approach. Fundam Clin Pharmacol. 1999, 13 (5): 595-604.View ArticlePubMedGoogle Scholar
- Pignon T: Pharmacokinetics of high-dose methotrexate in adult osteogenic sarcoma. Cancer Chemother Pharmacol. 1994, 33 (5): 420-4. 10.1007/BF00686272.View ArticlePubMedGoogle Scholar
- Bruno R: Dosage predictions in high-dose methotrexate infusions. Part 2: Bayesian estimation of methotrexate clearance. Cancer Drug Deliv. 1985, 2 (4): 277-83.View ArticlePubMedGoogle Scholar
- Monjanel-Mouterde S: Bayesian population model of methotrexate to guide dosage adjustments for folate rescue in patients with breast cancer. J Clin Pharm Ther. 2002, 27 (3): 189-95. 10.1046/j.1365-2710.2002.00402.x.View ArticlePubMedGoogle Scholar
- Bauer RJ, Guzy S, Ng C: A survey of population analysis methods and software for complex pharmacokinetic and pharmacodynamic models with examples. Aaps J. 2007, 9 (1): E60-83. 10.1208/aapsj0901007.View ArticlePubMedPubMed CentralGoogle Scholar
- Jelliffe RW: Adaptive control of drug dosage regimens: basic foundations, relevant issues, and clinical examples. Int J Biomed Comput. 1994, 36 (1–2): 1-23. 10.1016/0020-7101(94)90091-4.View ArticlePubMedGoogle Scholar
- Pestotnik SL: Expert clinical decision support systems to enhance antimicrobial stewardship programs: insights from the society of infectious diseases pharmacists. Pharmacotherapy. 2005, 25 (8): 1116-25. 10.1592/phco.2005.25.8.1116.View ArticlePubMedGoogle Scholar
- Staes CJ: A case for manual entry of structured, coded laboratory data from multiple sources into an ambulatory electronic health record. J Am Med Inform Assoc. 2006, 13 (1): 12-5. 10.1197/jamia.M1813.View ArticlePubMedPubMed CentralGoogle Scholar
- TheraDoc. [http://www.theradoc.com]
- Cereplex. [http://www.cereplex.com]
- Brossette SE: A laboratory-based, hospital-wide, electronic marker for nosocomial infection: the future of infection control surveillance?. Am J Clin Pathol. 2006, 125 (1): 34-9. 10.1309/502A-UPR8-VE67-MBDE.View ArticlePubMedGoogle Scholar
- MedMined. [http://www.medmined.com]
- Evans RS, Pestotnik SL: Applications of medical informatics in antibiotic therapy. Adv Exp Med Biol. 1994, 349: 87-96.View ArticlePubMedGoogle Scholar
- Evans RS: Improving empiric antibiotic selection using computer decision support. Arch Intern Med. 1994, 154 (8): 878-84. 10.1001/archinte.154.8.878.View ArticlePubMedGoogle Scholar
- Pestotnik SL: Implementing antibiotic practice guidelines through computer-assisted decision support: clinical and financial outcomes. Ann Intern Med. 1996, 124 (10): 884-90.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6947/8/6/prepub

### Pre-publication history

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.