二维排列问题

Ysu2008

Ysu2008

2008年四川汶川大地震发生后，有人在网上搜索到一篇2006年发表于正式刊物的论文，该论文提前预报了08年的大地震。这让有关方面感到压力巨大。最后，这篇论文遭到主流学术界广泛批判。

因为论文采用的方法是不太为人所知的“可公度法”，由已故院士翁文波于上世纪将可公度现象应用于天灾预报而开创此法。



20世纪以来川滇地区发生 6.7 级以上的地震年份总计 25 年，作者选取了其中部分年份，排列出上面这张图。
这些年份构成几组等公差数列的规则结构。根据图中规律，不难算出问号年份正是 2008 年 。

这张图成了很多人的批评对象。批评者的理由是，作者不能从地震成因的角度解释为什么会有这么一个规则结构。

我认为，能不能从地震成因的角度解释并不是问题的关键。问题的关键在于：出现这么一个规则结构的概率有多大？

如果在100年中，随机选取 16 个年份，出现规则结构的概率很大，那么我们可以说论文作者的研究是没有意义的。

100年中随机选取16年构成规则结构的概率究竟有多大呢？我们可以模拟估算。

从 [1910,2010] 连续自然数区间中，随机抽取 16 个年份作为一个样本，每个样本有 24024 种升序排列，其中一些排列有可能是规则的（构成等公差数列），只要有一种排列是规则的，我们就说该样本是规则的。

随机抽取 50 万个样本，规则样本仅发现 2602 个，占比 0.005204 ，比例的 99% 置信区间估计为 (0.0049 , 0.0055)。

这就是说，论文中的规则结构出现的概率不到 1% ；或者说，论文作者随意拼凑出这个规则结构的概率不到 1% 。作者的发现显然很有意义。

Ysu2008

随机抽取 50 万个样本，规则样本仅发现 2602 个。
——有可能是程序设计错误导致的错误结果。如果对该问题有兴趣，欢迎重复我的模拟实验。