ECE3340 Numerical Methods

Classwork 3/26 - solutions
US Covid-19 data analysis

ECE 3340  & generic    
Han Q. Le (c)

Q. 1 Linear regression of simulated data

In[3]:=

Classwork_3_26_solutions_1.png

Out[3]=

Classwork_3_26_solutions_2.png

We want to create a variable y1= a x1 +b

In[4]:=

Classwork_3_26_solutions_3.gif

Out[4]=

Classwork_3_26_solutions_4.png

Out[5]=

Classwork_3_26_solutions_5.png

Out[6]=

Classwork_3_26_solutions_6.png

Classwork_3_26_solutions_7.png

Out[5]=

Classwork_3_26_solutions_8.png

Classwork_3_26_solutions_9.png

Out[6]=

Classwork_3_26_solutions_10.gif

In[7]:=

Classwork_3_26_solutions_11.gif

Out[8]=

Classwork_3_26_solutions_12.png

In[9]:=

Classwork_3_26_solutions_13.png

Out[9]=

Classwork_3_26_solutions_14.png

In[10]:=

Classwork_3_26_solutions_15.gif

Out[10]=

Classwork_3_26_solutions_16.png

Out[11]=

Classwork_3_26_solutions_17.gif

Out[12]=

Classwork_3_26_solutions_18.png

Out[13]=

Classwork_3_26_solutions_19.png

Out[14]=

Classwork_3_26_solutions_20.png

Out[15]=

Classwork_3_26_solutions_21.gif

Q. 2 Estimates of coefficients

In[7]:=

Classwork_3_26_solutions_22.gif

In[45]:=

Classwork_3_26_solutions_23.png

In[21]:=

Classwork_3_26_solutions_24.png

Out[21]=

Classwork_3_26_solutions_25.gif

Q. 3 Estimate coefficients: descriptive statistics

In[10]:=

Classwork_3_26_solutions_26.gif

Out[10]=

Classwork_3_26_solutions_27.png

Out[11]=

Classwork_3_26_solutions_28.png

In[12]:=

Classwork_3_26_solutions_29.png

Out[12]=

Classwork_3_26_solutions_30.gif

In[13]:=

Classwork_3_26_solutions_31.png

Out[13]=

Classwork_3_26_solutions_32.gif

In[14]:=

Classwork_3_26_solutions_33.png

Out[14]=

Classwork_3_26_solutions_34.gif

Q. 4 Estimate coefficients: Histogram and model

In[15]:=

Classwork_3_26_solutions_35.gif

Out[19]=

Classwork_3_26_solutions_36.gif

In[22]:=

Classwork_3_26_solutions_37.png

Out[22]=

Classwork_3_26_solutions_38.gif

Data covid

In[23]:=

Classwork_3_26_solutions_39.png

In[26]:=

Classwork_3_26_solutions_40.png

Out[26]=

Classwork_3_26_solutions_41.png

In[24]:=

Classwork_3_26_solutions_42.gif

In[26]:=

Classwork_3_26_solutions_43.gif

Q.5  US Covid data

In[31]:=

Classwork_3_26_solutions_44.gif

Out[31]=

Classwork_3_26_solutions_45.png

In[54]:=

Classwork_3_26_solutions_46.gif

Out[55]=

Classwork_3_26_solutions_47.gif

Out[56]=

Classwork_3_26_solutions_48.gif

In[36]:=

Classwork_3_26_solutions_49.png

Out[36]=

Classwork_3_26_solutions_50.png

In[37]:=

Classwork_3_26_solutions_51.png

Out[37]=

Classwork_3_26_solutions_52.png

In[47]:=

Classwork_3_26_solutions_53.gif

In[49]:=

Classwork_3_26_solutions_54.png

Out[49]=

Classwork_3_26_solutions_55.gif

In[50]:=

Classwork_3_26_solutions_56.png

Out[50]=

Classwork_3_26_solutions_57.png

In[51]:=

Classwork_3_26_solutions_58.png

Out[51]=

Classwork_3_26_solutions_59.png

In[57]:=

