In order to estimate parameter values for our model we explored both Markov Chain Monte Carlo and Maximum Likelihood Estimation techniques. Here we present the results for the five parameter model as we believe it to be the most accurate.
Markov Chain Monte Carlo:
Histograms and Trace Plots:
Maximum Likelihood Estimation:
In the model fits as a geographical parameter in the spatial kernel: this is used to weight the mixing component of the model to act heterogeneously.
As is an exponential function with range in (0,1) it is easy to measure the rate of decay:
The diameter of the Birmingham district we are working over is is 55km so the above data shows us the disease is very well localized spatially.
As stated in our model the probability a school i is infected at time t is given by This form allows us to study the effects of on the model as the probabilities factor implying independence. As smaller implies a greater chance of infection. We can then plot this drop off as a function of school size:
As the average school size in Birmingham is 227.5 we can attribute between 80% to 90% of infections occurring due to school mixing as oppose to external pressure.
The and parameter:
As these parameters are used in a negative binomial they are fairly easy to interpret.
Mean time is school is infectious: 4.3225 Weeks
Standard deviation: 2.3976 Weeks
As the time scale is of 9 weeks this shows us that once a school gets infected it generally stays that way for the majority of the epidemic.