본문 바로가기

Python/통계4

정규분포 with python 정규분포 확률밀도함수 정규분포의 확률밀도함수는 다음과 같습니다. $n(x;\mu, \sigma) = \frac {1}{\sqrt {2\pi}\sigma} e^{-\frac {1}{2\sigma^{2}}(x-\mu)^{2}}$, $-\infty 정규분포 평균 변형 legend = [] def norm (mu, var) : legend.append("N(" + str(mu) + ", " + str(var) + ")") return stats.norm(mu, var).pdf(x) plt.figure(figsize=(8, 6)) # 플롯 사이즈 지정 for i in np.arange(4) + 1: plt.plot(x, norm(i, 1)) # plot 추가 plt.xlabel("x") # x축 레이블 지정 pl.. 2019. 10. 15.
감마분포 with python 감마분포 확률밀도함수 감마분포의 확률밀도함수는 다음과 같습니다. $f(x; \alpha, \beta) = \frac {1}{\beta^{\alpha}\Gamma(\alpha)}x^{\alpha - 1} e^{-\frac {x}{\beta}}$ ($x > 0$) 감마분포 그래프 그리기(함수 사용) from numpy import linspace, exp # numpy.linspace(start, stop, num=50, endpoint=True, retstep=False, dtype=None, axis=0) # numpy.exp(x, /, out=None, *, where=True, casting='same_kind', order='K', dtype=None, subok=True[, signature, e.. 2019. 10. 14.
(2)정규분포의 누적분포함수를 python으로 구현해보자 def normal_cdf(x, mu=0,sigma=1): return (1 + math.erf((x - mu) / math.sqrt(2) / sigma)) / 2 정규분포의 누적분포함수는 오차함수를 사용하여 구현합니다. (erf: error function) 오차함수의 설명은 아래를 참고하세요(나중에 설명드리겠습니다) https://ko.wikipedia.org/wiki/%EC%98%A4%EC%B0%A8_%ED%95%A8%EC%88%98 오차 함수 - 위키백과, 우리 모두의 백과사전 위키백과, 우리 모두의 백과사전. ko.wikipedia.org 이 글의 출처는 data science from scratch o'reilly입니다. 너무나 좋은 책이니 꼭 한 번 사서 공부해보셨으면 좋겠습니다. 2019. 8. 22.
(1)정규분포 밀도함수 python 함수로 만들어보기 def normal_pdf(x, mu=0, sigma=1): sqrt_two_pi = math.sqrt(2 * math.pi) return (math.exp(-(x-mu) ** 2 / 2 / sigma ** 2) / (sqrt_two_pi * sigma)) ※ 정규분포의 밀도함수 $\mathit{f}(x | \mu,\sigma) = \frac{1}{\sqrt{2\pi \sigma }}exp(\frac{(x - \mu)^{2}}{2\sigma ^{2}})$ ($-\infty \leq x \leq \infty$) 이 글의 출처는 data science from scratch o'reilly입니다. 2019. 8. 20.