Copyright 2015 Quintiles

Size: px
Start display at page:

Download "Copyright 2015 Quintiles"

Transcription

1 fficiency of Ranomize Concentration-Controlle Trials Relative to Ranomize Dose-Controlle Trials, an Application to Personalize Dosing Trials Russell Reeve, PhD Quintiles, Inc. Copyright 2015 Quintiles

2 RCCT vs RDCT Nomenclature an setup Ranomize ose-controlle trial Patients are ranomize to ose groups Drug from ose istributes into the bloostream, yieling exposure (concentration) Drug in the bloostream fins its way to the site of action Drug at site of action Ranomize concentration-controlle trial Patients are ranomize to concentration groups Drug ose is calculate from the concentration target Drug from ose istributes into the bloostream, yieling exposure (concentration)» Hopefully this is close to the target Drug in the bloostream fins its way to the site of action Drug at site of action Concentration can be summarize by AUC, C max, C min, mean concentration, etc. 2

3 Assume linearity Base on power moel C = concentration (or summary statistic) Base on power moel, log C = log + is occasion-to-occasion (but within subject) variability, with Var = 2 0 an 1 may vary subject to subject, but constant within subject We will assume Var 0 = 2, an Var 1 = 0 For ease of exposition, also assume ose proportionality (i.e., 1 = 1) Challenge is to estimate 0 for each patient Algorithm First ose chosen at ranom from oses that achieve expecte exposure on average Observe exposure Use that to ajust ose to bring onto target 3

4 ffect of Reucing PK Variability on PD Separation What we are trying to achieve 4

5 ffect of Reucing PK Variability on PD Separation What we are trying to achieve 5

6 Bayesian Dose Ajustment Assume prior 0 ~ N(B, 2 ) The posterior at step i then has mean an variance 6

7 Convergence Converges with quaratic spee 7

8 Definitions Moel for effect Hill moel Y = a + q, where q = {1 + exp[( m)]} 1 Hypotheses: H 0 : {Y} 0 vs H A : {Y} follow Hill moel Let test T 1 require n 1 to achieve power p, an test T 2 require n 2 to achieve the same power p. Then we efine efficiency to be eff(t 1, T 2 ) = n 2 /n 1 (Lehmann 1986, p. 321) 8

9 RCCT Theorem Unerlies justification Let PK parameter x be efine as x = () + + form some function of the ose, patient-specific intercept, an resiual variability ~ N(0, 2 ) The variance 2 will vary epening on the esign, with RDCT being less variable than RCCT. We note that RCCT RDCT. Theorem 1 (Reeve 1994): Funamental RCCT Theorem eff(rcct, RDCT) 1 with equality if an only if RCCT = RDCT. Theorem 2: Consistency of Bayesian estimates (cf. Ghosal Theorem 4 quoting Ghosal, Ghosal, Samanta 1994 an Ghosal, Ghosh, Samanta 1995) lim t it = i 9

10 Commentary Improvement is a function of variability of within to between subjects (/) As long as ajustment is ampene enough, then result is asymptotically true Danger is when you o not ampen enough; see Deming s concern about ajustment in quality control setting Bayesian ajustment ampens enough Aaptive esigns have also been shown to be more efficient in this setting 10

11 ffect ffect Forms of Hill Moel Many forms available, all equivalent Hill moel was first evelope aroun 1903 = a + max h /(D 50 h + h ) = + ( range ½)/(1 + exp(h(log x log D 50 ))) Aliases Michaelis-Menten max 4-parameter Logistic Hill Moel Hill Moel (Log X) Dose Dose 11

