怎样用R语言求泊松分布参数

Python014

怎样用R语言求泊松分布参数,第1张

本帖最后由 jmpamao 于 2012-12-16 15:27 编辑

<9 记为 a

>17 记为 b

data <- read.table(text="

x Freq

9 282

10 547

11 704

12 682

13 527

14 340

15 188

16 91

17 39

",header=T)

pois.lambda <- function(a.min=0,a.max=8,a.Freq=73,b.min=18,b.max=50,b.Freq=23,data=data,p=0.95){

lambda <-c((sum(data[,1]*data[,2])+a.min*a.Freq+b.min*b.Freq)/(sum(data[,2])+a.Freq+b.Freq),

(sum(data[,1]*data[,2])+a.max*a.Freq+b.max*b.Freq)/(sum(data[,2])+a.Freq+b.Freq))

cat("lambda的取值范围",lambda,"\n")

lambda=round((lambda[1]*a.Freq+lambda[2]*b.Freq)/(a.Freq+b.Freq),digits=0)

pois=qpois(p,lambda)

cat("建议取值为",lambda,"\n",p,"的保证接通时,接话员数量",pois,"\n")

}

pois.lambda(0,8,73,18,50,23,data,p=0.95)

R语言泊松分布产生随机数函数rpois(n, lambda)

例子rpois(100,2) 按lambda=2的泊松分布产生了100个随机数 。

没有遇到过缺失值,你是不是哪里弄错了!

分析——非参数检验——K-S检验,在出现的对话框里,选中复选框组泊松分布即可。

泊松分布(Poisson distribution),台译卜瓦松分布,是一种统计与概率学里常见到的离散机率分布(discrete probability distribution)。泊松分布是以18~19 世纪的法国数学家西莫恩·德尼·泊松(Siméon-Denis Poisson)命名的,他在1838年时发表。但是这个分布却在更早些时候由贝努里家族的一个人描述过。就像当代科学史专家斯蒂芬·施蒂格勒(Stephen Stigler)所说的误称定律(the Law of Misonomy),数学中根本没有以其发明者命名的东西。