虹橋飛機(jī)場(chǎng)出租車大數(shù)據(jù)分析報(bào)告報(bào)告材料.

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

2、D從圖形可以猜想，其圖形是普松分布 普松分布的一階矩估計(jì)為myyy=myy-minyy;Iambda1=myyy ;3.28571普松分布的二階矩估計(jì)為Iambda2=ssyy* 2; 6.37372普松分布的區(qū)間估計(jì)為（見茆詩(shī)松和周紀(jì)薌，概率論與數(shù)理統(tǒng)計(jì)，中國(guó)統(tǒng)計(jì)出版社，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ù)個(gè)數(shù)，I為估計(jì)參數(shù)個(gè)數(shù)。而Kf為經(jīng)驗(yàn)離散頻率分布 函數(shù)和擬合分布的差的平方的加權(quán)和， 近視為KF統(tǒng)計(jì)量，自由度為(n-I-1 )，272.836是統(tǒng)計(jì)量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ì)和普松分布的二階矩估計(jì)不一致，這可能是估計(jì)不準(zhǔn)確的原因。為了消除這個(gè)原因，我們采用數(shù)據(jù)變換的方 法解決。對(duì)數(shù)據(jù)yy做變換：myyy=myy-minyy;c=(myyy)/(ssyy* 2);0.51551d=0;yy=i nt(myyy+d)*c);(

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

6、統(tǒng)計(jì)量 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平時(shí)柏松分布 2.63262 5.75277 0其中普松分布的一階矩估計(jì) (lambda1 =1.35714) 和普松分布的二階矩估計(jì) (lambda2 = 1.43731 ) 基本一致 . 這可能就是估計(jì)準(zhǔn)確的原因。 數(shù)據(jù)變換方法使得兩個(gè)估計(jì)接近相等。5. 經(jīng)濟(jì)性分析：變換控制圖根據(jù)變換后的估計(jì)分布的控制圖，變回到原來(lái)分布的控制圖的控制限 為UCL二 5.75277CL= 2.63262LCL=0要求每五個(gè)數(shù)據(jù)進(jìn)行平均得到一個(gè)點(diǎn) （低于五個(gè)數(shù)據(jù)平均時(shí)，由于隨 機(jī)性太強(qiáng)，控制圖波動(dòng)性太大，看不出規(guī)律性）。由此得到控制圖為：從圖中