Integration of SAS and NONMEM for Automation of Population Pharmacokinetic/Pharmacodynamic Modeling on UNIX systems
|
|
- Eustace O’Connor’
- 5 years ago
- Views:
Transcription
1 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, Buffalo NY ABSTRACT In drug development and other pharmaceutical applications, SAS is a powerful tool for statistical analysis, graph creation, text file processing, and system command submission, while NONMEM is a powerful software tool for population pharmacokinetic/ pharmacodynamic (PK/PD) modeling. The integration of SAS and NONMEM is therefore a natural solution for automation or semiautomation of population PK/PD modeling. Population PK/PD modeling, especially the process of covariate evaluation (including forward selection of the most significant covariate and backward elimination of the most non-significant covariate), can be automated using SAS programs as a recycling process of running NONMEM programs, comparing results, updating programs, and rerunning. For example, the results comparison step can be implemented via SAS macros that read NONMEM output files, list and sort the minimum objective function values, and aid in the identification of the next base model. The flowchart of the population PK/PD modeling, algorithm of the automation, and SAS and NONMEM code segments will be demonstrated in this paper. The automation is especially efficient and time-saving for population PK/PD analyses with many covariates to be evaluated. An additional benefit of this automated solution is the potential for reduction in error opportunities and the need for extensive quality assurance. INTRODUCTION Population pharmacokinetic/pharmacodynamic (PK/PD) modeling is becoming more and more important in the healthcare/pharmaceutical industry because necessary information for decision-making, drug development optimization, and clinical dosage individualization can be effectively obtained using population PK/PD analysis. Different methods for population PK/PD modeling have been proposed 1-4. Whatever method is used, the flowchart of the population PK/PD modeling process is similar to the one as illustrated in Figure 1. As is shown in the flowchart, one of the important processes in population PK/PD model development is covariate evaluation, identifying the relationship between the model parameters and covariates, such as demographic factors, clinical laboratory biochemical measurements, organ functionality, genotype/phenotype, and drug-drug interactions. This process is time-consuming and could account for up to 90% of the workload and total time for the whole population PK/PD analysis, depending on the complexity of the population PK/PD model structure, number of significant covariates, and methods used for covariate evaluation. The complexity of the population PK/PD model structure and number of covariates are generally determined by the complexity of the drug s metabolism, population PK/PD modeling strategy (such as physiologically-based or mechanism-based population PK/PD modeling), and availability of clinical trial data. Given a drug, the time required for covariate evaluation is mainly determined by the method used for this evaluation, which is commonly implemented by NONMEM programs 5. However, implementation of some methods might be very difficult in practice under certain circumstances. For example, the Likelihood Ratio Test (LRT) method with Wald s approximation proposed by Kowalski and Hutmacher 1 starts from a full model all relevant covariates are initially included, followed by a process of backward elimination of the most non- significant covariates from the model. Population PK/PD modeling with this method could never proceed if too many covariates are included in the full model or initial guesses for the model parameters are very different from the true values since runs might never converge. Therefore, selection of methodology for covariate evaluation is important and critical in determining the computation workload and time. To avoid the problems described above, this paper will discuss a commonly used method for covariate evaluation: stepwise forward selection followed by stepwise backward elimination. Data characterization SAS/NONMEM dataset creation Covariate distribution/statistics Drug concentration-time profiles Model Development Structural model Covariate evaluation Forward selection Backward elimination Random effect evaluation Model refinement Model validation/test/simulation Figure 1. Flowchart of population PK/PD modeling
2 The total workload and time for population PK/PD modeling can be split into two parts: machine (computation) workload/time and scientist workload/time. The description in the previous paragraph refers to machine workload and time. This paper will only focus on how to reduce the scientist s workload and time for population PK/PD modeling through automation using SAS programs. Specifically, automation of the process of covariate evaluation will be illustrated. FEASIBILITY OF AUTOMATING COVARIATE EVALUATION The issue of reducing a scientist s workload and time becomes important when the number of potential significant covariates is large and/or the number of model parameters on which covariate effects will be evaluated is large. A typical process of covariate evaluation via forward selection and backward elimination can be expressed as a recycling process as illustrated in Figure 2. Base model (base NONMEM program) Derived models (derived NONMEM programs) Running derived NONMEM programs Comparing running results of the derived programs More covariates? Yes Update base model (updated base NONMEM program) No Covariate Evaluation done Figure 2. Flowchart of forward selection or backward elimination Each derived/updated NONMEM program has only one more covariate included/removed from the previous base program. Therefore, the whole process of covariate evaluation can be easily automated by SAS macros. The structure of a NONMEM program which is composed of different functional blocks also makes automation feasible. Forward Selection During forward selection, the derived NONMEM programs are different from the base model by only one covariate term which is added on some model parameter and one more line in the $THETA block for the initial value of that covariate, as illustrated below. Segment of the base NONMEM program: $PK TVCL = THETA(1)+THETA(4)*(SGOT/20-1) CL = TVCL*EXP(ETA(1)) TV21 = THETA(2)+THETA(5)*(SEX- 0.53)+THETA(6)*(AGE/52-1) TVV = TV21+THETA(7)*(RAC1-0.69) V = TVV*EXP(ETA(2)) TVKA = THETA(3)+THETA(8)*(SGOT/20-1) KA = TVKA*EXP(ETA(3)) $THETA (0, 30); THETA(1) CLEARANCE (L/H) (0, 150); THETA(2) CENTRAL VOLUME (L) (0, 0.7); THETA(3) KA (PER H) (-20) (20) (-0.02) One derived NONMEM program: $PK TVCL = THETA(1)+THETA(4)*(SGOT/20-1)+THETA(9)*(RAC2-0.43) CL = TVCL*EXP(ETA(1)) TV21 = THETA(2)+THETA(5)*(SEX- 0.53)+THETA(6)*(AGE/52-1) TVV = TV21+THETA(7)*(RAC1-0.69) V = TVV*EXP(ETA(2)) TVKA = THETA(3)+THETA(8)*(SGOT/20-1) KA = TVKA*EXP(ETA(3)) $THETA (0, 30); THETA(1) CLEARANCE (L/H) (0, 150); THETA(2) CENTRAL VOLUME (L) (0, 0.7); THETA(3) KA (PER H) (-20) (20) (-0.02) (2) 2
3 . Compared to the base program, one term for the effect of a covariate has been added on a selected PK parameter and one initial guess for the coefficients in the added covariate term are appended at the end of the $THETA block (The number of the appended initial guess in the $THETA block should be equal to the number of the coefficients contained in the covariate term). The above updates can be implemented with a SAS macro like: %trnsctlm(path, baseprogramfilename, newprogramfilename, covariatetobeevaluated, covariatetermtobeadded, ParameterOnWhichTheCovariateTermAdded, initialguessforthecovariatecoefficient); This macro %trnsctlm implements the following tasks: 1. Read variable Xline ($80.) from a data file &path&baseprogramfilename line by line; 2. If index(xline, &ParameterOnWhichTheCovariateTermAdded ) = 1 and length(&covariatetermtobeadded) <79 - length(trim(xline)) <80 then a. Xline = Xline &covariatetermtobeadded ; Else do; Xline = newaccumulatetermforparameter substring(xline, length(¶meteronwhichthecovariateter madded)); Insert a line with: &ParameterOnWhichTheCovariateTermAdded = newaccumulatetermforparameter &covariatetermtobeadded ; 3. Append Xline = &initialguessforthecovariatecoefficient at the end of the block $THETA ; 4. Output the derived NONMEM program with a new file name &newprogramfilename. The derived NONMEM program in the previous example can be achieved by calling this macro with the macro variables being replaced with appropriate values, such as: &covariatetobeevaluated = RAC2 ; &covariatetermtobeadded = +THETA(9)*(RAC2-0.43) ; &ParameterOnWhichTheCovariateTermAdded = TVCL ; &initialguessforthecovariatecoefficient = (2) ; and so on. Running Derived NONMEM Programs Once a batch of new NONMEM programs are created, they can be started to run in parallel and/or in a queue simply by implementing the following SAS command: x cd /dir0/dir1/dir2/dir3/dir4/dir5/ x./nonmemcommand >& nonmemcommand.o & The script file nonmemcommand contains the command to run the derived NONMEM programs, such as # cd /dir0/dir1/dir2/dir3/dir4/dir5/dir6/ nonmem derivednonmemprogram1 nonmem derivednonmemprogram2 nonmem derivednonmemprogram3 nonmem derivednonmemprogram4 nonmem derivednonmemprogram5 # Backward Elimination It is easier to update NONMEM programs during backward elimination since the actual term which expresses the effect of a covariate is not necessarily physically removed from the program. The elimination of a covariate term from the program can be implemented just by fixing the corresponding coefficient(s) to zero(s), as expressed in this algorithm: 1. Read variable Xline ($80) from a data file: &path&baseprogramfilename line by line 2. If a covariate is selected to be removed, then set the corresponding initial guesses in the $THETA block to (0 FIXED). 3. Output the derived NONMEM program with a new file name &newprogramfilename. Results Comparison The purpose of the covariate evaluation process is to identify the most significant covariates for the PK of a drug. During the forward selection process, significant covariates are added to the PK model while, during the process of backward elimination, insignificant covariates are removed from the model. Therefore, to reach this goal, comparisons between the minimum objective function values of different NONMEM programs are required, which can be implemented by a SAS macro program that reads the objective function values from the output of each NONMEM program. The algorithm can be expressed as: 1. Read the minimum objective function value of the base NONMEM program from its NONMEM report file 2. Read the minimum objective function value of each derived NONMEM program and find the difference in the degrees of freedom between the derived program and the base program 3
4 3. Calculate p-values, based on the difference in minimum objective function values and the difference in degrees of freedom using the likelihood ratio test approach based on an asymptotic χ 2 distribution. 4. Identify the most significant covariate during forward selection of covariates (with the smallest p-value) or the most non-significant covariate during backward elimination (with the largest p-value) 5. The identified program is used as a new base program for further evaluation of covariate effects CASE STUDY Case Study 1: The data was from a sparse sampling Phase II clinical trial for drug XX1 with 331 patients enrolled. The drug was sequentially converted to two active/toxic metabolites in the human body. A five-compartment model was identified to best describe the drug s PK. Covariate effects were estimable on six PK parameters in this model. It took about 3 hours to finish running the structural model on a high speed UNIX system, with longer time to run when covariates were included in the PK model. According to the exploratory investigation and prior information, these six PK parameters are potentially associated with some or all of the 15 covariates. Therefore, for a full spectrum of covariate evaluation (main effect), there are 15*6 = 90 potential programs for the first round, not including the potential interactions between covariates. Even after reduction with physiological rationale, the total number of potential programs was still up to 50. Manual editing of the base program to generate these 50 programs would have been extremely time-consuming and prone to error. Traditionally, a series of graphs to evaluate the potential correlations between each PK parameter and each covariate are created to visually judge the significance of the covariate and the trend of the PK parameter with the covariate values so as to further largely reduce the number of covariates on each PK parameter. However, this method is disadvantageous since: 1). There are many graphs to be visually compared; 2). It is difficult to compare the significance of a covariate on a parent PK parameter with that on a metabolite PK parameter or the significance of a covariate on two different metabolite PK parameters; 3). It is difficult to compare different covariates on the same PK parameter when they are at different levels of magnitude, such as bilirubin at a level of 10-1 mg/dl while body weight is at a level of 10 2 kg. 4). It is almost impossible to visually decide which is more significant across covariates and species, such as bilirubin on the clearance of the parent drug and level of one enzyme on the metabolic rate constant. 5). Even with the number of total significant covariates largely reduced, whether or not those unselected covariates might be significant can not be confirmed. With the long run time for each NONMEM program, a scientist may have enough time to manually update NONMEM programs so that his computer does not have to wait for available programs to run. That is to say, the limiting process of the covariate evaluation is computation time. However, with manual update of NONMEM programs, the scientist has to be tied to the program updating and results comparison all the time. Furthermore, errors could frequently occur during the manual update of programs and retrieving/comparison of results. In contrast, with automation, the time that the scientist spends on updating programs and retrieving/comparing results is negligible. In addition, error occurrence could be largely reduced. As a result, the scientist could save time to perform other tasks. Case 2: The data was from three rich data Phase I clinical trials of drug XX2 (oral administration in different formulation and age groups) with 100 young patients enrolled. One active metabolite was identified for this drug. A two-compartment model was identified to best describe the drug PK. Altogether 21 covariates are potentially associated with 6 PK parameters. With physiological rationale, the total number of possible evaluations on effects of covariates on PK parameters could be reduced from 21*6=126 to 77. The computation time for the structural model was about 4 minutes. Obviously, program updating would be a rate limiting process for the evaluation of covariate effects if manual updating were conducted. With automation, computation time becomes the rate limiting process while the scientist is liberated from burden workload. Case 3: The data was from a sparse sampling Phase II clinical trial of drug XX3 (oral administration) with 200 patients enrolled. No metabolite was of concern. A onecompartmental model was identified to best describe the drug s PK. Altogether 14 covariates were potentially associated with 3 PK parameters. The total number of possible evaluations of the effect of covariates on PK parameters of the structural model was 30. It took about 20 seconds to run the structural model. In this case, program updating and results comparison were the rate limiting process for the covariate evaluation. Manual updating of programs and comparisons of results would require much more time than computation. That is, the computer would have to wait for the scientist to update programs for running. Even though the total number of programs could be further reduced via visual assessment, the total time for covariate evaluation would be still determined by manual update of programs. However, with automation, even though more programs were actually run, the total time for the covariate evaluation was still much less. SUMMARY Automation of program updating and the comparisons of results can be implemented with SAS macros and can save scientists a tremendous amount of time and dramatically decrease the opportunities for error occurrence. This automation does not change the total computation time for population PK/PD modeling, but largely reduces the 4
5 workload and the total time that a scientist spends on covariate evaluation and therefore reduces the total time for a project. Together with future automation of identifying structural models, integration of SAS and NONMEM provides a possibility for the complete automation of the whole life cycle of a population PK/PD analysis. Phone: ext. 228 Fax: jbf@cognigencorp.com Web: REFERENCE 1. K. G. Kowalski and M. M. Hutmacher. Efficient screening of covariates in population models using Wald's approximation to the likelihood ratio test. J Pharmacokinet Pharmacodyn 2001 Jun; 28(3): E. N. Jonsson and M. O. Karlson. Automated Covariate Model Building within NONMEM. Pharm Res 1998 Sep; 15 (9): J. W. Mandema, D. Verotta and L. B. Sheiner. Building Population Pharmacokinetic-pharmacodynamic Models. I. Models for Covariate Effects. J. Pharmacokin. Biopharm. 20: (1992) 4. J. F. Lawless and K. Singhal. Efficient Screening of Nonnormal Regression Models. Biometrics 34: (1978) 5. S. L. Beal and L. B. Sheiner (eds.). NONMEM users guides, NONMEM Project Group, University of California, San Francisco, CA, SAS and all other SAS Institute Inc. product r service names are registered trademarks or trade marks of SAS Institute Inc. in the USA and other countries. indicates USA registration. Other brand and product names are registered trademarks or trademarks of their respective companies. CONTACT INFORMATION The authors can be contacted at: Alan J Xiao, Ph.D. Population PK/PD Scientist Cognigen Corporation 395 Youngs Road Buffalo, NY Phone: ext. 265 Fax: alan.xiao@cognigencorp.com Web: Jill B Fiedler-Kelly, M.S. Vice President, Population PK/PD Cognigen Corporation 395 Youngs Road Buffalo, NY
Practice 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 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 informationMetabolite Identification and Characterization by Mining Mass Spectrometry Data with SAS and Python
PharmaSUG 2018 - Paper AD34 Metabolite Identification and Characterization by Mining Mass Spectrometry Data with SAS and Python Kristen Cardinal, Colorado Springs, Colorado, United States Hao Sun, Sun
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 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 informationDosing In NONMEM Data Sets an Enigma
PharmaSUG 2018 - Paper BB-02 ABSTRACT Dosing In NONMEM Data Sets an Enigma Sree Harsha Sreerama Reddy, Certara, Clarksburg, MD; Vishak Subramoney, Certara, Toronto, ON, Canada; Dosing data plays an integral
More informationA SAS/AF Application For Sample Size And Power Determination
A SAS/AF Application For Sample Size And Power Determination Fiona Portwood, Software Product Services Ltd. Abstract When planning a study, such as a clinical trial or toxicology experiment, the choice
More informationPerforming Deconvolution Operations in SAS for Input Rate Computation
ABSTRACT Paper PO11 Performing Deconvolution Operations in SAS for Input Rate Computation Steven Hege, Rho, Inc., Chapel Hill, NC. Deconvolution is a non-parametric technique used in the quantitative analysis
More informationPower calculation for non-inferiority trials comparing two Poisson distributions
Paper PK01 Power calculation for non-inferiority trials comparing two Poisson distributions Corinna Miede, Accovion GmbH, Marburg, Germany Jochen Mueller-Cohrs, Accovion GmbH, Marburg, Germany Abstract
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 informationThe SEQDESIGN Procedure
SAS/STAT 9.2 User s Guide, Second Edition The SEQDESIGN Procedure (Book Excerpt) This document is an individual chapter from the SAS/STAT 9.2 User s Guide, Second Edition. The correct bibliographic citation
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 informationPaper The syntax of the FILENAME is:
Paper 44-2010 Use of FILENAME statement to create automated reports with severe weather events data from the National Weather Service Valentin Todorov, Assurant Specialty Property, Atlanta, GA Abstract
More informationMachine search vs. humans in modeling or I m sorry Dave, I m afraid I. can t do that. Mark Sale M.D. Next Level Solutions.
Machine search vs. humans in modeling or I m sorry Dave, I m afraid I Mark Sale M.D. Next Level Solutions Nov 2, 2011 can t do that Current (local) search algorithm assumption is path/starting point independence
More informationPaper: ST-161. Techniques for Evidence-Based Decision Making Using SAS Ian Stockwell, The Hilltop UMBC, Baltimore, MD
Paper: ST-161 Techniques for Evidence-Based Decision Making Using SAS Ian Stockwell, The Hilltop Institute @ UMBC, Baltimore, MD ABSTRACT SAS has many tools that can be used for data analysis. From Freqs
More informationReaxys Improved synthetic planning with Reaxys
R&D SOLUTIONS Reaxys Improved synthetic planning with Reaxys Integration of custom inventory and supplier catalogs in Reaxys helps chemists overcome challenges in planning chemical synthesis due to the
More informationIntegrated Electricity Demand and Price Forecasting
Integrated Electricity Demand and Price Forecasting Create and Evaluate Forecasting Models The many interrelated factors which influence demand for electricity cannot be directly modeled by closed-form
More informationA GUI FOR EVOLVE ZAMS
A GUI FOR EVOLVE ZAMS D. R. Schlegel Computer Science Department Here the early work on a new user interface for the Evolve ZAMS stellar evolution code is presented. The initial goal of this project is
More informationStephen McDonald and Mark D. Wrona Waters Corporation, Milford, MA, USA INT RO DU C T ION. Batch-from-Structure Analysis WAT E R S SO LU T IONS
Chemical Intelligence in UNIFI: Enhancing in-silico Decision Making to Provide Confident and Accurate Results for Metabolite Identification Experiments Stephen McDonald and Mark D. Wrona Waters Corporation,
More informationMember Level Forecasting Using SAS Enterprise Guide and SAS Forecast Studio
Paper 2240 Member Level Forecasting Using SAS Enterprise Guide and SAS Forecast Studio Prudhvidhar R Perati, Leigh McCormack, BlueCross BlueShield of Tennessee, Inc., an Independent Company of BlueCross
More informationTraffic accidents and the road network in SAS/GIS
Traffic accidents and the road network in SAS/GIS Frank Poppe SWOV Institute for Road Safety Research, the Netherlands Introduction The first figure shows a screen snapshot of SAS/GIS with part of the
More informationChallenges in modelling the pharmacokinetics of isoniazid in South African tuberculosis patients
PAGE 2005 Pamplona, Spain Challenges in modelling the pharmacokinetics of isoniazid in South African tuberculosis patients Justin J Wilkins 1, Grant Langdon 1, Helen McIlleron 1, Goonaseelan Pillai 2,
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 informationQinlei Huang, St. Jude Children s Research Hospital, Memphis, TN Liang Zhu, St. Jude Children s Research Hospital, Memphis, TN
PharmaSUG 2014 - Paper SP04 %IC_LOGISTIC: A SAS Macro to Produce Sorted Information Criteria (AIC/BIC) List for PROC LOGISTIC for Model Selection ABSTRACT Qinlei Huang, St. Jude Children s Research Hospital,
More informationAn Approach for Identifiability of Population Pharmacokinetic-Pharmacodynamic Models
An Approach for Identifiability of Population Pharmacokinetic-Pharmacodynamic Models Vittal Shivva 1 * Julia Korell 12 Ian Tucker 1 Stephen Duffull 1 1 School of Pharmacy University of Otago Dunedin New
More informationGenerating Half-normal Plot for Zero-inflated Binomial Regression
Paper SP05 Generating Half-normal Plot for Zero-inflated Binomial Regression Zhao Yang, Xuezheng Sun Department of Epidemiology & Biostatistics University of South Carolina, Columbia, SC 29208 SUMMARY
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 informationAn Introduction to GLIF
An Introduction to GLIF Mor Peleg, Ph.D. Post-doctoral Fellow, SMI, Stanford Medical School, Stanford University, Stanford, CA Aziz A. Boxwala, M.B.B.S, Ph.D. Research Scientist and Instructor DSG, Harvard
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 informationDynamic Determination of Mixed Model Covariance Structures. in Double-blind Clinical Trials. Matthew Davis - Omnicare Clinical Research
PharmaSUG2010 - Paper SP12 Dynamic Determination of Mixed Model Covariance Structures in Double-blind Clinical Trials Matthew Davis - Omnicare Clinical Research Abstract With the computing power of SAS
More informationA New Approach to Modeling Covariate Effects and Individualization in Population Pharmacokinetics-Pharmacodynamics
Journal of Pharmacokinetics and Pharmacodynamics, Vol. 33, No. 1, February 2006 ( 2006) DOI: 10.1007/s10928-005-9000-2 A New Approach to Modeling Covariate Effects and Individualization in Population Pharmacokinetics-Pharmacodynamics
More informationPreparing Spatial Data
13 CHAPTER 2 Preparing Spatial Data Assessing Your Spatial Data Needs 13 Assessing Your Attribute Data 13 Determining Your Spatial Data Requirements 14 Locating a Source of Spatial Data 14 Performing Common
More informationEstimation and Model Selection in Mixed Effects Models Part I. Adeline Samson 1
Estimation and Model Selection in Mixed Effects Models Part I Adeline Samson 1 1 University Paris Descartes Summer school 2009 - Lipari, Italy These slides are based on Marc Lavielle s slides Outline 1
More informationwww.goldensoftware.com Why Create a Thematic Map? A thematic map visually represents the geographic distribution of data. MapViewer will help you to: understand demographics define sales or insurance territories
More informationHierarchical Anomaly Detection in Load Testing with StormRunner Load
White Paper Application Development, Test & Delivery Hierarchical Anomaly Detection in Load Testing with StormRunner Load A fresh approach to cloud-based website load testing is proving more effective
More informationHierarchical expectation propagation for Bayesian aggregation of average data
Hierarchical expectation propagation for Bayesian aggregation of average data Andrew Gelman, Columbia University Sebastian Weber, Novartis also Bob Carpenter, Daniel Lee, Frédéric Bois, Aki Vehtari, and
More informationAdministering your Enterprise Geodatabase using Python. Jill Penney
Administering your Enterprise Geodatabase using Python Jill Penney Assumptions Basic knowledge of python Basic knowledge enterprise geodatabases and workflows You want code Please turn off or silence cell
More informationReaxys Pipeline Pilot Components Installation and User Guide
1 1 Reaxys Pipeline Pilot components for Pipeline Pilot 9.5 Reaxys Pipeline Pilot Components Installation and User Guide Version 1.0 2 Introduction The Reaxys and Reaxys Medicinal Chemistry Application
More informationRANDOM and REPEATED statements - How to Use Them to Model the Covariance Structure in Proc Mixed. Charlie Liu, Dachuang Cao, Peiqi Chen, Tony Zagar
Paper S02-2007 RANDOM and REPEATED statements - How to Use Them to Model the Covariance Structure in Proc Mixed Charlie Liu, Dachuang Cao, Peiqi Chen, Tony Zagar Eli Lilly & Company, Indianapolis, IN 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 18, 2012
More informationENHANCING SOLID-PHASE DISK EXTRACTION PERFORMANCE WITH DESIGN OF EXPERIMENTS
ENHANCING SOLID-PHASE DISK EXTRACTION PERFORMANCE WITH DESIGN OF EXPERIMENTS BY CRAIG A. PERMAN Product Development Specialist 3M Company Filtration Products Laboratory Building 60-1W-17 St. Paul, MN 55144
More informationWeatherHawk Weather Station Protocol
WeatherHawk Weather Station Protocol Purpose To log atmosphere data using a WeatherHawk TM weather station Overview A weather station is setup to measure and record atmospheric measurements at 15 minute
More informationThe Comprehensive Report
High-Throughput Screening 2002: New Strategies and Technologies The Comprehensive Report Presented by HighTech Business Decisions 346 Rheem Blvd., Suite 208, Moraga, CA 94556 Tel: (925) 631-0920 Fax: (925)
More informationDISPLAYING THE POISSON REGRESSION ANALYSIS
Chapter 17 Poisson Regression Chapter Table of Contents DISPLAYING THE POISSON REGRESSION ANALYSIS...264 ModelInformation...269 SummaryofFit...269 AnalysisofDeviance...269 TypeIII(Wald)Tests...269 MODIFYING
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 informationECS8020 ORGANIC ELEMENTAL ANALYZER CHNS-O Analyzer
ECS8020 ORGANIC ELEMENTAL ANALYZER CHNS-O Analyzer ECS 8020 CHNS-O Analyzer Organic Elemental Analysis ECS 8020 is a C-H-N-S-O Elemental Analyzer Model based on the Dumas combustion method. ECS 8020 is
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 informationPerforming response surface analysis using the SAS RSREG procedure
Paper DV02-2012 Performing response surface analysis using the SAS RSREG procedure Zhiwu Li, National Database Nursing Quality Indicator and the Department of Biostatistics, University of Kansas Medical
More informationExtensions of Cox Model for Non-Proportional Hazards Purpose
PhUSE 2013 Paper SP07 Extensions of Cox Model for Non-Proportional Hazards Purpose Jadwiga Borucka, PAREXEL, Warsaw, Poland ABSTRACT Cox proportional hazard model is one of the most common methods used
More informationPrincipal Component Analysis, A Powerful Scoring Technique
Principal Component Analysis, A Powerful Scoring Technique George C. J. Fernandez, University of Nevada - Reno, Reno NV 89557 ABSTRACT Data mining is a collection of analytical techniques to uncover new
More informationSAS/STAT 15.1 User s Guide The SEQDESIGN Procedure
SAS/STAT 15.1 User s Guide The SEQDESIGN Procedure This document is an individual chapter from SAS/STAT 15.1 User s Guide. The correct bibliographic citation for this manual is as follows: SAS Institute
More informationHow the mean changes depends on the other variable. Plots can show what s happening...
Chapter 8 (continued) Section 8.2: Interaction models An interaction model includes one or several cross-product terms. Example: two predictors Y i = β 0 + β 1 x i1 + β 2 x i2 + β 12 x i1 x i2 + ɛ i. How
More informationPopulation Approach. Pharmacokinetics. Analysis Approaches. Describing The Signal and the Noise
1 Population Approach Describing The Signal and the Noise 2 Pharmacometrics Science Pharmacology Models Pharmacokinetics Pharmacodynamics Experiment Dose Concentration Effect 3 Analysis Approaches ANOVA
More informationPharmacometrics : Nonlinear mixed effect models in Statistics. Department of Statistics Ewha Womans University Eun-Kyung Lee
Pharmacometrics : Nonlinear mixed effect models in Statistics Department of Statistics Ewha Womans University Eun-Kyung Lee 1 Clinical Trials Preclinical trial: animal study Phase I trial: First in human,
More informationQuestion. Hypothesis testing. Example. Answer: hypothesis. Test: true or not? Question. Average is not the mean! μ average. Random deviation or not?
Hypothesis testing Question Very frequently: what is the possible value of μ? Sample: we know only the average! μ average. Random deviation or not? Standard error: the measure of the random deviation.
More informationMETHODS FOR POPULATION PHARMACOKINETICS AND PHARMACODYNAMICS
METHODS FOR POPULATION PHARMACOKINETICS AND PHARMACODYNAMICS Emily Anne Colby A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment of the requirements
More informationGrouping and Annual Reporting. Agenda. Grouping of extension and other variations
Grouping and Annual Reporting Zahra Hanaizi Scientific Administrator Safety & Efficacy of Medicines 1 Agenda Grouping of type IA & Periodic Report Rules and examples Procedural aspects Grouping of extension
More informationParkServe. The Trust for Public Land s. Emmalee Dolfi Gabriel Patterson-King ESRI User Conference 2017
The Trust for Public Land s ParkServe Automating Large-Scale Data Transformation, QA/QC and Analysis Emmalee Dolfi Gabriel Patterson-King ESRI User Conference 2017 Our mission The Trust for Public Land
More informationDesign of discrete-event simulations
Design of discrete-event simulations Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/kurssit/tlt-2707/ OUTLINE: Discrete event simulation; Event advance design; Unit-time
More informationTutorial 4: Power and Sample Size for the Two-sample t-test with Unequal Variances
Tutorial 4: Power and Sample Size for the Two-sample t-test with Unequal Variances Preface Power is the probability that a study will reject the null hypothesis. The estimated probability is a function
More informationPaper Equivalence Tests. Fei Wang and John Amrhein, McDougall Scientific Ltd.
Paper 11683-2016 Equivalence Tests Fei Wang and John Amrhein, McDougall Scientific Ltd. ABSTRACT Motivated by the frequent need for equivalence tests in clinical trials, this paper provides insights into
More informationABSTRACT INTRODUCTION SUMMARY OF ANALYSIS. Paper
Paper 1891-2014 Using SAS Enterprise Miner to predict the Injury Risk involved in Car Accidents Prateek Khare, Oklahoma State University; Vandana Reddy, Oklahoma State University; Goutam Chakraborty, Oklahoma
More informationAnalysis of Longitudinal Data: Comparison between PROC GLM and PROC MIXED.
Analysis of Longitudinal Data: Comparison between PROC GLM and PROC MIXED. Maribeth Johnson, Medical College of Georgia, Augusta, GA ABSTRACT Longitudinal data refers to datasets with multiple measurements
More informationarxiv: v1 [stat.me] 14 Jan 2019
arxiv:1901.04443v1 [stat.me] 14 Jan 2019 An Approach to Statistical Process Control that is New, Nonparametric, Simple, and Powerful W.J. Conover, Texas Tech University, Lubbock, Texas V. G. Tercero-Gómez,Tecnológico
More informationLeverage Sparse Information in Predictive Modeling
Leverage Sparse Information in Predictive Modeling Liang Xie Countrywide Home Loans, Countrywide Bank, FSB August 29, 2008 Abstract This paper examines an innovative method to leverage information from
More informationTMDD Model Translated from NONMEM (NM- TRAN) to Phoenix NLME (PML)
TMDD Model Translated from NONMEM (NM- TRAN) to Phoenix NLME (PML) Phoenix Modeling Language School March 22, 2018 Loan Pham, Ph.D. Senior Pharmacokinetic Scientist Camargo Pharmaceutical Services 9825
More informationLecture 7 Time-dependent Covariates in Cox Regression
Lecture 7 Time-dependent Covariates in Cox Regression So far, we ve been considering the following Cox PH model: λ(t Z) = λ 0 (t) exp(β Z) = λ 0 (t) exp( β j Z j ) where β j is the parameter for the the
More informationYun W. Alelyunas, Mark D. Wrona, Russell J. Mortishire-Smith, Nick Tomczyk, and Paul D. Rainville Waters Corporation, Milford, MA, USA INTRODUCTION
Quantitation by High Resolution Mass Spectrometry: Using Target Enhancement and Tof-MRM to Achieve Femtogram-level On-column Sensitivity for Quantitation of Drugs in Human Plasma Yun W. Alelyunas, Mark
More informationGMM Logistic Regression with Time-Dependent Covariates and Feedback Processes in SAS TM
Paper 1025-2017 GMM Logistic Regression with Time-Dependent Covariates and Feedback Processes in SAS TM Kyle M. Irimata, Arizona State University; Jeffrey R. Wilson, Arizona State University ABSTRACT The
More informationIrreversible Inhibition Kinetics
1 Irreversible Inhibition Kinetics Automation and Simulation Petr Kuzmič, Ph.D. BioKin, Ltd. 1. Automate the determination of biochemical parameters 2. PK/PD simulations with multiple injections Irreversible
More informationNonlinear random effects mixture models: Maximum likelihood estimation via the EM algorithm
Computational Statistics & Data Analysis 51 (2007) 6614 6623 www.elsevier.com/locate/csda Nonlinear random effects mixture models: Maximum lielihood estimation via the EM algorithm Xiaoning Wang a, Alan
More informationTutorial 3: Power and Sample Size for the Two-sample t-test with Equal Variances. Acknowledgements:
Tutorial 3: Power and Sample Size for the Two-sample t-test with Equal Variances Anna E. Barón, Keith E. Muller, Sarah M. Kreidler, and Deborah H. Glueck Acknowledgements: The project was supported in
More informationA MACRO-DRIVEN FORECASTING SYSTEM FOR EVALUATING FORECAST MODEL PERFORMANCE
A MACRO-DRIVEN ING SYSTEM FOR EVALUATING MODEL PERFORMANCE Bryan Sellers Ross Laboratories INTRODUCTION A major problem of forecasting aside from obtaining accurate forecasts is choosing among a wide range
More informationNWS/AFWA/Navy Office: JAN NWS (primary) and other NWS (see report) Name of NWS/AFWA/Navy Researcher Preparing Report: Jeff Craven (Alan Gerard)
University of Louisiana at Monroe Name of University Researcher Preparing Report: Dr. Paul J. Croft NWS/AFWA/Navy Office: JAN NWS (primary) and other NWS (see report) Name of NWS/AFWA/Navy Researcher Preparing
More informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY (formerly the Examinations of the Institute of Statisticians) GRADUATE DIPLOMA, 2007
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY (formerly the Examinations of the Institute of Statisticians) GRADUATE DIPLOMA, 2007 Applied Statistics I Time Allowed: Three Hours Candidates should answer
More informationA Rapid Approach to the Confirmation of Drug Metabolites in Preclinical and Clinical Bioanalysis Studies
A Rapid Approach to the Confirmation of Drug Metabolites in Preclinical and Clinical Bioanalysis Studies APPLICATION BENEFITS Regulatory guidelines and recommendations place a greater emphasis on the detection
More informationDesktop GIS for Geotechnical Engineering
Desktop GIS for Geotechnical Engineering Satya Priya Deputy General Manager (Software) RMSI, A-7, Sector 16 NOIDA 201 301, UP, INDIA Tel: +91-120-2511102 Fax: +91-120-2510963 Email: Satya.Priya@rmsi.com
More informationAutomated, intelligent sample preparation: Integration of the ESI prepfast Auto-dilution System with the Thermo Scientific icap 7400 ICP-OES
TECHNICAL NOTE 43252 Automated, intelligent sample preparation: Integration of the ESI prepfast Auto-dilution System with the Thermo Scientific icap 7400 ICP-OES Keywords Auto-dilution, Intelligent dilution,
More informationNINE CHOICE SERIAL REACTION TIME TASK
instrumentation and software for research NINE CHOICE SERIAL REACTION TIME TASK MED-STATE NOTATION PROCEDURE SOF-700RA-8 USER S MANUAL DOC-025 Rev. 1.3 Copyright 2013 All Rights Reserved MED Associates
More informationIn Search of Desirable Compounds
In Search of Desirable Compounds Adrijo Chakraborty University of Georgia Email: adrijoc@uga.edu Abhyuday Mandal University of Georgia Email: amandal@stat.uga.edu Kjell Johnson Arbor Analytics, LLC Email:
More informationCharity Quick, Rho, Inc, Chapel Hill, NC Paul Nguyen, Rho, Inc, Chapel Hill, NC
PharmaSUG 2016 - Paper DS09 Prepare for Re-entry: Challenges and Solutions for Handling Re-screened Subjects in SDTM ABSTRACT Charity Quick, Rho, Inc, Chapel Hill, NC Paul Nguyen, Rho, Inc, Chapel Hill,
More informationPharmacokinetics Introduction to
Wikimedia Commons Korinna Creative Commons Attribution-ShareAlike 3. Unported Introduction to Population PK Pharmacokinetics Introduction to Introduction to Population PK Helmut Schütz BEBAC Workshop Bucarest,
More informationSTAT 5500/6500 Conditional Logistic Regression for Matched Pairs
STAT 5500/6500 Conditional Logistic Regression for Matched Pairs The data for the tutorial came from support.sas.com, The LOGISTIC Procedure: Conditional Logistic Regression for Matched Pairs Data :: SAS/STAT(R)
More informationModeling Effect Modification and Higher-Order Interactions: Novel Approach for Repeated Measures Design using the LSMESTIMATE Statement in SAS 9.
Paper 400-015 Modeling Effect Modification and Higher-Order Interactions: Novel Approach for Repeated Measures Design using the LSMESTIMATE Statement in SAS 9.4 Pronabesh DasMahapatra, MD, MPH, PatientsLikeMe
More informationLiang Li, PhD. MD Anderson
Liang Li, PhD Biostatistics @ MD Anderson Genetic Epidemiology Work Group Meeting, December 5, 2013 The Multiphase Optimization Strategy (MOST) An increasingly popular study design strategy in behavioral
More informationASSESSMENT. Industry Solutions Harness the Power of GIS for Property Assessment
ASSESSMENT Industry Solutions Harness the Power of GIS for Property Assessment Esri Canada has thousands of customers worldwide who are using the transforming power of GIS technology to collect, maintain,
More informationEvaluating Corridors. Jeff Jenness Dan Majka Paul Beier. CorridorDesigner Workshop Fall 2007
Evaluating Corridors Jeff Jenness Dan Majka Paul Beier CorridorDesigner Workshop Fall 2007 1 1 Evaluating Corridors Test Alternative Corridors Least cost methods always provide a best solution, even when
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 informationChapter 3 Multiple Regression Complete Example
Department of Quantitative Methods & Information Systems ECON 504 Chapter 3 Multiple Regression Complete Example Spring 2013 Dr. Mohammad Zainal Review Goals After completing this lecture, you should be
More informationBefore and After Models in Observational Research Using Random Slopes and Intercepts
Paper 3643-2015 Before and After Models in Observational Research Using Random Slopes and Intercepts David J. Pasta, ICON Clinical Research, San Francisco, CA ABSTRACT In observational data analyses, it
More informationRobust design in model-based analysis of longitudinal clinical data
Robust design in model-based analysis of longitudinal clinical data Giulia Lestini, Sebastian Ueckert, France Mentré IAME UMR 1137, INSERM, University Paris Diderot, France PODE, June 0 016 Context Optimal
More informationMS-C1620 Statistical inference
MS-C1620 Statistical inference 10 Linear regression III Joni Virta Department of Mathematics and Systems Analysis School of Science Aalto University Academic year 2018 2019 Period III - IV 1 / 32 Contents
More informationTutorial 2: Power and Sample Size for the Paired Sample t-test
Tutorial 2: Power and Sample Size for the Paired Sample t-test Preface Power is the probability that a study will reject the null hypothesis. The estimated probability is a function of sample size, variability,
More informationReal-time Systems: Scheduling Periodic Tasks
Real-time Systems: Scheduling Periodic Tasks Advanced Operating Systems Lecture 15 This work is licensed under the Creative Commons Attribution-NoDerivatives 4.0 International License. To view a copy of
More informationThe Cause-Effect Diagram can be used under these Circumstances:
What is Cause and Effect Graph Testing Technique - How to Design Test Cases With Example? Cause-Effect Graph graphically shows the connection between a given outcome and all issues that manipulate the
More informationResearch Article. Multiple Imputation of Missing Covariates in NONMEM and Evaluation of the Method s Sensitivity to η-shrinkage
The AAPS Journal, Vol. 15, No. 4, October 2013 ( # 2013) DOI: 10.1208/s12248-013-9508-0 Research Article Multiple Imputation of Missing Covariates in NONMEM and Evaluation of the Method s Sensitivity to
More informationNormalization of Peak Demand for an Electric Utility using PROC MODEL
Normalization of Peak Demand for an Electric Utility using PROC MODEL Mark Harris, Jeff Brown, and Mark Gilbert* Central and South West Services, Inc. Tulsa, Oklahoma Abstract This paper discusses the
More informationStatistical Machine Learning Theory. From Multi-class Classification to Structured Output Prediction. Hisashi Kashima.
http://goo.gl/xilnmn Course website KYOTO UNIVERSITY Statistical Machine Learning Theory From Multi-class Classification to Structured Output Prediction Hisashi Kashima kashima@i.kyoto-u.ac.jp DEPARTMENT
More informationNonlinear mixed-effects models using Stata
Nonlinear mixed-effects models using Stata Yulia Marchenko Executive Director of Statistics StataCorp LP 2017 German Stata Users Group meeting Yulia Marchenko (StataCorp) 1 / 48 Outline What is NLMEM?
More information