R에서 이항 분포를 그리는 방법

R의 이항 분포 에 대한 확률 질량 함수를 플롯하려면 다음 함수를 사용할 수 있습니다.

• 확률질량함수를 생성하기 위한 dbinom(x, size, prob)
• 플롯(x, y, type = ‘h’) – 확률 질량 함수를 플롯하고 플롯이 히스토그램(type=’h’)임을 지정합니다.

확률 질량 함수를 플롯하려면 크기를 지정하기만 하면 됩니다.   dbinom() 함수의 (예: 시행 횟수) 및 prob (예: 주어진 시행의 성공 확률).

예를 들어, 다음 코드는 크기 = 20이고 prob = 0.3인 이항 분포에 대한 확률 질량 함수를 그리는 방법을 보여줍니다.

 success <- 0:20

plot(success, dbinom(success, size=20, prob=.3), type='h')

x축은 성공 횟수를 나타내고 y축은 20번의 시도에서 해당 횟수의 성공 확률을 나타냅니다.

제목을 추가하고, 축 레이블을 변경하고, 선의 너비를 늘려 플롯을 더 미적으로 보기 좋게 만들 수 있습니다.

 success <- 0:20

plot(success,dbinom(success,size=20,prob=.3),
     type='h',
     main='Binomial Distribution (n=20, p=0.3)',
     ylab='Probability',
     xlab ='# Successes',
     lwd=3)

다음 코드를 사용하여 차트에 표시된 각 성공 횟수에 대한 실제 확률을 얻을 수 있습니다.

 #prevent R from displaying numbers in scientific notation
options(scipen=999) 

#define range of successes
success <- 0:20

#display probability of success for each number of trials
dbinom(success, size=20, prob=.3)

[1] 0.00079792266297612 0.00683933711122388 0.02784587252426865
[4] 0.07160367220526231 0.13042097437387065 0.17886305056987975
[7] 0.19163898275344257 0.16426198521723651 0.11439673970486122
[10] 0.06536956554563482 0.03081708090008504 0.01200665489613703
[13] 0.00385928193090119 0.00101783259716075 0.00021810698510587
[16] 0.00003738976887529 0.00000500755833151 0.00000050496386536
[19] 0.00000003606884753 0.00000000162716605 0.00000000003486784

추가 리소스

이항분포 소개
이항 분포의 형태 이해