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

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

2、QUENCY60-FREQUENCYofyy從圖形可以猜想，其圖形是普松分布。普松分布的一階矩估計(jì)為myyy=myy-minyy;|lambda1=myyy;3.28571普松分布的二階矩估計(jì)為lambda2=ssyy*2;6.37372普松分布的區(qū)間估計(jì)為(見其詩松和周紀(jì)英，概率論與數(shù)理統(tǒng)計(jì)，中國統(tǒng)計(jì)出版社，2007,p334-337)lambdaxx=mn*myyy;arraykk1000kk1-kk1000;dok=1to1000;kkk=1-poisson(lambdaxx,k)-probchi(2*lambdaxx,2*(k+1);ifk>2andkkk<0thengot

3、ook;|end;ok:;lambdaxx=cinv(0.05,2*k)/(2*mn);3.11526lambdal=cinv(lambdau=cinv(0.025,0.97552*k)/(,(2*k+2*mn);1)/(3.079702*(mn);3.54494lambdasx=cinv(0.95,2*k)/(2*mn);3.50296lambdazwqj=cinv(0.5,2*k)/(2*mn);3.305322. 穩(wěn)健性分析：比較各種估計(jì)的差別普松分布的一階矩估計(jì)和普松分布的區(qū)間估計(jì)兩種估計(jì)比較接近，但普松分布的二階矩估計(jì)與前兩者差別很大。用三種估計(jì)的平均值作為參數(shù)lambda的估計(jì)應(yīng)是

4、比較穩(wěn)健的lambdaqz=(lambda1+lambda2+lambdazwqj)/3;4.32159即使不考慮二階矩估計(jì)，用其他兩種估計(jì)的平均值作為參數(shù)lambd的估計(jì)也應(yīng)是比較穩(wěn)健的lambdaqz=(lambda1+lambdazwqj)/2;3.295523. 協(xié)調(diào)性分析：比較估計(jì)分布與經(jīng)驗(yàn)離散頻率分布32診斷估計(jì)分布與經(jīng)驗(yàn)離散頻率分布圖形的協(xié)調(diào)性以上述參數(shù)作為普松分布的估計(jì)與經(jīng)驗(yàn)離散頻率分布函數(shù)的比較圖從圖形看出擬合不好.即使不考慮二階矩估計(jì)，用其他兩種估計(jì)的平均值作為參數(shù)lambd的估計(jì),以此參數(shù)作為普松分布的估計(jì)與經(jīng)驗(yàn)離散頻率分布函數(shù)的比較從圖形看出擬合也不好.盡管它看起來比前

5、一種好一點(diǎn)，即圖形的接近程度高一些，但兩種分布的差別還是較大。3.2:估計(jì)分布與經(jīng)驗(yàn)離散頻率分布理論的協(xié)調(diào)性用皮爾遜的Kf統(tǒng)計(jì)量進(jìn)行擬合檢驗(yàn)兩種情況都沒有通過.N=238,l=1,Kf=728,110>272,836（顯著性水平0.05）N=238,l=1,Kf=325,637>272,836（顯著性水平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),272.836是統(tǒng)計(jì)量Kf顯著性水平0.05臨界值。4.功能性分析：變換原始數(shù)據(jù)上述分布的各種參數(shù)為mnmyycdssyymaxyymi

6、nyypzyyp50mumyyylambdallambda2lambdaxxklambdal2383.28571102.52462803303.285713.285716.373723.1152693.07970lambdaulambdasxlambdazwqjlambdaqzxitazdddkkkpixitaLxitaUfun3.544943.502963.305324.321594.32159130.8967881.321597.32159平時(shí)柏松分布其中普松分布的一階矩估計(jì)和普松分布的二階矩估計(jì)不一致，這可能是估計(jì)不準(zhǔn)確的原因。為了消除這個(gè)原因，我們采用數(shù)據(jù)變換的方法解決。對(duì)數(shù)據(jù)yy做變

7、換：myyy=myy-minyy;c=(myyy)/(ssyy*2);0.51551d=0;yy=int(myyy+d)*c);(此時(shí)變換后的均值為變換前的c倍，而變換后的方差為變換前的cA2倍，令兩者相等，得到上述的變換公式，由于取整數(shù)函數(shù)int(*)有低估現(xiàn)象，可以用d來做小的調(diào)整)。用這些數(shù)據(jù)進(jìn)行上述同樣的分析得到圖形如下FREQUENCY0.0OS1.01.52.02.53.03.54.0yyMIDPOINTFREQUENCYofyy1.0804.0122.039從圖形看出效果很好。用皮爾遜的Kf統(tǒng)計(jì)量進(jìn)行擬合檢驗(yàn)通過.N=238,l=3,Kf=13.8675<270.684(顯

8、著性水平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是統(tǒng)計(jì)量Kf顯著性水平0.05臨界值。得到的分布參數(shù)為mnmyycdssyymaxyyminyypzyyp50mumyyylambdal2381.357140.5155101.19888401101.357141.35714lambda2lambdaxxklambdallambdaulambdasxlambdazwqjlambdaqzxitaz1.437311.2554241.233031.534151.505661.376751.390402.69714xitadddkkkpiLxitaUfunCLUCLLCL130.94735605.69714平時(shí)柏松分布2.632625.752770其中普松分布的一階矩估計(jì)（lambda1=1.35714和普松分布的二階矩估計(jì)（lambda2=1.43731屋本一致.這可能就是估計(jì)準(zhǔn)確的原因。數(shù)據(jù)變換方法使得兩個(gè)估計(jì)接近相等。5.經(jīng)濟(jì)性分析：變換控制圖根據(jù)變換后的估計(jì)分布的控制圖，變回到原來分布的控制圖的控制限為UCL=5.75277CL=2.63262LCL=0要求每五個(gè)數(shù)據(jù)進(jìn)行平均得到一個(gè)點(diǎn)（低于五個(gè)數(shù)據(jù)平均時(shí)，