《虹桥飞机场出租车大数据分析报告报告材料..doc》由会员分享，可在线阅读，更多相关《虹桥飞机场出租车大数据分析报告报告材料..doc（16页珍藏版）》请在三一文库上搜索。
1、实用标准文档虹桥飞机场出租车数据分析报告- 排队长度应山华东师大学金融与统计学院从数据我们可以得到五种数据：label yy1= '逗留时间 = 等待时间 + 服务时间 ' f10label yy2= '排队长度 ' ;fi 中非零数据个数label yy3= '服务时间 ' ;做差 fi-f(i-1)label yy4='服务流量 ' ;单位时间（一小时） fi 中非零数据个数label yy5='顾客流量 ' ;一天 fi 中非零数据个数下面对排队长度 数据进行分析。 服务流量、顾客流量 数据分析类似。1识别性
2、分析：识别原始数据对原始数据进行原始识别处理，先画柱状图形（直方图）和文案大全实用标准文档饼状图形如下文案大全实用标准文档从图形可以猜想，其图形是普松分布。普松分布的一阶矩估计为myyy=myy-minyy;lambda1=myyy ;3.28571普松分布的二阶矩估计为lambda2=ssyy*2 ; 6.37372文案大全实用标准文档普松分布的区间估计为（见茆诗松和周纪芗，概率论与数理统计，中国统计， 2007 ，p334-337 ）lambdaxx=mn*myyy;array kk1000 kk1-kk1000;do k= 1 to 1000 ;kkk=1-poisson(lambdax
3、x,k)-probchi(2 *lambdaxx,2*(k+ 1 );if k> 2 and kkk<0 then goto ok;end;ok:;lambdaxx=cinv(0.05, 2 *k)/(2 *mn);3.11526lambdal=cinv(0.025,2 *k)/(2*mn);3.07970lambdau=cinv(0.9755, (2*k+1 )/( 2*(mn); 3.54494lambdasx=cinv(0.95, 2 *k)/(2 *mn);3.50296lambdazwqj=cinv(0.5 , 2*k)/( 2*mn); 3.30532文案大全实用标准文
4、档2稳健性分析：比较各种估计的差别普松分布的一阶矩估计和普松分布的区间估计两种估计比较接近，但普松分布的二阶矩估计与前两者差别很大。用三种估计的平均值作为参数 lambda的估计应是比较稳健的lambdaqz=(lambda1+lambda2+lambdazwqj)/3;4.32159即使不考虑二阶矩估计，用其他两种估计的平均值作为参数lambd的估计也应是比较稳健的lambdaqz=(lambda1+lambdazwqj)/2 ; 3.295523 协调性分析：比较估计分布与经验离散频率分布3.2 ：诊断估计分布与经验离散频率分布图形的协调性以上述参数作为普松分布的估计与经验离散频率 分布函
5、数的比较图形为文案大全实用标准文档从图形看出拟合不好 .即使不考虑二阶矩估计，用其他两种估计的平均值作为参数lambd的估计 ,以此参数作为普松分布的估计与经验离散频率 分布函数的比较图形为文案大全实用标准文档从图形看出拟合也不好 . 尽管它看起来比前一种好一点，即图形的接近程度高一些 , 但两种分布的差别还是较大。3.2 ：估计分布与经验离散频率分布理论的协调性用皮尔逊的 Kf 统计量进行拟合检验两种情况都没有通过.N=238,l=1,Kf=728.110>272.836(显著性水平0.05)N=238,l=1,Kf=325.637>272.836(显著性水平文案大全实用标准文档
6、0.05)其中 n 为数据个数， l 为估计参数个数。而Kf 为经验离散频率 分布函数和拟合分布的差的平方的加权和，近视为KF 统计量，自由度为（ n-l-1 ），272.836 是统计量 Kf 显著性水平 0.05 临界值。4功能性分析：变换原始数据上述分布的各种参数为mnmyyc dssyymaxyy minyy pz yyp50 mumyyylambda1 lambda2 lambdaxx klambdal238 3.28571 1 0 2.5246280330 3.28571 3.28571 6.373723.11526 9 3.07970lambdau lambdasx lambda
7、zwqj lambdaqzxitazddd kkkpixitaLxitaUfun3.544943.502963.305324.321594.32159130.896788 1.321597.32159平时柏松分布其中普松分布的一阶矩估计和普松分布的二阶矩估计不一致, 这可能是估计不准确的原因。 为了消除这个原因， 我们采用数据变换的方法解决。对数据 yy 做变换：myyy=myy-minyy;文案大全实用标准文档c=(myyy)/(ssyy*2); 0.51551d=0;yy=int(myyy+d)*c);（此时变换后的均值为变换前的c 倍，而变换后的方差为变换前的c2倍,令两者相等，得到上述
8、的变换公式，由于取整数函数int(*)有低估现象，可以用d 来做小的调整）。用这些数据进行上述同样的分析得到图形如下文案大全实用标准文档文案大全实用标准文档从图形看出效果很好。用皮尔逊的 Kf 统计量进行拟合检验通过.N=238,l=3,Kf=13.8675<270.684( 显著性水平0.05)其中 n 为数据个数， l 为估计参数个数。而Kf 为经验离散频率 分布函数和拟合分布的差的平方的加权和，近视为KF 统计量，自由度为（n-l-1 ），270.684 是统计量 Kf 显著性水平 0.05 临界值。文案大全实用标准文档得到的分布参数为mnmyycdssyymaxyyminyypz
9、yyp50mumyyylambda12381.357140.5155101.19888401101.357141.35714lambda2lambdaxxklambdallambdaulambdasxlambdazwqjlambdaqzxitaz1.437311.2554241.233031.534151.505661.376751.390402.69714xitadddkkkpiLxitaUfunCLUCLLCL130.94735605.69714平时柏松分布2.632625.752770其中普松分布的一阶矩估计 (lambda1 =1.35714) 和普松分布的二阶矩估计 (lambda2
10、 = 1.43731 )基本一致 .这可能就是估计准确的原因。数据变换方法使得两个估计接近相等。5经济性分析：变换控制图根据变换后的估计分布的控制图， 变回到原来分布的控制图的控制限为UCL= 5.75277CL= 2.63262文案大全实用标准文档LCL=0要求每五个数据进行平均得到一个点 （低于五个数据平均时， 由于随机性太强，控制图波动性太大，看不出规律性） 。由此得到控制图为：从图中看到：排队长度分布具有高峰期发生， 高峰期发生的条件可以判定为队长大于五。用这个条件我们可以分析各种数据的高峰现象。对应的底峰期发生的条件也可以判定为队长不大于六 , 用这个条件我们也可以分析各种数据的低峰现象。文案大全实用标准文档文案大全