如何在箱中找到预期值

Question

我用随机变量生成了200个随机数，Y =指数（0.8）使用：

n=rexp(200,0.8)
bins=cut(n,breaks=c(0,1,2,3,max(n)))

我正在尝试计算每个bin的预期计数器数量，知道Y~Exponential（0.8）。我怎么能这样做？

Answer 1

您可以使用?rexp所用的同一组函数来计算每个范围内的比例，并使用它来计算预期数字：

num  <- 200
brks=c(0,1,2,3,Inf)

set.seed(12345)
n    <- rexp(num,0.8)
bins <- cut(n,breaks=brks)

expect.prop <- diff(pexp(brks, rate=0.8))
#[1] 0.55067104 0.24743245 0.11117856 0.09071795

rbind(
  observed=table(bins),
  expected=num*expect.prop
)
#            (0,1]    (1,2]    (2,3]  (3,Inf]
#observed 115.0000 42.00000 24.00000 19.00000
#expected 110.1342 49.48649 22.23571 18.14359