我怎样才能得到一辆汽车在 R 中其余汽车的泊松模拟？

Question

我有 mtcars 数据集我想要的是使用泊松分布进行模拟基本思想是使用 rpois() 模拟 hp，并将其与其余汽车进行比较，以查看一辆车的 hp 比另一辆低多少次。这样我就可以确定一辆汽车的马力低于其他汽车的可能性。为了简单起见，我在现实中使用这个例子我正在使用高尔夫球手并根据他们的平均得分我想知道赢得比赛的概率是多少

我可以用他们的分数测试一辆车的马力比另一辆车小的概率

n1 = 93
n2 = 100
sum(rpois(1000, n1) < rpois(1000, n2))/1000

在这种情况下，我正在测试马自达 RX4 和 Datsun 710 我想要一种方法来找到这个概率但是使用整个样本，换句话说，多少次马自达 RX4 的泊松模拟马力低于样本中的其他汽车

Answer 1

这是使用 outer 和矢量化比较函数 g.

比较所有汽车与所有汽车的方法

f <- function(n1, n2){
  mean(rpois(1000, n1) < rpois(1000, n2))
}
g <- Vectorize(f, c("n1", "n2"))
res <- outer(mtcars$hp, mtcars$hp, g)
dimnames(res) <- list(row.names(mtcars), row.names(mtcars))

我怎样才能得到一辆汽车在 R 中其余汽车的泊松模拟？

How can I get poisson simulation of one car over the rest of cars in R?

simulation

r

poisson