ENGI 4421 Propagation of Error Page 8-01

Similar documents
ENGI 4421 Joint Probability Distributions Page Joint Probability Distributions [Navidi sections 2.5 and 2.6; Devore sections

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

Econometric Methods. Review of Estimation

ENGI 3423 Simple Linear Regression Page 12-01

STATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS. x, where. = y - ˆ " 1

Functions of Random Variables

Simple Linear Regression

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Multiple Choice Test. Chapter Adequacy of Models for Regression

ESS Line Fitting

Summary of the lecture in Biostatistics

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

Mean is only appropriate for interval or ratio scales, not ordinal or nominal.

Objectives of Multiple Regression

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Maximum Likelihood Estimation

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE MODULE 5

Chapter Business Statistics: A First Course Fifth Edition. Learning Objectives. Correlation vs. Regression. In this chapter, you learn:

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Special Instructions / Useful Data

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Chapter 5 Properties of a Random Sample

Example: Multiple linear regression. Least squares regression. Repetition: Simple linear regression. Tron Anders Moger

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

STA302/1001-Fall 2008 Midterm Test October 21, 2008

Simple Linear Regression

STK4011 and STK9011 Autumn 2016

ε. Therefore, the estimate

Lecture 3. Sampling, sampling distributions, and parameter estimation

Third handout: On the Gini Index

Chapter 8. Inferences about More Than Two Population Central Values

Multivariate Transformation of Variables and Maximum Likelihood Estimation

X ε ) = 0, or equivalently, lim

b. There appears to be a positive relationship between X and Y; that is, as X increases, so does Y.

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy

Arithmetic Mean Suppose there is only a finite number N of items in the system of interest. Then the population arithmetic mean is

1. The weight of six Golden Retrievers is 66, 61, 70, 67, 92 and 66 pounds. The weight of six Labrador Retrievers is 54, 60, 72, 78, 84 and 67.

Statistics MINITAB - Lab 5

Sampling Theory MODULE V LECTURE - 14 RATIO AND PRODUCT METHODS OF ESTIMATION

1 Solution to Problem 6.40

Lecture 2: Linear Least Squares Regression

MEASURES OF DISPERSION

Chapter 8: Statistical Analysis of Simulated Data

MOLECULAR VIBRATIONS

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades

