nonlinear regression: untransformed data

17
1 Nonlinear Regression: Untransformed Data http://nm.MathForCollege.com Transforming Numerical Methods Education for STEM Undergraduates

Upload: others

Post on 18-Oct-2021

6 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Nonlinear Regression: Untransformed Data

1

Nonlinear Regression: Untransformed Data

http://nm.MathForCollege.comTransforming Numerical Methods Education for STEM Undergraduates

Page 2: Nonlinear Regression: Untransformed Data

Nonlinear Regression

𝑦 π‘Žπ‘’π‘¦ π‘Žπ‘₯

π‘¦π‘Žπ‘₯𝑏 π‘₯

Some popular nonlinear regression models:

1. Exponential model:

2. Power model:

3. Saturation growth model:

4. Polynomial model: 𝑦 π‘Ž0 π‘Ž1π‘₯ . . . π‘Žπ‘šπ‘₯

http://numericalmethods.eng.usf.edu2

Page 3: Nonlinear Regression: Untransformed Data

Nonlinear Regression

Given n data points (x1,y1), …… (xn,yn), best fit 𝑦 𝑓 π‘₯

to the data, where 𝑓 π‘₯ is a nonlinear function of π‘₯

Figure. Nonlinear regression model for discrete y vs. x data

𝑦 𝑓 π‘₯

π‘₯ , 𝑦

π‘₯ , 𝑦

π‘₯ ,𝑦

π‘₯ , 𝑦𝑦 𝑓 π‘₯

http://numericalmethods.eng.usf.edu3

Page 4: Nonlinear Regression: Untransformed Data

Exponential Model

π‘₯ , 𝑦 , π‘₯ , 𝑦 , ... , π‘₯ ,𝑦Given best fit 𝑦 π‘Žπ‘’ to the data.

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

𝑦 π‘Žπ‘’

π‘₯ , 𝑦

π‘₯ ,𝑦

π‘₯ ,𝑦

π‘₯ ,𝑦 𝑦 π‘Žπ‘’

http://numericalmethods.eng.usf.edu4

Page 5: Nonlinear Regression: Untransformed Data

Finding Constants of Exponential Model

π‘†π‘Ÿ 𝑦 π‘Žπ‘’

The sum of the square of the residuals is defined as

Differentiate with respect to a and b

πœ•π‘†πœ•π‘Ž 2 𝑦 π‘Žπ‘’ 𝑒 0

πœ•π‘†πœ•π‘ 2 𝑦 π‘Žπ‘’ π‘Žπ‘₯ 𝑒 0

http://numericalmethods.eng.usf.edu5

Page 6: Nonlinear Regression: Untransformed Data

Finding Constants of Exponential Model

Rewriting the equations, we obtain

𝑦 𝑒 π‘Ž 𝑒 0

𝑦 π‘₯ 𝑒 π‘Ž π‘₯ 𝑒 0

http://numericalmethods.eng.usf.edu6

Page 7: Nonlinear Regression: Untransformed Data

Finding constants of Exponential Model

𝑦 π‘₯ π‘’βˆ‘ 𝑦 π‘’βˆ‘ 𝑒 π‘₯ 𝑒 0

π‘Žβˆ‘ 𝑦 π‘’βˆ‘ 𝑒

Solving the first equation for a yields

http://numericalmethods.eng.usf.edu7

𝑦 π‘₯ 𝑒 π‘Ž π‘₯ 𝑒 0

Page 8: Nonlinear Regression: Untransformed Data

Example 1-Exponential Model

t(hrs) 0 1 3 5 7 91.000 0.891 0.708 0.562 0.447 0.355

Many patients get concerned when a test involves injection of a radioactive material. For example for scanning a gallbladder, a few drops of Technetium-99m isotope is used. Half of the Technetium-99m 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.

𝛾

http://numericalmethods.eng.usf.edu8

Page 9: Nonlinear Regression: Untransformed Data

Example 1-Exponential Model cont.

Find: a) The value of the regression constants 𝐴 πœ†andb) The half-life of Technetium-99mc) Radiation intensity after 24 hours

The relative intensity is related to time by the equation𝛾 𝐴𝑒

http://numericalmethods.eng.usf.edu9

Page 10: Nonlinear Regression: Untransformed Data

Plot of data

http://numericalmethods.eng.usf.edu10

Page 11: Nonlinear Regression: Untransformed Data

Constants of the Model

The value of Ξ» is found by solving the nonlinear equation

𝑓 πœ† 𝛾 𝑑 π‘’βˆ‘ 𝛾 π‘’βˆ‘ 𝑒

𝑑 𝑒 0

π΄βˆ‘ 𝛾 π‘’βˆ‘ 𝑒

𝛾 𝐴𝑒

http://numericalmethods.eng.usf.edu11

Page 12: Nonlinear Regression: Untransformed Data

Plotting the Equation in MATLAB

t (hrs) 0 1 3 5 7 9Ξ³ 1.000 0.891 0.708 0.562 0.447 0.355

http://numericalmethods.eng.usf.edu12

Page 13: Nonlinear Regression: Untransformed Data

Setting up the Equation in MATLAB

t=[0 1 3 5 7 9]gamma=[1 0.891 0.708 0.562 0.447 0.355]syms lamdasum1=sum(gamma.*t.*exp(lamda*t));sum2=sum(gamma.*exp(lamda*t));sum3=sum(exp(2*lamda*t));sum4=sum(t.*exp(2*lamda*t));f=sum1-sum2/sum3*sum4;

http://numericalmethods.eng.usf.edu13

Page 14: Nonlinear Regression: Untransformed Data

Calculating the Other Constant

The value of A can now be calculated

π΄βˆ‘ 𝛾 π‘’βˆ‘ 𝑒

0.9998

The exponential regression model then is𝛾 0.9998 𝑒 .

http://numericalmethods.eng.usf.edu14

Page 15: Nonlinear Regression: Untransformed Data

Plot of data and regression curve

.

http://numericalmethods.eng.usf.edu15

Page 16: Nonlinear Regression: Untransformed Data

Relative Intensity After 24 hrs

The relative intensity of radiation after 24 hours 𝛾 0.9998 𝑒 .

6.3160 10

This result implies that only6.316 10

0.9998 100 6.317%

radioactive intensity is left after 24 hours.

http://numericalmethods.eng.usf.edu16

Page 17: Nonlinear Regression: Untransformed Data

Homework

β€’ What is the half-life of Technetium-99m isotope?

β€’ Write a program in the language of your choice to find the constants of the model.

β€’ Compare the constants of this regression model with the one where the data is transformed.

β€’ What if the model was

http://numericalmethods.eng.usf.edu17

𝛾 𝑒