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?

questionAnswers(2)

yourAnswerToTheQuestion