residual. (Note that usually in descriptions of regression analysis, upper-case

ECON 482 / WH Hong The Simple Regression Model 1. Definition of the Simple Regression Model

hp calculators HP 30S Statistics Averages and Standard Deviations Average and Standard Deviation Practice Finding Averages and Standard Deviations

Chapter 4 Multiple Random Variables

Class 13,14 June 17, 19, 2015

Investigation of Partially Conditional RP Model with Response Error. Ed Stanek

THE ROYAL STATISTICAL SOCIETY HIGHER CERTIFICATE

Chapter 13 Student Lecture Notes 13-1

( ) = ( ) ( ) Chapter 13 Asymptotic Theory and Stochastic Regressors. Stochastic regressors model

Lecture 1: Introduction to Regression

ρ < 1 be five real numbers. The

TESTS BASED ON MAXIMUM LIKELIHOOD

Lecture 8: Linear Regression

The expected value of a sum of random variables,, is the sum of the expected values:

Evaluation of uncertainty in measurements

Multiple Linear Regression Analysis

2.28 The Wall Street Journal is probably referring to the average number of cubes used per glass measured for some population that they have chosen.

2SLS Estimates ECON In this case, begin with the assumption that E[ i

THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA

Parameter, Statistic and Random Samples

STA 105-M BASIC STATISTICS (This is a multiple choice paper.)

Midterm Exam 1, section 2 (Solution) Thursday, February hour, 15 minutes

STA 108 Applied Linear Models: Regression Analysis Spring Solution for Homework #1

COV. Violation of constant variance of ε i s but they are still independent. The error term (ε) is said to be heteroscedastic.

THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA

CLASS NOTES. for. PBAF 528: Quantitative Methods II SPRING Instructor: Jean Swanson. Daniel J. Evans School of Public Affairs

Lecture Notes Types of economic variables

Module 7. Lecture 7: Statistical parameter estimation

Correlation and Simple Linear Regression

The equation is sometimes presented in form Y = a + b x. This is reasonable, but it s not the notation we use.

CHAPTER VI Statistical Analysis of Experimental Data

Wu-Hausman Test: But if X and ε are independent, βˆ. ECON 324 Page 1

Chapter 2 Supplemental Text Material

Point Estimation: definition of estimators

X X X E[ ] E X E X. is the ()m n where the ( i,)th. j element is the mean of the ( i,)th., then

Chapter Statistics Background of Regression Analysis

Chapter 14 Logistic Regression Models

Lecture 1 Review of Fundamental Statistical Concepts

Lecture 07: Poles and Zeros

Midterm Exam 1, section 1 (Solution) Thursday, February hour, 15 minutes

6.867 Machine Learning

THE ROYAL STATISTICAL SOCIETY 2010 EXAMINATIONS SOLUTIONS GRADUATE DIPLOMA MODULE 2 STATISTICAL INFERENCE

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

1. A real number x is represented approximately by , and we are told that the relative error is 0.1 %. What is x? Note: There are two answers.

Lecture 9: Tolerant Testing

Chapter 11 Systematic Sampling

Laboratory I.10 It All Adds Up

THE ROYAL STATISTICAL SOCIETY 2016 EXAMINATIONS SOLUTIONS GRADUATE DIPLOMA MODULE 2

9 U-STATISTICS. Eh =(m!) 1 Eh(X (1),..., X (m ) ) i.i.d

Chapter 13, Part A Analysis of Variance and Experimental Design. Introduction to Analysis of Variance. Introduction to Analysis of Variance

Chapter 3 Sampling For Proportions and Percentages

Measures of Dispersion

Lecture Note to Rice Chapter 8

Lecture 1: Introduction to Regression

Using Statistics To Make Inferences 9

Transcription:

ENGI 441 Propagato of Error Page 8-01 Propagato of Error [Navd Chapter 3; ot Devore] Ay realstc measuremet procedure cotas error. Ay calculatos based o that measuremet wll therefore also cota a error. The errors are propagated from measuremets to the calculated value. From Chapter 1, for ay costats a, b ay rom quatty X we kow that [ ] E[ ] E ax + b = a X + b V[ ax + b] = a V[ X] = a X From Chapter 7, for ay costats { c1, c,, c } { X1, X,, X }, we kow that ax+ b ay set of depedet rom quattes [ ] E cx 1 1 + cx + + cx = c1e X 1 + ce X + + ce X V cx 1 1 + cx + + cx = c1 V X 1 + c V[ X] + + c V X A rom sample of observatos draw from a populato of mea µ varace a sample mea X for whch has E X = µ X V X = = X Example 8.01 The dameter of a crcle s measured to be 6.0 cm, wth the ucertaty represeted by the stard devato of 0.1 cm. Fd the ucertaty the crcumferece of the crcle. Let D = measured dameter of the crcle C = resultg calculated crcumferece of the crcle the

ENGI 441 Propagato of Error Page 8-0 Example 8.0 A wegh scale s used to estmate masses of objects up to 100 kg wth a ucertaty of 0.5 kg. The error each measuremet s depedet of the errors all other measuremets. Four compoets of a mache are measured masses of (8.7 ± 0.5) kg, (10. ± 0.5) kg, (1.7 ± 0.5) kg (15.1 ± 0.5) kg are reported. What s the total mass of the four compoets? Let M = total mass M = mass of the th compoet, the Note o quoted ucertaty: The ucertaty b the expresso (a ± b) s, by default, the approprate stard devato. However, b s sometmes take to be the maxmum possble error. O other occasos t s some multple of the stard devato (ofte or 3). I chapters 11 1 we wll express cofdece tervals the form (a ± b). Be sure whch of these terpretatos s teded whe you ecouter (a ± b).

ENGI 441 Propagato of Error Page 8-03 Example 8.03 The mass of a compoet s measured sx tmes o a wegh scale whose ucertaty s ukow. The sx measuremets ( kg) are 5.09, 5.16, 5.08, 5.10, 5.14 5.1. Estmate the mass of the compoet fd the ucertaty ths estmate. Let M = sample mea of the measured masses S = sample stard devato. The true ucertaty s ukow, but t may be estmated by the observed s. m = 30.69 6 = 5.115 Repeated Measuremets wth Dfferet certates Example 8.04 Two depedet estmates of the durato of a chemcal process, usg dfferet clocks, are (4.0 ± 0.) s (4.1 ± 0.1) s. Fd the ucertaty the average tme of 4.05 s. Let T = average of the tmes o the two clocks. 4.0 + 4.1 t = = 4.05

ENGI 441 Propagato of Error Page 8-04 Example 8.04 (cotued) However, oe clock s more precse tha the other. It makes sese that greater weght should be placed o the value provded by the more precse clock. Let us use a weghted average T W place of the smple average T. W 1 V V 1 V T ct c T T c T c T 1 W 1 where c s a weght the terval [0, 1]. We eed the value of c that mmses ( therefore geerates the best possble precso). I geeral, whe two depedet estmates for the same value are reported as x x, the the most precse average value wll be x, where 1 1 x c x c x c c c c 1 1 1, 1 1, 1 1 1

ENGI 441 Propagato of Error Page 8-05 Lear Combato of Depedet Measuremets The geeral expresso for the varace of a lear combato of rom quattes Y a X s 1 VY V a X a a Cov X, X Ths ca be re-wrtte as j j 1 1 j 1 V Y a V X a a Cov X, X where j j 1 1 j 1 a aa jj j 1 1 j 1 j s the correlato coeffcet betwee X X j V X The maxmum possble value of V[Y] occurs whe sg a a, j We the obta V. j j. Y a a a a Y 1 1 1 I all other cases, a coservatve estmate of the ucertaty Y a a a Y 1 1 a X s 1 Example 8.05 Suppose that we are ot sure whether or ot the two estmates, usg dfferet clocks, of the durato of a chemcal process Example 8.04 are depedet. The estmates are (4.0 ± 0.) s (4.1 ± 0.1) s. Fd a coservatve overall estmate. T T 1 4.0 4.1 T T1 T t 1 As before, 4.05

ENGI 441 Propagato of Error Page 8-06 No-Lear Fuctos of Oe Measuremet Suppose that (X) s a o-lear fucto of a rom quatty X whose true mea s μ. If X s a ubased estmator of μ, the, geeral, (X) wll be a based estmator of ( μ ). If X s close to μ, the the Taylor seres approxmato for (X), d 1 d ( X) = ( μ) + ( X μ) + ( X μ ) + dx! dx X = μ X = μ may be trucated to just the frst order, d d d ( X) ( μ) + ( X μ) = ( μ) μ + X dx X = μ dx X = μ dx X = μ Ths s a lear fucto of X, whose stard devato follows from the geeral formula ax+ b = a X : d X dx X = μ Ths s the formula for propagato of error. I practce we do ot kow the true value of μ, so we evaluate the dervatve at the observed value of X stead. Example 8.06 I Example 8.01, the dameter D of the crcle was measured to be (6.0 ± 0.1) cm. Estmate the area A of the crcle.

ENGI 441 Propagato of Error Page 8-07 Relatve certaty [for bous questos oly] If s a ubased measuremet whose true value s µ whose absolute ucertaty s, the the relatve ucertaty s the dmesoless quatty µ (also kow as the coeffcet of varato). I practce, the relatve ucertaty s calculated as d Note that l ( l ) =. d Therefore the absolute ucertaty (l ) s also the relatve ucertaty. Example 8.06 (aga) [for bous questos oly] Fd the relatve ucertaty the area of the crcle of dameter (6.0 ± 0.1) cm hece fd the absolute ucertaty the area of the crcle D A= π = π D l A= l π D = l π + l D = l π + l D ( ) ( ) 4 4 4 4 d 1 ( l A) 0 dd = + D = 6 = 3 A 1 1 1 D = 0.1 = = 3.3% A 3 3 30 1 1 6 9π 3π A A = π = = 0.94 cm 30 30 30 10

ENGI 441 Propagato of Error Page 8-08 Example 8.07 (Navd textbook, exercses 3.3, page 185, questo 8) [examable] 1 The refractve dex η of a pece of glass s related to the crtcal agle θ by η =. sθ The crtcal agle s reported to be (0.70 ± 0.0) rad. Estmate the refractve dex fd the ucertaty the estmate. No-Lear Fuctos of Several Measuremets [for bous questos oly] If { X1, X,, X } are depedet measuremets whose ucertates { 1,,, } small f = ( X1, X,, X ) s a o-lear fucto of { X1, X,, X }, the are + + + X X X 1 1 where the partal dervatves are evaluated at the observed value of ( X1, X,, X ) Ths s the multvarate propagato of error formula..

ENGI 441 Propagato of Error Page 8-09 If the { X1, X,, X } are ot depedet, the a coservatve estmate s + + + X X X 1 1 Example 8.08 (Navd textbook, exercses 3.4, page 193, questo 8(b)) [for bous questos oly] The pressure P ( kpa), temperature T ( K) volume V ( ltres) of oe mole of a deal gas are related by the equato PV = 8.31T Gve that P = (4.5 ± 0.03) kpa T = (90.11 ± 0.0) K, estmate V. V 8.31T 8.31 90.11 = = = 9.94068 P 4.5 1 V 1 8.31 V = 8.31TP = 8.31P = = 0.034 65 T 4.5 V 8.31 90.11 = 8.31TP = = 0.0409891 P ( 4.5) Assumg that the measuremets of pressure temperature are depedet, V V V T P ( ) ( ) 0.034 0.0 0.040 0.03 0.001408 T + P = + Therefore V = ( 9.9407 ± 0.0014 ) ltre If we caot assume depedece, the a coservatve estmate s gve by V V V T P 0.034 0.0 0.040 0.03 0.001915 T + P = + Therefore V = ( 9.9407 ± 0.0019 ) ltre [The aswer s gve to the ffth sgfcat fgure V, because P T were both quoted to that level of precso. However, the costat s quoted to oly 3 s.f., so the aswer may be urelable beyod the thrd sgfcat fgure.] [Ed of Chapter 8]

ENGI 441 Propagato of Error Page 8-10 [Space for Addtoal Notes]