科学发现年龄定律的再探讨

周爱民

（郑州大学图书馆 河南 郑州 450002）

梁立明、 赵红洲已探讨过科学发现年龄定律。 他们探讨的科学发现年龄定律分析方法已作为经典方法进入了许多科学计量学书中。 为了探索出科学发现年龄符合威布尔分布，我们希望不作任何参数假定的情况下，得到威布尔分布的参数，为此我们从威布尔分布基本概念谈起。

1 威布尔分布基本概念

定义：若非负随机变量X 有分布函数

其中a、γ、k 均大于0，则称X（或F（x））有威布尔分布。设随机变量x 服从三参数威布尔分布，则其概率密度函数和累积概率分布函数为

其中，K＞0 为形状参数，a＞0 为尺度参数，γ≥0 为位置参数。

表1 1500-1960 年按世纪分期重大科学发现与年龄的统计数据

2 重大科学发现的年龄数据（表1）

3 数据处理

威布尔分布一般只给出时间， 即第一个数值产生的时间、第二个数值产生的时间、第三个数值产生的时间、…、第n个数值产生的时间。 但1500-1960 年按世纪分期重大科学发现与年龄的统计数据却是分组型数据，即每十岁间产生了多少重大科学发现，因此我们不能直接利用威布尔模型。 我们必须给出新的方法。 我们仅分析1601-1700 年期间的数据。由于此数据表给出的是岁数与发生的频率，而不是重大科学发现产生时的具体岁数，所以我们下面用威布尔累积概率分布函数来作为拟合模型。

威布尔模型中，F（x）表示累积概率。 显然我们可以把问题转化成：

F（20）-F（10）=0.05

F（30）-F（20）=0.28

F（40）-F（30）=0.26

F（50）-F（40）=0.24

F（60）-F（50）=0.09

F（70）-F（60）=0.08

为此我们定义一个新差函数

Φ（x）的值如表2。

表2 Φ（x）的值

显然这是一个一元非线性回归问题， 我们利用dps 软件中的一元非线性回归程序计算得：