人口普查資料的一致性
在使用 Benford 定律之前,先確認你的資料集是否符合這個定律非常重要。資料集 census.2009 包含美國 19,509 個城鎮與城市的人口數(2009 年 7 月),此資料曾用於 Nigrini, M. J. (2012). 《Benford's Law: Application for Forensic Accounting, Auditing and Fraud Detection》, Wiley and Sons: New Jersey。現在來測試 2009 年的人口數是否滿足 Benford 定律。
繪圖會將圖形與圖例分開顯示,點選「Previous Plot」即可查看各張圖。
本練習屬於課程
R 的詐欺偵測
練習說明
- 載入
benford.analysis套件。 - 檢查變數
pop.2009是否符合第一位數的 Benford 定律分佈。注意,我們只看產生的第一張圖。其他圖超出本課程範圍。 - 將變數
pop.2009乘以 3,並再次檢查是否符合。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Load package benford.analysis
library(___)
data(census.2009)
# Check conformity
bfd.cen <- benford(___$___, number.of.digits = ___)
plot(bfd.cen, except = c("second order", "summation", "mantissa", "chi squared","abs diff", "ex summation", "Legend"), multiple = F)
# Multiply the data by 3 and check conformity again
data <- census.2009$pop.2009 * ___
bfd.cen3 <- benford(___, number.of.digits=___)
plot(bfd.cen3, except = c("second order", "summation", "mantissa", "chi squared","abs diff", "ex summation", "Legend"), multiple = F)