如何在 r 中绘制二项式分布
要绘制 R 中二项式分布的概率质量函数,我们可以使用以下函数:
- dbinom(x, size, prob)创建概率质量函数
- plot(x, y, type = ‘h’)绘制概率质量函数,指定该图是直方图 (type=’h’)
要绘制概率质量函数,只需指定大小 dbinom()函数中的(例如试验次数)和prob (例如给定试验的成功概率)。
例如,以下代码说明了如何绘制大小 = 20 且概率 = 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