7 The Analysis of Response Surfaces
|
|
- Benedict Newman
- 6 years ago
- Views:
Transcription
1 7 The Analysis of Response Surfaces Goal: the researcher is seeking the experimental conditions which are most desirable, that is, determine optimum design variable levels. Once the researcher has determined the experimental design region by selecting ranges of the experimental variables, an experimental design is chosen and the experimental data is collected. The coefficients in the model are then estimated. We are considering the fitted second-order model or the fitted response surface: ŷ = b + k b i x i i= + k b ij x i x j + i<j k b ii x 2 i (9) i= where the b s are the estimated coefficients. Eq (9) is used to predict the response y for given x, x 2,..., x k and enables the researcher to conduct an analysis of the fitted response surface. In matrix notation Eq (9) is given by where x = x x 2 x k b = Note: B is a symmetric matrix. ŷ = () b b b 2 /2 b k /2 b 2 b 22 b 2k /2 B = b b k,k /2 k b kk Suppose the goal of the experimenter is to estimate the conditions on x, x 2,..., x k which maximize the response y. Then, the maximum, if it exists, will be a set of conditions on (x, x 2,..., x k ) such that the partial derivatives are simultaneously zero. ŷ/ x, ŷ/ x 2,..., ŷ/ x k This point, say x s = [x,, x 2,,..., x k, ], is called the stationary point of the fitted response surface. We apply the rules for differentiating with respect to a column vector and then equate the result to zero: ŷ = ] [b + x b + x Bx = x x Solving for x yields the stationary point x s = This point may or may not be the point which maximizes the response. The point falls into one of three categories:. x s is the point at which the response surface attains a maximum. 2. x s is the point at which the response surface attains a minimum. 3. x s is a saddle point. That is, the response could either increase or decrease as you move away from x s (depending on what direction you move). 4
2 7. Canonical Analysis The goal is to determine the nature of the stationary point and the entire response surface. The nature of the stationary point is determined by the signs of the eigenvalues of B. The type of response surface is determined by examining the magnitudes of the eigenvalues. Let P be the k k matrix whose columns are the normalized eigenvectors associated with the eigenvalues λ, λ 2,..., λ k of B. Then where Λ = Diag(λ, λ 2,..., λ k ) is a diagonal matrix of the eigenvalues of B. To understand a response surface, it is necessary to translate the matrix form in () into another form which is easier to interpret. The analysis begins with a translation of response function in (9) from the origin (x, x 2,..., x k ) = (,,..., ) to a new origin at the stationary point x s. Let z = (z, z 2,..., z k ) be the vector of re-centered x coordinates: z = x x s (Shift) Now we want to rotate the axes to align with the contours of the response surface. We do this by use of the eigenvalues λ, λ 2,..., λ k and the matrix of eigenvectors P. The response function is expressed in terms of new variables w, w 2,..., w k which correspond to the principal axes. That is, rotate the original axes about x, x 2,..., x k by transforming z to w = ( w, w 2,..., w k ) via the matrix of normalized eigenvectors: w = P z (Rotation) 5
3 In terms of predicted values, this translation yields ŷ = b + x b + x Bx where ŷ s is the predicted response at the stationary point x s. Because PP = I, we see P = P, and because w = P w, we have z = Pw. Therefore, ŷ = ŷ s + z Bz The w-axes are the principal axes of the contour system. The variables w, w 2,..., w k are called the canonical variables. The canonical form of the response surface is obtained by rewriting Eq () as: () ŷ = ŷ s + k λ i wi 2 (2) i= If λ i >, then ŷ will increase as you move along w i in either direction from x s. If λ i <, then ŷ will decrease as you move along w i in either direction from x s. The analysis after reduction of the response surface to canonical form is called a canonical analysis. 6
4 7
5 7.2 Interpretation of the Response Surface The canonical form in (2) describes the nature of the stationary point and the nature of the response surface about the stationary point:. If λ, λ 2,..., λ k are all negative, then the stationary point is a point of. k Note that if all λ i < then λ i wi 2 < for any nonzero w. Thus, if we move away from x s in any direction, ŷ decreases. i= 2. If λ, λ 2,..., λ k are all positive, then the stationary point is a point of. k Note that if all λ i > then λ i wi 2 > for any nonzero w. Thus, if we move away from x s in any direction, ŷ increases. i= 3. If λ, λ 2,..., λ k are mixed in sign, then the stationary point is a. Note k that if the λ i are mixed in sign then λ i wi 2 can be positive or negative. Thus, an increase or decrease in ŷ depends upon the direction traveled away from x s. We now consider three other types of response surfaces: i=. The stationary ridge system indicates that there is a region, rather than a single stationary point, where there is approximately a maximum or minimum response. Figure a shows a maximum stationary ridge along w and Figure b shows a minimum stationary ridge along w. 8
6 2. The rising ridge system has its stationary point remote from the experimental region and the estimated response increases as you move along the axis toward the stationary point. See the dotted line in the Figure 2a. 3. The falling ridge system has its stationary point remote from the experimental region and the estimated response decreases as you move along the axis toward the stationary point. See the dotted line in the Figure 2b. To determine whether a stationary, rising, or falling ridge describes the response surface in the experimental region, we need to simultaneously study the location of the stationary point and the magnitude of the eigenvalues. The magnitude of the eigenvalues provides information about the shape of the response surface. This is best exemplified in response surfaces involving k = 2 variables (and, hence, two eigenvalues): Case I: Suppose λ and λ 2 are both negative and λ 2 is considerably greater than λ. See Figure 3a. The stationary point is a maximum but note the sensitivity of ŷ with respect to the two canonical variables w and w 2. As you move along the w axis away from the stationary point in either direction, there is less of a change (decrease) in the response relative to the change (decrease) that occurs when moving a similar distance along the w 2 axis away from the stationary point in either direction. Case II: Suppose λ and λ 2 are both positive and λ 2 is considerably greater than λ. See Figure 3b. The stationary point is a minimum but note the sensitivity of ŷ with respect to the two canonical variables w and w 2. As you move along the w axis away from the stationary point in either direction, there is less of a change (increase) in the response relative to the change (increase) that occurs when moving a similar distance along the w 2 axis away from the stationary point in either direction. In other words, for Cases I & II, the response surface is elongated in the w direction and is insensitive to small changes in movement along the w axis. 9
7 Case III: Suppose λ and λ 2 have different signs such that λ < and λ 2 >, and λ 2 is considerably greater than λ. See Figure 4a. The stationary point is a saddle point but note the sensitivity of ŷ with respect to the two canonical variables w and w 2. As you move along the w axis away from the stationary point in either direction, there is less of a change (decrease) in the response relative to the change (increase) that occurs when moving a similar distance along the w 2 axis away from the stationary point in either direction. Case IV: Suppose λ and λ 2 have different signs such that λ <, and λ 2 > and λ 2 is considerably smaller than λ. See Figure 4b. The stationary point is a saddle point but note the sensitivity of ŷ with respect to the two canonical variables w and w 2. As you move along the w 2 axis away from the stationary point in either direction, there is less of a change (increase) in the response relative to the change (decrease) that occurs when moving a similar distance along the w axis away from the stationary point in either direction. In general, for any response surface with k eigenvalues, if λ j is considerably greater than λ i for two eigenvalues λ i and λ j, the response surface contours are elongated along w i and are narrower along w j. In general, for any response surface with k eigenvalues, if λ i λ j for two eigenvalues λ i and λ j, the response surface contours are nearly evenly spaced along w i and w j. That is, change in ŷ along w change in ŷ along w 2 when moving similar distances along the axes away from the stationary point. Near-Stationary Ridges The most extreme case with two variables is the stationary ridge. An exact stationary ridge occurs when one of the eigenvalues is zero. For the stationary ridges shown in Figures a and b, we have λ =. Consider the stationary ridge systems in Figures A and B. If λ, then ŷ = ŷ s + λ w 2 + λ 2 w 2 2 ŷ s + λ w 2 2. Therefore, changing w will have little effect on the predicted value ŷ. In practice, if one of the eigenvalues is very small, then the response surface approximates a stationary ridge. 2
8 Mathematically, the stationary ridge condition arises as a limiting eigenvalue case from a maximum, minimum, or saddle point condition. Thus, it is unlikely that an exact stationary ridge condition (some λ i = ) would occur in practice. The practical implications of a near-stationary ridge (some λ i researcher. ) are important to the For example, if the stationary point is a maximum and λ is near zero, then for all practical purposes the stationary point does not clearly define optimum operating conditions. In such circumstances, the researcher has a range of potential operating conditions along the w axis, all of which yield a near-maximum response. Approximate Rising and Falling Ridges Suppose that for k = 2 variables that the stationary point is not near the experimental region and that λ < and λ 2 is near zero. The response surface approaches a rising ridge condition. See Figure 5a. The boxed area is the experimental design region. Figure 5a Suppose that for k = 2 variables that the stationary point is not near the experimental region and that λ > and λ 2 is near zero. The response surface approaches a falling ridge condition. See Figure 5b. The boxed area is the experimental design region. Figure 5b 2
9 7.2. Three Examples of Data Analysis EXAMPLE : Two types of fertilizers, a standard N-P-K combination and the other a nutritional supplement, were applied to experimental plots to assess the effects on the yield of peanuts measured in pounds per plot. The level of the amount (lb/plot) of each fertilizer applied to a plot was determined by the coordinate settings of a central composite rotatable design. The data below represent the harvested yields of peanuts taken from two replicates of each experimental combination. (Note: the uncoded fertilizer levels corresponding to coded 2 levels are rounded to decimal place.) Fertilizer Fertilizer Coded Coded Yields (lb/plot) 2 x x 2 Rep Rep DM LOG; CLEAR; OUT; CLEAR; ; ODS PRINTER PDF file= C:\COURSES\ST578\SAS\canon.pdf ; ODS LISTING; OPTIONS LS=72 NODATE NONUMBER; **********************; ***** EXAMPLE # *****; **********************; DATA IN; INPUT F F2 LINES; ; PROC RSREG DATA=IN PLOTS=ALL ; MODEL YIELD = F F2 / LACKFIT ; TITLE EXAMPLE -- USING UNCODED DATA IN PROC RSREG ; RUN; 22
10 23
11 EXAMPLE -- USING UNCODED DATA IN PROC RSREG..5 Fit Diagnostics for YIELD RStudent - RStudent Predicted Value Predicted Value Leverage - YIELD 2 5 Cook's D Quantile Predicted Value 5 5 Observation 3 Fit Mean Percent Observations Parameters Error DF MSE R-Square Adj R-Square Proportion Less 24
12 . Plots for YIELD EXAMPLE -- USING UNCODED DATA IN PROC RSREG.5 Plots for YIELD F F F F2 25
13 EXAMPLE 2: Data taken from Response Surface Methodology by Raymond Myers, pages In this experiment, the researcher wants to gain an insight into the influence of sealing temperature (x ), cooling bar temperature (x 2 ), and % polyethylene additive (x 3 ) on the seal strength in grams per inch of a breadwrapper stock. The actual levels of the variables are coded as follows: x = seal temp x 2 = cooling temp 55 9 On the uncoded scale, five levels of each variable were used: x 3 = % polyethylene..6 The design variables and responses are: x x x x x 2 x 3 y x x 2 x 3 y From the canonical analysis, the stationary point is x s = (x, x 2, x 3 ) = (.,.26,.68). The uncoded levels are ξ = (3)(.) = Seal Temperature ξ 2 = 55 + (9)(.26) = 57.3 Cool Temperature ξ 3 =. + (.6)(.68) =.5 Polyethelene Additive SAS Code for Example 2 ******************************************************; *** EXAMPLE # FACTOR CCD -- MYERS PP ***; ******************************************************; DATA DSGN; INPUT X X2 X3 LINES; ; PROC RSREG DATA=DSGN PLOTS=(ALL SURFACE(AT(X= ))); MODEL STRENGTH = X X2 X3 / LACKFIT ; TITLE EXAMPLE 2 -- USING CODED DATA IN PROC RSREG ; RUN; 26
14 27
15 EXAMPLE 2 -- USING CODED DATA IN PROC RSREG Fit Diagnostics for STRENGTH RStudent - RStudent Predicted Value Predicted Value Leverage 2 - STRENGTH Cook's D Quantile Predicted Value Observation 4 Fit Mean Percent Observations Parameters Error DF MSE R-Square Adj R-Square Proportion Less 28
16 EXAMPLE 2 -- USING CODED DATA IN PROC RSREG EXAMPLE 2 -- USING CODED DATA IN PROC RSREG 29
17 EXAMPLE 3: From Response Surfaces by Andre Khuri and John Cornell, pages In this experiment, the researcher wants to investigate the effects of three fertilizer ingredients on the yield of snap beans under field conditions. The fertilizer ingredients and actual amounts applied were (i) nitrogen (N) from.94 to 6.29 lb/plt, (ii) phosphoric acid (P 2 O 5 ) from.59 to 2.97 lb/plot, and (iii) potash (K 2 O) from.6 to 4.22 lb/plot. The response of interest is the average yield in pounds per plot of snap beans. When coded the levels of nitrogen, phosphoric acid, and potash are: x = N 3.62 x 2 = P 2O 5.78 x 3 = K 2O On the uncoded scale, five levels of each variable were used: The design variables and responses are: x x x x x 2 x 3 y x x 2 x 3 y x x 2 x 3 y SAS Code for Example 3: DM LOG; CLEAR; OUT; CLEAR; ; * ODS PRINTER PDF file= C:\COURSES\ST578\SAS\canon3p.pdf ; ODS LISTING; OPTIONS PS=54 LS=72 NODATE NONUMBER; TITLE 3 FACTOR CCD -- KHURI AND CORNELL (PAGE 9) ; DATA DSGN; INPUT X X2 X3 LINES; ; PROC RSREG DATA=DSGN PLOTS=(ALL SURFACE(AT(x2= ))); MODEL YIELD = X X2 X3 / LACKFIT ; RUN; 3
18 3 FACTOR CCD -- KHURI AND CORNELL (PAGE 9) Coding Coefficients for the Independent Variables Factor Subtracted off Divided by X.682 X2.682 X3.682 Response Surface for Variable YIELD Response Mean.98 3 FACTOR CCD -- KHURI AND CORNELL (PAGE 9) Root MSE R-Square.786 Coefficient of Variation Parameter Estimate from Type Standard I Sum Coded Parameter Regression DF Estimate DF of Squares ErrorR-Square t Value F Pr Value > t Pr > Data F Intercept Linear < X Quadratic X2 Crossproduct X3 Total Model X*X X2*X Sum of X2*X2 DF Squares.2623 Mean Square 2.4 F.576 Value Pr.5955 > F X3*X Lack of Fit X3*X2 Pure Error X3*X3 Total Error Factor DF Sum of Squares Mean Square F Value Pr > F X X X
19 3 FACTOR CCD -- KHURI AND CORNELL (PAGE 9) Parameter DF Estimate Standard Error t Value Pr > t Parameter Estimate from Coded Data Intercept < X X X X*X X2*X X2*X X3*X X3*X X3*X FACTOR CCD -- KHURI AND CORNELL (PAGE 9) Canonical Analysis Sum of Response of Surface Based on Coded Data Factor DF Squares Mean Square F Value Pr > F X Critical Value X2 4 Factor Coded Uncoded X X X X Predicted value at stationary point: Eigenvectors Eigenvalues X X2 X Stationary point is a saddle point. 32
20 3 FACTOR CCD -- KHURI AND CORNELL (PAGE 9) Fit Diagnostics for YIELD RStudent - RStudent Predicted Value Predicted Value Leverage YIELD 2 8 Cook's D Quantile 8 2 Predicted Value Observation 3 2 Fit Mean Percent 2-2 Observations Parameters Error DF MSE R-Square Adj R-Square Proportion Less 33
21 3 FACTOR CCD -- KHURI AND CORNELL (PAGE 9) 3 FACTOR CCD -- KHURI AND CORNELL (PAGE 9) 34
7.3 Ridge Analysis of the Response Surface
7.3 Ridge Analysis of the Response Surface When analyzing a fitted response surface, the researcher may find that the stationary point is outside of the experimental design region, but the researcher wants
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 information5.3 Three-Stage Nested Design Example
5.3 Three-Stage Nested Design Example A researcher designs an experiment to study the of a metal alloy. A three-stage nested design was conducted that included Two alloy chemistry compositions. Three ovens
More informationPerforming response surface analysis using the SAS RSREG procedure
Paper DV02-2012 Performing response surface analysis using the SAS RSREG procedure Zhiwu Li, National Database Nursing Quality Indicator and the Department of Biostatistics, University of Kansas Medical
More 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 information6 Designs with Split Plots
6 Designs with Split Plots Many factorial experimental designs are incorrectly analyzed because the assumption of complete randomization is not true. Many factorial experiments have one or more restrictions
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 information4.8 Alternate Analysis as a Oneway ANOVA
4.8 Alternate Analysis as a Oneway ANOVA Suppose we have data from a two-factor factorial design. The following method can be used to perform a multiple comparison test to compare treatment means as well
More informationSTATISTICS 479 Exam II (100 points)
Name STATISTICS 79 Exam II (1 points) 1. A SAS data set was created using the following input statement: Answer parts(a) to (e) below. input State $ City $ Pop199 Income Housing Electric; (a) () Give the
More information1) Answer the following questions as true (T) or false (F) by circling the appropriate letter.
1) Answer the following questions as true (T) or false (F) by circling the appropriate letter. T F T F T F a) Variance estimates should always be positive, but covariance estimates can be either positive
More informationEXST7015: Estimating tree weights from other morphometric variables Raw data print
Simple Linear Regression SAS example Page 1 1 ********************************************; 2 *** Data from Freund & Wilson (1993) ***; 3 *** TABLE 8.24 : ESTIMATING TREE WEIGHTS ***; 4 ********************************************;
More informationJ. Response Surface Methodology
J. Response Surface Methodology Response Surface Methodology. Chemical Example (Box, Hunter & Hunter) Objective: Find settings of R, Reaction Time and T, Temperature that produced maximum yield subject
More informationLecture notes on Regression & SAS example demonstration
Regression & Correlation (p. 215) When two variables are measured on a single experimental unit, the resulting data are called bivariate data. You can describe each variable individually, and you can also
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 informationCHAPTER II RESPONSE SURFACE METHODOLOGY
CHAPTER II RESPOSE SURFACE METHODOLOGY.1. Introduction Response surface methodology is useful for factorial eperiments with v quantitative factors that are under taken so as to determine the combination
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 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 informationCOMPREHENSIVE WRITTEN EXAMINATION, PAPER III FRIDAY AUGUST 26, 2005, 9:00 A.M. 1:00 P.M. STATISTICS 174 QUESTION
COMPREHENSIVE WRITTEN EXAMINATION, PAPER III FRIDAY AUGUST 26, 2005, 9:00 A.M. 1:00 P.M. STATISTICS 174 QUESTION Answer all parts. Closed book, calculators allowed. It is important to show all working,
More informationEXST Regression Techniques Page 1. We can also test the hypothesis H :" œ 0 versus H :"
EXST704 - Regression Techniques Page 1 Using F tests instead of t-tests We can also test the hypothesis H :" œ 0 versus H :" Á 0 with an F test.! " " " F œ MSRegression MSError This test is mathematically
More informationLINEAR REGRESSION. Copyright 2013, SAS Institute Inc. All rights reserved.
LINEAR REGRESSION LINEAR REGRESSION REGRESSION AND OTHER MODELS Type of Response Type of Predictors Categorical Continuous Continuous and Categorical Continuous Analysis of Variance (ANOVA) Ordinary Least
More informationPhysics 403. Segev BenZvi. Parameter Estimation, Correlations, and Error Bars. Department of Physics and Astronomy University of Rochester
Physics 403 Parameter Estimation, Correlations, and Error Bars Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Review of Last Class Best Estimates and Reliability
More informationSAS code for Levene s Test and Weighted Least Squares DM LOG; CLEAR; OUT; CLEAR; ; ODS GRAPHICS ON; OPTIONS NODATE NONUMBER;
4.5 4.9 Weighted WeightedLeast Least Squares Squares and and Transformation Transformation Example Example 97 8 SAS code for Levene s Test and Weighted Least Squares M LOG; LEAR; OUT; LEAR; ; OS GRAPHIS
More informationCircle a single answer for each multiple choice question. Your choice should be made clearly.
TEST #1 STA 4853 March 4, 215 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 31 questions. Circle
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 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 informationA A x i x j i j (i, j) (j, i) Let. Compute the value of for and
7.2 - Quadratic Forms quadratic form on is a function defined on whose value at a vector in can be computed by an expression of the form, where is an symmetric matrix. The matrix R n Q R n x R n Q(x) =
More informationOutline. Review regression diagnostics Remedial measures Weighted regression Ridge regression Robust regression Bootstrapping
Topic 19: Remedies Outline Review regression diagnostics Remedial measures Weighted regression Ridge regression Robust regression Bootstrapping Regression Diagnostics Summary Check normality of the residuals
More informationChapter 8 (More on Assumptions for the Simple Linear Regression)
EXST3201 Chapter 8b Geaghan Fall 2005: Page 1 Chapter 8 (More on Assumptions for the Simple Linear Regression) Your textbook considers the following assumptions: Linearity This is not something I usually
More informationCh 2: Simple Linear Regression
Ch 2: Simple Linear Regression 1. Simple Linear Regression Model A simple regression model with a single regressor x is y = β 0 + β 1 x + ɛ, where we assume that the error ɛ is independent random component
More informationMulticollinearity Exercise
Multicollinearity Exercise Use the attached SAS output to answer the questions. [OPTIONAL: Copy the SAS program below into the SAS editor window and run it.] You do not need to submit any output, so there
More informationVersion 1: Equality of Distributions. 3. F (x) and G(x) represent the distribution functions corresponding to the Xs and Y s, respectively.
4 Two-Sample Methods 4.1 The (Mann-Whitney) Wilcoxon Rank Sum Test Version 1: Equality of Distributions Assumptions: Given two independent random samples X 1, X 2,..., X n and Y 1, Y 2,..., Y m : 1. The
More information6 Single Sample Methods for a Location Parameter
6 Single Sample Methods for a Location Parameter If there are serious departures from parametric test assumptions (e.g., normality or symmetry), nonparametric tests on a measure of central tendency (usually
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 informationSimple Linear Regression
Simple Linear Regression September 24, 2008 Reading HH 8, GIll 4 Simple Linear Regression p.1/20 Problem Data: Observe pairs (Y i,x i ),i = 1,...n Response or dependent variable Y Predictor or independent
More informationAN ITERATION. In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b
AN ITERATION In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b In this, A is an n n matrix and b R n.systemsof this form arise
More informationAlgebra 1 Standards Curriculum Map Bourbon County Schools. Days Unit/Topic Standards Activities Learning Targets ( I Can Statements) 1-19 Unit 1
Algebra 1 Standards Curriculum Map Bourbon County Schools Level: Grade and/or Course: Updated: e.g. = Example only Days Unit/Topic Standards Activities Learning Targets ( I 1-19 Unit 1 A.SSE.1 Interpret
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More information17. Example SAS Commands for Analysis of a Classic Split-Plot Experiment 17. 1
17 Example SAS Commands for Analysis of a Classic SplitPlot Experiment 17 1 DELIMITED options nocenter nonumber nodate ls80; Format SCREEN OUTPUT proc import datafile"c:\data\simulatedsplitplotdatatxt"
More informationChapter 8 Quantitative and Qualitative Predictors
STAT 525 FALL 2017 Chapter 8 Quantitative and Qualitative Predictors Professor Dabao Zhang Polynomial Regression Multiple regression using X 2 i, X3 i, etc as additional predictors Generates quadratic,
More informationLecture 2 - Unconstrained Optimization Definition[Global Minimum and Maximum]Let f : S R be defined on a set S R n. Then
Lecture 2 - Unconstrained Optimization Definition[Global Minimum and Maximum]Let f : S R be defined on a set S R n. Then 1. x S is a global minimum point of f over S if f (x) f (x ) for any x S. 2. x S
More informationSOLUTION OF QUADRATIC EQUATIONS LESSON PLAN. A3 Topic Overview ALGEBRA
ALGEBRA A Topic Overview A SOLUTION OF QUADRATIC EQUATIONS This topic describes three methods of solving Quadratic equations. assumes you understand and have practised using the algebraic methods described
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 informationThe use of response surface analysis in obtaining maximum profit in oil palm industry
ORIGINAL ARTICLE The use of response surface analysis in obtaining maximum profit in oil palm industry Azme Khamis 1, Zuhaimy Ismail 2, Khalid Haron 3 and Ahmad Tarmizi Mohammed 4 Abstract Khamis, A.,
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More informationThis is a Randomized Block Design (RBD) with a single factor treatment arrangement (2 levels) which are fixed.
EXST3201 Chapter 13c Geaghan Fall 2005: Page 1 Linear Models Y ij = µ + βi + τ j + βτij + εijk This is a Randomized Block Design (RBD) with a single factor treatment arrangement (2 levels) which are fixed.
More informationFailure Time of System due to the Hot Electron Effect
of System due to the Hot Electron Effect 1 * exresist; 2 option ls=120 ps=75 nocenter nodate; 3 title of System due to the Hot Electron Effect ; 4 * TIME = failure time (hours) of a system due to drift
More informationCircle the single best answer for each multiple choice question. Your choice should be made clearly.
TEST #1 STA 4853 March 6, 2017 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 32 multiple choice
More informationMATH 5720: Unconstrained Optimization Hung Phan, UMass Lowell September 13, 2018
MATH 57: Unconstrained Optimization Hung Phan, UMass Lowell September 13, 18 1 Global and Local Optima Let a function f : S R be defined on a set S R n Definition 1 (minimizers and maximizers) (i) x S
More informationStatistics for exp. medical researchers Regression and Correlation
Faculty of Health Sciences Regression analysis Statistics for exp. medical researchers Regression and Correlation Lene Theil Skovgaard Sept. 28, 2015 Linear regression, Estimation and Testing Confidence
More information= (Type exponential notation with positive exponents)
1. Subtract. Simplify by collecting like radical terms if possible. 2 2 = (Simplify your answer) 2. Add. Simplify if possible. = (Simplify your answer) 3. Divide and simplify. = (Type exponential notation
More informationQuadratic Formula: - another method for solving quadratic equations (ax 2 + bx + c = 0)
In the previous lesson we showed how to solve quadratic equations that were not factorable and were not perfect squares by making perfect square trinomials using a process called completing the square.
More informationAssignment 9 Answer Keys
Assignment 9 Answer Keys Problem 1 (a) First, the respective means for the 8 level combinations are listed in the following table A B C Mean 26.00 + 34.67 + 39.67 + + 49.33 + 42.33 + + 37.67 + + 54.67
More information1. The dosage in milligrams D of a heartworm preventive for a dog who weighs X pounds is given by D x. Substitute 28 in place of x to get:
1. The dosage in milligrams D of a heartworm preventive for a dog who weighs X pounds is given by D x 28 pounds. ( ) = 136 ( ). Find the proper dosage for a dog that weighs 25 x Substitute 28 in place
More informationSTOR 455 STATISTICAL METHODS I
STOR 455 STATISTICAL METHODS I Jan Hannig Mul9variate Regression Y=X β + ε X is a regression matrix, β is a vector of parameters and ε are independent N(0,σ) Es9mated parameters b=(x X) - 1 X Y Predicted
More information14 Multiple Linear Regression
B.Sc./Cert./M.Sc. Qualif. - Statistics: Theory and Practice 14 Multiple Linear Regression 14.1 The multiple linear regression model In simple linear regression, the response variable y is expressed in
More informationExst7037 Multivariate Analysis Cancorr interpretation Page 1
Exst7037 Multivariate Analysis Cancorr interpretation Page 1 1 *** C03S3D1 ***; 2 ****************************************************************************; 3 *** The data set insulin shows data from
More informationLecture 5: Comparing Treatment Means Montgomery: Section 3-5
Lecture 5: Comparing Treatment Means Montgomery: Section 3-5 Page 1 Linear Combination of Means ANOVA: y ij = µ + τ i + ɛ ij = µ i + ɛ ij Linear combination: L = c 1 µ 1 + c 1 µ 2 +...+ c a µ a = a i=1
More informationIES 612/STA 4-573/STA Winter 2008 Week 1--IES 612-STA STA doc
IES 612/STA 4-573/STA 4-576 Winter 2008 Week 1--IES 612-STA 4-573-STA 4-576.doc Review Notes: [OL] = Ott & Longnecker Statistical Methods and Data Analysis, 5 th edition. [Handouts based on notes prepared
More information7.4 The c-chart (fixed sample size)
7.4 The c-chart (fixed sample size) Often there is greater concern for the total number of nonconformities (or defects) (C) than the fraction of nonconforming (defective) units (p). An inspection unit
More informationICS 6N Computational Linear Algebra Symmetric Matrices and Orthogonal Diagonalization
ICS 6N Computational Linear Algebra Symmetric Matrices and Orthogonal Diagonalization Xiaohui Xie University of California, Irvine xhx@uci.edu Xiaohui Xie (UCI) ICS 6N 1 / 21 Symmetric matrices An n n
More informationUNIT 5 QUADRATIC FUNCTIONS Lesson 1: Interpreting Structure in Expressions Instruction
Prerequisite Skills This lesson requires the use of the following skills: evaluating expressions using the order of operations evaluating expressions for a given value identifying parts of an expression
More informationRidge Regression. Chapter 335. Introduction. Multicollinearity. Effects of Multicollinearity. Sources of Multicollinearity
Chapter 335 Introduction is a technique for analyzing multiple regression data that suffer from multicollinearity. When multicollinearity occurs, least squares estimates are unbiased, but their variances
More informationCh 3: Multiple Linear Regression
Ch 3: Multiple Linear Regression 1. Multiple Linear Regression Model Multiple regression model has more than one regressor. For example, we have one response variable and two regressor variables: 1. delivery
More informationAdaptive Beamforming Algorithms
S. R. Zinka srinivasa_zinka@daiict.ac.in October 29, 2014 Outline 1 Least Mean Squares 2 Sample Matrix Inversion 3 Recursive Least Squares 4 Accelerated Gradient Approach 5 Conjugate Gradient Method Outline
More informationAlgebra III: Blizzard Bag #1 Exponential and Logarithm Functions
NAME : DATE: PERIOD: Algebra III: Blizzard Bag #1 Exponential and Logarithm Functions Students need to complete the following assignment, which will aid in review for the end of course exam. Look back
More informationIntroduction to Linear Regression Rebecca C. Steorts September 15, 2015
Introduction to Linear Regression Rebecca C. Steorts September 15, 2015 Today (Re-)Introduction to linear models and the model space What is linear regression Basic properties of linear regression Using
More informationSpring Frost Prediction Models in Cranberry
University of Massachusetts Amherst ScholarWorks@UMass Amherst North American Cranberry Researcher and Extension Workers Conference NACREW 2017 Aug 28th, 9:30 AM - 9:45 AM Spring Frost Prediction Models
More informationThe integration of response surface method in microsoft excel with visual basic application
Journal of Physics: Conference Series PAPER OPEN ACCESS The integration of response surface method in microsoft excel with visual basic application To cite this article: H Sofyan et al 2018 J. Phys.: Conf.
More informationMathematics High School Algebra I
Mathematics High School Algebra I All West Virginia teachers are responsible for classroom instruction that integrates content standards and mathematical habits of mind. Students in this course will focus
More informationAssignment #1 MAT121 Summer 2015 NAME:
Assignment #1 MAT11 Summer 015 NAME: Directions: Do ALL of your work on THIS handout in the space provided! Circle your final answer! On problems that your teacher would show work on be sure that you also
More informationREGRESSION DIAGNOSTICS AND REMEDIAL MEASURES
REGRESSION DIAGNOSTICS AND REMEDIAL MEASURES Lalmohan Bhar I.A.S.R.I., Library Avenue, Pusa, New Delhi 110 01 lmbhar@iasri.res.in 1. Introduction Regression analysis is a statistical methodology that utilizes
More informationPerformance Surfaces and Optimum Points
CSC 302 1.5 Neural Networks Performance Surfaces and Optimum Points 1 Entrance Performance learning is another important class of learning law. Network parameters are adjusted to optimize the performance
More informationChapter 11 : State SAT scores for 1982 Data Listing
EXST3201 Chapter 12a Geaghan Fall 2005: Page 1 Chapter 12 : Variable selection An example: State SAT scores In 1982 there was concern for scores of the Scholastic Aptitude Test (SAT) scores that varied
More informationdm'log;clear;output;clear'; options ps=512 ls=99 nocenter nodate nonumber nolabel FORMCHAR=" = -/\<>*"; ODS LISTING;
dm'log;clear;output;clear'; options ps=512 ls=99 nocenter nodate nonumber nolabel FORMCHAR=" ---- + ---+= -/\*"; ODS LISTING; *** Table 23.2 ********************************************; *** Moore, David
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 informationSimultaneous Optimization of Incomplete Multi-Response Experiments
Open Journal of Statistics, 05, 5, 430-444 Published Online August 05 in SciRes. http://www.scirp.org/journal/ojs http://dx.doi.org/0.436/ojs.05.55045 Simultaneous Optimization of Incomplete Multi-Response
More informationResponse-surface illustration
Response-surface illustration Russ Lenth October 23, 2017 Abstract In this vignette, we give an illustration, using simulated data, of a sequential-experimentation process to optimize a response surface.
More informationTopic 14: Inference in Multiple Regression
Topic 14: Inference in Multiple Regression Outline Review multiple linear regression Inference of regression coefficients Application to book example Inference of mean Application to book example Inference
More informationPY 351 Modern Physics Short assignment 4, Nov. 9, 2018, to be returned in class on Nov. 15.
PY 351 Modern Physics Short assignment 4, Nov. 9, 2018, to be returned in class on Nov. 15. You may write your answers on this sheet or on a separate paper. Remember to write your name on top. Please note:
More informationTopic 28: Unequal Replication in Two-Way ANOVA
Topic 28: Unequal Replication in Two-Way ANOVA Outline Two-way ANOVA with unequal numbers of observations in the cells Data and model Regression approach Parameter estimates Previous analyses with constant
More informationFunctions of Several Variables
Functions of Several Variables The Unconstrained Minimization Problem where In n dimensions the unconstrained problem is stated as f() x variables. minimize f()x x, is a scalar objective function of vector
More informationBE640 Intermediate Biostatistics 2. Regression and Correlation. Simple Linear Regression Software: SAS. Emergency Calls to the New York Auto Club
BE640 Intermediate Biostatistics 2. Regression and Correlation Simple Linear Regression Software: SAS Emergency Calls to the New York Auto Club Source: Chatterjee, S; Handcock MS and Simonoff JS A Casebook
More informationA quadratic expression is a mathematical expression that can be written in the form 2
118 CHAPTER Algebra.6 FACTORING AND THE QUADRATIC EQUATION Textbook Reference Section 5. CLAST OBJECTIVES Factor a quadratic expression Find the roots of a quadratic equation A quadratic expression is
More informationOdor attraction CRD Page 1
Odor attraction CRD Page 1 dm'log;clear;output;clear'; options ps=512 ls=99 nocenter nodate nonumber nolabel FORMCHAR=" ---- + ---+= -/\*"; ODS LISTING; *** Table 23.2 ********************************************;
More informationHONORS GEOMETRY Summer Skills Set
HONORS GEOMETRY Summer Skills Set Algebra Concepts Adding and Subtracting Rational Numbers To add or subtract fractions with the same denominator, add or subtract the numerators and write the sum or difference
More informationLinear Regression Linear Regression with Shrinkage
Linear Regression Linear Regression ith Shrinkage Introduction Regression means predicting a continuous (usually scalar) output y from a vector of continuous inputs (features) x. Example: Predicting vehicle
More informationStatistics 512: Solution to Homework#11. Problems 1-3 refer to the soybean sausage dataset of Problem 20.8 (ch21pr08.dat).
Statistics 512: Solution to Homework#11 Problems 1-3 refer to the soybean sausage dataset of Problem 20.8 (ch21pr08.dat). 1. Perform the two-way ANOVA without interaction for this model. Use the results
More informationProblems. Suppose both models are fitted to the same data. Show that SS Res, A SS Res, B
Simple Linear Regression 35 Problems 1 Consider a set of data (x i, y i ), i =1, 2,,n, and the following two regression models: y i = β 0 + β 1 x i + ε, (i =1, 2,,n), Model A y i = γ 0 + γ 1 x i + γ 2
More informationMaximum variance formulation
12.1. Principal Component Analysis 561 Figure 12.2 Principal component analysis seeks a space of lower dimensionality, known as the principal subspace and denoted by the magenta line, such that the orthogonal
More informationSIMPLE LINEAR REGRESSION
SIMPLE LINEAR REGRESSION In linear regreion, we conider the frequency ditribution of one variable (Y) at each of everal level of a econd variable (). Y i known a the dependent variable. The variable for
More informationBeginning Algebra. 1. Review of Pre-Algebra 1.1 Review of Integers 1.2 Review of Fractions
1. Review of Pre-Algebra 1.1 Review of Integers 1.2 Review of Fractions Beginning Algebra 1.3 Review of Decimal Numbers and Square Roots 1.4 Review of Percents 1.5 Real Number System 1.6 Translations:
More informationEXST 7015 Fall 2014 Lab 08: Polynomial Regression
EXST 7015 Fall 2014 Lab 08: Polynomial Regression OBJECTIVES Polynomial regression is a statistical modeling technique to fit the curvilinear data that either shows a maximum or a minimum in the curve,
More informationPRINCIPAL COMPONENTS ANALYSIS (PCA)
PRINCIPAL COMPONENTS ANALYSIS (PCA) Introduction PCA is considered an exploratory technique that can be used to gain a better understanding of the interrelationships between variables. PCA is performed
More informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More informationa. YOU MAY USE ONE 8.5 X11 TWO-SIDED CHEAT SHEET AND YOUR TEXTBOOK (OR COPY THEREOF).
STAT3503 Test 2 NOTE: a. YOU MAY USE ONE 8.5 X11 TWO-SIDED CHEAT SHEET AND YOUR TEXTBOOK (OR COPY THEREOF). b. YOU MAY USE ANY ELECTRONIC CALCULATOR. c. FOR FULL MARKS YOU MUST SHOW THE FORMULA YOU USE
More informationTopic 18: Model Selection and Diagnostics
Topic 18: Model Selection and Diagnostics Variable Selection We want to choose a best model that is a subset of the available explanatory variables Two separate problems 1. How many explanatory variables
More informationSolving Equations by Factoring. Solve the quadratic equation x 2 16 by factoring. We write the equation in standard form: x
11.1 E x a m p l e 1 714SECTION 11.1 OBJECTIVES 1. Solve quadratic equations by using the square root method 2. Solve quadratic equations by completing the square Here, we factor the quadratic member of
More informationTopic 13. Analysis of Covariance (ANCOVA) [ST&D chapter 17] 13.1 Introduction Review of regression concepts
Topic 13. Analysis of Covariance (ANCOVA) [ST&D chapter 17] 13.1 Introduction The analysis of covariance (ANCOVA) is a technique that is occasionally useful for improving the precision of an experiment.
More informationThe P/Q Mathematics Study Guide
The P/Q Mathematics Study Guide Copyright 007 by Lawrence Perez and Patrick Quigley All Rights Reserved Table of Contents Ch. Numerical Operations - Integers... - Fractions... - Proportion and Percent...
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 information