Nonlinear Regression Without Transformation of Data 1 Nonlinear

  • Slides: 19
Download presentation
Nonlinear Regression Without Transformation of Data 1

Nonlinear Regression Without Transformation of Data 1

Nonlinear Regression 2 http: //nm. Math. For. College. com

Nonlinear Regression 2 http: //nm. Math. For. College. com

Nonlinear Regression Figure. Nonlinear regression model for discrete y vs. x data 3 http:

Nonlinear Regression Figure. Nonlinear regression model for discrete y vs. x data 3 http: //nm. Math. For. College. com

Exponential Model Figure. Exponential model of nonlinear regression for y vs. x data 4

Exponential Model Figure. Exponential model of nonlinear regression for y vs. x data 4 http: //nm. Math. For. College. com

Finding Constants of Exponential Model The sum of the square of the residuals is

Finding Constants of Exponential Model The sum of the square of the residuals is defined as 5 http: //nm. Math. For. College. com

Finding Constants of Exponential Model Solving the first equation for a yields 6 http:

Finding Constants of Exponential Model Solving the first equation for a yields 6 http: //nm. Math. For. College. com

Example - Exponential Model Many patients get concerned when a test involves injection of

Example - Exponential Model Many patients get concerned when a test involves injection of a radioactive material. For example for scanning a gallbladder, a few drops of Technetium 99 m isotope is used. Half of the Technetium-99 m would be gone in about 6 hours. It, however, takes about 24 hours for the radiation levels to reach what we are exposed to in day-to-day activities. Below is given the relative intensity of radiation as a function of time. Table. Relative intensity of radiation as a function of time. t(hrs) 7 0 1 3 5 7 9 1. 000 0. 891 0. 708 0. 562 0. 447 0. 355 http: //nm. Math. For. College. com

Example - Exponential Model (contd) Table. Relative intensity of radiation as a function of

Example - Exponential Model (contd) Table. Relative intensity of radiation as a function of time. 8 0 1 3 5 7 9 1. 000 0. 891 0. 708 0. 562 0. 447 0. 355 http: //nm. Math. For. College. com

Plot of data 9 http: //nm. Math. For. College. com

Plot of data 9 http: //nm. Math. For. College. com

Constants of the Model The value of λ is found by solving the nonlinear

Constants of the Model The value of λ is found by solving the nonlinear equation 10 http: //nm. Math. For. College. com

Setting up the Equation in MATLAB 11 0 1 3 5 7 9 1.

Setting up the Equation in MATLAB 11 0 1 3 5 7 9 1. 000 0. 891 0. 708 0. 562 0. 447 0. 355 http: //nm. Math. For. College. com

Setting up the Equation in MATLAB t=[0 1 3 5 7 9] gamma=[1. 000

Setting up the Equation in MATLAB t=[0 1 3 5 7 9] gamma=[1. 000 0. 891 0. 708 0. 562 0. 447 0. 355] syms lamda sum 1=sum(gamma. *t. *exp(lamda*t)); sum 2=sum(gamma. *exp(lamda*t)); sum 3=sum(exp(2*lamda*t)); sum 4=sum(t. *exp(2*lamda*t)); f=sum 1 -sum 2/sum 3*sum 4; 12 http: //nm. Math. For. College. com

Calculating the Other Constant The value of �� can now be calculated The exponential

Calculating the Other Constant The value of �� can now be calculated The exponential regression model then is 13 http: //nm. Math. For. College. com

Plot of data and regression curve 14 http: //nm. Math. For. College. com

Plot of data and regression curve 14 http: //nm. Math. For. College. com

Relative Intensity After 24 hrs 15 http: //nm. Math. For. College. com

Relative Intensity After 24 hrs 15 http: //nm. Math. For. College. com

Homework 16 http: //nm. Math. For. College. com

Homework 16 http: //nm. Math. For. College. com

Transformed vs Untransformed Data 0 1 3 5 7 9 1. 000 0. 891

Transformed vs Untransformed Data 0 1 3 5 7 9 1. 000 0. 891 0. 708 0. 562 0. 447 0. 355

Transformed vs Untransformed Data x y 0 1. 0000 5 0. 8326 10 0.

Transformed vs Untransformed Data x y 0 1. 0000 5 0. 8326 10 0. 6738 15 0. 5837 20 0. 5150 25 0. 4163 40 0. 3219 60 0. 2466 90 0. 1803