Estimação dos parametro about modelagem-matematica-iv

[fccoelho]

@wellington36 , tira o self da função SIR. vai ficar:

def SIR(y,t,p):
    ...

[fccoelho]

@wellington36, este log retorna um vetor, ou seja, temos um mu para cada ponto da curva.

Acho que o seu erro é que vc está usando dados definidos como uma lista:

yobs = [0.15, 0.10, 0.04, 0.009, 0.001]

Quando a saida do odeint retorna um vetor coluna:

# Gerando curvas simuladas
y = odeint(SIR, t=times, y0=[0.99, 0.01], args=((beta, gamma),), rtol=1e-8)
# Simulando dados  Assumindo uma distribuição log-normal com média igual às séries simuladas
yobs = np.random.lognormal(mean=np.log(y[1::]), sigma=[0.2, 0.3])

acho que o erro vem disso, o shape da série não bate com o esperado.

[wellington36]

@fccoelho, uma outra questão que estou com dúvida é que no jupyter do curso: Aula 14 - Estimando Parâmetros.ipynb, mais expecificamente no seguinte ponto:

o senhor faz o log do que seria uma curva não entendi como essa curva pode retornar um número (no meu caso isso da um erro)

Qualquer coisa esté é meu jupyter com as imprementações até o momento: https://github.com/wellington36/Modelling-the-dynamics-of-language-death/blob/main/implementa%C3%A7%C3%A3o.ipynb

[wellington36]

De fato, resolveu o bug.

[wellington36]

@fccoelho, sem querer encher o saco do senhor acho que agora sim achei algo problematico, na celula, do mesmo jupyter mencionado acima temos:

Gerando este erro, que me parece uma questão de precisão da maquina:

[fccoelho]

vc não deve usar distribuição normal par parâmetros que são estritamente positivos (neste caso o a). Joga um pm.Bound, como foi feito nos outros, para truncar a distribuição.

[wellington36]

Obrigado, agora ta dando certo: