R에서 포아송 분포를 그리는 방법

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

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

확률질량함수를 그래프로 그리려면 간단히 람다를 지정하면 됩니다.   (예를 들어 이벤트 발생률) dpois() 함수에서.

예를 들어, 다음 코드는 람다 = 5인 포아송 분포에 대한 확률 질량 함수를 그리는 방법을 보여줍니다.

 #define range of "successes"
success <- 0:20

#create plot of probability mass function
plot(success, dois(success, lambda=5), type='h')

x축은 “성공” 횟수(예: 발생한 이벤트 수)를 표시하고 y축은 20번의 시도에서 해당 성공 횟수를 달성할 확률을 표시합니다.

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

 success <- 0:20

plot(success, dois(success, lambda=5),
     type='h',
     main='Fish Distribution (lambda = 5)',
     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
dpois(success, lambda=5)

[1] 0.0067379469991 0.0336897349954 0.0842243374886 0.1403738958143
[5] 0.1754673697679 0.1754673697679 0.1462228081399 0.1044448629571
[9] 0.0652780393482 0.0362655774156 0.0181327887078 0.0082421766854
[13] 0.0034342402856 0.0013208616483 0.0004717363030 0.0001572454343
[17] 0.0000491391982 0.0000144527054 0.0000040146404 0.0000010564843
[21] 0.0000002641211