R でポアソン分布をプロットする方法
R のポアソン分布の確率質量関数をプロットするには、次の関数を使用できます。
- dpois(x, lambda) は確率質量関数を作成します
- Lot(x, y, type = ‘h’) は、プロットがヒストグラム (type=’h’) であることを指定して、確率質量関数をプロットします。
確率質量関数をプロットするには、単純にlambdaを指定します。 (イベントの発生率など) 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