Summary and Discussion �Lab 6: �MLE fitting of a dynamic model

​

Reshma Kassanjee

University of Cape Town �SACEMA

​

MMED 2024

​

​

Some slides by Juliet Pulliam

​

​

��

​

​

000

Goals

• Understand how to simulate cross-sectional prevalence data around a derived and assumed epidemic trajectory
• Calculate the likelihood of prevalence data through time, given a fully specified epidemic model
• Understand that the likelihood is a function of the hypothesized model parameters (and data-generating process) and that “fitting” the dynamic model involves maximizing this likelihood
• Understand why we transform parameters for fitting
• Be aware that the choice of optimization algorithms affects the outcome of the optimization
• Create 95% confidence regions for a multivariate model fit

000

The real world

000

4

The ‘true’ prevalence over time

000

5

Simulate prevalence data

000

6

No longer know the true prevalence

🡪 Fit the model to data by estimating unknown parameters by � maximizing the likelihood function (using optimization algorithms)

000

The model world = real world

000

β × exp(-α×I/N)

The model world = real world

000

The ‘fitted’ prevalence over time

9

000

Parameter transformations

​

• Used log of the parameters when optimizing

​

Why?

​

• Most optimization algorithms assume the inputs are defined on a scale from –∞ to ∞

​

• More efficiently explore parameter space

​

000

Parameter transformations

000

Parameter transformations

000

Optimisation algorithms

​

• Carefully choose and modify the settings of your optimization algorithms

​

• Sometimes you may want to use multiple algorithms (and multiple starting values)

​

000

Parameter estimates

14

000

Parameter estimates

15

000

Parameter estimates

16

Uncertainties??

000

Normal distribution

17

000

Likelihood ratio

18

000

Take-away points

• Always fit fake data (e.g., simulated by your model) before attempting to fit real data: to validate the fitting approach and identify any errors

​

• Think about your optimization algorithms and what they are doing, and choose an approach that will be both efficient and accurate – this may include transforming parameters to scales that span the real numbers

​

000

This presentation is made available through a Creative Commons Attribution-Noncommercial license. Details of the license and permitted uses are available at� http://creativecommons.org/licenses/by-nc/3.0/

Title: Lab 6 Summary: MLE fitting of a dynamic model to prevalence data

Attribution: Juliet R.C. Pulliam & Reshma Kassanjee & Steve E. Bellan, Clinic on the Meaningful Modeling of Epidemiological Data

Source URL: http://www.ici3d.org/MMED/tutorials/Lab6_summary.pdf

For further information please contact admin@ici3d.org.

© 2024 International Clinics on Infectious Disease Dynamics and Data

000