pymc

Dopasowanie funkcji nieliniowej do danych/obserwacji za pomocą pyMCMC/pyMC

Próbuję dopasować niektóre dane za pomocą funkcji Gaussa (i bardziej złożonych). Poniżej stworzyłem mały przykład. Moje pie ... duje się poniżej. Pymc.Matplot.wykresy (MDL) wyglądają tak, pokazując ładnie rozłożone rozkłady. To jest dobre, prawda?

pymc3: wiele obserwowanych wartości

Mam kilka danych obserwacyjnych, dla których chciałbym oszacować parametry, i pomyślałem, że będzie to dobra okazja, aby wypr ... opolis() success_samples = sample(20000, step) #, start) plt.figure(figsize=(10, 10)) _ = traceplot(success_samples)

Region o najwyższym zagęszczeniu tylnym i Regionie Centralnym

Biorąc pod uwagę tylną p (Θ / D) nad niektórymi parametrami Θ, można zdefiniować następująco: Najwyższy Tylny Obszar Gęst ... trycznych(np. Beta, Gaussa itp.) czy są jakieś wbudowane lub biblioteki do obliczenia tego za pomocą SciPylub statsmodels?