<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)

太感谢了

建议值的部分 不是很完善
就叫 近似吧 额呵呵

这边书看着不错, 请问叫啥名字

记得作分布拟合检验适当调整数据分组