доверительный интервал для самих данных с использованием lmfit в python

введите здесь описание изображения

Вот ссылка на реализацию LMFIT доверительных интервалов параметров: http://lmfit.github.io/lmfit-py/confidence.html Вот код, который я использую:

import lmfit
import numpy as np

# x = np.linspace(1, 10, 250)
# np.random.seed(0)
# y = 1. -np.exp(-(x)/10.) +  0.1*np.random.randn(len(x))

pars = lmfit.Parameters()
pars.add_many(('n', 1.),  ('tau', 3.))

# def residual(pars,data=None):
def residual(pars):
    v = pars.valuesdict()
    # if data is None:
    #     return 1.0 - np.exp(-(x**v['n'])/v['tau'])
    return 1.0 - np.exp(-(x**v['n'])/v['tau'])-y

# create Minimizer
mini = lmfit.Minimizer(residual, pars)

# first solve with Nelder-Mead
out1 = mini.minimize(method='Nelder')

out2 = mini.minimize(method='leastsq', params=out1.params)

lmfit.report_fit(out2.params, min_correl=0.5)

ci, trace = lmfit.conf_interval(mini, out2, sigmas=[0.95],
                                trace=True, verbose=False)
lmfit.printfuncs.report_ci(ci)

person bazilevs31    schedule 21.10.2015    source источник
comment
Не могли бы вы отредактировать свой вопрос? Первая часть текста вырезана; вторую часть невозможно читать, спасибо!   -  person Cleb    schedule 22.10.2015


Ответы (1)


Немного сложно понять заголовок Доверительный интервал для самих данных с использованием lmfit в python (нет данных) или первое предложение Я выполняю подбор кривой с помощью пакета lmfit< /em> (вам нужны данные, чтобы соответствовать).

Я думаю, что вы просите, чтобы получить экстремальные значения для функции модели, которая лучше всего соответствует вашим данным. Если да, получится ли оценить вашу функцию со всеми комбинациями значений параметров best +/- delta (где delta может быть любым уровнем неопределенности, который вам нравится) и взять экстремальные значения функции модели? Это не очень автоматизировано, но не должно быть слишком сложно.

person M Newville    schedule 24.10.2015
comment
Он ищет полосы уверенности: en.wikipedia.org/wiki/Confidence_and_prediction_bands. - person Christian K.; 26.10.2015