如何在 R 中对其他汽车进行泊松模拟？答案

【问题标题】：How can I get poisson simulation of one car over the rest of cars in R?如何在 R 中对其他汽车进行泊松模拟？
【发布时间】：2021-06-08 14:40:07
【问题描述】：

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

我可以通过使用他们的分数来测试一辆汽车的马力低于另一辆的概率

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

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

【问题讨论】：

泊松分布似乎不适合这种情况。这个答案建议您应该研究更有意义的贝叶斯推理/网络：stackoverflow.com/questions/26769518/…

标签： r simulation poisson

【解决方案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))

【讨论】：

第一行 (Mazda Rx4) 的第三列 (datsun710) 给出了 Datsun 的 hp 低于 Mazda rx4 的概率？
@PastorSoto 是的，就是这样。
我正在检查，似乎情况正好相反。与列相比，一辆车的 mpg 较低的概率在行中。在这个例子中，我认为应该是马自达RX4的马力低于达特桑710的概率。您可以检查并告诉我。另外，使用这个代码我能找到一辆车的马力比其他车最低的概率吗？