Plot of Laser Operating Current as a Function of Time Laser Test Data
|
|
- Linda Chandler
- 5 years ago
- Views:
Transcription
1 Chapter Repeated Measures Data and Random Parameter Models Repeated Measures Data and Random Parameter Models Chapter Objectives Understand applications of growth curve models to describe the results of statistical studies in which repeated measures are made on a sample of units from some population or process. Part of the Iowa State University NSF/ILI project Beyond Traditional Statistical Methods Copyright 999 D. Cook, W. M. Duckworth, M. S. Kaiser, W. Q. Meeker, and W. R. Stephenson. Developed as part of NSF/ILI grant DUE97. January 3, 7h min - Explore different kinds of models for repeated measure data. Empirical Mechanistic (e.g., from systems of differential equations) Understand how to fit nonlinear regression models. Understand the different sources of variability in repeated measures data. Become familiar with methods for data analysis and inference for repeated measure data. - Plot of Laser Operating Current as a Function of Time Laser Test Data Percent Increase in Operating Current 3 Hours Measured percentage increase in operating current over time for GaAs lasers tested at C. Fifteen () devices each measured every hours up to hours of operation. For these device and the corresponding application, a D f = % increase in current was the specified failure level. Engineers wanted to predict life of lasers operating at C Plasma Concentrations of Indomethicin Following Bolus Intravenous Injection Constant-Rate Growth Model For some types of growth, rate will be approximately constant: dy(t) = β With y() = β, this simple differential equation has the solution y(t) =β + β t where β is the intercept and β is the slope (growth rate). This linear model is useful describe certain growth patterns. Concentration (mg/l)..... Hours - -
2 Plasma Concentrations of Indomethicin Following Bolus Intravenous Injection Plasma Concentrations of Indomethicin Following Intravenous Injection Blood plasma concentrations of indomethicin following bolus intravenous injection. Data from Kwan et al. (97) Kinetics of indomethicin absorption, elimination, and enterohepatic circulation in man. Also, page of Davidian and Giltinan (99) Data collected on six human volunteers at time points from minutes to to hours. Pharmacokinetic researchers need to know the rate of passing to elimination and the time (or expected time) until a specified level of indomethicin is left. log(concentration) (log(mcg/ml)) Interest centers on the identification of demographic and physiological characteristics explaining the variability in drug response information to be used in choosing effective dosage regimes Time (hr) -7 - Compartmental Model for Indomethicin Absorption and Elimination Two-Compartment Model for Blood Plasma Concentrations of Indomethicin Following Bolus Intravenous Injection A two-compartment model provides useful model to describe the concentration in the blood over time. The following differential equations describe the two compartment model da (t) = k A (t), and da (t) = k A (t) k A (t). Injection Site Blood k A A k Elimination Then, using standard methods for solving systems of linear differential equations the concentration in the blood at time t is A (t) = θ exp( k t) θ exp( k t) θ = A ()k /(k k ) θ = θ + A () Splus has the biexponential function for this model Biexponential Function plot.function(biexp,:,a=.9, A=3., lrc=-.7, lrc=.) Serum Concentrations of Theophylline Following Oral Administration response Concentration (mg/l) time Time since drug administration (Hours) - -
3 Trellis Plot of Serum Concentrations of Theophylline Following Oral Administration Serum Concentrations of Theophylline Following Oral Administration 9 Data from Boeckmann et al. (99). Theophylline is an anti-asthmatic agent. Data collected on twelve subjects. Serum concentrations were measured at time points from minutes to to hours. Theophylline concentration (mg/l) 7 3 Dosage was varied in proportion to body weight. Time since drug administration (hr) Two-Compartment Model for Theophylline Absorption and Elimination Following Oral Administration Two-Compartment Model for Theophylline Absorption and Elimination Following Oral Administration Stomach k Blood A A k =V/Cl Elimination - The following two-compartment (gut and blood) model is used to describe the concentration in the blood as a function of time. da (t) = k A (t), and where da (t) = k A (t) Cl V A (t) A (t) is the amount in the stomach at time t. A (t) is the concentration in the blood at time t. The rate k has units /hour. Clearance Cl has units L/hr/kg. Volume of distribution V has units L/kg. The elimination rate k = Cl/V has units /hour. - Two-Compartment Model for Theophylline Absorption and Elimination With initial conditions A () = Dose and A () =, the concentration in the blood as a function of time is obtained by solving the system of differential equations to give: A (t) = Dose k V (k Cl/V ) [ exp ( Cl V t The amount left in the stomach at time t is A (t) =Dose exp( k t) ) ] exp( k t) The model fit to the blood concentration data was reparameterized as A (t)= Dose exp(θ ) exp(θ ) exp(θ 3 )[exp(θ ) exp(θ )] {exp[ exp(θ )t] exp[ exp(θ )t]} where θ =log(k ), θ =log(cl/v ), θ 3 = log(cl), are unrestricted in sign. response First Order Log Function plot.function(first.order.log, time=:, Dose=,lCl=-3.7,lka= -., lke=-.) Splus has function first.order.log function for this model. -7 time -
4 Device-B Power Drop Accelerated Degradation Test Results at C, 9 C,and37 C (Use conditions C) Device-B Power Drop Accelerated Test Results Power drop in db Degrees C 9 Degrees C 37 Degrees C RF amplifier IC device, used in a satellite communications system. System life was designed to be between and years. A six-month test at high temperature was to be used to predict life at C junction temperature. The physics/chemistry of the failure mechanism was wellunderstood. 3 Hours -9 - Device-B Power Drop Simple One-Step Chemical Reaction Leading to Failure Simple One-Step Chemical Reaction A (t) is the amount of harmful material available for reaction at time t A (t) is proportional to the amount of failure-causing compounds at time t. k A A Simple one-step chemical reaction: A k A Power drop proportional to A (t) The rate equations for this reaction are da = k A and da = k A - - Device-B Power Drop Simple One-Step Chemical Reaction Leading to Failure (contd.) Solution to system of differential equations: A (t) = A () exp( k t) A (t) = A () + A ()[ exp( k t)] where A () and A () are initial conditions. If A () =, then lim t A (t) =A () and the solution for A (t) (the function of primary interest) can be reexpressed as A (t) = A ()[ exp( k t)] A simple -step diffusion process has the same solution. Acceleration of Degradation The Arrhenius model describing the effect that temperature has on the rate of a simple one-step chemical reaction is ( ) E a R(temp) =γ exp k B (temp + 73.) where temp is temperature in C and k B =. is Boltzmann s constant in units of electron volts per C. The pre-exponential factor γ and the reaction activation energy E a are characteristics of the product or material. The Acceleration Factor between temp and temp U is AF(temp) =AF(temp, temp U,E R(temp) a)= R(temp U ) When temp > temp U, AF(temp, temp U,E a) >. -3 -
5 Simple Reversible Chemical Reaction Simple Reversible Chemical Reaction The rate equations for this reaction are da = k A + k A and da = k A k A k A A k In this model the failure causing reaction is reversible. With A () =, the solution of the system of differential equations for this model gives, for A k A (t) =A () { exp[ (k + k )t]}. k + k - - Orange Tree Circumference Growth Orange Tree Circumference Growth Data collected by measuring trunk circumference on five different trees grown in Riverside, Calfornia, between 99 and 973. Measurements taken at,,,, 3, 37, and days. Researchers are interested in characterizing the tree-to-tree variability and seeking explanations for why some trees grow faster than others. Trunk circumference (mm) Data from Draper and Smith (9), page. Time since December 3, 9 (days) -7 - Trellis Plot of Orange Tree Circumference Growth Logistic (Autocatalytic) Model For some types of growth, growth rate will increase to a point and then decline. Trunk circumference (mm) 3 The following differential equation has rate depending on has a basic rate, κ, current size f, and a limiting size α. dy(t) = Asym scale y(t) [Asym y(t)]. This differential equation has the general solution y(t) =Asym/{ + exp[ (t T )/scale]} where Asym is the asymptote, T is the time at which % growth has been achieved, and scale is a scale parameter describing the steepness of the growth curve. This S-shaped logistic curve is a useful curve to describe certain growth patterns. Time since December 3, 9 (days). -9 Splus has the logistic function for this model. -3
6 Orange Tree Circumference Growth Data and Fitted Logistic Function Logistic Function plot.function(logistic,:,asym=7,t=,scal=379) Orange Tree Circumference Growth response inches time Days -3-3 Fatigue Crack Size Observations for Alloy-A (Bogdanoff & Kozin 9) Bogdonoff-Kozin Fatigue Crack Growth Data Notched Compact Fatigue Test Specimen Inches Mcycles Alloy-A Fatigue Crack-Size Data Fatigue Crack Size Observations for Alloy-A (Bogdanoff & Kozin 9) Data from Hudak, Saxena, Bucci, and Malcolm (97) Bogdanoff and Kozin (9, page ), and Lu and Meeker (993). Bogdonoff-Kozin Fatigue Crack Growth Data Suppose investigators wanted to:... Estimate materials parameters related to crack growth. Inches Estimate time (measured in number of cycles) at which % of the cracks would reach. inches Assess adequacy of the paris model to describe crack growth rate Mcycles -3-3
7 Fatigue Crack Size Observations for Alloy-A and Nonlinear Least Squares Paris Law Fit Bogdonoff-Kozin Fatigue Crack Growth Data Degradation Data in Reliability Analysis Sometimes possible to measure degradation directly over time Inches.... Continuously. At specific points in time. Degradation is natural response for some tests. Degradation data can provide considerably more reliability information than censored failure-time data (especially with few or no failures) Mcycles -37 Direct observation of the degradation process allows direct modeling of the failure-causing mechanism. -3 Paris Crack Growth Model Paris Model with no Variability The paris model is da(t) = C [ K(a)] m is a commonly used empirical model to describe the growth of fatigue cracks over some range of size. C>andm> are materials properties K(a) is the stress intensity function. Form of K(a) depends on applied stress, part dimensions, and geometry. To model a two-dimensional edge-crack in a plate with a crack that is small relative to the wih of the plate (say less than 3%), K(a) =Stress πa and the solution to the resulting differential equation is [ {a()} m/ +( m/) C (Stress π) m t ] m, m a(t) = a() exp [ C (Stress π) t ], m = Crack Size (mm)..... Cycles Paris Model with Unit-to-Unit Variability in Initial Crack Size but with Fixed Materials Parameters and Constant Stress. Models for Variation in Degradation and Failure Time If all manufactured units were identical, operated at exactly the same time, under exactly the same conditions, and in exactly the same environment, and if every unit failed as it reached a particular critical level of degradation, then all units would fail at exactly the same time. Crack Size (mm)... Need to identify and model important sources of variability in the degradation process. Quantities that might be modeled as random include: Initial conditions (flaw size, amount of material). Materials parameters (related to degradation rate). Level of degradation at which unit will fail.. Cycles Stochastic process variability (e.g., stress of other environmental variables changing over time). - -
8 Paris Model with Unit-to-Unit Variability in the Initial Crack Size and Materials Parameters but Constant Stress Paris Model with Unit-to-Unit Variability in the Initial Crack Size and Materials Parameters and Stochastic Stress.... Crack Size (mm). Crack Size (mm)..... Cycles -3 Cycles - General Response Path Model Estimation of Degradation Model Parameters D ij = D(t ij,β i,...,β ki ) is the degradation path for unit i at time t (measured in hours, cycles, etc.). Observed sample path of unit i at time t j is y ij = D ij + ɛ ij, i =,...,n, j =,...,m i Residuals ɛ ij NOR (,σ ɛ) describe a combination of measurement error and model error. For unit i, β i,...,β ki is a vector of k unknown parameters. Some of the β i,...,β ki are random from unit to unit. Model appropriate function of β i,...,β ki with multivariate normal distribution (MVN) with parameters µ β and Σ β. The likelihood for the random-parameter path model is L(µ β, Σ β,σ ɛ DATA) n mi =... φ nor(ζ ij) f β (β i,...,β ki ; µ β, Σ β )dβ i,...,dβ ki σ ɛ i= j= where ζ ij =[y ij D(t ij,β i,...,β ki )]/σ ɛ and f β (β i,...,β ki ; µ β, Σ β ) is the multivariate normal distribution density function. Each evaluation of L(µ β, Σ β,σ ɛ DATA) will, in general, require numerical approximation of n integrals of dimension k. Maximization of L(µ β, Σ β,σ ɛ DATA) computationally difficult. - - Alloy-A Fatigue Crack Size Observations and Fitted Paris-Rule Model Crack Length.... Estimates of Fatigue Data Model Parameters for Alloy-A Splus gives the following approximate ML estimates. ( ) ( ).7..9 ˆµ β =, ˆΣ 3.73 β =.9.9 and σ ɛ =.3. Here β = C, β = m Millions of Cycles -7 -
9 Estimates β i Versus β i, i =,..., and Contours for the Fitted Bivariate Normal Distribution 7 Models Relating Degradation and Failure Beta Soft failures are defined to occur at a specified degradation level. In some products there is a gradual loss of performance (e.g., decreasing light output from a fluorescent light bulb). Use fixed D f to denote the critical level defining failure for a degradation path. 3 Beta -9 - Alloy-A Fatigue Crack Size Observations and Fitted Paris-Rule Model Evaluation of F (t) Distribution of Time to First Crossing Crack Length Millions of Cycles Direct evaluation of F (t): Closed forms available for simple problems (e.g., a single random variable and other special cases). Numerical integration: Useful for a small number of random variables (e.g., or 3). FORM (first order) approximation: Rapid computation, but uncertain approximation. Used frequently in engineering problems with large number of random variables. Monte Carlo simulation: General method. Needs much computer time to evaluate small probabilities. Can use importance sampling. Estimate failure probabilities by evaluating at ML estimates. - - Evaluation of F (t) by Numerical Integration Bootstrap Estimates of F (t) =Pr(T t) The failure (crossing probability) can be expressed as Pr(T t) =F (t) =F (t; θ β )=Pr[D(t, β,...,β k ) > D f ]. If (β,β ) follows a bivariate normal distribution with parameters µ β,µ β σ β,σ β,ρ β,β,thenp (T t) = Φnor g(d f,t,β ) µ β β σ β β φ nor σ β ( ) β µ β dβ σ β where g(d f,t,β )isthevalueofβ for given β, that gives D(t) =D f. Cumulative Probability.... Method generalizes to multivariate normal, but requires correspondingly higher-order integration Millions of Cycles -3 -
10 Degradation Estimate of F (t) with Pointwise Two-Sided 9% and % Bootstrap Bias-Corrected Percentile Confidence Intervals, Based on the Crack-Size Data Censored at t c =.. The Nonparametric Estimate of F (t) Indicated by Dots Probability <- Censoring time Millions of Cycles cdf Estimate % Confidence Intervals 9% Confidence Intervals - > Indometh.lis <- nlslist(data=indometh, + conc ~ biexp(time, A, A,lrc,lrc), cluster= ~ Subject) A A lrc lrc A A lrc lrc A A lrc lrc A A lrc lrc A A lrc lrc A A lrc lrc > Indometh.nlme <- nlme(indometh.lis) > plot(indometh.nlme) - Fitted Values of Plasma Concentrations of Indomethicin Following Intravenous Injection Plasma Concentrations of Indomethicin Estimates of Random Parameters A -... Concentration (mcg/ml).... Time since drug administration (Hours) A lrc lrc > anova(indometh.nlme) Value Std.Error z ratio A A lrc lrc > > summary(indometh.nlme) Call: Model: conc ~ biexp(time, A, A, lrc, lrc) Fixed: list(a ~., A ~., lrc ~., lrc ~.) Random: list(a ~., A ~., lrc ~., lrc ~.) Cluster: ~ Subject Data: Indometh Estimation Method: ML Convergence at iteration: Approximate Loglikelihood:.3 AIC: -9.7 BIC: Variance/Covariance Components Estimate(s): Structure: matrixlog Standard Deviation(s) of Random Effect(s) A A lrc lrc Correlation of Random Effects A A lrc A.993 lrc.9.33 lrc Cluster Residual Variance:.9 Fixed Effects Estimate(s): Value Approx. Std.Error z ratio(c) A A lrc lrc
11 Conditional Correlation(s) of Fixed Effects Estimates A A lrc A.39 lrc lrc Random Effects (Conditional Modes): A A lrc lrc Standardized Population-Average Residuals: Min Q Med Q3 Max Number of Observations: Number of Clusters: - > predict(indometh.nlme,data=indometh,cluster= ~ Subject) cluster fit.cluster fit.population Fitted Values of Serum Concentrations of Theophylline Following Oral Administration Fitted Values of Orange Tree Circumference Growth Concentration (mg/l) Time since drug administration (Hours) Trunk circumference (mm) Time since December 3, 9 (days) -3 - Trellis Plot of the Results of an Experiment to Compare Two Genotypes of Soybeans Summary Repeated measure data occur in many different areas of application from biological to engineering F 99 P For many purposes, it is better to have a mechanistic (e.g., one derived from kinetics) but empirical models can also be useful Leaf weight/plant (g) 99 F 99 P 3 Modeling the unit-to-unit variability in model parameters is useful for answering many important questions about variability in the population or process being studied. 3 9 F 9 P Modern software (e.g., the Splus function nlme) make it easy to fit nonlinear regression models with random parameters to reflect unit-to-unit variability. Days after planting - It is also possible to fit models with autocorrelated residual terms. -
nonlinear mixed effect model fitting with nlme
nonlinear mixed effect model fitting with nlme David Lamparter March 29, 2010 David Lamparter nonlinear mixed effect model fitting with nlme Purpose of nonlinear mixed effects modeling nonlinearity fitting
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 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 informationAccelerated Testing Obtaining Reliability Information Quickly
Accelerated Testing Background Accelerated Testing Obtaining Reliability Information Quickly William Q. Meeker Department of Statistics and Center for Nondestructive Evaluation Iowa State University Ames,
More informationChapter 18. Accelerated Test Models. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 18 Accelerated Test Models William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors text
More informationAccelerated Degradation Tests: Modeling and Analysis
Statistics Preprints Statistics 9-20-1999 Accelerated Degradation Tests: Modeling and Analysis William Q. Meeker Iowa State University, wqmeeker@iastate.edu Luis A. Escobar Louisiana State University C.
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 informationChapter 17. Failure-Time Regression Analysis. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 17 Failure-Time Regression Analysis William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors
More informationDegradation data analysis for samples under unequal operating conditions: a case study on train wheels
Degradation data analysis for samples under unequal operating conditions: a case study on train wheels Marta AFONSO FREITAS UNIVERSIDADE FEDERAL DE MINAS GERAIS DEPARTMENT OF PRODUCTION ENGINEERING marta.afreitas@gmail.com
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 informationn =10,220 observations. Smaller samples analyzed here to illustrate sample size effect.
Chapter 7 Parametric Likelihood Fitting Concepts: Chapter 7 Parametric Likelihood Fitting Concepts: Objectives Show how to compute a likelihood for a parametric model using discrete data. Show how to compute
More informationData: Singly right censored observations from a temperatureaccelerated
Chapter 19 Analyzing Accelerated Life Test Data William Q Meeker and Luis A Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W Q Meeker and L A Escobar Based on the authors
More informationArrhenius Plot. Sample StatFolio: arrhenius.sgp
Summary The procedure is designed to plot data from an accelerated life test in which failure times have been recorded and percentiles estimated at a number of different temperatures. The percentiles P
More informationSemiconductor Reliability
Semiconductor Reliability. Semiconductor Device Failure Region Below figure shows the time-dependent change in the semiconductor device failure rate. Discussions on failure rate change in time often classify
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 informationACCELERATED DESTRUCTIVE DEGRADATION TEST PLANNING. Presented by Luis A. Escobar Experimental Statistics LSU, Baton Rouge LA 70803
ACCELERATED DESTRUCTIVE DEGRADATION TEST PLANNING Presented by Luis A. Escobar Experimental Statistics LSU, Baton Rouge LA 70803 This is jointly work with Ying Shi and William Q. Meeker both from Iowa
More informationChapter 15. System Reliability Concepts and Methods. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 15 System Reliability Concepts and Methods William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on
More informationAccelerated Destructive Degradation Tests: Data, Models, and Analysis
Statistics Preprints Statistics 03 Accelerated Destructive Degradation Tests: Data, Models, and Analysis Luis A. Escobar Louisiana State University William Q. Meeker Iowa State University, wqmeeker@iastate.edu
More informationON THE CONSEQUENCES OF MISSPECIFING ASSUMPTIONS CONCERNING RESIDUALS DISTRIBUTION IN A REPEATED MEASURES AND NONLINEAR MIXED MODELLING CONTEXT
ON THE CONSEQUENCES OF MISSPECIFING ASSUMPTIONS CONCERNING RESIDUALS DISTRIBUTION IN A REPEATED MEASURES AND NONLINEAR MIXED MODELLING CONTEXT Rachid el Halimi and Jordi Ocaña Departament d Estadística
More informationBivariate Degradation Modeling Based on Gamma Process
Bivariate Degradation Modeling Based on Gamma Process Jinglun Zhou Zhengqiang Pan Member IAENG and Quan Sun Abstract Many highly reliable products have two or more performance characteristics (PCs). The
More informationMultistate Modeling and Applications
Multistate Modeling and Applications Yang Yang Department of Statistics University of Michigan, Ann Arbor IBM Research Graduate Student Workshop: Statistics for a Smarter Planet Yang Yang (UM, Ann Arbor)
More informationLifetime prediction and confidence bounds in accelerated degradation testing for lognormal response distributions with an Arrhenius rate relationship
Scholars' Mine Doctoral Dissertations Student Research & Creative Works Spring 01 Lifetime prediction and confidence bounds in accelerated degradation testing for lognormal response distributions with
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 informationChapter 9. Bootstrap Confidence Intervals. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 9 Bootstrap Confidence Intervals William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors
More informationComputationally Efficient Estimation of Multilevel High-Dimensional Latent Variable Models
Computationally Efficient Estimation of Multilevel High-Dimensional Latent Variable Models Tihomir Asparouhov 1, Bengt Muthen 2 Muthen & Muthen 1 UCLA 2 Abstract Multilevel analysis often leads to modeling
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 informationSTATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL. A Thesis. Presented to the. Faculty of. San Diego State University
STATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL A Thesis Presented to the Faculty of San Diego State University In Partial Fulfillment of the Requirements for the Degree
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 informationSTAT 6350 Analysis of Lifetime Data. Probability Plotting
STAT 6350 Analysis of Lifetime Data Probability Plotting Purpose of Probability Plots Probability plots are an important tool for analyzing data and have been particular popular in the analysis of life
More informationPrediction for Two-Phase Degradation. A thesis presented to. the faculty of. In partial fulfillment. of the requirements for the degree
Bayesian Degradation Analysis Considering Competing Risks and Residual-Life Prediction for Two-Phase Degradation A thesis presented to the faculty of the Russ College of Engineering and Technology of Ohio
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 informationParameters Estimation for a Linear Exponential Distribution Based on Grouped Data
International Mathematical Forum, 3, 2008, no. 33, 1643-1654 Parameters Estimation for a Linear Exponential Distribution Based on Grouped Data A. Al-khedhairi Department of Statistics and O.R. Faculty
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationLinear Mixed Models. One-way layout REML. Likelihood. Another perspective. Relationship to classical ideas. Drawbacks.
Linear Mixed Models One-way layout Y = Xβ + Zb + ɛ where X and Z are specified design matrices, β is a vector of fixed effect coefficients, b and ɛ are random, mean zero, Gaussian if needed. Usually think
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 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 informationOptimum Life Test Plans of Electrical Insulation for Thermal Stress. Hideo Hirose, Takenori Sakumura, Naoki Tabuchi and Takeru Kiyosue
Proceedings of the International MultiConference of Engineers and Computer Scientists 204 Vol II, IMECS 204, March 2-4, 204, Hong Kong Optimum Life Test Plans of Electrical Insulation for Thermal Stress
More informationBayesian Modeling of Accelerated Life Tests with Random Effects
Bayesian Modeling of Accelerated Life Tests with Random Effects Ramón V. León Avery J. Ashby Jayanth Thyagarajan Joint Statistical Meeting August, 00 Toronto, Canada Abstract We show how to use Bayesian
More informationPHA 4123 First Exam Fall On my honor, I have neither given nor received unauthorized aid in doing this assignment.
PHA 4123 First Exam Fall 1998 On my honor, I have neither given nor received unauthorized aid in doing this assignment. TYPED KEY Name Question 1. / pts 2. /2 pts 3. / pts 4. / pts. / pts 6. /10 pts 7.
More informationStandard Errors & Confidence Intervals. N(0, I( β) 1 ), I( β) = [ 2 l(β, φ; y) β i β β= β j
Standard Errors & Confidence Intervals β β asy N(0, I( β) 1 ), where I( β) = [ 2 l(β, φ; y) ] β i β β= β j We can obtain asymptotic 100(1 α)% confidence intervals for β j using: β j ± Z 1 α/2 se( β j )
More informationSimultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions
Statistics Preprints Statistics 10-2014 Simultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions Yimeng Xie Virginia Tech Yili Hong Virginia Tech Luis A. Escobar Louisiana
More informationBootstrapping the Confidence Intervals of R 2 MAD for Samples from Contaminated Standard Logistic Distribution
Pertanika J. Sci. & Technol. 18 (1): 209 221 (2010) ISSN: 0128-7680 Universiti Putra Malaysia Press Bootstrapping the Confidence Intervals of R 2 MAD for Samples from Contaminated Standard Logistic Distribution
More informationBayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features. Yangxin Huang
Bayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features Yangxin Huang Department of Epidemiology and Biostatistics, COPH, USF, Tampa, FL yhuang@health.usf.edu January
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 informationLife Prediction of Structural Components
Life Prediction of Structural Components Steve Duffy Cleveland State University, Connecticut Reserve Technologies, LLC Lesley A. Janosik NASA Glenn Research Center Andy Wereszczak Oak Ridge National Laboratory
More informationBeka 2 Cpt: Two Compartment Model - Loading Dose And Maintenance Infusion
MEDSCI 719 Pharmacometrics Beka 2 Cpt: Two Compartment Model - Loading Dose And Maintenance Infusion Objective 1. To observe the time course of drug concentration in the central and peripheral compartments
More informationLatent Variable Models for Binary Data. Suppose that for a given vector of explanatory variables x, the latent
Latent Variable Models for Binary Data Suppose that for a given vector of explanatory variables x, the latent variable, U, has a continuous cumulative distribution function F (u; x) and that the binary
More informationECON 5350 Class Notes Functional Form and Structural Change
ECON 5350 Class Notes Functional Form and Structural Change 1 Introduction Although OLS is considered a linear estimator, it does not mean that the relationship between Y and X needs to be linear. In this
More informationAccelerated Destructive Degradation Test Planning
Accelerated Destructive Degradation Test Planning Ying Shi Dept. of Statistics Iowa State University Ames, IA 50011 yshi@iastate.edu Luis A. Escobar Dept. of Experimental Statistics Louisiana State University
More informationMODIFIED MONTE CARLO WITH LATIN HYPERCUBE METHOD
MODIFIED MONTE CARLO WITH LATIN HYPERCUBE METHOD Latin hypercube sampling (LHS) was introduced by McKay, Conover and Beckman as a solution to increase the efficiency of computer simulations. This technique
More informationNon-maximum likelihood estimation and statistical inference for linear and nonlinear mixed models
Optimum Design for Mixed Effects Non-Linear and generalized Linear Models Cambridge, August 9-12, 2011 Non-maximum likelihood estimation and statistical inference for linear and nonlinear mixed models
More informationNotes on the Multivariate Normal and Related Topics
Version: July 10, 2013 Notes on the Multivariate Normal and Related Topics Let me refresh your memory about the distinctions between population and sample; parameters and statistics; population distributions
More informationQuantile POD for Hit-Miss Data
Quantile POD for Hit-Miss Data Yew-Meng Koh a and William Q. Meeker a a Center for Nondestructive Evaluation, Department of Statistics, Iowa State niversity, Ames, Iowa 50010 Abstract. Probability of detection
More informationGeneralized, Linear, and Mixed Models
Generalized, Linear, and Mixed Models CHARLES E. McCULLOCH SHAYLER.SEARLE Departments of Statistical Science and Biometrics Cornell University A WILEY-INTERSCIENCE PUBLICATION JOHN WILEY & SONS, INC. New
More informationUsing Non-Linear Mixed Models for Agricultural Data
Using Non-Linear Mixed Models for Agricultural Data Fernando E. Miguez Energy Biosciences Institute Crop Sciences University of Illinois, Urbana-Champaign miguez@illinois.edu Oct 8 th, 8 Outline 1 Introduction
More informationModel-Assisted Probability of Detection for Ultrasonic Structural Health Monitoring
4th European-American Workshop on Reliability of NDE - Th.2.A.2 Model-Assisted Probability of Detection for Ultrasonic Structural Health Monitoring Adam C. COBB and Jay FISHER, Southwest Research Institute,
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 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 informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationImproving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates
Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates Anastasios (Butch) Tsiatis Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
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 informationConstrained estimation for binary and survival data
Constrained estimation for binary and survival data Jeremy M. G. Taylor Yong Seok Park John D. Kalbfleisch Biostatistics, University of Michigan May, 2010 () Constrained estimation May, 2010 1 / 43 Outline
More informationSTAT Section 2.1: Basic Inference. Basic Definitions
STAT 518 --- Section 2.1: Basic Inference Basic Definitions Population: The collection of all the individuals of interest. This collection may be or even. Sample: A collection of elements of the population.
More informationBayesian inference for sample surveys. Roderick Little Module 2: Bayesian models for simple random samples
Bayesian inference for sample surveys Roderick Little Module : Bayesian models for simple random samples Superpopulation Modeling: Estimating parameters Various principles: least squares, method of moments,
More informationEstimation of AUC from 0 to Infinity in Serial Sacrifice Designs
Estimation of AUC from 0 to Infinity in Serial Sacrifice Designs Martin J. Wolfsegger Department of Biostatistics, Baxter AG, Vienna, Austria Thomas Jaki Department of Statistics, University of South Carolina,
More informationFinancial Econometrics
Financial Econometrics Estimation and Inference Gerald P. Dwyer Trinity College, Dublin January 2013 Who am I? Visiting Professor and BB&T Scholar at Clemson University Federal Reserve Bank of Atlanta
More informationSTAT 510 Final Exam Spring 2015
STAT 510 Final Exam Spring 2015 Instructions: The is a closed-notes, closed-book exam No calculator or electronic device of any kind may be used Use nothing but a pen or pencil Please write your name and
More informationReliability Engineering I
Happiness is taking the reliability final exam. Reliability Engineering I ENM/MSC 565 Review for the Final Exam Vital Statistics What R&M concepts covered in the course When Monday April 29 from 4:30 6:00
More informationBayesian Methods for Accelerated Destructive Degradation Test Planning
Statistics Preprints Statistics 11-2010 Bayesian Methods for Accelerated Destructive Degradation Test Planning Ying Shi Iowa State University William Q. Meeker Iowa State University, wqmeeker@iastate.edu
More informationRecent Advances in Reliability Estimation of Time-Dependent Problems Using the Concept of Composite Limit State
Automotive Research Center A U.S. Army Center of Excellence for Modeling and Simulation of Ground Vehicles Recent Advances in Reliability Estimation of Time-Dependent Problems Using the Concept of Composite
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 informationOpen Problems in Mixed Models
xxiii Determining how to deal with a not positive definite covariance matrix of random effects, D during maximum likelihood estimation algorithms. Several strategies are discussed in Section 2.15. For
More informationThe estimation methods in Nonlinear Mixed-Effects Models (NLMM) still largely rely on numerical approximation of the likelihood function
ABSTRACT Title of Dissertation: Diagnostics for Nonlinear Mixed-Effects Models Mohamed O. Nagem, Doctor of Philosophy, 2009 Dissertation directed by: Professor Benjamin Kedem Department of Mathematics
More informationA Simulation Study on Confidence Interval Procedures of Some Mean Cumulative Function Estimators
Statistics Preprints Statistics -00 A Simulation Study on Confidence Interval Procedures of Some Mean Cumulative Function Estimators Jianying Zuo Iowa State University, jiyizu@iastate.edu William Q. Meeker
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationUsing Accelerated Life Tests Results to Predict Product Field Reliability
Statistics Preprints Statistics 6-2008 Using Accelerated Life Tests Results to Predict Product Field Reliability William Q. Meeker Iowa State University, wqmeeker@iastate.edu Luis A. Escobar Louisiana
More informationMonte Carlo Studies. The response in a Monte Carlo study is a random variable.
Monte Carlo Studies The response in a Monte Carlo study is a random variable. The response in a Monte Carlo study has a variance that comes from the variance of the stochastic elements in the data-generating
More informationSubject CS1 Actuarial Statistics 1 Core Principles
Institute of Actuaries of India Subject CS1 Actuarial Statistics 1 Core Principles For 2019 Examinations Aim The aim of the Actuarial Statistics 1 subject is to provide a grounding in mathematical and
More informationAdvanced Statistics II: Non Parametric Tests
Advanced Statistics II: Non Parametric Tests Aurélien Garivier ParisTech February 27, 2011 Outline Fitting a distribution Rank Tests for the comparison of two samples Two unrelated samples: Mann-Whitney
More informationInstitute of Actuaries of India
Institute of Actuaries of India Subject CT3 Probability and Mathematical Statistics For 2018 Examinations Subject CT3 Probability and Mathematical Statistics Core Technical Syllabus 1 June 2017 Aim The
More informationResearch Article Applying Hierarchical Bayesian Neural Network in Failure Time Prediction
Mathematical Problems in Engineering Volume 2012, Article ID 953848, 11 pages doi:10.1155/2012/953848 Research Article Applying Hierarchical Bayesian Neural Network in Failure Time Prediction Ling-Jing
More informationOptimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests
International Journal of Performability Engineering, Vol., No., January 24, pp.3-4. RAMS Consultants Printed in India Optimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests N. CHANDRA *, MASHROOR
More informationSTA 216: GENERALIZED LINEAR MODELS. Lecture 1. Review and Introduction. Much of statistics is based on the assumption that random
STA 216: GENERALIZED LINEAR MODELS Lecture 1. Review and Introduction Much of statistics is based on the assumption that random variables are continuous & normally distributed. Normal linear regression
More informationA Tool for Evaluating Time-Varying-Stress Accelerated Life Test Plans with Log-Location- Scale Distributions
Statistics Preprints Statistics 6-2010 A Tool for Evaluating Time-Varying-Stress Accelerated Life Test Plans with Log-Location- Scale Distributions Yili Hong Virginia Tech Haiming Ma Iowa State University,
More informationSTAT 6350 Analysis of Lifetime Data. Failure-time Regression Analysis
STAT 6350 Analysis of Lifetime Data Failure-time Regression Analysis Explanatory Variables for Failure Times Usually explanatory variables explain/predict why some units fail quickly and some units survive
More informationOnline appendix to accompany:
Online appendix to accompany: Preacher, K. J., & Hancock, G. R. (submitted). Meaningful aspects of change as novel random coefficients: A general method for reparameterizing longitudinal models. Contents
More informationChapter 17: Undirected Graphical Models
Chapter 17: Undirected Graphical Models The Elements of Statistical Learning Biaobin Jiang Department of Biological Sciences Purdue University bjiang@purdue.edu October 30, 2014 Biaobin Jiang (Purdue)
More informationBayesian Inference in GLMs. Frequentists typically base inferences on MLEs, asymptotic confidence
Bayesian Inference in GLMs Frequentists typically base inferences on MLEs, asymptotic confidence limits, and log-likelihood ratio tests Bayesians base inferences on the posterior distribution of the unknowns
More informationLOGISTIC REGRESSION Joseph M. Hilbe
LOGISTIC REGRESSION Joseph M. Hilbe Arizona State University Logistic regression is the most common method used to model binary response data. When the response is binary, it typically takes the form of
More informationGauge Plots. Gauge Plots JAPANESE BEETLE DATA MAXIMUM LIKELIHOOD FOR SPATIALLY CORRELATED DISCRETE DATA JAPANESE BEETLE DATA
JAPANESE BEETLE DATA 6 MAXIMUM LIKELIHOOD FOR SPATIALLY CORRELATED DISCRETE DATA Gauge Plots TuscaroraLisa Central Madsen Fairways, 996 January 9, 7 Grubs Adult Activity Grub Counts 6 8 Organic Matter
More informationMA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems
MA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems Hierarchical Statistical Models for Complex Data Structures Motivation Basic hierarchical model and assumptions Statistical inference
More informationApproximation of Survival Function by Taylor Series for General Partly Interval Censored Data
Malaysian Journal of Mathematical Sciences 11(3): 33 315 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Approximation of Survival Function by Taylor
More informationWarranty Prediction Based on Auxiliary Use-rate Information
Statistics Preprints Statistics 2009 Warranty Prediction Based on Auxiliary Use-rate Information Yili Hong Virginia Polytechnic Institute and State University William Q. Meeker Iowa State University, wqmeeker@iastate.edu
More information20. REML Estimation of Variance Components. Copyright c 2018 (Iowa State University) 20. Statistics / 36
20. REML Estimation of Variance Components Copyright c 2018 (Iowa State University) 20. Statistics 510 1 / 36 Consider the General Linear Model y = Xβ + ɛ, where ɛ N(0, Σ) and Σ is an n n positive definite
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS330 / MAS83 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-0 8 Parametric models 8. Introduction In the last few sections (the KM
More informationE(Y ij b i ) = f(x ijβ i ), (13.1) β i = A i β + B i b i. (13.2)
1 Advanced topics 1.1 Introduction In this chapter, we conclude with brief overviews of several advanced topics. Each of these topics could realistically be the subject of an entire course! 1. Generalized
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS3301 / MAS8311 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-10 1 13 The Cox proportional hazards model 13.1 Introduction In the
More informationStatistical Prediction Based on Censored Life Data. Luis A. Escobar Department of Experimental Statistics Louisiana State University.
Statistical Prediction Based on Censored Life Data Overview Luis A. Escobar Department of Experimental Statistics Louisiana State University and William Q. Meeker Department of Statistics Iowa State University
More informationA Nonlinear Random Coefficients Model. for Degradation Testing
A Nonlinear Random Coefficients Model for Degradation Testing Suk Joo Bae Paul H. Kvam Georgia Institute of Technology February 8, 2004 Abstract As an alternative to traditional life testing, degradation
More information