R中的本福德定律

Question

我正在嘗試在所有供應商的發票中使用 R 中的 benford.analysis 包來實施本福德定律。 在整個數據集中，數據證實了這一點。 我正在嘗試找到一種按供應商分組的方法，以確定是否有個別供應商通過不符合標准來顯示欺詐指標。 有沒有辦法按組打破不合格？

Answer 1

這是一種使用group_by和group_map為每個組創建benford.analysis圖的方法。 在此示例中，按物種對鳶尾花數據進行分組並對萼片長度變量進行分析。

在group_map()中， .x表示分組的子集數據， .y表示組的名稱。

library(dplyr)
library(benford.analysis)
iris %>% 
  group_by(Species) %>% 
  group_map(.f = ~ plot(benford(.x$Sepal.Length)))