R語言中提供了四類有關統計分布的函數(密度函數,累計分布函數,分位函數,隨機數函數)。分別在代表該分布的R函數前加上相應前綴獲得(d,p,q,r)。如:
1)正態分布的函數是norm,命令dnorm(0)就可以獲得正態分布的密度函數在0處的值(0.3989)(默認為標准正態分布)。
2)同理,pnorm(0)是0.5就是正態分布的累計密度函數在0處的值。
3)而qnorm(0.5)則得到的是0,即標准正態分布在0.5處的分位數是0(在來個比較常用的:qnorm(0.975)就是那個估計中經常用到的1.96了)。
4)最后一個rnorm(n)則是按正態分布隨機產生n個數據。
上面正態分布的參數平均值和方差都是默認的0和1,你可以通過在函數里顯示指定這些參數對其進行更改。如dnorm(0,1,2)則得出的是均值為1,標准差為2的正態分布在0處的概率值。
要注意的是()內的順序不能顛倒。
關於二項分布的有關函數為:
The Binomial Distribution Description
Density, distribution function, quantile function and random generation for the binomial distribution with parameters size and prob.
dbinom(x, size, prob, log = FALSE)# 可用於計算二項分布的概率。
pbinom(q, size, prob, lower.tail = TRUE, log.p = FALSE)
qbinom(p, size, prob, lower.tail = TRUE, log.p = FALSE)
rbinom(n, size, prob)
例子: 已知某批雞蛋的孵出率prob為0.9,抽取size為5個雞蛋檢查其孵化情況,發現最終x=3個雞蛋孵化,求二項分布的概率。
R中的求解如下:
dbinom(3,5,0.9,log=F)
[1] 0.0729
#結果得到0.0729。