8 RESPONSE SURFACE DESIGNS
|
|
- Caroline Dawson
- 6 years ago
- Views:
Transcription
1 8 RESPONSE SURFACE DESIGNS Desirable Properties of a Response Surface Design 1. It should generate a satisfactory distribution of information throughout the design region. 2. It should ensure that the fitted value ŷ at x be as close to the true value at x. 3. It should provide an estimate of the pure experimental error. 4. It should give good detectability of model lack of fit. 5. It should allow experiments to be performed in blocks. 6. It should allow for designs of increasing order to be constructed sequentially. 7. It should be efficient with respect to the number of experimental runs. 8. It should be cost effective. 9. It should provide a check of the homogeneity of variance assumption. 10. It should provide a good distribution of the scaled prediction variance Var[ŷ(x)]/σ 2 throughout the design region. 11. It should be robust to the presence of outliers in the data. 8.1 Prediction Variance Recall that the prediction variance at any point x in the design region is given by Var[ŷ(x)] = (13) where X is the model matrix corresponding to x (m) which is a vector corresponding to the model terms. For example: For the first-order model, x (m) = [1, x 1, x 2,..., x k ] For the interactive model, x (m) = [1, x 1, x 2,..., x k, x 1 x 2,..., x k 1 x k ] For the second-order model, x (m) = [1, x 1, x 2,..., x k, x 2 1, x 2 2,..., x 2 k, x 1 x 2,..., x k 1 x k ] The definition of Var[ŷ(x)] in (13) indicates that 1. Var[ŷ(x)] varies based on the location of x. 2. Var[ŷ(x)] is dependent on the choice of experimental design (because of (X X) 1 ). In design comparison studies, the scaled prediction variance function NVar[ŷ(x)] σ 2 = is often used because the division by σ 2 makes the quantity scale-free and the multiplication by the design size N allows this quantity to reflect variance on a per observation basis. 148
2 That is, when two designs are being compared, scaling by N penalizes the design with the larger sample size. Thus, emphasis is also placed on design size efficiency. EXAMPLE: Consider a k-factor experiment having N experimental runs. The design matrix D is the N k matrix whose rows correspond to the N experimental runs. For any p-parameter polynomial response surface model, the model matrix X associated with D is the N p matrix whose columns correspond to the p terms in the model. For the 3 2 design, (x 1, x 2 ) { 1, 0, 1} { 1, 0, 1}: Design Matrix D Model Matrix X (p = 6) k = 2 y = β 0 + β 1 x 1 + β 2 x 2 + β 12 x 1 x 2 N = 9 +β 11 x β 22 x ɛ x 1 x x 1 x 2 x 1 x 2 x 2 1 x The problem of choosing a best design D for fitting the parameters of a linear model can be interpreted in more than one way. Examples: D produces model coefficient estimates with smallest variance. X should have orthogonal columns. For second or higher-order polynomial response surface models, there is not a unique class of best designs (Box and Hunter 1957). Coefficient estimates should be studied simultaneously. Therefore, one desirable design property is to produce predicted values Ŷ (x) with small variance, i.e., small Var(Ŷ (x)). For example, consider two 6-point designs with 1 factor. The experiment consists of collecting data at temperatures between 80 and 100. Temperature is coded as
3 Design D1 X for D1 Design D2 X for D2 N = 6 y = β 0 + β 1 x + ɛ N = 6 y = β 0 + β 1 x + ɛ For D1 : [ ] [ ] 6 0 1/6 0 (X X) = (X X) 1 = /4 For D2 : [ ] [ ] 6 0 1/6 0 (X X) = (X X) 1 = 0 5/2 0 2/5 [ ] [ ] 1/6 0 1 For D1 : V (x) = 6 [ 1 x ] 0 1/4 x = [ ] [ ] 1/6 0 1 For D2 : V (x) = 6 [ 1 x ] 0 2/5 x = If we consider the second-order model then we can show that For D1 : V (x) = For D2 : V (x) = y = β 0 + β 1 x 1 + β 11 x ɛ : For k = 1 factor, V (x) can be plotted across the interval design space. For the first-order response surface model depicted in Figure 1a, V (x) is uniformly better for D1 than for D2. This is not true, however, for V (x) for the second-order response surface model shown in Figure 1b. V (x) is smaller for D1 than D2 near ±1 because D1 has replicated endpoints at ±1. V (x) is smaller for D2 than D1 near 0 because D2 has replicated points at 0 and at ±.5. Each design has its strengths and weaknesses with respect to the prediction variance of the second order response surface model. It is important to remember that the prediction variance function is dependent on the response surface design and the response surface model. 150
4 Figure 1a: First Order Model Figure 1b: Second Order Model Design D1 Design D2 For designs with k = 2 factors (x 1, x 2 ), V (x) can be displayed using a contour plot or a 3-dimensional surface plot. Example: For the N = 9 point 3 2 design with (x 1, x 2 ) { 1, 0, 1} { 1, 0, 1}: V (x) = A contour plot and a 3-dimensional surface plot of V (x) for the 3 2 design are shown in Figure 2a and Figure 2b. Figure 2a: Contour plot Figure 2b: 3-dimensional plot Graphical techniques for evaluating prediction variance properties throughout the experimental region have been developed for studying k 3 factors. These will be discussed later in the course. 151
5 152
6 8.2 Design of Experiments for First-Order Models Situation: An experimental design consisting of N runs is to be conducted on x 1, x 2,..., x k, a single response y is to be recorded, and a first-order model y i β 0 + β 1 x 1 + β 2 x β k x k is considered adequate. We will use coded design variables such that x j [ 1, 1] j. That is, the design region R(x) is a k-dimensional hypercube. Once the data is collected, we use the method of least squares to fit the first-order model ŷ = β 0 + k β i x i i=1 An orthogonal design or orthogonal array is an experimental design X = [1, x 1, x 2,..., x k ] where x j is the j th column of X is x i x j = 0 for all i j and 1 x j = 0 for all j. As a result, if two columns are orthogonal, then the two corresponding variables are linearly independent and the corresponding estimates are independent. The goal is to simultaneously minimize the variances of the b j s and, if possible, retain orthogonality among the columns of X. In other words, we want (X X) 1 to be a diagonal matrix whose diagonal entries are minimized. (Or equivalently, we want (X X) to be a diagonal matrix whose diagonal entries are maximized.) Variance Optimality Theorem: Suppose the first-order model ŷ = β 0 + k β i x i be fitted and the design size N is fixed. If x j [ 1, 1] for j = 1, 2,..., k, then var(b j ) is minimized if the design is orthogonal and all x i levels are ±1 for j = 1, 2,..., k. Outline of proof: i=1 is to The i th diagonal element of X X = N j=1 x2 ij. This will be will be maximized when each x ij = ±1, and will yield values on the diagonal = N. If the columns of X are also orthogonal, then X X is a diagonal matrix NI k+1. Thus, (X X) 1 = (1/N)I k+1, and the variance/covariance matrix σ 2 (X X) 1 = (σ 2 /N)I k+1. Thus, the diagonal entries will be minimized. Recall that the variance of any parameter estimate var(b j ) is a diagonal entry of σ 2 (X X) 1. Thus, for each j, the var(b j ) is minimized by designs satisfying the stated conditions. The two-level full-factorial 2 k designs, fractional-factorial 2 k p designs of at least Resolution III, and Plackett-Burman designs satisfy this theorem. A saturated design for a model with P parameters (excluding the intercept) is a design which allows estimation of all P model parameters in N = points. This implies there are no degrees of freedom for error. Therefore, saturated designs are used primarily for parameter estimation. 153
7 Another type of orthogonal first-order design is the simplex design. A simplex design is an orthogonal saturated design with N = k + 1 whose points represent the vertices of a regular-sided figure. In general, the design points are given by the rows of the matrix D = x 11 x 21 x k1 x 12 x 22 x k2 x 1,N x 2,N x k,n such that the angle that any two points makes with the origin is θ where cos(θ) = For k = 2, the points form an equilateral triangle with cos(θ) = 1/2 (i.e., θ = 120 ). For k = 3, the points form a regular tetrahedron (pyramid) with cos(θ) = 1/3. Rotation of a k-variable simplex design will yield another k-variable simplex design. To construct an k-variable simplex design with N = k + 1-points, you begin with any N N orthogonal matrix P with equal first column elements. The last N 1 columns of the matrix NP yield a simplex design. A matrix P can be obtained by selecting any nonsingular N N matrix Q with equal first column elements and then applying the Gram-Schmidt orthonormalization techniques to the columns of Q. A simpler method is to use the following table (from Montgomery s Design and Analysis of Experiments text). The table contains the coefficients of orthogonal polynomials to form Q and then apply the Gram-Schmidt orthonormalization technique. 154
8 155
9 156
10 A Three-Factor Central Composite Design 8.3 Central Composite Designs One of the most popular and commonly used classes of experimental designs for fitting the second order model are the central composite designs introduced by Box and Wilson (1951). Assuming k 2 design variables, the central composite design (CCD) consists of: (i) An f = 2 k p full (p = 0) or fractional (p > 0) factorial design of at least Resolution V. Each point is of the form (x 1,..., x k ) = (±1, ±1,..., ±1). (ii) 2k axial or star points of the form (x 1,..., x i,..., x k ) = (0,.., 0, ±α, 0,..., 0) for 1 i k. (iii) n c center points (x 1,..., x k ) = (0, 0,..., 0). If α = 1 for the axial points, then the design is referred to as a face-centered cube design. The CCD contains N = f + 2k + n c points to estimate the ( ) k+2 2 parameters in the second-order model. Each of the three types of points in a central composite design play different roles: The factorial points allow estimation of the first-order and interaction terms. The axial points allow estimation of the squared terms. The center points provide an internal estimate of pure error used to test for lack of fit and also contribute toward estimation of the squared terms. 157
11 Table of Central Composite Designs k Fraction Resolution Defining Factorial Axial Center Total Relation Points Points Points Points n c 8 + n c n c 14 + n c n c 24 + n c V E=ABCD n c 26 + n c VI F=ABCDE n c 44 + n c VII G=ABCDEF n c 78 + n c V G=CDEF H=ABEF n c 80 + n c VI H=CDEFG J=ABEFG n c n c Central Composite Designs (Face-Centered Cube) ============================================== k=3 k=4 k=5 ======== ============== ============== : : : : : : : : : : : : For a spherical design region, replace the (±1, 0, 0,..., 0), (0, ±1, 0,..., 0),..., (0, 0,..., ±1) axial points with (±α, 0, 0,..., 0), (0, ±α, 0,..., 0),..., (0, 0,..., ±α). The optimal choice of α will be discussed later. 158
12 8.4 Box-Behnken Designs A second class of experimental designs for quadratic regression is the class of Box-Behnken designs introduced by Box and Behnken (1960). Assuming k 3, most of the Box-Behnken designs (BBD) are constructed by combining two-level factorial designs with balanced incomplete block designs (BIBD). Associated with every BIBD, and hence, every BBD considered, are the following parameters: k = the number of design variables. b = the number of blocks in the BIBD. t = the number of design variables per block. r = the number of blocks in which a design variable appears. λ = the number of times that each pair of design variables appears in the same block. It must hold that λ =. When constructing a BBD: 1. The t columns defining a 2 t factorial design (or 2 t p fractional-factorial design for larger t) with levels ±1 replace the t design variables appearing in each block in the BIBD. 2. The remaining k t columns are set at mid-level n c mid-level center points (0,..., 0) are included in the design. The total BBD size is N = where f = 2 t. It should be noted that several of the designs (e.g., when k = 6) proposed by Box and Behnken (1960) are based on designs with partially balanced incomplete block designs (PBIBDs). Converting a 3-Factor BIBD into a 3-Variable BBD BIBD Factors Blocks A B C BBD Variables Blocks x 1 x 2 x CP
13 160
14 161
15 Box-Behnken Designs for k 7 162
Design Moments. Many properties of experimental designs are quantified by design moments. Given a model matrix 1 x 11 x 21 x k1 1 x 12 x 22 x k2 X =
8.5 Rotatability Recall: Var[ŷ(x)] = σ 2 x (m) () 1 x (m) is the prediction variance and NVar[ŷ(x)]/σ 2 = Nx (m) () 1 x (m) is the scaled prediction variance. A design is rotatable if the prediction variance
More informationAnalysis of Variance and Design of Experiments-II
Analysis of Variance and Design of Experiments-II MODULE VIII LECTURE - 36 RESPONSE SURFACE DESIGNS Dr. Shalabh Department of Mathematics & Statistics Indian Institute of Technology Kanpur 2 Design for
More informationResponse Surface Methodology
Response Surface Methodology Process and Product Optimization Using Designed Experiments Second Edition RAYMOND H. MYERS Virginia Polytechnic Institute and State University DOUGLAS C. MONTGOMERY Arizona
More informationResponse Surface Methodology:
Response Surface Methodology: Process and Product Optimization Using Designed Experiments RAYMOND H. MYERS Virginia Polytechnic Institute and State University DOUGLAS C. MONTGOMERY Arizona State University
More informationOptimal design of experiments
Optimal design of experiments Session 4: Some theory Peter Goos / 40 Optimal design theory continuous or approximate optimal designs implicitly assume an infinitely large number of observations are available
More informationExperimental designs for multiple responses with different models
Graduate Theses and Dissertations Graduate College 2015 Experimental designs for multiple responses with different models Wilmina Mary Marget Iowa State University Follow this and additional works at:
More informationResponse surface methodology
Advanced Review André I. Khuri 1 and Siuli Mukhopadhyay 2 The purpose of this article is to provide a survey of the various stages in the development of response surface methodology (RSM). The coverage
More informationUseful Numerical Statistics of Some Response Surface Methodology Designs
Journal of Mathematics Research; Vol. 8, No. 4; August 20 ISSN 19-9795 E-ISSN 19-9809 Published by Canadian Center of Science and Education Useful Numerical Statistics of Some Response Surface Methodology
More information22m:033 Notes: 7.1 Diagonalization of Symmetric Matrices
m:33 Notes: 7. Diagonalization of Symmetric Matrices Dennis Roseman University of Iowa Iowa City, IA http://www.math.uiowa.edu/ roseman May 3, Symmetric matrices Definition. A symmetric matrix is a matrix
More information2 Introduction to Response Surface Methodology
2 Introduction to Response Surface Methodology 2.1 Goals of Response Surface Methods The experimenter is often interested in 1. Finding a suitable approximating function for the purpose of predicting a
More informationResponse Surface Methodology IV
LECTURE 8 Response Surface Methodology IV 1. Bias and Variance If y x is the response of the system at the point x, or in short hand, y x = f (x), then we can write η x = E(y x ). This is the true, and
More informationA Study on Minimal-point Composite Designs
A Study on Minimal-point Composite Designs by Yin-Jie Huang Advisor Ray-Bing Chen Institute of Statistics, National University of Kaohsiung Kaohsiung, Taiwan 8 R.O.C. July 007 Contents Introduction Construction
More informationRESPONSE SURFACE MODELLING, RSM
CHEM-E3205 BIOPROCESS OPTIMIZATION AND SIMULATION LECTURE 3 RESPONSE SURFACE MODELLING, RSM Tool for process optimization HISTORY Statistical experimental design pioneering work R.A. Fisher in 1925: Statistical
More informationSIX SIGMA IMPROVE
SIX SIGMA IMPROVE 1. For a simplex-lattice design the following formula or equation determines: A. The canonical formula for linear coefficients B. The portion of each polynomial in the experimental model
More information7. Response Surface Methodology (Ch.10. Regression Modeling Ch. 11. Response Surface Methodology)
7. Response Surface Methodology (Ch.10. Regression Modeling Ch. 11. Response Surface Methodology) Hae-Jin Choi School of Mechanical Engineering, Chung-Ang University 1 Introduction Response surface methodology,
More informationInverses. Stephen Boyd. EE103 Stanford University. October 28, 2017
Inverses Stephen Boyd EE103 Stanford University October 28, 2017 Outline Left and right inverses Inverse Solving linear equations Examples Pseudo-inverse Left and right inverses 2 Left inverses a number
More informationFIVE-FACTOR CENTRAL COMPOSITE DESIGNS ROBUST TO A PAIR OF MISSING OBSERVATIONS. Munir Akhtar Vice-Chanselor, Islamia University Bahawalpur, Pakistan.
Journal of Research (Science), Bahauddin Zakariya University, Multan, Pakistan. Vol.2, No.2, December 2, pp. 5-5 ISSN 2-2 FIVE-FACTOR CENTRAL COMPOSITE DESIGNS ROBUST TO A PAIR OF MISSING OBSERVATIONS
More informationDESIGN OF EXPERIMENT ERT 427 Response Surface Methodology (RSM) Miss Hanna Ilyani Zulhaimi
+ DESIGN OF EXPERIMENT ERT 427 Response Surface Methodology (RSM) Miss Hanna Ilyani Zulhaimi + Outline n Definition of Response Surface Methodology n Method of Steepest Ascent n Second-Order Response Surface
More informationPractical Statistics for the Analytical Scientist Table of Contents
Practical Statistics for the Analytical Scientist Table of Contents Chapter 1 Introduction - Choosing the Correct Statistics 1.1 Introduction 1.2 Choosing the Right Statistical Procedures 1.2.1 Planning
More informationDesign & Analysis of Experiments 7E 2009 Montgomery
Chapter 5 1 Introduction to Factorial Design Study the effects of 2 or more factors All possible combinations of factor levels are investigated For example, if there are a levels of factor A and b levels
More informationProcess/product optimization using design of experiments and response surface methodology
Process/product optimization using design of experiments and response surface methodology Mikko Mäkelä Sveriges landbruksuniversitet Swedish University of Agricultural Sciences Department of Forest Biomaterials
More informationChemometrics Unit 4 Response Surface Methodology
Chemometrics Unit 4 Response Surface Methodology Chemometrics Unit 4. Response Surface Methodology In Unit 3 the first two phases of experimental design - definition and screening - were discussed. In
More informationLinear Models Review
Linear Models Review Vectors in IR n will be written as ordered n-tuples which are understood to be column vectors, or n 1 matrices. A vector variable will be indicted with bold face, and the prime sign
More informationUnit 12: Response Surface Methodology and Optimality Criteria
Unit 12: Response Surface Methodology and Optimality Criteria STA 643: Advanced Experimental Design Derek S. Young 1 Learning Objectives Revisit your knowledge of polynomial regression Know how to use
More informationBayesian Regressions in Experimental Design
Bayesian Regressions in Experimental Design Stanley Sawyer Washington University Vs. April 9, 8. Introduction. The purpose here is to derive a formula for the posterior probabilities given observed data
More informationWeighted Least Squares
Weighted Least Squares The standard linear model assumes that Var(ε i ) = σ 2 for i = 1,..., n. As we have seen, however, there are instances where Var(Y X = x i ) = Var(ε i ) = σ2 w i. Here w 1,..., w
More informationLINEAR ALGEBRA 1, 2012-I PARTIAL EXAM 3 SOLUTIONS TO PRACTICE PROBLEMS
LINEAR ALGEBRA, -I PARTIAL EXAM SOLUTIONS TO PRACTICE PROBLEMS Problem (a) For each of the two matrices below, (i) determine whether it is diagonalizable, (ii) determine whether it is orthogonally diagonalizable,
More informationPrinciple Components Analysis (PCA) Relationship Between a Linear Combination of Variables and Axes Rotation for PCA
Principle Components Analysis (PCA) Relationship Between a Linear Combination of Variables and Axes Rotation for PCA Principle Components Analysis: Uses one group of variables (we will call this X) In
More informationLecture 7. Econ August 18
Lecture 7 Econ 2001 2015 August 18 Lecture 7 Outline First, the theorem of the maximum, an amazing result about continuity in optimization problems. Then, we start linear algebra, mostly looking at familiar
More informationMA 575 Linear Models: Cedric E. Ginestet, Boston University Mixed Effects Estimation, Residuals Diagnostics Week 11, Lecture 1
MA 575 Linear Models: Cedric E Ginestet, Boston University Mixed Effects Estimation, Residuals Diagnostics Week 11, Lecture 1 1 Within-group Correlation Let us recall the simple two-level hierarchical
More informationExperimental Design and Optimization
. Experimental Design Stages a) Identifying the factors which may affect the results of an experiment; b) Designing the experiment so that the effects of uncontrolled factors are minimized; c) Using statistical
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationSolutions to Review Problems for Chapter 6 ( ), 7.1
Solutions to Review Problems for Chapter (-, 7 The Final Exam is on Thursday, June,, : AM : AM at NESBITT Final Exam Breakdown Sections % -,7-9,- - % -9,-,7,-,-7 - % -, 7 - % Let u u and v Let x x x x,
More informationLeverage. the response is in line with the other values, or the high leverage has caused the fitted model to be pulled toward the observed response.
Leverage Some cases have high leverage, the potential to greatly affect the fit. These cases are outliers in the space of predictors. Often the residuals for these cases are not large because the response
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationWeighted Least Squares
Weighted Least Squares The standard linear model assumes that Var(ε i ) = σ 2 for i = 1,..., n. As we have seen, however, there are instances where Var(Y X = x i ) = Var(ε i ) = σ2 w i. Here w 1,..., w
More informationResponse Surface Methodology: Process and Product Optimization Using Designed Experiments, 3rd Edition
Brochure More information from http://www.researchandmarkets.com/reports/705963/ Response Surface Methodology: Process and Product Optimization Using Designed Experiments, 3rd Edition Description: Identifying
More informationA COMPARISON STUDY OF SECOND-ORDER SCREENING DESIGNS AND THEIR EXTENSION DISSERTATION. Shane A. Dougherty, Lt Col AFIT-ENS-DS-13-D-01
A COMPARISON STUDY OF SECOND-ORDER SCREENING DESIGNS AND THEIR EXTENSION DISSERTATION Shane A. Dougherty, Lt Col AFIT-ENS-DS-13-D-01 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY AIR FORCE INSTITUTE OF TECHNOLOGY
More informationVectors and Matrices Statistics with Vectors and Matrices
Vectors and Matrices Statistics with Vectors and Matrices Lecture 3 September 7, 005 Analysis Lecture #3-9/7/005 Slide 1 of 55 Today s Lecture Vectors and Matrices (Supplement A - augmented with SAS proc
More informationREVIEW OF EXPERIMENTAL DESIGN IN ANALYTICAL CHEMISTRY
Page55 Indo American Journal of Pharmaceutical Research, 217 ISSN NO: 2231-6876 REVIEW OF EXPERIMENTAL DESIGN IN ANALYTICAL CHEMISTRY T.Sudha *, G.Divya, J.Sujaritha, P. Duraimurugan Department of Pharmaceutical
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationExplaining Correlations by Plotting Orthogonal Contrasts
Explaining Correlations by Plotting Orthogonal Contrasts Øyvind Langsrud MATFORSK, Norwegian Food Research Institute. www.matforsk.no/ola/ To appear in The American Statistician www.amstat.org/publications/tas/
More informationAppendix IV Experimental Design
Experimental Design The aim of pharmaceutical formulation and development is to develop an acceptable pharmaceutical formulation in the shortest possible time, using minimum number of working hours and
More informationSection 6.4. The Gram Schmidt Process
Section 6.4 The Gram Schmidt Process Motivation The procedures in 6 start with an orthogonal basis {u, u,..., u m}. Find the B-coordinates of a vector x using dot products: x = m i= x u i u i u i u i Find
More informationLinear Analysis Lecture 16
Linear Analysis Lecture 16 The QR Factorization Recall the Gram-Schmidt orthogonalization process. Let V be an inner product space, and suppose a 1,..., a n V are linearly independent. Define q 1,...,
More informationTWO-LEVEL FACTORIAL EXPERIMENTS: IRREGULAR FRACTIONS
STAT 512 2-Level Factorial Experiments: Irregular Fractions 1 TWO-LEVEL FACTORIAL EXPERIMENTS: IRREGULAR FRACTIONS A major practical weakness of regular fractional factorial designs is that N must be a
More informationExperimental design. Matti Hotokka Department of Physical Chemistry Åbo Akademi University
Experimental design Matti Hotokka Department of Physical Chemistry Åbo Akademi University Contents Elementary concepts Regression Validation Design of Experiments Definitions Random sampling Factorial
More informationStatistical View of Least Squares
May 23, 2006 Purpose of Regression Some Examples Least Squares Purpose of Regression Purpose of Regression Some Examples Least Squares Suppose we have two variables x and y Purpose of Regression Some Examples
More informationPreface Introduction to Statistics and Data Analysis Overview: Statistical Inference, Samples, Populations, and Experimental Design The Role of
Preface Introduction to Statistics and Data Analysis Overview: Statistical Inference, Samples, Populations, and Experimental Design The Role of Probability Sampling Procedures Collection of Data Measures
More information1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal
. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal 3 9 matrix D such that A = P DP, for A =. 3 4 3 (a) P = 4, D =. 3 (b) P = 4, D =. (c) P = 4 8 4, D =. 3 (d) P
More informationDefinitive Screening Designs
Definitive Screening Designs Bradley Jones September 2011 Copyright 2008, SAS Institute Inc. All rights reserved. Joint work with Chris Nachtsheim Outline 1. Motivation 2. Design Structure 3. Design Construction
More informationLecture 3: QR-Factorization
Lecture 3: QR-Factorization This lecture introduces the Gram Schmidt orthonormalization process and the associated QR-factorization of matrices It also outlines some applications of this factorization
More informationClasses of Second-Order Split-Plot Designs
Classes of Second-Order Split-Plot Designs DATAWorks 2018 Springfield, VA Luis A. Cortés, Ph.D. The MITRE Corporation James R. Simpson, Ph.D. JK Analytics, Inc. Peter Parker, Ph.D. NASA 22 March 2018 Outline
More information235 Final exam review questions
5 Final exam review questions Paul Hacking December 4, 0 () Let A be an n n matrix and T : R n R n, T (x) = Ax the linear transformation with matrix A. What does it mean to say that a vector v R n is an
More informationMATRIX ALGEBRA. or x = (x 1,..., x n ) R n. y 1 y 2. x 2. x m. y m. y = cos θ 1 = x 1 L x. sin θ 1 = x 2. cos θ 2 = y 1 L y.
as Basics Vectors MATRIX ALGEBRA An array of n real numbers x, x,, x n is called a vector and it is written x = x x n or x = x,, x n R n prime operation=transposing a column to a row Basic vector operations
More informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More information2. Matrix Algebra and Random Vectors
2. Matrix Algebra and Random Vectors 2.1 Introduction Multivariate data can be conveniently display as array of numbers. In general, a rectangular array of numbers with, for instance, n rows and p columns
More informationMA 265 FINAL EXAM Fall 2012
MA 265 FINAL EXAM Fall 22 NAME: INSTRUCTOR S NAME:. There are a total of 25 problems. You should show work on the exam sheet, and pencil in the correct answer on the scantron. 2. No books, notes, or calculators
More informationResponse Surface Methodology
Response Surface Methodology Bruce A Craig Department of Statistics Purdue University STAT 514 Topic 27 1 Response Surface Methodology Interested in response y in relation to numeric factors x Relationship
More information2. Review of Linear Algebra
2. Review of Linear Algebra ECE 83, Spring 217 In this course we will represent signals as vectors and operators (e.g., filters, transforms, etc) as matrices. This lecture reviews basic concepts from linear
More informationDesign of Engineering Experiments Chapter 5 Introduction to Factorials
Design of Engineering Experiments Chapter 5 Introduction to Factorials Text reference, Chapter 5 page 170 General principles of factorial experiments The two-factor factorial with fixed effects The ANOVA
More informationPharmaceutical Experimental Design and Interpretation
Pharmaceutical Experimental Design and Interpretation N. ANTHONY ARMSTRONG, B. Pharm., Ph.D., F.R.Pharm.S., MCPP. KENNETH C. JAMES, M. Pharm., Ph.D., D.Sc, FRSC, F.R.Pharm.S., C.Chem. Welsh School of Pharmacy,
More informationENGG5781 Matrix Analysis and Computations Lecture 8: QR Decomposition
ENGG5781 Matrix Analysis and Computations Lecture 8: QR Decomposition Wing-Kin (Ken) Ma 2017 2018 Term 2 Department of Electronic Engineering The Chinese University of Hong Kong Lecture 8: QR Decomposition
More informationTwo-Stage Computing Budget Allocation. Approach for Response Surface Method PENG JI
Two-Stage Computing Budget Allocation Approach for Response Surface Method PENG JI NATIONAL UNIVERSITY OF SINGAPORE 2005 Two-Stage Computing Budget Allocation Approach for Response Surface Method PENG
More informationECE521 week 3: 23/26 January 2017
ECE521 week 3: 23/26 January 2017 Outline Probabilistic interpretation of linear regression - Maximum likelihood estimation (MLE) - Maximum a posteriori (MAP) estimation Bias-variance trade-off Linear
More informationk is a product of elementary matrices.
Mathematics, Spring Lecture (Wilson) Final Eam May, ANSWERS Problem (5 points) (a) There are three kinds of elementary row operations and associated elementary matrices. Describe what each kind of operation
More informationLinear Methods for Regression. Lijun Zhang
Linear Methods for Regression Lijun Zhang zlj@nju.edu.cn http://cs.nju.edu.cn/zlj Outline Introduction Linear Regression Models and Least Squares Subset Selection Shrinkage Methods Methods Using Derived
More informationResponse Surface Methodology III
LECTURE 7 Response Surface Methodology III 1. Canonical Form of Response Surface Models To examine the estimated regression model we have several choices. First, we could plot response contours. Remember
More informationOrthogonal Transformations
Orthogonal Transformations Tom Lyche University of Oslo Norway Orthogonal Transformations p. 1/3 Applications of Qx with Q T Q = I 1. solving least squares problems (today) 2. solving linear equations
More informationDesign of Engineering Experiments Part 5 The 2 k Factorial Design
Design of Engineering Experiments Part 5 The 2 k Factorial Design Text reference, Special case of the general factorial design; k factors, all at two levels The two levels are usually called low and high
More informationObjective Experiments Glossary of Statistical Terms
Objective Experiments Glossary of Statistical Terms This glossary is intended to provide friendly definitions for terms used commonly in engineering and science. It is not intended to be absolutely precise.
More informationDesign and Analysis of Experiments
Design and Analysis of Experiments Part IX: Response Surface Methodology Prof. Dr. Anselmo E de Oliveira anselmo.quimica.ufg.br anselmo.disciplinas@gmail.com Methods Math Statistics Models/Analyses Response
More informationMATH 1120 (LINEAR ALGEBRA 1), FINAL EXAM FALL 2011 SOLUTIONS TO PRACTICE VERSION
MATH (LINEAR ALGEBRA ) FINAL EXAM FALL SOLUTIONS TO PRACTICE VERSION Problem (a) For each matrix below (i) find a basis for its column space (ii) find a basis for its row space (iii) determine whether
More informationContents. Response Surface Designs. Contents. References.
Response Surface Designs Designs for continuous variables Frédéric Bertrand 1 1 IRMA, Université de Strasbourg Strasbourg, France ENSAI 3 e Année 2017-2018 Setting Visualizing the Response First-Order
More information5. Orthogonal matrices
L Vandenberghe EE133A (Spring 2017) 5 Orthogonal matrices matrices with orthonormal columns orthogonal matrices tall matrices with orthonormal columns complex matrices with orthonormal columns 5-1 Orthonormal
More informationProbability and Statistics Notes
Probability and Statistics Notes Chapter Seven Jesse Crawford Department of Mathematics Tarleton State University Spring 2011 (Tarleton State University) Chapter Seven Notes Spring 2011 1 / 42 Outline
More informationDiagonalizing Matrices
Diagonalizing Matrices Massoud Malek A A Let A = A k be an n n non-singular matrix and let B = A = [B, B,, B k,, B n ] Then A n A B = A A 0 0 A k [B, B,, B k,, B n ] = 0 0 = I n 0 A n Notice that A i B
More information1 Last time: least-squares problems
MATH Linear algebra (Fall 07) Lecture Last time: least-squares problems Definition. If A is an m n matrix and b R m, then a least-squares solution to the linear system Ax = b is a vector x R n such that
More informationFACTOR SCREENING AND RESPONSE SURFACE EXPLORATION
Statistica Sinica 11(2001), 553-604 FACTOR SCREENING AND RESPONSE SURFACE EXPLORATION Shao-Wei Cheng and C. F. J. Wu Academia Sinica and University of Michigan Abstract: Standard practice in response surface
More informationMATH 304 Linear Algebra Lecture 34: Review for Test 2.
MATH 304 Linear Algebra Lecture 34: Review for Test 2. Topics for Test 2 Linear transformations (Leon 4.1 4.3) Matrix transformations Matrix of a linear mapping Similar matrices Orthogonality (Leon 5.1
More informationLinear Regression. In this problem sheet, we consider the problem of linear regression with p predictors and one intercept,
Linear Regression In this problem sheet, we consider the problem of linear regression with p predictors and one intercept, y = Xβ + ɛ, where y t = (y 1,..., y n ) is the column vector of target values,
More information1. Vectors.
1. Vectors 1.1 Vectors and Matrices Linear algebra is concerned with two basic kinds of quantities: vectors and matrices. 1.1 Vectors and Matrices Scalars and Vectors - Scalar: a numerical value denoted
More informationChapter 6: Orthogonality
Chapter 6: Orthogonality (Last Updated: November 7, 7) These notes are derived primarily from Linear Algebra and its applications by David Lay (4ed). A few theorems have been moved around.. Inner products
More informationThere are two things that are particularly nice about the first basis
Orthogonality and the Gram-Schmidt Process In Chapter 4, we spent a great deal of time studying the problem of finding a basis for a vector space We know that a basis for a vector space can potentially
More informationMath 423/533: The Main Theoretical Topics
Math 423/533: The Main Theoretical Topics Notation sample size n, data index i number of predictors, p (p = 2 for simple linear regression) y i : response for individual i x i = (x i1,..., x ip ) (1 p)
More informationLinear Least Squares Problems
Linear Least Squares Problems Introduction We have N data points (x 1,y 1 ),...(x N,y N ). We assume that the data values are given by y j = g(x j ) + e j, j = 1,...,N where g(x) = c 1 g 1 (x) + + c n
More informationDisegno Sperimentale (DoE) come strumento per QbD
Giornata di Studio Disegno Sperimentale (DoE) come strumento per QbD Università degli Studi di Milano Dipartimento di Scienze Farmaceutiche Milano, 22 aprile 2013 Dr. Lorenza Broccardo Introduction Nowadays,
More informationFactor Analysis (10/2/13)
STA561: Probabilistic machine learning Factor Analysis (10/2/13) Lecturer: Barbara Engelhardt Scribes: Li Zhu, Fan Li, Ni Guan Factor Analysis Factor analysis is related to the mixture models we have studied.
More informationMATH602: APPLIED STATISTICS
MATH602: APPLIED STATISTICS Dr. Srinivas R. Chakravarthy Department of Science and Mathematics KETTERING UNIVERSITY Flint, MI 48504-4898 Lecture 10 1 FRACTIONAL FACTORIAL DESIGNS Complete factorial designs
More informationDesign and Analysis of Experiments Prof. Jhareshwar Maiti Department of Industrial and Systems Engineering Indian Institute of Technology, Kharagpur
Design and Analysis of Experiments Prof. Jhareshwar Maiti Department of Industrial and Systems Engineering Indian Institute of Technology, Kharagpur Lecture 51 Plackett Burman Designs Hello, welcome. We
More informationMTH 2032 SemesterII
MTH 202 SemesterII 2010-11 Linear Algebra Worked Examples Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education December 28, 2011 ii Contents Table of Contents
More information1. Background: The SVD and the best basis (questions selected from Ch. 6- Can you fill in the exercises?)
Math 35 Exam Review SOLUTIONS Overview In this third of the course we focused on linear learning algorithms to model data. summarize: To. Background: The SVD and the best basis (questions selected from
More informationCHAPTER 4 EXPERIMENTAL DESIGN. 4.1 Introduction. Experimentation plays an important role in new product design, manufacturing
CHAPTER 4 EXPERIMENTAL DESIGN 4.1 Introduction Experimentation plays an important role in new product design, manufacturing process development and process improvement. The objective in all cases may be
More information[y i α βx i ] 2 (2) Q = i=1
Least squares fits This section has no probability in it. There are no random variables. We are given n points (x i, y i ) and want to find the equation of the line that best fits them. We take the equation
More informationThe value of a problem is not so much coming up with the answer as in the ideas and attempted ideas it forces on the would be solver I.N.
Math 410 Homework Problems In the following pages you will find all of the homework problems for the semester. Homework should be written out neatly and stapled and turned in at the beginning of class
More information18.S096 Problem Set 3 Fall 2013 Regression Analysis Due Date: 10/8/2013
18.S096 Problem Set 3 Fall 013 Regression Analysis Due Date: 10/8/013 he Projection( Hat ) Matrix and Case Influence/Leverage Recall the setup for a linear regression model y = Xβ + ɛ where y and ɛ are
More informationSparse Legendre expansions via l 1 minimization
Sparse Legendre expansions via l 1 minimization Rachel Ward, Courant Institute, NYU Joint work with Holger Rauhut, Hausdorff Center for Mathematics, Bonn, Germany. June 8, 2010 Outline Sparse recovery
More information11 Hypothesis Testing
28 11 Hypothesis Testing 111 Introduction Suppose we want to test the hypothesis: H : A q p β p 1 q 1 In terms of the rows of A this can be written as a 1 a q β, ie a i β for each row of A (here a i denotes
More informationConstruction and analysis of Es 2 efficient supersaturated designs
Construction and analysis of Es 2 efficient supersaturated designs Yufeng Liu a Shiling Ruan b Angela M. Dean b, a Department of Statistics and Operations Research, Carolina Center for Genome Sciences,
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 information