An effective approach for obtaining optimal sampling windows for population pharmacokinetic experiments
|
|
- Louise Barnett
- 5 years ago
- Views:
Transcription
1 An effective approach for obtaining optimal sampling windows for population pharmacokinetic experiments Kayode Ogungbenro and Leon Aarons Centre for Applied Pharmacokinetic Research School of Pharmacy and Pharmaceutical Sciences The University of Manchester, Manchester United Kingdom
2 Introduction Population pharmacokinetics involves collection of blood samples Sample collection at specific times may not be feasible less informative Size of the trial (Phase III elays in seeing the medical personnel Poor patient compliance with respect to dosing times More immediate medical procedure Sampling windows controlling sampling times Sampling within some time intervals Gives flexibility, informative data and satisfactory parameter estimation
3 Population PK modelling Level Individual level yi f ( θi, ti + ε i,..., ni, i,..., N Level Population level θ i g( θ, bi ~ N(0, Ω Parameters Ψ ( 0, σ σ f ( θ, t ε ~ N + i i i b i ω Ω ω. [ θ,..., θ, ω,..., ω, ω,..., ω,..., ω, σ σ ] p p p pp, ω ω... ω pp Sampling Times ξ i [ t,..., t ] i in i 3
4 4 Population Fisher information Matrix Log-likelihood of observations Expressions for PFIM based on mixed effects modelling by maximum likelihood method Ψ Ψ Ψ Ψ T i i y l E F ; (, ( ξ Approximations based on linearization of the model Ψ Ψ + Ψ, ( V V V V tr J V J F z a T z a az ξ i Population design Ξ Q Q N N,...,,..., ξ ξ Ψ Ψ Ψ Ξ Q q q q Q q N i i F N F F q q q, (, (, ( ξ ξ
5 5 Optimal esign -Efficiency, ( max arg Ξ Ψ Ξ F χ / dim(, (, (,, ( Ψ Ψ Ξ Ψ Ξ Ξ Ξ Ψ F F eff Optimality criteria - -optimality Population -optimal design
6 Sampling windows determination - approaches uffull et al. (Pharmaceutical Research 00 Sampling windows design obtained as marginal windows associated with each time points varying one sampling time at a time until determinant is reduced by 5% Graham and Aarons (Statistics in Medicine 006 Two stage approach and a quadratic loss function Sampling windows design result in specified loss of efficiency compared to fixed -optimal time points Patan and Bogacka (Advances in Model-Oriented esign and analysis 007 Based on equivalence theorem for -optimal continuous designs using the variance function 6
7 Sampling windows determination - approaches Summary of attributes Attribute uffull G&A P&B Parameter Sensitivities Yes No Yes Assessment of Sampling Windows Efficiency No Yes Yes Exact Yes Yes No esign Continuous Yes Yes Yes 7
8 Sampling windows determination new approach Main features Very efficient and effective Can be applied to different types of design Exact Continuous Reflects parameter sensitivities Less flexibility (narrow window for high parameter sensitivities (important to sample close to optimal time point More flexibility (wide window for low parameter sensitivities (less important to sample close to optimal time point The efficiency of the sampling windows design can be assessed ointly 8
9 9 Sampling windows determination Sampling windows population design Fixed time population -optimal design Ξ L n L U n U W t t t t,...,,..., L U t t t t δ δ + [ ] δ,...,δ n Δ [ ] n t t N,..., Ξ Ξ ξ
10 Sampling windows determination Efficiency functions conditional and oint Uniform or loguniform distribution eff W / dim( Ψ E F( Ψ, Ξ ( δ W ( Ψ, Ξ ( δ / dim( Ψ F( Ψ, Ξ eff E F( Ψ, Ξ / dim( Ψ W ( Ψ, Ξ ( Δ / dim( Ψ W F( Ψ, Ξ ( Δ 0
11 Sampling windows determination 3 stage approach Stage : Optimisation of fixed -optimal time points Stage : Assuming a distribution uniform or loguniform efine a target mean efficiency level, eff 0 Optimise one window length at a time using a quadratic function δ arg min δ δ Δ [( ] W eff ( Ψ, Ξ ( eff 0
12 Sampling windows determination Stage 3: Evaluate the efficiency of the oint sampling windows eff ( Ψ, Ξ W ( Δ Check if eff is greater than eff 0 if not (ideally reduce the window lengths by equal percentage, % and obtain a new vector for the lengths of the sampling windows δ δ ( new ( old ( old ( δ * 0.0 Repeat Stage 3 until the required efficiency level is obtained Take the last vector of sampling windows lengths as the optimal sampling windows lengths
13 Sampling windows determination Apart from mean efficiency level, percentiles can also be used - mean efficiency can produce variation in the realized design depending on distribution of samples within the windows Sampling at or close to boundaries less efficient design Sampling at or close to fixed -optimal time points more efficient design Constraints during optimisation especially if any of the fixed -optimal time points is at the boundary or near the boundary Ensure sampling windows do not extend outside the design space (sampling at negative time points 3
14 Sampling windows determination Example Individual continuous design One compartment IV bolus model (proportional residual Ψ y ose V ( Cl / V t esign region 0 and 4 hrs e [ Cl, V, σ ] [ 3,30,0.04] ose 450mg Sampling windows P&B and new approach 95% and 90% efficiency levels Mean, 5th and 0th percentile efficiency criteria Uniform and loguniform sample distributions 4
15 Sampling windows determination Example Fixed -optimal time points 0 and 4 hr (0.5,0.5 Efficiency level Criteria P&B Uniform New Approach P&B Loguniform New Approach Mean % 0 th th Mean % 0 th th
16 Sampling windows determination Example Individual continuous design One compartment first order absorption model at steady state (the data and the model are log transformed log( y Ψ Fka log Vka Cl F, mg, τ hr e e / V kat e e Clτ / V kaτ [ Cl, V, ka, σ ] [.55,00,.08,0.05] Clt + ε ( σ esign region 0. and hr Sampling windows P&B and new approach 95% and 90% efficiency levels Mean, 5th and 0th percentile efficiency criteria Uniform and loguniform sample distributions 6
17 Sampling windows determination Example Fixed -optimal time points 0.9,.46, hr (0.33,0.33,0.33 Efficiency level Criteria P&B Uniform New Approach P&B Loguniform New Approach Mean % 0 th th Mean % 0 th th
18 Sampling windows determination Example (95% Efficiency, Mean criteria and Uniform distribution Fixed -optimal time points 0.9,.46, hr (0.33,0.33,0.33 sensitivity Equivalence theorem - sensitivity time (min 8
19 Sampling windows determination Example 3 Population exact design One compartment first order absorption model at steady state (the data and the model are log transformed Ψ Cliti / Vi kaiti Fka i e e log( y i log ε ( σ Cl τ / V ka τ + i i i i Vikai Cl i e e [ Cl, V, ka, ω, ω, ω, σ ] [.55,00,.08,0.09,0.09,0.09,0.05] Cl F, mg, hr V ka τ, b ] [ 0,] [a subects, group and 3 time points Sampling windows G&A, uffull and new approach 95% and 90% efficiency levels Mean and 0th percentile efficiency criteria Uniform and loguniform sample distributions 9
20 Sampling windows determination Example (Uniform distribution Optimal design (fixed time 0.30,.9 and hrs (det Efficiency Criterion G&A uffull New approach Mean (0.8*.6-.7 ( ( ( ( ( ( ( (4.0 det** [0.90] det 4.3 [0.74] det 7.9 [0.9] 90% 0% ( ( ( ( ( ( ( ( (3.9 det [0.90] det 48.6 [0.78] det 7.36 [0.90] Mean ( ( ( ( ( ( ( ( (.78 det [0.95] det [0.84] det [0.95] 95% 0% ( ( ( ( ( ( ( ( (.0 det [0.95] det 70.7 [0.90] det 8.7 [0.95] *Window half lengths in parentheses, **det is the normalised determinant 0
21 Sampling windows determination Example (Loguniform distribution Optimal design (fixed time 0.30,.9 and hrs (det Efficiency Criterion G&A uffull New approach 90% Mean 0% (0.54* ( (0.54 det** 7.43 [0.90] ( ( (0.40 det 70.4 [0.90] ( ( (0.7 det 9.9 [0.68] ( ( (0.44 det 45. [ ( ( (0.40 det 7.5 [0.9] ( ( (0.30 det [0.9] Mean ( ( ( ( ( ( ( ( (0.5 det 80.8 [0.95] det [0.83] det 80.5 [0.95] 95% 0% ( ( ( ( ( ( ( ( (0.7 det 80. [0.95] det [0.89] det [0.95] *Window half lengths (on the logarithmic scale in parentheses, **det is the normalised determinant
22 Sampling windows determination Example (uniform distribution, 90% efficiency level and mean criteria G&A uffull New ( ( ( ( ( ( ( ( (4.0
23 Sampling windows determination Example (simulations, RE 3 Fixed Time 3 G&A RE RE Cl V ka λ Cl λ V λ ka σ Cl V ka λ Cl λ V λ ka σ 3 uffull 3 New RE RE Cl V ka λ Cl λ V λ ka σ Cl V ka λ Cl λ V λ ka σ 3
24 Sampling windows determination Conclusion Sampling windows provide adequate flexibility (controlled flexibility for sample collection The new approach is efficient and reflects parameter sensitivities This approach can be applied to both exact and continuous design as well as multiresponse designs Choice of efficiency level, criteria function and parameter distribution must be balanced against other design properties Efficient population PK experiment can provide improved parameter estimates and can help to reduce cost and time 4
Robust 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 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 informationDESIGN EVALUATION AND OPTIMISATION IN CROSSOVER PHARMACOKINETIC STUDIES ANALYSED BY NONLINEAR MIXED EFFECTS MODELS
DESIGN EVALUATION AND OPTIMISATION IN CROSSOVER PHARMACOKINETIC STUDIES ANALYSED BY NONLINEAR MIXED EFFECTS MODELS Thu Thuy Nguyen, Caroline Bazzoli, France Mentré UMR 738 INSERM - University Paris Diderot,
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 informationADAPTIVE EXPERIMENTAL DESIGNS. Maciej Patan and Barbara Bogacka. University of Zielona Góra, Poland and Queen Mary, University of London
ADAPTIVE EXPERIMENTAL DESIGNS FOR SIMULTANEOUS PK AND DOSE-SELECTION STUDIES IN PHASE I CLINICAL TRIALS Maciej Patan and Barbara Bogacka University of Zielona Góra, Poland and Queen Mary, University of
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 informationEvaluation of the Fisher information matrix in nonlinear mixed eect models without linearization
Evaluation of the Fisher information matrix in nonlinear mixed eect models without linearization Sebastian Ueckert, Marie-Karelle Riviere, and France Mentré INSERM, IAME, UMR 1137, F-75018 Paris, France;
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 informationHandling parameter constraints in complex/mechanistic population pharmacokinetic models
Handling parameter constraints in complex/mechanistic population pharmacokinetic models An application of the multivariate logistic-normal distribution Nikolaos Tsamandouras Centre for Applied Pharmacokinetic
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 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 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 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 informationMinimum Message Length Analysis of the Behrens Fisher Problem
Analysis of the Behrens Fisher Problem Enes Makalic and Daniel F Schmidt Centre for MEGA Epidemiology The University of Melbourne Solomonoff 85th Memorial Conference, 2011 Outline Introduction 1 Introduction
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 informationAP-Optimum Designs for Minimizing the Average Variance and Probability-Based Optimality
AP-Optimum Designs for Minimizing the Average Variance and Probability-Based Optimality Authors: N. M. Kilany Faculty of Science, Menoufia University Menoufia, Egypt. (neveenkilany@hotmail.com) W. A. Hassanein
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 informationLecture Stat Information Criterion
Lecture Stat 461-561 Information Criterion Arnaud Doucet February 2008 Arnaud Doucet () February 2008 1 / 34 Review of Maximum Likelihood Approach We have data X i i.i.d. g (x). We model the distribution
More informationDose-finding for Multi-drug Combinations
Outline Background Methods Results Conclusions Multiple-agent Trials In trials combining more than one drug, monotonicity assumption may not hold for every dose The ordering between toxicity probabilities
More informationDose-finding for Multi-drug Combinations
September 16, 2011 Outline Background Methods Results Conclusions Multiple-agent Trials In trials combining more than one drug, monotonicity assumption may not hold for every dose The ordering between
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 informationTime-Dependent Conduction :
Time-Dependent Conduction : The Lumped Capacitance Method Chapter Five Sections 5.1 thru 5.3 Transient Conduction A heat transfer process for which the temperature varies with time, as well as location
More informationPenalized D-optimal design for dose finding
Penalized for dose finding Luc Pronzato Laboratoire I3S, CNRS-Univ. Nice Sophia Antipolis, France 2/ 35 Outline 3/ 35 Bernoulli-type experiments: Y i {0,1} (success or failure) η(x,θ) = Prob(Y i = 1 x
More informationGARCH Models Estimation and Inference. Eduardo Rossi University of Pavia
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationOptimum Designs for the Equality of Parameters in Enzyme Inhibition Kinetic Models
Optimum Designs for the Equality of Parameters in Enzyme Inhibition Kinetic Models Anthony C. Atkinson, Department of Statistics, London School of Economics, London WC2A 2AE, UK and Barbara Bogacka, School
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 informationCentral Limit Theorem ( 5.3)
Central Limit Theorem ( 5.3) Let X 1, X 2,... be a sequence of independent random variables, each having n mean µ and variance σ 2. Then the distribution of the partial sum S n = X i i=1 becomes approximately
More informationOn construction of constrained optimum designs
On construction of constrained optimum designs Institute of Control and Computation Engineering University of Zielona Góra, Poland DEMA2008, Cambridge, 15 August 2008 Numerical algorithms to construct
More informationarxiv: v1 [math.st] 22 Dec 2018
Optimal Designs for Prediction in Two Treatment Groups Rom Coefficient Regression Models Maryna Prus Otto-von-Guericke University Magdeburg, Institute for Mathematical Stochastics, PF 4, D-396 Magdeburg,
More informationImproving a safety of the Continual Reassessment Method via a modified allocation rule
Improving a safety of the Continual Reassessment Method via a modified allocation rule Pavel Mozgunov, Thomas Jaki Medical and Pharmaceutical Statistics Research Unit, Department of Mathematics and Statistics,
More informationDesign of HIV Dynamic Experiments: A Case Study
Design of HIV Dynamic Experiments: A Case Study Cong Han Department of Biostatistics University of Washington Kathryn Chaloner Department of Biostatistics University of Iowa Nonlinear Mixed-Effects Models
More informationThe Hybrid Likelihood: Combining Parametric and Empirical Likelihoods
1/25 The Hybrid Likelihood: Combining Parametric and Empirical Likelihoods Nils Lid Hjort (with Ingrid Van Keilegom and Ian McKeague) Department of Mathematics, University of Oslo Building Bridges (at
More informationOptimum designs for model. discrimination and estimation. in Binary Response Models
Optimum designs for model discrimination and estimation in Binary Response Models by Wei-Shan Hsieh Advisor Mong-Na Lo Huang Department of Applied Mathematics National Sun Yat-sen University Kaohsiung,
More informationA comparison of estimation methods in nonlinear mixed effects models using a blind analysis
A comparison of estimation methods in nonlinear mixed effects models using a blind analysis Pascal Girard, PhD EA 3738 CTO, INSERM, University Claude Bernard Lyon I, Lyon France Mentré, PhD, MD Dept Biostatistics,
More informationAccurate Maximum Likelihood Estimation for Parametric Population Analysis. Bob Leary UCSD/SDSC and LAPK, USC School of Medicine
Accurate Maximum Likelihood Estimation for Parametric Population Analysis Bob Leary UCSD/SDSC and LAPK, USC School of Medicine Why use parametric maximum likelihood estimators? Consistency: θˆ θ as N ML
More informationNONLINEAR MODELS IN MULTIVARIATE POPULATION BIOEQUIVALENCE TESTING
Virginia Commonwealth University VCU Scholars Compass Theses and Dissertations Graduate School NONLINEAR MODELS IN MULTIVARIATE POPULATION BIOEQUIVALENCE TESTING Bassam Dahman Virginia Commonwealth University
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 informationInformation geometry for bivariate distribution control
Information geometry for bivariate distribution control C.T.J.Dodson + Hong Wang Mathematics + Control Systems Centre, University of Manchester Institute of Science and Technology Optimal control of stochastic
More informationWald s theorem and the Asimov data set
Wald s theorem and the Asimov data set Eilam Gross & Ofer Vitells ATLAS statistics forum, Dec. 009 1 Outline We have previously guessed that the median significance of many toy MC experiments could be
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 informationOn the efficiency of two-stage adaptive designs
On the efficiency of two-stage adaptive designs Björn Bornkamp (Novartis Pharma AG) Based on: Dette, H., Bornkamp, B. and Bretz F. (2010): On the efficiency of adaptive designs www.statistik.tu-dortmund.de/sfb823-dp2010.html
More informationAccounting for Population Uncertainty in Covariance Structure Analysis
Accounting for Population Uncertainty in Structure Analysis Boston College May 21, 2013 Joint work with: Michael W. Browne The Ohio State University matrix among observed variables are usually implied
More informationUNIVERSITY OF TORONTO Faculty of Arts and Science
UNIVERSITY OF TORONTO Faculty of Arts and Science December 2013 Final Examination STA442H1F/2101HF Methods of Applied Statistics Jerry Brunner Duration - 3 hours Aids: Calculator Model(s): Any calculator
More informationDuke University. Duke Biostatistics and Bioinformatics (B&B) Working Paper Series. Randomized Phase II Clinical Trials using Fisher s Exact Test
Duke University Duke Biostatistics and Bioinformatics (B&B) Working Paper Series Year 2010 Paper 7 Randomized Phase II Clinical Trials using Fisher s Exact Test Sin-Ho Jung sinho.jung@duke.edu This working
More informationDesigns for Clinical Trials
Designs for Clinical Trials Chapter 5 Reading Instructions 5.: Introduction 5.: Parallel Group Designs (read) 5.3: Cluster Randomized Designs (less important) 5.4: Crossover Designs (read+copies) 5.5:
More informationSome Assorted Formulae. Some confidence intervals: σ n. x ± z α/2. x ± t n 1;α/2 n. ˆp(1 ˆp) ˆp ± z α/2 n. χ 2 n 1;1 α/2. n 1;α/2
STA 248 H1S MIDTERM TEST February 26, 2008 SURNAME: SOLUTIONS GIVEN NAME: STUDENT NUMBER: INSTRUCTIONS: Time: 1 hour and 50 minutes Aids allowed: calculator Tables of the standard normal, t and chi-square
More informationGARCH Models Estimation and Inference
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 1 Likelihood function The procedure most often used in estimating θ 0 in
More informationY1 Double Maths Assignment λ (lambda) Exam Paper to do Core 1 Solomon C on the VLE. Drill
α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Nature is an infinite sphere of which the centre is everywhere and the circumference nowhere Blaise Pascal Y Double Maths Assignment λ (lambda) Tracking
More informationA HEURISTIC ALGORITHM FOR CAPACITY SIZING OF FIXED ROUTING MULTISERVICE ERLANG LOSS NETWORKS
A HEURISTIC ALGORITHM FOR CAPACITY SIZING OF FIXED ROUTING MULTISERVICE ERLANG LOSS NETWORKS Zdenko Vrdolak and Mladen Kos University of Zagreb, Faculty of Electrical Engineering and Computing, Unska 3,
More informationSGN Advanced Signal Processing: Lecture 8 Parameter estimation for AR and MA models. Model order selection
SG 21006 Advanced Signal Processing: Lecture 8 Parameter estimation for AR and MA models. Model order selection Ioan Tabus Department of Signal Processing Tampere University of Technology Finland 1 / 28
More informationLecture 32: Asymptotic confidence sets and likelihoods
Lecture 32: Asymptotic confidence sets and likelihoods Asymptotic criterion In some problems, especially in nonparametric problems, it is difficult to find a reasonable confidence set with a given confidence
More informationEfficient computation of Bayesian optimal discriminating designs
Efficient computation of Bayesian optimal discriminating designs Holger Dette Ruhr-Universität Bochum Fakultät für Mathematik 44780 Bochum, Germany e-mail: holger.dette@rub.de Roman Guchenko, Viatcheslav
More informationFace Recognition and Biometric Systems
The Eigenfaces method Plan of the lecture Principal Components Analysis main idea Feature extraction by PCA face recognition Eigenfaces training feature extraction Literature M.A.Turk, A.P.Pentland Face
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 informationRerandomization to Balance Covariates
Rerandomization to Balance Covariates Kari Lock Morgan Department of Statistics Penn State University Joint work with Don Rubin University of Minnesota Biostatistics 4/27/16 The Gold Standard Randomized
More informationSwitching Regime Estimation
Switching Regime Estimation Series de Tiempo BIrkbeck March 2013 Martin Sola (FE) Markov Switching models 01/13 1 / 52 The economy (the time series) often behaves very different in periods such as booms
More informationDiagnostics can identify two possible areas of failure of assumptions when fitting linear models.
1 Transformations 1.1 Introduction Diagnostics can identify two possible areas of failure of assumptions when fitting linear models. (i) lack of Normality (ii) heterogeneity of variances It is important
More informationGroup Sequential Tests for Delayed Responses. Christopher Jennison. Lisa Hampson. Workshop on Special Topics on Sequential Methodology
Group Sequential Tests for Delayed Responses Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj Lisa Hampson Department of Mathematics and Statistics,
More informationMAS3301 Bayesian Statistics Problems 5 and Solutions
MAS3301 Bayesian Statistics Problems 5 and Solutions Semester 008-9 Problems 5 1. (Some of this question is also in Problems 4). I recorded the attendance of students at tutorials for a module. Suppose
More informationRecursive Estimation and Order Determination of Space-Time Autoregressive Processes
Recursive Estimation and Order Determination of Space-Time Autoregressive Processes Ana Monica C. Antunes, Barry G. Quinn & T. Subba Rao First version: 23 September 2005 Research Report No. 4, 2005, Probability
More informationGeneralizing the MCPMod methodology beyond normal, independent data
Generalizing the MCPMod methodology beyond normal, independent data José Pinheiro Joint work with Frank Bretz and Björn Bornkamp Novartis AG Trends and Innovations in Clinical Trial Statistics Conference
More informationNonparametric Identification of a Binary Random Factor in Cross Section Data - Supplemental Appendix
Nonparametric Identification of a Binary Random Factor in Cross Section Data - Supplemental Appendix Yingying Dong and Arthur Lewbel California State University Fullerton and Boston College July 2010 Abstract
More informationGroup Sequential Designs: Theory, Computation and Optimisation
Group Sequential Designs: Theory, Computation and Optimisation Christopher Jennison Department of Mathematical Sciences, University of Bath, UK http://people.bath.ac.uk/mascj 8th International Conference
More informationUnit 9: Inferences for Proportions and Count Data
Unit 9: Inferences for Proportions and Count Data Statistics 571: Statistical Methods Ramón V. León 1/15/008 Unit 9 - Stat 571 - Ramón V. León 1 Large Sample Confidence Interval for Proportion ( pˆ p)
More information1. Describing patient variability, 2. Minimizing patient variability, 3. Describing and Minimizing Intraindividual
1. Describing patient variability, 2. Minimizing patient variability, 3. Describing and Minimizing Intraindividual Variability. 4. Improving the reporting of lab assay results. Roger Jelliffe MD, Alan
More informationD-Optimal Designs for Second-Order Response Surface Models with Qualitative Factors
Journal of Data Science 920), 39-53 D-Optimal Designs for Second-Order Response Surface Models with Qualitative Factors Chuan-Pin Lee and Mong-Na Lo Huang National Sun Yat-sen University Abstract: Central
More informationUsing cold vapor generation atomic absorption to determine mercury impurities in pharmaceutical products
APPLICATION NOTE Using cold vapor generation atomic absorption to determine mercury impurities in pharmaceutical products AN43381 Dr. Jianfeng Cui, Thermo Fisher Scientific, Bremen, Germany Key Words al
More informationMasters Comprehensive Examination Department of Statistics, University of Florida
Masters Comprehensive Examination Department of Statistics, University of Florida May 6, 003, 8:00 am - :00 noon Instructions: You have four hours to answer questions in this examination You must show
More informationBIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation
BIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation Yujin Chung November 29th, 2016 Fall 2016 Yujin Chung Lec13: MLE Fall 2016 1/24 Previous Parametric tests Mean comparisons (normality assumption)
More informationThreshold estimation in marginal modelling of spatially-dependent non-stationary extremes
Threshold estimation in marginal modelling of spatially-dependent non-stationary extremes Philip Jonathan Shell Technology Centre Thornton, Chester philip.jonathan@shell.com Paul Northrop University College
More informationOPTIMAL DESIGN INPUTS FOR EXPERIMENTAL CHAPTER 17. Organization of chapter in ISSO. Background. Linear models
CHAPTER 17 Slides for Introduction to Stochastic Search and Optimization (ISSO)by J. C. Spall OPTIMAL DESIGN FOR EXPERIMENTAL INPUTS Organization of chapter in ISSO Background Motivation Finite sample
More informationPh.D. Qualifying Exam Friday Saturday, January 6 7, 2017
Ph.D. Qualifying Exam Friday Saturday, January 6 7, 2017 Put your solution to each problem on a separate sheet of paper. Problem 1. (5106) Let X 1, X 2,, X n be a sequence of i.i.d. observations from a
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 informationWell-developed and understood properties
1 INTRODUCTION TO LINEAR MODELS 1 THE CLASSICAL LINEAR MODEL Most commonly used statistical models Flexible models Well-developed and understood properties Ease of interpretation Building block for more
More informationThe decomposition of inequality and poverty
The decomposition of inequality and poverty The decomposition of the FGT index The FGT poverty index for a population composed of K groups can be written as follows: P(z; a) f(k)p(k; z; a) K k where P(k;
More informationCTP-656 Tablet Confirmed Superiority Of Pharmacokinetic Profile Relative To Kalydeco in Phase I Clinical Studies
Tablet Confirmed Superiority Of Pharmacokinetic Profile Relative To Kalydeco in Phase I Clinical Studies 39 th European Cystic Fibrosis Conference 8-11 June 2016, Basel, Switzerland 2014 Concert Pharmaceuticals,
More informationACCOUNTING FOR INPUT-MODEL AND INPUT-PARAMETER UNCERTAINTIES IN SIMULATION. <www.ie.ncsu.edu/jwilson> May 22, 2006
ACCOUNTING FOR INPUT-MODEL AND INPUT-PARAMETER UNCERTAINTIES IN SIMULATION Slide 1 Faker Zouaoui Sabre Holdings James R. Wilson NC State University May, 006 Slide From American
More informationUsing SAS Proc Power to Perform Model-based Power Analysis for Clinical Pharmacology Studies Peng Sun, Merck & Co., Inc.
PharmaSUG00 - Paper SP05 Using SAS Proc Power to Perform Model-based Power Analysis for Clinical Pharmacology Studies Peng Sun, Merck & Co., Inc., North Wales, PA ABSRAC In this paper, we demonstrate that
More informationThe requirements of a plant may be expressed in terms of (a) settling time (b) damping ratio (c) peak overshoot --- in time domain
Compensators To improve the performance of a given plant or system G f(s) it may be necessary to use a compensator or controller G c(s). Compensator Plant G c (s) G f (s) The requirements of a plant may
More informationTesting for bioequivalence of highly variable drugs from TR-RT crossover designs with heterogeneous residual variances
Testing for bioequivalence of highly variable drugs from T-T crossover designs with heterogeneous residual variances Christopher I. Vahl, PhD Department of Statistics Kansas State University Qing Kang,
More informationGeneralized Method of Moments (GMM) Estimation
Econometrics 2 Fall 2004 Generalized Method of Moments (GMM) Estimation Heino Bohn Nielsen of29 Outline of the Lecture () Introduction. (2) Moment conditions and methods of moments (MM) estimation. Ordinary
More informationInterval Estimation for the Ratio and Difference of Two Lognormal Means
UW Biostatistics Working Paper Series 12-7-2005 Interval Estimation for the Ratio and Difference of Two Lognormal Means Yea-Hung Chen University of Washington, yeahung@u.washington.edu Xiao-Hua Zhou University
More informationA Note on Cochran Test for Homogeneity in Two Ways ANOVA and Meta-Analysis
Open Journal of Statistics, 05, 5, 787-796 Published Online December 05 in SciRes. http://www.scirp.org/journal/ojs http://dx.doi.org/0.436/ojs.05.57078 A Note on Cochran Test for Homogeneity in Two Ways
More informationInternational Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: , ISSN(Online): Vol.9, No.9, pp , 2016
International Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: 0974-4304, ISSN(Online): 2455-9563 Vol.9, No.9, pp 477-487, 2016 Modeling of Multi-Response Longitudinal DataUsing Linear Mixed Modelin
More informationOn the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models
On the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models Thomas Kneib Department of Mathematics Carl von Ossietzky University Oldenburg Sonja Greven Department of
More informationMultivariate Statistics
Multivariate Statistics Chapter 2: Multivariate distributions and inference Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2016/2017 Master in Mathematical
More informationE-optimal approximate block designs for treatment-control comparisons
E-optimal approximate block designs for treatment-control comparisons Samuel Rosa 1 1 Faculty of Mathematics, Physics and Informatics, Comenius University in Bratislava, Slovakia We study E-optimal block
More informationNeighbourhoods of Randomness and Independence
Neighbourhoods of Randomness and Independence C.T.J. Dodson School of Mathematics, Manchester University Augment information geometric measures in spaces of distributions, via explicit geometric representations
More informationSTAT Advanced Bayesian Inference
1 / 32 STAT 625 - Advanced Bayesian Inference Meng Li Department of Statistics Jan 23, 218 The Dirichlet distribution 2 / 32 θ Dirichlet(a 1,...,a k ) with density p(θ 1,θ 2,...,θ k ) = k j=1 Γ(a j) Γ(
More information0, otherwise. U = Y 1 Y 2 Hint: Use either the method of distribution functions or a bivariate transformation. (b) Find E(U).
1. Suppose Y U(0, 2) so that the probability density function (pdf) of Y is 1 2, 0 < y < 2 (a) Find the pdf of U = Y 4 + 1. Make sure to note the support. (c) Suppose Y 1, Y 2,..., Y n is an iid sample
More informationModularity of Galois representations. imaginary quadratic fields
over imaginary quadratic fields Krzysztof Klosin (joint with T. Berger) City University of New York December 3, 2011 Notation F =im. quadr. field, p # Cl F d K, fix p p Notation F =im. quadr. field, p
More informationSupport Vector Machines
Support Vector Machines Sridhar Mahadevan mahadeva@cs.umass.edu University of Massachusetts Sridhar Mahadevan: CMPSCI 689 p. 1/32 Margin Classifiers margin b = 0 Sridhar Mahadevan: CMPSCI 689 p.
More informationOptimal input design for nonlinear dynamical systems: a graph-theory approach
Optimal input design for nonlinear dynamical systems: a graph-theory approach Patricio E. Valenzuela Department of Automatic Control and ACCESS Linnaeus Centre KTH Royal Institute of Technology, Stockholm,
More informationModel comparison: Deviance-based approaches
Model comparison: Deviance-based approaches Patrick Breheny February 19 Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/23 Model comparison Thus far, we have looked at residuals in a fairly
More informationPSI Journal Club March 10 th, 2016
PSI Journal Club March 1 th, 1 The analysis of incontinence episodes and other count data in patients with Overactive Bladder (OAB) by Poisson and negative binomial regression Martina R, Kay R, van Maanen
More informationEngineering Part IIB: Module 4F10 Statistical Pattern Processing Lecture 5: Single Layer Perceptrons & Estimating Linear Classifiers
Engineering Part IIB: Module 4F0 Statistical Pattern Processing Lecture 5: Single Layer Perceptrons & Estimating Linear Classifiers Phil Woodland: pcw@eng.cam.ac.uk Michaelmas 202 Engineering Part IIB:
More informationThis midterm covers Chapters 6 and 7 in WMS (and the notes). The following problems are stratified by chapter.
This midterm covers Chapters 6 and 7 in WMS (and the notes). The following problems are stratified by chapter. Chapter 6 Problems 1. Suppose that Y U(0, 2) so that the probability density function (pdf)
More informationWeb-based Supplementary Material for A Two-Part Joint. Model for the Analysis of Survival and Longitudinal Binary. Data with excess Zeros
Web-based Supplementary Material for A Two-Part Joint Model for the Analysis of Survival and Longitudinal Binary Data with excess Zeros Dimitris Rizopoulos, 1 Geert Verbeke, 1 Emmanuel Lesaffre 1 and Yves
More informationOptimal Experiment Design for the estimation of moisture material properties
Optimal Experiment Design for the estimation of moisture material properties Julien Berger in collaboration with: T. Busser, D. Dutykh, N. Mendes, H. Le Meur Journée SFT, Paris, Mai 2 nd, 218 Context Comparison
More informationLecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 6 Classical Control Overview IV Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lead Lag Compensator Design Dr. Radhakant Padhi Asst.
More information