[发明专利]一种应用于区域路网的交通事故严重度预测方法

权利要求书

1.一种应用于区域路网的交通事故严重度预测方法，其特征是按如下步骤进行：

步骤一、区域路网道路交通事故数据的采集与预处理；

从道路交通事故数据库中获取N起事故数据作为事故数据集D，并从任意第i起事故数据中选取K个分类变量组成集合X＝{x₁,x₂,…,x_k,…,x_K}来表征第i起事故，其中，x_k表示第k个分类变量，且第k个分类变量x_k包含C_k种类别，第k个分类变量x_k在C_k种类别中的取值记为s_k，令s_ik表示第i起事故的第k个分类变量的取值，则第i起事故中所有K个分类变量的取值所组成的分类变量取值集合记为S_i＝{s_i1,s_i2,...,s_ik,...,s_iK}；令表示第i起事故的K个分类变量的所有可能取值中的任意一种取值集合；k＝1,2,3,...,K；i＝1,2,3,...,N；

将第i起事故的严重度作为预测变量，记为y_i，且y_i的取值为“0”或“1”分别表示非死亡事故和死亡事故；

步骤二、根据区域路网道路交通事故数据，建立潜在类别分析模型；

步骤2.1、定义所述潜在类别分析模型中存在一个潜在类别变量V，V包含T种类别，且任意一种类别记为t，t＝1,2,...,T；令第i起事故中潜在类别变量V的取值记为V_i；

步骤2.1.1、定义外循环次数为τ、最大外循环迭代次数为τ_max；令第τ次所设置的类别数目为T_τ；初始化τ＝1；

步骤2.1.2、初始化t＝1；

步骤2.1.3、初利用式(1)得到第i起事故V_i取值为t，即属于第t种潜在类别时，第i起事故在K个分类变量上的取值集合为的条件概率

式(1)中，P(s_ik＝s_k|V_i＝t)表示第i起事故属于第t个潜在类别时，第k个分类变量上取值为s_k的条件概率；

步骤2.1.4、利用式(2)得到第i起事故中K个分类变量取值集合为的非条件概率即潜在类别分析模型的联合概率

式(2)中，P(V_i＝t)是第i起事故属于第t个潜在类别的概率，潜在类别t占总体的比率；

步骤2.2、采用极大似然法进行模型参数估计，得到潜在类别概率和分类变量条件概率的估计值以及潜在类别分析模型的第τ次极大似然函数值L_τ；

步骤2.3、利用式(3)计算第i起事故被分类到第t个潜在类别的后验概率

步骤2.4、令t+1赋值给t，并判断t＞T_τ是否成立，若成立，则执行步骤2.5；否则，返回步骤2.1.3执行；

步骤2.5、利用式(4)、式(5)、式(6)和式(7)得到模型拟合评价指标，包括：第τ次信息评价指标AIC_τ、第τ次贝叶斯信息准则BIC_τ、第τ次样本校正的贝叶斯信息准则aBIC_τ、第τ次熵值

AIC_τ＝-2ln(L_τ)+2M (4)

BIC_τ＝-2ln(L_τ)+ln(N)×M (5)

aBIC_τ＝-2ln(L_τ)+ln(n^*)×M (6)

式(4)、式(5)、式(6)和式(7)中，M为潜在类别分析模型中未知参数的个数；n^*是调整后的样本量，且n^*＝(N+2)/24；

步骤2.6、将τ+1赋值给后τ，判断τ＞τ_max是否成立，若成立，则执行步骤2.7；否则，返回步骤2.1.3执行；

步骤2.7、从τ_max次信息评价指标AIC、贝叶斯信息准则BIC、样本校正的贝叶斯信息准则aBIC和熵值R²中选出各个模型拟合评价指标均取到最优值时所对应的潜在类别个数，记为T^*；将所述事故数据集D划分为T^*个事故子类别，记为表示第t^*个事故子类别中的事故数据，t^*＝1,2,…,T^*；

步骤三、根据潜在类别分析模型结果，对T^*个事故子类别分别建立CART决策树模型；

步骤3.1、令所述第t^*个事故子类别中的事故数据作为训练样本集，令K个分类变量所组成的集合X为所述CART决策树模型中的特征集；令结点样本阈值为σ、特征值切分点为α、Gini指数阈值为ε；

步骤3.2、初始化t^*＝1；

步骤3.3、将所述训练样本集特征集X、定义结点样本阈值σ和Gini指数阈值ε输入所述CART决策树模型；

步骤3.4、令t^*+1赋值给t^*，并判断t^*＞T^*是否成立，若成立，则表示得到T^*个二叉决策树，并执行步骤3.5；否则，返回步骤3.3执行；

步骤3.5、根据所述T^*个二叉决策树的树形图，确定分类变量间的交互作用项，其中，第t^*个事故子类别对应的二叉决策树所确定的交互作用项；

步骤四、对T^*个事故子类别分别建立基于二元logistic回归的事故严重度模型；

步骤4.1、将所述第t^*个事故子类别中的事故数据作为事故严重度模型的拟合数据，以K个分类变量所组成集合X和第t^*个事故子类别的交互作用项共同作为所述事故严重度模型的自变量X^*；定义第t^*个事故子类别包含J个事故数据，J的值为第j起事故的预测变量记为y_j；

步骤4.2、初始化t^*＝1；

步骤4.3、利用式(11)得到基于二元logistic回归在自变量X^*条件下死亡事故即y_j＝1的发生概率P(y＝1|X^*)：

式(11)中，w^*为自变量X^*的回归系数；

步骤4.4、利用极大似然法估计所述二元logistic回归的事故严重度模型的参数w^*：

对于第j起事故，为给定自变量条件下y_j＝1的概率，则给定自变量条件下y_j＝0的概率为1-P_j；并利用式(12)得到似然函数L(w^*)：

利用极大似然估计，求出使得L(w^*)取得最大值时的估计参数w′；

根据估计参数w′得到第j起事故在自变量条件下y_j＝1的预测概率从而得到J起事故的预测概率并进行升序排序，得到排序后的预测概率集合记为{P′₁,...,P′_j,...,P′_J}；

步骤4.5、调整事故严重度模型的预测分类阈值；

步骤4.6、令t^*+1赋值给t^*，并判断t^*＞T^*是否成立，若成立，则表示获得T^*个事故严重度预测模型，否则，返回步骤4.3执行。