silpara / simulators

Simulators for Compartmental Models in Epidemiology
37 stars 38 forks source link

SEIRD Simulator in Python #1

Closed pinkfloydisch closed 4 years ago

pinkfloydisch commented 4 years ago

Hi, Thanks a lot for very nice, clean and useful code. surprisingly it is exactly the mode which I am working on it. I am studying mathematics and it is very important for me to estimate the parameters, especially beta, May I ask you to write additional code for SEIRD model regarding the least square error and minimizing it which leads to estimate the parameters. I am working on it but I have problem in coding part of it. the problem that I have is that maybe separately I or D are fitted to the data but all together usually no. my first guess for beta is something between 0.5 to 1 , sigma 1 / 5.6 , gamma = 1/14 and mu = 0.02

the problem is that it is difficult for me how to write a code which I can use real data of infected and death people (even recovered) at the same time into my code to minimize the least square in several aspects and estimate the parameters better.

Sorry for the long message. I thought maybe you can help me. have a nice and safe days Best regards, Nik

silpara commented 4 years ago

Hi Nik,

It was in the plan. Have published a post: Estimating Parameters of Compartmental Models from Observed Data

Here are the notebooks demonstrating the parameter estimation Simulation & Parameter Estimation of SEIR Model Simulation & Parameter Estimation of SEIRD Model

I have no background in medicine or epidemiology, neither this code has been vetted by any expert, so be careful for what purpose you use it. Ask experts to you know to validate the results if you are planning to use them one way or the other.

Cheers, Abhay