我怎样才能得到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))