Classwork_3_26_solutions_60.png

Out[57]=

Classwork_3_26_solutions_61.gif

In[53]:=

Classwork_3_26_solutions_62.png

Out[53]=

Classwork_3_26_solutions_63.png

In[58]:=

Classwork_3_26_solutions_64.png

Out[58]=

Classwork_3_26_solutions_65.gif

In[59]:=

Classwork_3_26_solutions_66.png

Out[59]=

Classwork_3_26_solutions_67.gif

Q.6  Confidence of the estimates

In[33]:=

Classwork_3_26_solutions_68.png

Out[33]=

Classwork_3_26_solutions_69.png

In[41]:=

Classwork_3_26_solutions_70.png

Out[41]=


Estimate Standard Error t-Statistic P-Value
1 11.2092 0.0416336 269.235 9.58726*10^^-42
x 0.286483 0.00280058 102.294 4.35642*10^^-32

In[60]:=

Classwork_3_26_solutions_71.png

Out[60]=

Classwork_3_26_solutions_72.png

In[61]:=

Classwork_3_26_solutions_73.png

Out[61]=

Classwork_3_26_solutions_74.png

In[62]:=

Classwork_3_26_solutions_75.png

Out[62]=

Classwork_3_26_solutions_76.png

In[63]:=

Classwork_3_26_solutions_77.png

Out[63]=

Classwork_3_26_solutions_78.gif

Q.7  Projection

In[64]:=

Classwork_3_26_solutions_79.png

Out[64]=

Classwork_3_26_solutions_80.png

In[65]:=

Classwork_3_26_solutions_81.png

Out[65]=

Classwork_3_26_solutions_82.gif

In[83]:=

Classwork_3_26_solutions_83.png

Out[83]=

Estimates of infected population in the US
(if no preventive measure and same infection rate)

date low estimate high estimate
Fri Mar, 27 89567 105742
Sat Mar, 28 118650 141412
Sun Mar, 29 157166 189126
Mon Mar, 30 208173 252955
Tue Mar, 31 275719 338343
Wed Apr, 01 365165 452574
Thu Apr, 02 483612 605396

Q.8  Trend

Daily trend

In[84]:=

Classwork_3_26_solutions_84.png

Out[84]=

Classwork_3_26_solutions_85.png

In[85]:=

Classwork_3_26_solutions_86.png

Out[85]=

Classwork_3_26_solutions_87.png

In[86]:=

Classwork_3_26_solutions_88.png

Out[86]=

Classwork_3_26_solutions_89.gif

In[87]:=

Classwork_3_26_solutions_90.png

Out[87]=

Classwork_3_26_solutions_91.gif

In[88]:=

Classwork_3_26_solutions_92.gif

Out[89]=

Classwork_3_26_solutions_93.gif

In[90]:=

Classwork_3_26_solutions_94.gif

Out[91]=

Classwork_3_26_solutions_95.gif

For HW - Logistic curve fit for most recent total cases

Out[91]=

Classwork_3_26_solutions_96.gif

Is the slight curvature significant? Consider model:  Classwork_3_26_solutions_97.png
We can fit to the most recent trend to see if the bend over is real

In[114]:=

Classwork_3_26_solutions_98.gif

In[207]:=

Classwork_3_26_solutions_99.png

Out[207]=

Classwork_3_26_solutions_100.gif

In[236]:=

Classwork_3_26_solutions_101.gif

Out[244]=

the turning point in the US:   Mar, 25
estimate maximum infections eventually 129084.

In[216]:=

Classwork_3_26_solutions_102.png

Out[216]=

Classwork_3_26_solutions_103.gif

In[218]:=

Classwork_3_26_solutions_104.png

Out[218]=

Classwork_3_26_solutions_105.png

For the best case scenario, the maximum infection in the US will max out between 108 K and 150 K if the bending trend in recent data is true. Below is the equivalent graph on the App (Log scale)

Classwork_3_26_solutions_106.gif

Q.9  Spatial Laplace  distribution

Link to HW 5





Created with the Wolfram Language