数学建模代写  IMM Project math modeling

 IMM Project math modeling

数学建模代写 What are the major weaknesses of the SIR model for modeling Covid, or any epidemic? In what sort of situations is it particularly

In this project, you will analyze the SIR epidemic model, and apply it to real Covid-19 epidemiological data, for at least two regions. For the analysis of data for two regions, the Python notebook CovidDataExample.ipynb demonstrates how to read in data, extract time series for various regions, and make plots of the variables that enter the SIR model. The notebook SIR_data_fit.ipynb provides a lengthy discussion of the SIR model parameter and analysis, as well as code to find numerical solutions its deterministic form, as well as a fitting example for data for Japan.

Analysis of the SIR model  数学建模代写

For the analysis of the SIR model, provide the following:

1. What are the major weaknesses of the SIR model for modeling Covid, or any epidemic? In what sort of situations is it particularly well-suited?




2. How could the SIR model be modified to include separate compartments for dead and recovered populations? How about for vaccinated people? Write down the deterministic version of such a model and explain the new parameters.




3. News and government agencies often show graphs of “daily new cases”, which is the number of new positive cases per day (I’ll call it C(t)). Using variables from the model (S, I and R), what is this quantity?




4. How might you compute the replacement number rt without knowing all of the underlying model parameters? For example, there is a form that involves only the removal rate b. Alternatively, how might you simply calculate the exponential growth rate of I for some short period of time during which an outbreak is starting?




5. Draw the phase plane of the solution space for the SIR model, using just the equations for S and I. You may want to use the nondimensionalized population fractions s = S/N and ı = I/N. 6. Derive the peak number of infectives Imax, and if possible, estimate the time at which this peak occur. This may be done, for example, using the approximate solution derived in chapter 10 of Murray (2002, provided in the Resources folder).




6. Show that the peak in the daily new cases C(t) occurs before the peak in infections.

Analysis of data for two regions

For the two regions, you will choose:

1. A region that controlled the virus and exhibits a classic SIR-like curve for number of infected individuals (e.g. places like Hubei, or a cruise ship like the Diamond Princess)




2. A region with an uneven curve or recent outbreak, and hence an unresolved near future (e.g. Thailand).

For the simpler region (1), or a time period over which the infections reach a peak and then decline in the more challenging region (2), estimate at a minimum the following:  数学建模代写

1. The parameters a (average number of adequate contacts per person, per unit time) and b (the inverse average removal time). These can be found using the Scipy optimization/parameter fitting method curve_fit, as I I recommend that you leave N as a free variable in the fit.




2. The basic reproduction number r0 (compare this to the estimate r0 = 2.2 (e.g. here)




3. The population fraction that needs to be immuned in order to have herd immunity




4. If possible, consider including the timeseries of the vaccinated population in the model, based on your answer to Analysis question 2 above.

For the more challenging region (2), in a period during which there is a rapidly growing outbreak (e.g. for India in the last few weeks), estimate the exponential growth rate for I(t), and thereby estimate the time period over which the number of infected people is expected to double.

The form of the finished product should be a Python notebook, including both the analysis of the SIR model, as well as the analysis of the two regions. [However, if you prefer, you may provide the math of the analysis of the model in a separate PDF scan of your very neatly handwritten work. The text-based questions should be included in the notebook.]

Remember that JupyterHub cells can be used for text by selecting Markdown instead of Code; LaTeX can and should be used in these blocks for math. For helpusing Markdown in Jupyter, see this page, and for help with using LaTeX in general, here is a good cheatsheet. And of course, look at the many examples provided in all the shared notebooks.  数学建模代写

The final product should include:

• A title and project description: regions you are analyzing and why, a description of the data, caveats about this particular data (e.g. is the disease on the rise in these regions, or has the curve been flattened?), etc...
• Analysis of the SIR model, listed above
• Read-in of the data, with plots of I and R (and S, possibly on a separate plot, if it is much larger than the former), with labels and description.
• Perform fits for your data, using either the Gauss-Newton method, or another optimization package. Provide comparison plots for the fit and the data, for each of the regions you’ve analyzed.
• Compute the quanatities and analysis listed above.
• Conclude your work by interpreting the results. Why did different regions give different r0 alues? What does this say about containment strategies? Do you believe the results? Provide any references (including basic news articles) that you use to understand and interpret your results.

更多代写：计算机作业代写    经济代考   essay代写      AI作业代写     统计代考   经济paper代写