虹橋飛機場出租車大數(shù)據(jù)分析報告報告材料.

1、虹橋飛機場出租車數(shù)據(jù)分析報告-排隊長度張應(yīng)山華東師范大學(xué)金融與統(tǒng)計學(xué)院從數(shù)據(jù)我們可以得到五種數(shù)據(jù)：label yy仁逗留時間二等待時間+ 服務(wù)時間;f10label yy2=排隊長度;fi中非零數(shù)據(jù)個數(shù)label yy3=服務(wù)時間; 做差 fi-f（i-1）label yy4=服務(wù)流量;單位時間（一小時）內(nèi)fi中非零數(shù)據(jù)個數(shù)label yy5=顧客流量;一天內(nèi)fi中非零數(shù)據(jù)個數(shù)F面對排隊長度數(shù)據(jù)進行分析。服務(wù)流量、顧客流量 數(shù)據(jù)分析類似。1. 識別性分析：識別原始數(shù)據(jù)對原始數(shù)據(jù)進行原始識別處理，先畫柱狀圖形（直方圖）和餅狀圖形如下FREQUENCYFREQUENCY of yy11?S7303

2、D從圖形可以猜想，其圖形是普松分布 普松分布的一階矩估計為myyy=myy-minyy;Iambda1=myyy ;3.28571普松分布的二階矩估計為Iambda2=ssyy* 2; 6.37372普松分布的區(qū)間估計為（見茆詩松和周紀薌，概率論與數(shù)理統(tǒng)計，中國統(tǒng)計出版社，2007，p334-337）lambdaxx=mn *myyy;array kk 1000 kk1-kk1000;do k= 1 to 1000;kkk= 1-poisson(lambdaxx,k)-probchi( 2*lambdaxx, 2*(k+ 1);if k 2 and kkk272.836 （顯著性水平0.05）

3、N=238, l=1, Kf= 325.637 272.836 （顯著性水平0.05) 其中n為數(shù)據(jù)個數(shù)，I為估計參數(shù)個數(shù)。而Kf為經(jīng)驗離散頻率分布 函數(shù)和擬合分布的差的平方的加權(quán)和， 近視為KF統(tǒng)計量，自由度為(n-I-1 )，272.836是統(tǒng)計量Kf顯著性水平0.05臨界值。4. 功能性分析：變換原始數(shù)據(jù)上述分布的各種參數(shù)為mn myy c d ssyy maxyy minyy pz yyp50 mu myyy lambdal Iambda2 lambdaxx k lambdal238 3.28571 1 0 2.52462 803 30 3.28571 3.28571 6.37372

4、3.11526 9 3.07970lambdau lambdasx lambdazwqj lambdaqz xitaz ddd kkk p i xitaL xitaU fun3.54494 3.50296 3.30532 4.32159 4.32159 1 3 0.89678 8 1.32159 7.32159平時柏松分布其中普松分布的一階矩估計和普松分布的二階矩估計不一致，這可能是估計不準確的原因。為了消除這個原因，我們采用數(shù)據(jù)變換的方 法解決。對數(shù)據(jù)yy做變換：myyy=myy-minyy;c=(myyy)/(ssyy* 2);0.51551d=0;yy=i nt(myyy+d)*c);(

5、此時變換后的均值為變換前的 c倍，而變換后的方差為變換前的 cA2倍,令兩者相等，得到上述的變換公式，由于取整數(shù)函數(shù)int(*)有低估現(xiàn)象，可以用d來做小的調(diào)整)。用這些數(shù)據(jù)進行上述同樣的分析得到圖形如下貯 MIDPOINTFREQUENCY o yy2 . 0VV從圖形看出效果很好。用皮爾遜的 Kf 統(tǒng)計量進行擬合檢驗通過 .N=238, l=3, Kf= 13.8675 270.684 （ 顯著性水平 0.05） 其中 n 為數(shù)據(jù)個數(shù)， l 為估計參數(shù)個數(shù)。而 Kf 為經(jīng)驗 離散頻率 分布 函數(shù)和擬合分布的差的平方的加權(quán)和，近視為KF統(tǒng)計量，自由度為（ n-l-1 ）， 270.684 是

6、統(tǒng)計量 Kf 顯著性水平 0.05 臨界值。得到的分布參數(shù)為mn myy c d ssyy maxyy minyy pz yyp50 mu myyy lambda1238 1.35714 0.51551 0 1.19888 4 0 1 1 0 1.35714 1.35714lambda2 lambdaxx k lambdal lambdau lambdasx lambdazwqj lambdaqz xitaz1.43731 1.25542 4 1.23303 1.53415 1.50566 1.37675 1.39040 2.69714xitaddd kkk p i L xitaU fun CL UCL LCL1 3 0.94735 6 0 5.69714平時柏松分布 2.63262 5.75277 0其中普松分布的一階矩估計 (lambda1 =1.35714) 和普松分布的二階矩估計 (lambda2 = 1.43731 ) 基本一致 . 這可能就是估計準確的原因。 數(shù)據(jù)變換方法使得兩個估計接近相等。5. 經(jīng)濟性分析：變換控制圖根據(jù)變換后的估計分布的控制圖，變回到原來分布的控制圖的控制限 為UCL二 5.75277CL= 2.63262LCL=0要求每五個數(shù)據(jù)進行平均得到一個點 （低于五個數(shù)據(jù)平均時，由于隨 機性太強，控制圖波動性太大，看不出規(guī)律性）。由此得到控制圖為：從圖中