12 12 Hill Moel Forms All are equivalent functions of the ose h h h K max 0 ) ( h h h K ) ( ) ( 0 0 ) ) / ( (1 ) ( max h K ) ) / ( (1 1 1 ) ( 0 h K P ]} [log exp{ 1 ) ( max 0 c h ) ) ( (1 1 1 ) ( 0 h P 2 1 ]} [log exp{ 1 1 ) ( max 0 c h

13 D-Optimal Design Let = ( 0, max, c, h) be the parameters Let y(; ) be the Hill moel Z = y/ Determinant is D = Z T Z Want to fin esign that maximizes D. 13

14 v 1 1 exp D 22 Optimal Define () log D 22 log max 2 Then () = exp()/[1 + exp()] 4 = 0. This can be solve approximately as opt 4(3 + 2e)/(3 + e) 2. This implies v 0.23 or

15 Robustness of Optimal Location D D22 Hill

16 Aaptive Design Start with optimal esign for first K subjects Upate istributions of parameters Use D-Bayes esign to pick next set of M subjects Repeat process Shoul beat fixe esign 16

17 Simulation Stuy Truth 0 = 40, max = -5, h = 1.5, = 1, c {1, 1.1,, 1.9, 2.0} This moel base on actual experience in clinical trial Sample size N = 300 Fixe esign Doses = 0, 20/3, 40/3, an 4 (equal allocation of patients) D-optimal esign 4 oses (equal allocation of patients) D-Bayes optimal esign Use uniform prior for K = 120, M = 20 Look at ability to accurate an precisely estimate log D 50, log D 70, an log D

18 Comparison of Performance Bias SD D 50 D 70 Re = fixe esign Blue = Bayes aaptive Green = D optimal D 90 18

19 ACR Treatment ffect Within Patient Aapting Back to RA motivating example We have several options Fixe esign A hoc D-optimal Bayesian aaptive esign Within patient aaptive npoint is ACR20 Binary If a patient achieve 20% reuction in symptoms, then set to 1; otherwise set to Dose 19

20 ACR20 proportion Personalize Design ACR20 is binary P(ACR20 = 1) follows Hill moel Hill Moels for Personalize Trial Question: Can we estimate D70 better than with Bayesian esign or fixe esign? Patients are seen at baseline, weeks 4, 8, 12, 16, an 24 If ACR20 = 0 for 3 consecutive visits, increase ose If ACR20 = 1 for 3 consecutive visits, ecrease ose Dose 20

21 Bias Stanar Deviation Simulation Results of Personalize Dosing Bias Stanar Deviation Fixe Bayes Pesonalize Fixe Bayes Pesonalize D50 Possibility D50 Possibility 21

22 More Challenging: Personalize Design on Binary npoint xten to PD enpoints (Reeve an Ferguson 2013) Objective is to estimate minimum effective ose (MD) MD efine as ose that yiel some p% response Consier binary enpoints (e.g., rheumatoi arthritis trial) Measure response repeately Assume P(therapeutic success) = log ose Upate preictors of 0 an 1 via Bayesian upates for each patient Upate patient s ose No TS for 3 observations, increase ose TS for 3 observations, ecrease ose Otherwise, keep at current ose 22

23 Personalize Design is More fficient Simulation results 23

24 Discussion of Themes Tie together the 3 types of esigns We have iscusse three seemingly isparate esigns, but they are linke by a common theme RCCT Bayesian aaptive Personalize All are using ata (Bayesian upating in all 3) to correct for mis-specification of our moel RCCT corrects for misspecifie exposure Bayesian aaptive an Personalize misspecifie ose-response» Coul be even better if we correct for misspecifie ose-response an exposure simultaneously By correcting for misspecifications, we achieve efficiency in esign 24

25 ACR Treatment ffect Conclusions an Wrapup RCCT an other Bayesian aaptations can significantly improve the efficiency of esigns ncourages us to think in terms of exposure-response moels These trials provie us with better estimates of the ose response Dose 25

Optimal Design for Hill Model

Optimal Design for Hill Model Optimal Design for Hill Moel DA 01 D-optimal; aaptive esign; personalize esign Russell Reeve, PhD Quintiles Durham, NC, USA Copyright 01 Quintiles Outline Motivating problem, from rheumatoi arthritis Hill

More information

The Role of Models in Model-Assisted and Model- Dependent Estimation for Domains and Small Areas

The Role of Models in Model-Assisted and Model- Dependent Estimation for Domains and Small Areas The Role of Moels in Moel-Assiste an Moel- Depenent Estimation for Domains an Small Areas Risto Lehtonen University of Helsini Mio Myrsylä University of Pennsylvania Carl-Eri Särnal University of Montreal

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Survey-weighted Unit-Level Small Area Estimation

Survey-weighted Unit-Level Small Area Estimation Survey-weighte Unit-Level Small Area Estimation Jan Pablo Burgar an Patricia Dörr Abstract For evience-base regional policy making, geographically ifferentiate estimates of socio-economic inicators are

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

Topic 7: Convergence of Random Variables

Topic 7: Convergence of Random Variables Topic 7: Convergence of Ranom Variables Course 003, 2016 Page 0 The Inference Problem So far, our starting point has been a given probability space (S, F, P). We now look at how to generate information

More information

A COMPARISON OF SMALL AREA AND CALIBRATION ESTIMATORS VIA SIMULATION

A COMPARISON OF SMALL AREA AND CALIBRATION ESTIMATORS VIA SIMULATION SAISICS IN RANSIION new series an SURVEY MEHODOLOGY 133 SAISICS IN RANSIION new series an SURVEY MEHODOLOGY Joint Issue: Small Area Estimation 014 Vol. 17, No. 1, pp. 133 154 A COMPARISON OF SMALL AREA

More information

Cascaded redundancy reduction

Cascaded redundancy reduction Network: Comput. Neural Syst. 9 (1998) 73 84. Printe in the UK PII: S0954-898X(98)88342-5 Cascae reunancy reuction Virginia R e Sa an Geoffrey E Hinton Department of Computer Science, University of Toronto,

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Test of Hypotheses in a Time Trend Panel Data Model with Serially Correlated Error Component Disturbances

Test of Hypotheses in a Time Trend Panel Data Model with Serially Correlated Error Component Disturbances HE UNIVERSIY OF EXAS A SAN ANONIO, COLLEGE OF BUSINESS Working Paper SERIES Date September 25, 205 WP # 000ECO-66-205 est of Hypotheses in a ime ren Panel Data Moel with Serially Correlate Error Component

More information

Tutorial on Maximum Likelyhood Estimation: Parametric Density Estimation

Tutorial on Maximum Likelyhood Estimation: Parametric Density Estimation Tutorial on Maximum Likelyhoo Estimation: Parametric Density Estimation Suhir B Kylasa 03/13/2014 1 Motivation Suppose one wishes to etermine just how biase an unfair coin is. Call the probability of tossing

More information

VI. Linking and Equating: Getting from A to B Unleashing the full power of Rasch models means identifying, perhaps conceiving an important aspect,

VI. Linking and Equating: Getting from A to B Unleashing the full power of Rasch models means identifying, perhaps conceiving an important aspect, VI. Linking an Equating: Getting from A to B Unleashing the full power of Rasch moels means ientifying, perhaps conceiving an important aspect, efining a useful construct, an calibrating a pool of relevant

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: EBLUP Unit Level for Small Area Estimation Contents General section... 3 1. Summary... 3 2. General

More information

CONTROL CHARTS FOR VARIABLES

CONTROL CHARTS FOR VARIABLES UNIT CONTOL CHATS FO VAIABLES Structure.1 Introuction Objectives. Control Chart Technique.3 Control Charts for Variables.4 Control Chart for Mean(-Chart).5 ange Chart (-Chart).6 Stanar Deviation Chart

More information

Image Denoising Using Spatial Adaptive Thresholding

Image Denoising Using Spatial Adaptive Thresholding International Journal of Engineering Technology, Management an Applie Sciences Image Denoising Using Spatial Aaptive Thresholing Raneesh Mishra M. Tech Stuent, Department of Electronics & Communication,

More information

State observers and recursive filters in classical feedback control theory

State observers and recursive filters in classical feedback control theory State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent

More information

Bayesian Updating: Discrete Priors: Spring

Bayesian Updating: Discrete Priors: Spring Bayesian Updating: Discrete Priors: 18.05 Spring 2017 http://xkcd.com/1236/ Learning from experience Which treatment would you choose? 1. Treatment 1: cured 100% of patients in a trial. 2. Treatment 2:

More information

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask 5.4 FUNDAMENTAL THEOREM OF CALCULUS Do you remember the Funamental Theorem of Algebra? Just thought I' ask The Funamental Theorem of Calculus has two parts. These two parts tie together the concept of

More information

PD Controller for Car-Following Models Based on Real Data

PD Controller for Car-Following Models Based on Real Data PD Controller for Car-Following Moels Base on Real Data Xiaopeng Fang, Hung A. Pham an Minoru Kobayashi Department of Mechanical Engineering Iowa State University, Ames, IA 5 Hona R&D The car following

More information

Neural Network Controller for Robotic Manipulator

Neural Network Controller for Robotic Manipulator MMAE54 Robotics- Class Project Paper Neural Network Controller for Robotic Manipulator Kai Qian Department of Biomeical Engineering, Illinois Institute of echnology, Chicago, IL 666 USA. Introuction Artificial

More information

Exam 2 Answers Math , Fall log x dx = x log x x + C. log u du = 1 3

Exam 2 Answers Math , Fall log x dx = x log x x + C. log u du = 1 3 Exam Answers Math -, Fall 7. Show, using any metho you like, that log x = x log x x + C. Answer: (x log x x+c) = x x + log x + = log x. Thus log x = x log x x+c.. Compute these. Remember to put boxes aroun

More information

A Course in Machine Learning

A Course in Machine Learning A Course in Machine Learning Hal Daumé III 12 EFFICIENT LEARNING So far, our focus has been on moels of learning an basic algorithms for those moels. We have not place much emphasis on how to learn quickly.

More information

Wiener Deconvolution: Theoretical Basis

Wiener Deconvolution: Theoretical Basis Wiener Deconvolution: Theoretical Basis The Wiener Deconvolution is a technique use to obtain the phase-velocity ispersion curve an the attenuation coefficients, by a two-stations metho, from two pre-processe

More information

Equilibrium in Queues Under Unknown Service Times and Service Value

Equilibrium in Queues Under Unknown Service Times and Service Value University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University

More information

Spurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009

Spurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009 Spurious Significance of reatment Effects in Overfitte Fixe Effect Moels Albrecht Ritschl LSE an CEPR March 2009 Introuction Evaluating subsample means across groups an time perios is common in panel stuies

More information

Maximal Causes for Non-linear Component Extraction

Maximal Causes for Non-linear Component Extraction Journal of Machine Learning Research 9 (2008) 1227-1267 Submitte 5/07; Revise 11/07; Publishe 6/08 Maximal Causes for Non-linear Component Extraction Jörg Lücke Maneesh Sahani Gatsby Computational Neuroscience

More information

by using the derivative rules. o Building blocks: d

by using the derivative rules. o Building blocks: d Calculus for Business an Social Sciences - Prof D Yuen Eam Review version /9/01 Check website for any poste typos an upates Eam is on Sections, 5, 6,, 1,, Derivatives Rules Know how to fin the formula

More information

Connections Between Duality in Control Theory and

Connections Between Duality in Control Theory and Connections Between Duality in Control heory an Convex Optimization V. Balakrishnan 1 an L. Vanenberghe 2 Abstract Several important problems in control theory can be reformulate as convex optimization

More information

PubH 7470: STATISTICS FOR TRANSLATIONAL & CLINICAL RESEARCH

PubH 7470: STATISTICS FOR TRANSLATIONAL & CLINICAL RESEARCH PubH 7470: STATISTICS FOR TRANSLATIONAL & CLINICAL RESEARCH The First Step: SAMPLE SIZE DETERMINATION THE ULTIMATE GOAL The most important, ultimate step of any of clinical research is to do draw inferences;

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Upsetting the contingency table: Causal induction over sequences of point events

Upsetting the contingency table: Causal induction over sequences of point events Upsetting the contingency table: Causal inuction over sequences of point events Michael D. Pacer (mpacer@berkeley.eu) Thomas L. Griffiths (tom griffiths@berkeley.eu) Department of Psychology, 321 Tolman

More information

A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL

A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL 6th European Signal Processing Conference EUSIPCO 28, Lausanne, Switzerlan, August 25-29, 28, copyright by EURASIP A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL Leonaro Tomazeli

More information

Bayesian Hierarchical Models for Data Extrapolation and Analysis in Rare and Pediatric Disease Clinical Trials

Bayesian Hierarchical Models for Data Extrapolation and Analysis in Rare and Pediatric Disease Clinical Trials Bayesian Hierarchical Moels for Data Extrapolation an Analysis in Rare an Peiatric Disease Clinical Trials A DISSERTATION SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA

More information

Optimal Signal Detection for False Track Discrimination

Optimal Signal Detection for False Track Discrimination Optimal Signal Detection for False Track Discrimination Thomas Hanselmann Darko Mušicki Dept. of Electrical an Electronic Eng. Dept. of Electrical an Electronic Eng. The University of Melbourne The University

More information

Personalized Treatment Selection Based on Randomized Clinical Trials. Tianxi Cai Department of Biostatistics Harvard School of Public Health

Personalized Treatment Selection Based on Randomized Clinical Trials. Tianxi Cai Department of Biostatistics Harvard School of Public Health Personalized Treatment Selection Based on Randomized Clinical Trials Tianxi Cai Department of Biostatistics Harvard School of Public Health Outline Motivation A systematic approach to separating subpopulations

More information

Collapsed Variational Inference for HDP

Collapsed Variational Inference for HDP Collapse Variational Inference for HDP Yee W. Teh Davi Newman an Max Welling Publishe on NIPS 2007 Discussion le by Iulian Pruteanu Outline Introuction Hierarchical Bayesian moel for LDA Collapse VB inference

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A

More information

Real-time arrival prediction models for light rail train systems EDOUARD NAYE

Real-time arrival prediction models for light rail train systems EDOUARD NAYE DEGREE PROJECT IN TRANSPORT AND LOCATION ANALYSIS STOCKHOLM, SWEDEN 14 Real-time arrival preiction moels for light rail train systems EDOUARD NAYE KTH ROYAL INSTITUTE OF TECHNOLOGY SCHOOL OF ARCHITECTURE

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

A comparison of small area estimators of counts aligned with direct higher level estimates

A comparison of small area estimators of counts aligned with direct higher level estimates A comparison of small area estimators of counts aligne with irect higher level estimates Giorgio E. Montanari, M. Giovanna Ranalli an Cecilia Vicarelli Abstract Inirect estimators for small areas use auxiliary

More information

How to Minimize Maximum Regret in Repeated Decision-Making

How to Minimize Maximum Regret in Repeated Decision-Making How to Minimize Maximum Regret in Repeate Decision-Making Karl H. Schlag July 3 2003 Economics Department, European University Institute, Via ella Piazzuola 43, 033 Florence, Italy, Tel: 0039-0-4689, email:

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

Constraint Reformulation and a Lagrangian Relaxation based Solution Algorithm for a Least Expected Time Path Problem Abstract 1.

Constraint Reformulation and a Lagrangian Relaxation based Solution Algorithm for a Least Expected Time Path Problem Abstract 1. Constraint Reformulation an a Lagrangian Relaation base Solution Algorithm for a Least Epecte Time Path Problem Liing Yang State Key Laboratory of Rail Traffic Control an Safety, Being Jiaotong University,

More information

Math 115 Section 018 Course Note

Math 115 Section 018 Course Note Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of

More information

Bayesian Updating: Discrete Priors: Spring

Bayesian Updating: Discrete Priors: Spring Bayesian Updating: Discrete Priors: 18.05 Spring 2017 http://xkcd.com/1236/ Learning from experience Which treatment would you choose? 1. Treatment 1: cured 100% of patients in a trial. 2. Treatment 2:

More information

β ˆ j, and the SD path uses the local gradient

β ˆ j, and the SD path uses the local gradient Proceeings of the 00 Winter Simulation Conference E. Yücesan, C.-H. Chen, J. L. Snowon, an J. M. Charnes, es. RESPONSE SURFACE METHODOLOGY REVISITED Ebru Angün Jack P.C. Kleijnen Department of Information

More information

CMSC 313 Preview Slides

CMSC 313 Preview Slides CMSC 33 Preview Slies These are raft slies. The actual slies presente in lecture may be ifferent ue to last minute changes, scheule slippage,... UMBC, CMSC33, Richar Chang CMSC 33 Lecture

More information

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Two-stage Adaptive Randomization for Delayed Response in Clinical Trials

Two-stage Adaptive Randomization for Delayed Response in Clinical Trials Two-stage Adaptive Randomization for Delayed Response in Clinical Trials Guosheng Yin Department of Statistics and Actuarial Science The University of Hong Kong Joint work with J. Xu PSI and RSS Journal

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

A Unified Theorem on SDP Rank Reduction

A Unified Theorem on SDP Rank Reduction A Unifie heorem on SDP Ran Reuction Anthony Man Cho So, Yinyu Ye, Jiawei Zhang November 9, 006 Abstract We consier the problem of fining a low ran approximate solution to a system of linear equations in

More information

Planar sheath and presheath

Planar sheath and presheath 5/11/1 Flui-Poisson System Planar sheath an presheath 1 Planar sheath an presheath A plasma between plane parallel walls evelops a positive potential which equalizes the rate of loss of electrons an ions.

More information

Introduction. A Dirichlet Form approach to MCMC Optimal Scaling. MCMC idea

Introduction. A Dirichlet Form approach to MCMC Optimal Scaling. MCMC idea Introuction A Dirichlet Form approach to MCMC Optimal Scaling Markov chain Monte Carlo (MCMC quotes: Metropolis et al. (1953, running coe on the Los Alamos MANIAC: a feasible approach to statistical mechanics

More information

Lecture 6 : Dimensionality Reduction

Lecture 6 : Dimensionality Reduction CPS290: Algorithmic Founations of Data Science February 3, 207 Lecture 6 : Dimensionality Reuction Lecturer: Kamesh Munagala Scribe: Kamesh Munagala In this lecture, we will consier the roblem of maing

More information

A Modification of the Jarque-Bera Test. for Normality

A Modification of the Jarque-Bera Test. for Normality Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam

More information

Leaving Randomness to Nature: d-dimensional Product Codes through the lens of Generalized-LDPC codes

Leaving Randomness to Nature: d-dimensional Product Codes through the lens of Generalized-LDPC codes Leaving Ranomness to Nature: -Dimensional Prouct Coes through the lens of Generalize-LDPC coes Tavor Baharav, Kannan Ramchanran Dept. of Electrical Engineering an Computer Sciences, U.C. Berkeley {tavorb,

More information

Transformations of Random Variables

Transformations of Random Variables Transformations of Ranom Variables September, 2009 We begin with a ranom variable an we want to start looking at the ranom variable Y = g() = g where the function g : R R. The inverse image of a set A,

More information

Influence of weight initialization on multilayer perceptron performance

Influence of weight initialization on multilayer perceptron performance Influence of weight initialization on multilayer perceptron performance M. Karouia (1,2) T. Denœux (1) R. Lengellé (1) (1) Université e Compiègne U.R.A. CNRS 817 Heuiasyc BP 649 - F-66 Compiègne ceex -

More information

Lecture 6: Calculus. In Song Kim. September 7, 2011

Lecture 6: Calculus. In Song Kim. September 7, 2011 Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear

More information

Event based Kalman filter observer for rotary high speed on/off valve

Event based Kalman filter observer for rotary high speed on/off valve 28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 28 WeC9.6 Event base Kalman filter observer for rotary high spee on/off valve Meng Wang, Perry Y. Li ERC for Compact

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Introduction to Machine Learning

Introduction to Machine Learning How o you estimate p(y x)? Outline Contents Introuction to Machine Learning Logistic Regression Varun Chanola April 9, 207 Generative vs. Discriminative Classifiers 2 Logistic Regression 2 3 Logistic Regression

More information

Subspace Estimation from Incomplete Observations: A High-Dimensional Analysis

Subspace Estimation from Incomplete Observations: A High-Dimensional Analysis Subspace Estimation from Incomplete Observations: A High-Dimensional Analysis Chuang Wang, Yonina C. Elar, Fellow, IEEE an Yue M. Lu, Senior Member, IEEE Abstract We present a high-imensional analysis

More information

A. Incorrect! The letter t does not appear in the expression of the given integral

A. Incorrect! The letter t does not appear in the expression of the given integral AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Sample Size and Power I: Binary Outcomes. James Ware, PhD Harvard School of Public Health Boston, MA

Sample Size and Power I: Binary Outcomes. James Ware, PhD Harvard School of Public Health Boston, MA Sample Size and Power I: Binary Outcomes James Ware, PhD Harvard School of Public Health Boston, MA Sample Size and Power Principles: Sample size calculations are an essential part of study design Consider

More information

Topic Modeling: Beyond Bag-of-Words

Topic Modeling: Beyond Bag-of-Words Hanna M. Wallach Cavenish Laboratory, University of Cambrige, Cambrige CB3 0HE, UK hmw26@cam.ac.u Abstract Some moels of textual corpora employ text generation methos involving n-gram statistics, while

More information

CONFIRMATORY FACTOR ANALYSIS

CONFIRMATORY FACTOR ANALYSIS 1 CONFIRMATORY FACTOR ANALYSIS The purpose of confirmatory factor analysis (CFA) is to explain the pattern of associations among a set of observe variables in terms of a smaller number of unerlying latent

More information

AP Calculus. Derivatives and Their Applications. Presenter Notes

AP Calculus. Derivatives and Their Applications. Presenter Notes AP Calculus Derivatives an Their Applications Presenter Notes 2017 2018 EDITION Copyright 2017 National Math + Science Initiative, Dallas, Texas. All rights reserve. Visit us online at www.nms.org Copyright

More information

Fall 2016: Calculus I Final

Fall 2016: Calculus I Final Answer the questions in the spaces provie on the question sheets. If you run out of room for an answer, continue on the back of the page. NO calculators or other electronic evices, books or notes are allowe

More information

Linear Regression with Limited Observation

Linear Regression with Limited Observation Ela Hazan Tomer Koren Technion Israel Institute of Technology, Technion City 32000, Haifa, Israel ehazan@ie.technion.ac.il tomerk@cs.technion.ac.il Abstract We consier the most common variants of linear

More information

The new concepts of measurement error s regularities and effect characteristics

The new concepts of measurement error s regularities and effect characteristics The new concepts of measurement error s regularities an effect characteristics Ye Xiaoming[1,] Liu Haibo [3,,] Ling Mo[3] Xiao Xuebin [5] [1] School of Geoesy an Geomatics, Wuhan University, Wuhan, Hubei,

More information

Simple deconvolution of time varying signal of gas phase variation in a 12 MW CFB biomass boiler

Simple deconvolution of time varying signal of gas phase variation in a 12 MW CFB biomass boiler Simple econvolution of time varying signal of gas phase variation in a MW CFB biomass boiler Byungho Song*, Anreas Johansson, Lars-Erik Åman, Filip Johnsson, an Bo Leckner *Department of Chemical Engineering,

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

A Minimum Variance Method for Lidar Signal Inversion

A Minimum Variance Method for Lidar Signal Inversion 468 J O U R N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y VOLUME 31 A Minimum Variance Metho for Liar Signal Inversion ANDREJA SU SNIK Centre for Atmospheric Research, University

More information

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow

More information

Designing of Acceptance Double Sampling Plan for Life Test Based on Percentiles of Exponentiated Rayleigh Distribution

Designing of Acceptance Double Sampling Plan for Life Test Based on Percentiles of Exponentiated Rayleigh Distribution International Journal of Statistics an Systems ISSN 973-675 Volume, Number 3 (7), pp. 475-484 Research Inia Publications http://www.ripublication.com Designing of Acceptance Double Sampling Plan for Life

More information

Outline. Calculus for the Life Sciences. Average Height and Weight of Girls 1. Chain Rule. Lecture Notes Chain Rule

Outline. Calculus for the Life Sciences. Average Height and Weight of Girls 1. Chain Rule. Lecture Notes Chain Rule Outline Calculus for the Life Sciences Lecture Notes Joseph M. Mahaffy, jmahaffy@mail.ssu.eu Department of Mathematics an Statistics Dynamical Systems Group Computational Sciences Research Center San Diego

More information

ELEC E7210: Communication Theory. Lecture 4: Equalization

ELEC E7210: Communication Theory. Lecture 4: Equalization ELEC E7210: Communication Theory Lecture 4: Equalization Equalization Delay sprea ISI irreucible error floor if the symbol time is on the same orer as the rms elay sprea. DF: Equalization a receiver signal

More information

Non-Linear Bayesian CBRN Source Term Estimation

Non-Linear Bayesian CBRN Source Term Estimation Non-Linear Bayesian CBRN Source Term Estimation Peter Robins Hazar Assessment, Simulation an Preiction Group Dstl Porton Down, UK. probins@stl.gov.uk Paul Thomas Hazar Assessment, Simulation an Preiction

More information

Adaptive Adjustment of Noise Covariance in Kalman Filter for Dynamic State Estimation

Adaptive Adjustment of Noise Covariance in Kalman Filter for Dynamic State Estimation Aaptive Ajustment of Noise Covariance in Kalman Filter for Dynamic State Estimation Shahroh Ahlaghi, Stuent Member, IEEE Ning Zhou, Senior Member, IEEE Electrical an Computer Engineering Department, Binghamton

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

Learning Automata in Games with Memory with Application to Circuit-Switched Routing

Learning Automata in Games with Memory with Application to Circuit-Switched Routing Learning Automata in Games with Memory with Application to Circuit-Switche Routing Murat Alanyali Abstract A general setting is consiere in which autonomous users interact by means of a finite-state controlle

More information

American Society of Agricultural Engineers PAPER NO PRAIRIE RAINFALL,CHARACTERISTICS

American Society of Agricultural Engineers PAPER NO PRAIRIE RAINFALL,CHARACTERISTICS - PAPER NO. 79-2108 PRAIRIE RAINFALL,CHARACTERISTICS G.E. Dyck an D.M. Gray Research Engineer an Chairman Division of Hyrology University of Saskatchewan Saskatoon, Saskatchewan, Canaa For presentation

More information

1. Aufgabenblatt zur Vorlesung Probability Theory

1. Aufgabenblatt zur Vorlesung Probability Theory 24.10.17 1. Aufgabenblatt zur Vorlesung By (Ω, A, P ) we always enote the unerlying probability space, unless state otherwise. 1. Let r > 0, an efine f(x) = 1 [0, [ (x) exp( r x), x R. a) Show that p f

More information

APPPHYS 217 Thursday 8 April 2010

APPPHYS 217 Thursday 8 April 2010 APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q

More information

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device Close an Open Loop Optimal Control of Buffer an Energy of a Wireless Device V. S. Borkar School of Technology an Computer Science TIFR, umbai, Inia. borkar@tifr.res.in A. A. Kherani B. J. Prabhu INRIA

More information

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 29 Sep 1997

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 29 Sep 1997 Dopant Spatial Distributions: Sample Inepenent Response Function An Maximum Entropy Reconstruction D. P. Chu an M. G. Dowsett Department of Physics, University of Warwick, Coventry CV4 7AL, UK arxiv:con-mat/970933v

More information

arxiv: v1 [hep-lat] 19 Nov 2013

arxiv: v1 [hep-lat] 19 Nov 2013 HU-EP-13/69 SFB/CPP-13-98 DESY 13-225 Applicability of Quasi-Monte Carlo for lattice systems arxiv:1311.4726v1 [hep-lat] 19 ov 2013, a,b Tobias Hartung, c Karl Jansen, b Hernan Leovey, Anreas Griewank

More information

Monte Carlo Methods with Reduced Error

Monte Carlo Methods with Reduced Error Monte Carlo Methos with Reuce Error As has been shown, the probable error in Monte Carlo algorithms when no information about the smoothness of the function is use is Dξ r N = c N. It is important for

More information

MA119-A Applied Calculus for Business Fall Homework 5 Solutions Due 10/4/ :30AM

MA119-A Applied Calculus for Business Fall Homework 5 Solutions Due 10/4/ :30AM MA9-A Applie Calculus for Business 2006 Fall Homework 5 Solutions Due 0/4/2006 0:0AM. #0 Fin the erivative of the function by using the rules of i erentiation. r f (r) r. #4 Fin the erivative of the function

More information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.2 Shot peening - modeling 3 PROCEEDINGS 3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:

More information

Deliverable 2.2. Small Area Estimation of Indicators on Poverty and Social Exclusion

Deliverable 2.2. Small Area Estimation of Indicators on Poverty and Social Exclusion Deliverable 2.2 Small Area Estimation of Inicators on Poverty an Social Exclusion Version: 2011 Risto Lehtonen, Ari Veijanen, Mio Myrsylä an Maria Valaste The project FP7-SSH-2007-217322 AMELI is supporte

More information

Learning in Monopolies with Delayed Price Information

Learning in Monopolies with Delayed Price Information Learning in Monopolies with Delaye Price Information Akio Matsumoto y Chuo University Ferenc Sziarovszky z University of Pécs February 28, 2013 Abstract We call the intercept of the price function with

More information

Gaussian processes with monotonicity information

Gaussian processes with monotonicity information Gaussian processes with monotonicity information Anonymous Author Anonymous Author Unknown Institution Unknown Institution Abstract A metho for using monotonicity information in multivariate Gaussian process

More information

KNN Particle Filters for Dynamic Hybrid Bayesian Networks

KNN Particle Filters for Dynamic Hybrid Bayesian Networks KNN Particle Filters for Dynamic Hybri Bayesian Networs H. D. Chen an K. C. Chang Dept. of Systems Engineering an Operations Research George Mason University MS 4A6, 4400 University Dr. Fairfax, VA 22030

More information

Robustifying Trial-Derived Treatment Rules to a Target Population

Robustifying Trial-Derived Treatment Rules to a Target Population 1/ 39 Robustifying Trial-Derived Treatment Rules to a Target Population Yingqi Zhao Public Health Sciences Division Fred Hutchinson Cancer Research Center Workshop on Perspectives and Analysis for Personalized

More information

TCP throughput and timeout steady state and time-varying dynamics

TCP throughput and timeout steady state and time-varying dynamics TCP throughput an timeout steay state an time-varying ynamics Stephan Bohacek an Khushboo Shah Dept. of Electrical an Computer Engineering Dept. of Electrical Engineering University of Delaware University

More information