У меня есть 1-мерный массив данных:
a = np.array([1,2,3,4,4,4,5,5,5,5,4,4,4,6,7,8])
для которого я хочу получить доверительный интервал 68% (т.е. 1 сигма).
Первый комментарий в этом ответе гласит, что этого можно достичь с помощью scipy.stats.norm.interval
из scipy.stats.norm, используя:
from scipy import stats
import numpy as np
mean, sigma = np.mean(a), np.std(a)
conf_int = stats.norm.interval(0.68, loc=mean, scale=sigma)
Но комментарий в этот пост утверждает, что фактический правильный способ получения доверительного интервала:
conf_int = stats.norm.interval(0.68, loc=mean, scale=sigma / np.sqrt(len(a)))
т.е. на сигме используется коэффициент 1/np.sqrt(len(a))
.
Вопрос: какая версия правильная?