Hoe een poisson-verdeling in r te plotten

Om de waarschijnlijkheidsmassafunctie voor een Poisson-verdeling in R uit te zetten, kunnen we de volgende functies gebruiken:

• dpois(x, lambda) om de waarschijnlijkheidsmassafunctie te creëren
• plot(x, y, type = ‚h‘) om de waarschijnlijkheidsmassafunctie te plotten, waarbij wordt gespecificeerd dat de plot een histogram is (type=’h‘)

Om de waarschijnlijkheidsmassafunctie uit te zetten, geeft u eenvoudigweg lambda op   (bijvoorbeeld het aantal gebeurtenissen) in de functie dpois() .

De volgende code illustreert bijvoorbeeld hoe u een waarschijnlijkheidsmassafunctie plot voor een Poisson-verdeling met lambda = 5:

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

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

De x-as toont het aantal ‘successen’ – bijvoorbeeld het aantal gebeurtenissen dat heeft plaatsgevonden – en de y-as toont de waarschijnlijkheid dat dat aantal successen in twintig pogingen wordt behaald.

We kunnen een titel toevoegen, de aslabels wijzigen en de breedte van de lijnen vergroten om de plot esthetisch aantrekkelijker te maken:

 success <- 0:20

plot(success, dois(success, lambda=5),
     type='h',
     main='Fish Distribution (lambda = 5)',
     ylab='Probability',
     xlab ='# Successes',
     lwd=3)

We kunnen de volgende code gebruiken om de werkelijke kansen te krijgen voor elk aantal successen dat in de grafiek wordt weergegeven:

 #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