Como calcular a probabilidade de ajuste de curva no scipy?
Eu tenho um modelo não-linear que se parece com isso:
A linha sólida escura é adequada ao modelo e a parte cinza são os dados brutos.
Versão curta da pergunta: como obtenho a probabilidade desse modelo para que eu possa executar o teste da razão de verossimilhança de log? Suponha que o resíduo seja normalmente distribuído.
Sou relativamente novo em estatística, e meus pensamentos atuais são:
Obtenha o residual do ajuste da curva e calcule a variação do residual;
Use esta equação E conecte a variação do residual ao quadrado sigma, x_i como experimento e mu conforme o modelo;
Calcule a razão de log-verossimilhança.
Alguém poderia me ajudar, com essas duas perguntas da versão completa?
Meu método está correto? (Acho que sim, mas seria ótimo ter certeza!)
Existem funções prontas no python / scipy / statsmodels para fazer isso por mim?