COMMON ERRORS IN NUMERICAL PROBLEM SOLVING HOW CAN THEY BE DETECTED AND PREVENTED
|
|
- Barry Fisher
- 5 years ago
- Views:
Transcription
1 COMMON ERRORS IN NUMERICL PROBLEM SOLVING HOW CN THEY BE DETECTED ND PREVENTED Mordecha Shacham, Ben-Guron Unversty, Beer-Sheva, Israel Mchael B. Cutlp, Unversty of Connectcut, Storrs, CT 6269, US Mchael Elly, Intel Corporaton, Qryat Gat, Israel Introducton The ntroducton of mathematcal software packages such as Excel 1, MTLB 2 and Polymath 3 nto the feld of engneerng problem solvng has brought the benefts of hgher speed and precson n obtanng the results. However, numercal soluton technques have also ntroduced some new sources of errors. These errors may often pass undetected because of the lack of experence n the use of the new computatonal tools. few examples of errors orgnatng from numercal sources follow. 1. Regresson and nalyss of Data The most common errors n regresson of data orgnate from the use of regresson models wth too many or too few parameters [1], fttng polynomal models wthout proper scalng (standardzaton) of the temperature data, correlaton of data when the model equatons are mproperly lnearzed [1] or when expermental desgn for obtanng the data s not satsfactory [3]. 2. Systems of Nonlnear lgebrac Equatons (NLE) There are many examples n the lterature showng that dentfcaton of good ntal guesses and formulatng the problem correctly for enablng convergence to a soluton of an NLE system may represent a major challenge [4]. Even f the soluton s found, t may be nfeasble n the physcal sense (soluton wth negatve concentratons, for example) or a false soluton caused by mproper varable and/or functon scalng. 3. Ordnary Dfferental Equatons (ODE) Indscrmnate use of default error tolerances of the ODE solver tools s the most common source of errors n solvng ordnary dfferental equatons [5]. Falure to use the proper ntegraton algorthm (stff vs. non-stff), carelessly roundng numbers n the model equatons, usng the model outsde the doman of ts valdty and the use of nadequate resoluton n presentng the results have been documented [1],[2] as addtonal common sources of errors n solvng ODEs. Beng aware of potental causes to errors n numercal problem solvng s a prelmnary requrement to detecton and preventon of such errors. In ths paper three potental causes of errors are descrbed and demonstrated. Means to prevent these errors are proposed. No Consderaton Gven to Problem Parameter Uncertantes The mportance of the analyss of the senstvty of the problem soluton to small changes n the parameter values has been recognzed n the feld of optmzaton. Edgar et al. (21), for example, wrote "You should always be aware of the senstvty of the optmal answer Parameter values usually contan errors or uncertantes. Informaton concernng the senstvty of the optmum to changes or varatons n a parameter s therefore very mportant n optmal process desgn." However, senstvty analyss s rarely carred for other types of problems. In 1 Excel s a trademark of Mcrosoft Corporaton ( 2 MTLB s a trademark of The Math Works, Inc. ( 3 POLYMTH s a product of Polymath Software (
2 ths secton, the mportance of the senstvty analyss n solvng systems of nonlnear algebrac equatons wll be demonstrated. The example consdered was orgnally publshed by Froment and Bshoff (199), and ts soluton usng Mathcad was presented by Parulekar (26). The reacton consdered s the catalytc hydrogenaton of olefns (component ) n an sothermal CSTR. Dfferental mass balance on component yelds dc dt C = C τ r (1) where C s the nlet concentraton of (mol/l), τ = V / v s the space tme (s), V s the reactor volume (L), v s the nlet flow rate (L/s) and r s the reacton rate per unt reactor volume (mol/l s): r C = (2) ( 1 C ) 2 The concentraton of at steady-state s sought for the operatng condtons: C = 13 mol/l; V = 1 L; and v =.2 L/s. Fndng the steady-state soluton requres solvng the nonlnear algebrac equaton f(c ) = where f C C ( C ) = r (3) τ In Fgure 1, f(c ) s plotted versus C for the regon C 1. Observe that n ths regon there are three real solutons. The solutons as obtaned by the Polymath 6.1 software package are: C * = [ , , and ]. Polymath dsplays the results showng 7 sgnfcant dgts. Mathcad dsplays the same results wth 2 sgnfcant dgts (Parulecar, 26). The dscrepancy between the number of sgnfcant dgts n the CSTR parameters (one or two) and the number of those dgts reported n the soluton (up to twenty) cannot be mssed. The determnaton of the number of sgnfcant dgts to be shown n the documentaton of the problem soluton should be based, n ths partcular case, on the uncertanty assocated wth the parameter values. In Table 1, the values of C at the soluton are presented for varous levels of uncertanty n the flow rate v. ssumng an uncertanty of 1-4 L/s (thus flow rate of.21 L/s nstead of.2 L/s) yelds soluton values whch match the base soluton only n the frst two dgts. For uncertanty of 1-3 L/s there s only one correct dgt n the frst two roots and two correct dgts n the largest root. For an uncertanty of.1 L/s there s not even a sngle correct dgt n the frst two roots and one n the largest root. Thus, any realstc estmate for the uncertanty n the value of the flow rate leads to the concluson that the solutons for C should be rounded up to two sgnfcant dgts as all the addtonal dgts are uncertan. It should be emphaszed that ths analyss apples only to physcal quanttes whch obtaned as the fnal results of the computaton. There are other cases where roundng the numbers may lead to sgnfcant errors. Ths wll be demonstrated n one of the followng examples. s the value of the flow rate s specfed by one sgnfcant dgt, only one possble assumpton regardng the uncertanty could be that the flow rate was measured (or controlled)
3 wth such an accuracy, thus ts value could vary n the range of.15 L/s v.24 L/s. In the last two rows of Table 1 the solutons for the two extreme values of v are shown. For v =.24 L/s there s only one real root at C = mol/l whle for v =.15 L/s there s one real root at C =.3419 mol/l. Thus, n these cases there s only one root nstead of the three roots for the base value of v =.2 L/s, and the roots are completely dfferent. Use of Inapproprate Statstcs and Graphcal Representatons n nalyss of Regresson Models The use of napproprate statstcs and graphcal representaton may often lead to the acceptance unsatsfactory regresson models. Ths phenomenon wll be demonstrated by fttng the Clausus-Clapeyron equaton to vapor pressure (P V ) data of ethane. The vapor pressure data, contanng N = 17 ponts, are avalable [8] n the temperature range of T = 92 K (P V = 1.7 Pa) through T = 34 K (P V = 4.738*1 6 Pa). Ths range covers essentally the lqud phase regon between the trple pont temperature (= K) and the crtcal temperature (T C = K) of ethane. The Clausus-Clapeyron equaton s consdered napproprate for modelng vapor pressure for such a wde range of temperatures and we want to examne what knd of statstcs and/or graphcal representaton can reveal ths napproprateness. The Clausus-Clapeyron equaton s gven by B ln ( P V ) = + (4) T where and B are coeffcents calculated by regresson of expermental data. In Fgure 2 ln(p V ) exp s plotted versus 1/T. The straght lne ftted to the data s also shown. The lnear equaton obtaned s yˆ = x (5) where N 2 = 1 = N = 1 x = 1/ T and y ˆ = ln( P ) V calc. The correlaton coeffcent R2 s defned ( yˆ y) 2 R (6) ( ) 2 y y where y s the sample mean of the dependent varable. The value of R s bounded: R 1. If R s close to 1, there s a strong correlaton between the varables, whereas a value close to zero ndcates a weak or no correlaton. The R 2 value n ths case (as shown n Fgure 2) s.999, very close to 1. Thus, the plot of ln(p V ) exp versus 1/T and the statstcal R 2 both ndcate that the Clausus-Clapeyron equaton represents the vapor pressure data adequately over the entre range from the trple pont to T C. nother wdely used graphcal representaton for checkng adequacy of a correlaton s the plot of the calculated value of the dependent varable ŷ versus the expermental value y. In case of a good ft a straght lne should be obtaned wth the slope ~1.. In Fgure 3, ln(p V ) calc (as calculated by the Clausus-Clapeyron equaton) s plotted versus ln(p V ) exp. Observe that the data ponts are algned along a straght lne wth the slope of.9999 wth only small devatons. Thus, ths type of representaton ndcates also a good ft between the Clausus-Clapeyron equaton and the vapor pressure data over the entre range. However, comparng ln(p V ) calc and ln(p V ) exp at T = 92 K, for example, yelds ln(p V ) exp =.53; ln(p V ) calc =.95, whch s an ~ 8 %
4 dfference. Ths defntely cannot be defned as adequate representaton of the expermental data. The proper graphcal means to test the goodness of ft between a regresson model and the data s the "resdual" plot. The resdual εˆ s defned as ˆ ε = y yˆ (7) random dstrbuton of the resduals around zero ndcates that the regresson model represents the data correctly. defnte trend or pattern n the resdual plot may ndcate ether lack of ft of the model or that the assumed error dstrbuton for the data (.e. random error dstrbuton n y) s ncorrect. In Fgure 4 the resdual: ln(p V ) exp - ln(p V ) calc s plotted versus ln(p V ) exp. Observe that there s defnte trend (curvature) n the resduals, ndcatng that the Clausus-Clapeyron equaton cannot represent the data adequately. Careless Roundng of Model Parameters. In the DIPPR database [9], the Redel equaton s recommended for modelng vapor pressure data. The Redel equaton can be wrtten B E ln ( P V) = + + C lnt + DT (8) T where,b,c,d and E are the parameters of the regresson model. Usually the value of E s set at 2 to enable calculaton of the rest of the parameters by lnear regresson. Usng the same vapor pressure data that were dscussed n the prevous secton, the DIPPR staff obtaned the followng parameter values: = E+1; B = E+3; C = E+; and D = E-5. In Fgure 5, the resdual plot of ln(p V ) calculated by the Redel equaton usng the DIPPR parameter values s shown. Observe that the resduals n ths case are smaller by at least one order of magntude than the resduals obtaned usng the Clausus-Clapeyron equaton. There are two separate regons. In the frst one (up to ln(p V ) ~ 12) the resduals are n the range of -.4 to.2, not showng any partcular trend. In the second regon (ln(p V ) > 12) there s an order of magntude reducton n the resdual values (all smaller n absolute value than.3), however there s some trend n the resdual values n ths regon. To obtan perfectly random dstrbuton of the resduals, a dfferent correlaton should be used (the Wagner equaton, for example). Improvement of the correlaton by use of a dfferent type of equaton, however, s outsde of the scope of the present paper. The "copy and paste" faclty s a very convenent and tme-savng opton when transferrng numercal and other data between varous data-bases and software packages. For example, the Redel equaton constants for varous compounds can be coped nto an Excel spreadsheet. fter pastng the parameter values nto Excel, they are stored wth the exact same value as they appear n the DIPPR database, however they are dsplayed rounded to two three decmal dgts: = 5.19E+1; B = -2.6E+3; C = -5.13E+; and D = 1.49E-5 when the default number dsplay format s used. Copyng these values from Excel and pastng t nto a document or a dfferent software package wll results n a loss of the addtonal dgts whch were ncluded n the DIPPR publshed values. Calculatng the vapor pressure usng the rounded values of the parameters leads to sgnfcant errors n the calculated pressure values. The resdual plots of P V calculated by the Redel equaton usng exact and rounded parameter values are shown n Fgure 6. Observe that the resdual (error) n the hghest pressure range s of the order of 2.5% when rounded
5 parameter values are used n comparson to.3 % wth exact parameter values. Thus, careless roundng of parameter values causes the predcton error, n ths example, to ncrease by almost one order of magntude. Conclusons complete defnton of a problem should specfy the parameter values and the uncertanty levels for these parameters. Specfcaton of the uncertanty levels of the parameters enables a determnaton of the uncertanty of the soluton. Ths provdes justfcaton to number of sgnfcant dgts to be used n documentaton of the results and sometmes (as n the example presented) t may reveal that the proposed process cannot be mplemented n practce because of levels of uncertanty that are too hgh. The correlaton coeffcent ( R 2 ) and a smple plot of calculated values versus expermental values can be msleadng n determnng the goodness of the ft between the data and the regresson model. The resdual plot s the most relable means for examnng the goodness of the ft. transfer of numercal data between varous databases, documents and programs by electronc "copy and paste" may result n a loss of sgnfcant dgts of mportant parameter values. Beware of unntended roundng durng such transfer. Proper attenton to the potental causes of errors that are presented here and n prevous publcatons [1],[2],[3],[4],[5] can consderably mprove the relablty of the results n numercal problem solvng and n the correlaton of data by regresson. References 1. Brauner, N., M. Shacham and M. B. Cutlp, Computatonal Results: How Relable re They? - Systematc pproach to Modal Valdaton'', Chem. Engr. Educaton, 3(1), 2-25 (1996). 2. Shacham, M., N. Brauner and M. Pozn, "Potental Ptfalls n Usng General Purpose Software for Interactve Soluton of Ordnary Dfferental Equatons", cta Chmca Slovenca, 42(1), (1995). 3. Shacham, M. and N. Brauner, "Correlaton and Over-correlaton of Heterogeneous Reacton Rate Data", Chem. Engr. Educaton, 29(1) 22-25, 45 (1995). 4. Shacham, M., N. Brauner and M. B. Cutlp, Web-based Lbrary for Testng Performance of Numercal Software for Solvng Nonlnear lgebrac Equatons, Comput. Chem. Eng. 26(4-5), (22). 5. Shacham, M., N. Brauner, W. R. shurst, and M. B. Cutlp, "Can I Trust ths Software Package? n Exercse n Valdaton of Computatonal Results ", Chem. Engr. Educaton, 42(1), (28). 6. Edgar, T. F., Hmmelblau, D. M., and Lasdon, L. S., Optmzaton of Chemcal Processes, 2nd Ed., McGraw-Hll, New York, Froment, G.F., and Bshoff, K. B., Chemcal Reactor nalyss and Desgn, 2 nd Ed., Wley, New York, Parulekar, S. J., Numercal Problem Solvng Usng Mathcad n Undergraduate Reacton Engneerng, Chem. Engr. Educaton, 4 (1), 14 (26). 8. Ingham, H.; Frend, D.G.; Ely, J.F.; "Thermophyscal Propertes of Ethane"; J. Phys. Ref. Data, 2, 275(1991).
6 9. Rowley, R.L.; Wldng, W.V.; Oscarson, J.L.; Yang, Y.; Zundel, N.. DIPPR Data Complaton of Pure Chemcal Propertes Desgn Insttute for Physcal Propertes, (http//dppr.byu.edu), Brgham Young Unversty Provo Utah, 26. Table 1. Isothermal CSTR solutons for varous levels of uncertanty n the flow rate v. Feed flow rate (v ) L/s C, 1st Root C, 2nd Root C, 3rd Root Correct dgts are shown n bold f(c ) C Fgure 1. Plot of f(c ) versus C for v =.2 L/s 16 ln(pv )exp 12 8 y = x R 2 = /T (K -1 )
7 Fgure 2. Plot of ln(p V ) exp versus 1/T for ethane between the trple pont and T C ln(pv )Calc 8 y =.9999x ln(p V ) exp Fgure 3. Plot of ln(p V ) calc, calculated by the Clausus-Clapeyron equaton, versus ln(p V ) exp for ethane between the trple pont and T C.2.1 Resdual ln(p V ) exp Fgure 4. Resdual plot for ln(p V ) calculated by the Clausus - Clapeyron equaton for ethane between the trple pont T C Resdual ln(p V ) exp Fgure 5. Resdual plot for ln(p V ) calculated by the Redel equaton for ethane between the trple pont and T C
8 Exact Rounded 4.E+4.E+ Resdual -4.E+4-8.E+4-1.2E+5-1.6E+5.E+ 1.E+6 2.E+6 3.E+6 4.E+6 5.E+6 P Vexp (Pa) Fgure 6. Resdual plot for P V, calculated by the Redel equaton, for ethane between the trple pont and the T C
is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors
Multple Lnear and Polynomal Regresson wth Statstcal Analyss Gven a set of data of measured (or observed) values of a dependent varable: y versus n ndependent varables x 1, x, x n, multple lnear regresson
More informationChapter 11: Simple Linear Regression and Correlation
Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests
More informationComparison of Regression Lines
STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence
More informationConsidering Numerical Error Propagation in Modeling and Regression of Data
Sesson 2520 Consderng Numercal Error Propagaton n Modelng and Regresson of Data Nema Brauner, Mordecha Shacham Tel-Avv Unversty/Ben-Guron Unversty of the Negev Abstract The use of user-frendly nteractve
More informationNegative Binomial Regression
STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...
More informationStatistics for Economics & Business
Statstcs for Economcs & Busness Smple Lnear Regresson Learnng Objectves In ths chapter, you learn: How to use regresson analyss to predct the value of a dependent varable based on an ndependent varable
More informationStatistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation
Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 13 The Smple Lnear Regresson Model and Correlaton 1999 Prentce-Hall, Inc. Chap. 13-1 Chapter Topcs Types of Regresson Models Determnng the Smple Lnear
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationDepartment of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6
Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.
More informationThis column is a continuation of our previous column
Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard
More informationPolynomial Regression Models
LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationLecture 16 Statistical Analysis in Biomaterials Research (Part II)
3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan
More informationStatistics for Business and Economics
Statstcs for Busness and Economcs Chapter 11 Smple Regresson Copyrght 010 Pearson Educaton, Inc. Publshng as Prentce Hall Ch. 11-1 11.1 Overvew of Lnear Models n An equaton can be ft to show the best lnear
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationx = , so that calculated
Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to
More informationIntroduction to Regression
Introducton to Regresson Dr Tom Ilvento Department of Food and Resource Economcs Overvew The last part of the course wll focus on Regresson Analyss Ths s one of the more powerful statstcal technques Provdes
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More information10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization
10.34 Numercal Methods Appled to Chemcal Engneerng Fall 2015 Homework #3: Systems of Nonlnear Equatons and Optmzaton Problem 1 (30 ponts). A (homogeneous) azeotrope s a composton of a multcomponent mxture
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationStatistics MINITAB - Lab 2
Statstcs 20080 MINITAB - Lab 2 1. Smple Lnear Regresson In smple lnear regresson we attempt to model a lnear relatonshp between two varables wth a straght lne and make statstcal nferences concernng that
More informationBasic Business Statistics, 10/e
Chapter 13 13-1 Basc Busness Statstcs 11 th Edton Chapter 13 Smple Lnear Regresson Basc Busness Statstcs, 11e 009 Prentce-Hall, Inc. Chap 13-1 Learnng Objectves In ths chapter, you learn: How to use regresson
More informationPsychology 282 Lecture #24 Outline Regression Diagnostics: Outliers
Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationStatistics II Final Exam 26/6/18
Statstcs II Fnal Exam 26/6/18 Academc Year 2017/18 Solutons Exam duraton: 2 h 30 mn 1. (3 ponts) A town hall s conductng a study to determne the amount of leftover food produced by the restaurants n the
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
More information[The following data appear in Wooldridge Q2.3.] The table below contains the ACT score and college GPA for eight college students.
PPOL 59-3 Problem Set Exercses n Smple Regresson Due n class /8/7 In ths problem set, you are asked to compute varous statstcs by hand to gve you a better sense of the mechancs of the Pearson correlaton
More informationChapter 6. Supplemental Text Material
Chapter 6. Supplemental Text Materal S6-. actor Effect Estmates are Least Squares Estmates We have gven heurstc or ntutve explanatons of how the estmates of the factor effects are obtaned n the textboo.
More information2016 Wiley. Study Session 2: Ethical and Professional Standards Application
6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton
More informationChapter 15 - Multiple Regression
Chapter - Multple Regresson Chapter - Multple Regresson Multple Regresson Model The equaton that descrbes how the dependent varable y s related to the ndependent varables x, x,... x p and an error term
More informationCorrelation and Regression. Correlation 9.1. Correlation. Chapter 9
Chapter 9 Correlaton and Regresson 9. Correlaton Correlaton A correlaton s a relatonshp between two varables. The data can be represented b the ordered pars (, ) where s the ndependent (or eplanator) varable,
More informationAdiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram
Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationChapter 15 Student Lecture Notes 15-1
Chapter 15 Student Lecture Notes 15-1 Basc Busness Statstcs (9 th Edton) Chapter 15 Multple Regresson Model Buldng 004 Prentce-Hall, Inc. Chap 15-1 Chapter Topcs The Quadratc Regresson Model Usng Transformatons
More informationSTAT 3008 Applied Regression Analysis
STAT 3008 Appled Regresson Analyss Tutoral : Smple Lnear Regresson LAI Chun He Department of Statstcs, The Chnese Unversty of Hong Kong 1 Model Assumpton To quantfy the relatonshp between two factors,
More informationNumber of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k
ANOVA Model and Matrx Computatons Notaton The followng notaton s used throughout ths chapter unless otherwse stated: N F CN Y Z j w W Number of cases Number of factors Number of covarates Number of levels
More informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationLab 2e Thermal System Response and Effective Heat Transfer Coefficient
58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationDepartment of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution
Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable
More informationLecture 3 Stat102, Spring 2007
Lecture 3 Stat0, Sprng 007 Chapter 3. 3.: Introducton to regresson analyss Lnear regresson as a descrptve technque The least-squares equatons Chapter 3.3 Samplng dstrbuton of b 0, b. Contnued n net lecture
More informationA Robust Method for Calculating the Correlation Coefficient
A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal
More information1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands
Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of
More informationStatistics Chapter 4
Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment
More informationPHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University
PHYS 45 Sprng semester 7 Lecture : Dealng wth Expermental Uncertantes Ron Refenberger Brck anotechnology Center Purdue Unversty Lecture Introductory Comments Expermental errors (really expermental uncertantes)
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationDurban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications
Durban Watson for Testng the Lack-of-Ft of Polynomal Regresson Models wthout Replcatons Ruba A. Alyaf, Maha A. Omar, Abdullah A. Al-Shha ralyaf@ksu.edu.sa, maomar@ksu.edu.sa, aalshha@ksu.edu.sa Department
More informationThe Ordinary Least Squares (OLS) Estimator
The Ordnary Least Squares (OLS) Estmator 1 Regresson Analyss Regresson Analyss: a statstcal technque for nvestgatng and modelng the relatonshp between varables. Applcatons: Engneerng, the physcal and chemcal
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationMeasurement Uncertainties Reference
Measurement Uncertantes Reference Introducton We all ntutvely now that no epermental measurement can be perfect. It s possble to mae ths dea quanttatve. It can be stated ths way: the result of an ndvdual
More informationT E C O L O T E R E S E A R C H, I N C.
T E C O L O T E R E S E A R C H, I N C. B rdg n g En g neern g a nd Econo mcs S nce 1973 THE MINIMUM-UNBIASED-PERCENTAGE ERROR (MUPE) METHOD IN CER DEVELOPMENT Thrd Jont Annual ISPA/SCEA Internatonal Conference
More informationChapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.
Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where y + = β + β e for =,..., y and are observable varables e s a random error How can an estmaton rule be constructed for the
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationAnswers Problem Set 2 Chem 314A Williamsen Spring 2000
Answers Problem Set Chem 314A Wllamsen Sprng 000 1) Gve me the followng crtcal values from the statstcal tables. a) z-statstc,-sded test, 99.7% confdence lmt ±3 b) t-statstc (Case I), 1-sded test, 95%
More informationLecture 6: Introduction to Linear Regression
Lecture 6: Introducton to Lnear Regresson An Manchakul amancha@jhsph.edu 24 Aprl 27 Lnear regresson: man dea Lnear regresson can be used to study an outcome as a lnear functon of a predctor Example: 6
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationx i1 =1 for all i (the constant ).
Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationLearning Objectives for Chapter 11
Chapter : Lnear Regresson and Correlaton Methods Hldebrand, Ott and Gray Basc Statstcal Ideas for Managers Second Edton Learnng Objectves for Chapter Usng the scatterplot n regresson analyss Usng the method
More informationBOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu
BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS M. Krshna Reddy, B. Naveen Kumar and Y. Ramu Department of Statstcs, Osmana Unversty, Hyderabad -500 007, Inda. nanbyrozu@gmal.com, ramu0@gmal.com
More informationResource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis
Resource Allocaton and Decson Analss (ECON 800) Sprng 04 Foundatons of Regresson Analss Readng: Regresson Analss (ECON 800 Coursepak, Page 3) Defntons and Concepts: Regresson Analss statstcal technques
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationCopyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for P Charts. Dr. Wayne A. Taylor
Taylor Enterprses, Inc. Control Lmts for P Charts Copyrght 2017 by Taylor Enterprses, Inc., All Rghts Reserved. Control Lmts for P Charts Dr. Wayne A. Taylor Abstract: P charts are used for count data
More informationDevelopment of a Semi-Automated Approach for Regional Corrector Surface Modeling in GPS-Levelling
Development of a Sem-Automated Approach for Regonal Corrector Surface Modelng n GPS-Levellng G. Fotopoulos, C. Kotsaks, M.G. Sders, and N. El-Shemy Presented at the Annual Canadan Geophyscal Unon Meetng
More informationLecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding
Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 008 Recall: man dea of lnear regresson Lnear regresson can be used to study
More informationStatistics for Managers Using Microsoft Excel/SPSS Chapter 14 Multiple Regression Models
Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 14 Multple Regresson Models 1999 Prentce-Hall, Inc. Chap. 14-1 Chapter Topcs The Multple Regresson Model Contrbuton of Indvdual Independent Varables
More informationAn identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites
IOP Conference Seres: Materals Scence and Engneerng PAPER OPE ACCESS An dentfcaton algorthm of model knetc parameters of the nterfacal layer growth n fber compostes o cte ths artcle: V Zubov et al 216
More information18. SIMPLE LINEAR REGRESSION III
8. SIMPLE LINEAR REGRESSION III US Domestc Beers: Calores vs. % Alcohol Ftted Values and Resduals To each observed x, there corresponds a y-value on the ftted lne, y ˆ ˆ = α + x. The are called ftted values.
More informationCinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure
nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n
More informationLecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding
Recall: man dea of lnear regresson Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 8 Lnear regresson can be used to study an
More information28. SIMPLE LINEAR REGRESSION III
8. SIMPLE LINEAR REGRESSION III Ftted Values and Resduals US Domestc Beers: Calores vs. % Alcohol To each observed x, there corresponds a y-value on the ftted lne, y ˆ = βˆ + βˆ x. The are called ftted
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationIntroduction to Generalized Linear Models
INTRODUCTION TO STATISTICAL MODELLING TRINITY 00 Introducton to Generalzed Lnear Models I. Motvaton In ths lecture we extend the deas of lnear regresson to the more general dea of a generalzed lnear model
More informationOriginated from experimental optimization where measurements are very noisy Approximation can be actually more accurate than
Surrogate (approxmatons) Orgnated from expermental optmzaton where measurements are ver nos Approxmaton can be actuall more accurate than data! Great nterest now n applng these technques to computer smulatons
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationChapter 8 Indicator Variables
Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationSPANC -- SPlitpole ANalysis Code User Manual
Functonal Descrpton of Code SPANC -- SPltpole ANalyss Code User Manual Author: Dale Vsser Date: 14 January 00 Spanc s a code created by Dale Vsser for easer calbratons of poston spectra from magnetc spectrometer
More informationSIMPLE LINEAR REGRESSION
Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two
More informationChapter 6. Supplemental Text Material. Run, i X i1 X i2 X i1 X i2 Response total (1) a b ab
Chapter 6. Supplemental Text Materal 6-. actor Effect Estmates are Least Squares Estmates We have gven heurstc or ntutve explanatons of how the estmates of the factor effects are obtaned n the textboo.
More informationwhere I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X).
11.4.1 Estmaton of Multple Regresson Coeffcents In multple lnear regresson, we essentally solve n equatons for the p unnown parameters. hus n must e equal to or greater than p and n practce n should e
More informationLaboratory 1c: Method of Least Squares
Lab 1c, Least Squares Laboratory 1c: Method of Least Squares Introducton Consder the graph of expermental data n Fgure 1. In ths experment x s the ndependent varable and y the dependent varable. Clearly
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationScatter Plot x
Construct a scatter plot usng excel for the gven data. Determne whether there s a postve lnear correlaton, negatve lnear correlaton, or no lnear correlaton. Complete the table and fnd the correlaton coeffcent
More informationChapter 14 Simple Linear Regression
Chapter 4 Smple Lnear Regresson Chapter 4 - Smple Lnear Regresson Manageral decsons often are based on the relatonshp between two or more varables. Regresson analss can be used to develop an equaton showng
More informationUsing T.O.M to Estimate Parameter of distributions that have not Single Exponential Family
IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran
More informationChapter 5 Multilevel Models
Chapter 5 Multlevel Models 5.1 Cross-sectonal multlevel models 5.1.1 Two-level models 5.1.2 Multple level models 5.1.3 Multple level modelng n other felds 5.2 Longtudnal multlevel models 5.2.1 Two-level
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More information