发明名称 一种基于地理加权泊松回归的县级交通事故的预测方法
摘要 本发明公开了一种基于地理加权泊松回归的县级交通事故的预测方法;首先通过数据调查采集研究范围内各县的历史交通事故频次数据,然后设定县级事故频次和影响因素之间的模型关系和每个变量的影响系数的地理位置矩阵形式,以及进行参数估计过程中的地理权重函数,其后根据地理权重的泊松回归模型进行参数估计,克服了以往县级交通事故预测中没有考虑空间异质性而造成预测精准度不高的问题,为县级交通事故预测和安全防治政策指定提供科学、合理的指导。
申请公布号 CN103218668A 申请公布日期 2013.07.24
申请号 CN201310041718.0 申请日期 2013.02.01
申请人 东南大学 发明人 李志斌;刘攀;王炜;徐铖铖
分类号 G06Q10/04(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种基于地理加权泊松回归的县级交通事故预测方法,其特征是包括步骤:1)通过数据调查采集研究范围内各县的历史交通事故频次数据,以及各县域内人口统计学属性、社会经济属性、道路基础设施属性及交通量属性数据,提取出交通事故频次、交通量、道路属性、社会人口统计特征、车辆登记数量和持有驾照人数比例参数信息;2)采用泊松分布模型形式设定县级事故频次和影响因素之间的模型关系:ln(Y)=ln(β0(ui))+β1(ui)ln(DVMT)+β2(ui)X2+β3(ui)X3+…+βJ(ui)XJ+□    (1)其中,事故数Y为县级交通事故数量的预测值;βJ(ui)表示县城地理中心的坐标位置ui=(uxi,uyi);X1…XJ为第1J个自变量,□为随机误差变量;3)设置每个变量的影响系数的地理位置矩阵形式,并进行参数估计过程中的地理权重函数: <mrow> <mi>&beta;</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>&beta;</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>&beta;</mi> <mi>J</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>x</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>y</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&beta;</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>x</mi> <mn>2</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>y</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>x</mi> <mn>2</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>y</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>&beta;</mi> <mi>J</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mrow> <mi>x</mi> <mn>2</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>u</mi> <mrow> <mi>y</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&beta;</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>xn</mi> </msub> <mo>,</mo> <msub> <mi>u</mi> <mi>yn</mi> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>xn</mi> </msub> <mo>,</mo> <msub> <mi>u</mi> <mi>yn</mi> </msub> <mo>)</mo> </mrow> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>&beta;</mi> <mi>J</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>xn</mi> </msub> <mo>,</mo> <msub> <mi>u</mi> <mi>yn</mi> </msub> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>其中,n是县城数量; <mrow> <mover> <mi>&beta;</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>X</mi> <mi>T</mi> </msup> <mi>W</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>xi</mi> </msub> <mo>,</mo> <msub> <mi>u</mi> <mi>yi</mi> </msub> <mo>)</mo> </mrow> <mi>X</mi> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>X</mi> <mi>T</mi> </msup> <mi>W</mi> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>xi</mi> </msub> <mo>,</mo> <msub> <mi>u</mi> <mi>yi</mi> </msub> <mo>)</mo> </mrow> <mi>Y</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>W(uxi,uyi)通过n个空间加权矩阵表示一个n,通过W(i)表示为: <mrow> <mi>W</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>w</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>w</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> </mtd> <mtd> <msub> <mi>w</mi> <mi>in</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>其中wij(j=1,2,…,n)是i县城模型标定中给j县城的权重;运用高斯和毕方运算给出周围县城对其影响的权重:高斯: <mrow> <mi>wij</mi> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>&times;</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <mi>G</mi> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>毕方: <mrow> <mi>wij</mi> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mo>[</mo> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>/</mo> <msub> <mi>G</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msup> <mo>]</mo> <mn>2</mn> </msup> </mtd> <mtd> <mi>if</mi> <mo>|</mo> <mo>|</mo> <msub> <mi>u</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>&lt;</mo> <msub> <mi>G</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>otherwise</mi> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>其中,Gi为带宽变量,用于确定周围县城对目标县城产生影响的空间范围,Gi趋近于无穷大时,wij接近1,表明所有县城都对目标县城有影响;4)基于地理权重的泊松回归模型进行参数估计,得出每个变量系数随地理位置的变化值,以及各变量的显著性水平;5)将某年各县域内显著变量的属性值输入到模型中,对于各县级的交通事故频次进行预测;6)检验各县级事故频次预测值的残差是否有空间相关性,若无空间相关性,表明建立的基于地理权重的泊松回归模型可以用于有效预测县级交通事故数;如具有空间相关性,则应重新选择属性变量,从步骤4)开始重复,直至预测值的残差不具有空间相关性;所述预测值的残差的空间相关性检测方法为由莫兰开发的一种研究空间自相关性的方法,Moran’s I: <mrow> <msup> <mi>Moran</mi> <mo>&prime;</mo> </msup> <mi>sI</mi> <mo>=</mo> <mfrac> <mrow> <mi>N&Sigma;i&Sigma;j</mi> <msub> <mi>&omega;</mi> <mi>ij</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>Y</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>j</mi> </msub> <mo>-</mo> <mover> <mi>Y</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&Sigma;i&Sigma;j</mi> <msub> <mi>&omega;</mi> <mi>ij</mi> </msub> <mi>&Sigma;i</mi> <msup> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>Y</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>其中,N是县城数量;Yi和Yj是县城i和县城j残余的事故;Y是所有县城的平均残余,ωij是i县城模型标定中给j县城的权重,是空间权重矩阵的元素,呈现的空间相关性结构;一个正或负值的Moran’s I表示各县城的正或负的空间自相关性多;Moran’s I的范围是从‑1到+1;零值表示一个随机的空间格局;为了检验统计假设,将Moran’s I转化为Z‑值,该值大于1.96或小于‑1.96表示在统计上很显著的5%的置信水平。
地址 210096 江苏省南京市玄武区四牌楼2号