Acelerar run_SEIR_BAYES_model

[giulianonetto]

Em vez de fazer um for loop ao longo de runs iterações, geramos runs observações a partir dos random number generators e fazemos a extrapolação temporal com as matrizes/vetores diretamente.

3-4x speed up: https://github.com/giulianonetto/COVID-19/blob/master/models/speed%20up%20run_SEIR.ipynb

[dsevero]

@giulianonetto Consegues fazer uns plots só pra ver se os dois tão batendo? Parece que ta certo, mas sempre bom conferir.

[giulianonetto]

Parece que as matrizes geram medias iguais para cada tempo t:

Adicionei no notebook distribuicao dos quantis por tempo t em cada matriz gerada tbm, e.g.:

Se quiser dar uma olhada no notebook, tem o codigo todo (o mesmo que botei o link antes).