확률분포 (4) 썸네일형 리스트형 [통계학|Scipy] scipy 로 정규분포 그래프 + 구간확률 구하기. 문제 : 모평균이 50.2이고 모분산이 48인 어떤 모집단으로부터 크기가 100인 표본이 임의추출되었다고 할 때 표본평균이 47에서 51 사이에 있을 확률을 구해보자. scipy.stats 에 보면, norm 이라는 클래스가 있다. 정규분포를 다룰 수 있는 클래스를 만들어준다. 정규분포는 평균과 표준편차 두 수치로 결정된다. norm 클래스 생성자 인자도 이 두개. 문제에서 표분평균의 표본분산은 모분사/표본크기 이고, 중심극한정리에 따라 정규분포를 따른다. N(x, m=502, sig^2=48/100) 을 구하면 된다. 구간의 확률은 norm 클래스의 cdf 함수의 차로 구할 수 있다. import numpy as np import matplotlib.pyplot as plt import scipy as.. [Python] 이항분포 그래프 그리기 import sys from functools import lru_cache import math import numpy as np import matplotlib.pyplot as plt @lru_cache(None) def ncr(n, r): """조합. 재귀식을 이용함. n이 커지면 스택오버플로우 발생. """ if r in (0, n): return 1 return ncr(n - 1, r) + ncr(n - 1, r - 1) def bidist(n, p): """이항분포""" q = 1 - p dist = np.array([ncr(n, k) * (p ** k) * (q ** (n - k)) for k in range(n + 1)]) return dist def bicoeff(n): """이항계수... Poisson Distribution E(x)의 테일러전개와 관계있는 확률분포. Poisson Distribution P(x ; λ) = ( λx / x! ) / eλ 이 확률분포는 전체 확률의 합이 1이 될까? ea 의 테일러 전개를 보자. ea = a0 / 0! + a1 / 1! + a2 / 2! + ... + an / n! + ... 위 식에서 a 를 λ 로 치환하고 보면, 좌변은 eλ 이고, 우변의 각항은 P(0;λ), P(1;λ), ... P(n;λ)... (에 비례하는 값)임을 알 수 있다. 거꾸로 말하자면, 푸아송 분포는 eλ의 테일러전개의 n번째 항의 비율을 확률로 갖는 분포이다. 이렇게 식을 외우면 좀 더 기억이 잘 되겠는데, 무슨 직관(=>see also 2)이 더 얻어질 수 있을까? 테일러전개에서 전체에 기여하는 항은 n ~ x 인 항과 그 주변 항들이다? see.. [MFC|CPP] 사구모양의 포텐셜 만들기 다음 신지식에 재미있어 보이는 게 있길래 답변달다가 만들어 본 것. 2차원 평면상에 주어진 점 P를 중심으로 둥그런 모양의 산을 만들어 보란다. 3차원 그래픽까지 구현할라면 죽을 것 같아서, 2차원 평면의 각 점에서의 함수값은 색으로 표현했다. 단색으로 표현할라니까 분해능이 256가지 뿐이다. 왼쪽 버튼을 클릭하면, 분지의 중앙점이 되는 P가 바뀌고, 오른쪽 버튼을 클릭하면 산이 바라보는 타겟 T가 바뀐다. 휠을 돌리면 분지의 반지름이 넓어진다. 소스코드다. VC6.0 에서 만들었다. 이전 1 다음
