对恐怖袭击事件的多元统计分析

谢振东　帆　周子雅　欧辉摘 要 恐怖击威胁人类社会安全选取全球恐怖主义数据库（Glbl rrr bG）0507年世界上发生恐怖击事件记录根据相关性对0506年知作案组织或人对应事件数据进行整合简化使用二阶聚类得到分类数根据危害性从到选出前5犯罪嫌疑人利用判别分析对07年知作案组织或人对应事件进行概率预测得出嫌疑人嫌疑程选取影响恐怖击重要指标应用因子分析研究0507年恐怖击事件发生规律得到恐怖事件地域发展趋势关键词社会统计学;趋势预测;多元分析;恐怖击图分类 8献标识码 br rrr k r r l r r x l Glbl rrr b （Glbl rrr b G） 05～07 rrr k rl rr v r 05 06 k r rgz r vl rrg v gr b l rr lrg r g x br l rg r r bg ll v g r l r 07 k r rgz r vl rrg v rbbl r l " gr ; r r r g rrr k r l l r l rr l rrr k 05～07 rgl vl r rrr bK r l ; r r; lvr l; rrr k 引 言恐怖主义是人类共威胁对恐怖击事件相关数据深入分析有助加深人们对恐怖主义认识反恐防恐提供有价值信息支持（06）[]探了0多年全球恐怖主义地理和上集性与持续性说明恐怖击事件非随机具有可研究性（06）[]提出了对欧洲圣战主义十年预测0年圣战分子欧洲发动击频率可能会呈现波动曲线峰值将逐渐升高李益斌（08）[3]利用聚类分析法以G恐怖击频率和烈依据对0000年和006年两阶段欧洲恐怖击按国进行分类根据实验结异分析当前欧洲恐怖击新态势杨振柳等（09） []对G进行数据挖掘利用主成分分析方法建立了K算法聚类模型通该聚类模型将相似特征事件进行归纳处理侦破恐怖击事件提供信息支撑华雅伦和王奇（08） [5]通使用G数据对0005年欧洲恐怖主義活动进行统计分析发现呈现出以热地区心扩散辐射圈恐怖分子击日趋专业化和技术化特总结了欧洲国反恐验彭如香等（09） [6]以0006年G数据对全球恐怖主义活动进行定量分析与可视化分析得出近些年全球恐怖主义活动发展态势结合近几年全球恐怖事件特总结出全球恐怖主义活动特征和发展趋势是有价值数据G或公开发表及出版献其缺失部分采用期望值化方法即算法估计所得采用二阶聚类方法可对含有类别变量和连续变量型数据集进行分析并动确定终分类数采用贝叶斯判别法不仅能够考虑到各总体出现先验概率又能考虑到错判造成损失且其判别效相比其他几种判别方法更加理想 分析方法 统计分析步骤对0506年发生、尚有组织或人宣称责恐怖击事件进行分析出主导某恐怖击事件有可能嫌疑人步将相关事件根据相关性进行整合二步利用二阶聚类模型对0506年处理数据进行二阶聚类分析将数据进行分类并按照危害性从到选出前五犯罪嫌疑人三步依据判别分析法贝叶斯方法对数据进行建模与分析得出嫌疑人嫌疑程研究近三年恐怖击事件发生规律选取恐怖击发生城市心次数、疑似恐怖击次数、恐怖击成功次数、恐怖击造成死亡总人数、亡人事件分比、平每起事件死亡人数、恐怖击造成受伤总人数、人质遭受绑架次数等8指标以区域分类变量进行分析研究下年全球或某些重地区恐怖击态势恐怖主义指数受许多因素影响依据世界济与和平研究所发布56国恐怖主义风险指数报告选取了其总起数、总亡数、总受伤人数和财产损失等四指标通计算四指标分值再合并得到综合分值考虑到目前财产损失价格数值估计并不十分准确且获取较困难而心城市是国济政治心尤重要所以加入“恐怖击发生城市心次数”这指标人质绑架会给人們带巨心理压力造成社会定程动荡不安对国面对恐怖击风险既要考虑直接恐怖击风险还要考虑国面对恐怖击脆弱性脆弱性低味着类别恐怖击可能会造成更人员伤亡和财产损失所以加入“亡人事件分比、平每起事件死亡人数这两指标” 统计模型二阶聚类模型步准聚类程（rlr ）采用专门应用巨型数据BR（Bl rv Rg lrg g rr）聚类算法建构多水平结构聚类特征树（lr r rr）其特它是种高平衡树（gbl r）—般用Vr建立r3 依据事件特征发现恐怖击事件制造者对0506年发生、尚有组织或人宣称责恐怖击事件运用二阶聚类方法将可能是恐怖组织或人不、不地多次作案若干案件归类并选出前五犯罪嫌疑人依据贝叶斯判别分析方法对07年所发生事件犯罪嫌疑人可能性进行归类判别07年主导某恐怖击事件有可能嫌疑人由数据复杂性以及海量性先将相关事件进行整合整合规则由下列指标定v、x、rg、v、brr、、、k、rg 、、kll、、rr、k、_接着利用二阶聚类模型对0506年处理数据进行二阶聚类分析结如图所示考虑到恐怖组织与人种类复杂性与多样性满足分类数据差异性前提下将0506年处理数据分50类基贝叶斯判别分析模型进行判断组织或人嫌疑程根据指标变量提取特征变量以某组织或人对做了该恐怖击事件概率作衡量标准列出分类函数系数如表所示图可清晰看出危险排序前五位以及其组质心可以通概率判别07年各事件是这五危险高人或组织所可能性例如选取07009003、0700037表即事件07009003嫌疑人所概率是07900四嫌疑人所概率是009990700037嫌疑人所概率是00000二嫌疑人所概率是099997四嫌疑人所概率是00000 对恐怖击态势分析对近三年即0507年恐怖击事件分别分析并进行横向和纵向对比出规律研究下年全球或某些重地区恐怖击态势K和巴特利检验（表3）显示K取样适切性量数06305所以该数据适合做因子分析根据表采用主成分分析集了8原始变量信息9676%以各因子方差贡献率占两因子总方差贡献率比重作权重进行加权且汇总基因子分析模型得出各城市综合得分其表达式（6765×+856×）÷9676（0）由表可知公共因子“击城市心次数”、“疑似恐怖击次数”、“击成功次数”、“死亡总人数”、“受伤总人数”、“人质遭受绑架次数”上荷值都很主要反映整体遭受恐怖击风险命名基风险因子而公共因子“亡人事件分比”、“平每起事件死亡人数”荷值较高主要反映应对恐怖击能力命名脆弱性因子得出0507年（基风险因子）、（脆弱性因子）两主要因子对比如表5和表6所示基风险因子得分反映该地区遭受恐怖击风险程纵向对比可知区域0（东和北非）近三年基风险因子得分高是要加强防恐怖击事件地区区域（东亚）风险性低安全性高横向比较可知区域0（东和北非）尽管三年期风险性有所下降但依然位地区首位所以仍是反恐重地区区域6（南亚）风险性还提高所以该地区下年反恐重地区脆弱性因子得分反映该地区应对恐怖击事件能力纵向对比可知区域（撒哈拉以南非洲）近三年脆弱性因子得分高是要提高其应对能力地区相对说区域（澳利亚和洋洲）应对能力强横向比较可知区域（撒哈拉以南非洲）三年应对处理恐怖击能力还下降值得提是区域（东亚）应对能力不断提高结合基风险因子和脆弱性因子得到综合得分由图3可知区域0（东和北非）及区域6（南亚）恐怖击事态较严重不仅遭受恐怖击风险性高而且其应对处理恐怖击能力也较低另外区域5（东南亚）综合得分不断提高可知恐怖击事件有向东南亚区域移趋势5 结 论通灵活应用二阶聚类、判别分析、因子分析对恐怖击事件记录数据进行量化分析考虑到事件相关性对0506年知作案组织或人对应事件数据進行整合简化应用二阶聚类充分利用数据差异性选取分类数50类因恐怖组织与人种类繁多且具有多样性这50类每类所包含事件数多少从到选出其前5对恐怖事件嫌疑程基贝叶斯判别分析模型进行判断组织或人嫌疑程根据指标变量提取主要特征计算出某组织或人实施该恐怖击事件概率根据其作衡量标准根据主要原因、空特性、蔓延特性、级别分布等规律进而分析总结近三年恐怖击事件发生规律和并预测下年全球某些重地区恐怖击态势选取了8指标以区域分类变量进行因子分析得到0507年关地区两公共因子上得分和总得分横向和纵向对比可知区域0（东和北非）恐怖击事态较严重不仅体现遭受恐怖击风险性高而且其应对处理恐怖击能力也较低下年反恐重地区区域0（东和北非）与区域（撒哈拉以南非洲）也是要加强反恐地区另外根据区域5（东南亚）综合得分不断提高可知恐怖击事件有向东南亚区域移趋势总体说近三年恐怖事件发生区域及其风险性较稳定参考献[] K r r rrr xlr r glbl r rrr v vl vr 0 r[] l rrr0639（）0309[] r r v [] rv rrr060（6）5670[3] 李益斌欧洲恐怖主义新态势及原因分析——基聚类分析法[]情报杂志0837（3）5563[] 杨振柳李颖钟子森基数据挖掘恐怖击嫌疑对象判断研究[]情报探09（）350[5] 华雅伦王奇基G数据库欧洲反恐形势分析及对我国启示[]犯罪研究08（5）905[6] 彭如香张奥博杨涛孔华锋基G全球恐怖主义活动现状与发展趋势研究[]计算机应用与软件0936（）5+ 相关热词 统计分析恐怖击事件