반응형
/*******************************************************************************************************************
-- Title : [Stats] 이항분포 함수 - dbinom(), pbinom(), qbinom(), rbinom()
-- Reference : http://blog.naver.com/padosori60/220755436419
-- Key word : 통계 statistics dbinom rbinom qbinom pbinom 이항 분포 함수 이항분포 이항분포함수
*******************************************************************************************************************/
■ 이항 분포 확률
■ 이항 분포 함수
■ dbinom() 함수 : 이항분포의 확률값
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
|
# X가 이항분포 B(10, 0.2)일 경우, P[X=2]의 확률값은?
dbinom(2, 10, 0.2)
# [1] 0.3019899
a = seq(0, 50, 1)
b = dbinom(a, 50, 0.2)
plot(a,b)
b = dbinom(a, 50, 0.6)
plot(a,b)
a = seq(0, 100, 1)
b = dbinom(a, 100, 0.6)
plot(a, b)
|
cs |
■ pbinom() 함수 : 이항분포의 누적확률값
|
1
2
3
4
5
6
7
8
9
10
|
# X가 이항분포 B(10, 0.2)이고 P[X<=3]인 누적확률값을 구하면
pbinom(3, 10, 0.2)
# [1] 0.8791261
dbinom(0, 10, 0.2) # 0.1073742
dbinom(1, 10, 0.2) # 0.2684355
dbinom(2, 10, 0.2) # 0.3019899
dbinom(3, 10, 0.2) # 0.2013266
0.1073742 + 0.2684355 + 0.3019899 + 0.2013266 # 0.8791262 = pbinom(3, 10, 0.2)와 유사
|
cs |
■ qbinom() 함수 : 이항분포의 백분위수
|
1
2
|
# qbinom( )함수는 백분위수를 구하는 함수로, pbinom( )함수의 역의 값
qbinom(0.6, 10, 0.2) # [1] 2
|
cs |
■ rbinom() 함수 : 이항분포의 난수발생
|
1
2
|
rbinom(5, 100, 0.3) # [1] 33 29 32 35 39
rbinom(5, 100, 0.7) # [1] 70 65 63 62 70
|
cs |
반응형