爱吱声

标题: 本福特定律验证疫情数据 [打印本页]
作者: togo    时间: 2020-2-13 10:00
标题: 本福特定律验证疫情数据
YouTube有一期李永乐老师教大家用本福特定律验证淘宝双十一数据是否作假的视频，我看了以后有点心得。
  n: Q! m1 M: O* n觉得用来验证一下现在的疫情数据可能有些参考意义# n& s; P$ m9 e, O

1 E' R9 i" c' D9 P0 ~
! Z! ~) V  _4 d3 C5 W$ V什么是本福特定律
* n! v. x, V' G% B" n9 e**********************- ^/ b8 A; T. B5 [# ~1 M3 k
1935年的某一天，物理学家富兰克.本福特（Frank.Benford）在图书馆查阅资料，他在翻阅对数表时，发现对数表的头几页要比后面的页脏一些。$ y0 E" m3 i8 f  L4 X( c" p
话说聪明的脑袋是一样的，愚笨的脑袋各式各样。牛顿的脑袋被苹果砸中，于是发现了万有引力。本福特也是如此，他拍了一下脑袋，发现了“本福特定律”。$ a, O3 H, m# N" V" h) Z
对数表的前几页比后面的脏，这说明有更多的人查阅头几页，这说明以1、2、3开头的数据比7、8、9开头的数据多。
( _8 v3 b. V0 r7 v% R6 k, d本福特搜集了人口、地理、经济等许多统计数据进一步分析，发现自然数据源，只要样本足够多，数据中以1为打头的数字出现的频率并不是1/9，而是30.1%。以2开头的数字出现的频率是17.6%。往后出现频率依次减少，以9为首的数字出现的频率最低，只有4.6%。
7 o2 S# k- W0 j( \$ W
数字
出现频率
# o+ U) U8 \8 p0 f
1
0.301

0 j& r5 g- ?+ y5 Y2 W
2
0.176

& n/ _7 S+ C0 m/ Z" X) o
3
0.125
0 i/ U. T: F5 k2 F" o: Z
4
0.097
- X  H5 D0 `1 U; k
5
0.079

8 H& U' k1 ?) Q6 E
6
0.067
/ t5 k0 @9 W8 V
7
0.058

6 K+ ?) z7 V6 G2 z
8
0.051
' Z3 {+ e3 j) K# l1 P
9
0.046
8 ], M; F4 z$ q$ n
合计
1.000
本福特定律的应用条件是：
) a  D$ {- C( b( v6 v/ C
1.数据不能是规律排序的。( u! A: z+ P: D) E0 @( j
2.数据不能经过人为修饰。% M& p6 E/ S# e( k5 J

6 R2 W7 }5 j& x' [9 U% T% v, k0 a************************************: ?* ^, ?1 \- X8 ?2 Z

, q( Q: L2 j% E( |4 B2 g" c. @8 K" j, R! }8 U) \) D& c$ Z9 P

* s4 e8 l1 n& E% ?2 d我从网上找到了一个疫情数据的统计表，刨去了比例计算的列，统计了其他数据里面首位出现的次数，粗粗比较一下两个数据的柱状图，感觉拟合是不错的。结论是疫情数据看不出有造假的痕迹。
  H+ U, ?8 u% r2 I- O4 s7 g! y3 G4 u$ z2 |  H2 f; O
9 X/ C& f9 C9 v- e7 ]0 Z

7 D" |9 U# \3 O+ L原始数据表. M; d$ k' ^- ?& q
8 f2 i2 K+ k: w6 v0 H7 Y4 O
作者: 龙血树    时间: 2020-2-13 10:29
拟合如何应该有检验指标以及参照的critical value？
作者: togo    时间: 2020-2-13 11:51
龙血树 发表于 2020-2-12 18:29
4 s2 V  H  m: _# g( m9 U+ Z% J' m拟合如何应该有检验指标以及参照的critical value？
" [2 W# r7 M( E; |4 d+ z+ A$ z% N: l
等明天上班后算一下correlation coefficient 和 p



欢迎光临 爱吱声 (http://www.aswetalk.net/bbs/) Powered by Discuz! X3.2