• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于遷徙數(shù)據(jù)估計(jì)武漢感染新型冠狀病毒的人員數(shù)量

    2020-06-17 06:44:12原子霞賈祖瑤
    關(guān)鍵詞:估計(jì)值中位數(shù)比率

    楊 政,原子霞,賈祖瑤

    (1. 電子科技大學(xué)經(jīng)濟(jì)與管理學(xué)院 成都 611731;2. 電子科技大學(xué)數(shù)學(xué)科學(xué)學(xué)院 成都 611731)

    新型冠狀病毒肺炎已經(jīng)成為國(guó)際關(guān)注的重大緊急公共衛(wèi)生事件,給人民的生命和生活造成嚴(yán)重危害。因此,阻擊病毒傳染成了全國(guó)人民的共同戰(zhàn)役。從2020 年1 月23 日武漢開始“封城”后,各個(gè)省市采取多種防控措施。居民按照專家建議減少外出活動(dòng),在家隔離以降低被感染的風(fēng)險(xiǎn)。

    2020 年1 月23 日?2 月4 日,武漢市確診人數(shù)和疑似病例的數(shù)據(jù)不斷升高。國(guó)家在武漢投入更多的力量醫(yī)治確診病人,如緊急調(diào)撥物資、建立火神山和雷神山醫(yī)院、派出多批次的支援醫(yī)療隊(duì)等。這些措施給全國(guó)人民帶來(lái)了戰(zhàn)勝病毒的信心。此時(shí),明確武漢市感染者的數(shù)量對(duì)于防控、診斷和治療有重要意義。那么,武漢市目前受感染的人數(shù)有多少?這是本文擬研究的問(wèn)題。

    從Elsevier 數(shù)據(jù)庫(kù)查閱到新型冠狀病毒的相關(guān)論文大約有70 余篇,大致分為兩類。一類側(cè)重于從醫(yī)學(xué)方面探討新型冠狀病毒的來(lái)源[1]、發(fā)現(xiàn)和臨床診斷[2]、病毒基因分析[3]、公眾心理健康[4]以及如何控制病毒流行[5]等問(wèn)題。

    另一類論文利用大數(shù)據(jù)、傳播動(dòng)力學(xué)模型、統(tǒng)計(jì)計(jì)算方法等工具對(duì)疫情進(jìn)行了預(yù)測(cè)分析。文獻(xiàn)[6]基于包括“易感態(tài)-潛伏態(tài)-感染態(tài)-移除態(tài)”的SEIR 倉(cāng)室模型,對(duì)病毒的基本再生數(shù)進(jìn)行估計(jì)。以《人民日?qǐng)?bào)》新型冠狀病毒肺炎疫情實(shí)時(shí)動(dòng)態(tài)數(shù)據(jù)為基準(zhǔn),估計(jì)基本再生數(shù)在2.8~3.3 之間;以國(guó)外同行預(yù)測(cè)的感染人數(shù)為基準(zhǔn),基本再生數(shù)在3.2~3.9 之間。文獻(xiàn)[7]利用傳播動(dòng)力學(xué)模型,對(duì)新型冠狀病毒肺炎傳播風(fēng)險(xiǎn)進(jìn)行了預(yù)測(cè)分析。該文利用2020 年1 月10 日?1 月22 日的報(bào)告疫情數(shù)據(jù),采用動(dòng)力學(xué)模型和統(tǒng)計(jì)計(jì)算方法預(yù)測(cè)基本再生數(shù)為6.47(95%置信區(qū)間為 5.71~7.23),給出了疫情的達(dá)峰時(shí)間、峰值及最終感染規(guī)模,按照2020 年1 月22 日前的控制措施,疫情將在3 月10 日左右達(dá)到峰值。文獻(xiàn)[8]分析了分布在全國(guó)31 個(gè)省市自治區(qū)、552 家醫(yī)院的1 099 個(gè)確診病例的臨床特征、潛伏期、診斷情況、治療方式等要素,發(fā)現(xiàn)新型冠狀病毒感染的平均潛伏期為3 天。文獻(xiàn)[9]預(yù)測(cè)了新型冠狀病毒感染者的人數(shù),估計(jì)2020 年1 月25 日的感染人數(shù)約7.5 萬(wàn)人。文獻(xiàn)[10]根據(jù)自然增長(zhǎng)規(guī)律動(dòng)態(tài)提出數(shù)據(jù)驅(qū)動(dòng)的預(yù)測(cè)方法,跟蹤疫情發(fā)展并檢測(cè)干預(yù)措施的有效性。在2020 年2 月5 日預(yù)測(cè)約4 天后(2 月9 日)達(dá)到峰值,確診病例總數(shù)將在3.7 萬(wàn)~4.4 萬(wàn)之間。

    上述研究并沒(méi)有直接針對(duì)武漢市的感染者人數(shù)進(jìn)行預(yù)估。本文在疫情前期把武漢市所有民眾看作一個(gè)樣本總體,離開武漢和留在武漢是由同一個(gè)總體分布中抽取的兩組隨機(jī)樣本。本文把離開武漢的民眾視為實(shí)驗(yàn)組樣本,把留在武漢的民眾視為對(duì)照組樣本。武漢市受到醫(yī)護(hù)人員不足和醫(yī)療物資緊缺的約束,對(duì)照組樣本的確診人數(shù)低于實(shí)際被感染人數(shù)。所以,利用實(shí)驗(yàn)組樣本的確診數(shù)據(jù),分析其統(tǒng)計(jì)分布的數(shù)字特征,借鑒這些數(shù)字特征對(duì)武漢市的感染人數(shù)進(jìn)行估計(jì)。簡(jiǎn)單來(lái)說(shuō),就是利用實(shí)驗(yàn)組樣本的統(tǒng)計(jì)參數(shù),估計(jì)對(duì)照組樣本未受約束時(shí)的發(fā)展?fàn)顩r。

    1 數(shù)據(jù)和假設(shè)

    本文利用百度遷徙數(shù)據(jù)來(lái)估計(jì)武漢目前受感染的人數(shù)。統(tǒng)計(jì)從武漢遷入人員數(shù)量排名前50 位城市的人數(shù),遷徙時(shí)間從2020 年1 月10 日至1 月22 日。表1 給出了排名前50 位城市從武漢遷入的人數(shù)。

    表1 排名前50 城市從武漢遷入的人數(shù)

    為了估算武漢市感染新型冠狀病毒的人數(shù),本文做出如下假設(shè)。

    假設(shè)1:2020 年1 月10 日?1 月22 日離開武漢的遷徙數(shù)據(jù)是準(zhǔn)確的。

    由于并沒(méi)有一手的人員流動(dòng)數(shù)據(jù),故以網(wǎng)絡(luò)報(bào)道的百度遷徙數(shù)據(jù)為準(zhǔn)。百度遷徙數(shù)據(jù)是根據(jù)遷徙人員的手機(jī)在不同地點(diǎn)的定位統(tǒng)計(jì)的,準(zhǔn)確性較高。例如本文根據(jù)遷徙數(shù)據(jù)計(jì)算從武漢到信陽(yáng)的人數(shù)是7.19 萬(wàn)人(新聞報(bào)道的數(shù)據(jù)是8.046萬(wàn)人,包括了22 日之后的遷徙人數(shù))。因此,在2020 年1 月22 日之前的遷徙數(shù)據(jù)具有較高的可信度。在50 個(gè)城市中,遷徙人員是確診新冠病毒肺炎患者的直接來(lái)源,至少在潛伏期間遷徙數(shù)據(jù)是一個(gè)最主要的影響因子。因此,50 個(gè)城市的確診病例應(yīng)該是以武漢遷入人員為主。在潛伏期之后,武漢遷入人員經(jīng)過(guò)病毒潛伏期,確診比率應(yīng)該會(huì)下降。在病毒的潛伏期內(nèi),50 個(gè)城市本地居民受到武漢遷入人員傳染,確診比率應(yīng)該會(huì)升高。從武漢遷入50 個(gè)城市的人員數(shù)量仍然是一個(gè)基礎(chǔ)變量,它持續(xù)影響了后續(xù)感染患者的數(shù)量。

    假設(shè)2:從武漢遷徙到50 個(gè)城市的人數(shù)在未來(lái)一段時(shí)間內(nèi)保持不變。

    這個(gè)假設(shè)是為了對(duì)應(yīng)“封城”后武漢的人數(shù)保持不變。計(jì)算50 個(gè)城市感染患者人數(shù)的基數(shù)也應(yīng)該保持不變。事實(shí)上,這個(gè)假設(shè)對(duì)個(gè)別城市可能不成立,如溫州市,據(jù)報(bào)道2020 年1 月22 日之后從武漢到溫州的人數(shù)比之前的人數(shù)更多。這對(duì)研究結(jié)果有一定影響,因此在后面研究中做了一些修正。

    假設(shè)3:留在武漢和離開武漢人員感染病毒的概率是同一個(gè)分布。

    這個(gè)假設(shè)是估計(jì)武漢感染人數(shù)的一個(gè)重要前提。假設(shè)意味著離開和留在武漢的人員都是同一個(gè)概率分布的樣本。這個(gè)假設(shè)在2020 年1 月22 日之后幾天可能是成立的。由于本文整理數(shù)據(jù)是2020 年從1 月29 日開始,距離1 月22 日已經(jīng)有7 天的時(shí)間間隔。通過(guò)檢驗(yàn)湖北省內(nèi)城市和省外城市的確診比率,發(fā)現(xiàn)兩者并不是同一個(gè)分布。為此,本文把50 個(gè)城市分成兩組。一組是湖北省內(nèi)的15 個(gè)城市,另一組是省外的35 個(gè)城市。即使這兩組樣本的分布不同,仍然假設(shè)武漢屬于這兩組樣本中的某一種情況。如果武漢不屬于這兩組樣本,那會(huì)是一種最差的情況。

    假設(shè)4:樣本期內(nèi)留在武漢和離開武漢確診的人員都是源于相同的病毒感染模式。

    這個(gè)假設(shè)表示武漢市的確診人員和其他城市的確診人員是相同的感染模式。比如,其他城市前期的確診人員都是從武漢遷入被確診,后期的部分確診人員是受遷移人員傳染而被確診。武漢市早期的確診人員和其他城市早期的確診人員是同一批感染者,后期的部分人員受他人傳染被確診。

    除了收集50 個(gè)城市從武漢遷徙來(lái)的人員數(shù)量,利用百度新型冠狀病毒肺炎-疫情實(shí)時(shí)大數(shù)據(jù)報(bào)告,收集整理這些城市從2020 年1 月29 日至2 月9 日的確診人數(shù),并從萬(wàn)德(Wind)數(shù)據(jù)庫(kù)收集武漢市的戶籍人口數(shù)據(jù)。

    2 湖北省內(nèi)外城市確診比率的計(jì)算

    用i=1,2,···,50表示50 個(gè)城市中的第i 個(gè)城市 ,用 xt,i表 示第i 個(gè)城市在第t天的累積確診人數(shù),用 yi表示從武漢遷出到第i 個(gè)城市的人數(shù)。計(jì)算第i 個(gè)城市每天的確診比率 pt,i為:

    這樣得到50 個(gè)城市累計(jì)確診比率的時(shí)間序列數(shù)據(jù)。

    接下來(lái),用R 軟件對(duì)各城市每天的數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)分析。為了避免個(gè)別城市的特殊值影響整體分析,除了最大值和最小值,還在描述性統(tǒng)計(jì)中增加了第二大值和第二小值。表2 給出湖北省內(nèi)和省外城市在確診比率方面的描述性統(tǒng)計(jì)。由于確診人數(shù)是累計(jì)值而不是每天增量值。因此均值、中位數(shù)、最大值、最小值和標(biāo)準(zhǔn)差等隨著時(shí)間增加而變大。對(duì)比省內(nèi)外城市,省外城市在均值,中位數(shù)、標(biāo)準(zhǔn)差、最大值及第二大值這些統(tǒng)計(jì)量的數(shù)值上都大于當(dāng)日的省內(nèi)城市。省外城市確診比率的最小值和第二小值在5 日和6 日之前都小于省內(nèi)城市,之后都大于省內(nèi)城市。原因是省內(nèi)城市早期有確診人員,初始值較大,但是受限于醫(yī)務(wù)人員的不足,確診比率增長(zhǎng)較慢。省外城市的初始值小,隨著省外城市的醫(yī)療資源充足,潛在患者被迅速確診,確診比率的最小值和第二小值迅速超過(guò)省內(nèi)城市。

    表2 中雅克貝拉(Bera–Jarque, JB)檢驗(yàn)統(tǒng)計(jì)量用于檢驗(yàn)分布是否屬于正態(tài)分布。對(duì)于每日確診比率,檢驗(yàn)原假設(shè) H0:pt,服從標(biāo)準(zhǔn)正態(tài)分布。根據(jù)文獻(xiàn)[11],定義JB 統(tǒng)計(jì)量為:

    表2 湖北省內(nèi)和省外城市確診比率的描述性統(tǒng)計(jì)

    表3 每天確診比率極值對(duì)應(yīng)的城市

    省外城市在2020 年1 月29 日?2 月9 日期間的變化不大。前3 日確診比率最高的是浙江溫州和杭州。隨著浙江采取嚴(yán)格的防控措施,2 月1 日之后確診比率第二大值出現(xiàn)在深圳市。確診比率最小值一直由貴陽(yáng)市保持。第二小值在洛陽(yáng)市、廈門市和石家莊市中變換。

    對(duì)比湖北省內(nèi)城市和省外城市的表現(xiàn)??紤]到省內(nèi)和省外城市的樣本數(shù)和分布相同,采用文獻(xiàn)[12]提出的F 檢驗(yàn)做省內(nèi)外的均值相等性檢驗(yàn),采用文獻(xiàn)[13]的 χ2檢驗(yàn)做中位數(shù)相等性檢驗(yàn)。應(yīng)用R 軟件檢驗(yàn)均值和中位數(shù)的相等性,檢驗(yàn)結(jié)果見(jiàn)表4。

    從均值檢驗(yàn)結(jié)果看,表4 顯示省內(nèi)外均值相等的原假設(shè)在5%水平被顯著拒絕,說(shuō)明湖北省內(nèi)城市和省外城市的均值差異明顯。從中位數(shù)檢驗(yàn)看,2020 年1 月30 日?2 月1 日,中位數(shù)相等的原假設(shè)在1%水平上被顯著拒絕。2 月2 日的中位數(shù)檢驗(yàn)在10%水平上并不顯著。在2 月3 日?6 日的結(jié)果出現(xiàn)反轉(zhuǎn),顯示省內(nèi)中位數(shù)持續(xù)低于省外的中位數(shù)。一個(gè)好的信號(hào)出現(xiàn)在2 月7 日,p 值在10%水平上不拒絕省內(nèi)外中位數(shù)相等的原假設(shè)。這說(shuō)明湖北省內(nèi)和省外的確診狀況暫時(shí)進(jìn)入一個(gè)新階段。2 月8 日和9 日的中位數(shù)檢驗(yàn)結(jié)果強(qiáng)化了這一結(jié)論。

    3 武漢市感染病毒人數(shù)的估計(jì)

    根據(jù)式(1),估計(jì)每日武漢感染人數(shù)為:

    表4 均值和中位數(shù)檢驗(yàn)

    式中,y =883.73萬(wàn) 是武漢市的戶籍人數(shù); pt,j是表2的比率,表示第t 日第j 種情況下的數(shù)值; j 分別表示均值、中位數(shù)、最大值、第二大值、第二小值和最小值這6 種情況。常見(jiàn)的統(tǒng)計(jì)估計(jì)應(yīng)該包括某些置信水平下的區(qū)間估計(jì),比如估計(jì)武漢感染人數(shù)的95%區(qū)間,在本文中并沒(méi)有做區(qū)間估計(jì)。由于武漢市的情況很特殊,也可能不屬于省內(nèi)和省外的兩種分布,匯報(bào)區(qū)間估計(jì)的意義不大。直接采用4 種極值比率來(lái)估計(jì),這樣能夠看到極端情況下武漢市感染人數(shù)的估計(jì)值。

    表5 給出了6 種情況估計(jì)的武漢感染人數(shù)。表5 的最后一列是每日公布的武漢市確診人數(shù)。為了更好地理解估計(jì)結(jié)果,把估計(jì)值除以每日公布的確診人數(shù),得到估計(jì)值和公布確診人數(shù)的比值,結(jié)果見(jiàn)表6。表5 和6 的結(jié)果總結(jié)為以下3 點(diǎn)。

    表5 6 種情況下武漢感染人數(shù)的估計(jì)

    1) 按照確診比率的最小值 (即最小比率)來(lái)估計(jì)?;谑?nèi)城市確診比率最小值(即潛江市),估計(jì)武漢市受感染人數(shù)。除了前3 日的估計(jì)人數(shù)低于確診人數(shù)外,后面5 日的估計(jì)人數(shù)都高于確診人數(shù)。自2020 年2 月8 日開始,武漢市確診人數(shù)開始大于估計(jì)的感染人數(shù)。再按第二小值的比率估計(jì),1 月29 日的估計(jì)值是1 971 人,和確診的1 905 人較為接近。2 月9 日,根據(jù)第二小值(即恩施市)的確診比率估計(jì)武漢市感染人數(shù),估計(jì)值大約是確診人數(shù)的1.1 倍。

    從省外城市來(lái)看,按照確診比率的最小值(即貴陽(yáng)市)估計(jì)2020 年2 月1 日武漢市的受感染人數(shù),估計(jì)值和確診人數(shù)持平。用第二小值(即洛陽(yáng)市)估計(jì),在2 月1 日的估計(jì)值已經(jīng)開始大于確診人數(shù)。按照省外城市確診比率的最小值,即最樂(lè)觀的估計(jì),2 月2 日之后武漢市感染人數(shù)的估計(jì)值全部大于確診人數(shù)。

    2) 按照省內(nèi)城市的平均值估計(jì),表5 顯示武漢在2020 年1 月29 日?2 月9 日的感染人數(shù)分別是4 483人和31 531 人。表6 展示了比值的動(dòng)態(tài)變化,由1 月29 日的2.4 倍開始增加,到2 月1 日到達(dá)峰值即3.4 倍,之后比值開始持續(xù)減少,2 月9 日的比值是2.1 倍。

    在按照省外城市的平均值估計(jì),從2020 年1 月29 日?2 月9 日的感染人數(shù)分別是1.304 0 萬(wàn)人和5.324 8 萬(wàn)人。表6 的比值隨時(shí)間變化的動(dòng)態(tài)特征和省內(nèi)情況類似,從1 月29 日的6.8 倍到1 月31 日達(dá)到峰值即8.3 倍,比值從2 月1 日開始持續(xù)下降,到2 月9 日估計(jì)的感染人數(shù)是確診人數(shù)的3.6 倍。

    表6 武漢感染人數(shù)的估計(jì)和公布確診人數(shù)的比值

    按中位數(shù)估計(jì),省內(nèi)城市的比值從1 月31 日的3 倍減少到2 月9 日的2 倍。利用省外城市估計(jì)的感染人數(shù)和確診人數(shù)的比值,從1 月29 日的4.7 倍增加到1 月31 日的6 倍,再逐步減少到2 月9 日的2.6 倍。

    3) 從感染確診比率的最大值(即最大比率)來(lái)估計(jì)。按照省內(nèi)城市的最大值(即十堰市和隨州市)估計(jì),從2020 年1 月29 日?2 月9 日,估計(jì)的武漢感染人數(shù)是當(dāng)日確診人數(shù)的4.5 倍和3.9 倍,期間比值在2 月4 日達(dá)到最大的5.9 倍。按照省內(nèi)城市的第二大值(分別是荊門市和十堰市)來(lái)估計(jì),1 月29 日和2 月9 日的估計(jì)值分別是0.811 7 萬(wàn)人和4.598 3 萬(wàn)人。武漢市估計(jì)的感染人數(shù)是確診人數(shù)的4.3 倍和3.1 倍。

    按照省外城市的最大值(即溫州市)估計(jì),從2020 年1 月29 日?2 月9 日,估計(jì)的感染人數(shù)和確診人數(shù)比值從45.6 倍(1 月29 日)增加到65.5 倍(1 月31 日),再 逐 步 減 少 到22.8 倍(2 月9 日)。1 月29 日的估計(jì)值是8.685 萬(wàn)人,2 月9 日的估計(jì)值是34.130 萬(wàn)人。用省外城市的第二大值來(lái)估計(jì),1 月29 日(杭州市)和2 月9 日(深圳市)的估計(jì)值分別是2.666 9 萬(wàn)人和12.392 5 萬(wàn)人,估計(jì)的感染人數(shù)分別是當(dāng)日確診人數(shù)的14 倍和8.3 倍。

    用省外城市的最大比率(即溫州市)估計(jì)出武漢市的感染人數(shù)在2 月9 日是34.130 萬(wàn)人,這個(gè)結(jié)果令人吃驚。追查從武漢回到溫州的人數(shù),溫州市副市長(zhǎng)在2020 年1 月29 日采訪中提到:“武漢‘封城’后,1 月23 日至27 日5 天,仍然有1.88 萬(wàn)人從湖北特別是武漢到達(dá)溫州,平均每天有3 600 多人”。因此從武漢回到溫州是3.04 萬(wàn)人,大于百度遷徙數(shù)據(jù)計(jì)算的2020 年1 月10 日至22 日的1.16 萬(wàn)人,回到溫州的實(shí)際人數(shù)增加了1.6 倍。假設(shè)溫州的累計(jì)確診人數(shù)是從3.04 萬(wàn)人中得到的,那么估算武漢市感染人數(shù)大約為13.127 萬(wàn)人(34.130 萬(wàn)人/2.6),仍然高于由第二大值(即深圳市)確診比率估計(jì)的12.393 萬(wàn)人。調(diào)整后的估計(jì)感染人數(shù)與當(dāng)日確診人數(shù)的比值是8.7 倍,高于第二大值的8.3 倍。從表6 看到,省外城市感染比率第二大的城市在1 月29 日是杭州。從2 月1 日之后就是深圳市。顯然,從確診率高的溫州市和深圳市估計(jì)武漢市的感染人數(shù),結(jié)果較為一致。

    4 省內(nèi)外城市均值和中位數(shù)的Bootstrap估計(jì)

    由于研究樣本較少,尤其是省內(nèi)的15 個(gè)樣本屬于小樣本情形。對(duì)均值和中位數(shù)的估計(jì)可能會(huì)有一些影響。為此,采用Bootstrap 方法[14]重新估計(jì)每日的均值和中位數(shù)。具體步驟為:

    對(duì)湖北省內(nèi)外城市的每日樣本進(jìn)行Bootstrap抽樣B=10 000 次,計(jì)算得到Bootstrap 均值和中位數(shù),結(jié)果見(jiàn)表7。圖1 給出從2020 年1 月29 日到2 月9 日每天的直方圖。

    表7 基于bootstrap 抽樣計(jì)算的省內(nèi)外城市的均值和中位數(shù)

    把表7 的Bootstrap 均值和中位數(shù)與表2 的均值和中位數(shù)對(duì)比,二者數(shù)據(jù)非常接近,說(shuō)明均值和中位數(shù)的結(jié)果是魯棒的。

    圖1 的直方圖反映了湖北省內(nèi)城市和省外城市確診比率在均值上的差異。直方圖是經(jīng)驗(yàn)分布的直觀表現(xiàn)。圖1 表明省內(nèi)外均值的差異是多方面的。省內(nèi)確診比率的均值始終小于省外均值,峰值高表明確診比率在均值周圍的頻次非常高,表明省內(nèi)城市的確診比率非常集中,方差小同樣表明確診比率在均值周圍變化小。省外城市確診比率的特點(diǎn)是方差更大。

    圖1 的另一個(gè)重要特征是隨著時(shí)間變化,省內(nèi)城市的直方圖和省外城市的直方圖產(chǎn)生了交集,當(dāng)時(shí)間增加,交集重合的部分越來(lái)越多。從表7 看到,省內(nèi)城市的均值以更快的速度增加,1 月29 日是0.508,2 月9 日是3.570,大約增加了6 倍。省外城市的均值在1 月29 日是1.473,2 月9 日是6.029,大約增加了3 倍。從表7 還可以看到從1 月29 日到2 月9 日,省內(nèi)城市的中位數(shù)增加了6.5 倍,而省外城市的中位數(shù)增加了3.2 倍。這表明隨著湖北省內(nèi)城市醫(yī)療條件的改善,確診比率提高得越來(lái)越快,逐漸跟上省外城市的確診趨勢(shì)。只有當(dāng)省內(nèi)外的均值和中位數(shù)在統(tǒng)計(jì)檢驗(yàn)上不再有顯著差異時(shí),省內(nèi)外的疫情達(dá)到了相同的水平。

    圖1 基于10 000 次Bootstrap 抽樣估計(jì)均值的直方圖

    5 結(jié) 束 語(yǔ)

    5.1 研究結(jié)論

    利用2020 年1 月10 日?1 月22 日的百度遷徙數(shù)據(jù),本文統(tǒng)計(jì)從武漢市到全國(guó)50 個(gè)城市的遷移人數(shù)。同時(shí),收集2020 年1 月29 日?2 月9 日這12 天內(nèi)50 個(gè)城市感染新型冠狀病毒確診人數(shù)的數(shù)據(jù)。首先,利用統(tǒng)計(jì)方法計(jì)算了感染新型冠狀病毒人數(shù)占遷移人數(shù)的比率。其次,對(duì)省內(nèi)外每天的確診比率進(jìn)行描述性統(tǒng)計(jì),以及均值和中位數(shù)相等性檢驗(yàn)。接下來(lái),根據(jù)省內(nèi)外的統(tǒng)計(jì)結(jié)果,對(duì)武漢市的感染人數(shù)進(jìn)行估計(jì)。最后對(duì)均值和中位數(shù)進(jìn)行了Bootstrap 抽樣計(jì)算,均值和中位數(shù)結(jié)果具有穩(wěn)健性。本文研究得到以下結(jié)論。

    1) 通過(guò)對(duì)比發(fā)現(xiàn),湖北省內(nèi)城市確診人數(shù)的均值和中位數(shù)都低于省外城市的均值和中位數(shù)。原因是疫情初期湖北省內(nèi)的醫(yī)療資源不足,許多感染病人還未得到有效的診斷和治療。潛在病人尚未被發(fā)現(xiàn),這需要特別重視。隨州、十堰、襄陽(yáng)、宜昌市和荊門等城市,在樣本期內(nèi)的確診比率處于15 個(gè)省內(nèi)城市的前兩位。從均值檢驗(yàn)來(lái)看,湖北省內(nèi)城市和省外城市的差異是顯著的。從中位數(shù)檢驗(yàn)來(lái)看,湖北省內(nèi)城市和省外城市的差異在2 月7 日發(fā)生了改變,不拒絕在10%水平上中位數(shù)相等的假設(shè)。這表明湖北省內(nèi)城市感染者的確診逐漸趕上省外城市確診的速度。當(dāng)省內(nèi)和省外城市在均值檢驗(yàn)也無(wú)顯著差異時(shí),才能認(rèn)為省內(nèi)外疫情狀態(tài)達(dá)到同一個(gè)水平,疫情防控進(jìn)入一個(gè)新的階段。

    2) 從最近一天(2 月9 日)的情況來(lái)看,利用省內(nèi)城市確診比率的均值和中位數(shù)估計(jì)武漢市的感染人數(shù),估計(jì)值是確診人數(shù)的2.1 倍和2 倍。利用最大值和第二大值估計(jì),感染人數(shù)是確診人數(shù)的3.9 倍和3.1 倍。利用省外城市的均值和中位數(shù)估計(jì)武漢市的感染人數(shù),估計(jì)值是確診人數(shù)的3.6 倍和2.6 倍。用最大值和第二大值估計(jì)的武漢市感染人數(shù),是確診人數(shù)的8.7 倍(修正后)和8.3 倍。這些結(jié)果無(wú)不說(shuō)明武漢市內(nèi)有很多潛在的感染患者尚未得到診斷。根據(jù)通報(bào)信息,武漢市前期已經(jīng)征用和開辟了9 000 張床位,雷神山醫(yī)院的1 000~1 500 張床位,火神山醫(yī)院的700~1 000 張床位。2 月4 日武漢市征用11 家場(chǎng)館改造成“方艙醫(yī)院”,改造完成后,可提供萬(wàn)余張床位。這些床位數(shù)量加在一起仍然小于估計(jì)的感染人數(shù)。

    3) 利用Bootstrap 方法重新估計(jì)湖北省內(nèi)外城市確診比率的均值和中位數(shù)。穩(wěn)健性的結(jié)論進(jìn)一步支持對(duì)武漢市感染患者的預(yù)測(cè)結(jié)果。

    5.2 進(jìn)一步的討論

    首先,由于作者不具備醫(yī)學(xué)方面的專業(yè)知識(shí),無(wú)法從傳染病模型、病毒潛伏期、基本再生數(shù)及感染傳播機(jī)制等方面進(jìn)行分析。本文的假設(shè)也忽略了病毒二代傳播在不同地方的差異性,這些差異性對(duì)估計(jì)結(jié)果會(huì)有一定的影響,使得估計(jì)值和實(shí)際感染人數(shù)有一定偏差。

    其次,基礎(chǔ)數(shù)據(jù)的準(zhǔn)確性會(huì)影響估計(jì)結(jié)果。由于研究數(shù)據(jù)是根據(jù)網(wǎng)絡(luò)上百度遷移數(shù)據(jù)整理得到,而實(shí)際情況更復(fù)雜,整理的數(shù)據(jù)與實(shí)際數(shù)據(jù)有差異。50 個(gè)城市遷徙人員的基數(shù)變小使得計(jì)算的確診率偏高,導(dǎo)致武漢感染人數(shù)的估計(jì)值也偏高。

    最后,文中沒(méi)有考慮50 個(gè)城市每日增加確診人數(shù)的動(dòng)態(tài)特征,利用面板分析方法研究動(dòng)態(tài)數(shù)據(jù)會(huì)得到新的啟示,比如判斷疫情拐點(diǎn)的出現(xiàn)。利用更多的數(shù)據(jù)信息,更復(fù)雜的統(tǒng)計(jì)和大數(shù)據(jù)研究方法,研究結(jié)論將會(huì)更豐富。

    猜你喜歡
    估計(jì)值中位數(shù)比率
    一類具有時(shí)滯及反饋控制的非自治非線性比率依賴食物鏈模型
    一道樣本的數(shù)字特征與頻率分布直方圖的交匯問(wèn)題
    統(tǒng)計(jì)信息
    2018年4月世界粗鋼產(chǎn)量表(續(xù))萬(wàn)噸
    中位數(shù)計(jì)算公式及數(shù)學(xué)性質(zhì)的新認(rèn)識(shí)
    一種適用于微弱信號(hào)的新穎雙峰值比率捕獲策略
    2015年中考數(shù)學(xué)模擬試題(五)
    2015年中考數(shù)學(xué)模擬試題(二)
    導(dǎo)學(xué)案不能淪落為“習(xí)題單”:以“中位數(shù)和眾數(shù)”的導(dǎo)學(xué)案為例
    2014年5月世界粗鋼產(chǎn)量表萬(wàn)噸
    亚洲成av人片在线播放无| 香蕉av资源在线| 十八禁网站免费在线| 免费搜索国产男女视频| 99热这里只有精品一区| 日韩国内少妇激情av| 神马国产精品三级电影在线观看| 亚洲精品色激情综合| 两个人视频免费观看高清| 夜夜夜夜夜久久久久| 麻豆一二三区av精品| 三级国产精品欧美在线观看| 国产一区二区在线av高清观看| 亚洲最大成人手机在线| 啦啦啦韩国在线观看视频| av.在线天堂| 日韩欧美一区二区三区在线观看| 偷拍熟女少妇极品色| 国产黄色视频一区二区在线观看 | 成人无遮挡网站| 老司机福利观看| 欧美人与善性xxx| 亚洲激情五月婷婷啪啪| 一级毛片久久久久久久久女| 又爽又黄无遮挡网站| 赤兔流量卡办理| 欧美一区二区精品小视频在线| 国产高清视频在线观看网站| 日韩亚洲欧美综合| 可以在线观看毛片的网站| 寂寞人妻少妇视频99o| 免费电影在线观看免费观看| 男人和女人高潮做爰伦理| 午夜久久久久精精品| 国语自产精品视频在线第100页| 少妇丰满av| 成人午夜高清在线视频| 国产精品永久免费网站| 精品一区二区三区视频在线观看免费| 精品久久国产蜜桃| 99热这里只有是精品在线观看| 真人做人爱边吃奶动态| 午夜福利高清视频| av在线天堂中文字幕| 麻豆乱淫一区二区| 少妇熟女aⅴ在线视频| av中文乱码字幕在线| 国产一区二区三区在线臀色熟女| 91精品国产九色| 久久久久久久久久成人| 日本 av在线| 免费搜索国产男女视频| 日本一本二区三区精品| 亚洲国产高清在线一区二区三| 我的女老师完整版在线观看| 亚洲欧美清纯卡通| 九九在线视频观看精品| 亚州av有码| 国产亚洲av嫩草精品影院| 国产综合懂色| 久久久久久九九精品二区国产| 精品久久久久久久久av| 成人性生交大片免费视频hd| 亚洲无线观看免费| 成人三级黄色视频| 综合色丁香网| 亚洲精品乱码久久久v下载方式| 日韩精品中文字幕看吧| 亚州av有码| 1024手机看黄色片| 人妻丰满熟妇av一区二区三区| 精品一区二区三区人妻视频| 亚洲内射少妇av| 国产人妻一区二区三区在| 国产探花在线观看一区二区| 日产精品乱码卡一卡2卡三| 亚洲性夜色夜夜综合| a级毛色黄片| 久久久久免费精品人妻一区二区| 村上凉子中文字幕在线| 一个人观看的视频www高清免费观看| 成年av动漫网址| 国产成人freesex在线 | 少妇的逼水好多| 日韩精品青青久久久久久| 日韩强制内射视频| 夜夜夜夜夜久久久久| 一个人看视频在线观看www免费| 亚洲自拍偷在线| 亚洲七黄色美女视频| 亚洲美女搞黄在线观看 | 麻豆成人午夜福利视频| 日韩成人伦理影院| 国产精品av视频在线免费观看| 亚洲婷婷狠狠爱综合网| 中文字幕久久专区| 国产三级在线视频| 性插视频无遮挡在线免费观看| 欧美国产日韩亚洲一区| 久久久久久九九精品二区国产| a级毛片a级免费在线| 久久99热6这里只有精品| 亚洲无线观看免费| eeuss影院久久| 少妇猛男粗大的猛烈进出视频 | 欧美日韩一区二区视频在线观看视频在线 | 午夜福利在线在线| 日日摸夜夜添夜夜爱| 内地一区二区视频在线| 国产精品亚洲一级av第二区| 国产av不卡久久| 狂野欧美白嫩少妇大欣赏| 欧美日本视频| 国产高清激情床上av| 久久久久国内视频| 久久久欧美国产精品| 免费看a级黄色片| 最近中文字幕高清免费大全6| 九九爱精品视频在线观看| 夜夜看夜夜爽夜夜摸| 天堂√8在线中文| 久久精品夜色国产| 免费无遮挡裸体视频| 我的女老师完整版在线观看| 欧美色欧美亚洲另类二区| 特级一级黄色大片| 久久久久久久久大av| 老司机影院成人| 午夜视频国产福利| 最新中文字幕久久久久| 在线免费十八禁| 亚洲中文日韩欧美视频| 日韩精品有码人妻一区| 麻豆国产av国片精品| 亚洲人成网站在线观看播放| 一区福利在线观看| 国产真实乱freesex| 亚洲成av人片在线播放无| 日韩中字成人| 亚洲国产色片| 久久久久国产网址| 免费高清视频大片| 国产精品伦人一区二区| 午夜福利在线观看免费完整高清在 | 我要搜黄色片| 国产 一区 欧美 日韩| 1000部很黄的大片| 精品免费久久久久久久清纯| 中文字幕人妻熟人妻熟丝袜美| 人妻夜夜爽99麻豆av| 听说在线观看完整版免费高清| 蜜桃亚洲精品一区二区三区| 九九在线视频观看精品| 国产精品女同一区二区软件| 最近最新中文字幕大全电影3| 成人毛片a级毛片在线播放| 舔av片在线| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品久久久久久毛片| 国产av不卡久久| 欧美3d第一页| 悠悠久久av| 俄罗斯特黄特色一大片| 精品一区二区三区av网在线观看| 亚洲精品成人久久久久久| 天天躁日日操中文字幕| 午夜日韩欧美国产| 久久久久久伊人网av| 亚洲国产色片| 人人妻人人澡人人爽人人夜夜 | 国产大屁股一区二区在线视频| 国产高潮美女av| 一级毛片久久久久久久久女| 国产精品一区二区免费欧美| 日产精品乱码卡一卡2卡三| 卡戴珊不雅视频在线播放| 日韩大尺度精品在线看网址| 久久精品夜夜夜夜夜久久蜜豆| 毛片女人毛片| 色综合色国产| a级毛片a级免费在线| 蜜桃久久精品国产亚洲av| 麻豆久久精品国产亚洲av| 黄色视频,在线免费观看| 日本在线视频免费播放| 亚洲精品亚洲一区二区| 午夜老司机福利剧场| 亚洲人成网站在线观看播放| 久久久久久久久久黄片| 成人鲁丝片一二三区免费| 国产老妇女一区| 丰满乱子伦码专区| 欧美国产日韩亚洲一区| 国产高清不卡午夜福利| 超碰av人人做人人爽久久| 精品久久久久久久久久久久久| 又粗又爽又猛毛片免费看| 69av精品久久久久久| 天堂网av新在线| 欧美日韩在线观看h| 日本一二三区视频观看| 日本熟妇午夜| or卡值多少钱| 亚洲在线自拍视频| 成人国产麻豆网| 简卡轻食公司| 日本黄色视频三级网站网址| 91久久精品国产一区二区成人| 日韩大尺度精品在线看网址| avwww免费| 久99久视频精品免费| 大香蕉久久网| 国产精品精品国产色婷婷| av黄色大香蕉| 国产真实伦视频高清在线观看| 亚洲欧美清纯卡通| h日本视频在线播放| 一级毛片电影观看 | 嫩草影院入口| 伦精品一区二区三区| 亚洲av.av天堂| 国产视频内射| 日韩人妻高清精品专区| 精品久久久久久久末码| 伊人久久精品亚洲午夜| 露出奶头的视频| 日韩欧美国产在线观看| eeuss影院久久| 熟女人妻精品中文字幕| 永久网站在线| 69av精品久久久久久| 精品久久久久久久人妻蜜臀av| 国产蜜桃级精品一区二区三区| 99热全是精品| 亚洲一区二区三区色噜噜| 亚洲七黄色美女视频| 亚洲丝袜综合中文字幕| 深夜精品福利| 欧美国产日韩亚洲一区| eeuss影院久久| 国产精品美女特级片免费视频播放器| 免费观看在线日韩| 欧美成人a在线观看| 亚洲中文日韩欧美视频| 亚洲内射少妇av| 亚洲第一电影网av| 2021天堂中文幕一二区在线观| 夜夜爽天天搞| 久久久久国产精品人妻aⅴ院| 久久久精品94久久精品| 国产aⅴ精品一区二区三区波| 亚洲美女黄片视频| 99久国产av精品国产电影| 看免费成人av毛片| 日韩成人av中文字幕在线观看 | 日韩成人av中文字幕在线观看 | 国产高清有码在线观看视频| 中国美白少妇内射xxxbb| 精品午夜福利在线看| 如何舔出高潮| 国产真实伦视频高清在线观看| 国产淫片久久久久久久久| 免费av不卡在线播放| 超碰av人人做人人爽久久| 欧美最黄视频在线播放免费| 日韩欧美三级三区| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品日韩av片在线观看| 美女高潮的动态| 国内精品美女久久久久久| 色播亚洲综合网| 日日摸夜夜添夜夜添av毛片| 露出奶头的视频| 精品久久久噜噜| 国内揄拍国产精品人妻在线| 中出人妻视频一区二区| 精品人妻视频免费看| 最近2019中文字幕mv第一页| av在线蜜桃| 亚洲五月天丁香| av在线亚洲专区| 亚洲内射少妇av| 国产综合懂色| 久久精品国产清高在天天线| 久久久久久久久久久丰满| 成人三级黄色视频| 成人二区视频| 亚洲美女黄片视频| 最新在线观看一区二区三区| 神马国产精品三级电影在线观看| 国产黄色视频一区二区在线观看 | 黄色日韩在线| av卡一久久| 欧美绝顶高潮抽搐喷水| 欧美日韩乱码在线| 国产精品99久久久久久久久| 久久综合国产亚洲精品| 亚洲成av人片在线播放无| 哪里可以看免费的av片| 真实男女啪啪啪动态图| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 久久精品夜夜夜夜夜久久蜜豆| 女人被狂操c到高潮| 97碰自拍视频| 国产欧美日韩一区二区精品| 日本-黄色视频高清免费观看| 一个人观看的视频www高清免费观看| 久久久久久久久久成人| 日本免费a在线| 赤兔流量卡办理| 国产av不卡久久| 真实男女啪啪啪动态图| 亚州av有码| 99久久成人亚洲精品观看| eeuss影院久久| 免费看av在线观看网站| 欧美日韩在线观看h| 韩国av在线不卡| 午夜福利高清视频| 国产麻豆成人av免费视频| 色哟哟哟哟哟哟| 99riav亚洲国产免费| 国内精品久久久久精免费| 美女大奶头视频| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区色噜噜| 高清毛片免费观看视频网站| 亚洲国产高清在线一区二区三| 久久久久九九精品影院| 久久久久国产精品人妻aⅴ院| 综合色丁香网| 国产精品99久久久久久久久| 亚洲av第一区精品v没综合| 日日摸夜夜添夜夜添av毛片| 国产成人影院久久av| 欧美最新免费一区二区三区| 在线观看午夜福利视频| 插逼视频在线观看| 给我免费播放毛片高清在线观看| 日本成人三级电影网站| 一本久久中文字幕| 不卡一级毛片| av女优亚洲男人天堂| 久久精品91蜜桃| 亚洲精品国产av成人精品 | 国产精品一区www在线观看| 寂寞人妻少妇视频99o| 国产亚洲精品久久久com| 亚洲中文字幕日韩| 99久久精品国产国产毛片| 变态另类成人亚洲欧美熟女| 国模一区二区三区四区视频| 欧美区成人在线视频| 18禁在线无遮挡免费观看视频 | 麻豆成人午夜福利视频| 国产乱人视频| 非洲黑人性xxxx精品又粗又长| 日韩av不卡免费在线播放| 亚洲成人久久爱视频| 日韩欧美 国产精品| 国产高清视频在线观看网站| 国产精品电影一区二区三区| 九九热线精品视视频播放| 国产av麻豆久久久久久久| 国产成人福利小说| 欧美高清成人免费视频www| 久久久久久大精品| 中国美白少妇内射xxxbb| 成人漫画全彩无遮挡| 91在线观看av| 超碰av人人做人人爽久久| 伦精品一区二区三区| 一个人免费在线观看电影| 国产伦一二天堂av在线观看| 久久6这里有精品| 欧美不卡视频在线免费观看| 欧美成人一区二区免费高清观看| 欧美三级亚洲精品| 国产人妻一区二区三区在| 久久热精品热| 成年免费大片在线观看| 日韩中字成人| 99久久成人亚洲精品观看| 久久精品国产清高在天天线| 人妻丰满熟妇av一区二区三区| 哪里可以看免费的av片| 国产精品人妻久久久久久| 国产精品精品国产色婷婷| 亚洲性久久影院| 少妇裸体淫交视频免费看高清| 两性午夜刺激爽爽歪歪视频在线观看| 女人十人毛片免费观看3o分钟| 久久久久国产精品人妻aⅴ院| 亚洲成av人片在线播放无| 国产一区二区在线av高清观看| 又黄又爽又刺激的免费视频.| 一级a爱片免费观看的视频| 99热6这里只有精品| 久久精品人妻少妇| 精品欧美国产一区二区三| 一个人看视频在线观看www免费| 69人妻影院| 国产精品一区二区三区四区免费观看 | 欧美国产日韩亚洲一区| 99热只有精品国产| 欧美一级a爱片免费观看看| 亚洲内射少妇av| 午夜免费男女啪啪视频观看 | 欧美最新免费一区二区三区| 毛片女人毛片| 日韩精品青青久久久久久| a级毛色黄片| 97人妻精品一区二区三区麻豆| 成熟少妇高潮喷水视频| 香蕉av资源在线| 亚洲在线观看片| 女生性感内裤真人,穿戴方法视频| 97超碰精品成人国产| 天美传媒精品一区二区| 国产白丝娇喘喷水9色精品| 最后的刺客免费高清国语| 国产一级毛片七仙女欲春2| 亚洲成人久久爱视频| 蜜臀久久99精品久久宅男| 久久精品国产99精品国产亚洲性色| 五月伊人婷婷丁香| 老师上课跳d突然被开到最大视频| 国产高清视频在线观看网站| 18+在线观看网站| 国产黄色小视频在线观看| 亚洲国产精品合色在线| 亚洲精品国产av成人精品 | 中出人妻视频一区二区| 免费av毛片视频| 国产高清三级在线| 国产精品爽爽va在线观看网站| 欧美激情国产日韩精品一区| 美女内射精品一级片tv| 永久网站在线| 97超级碰碰碰精品色视频在线观看| 亚洲成人中文字幕在线播放| 天堂影院成人在线观看| 少妇的逼好多水| 日韩欧美在线乱码| 91久久精品国产一区二区成人| 99视频精品全部免费 在线| 久久午夜亚洲精品久久| 婷婷精品国产亚洲av在线| 极品教师在线视频| 啦啦啦观看免费观看视频高清| 国产私拍福利视频在线观看| 亚洲av中文字字幕乱码综合| 亚洲七黄色美女视频| 三级男女做爰猛烈吃奶摸视频| 国产亚洲av嫩草精品影院| 男女啪啪激烈高潮av片| 亚洲av不卡在线观看| 国产私拍福利视频在线观看| 成人精品一区二区免费| 69av精品久久久久久| av.在线天堂| 51国产日韩欧美| 一级a爱片免费观看的视频| 日韩高清综合在线| 色哟哟哟哟哟哟| 99热这里只有是精品在线观看| 少妇丰满av| 国产在线精品亚洲第一网站| 成人特级av手机在线观看| 亚洲最大成人中文| 毛片女人毛片| 2021天堂中文幕一二区在线观| 成年av动漫网址| 国产成年人精品一区二区| 麻豆一二三区av精品| 亚洲七黄色美女视频| 国产亚洲精品久久久久久毛片| 精品人妻熟女av久视频| 婷婷色综合大香蕉| 精品一区二区三区av网在线观看| 深爱激情五月婷婷| av福利片在线观看| 国产麻豆成人av免费视频| 18禁黄网站禁片免费观看直播| 天堂影院成人在线观看| a级毛片a级免费在线| 国产精品一区二区三区四区免费观看 | 不卡视频在线观看欧美| 亚洲精品国产av成人精品 | 亚洲图色成人| 精品久久国产蜜桃| 大香蕉久久网| 在线免费十八禁| 亚洲精品日韩av片在线观看| 亚洲av成人精品一区久久| 黄色视频,在线免费观看| 一个人看的www免费观看视频| 麻豆成人午夜福利视频| www日本黄色视频网| 精品久久久久久久久久免费视频| 亚洲性夜色夜夜综合| 亚洲精品一卡2卡三卡4卡5卡| 午夜爱爱视频在线播放| 国产不卡一卡二| 啦啦啦观看免费观看视频高清| 少妇裸体淫交视频免费看高清| 在线观看免费视频日本深夜| 免费在线观看成人毛片| 久久这里只有精品中国| 亚洲丝袜综合中文字幕| 欧美在线一区亚洲| 国产免费一级a男人的天堂| 久久精品国产99精品国产亚洲性色| 国产精品免费一区二区三区在线| 日日摸夜夜添夜夜添av毛片| 亚洲va在线va天堂va国产| 亚洲精品日韩在线中文字幕 | 在线a可以看的网站| 能在线免费观看的黄片| 在线观看一区二区三区| 国产精品野战在线观看| 天天躁日日操中文字幕| 日韩精品青青久久久久久| 国产久久久一区二区三区| 九色成人免费人妻av| 不卡视频在线观看欧美| 在现免费观看毛片| 成人三级黄色视频| 人妻丰满熟妇av一区二区三区| 校园人妻丝袜中文字幕| 精品一区二区三区视频在线观看免费| 久久精品国产99精品国产亚洲性色| 亚洲av美国av| 熟妇人妻久久中文字幕3abv| 男人的好看免费观看在线视频| 日本五十路高清| 国产在视频线在精品| 听说在线观看完整版免费高清| 国产亚洲精品久久久久久毛片| 日本黄色片子视频| 国产精品人妻久久久影院| 国产伦一二天堂av在线观看| 国内精品宾馆在线| 欧美日韩精品成人综合77777| 精品国内亚洲2022精品成人| 亚洲欧美日韩卡通动漫| 国产高清视频在线播放一区| 欧美极品一区二区三区四区| 国产精品久久视频播放| 男女啪啪激烈高潮av片| 又爽又黄a免费视频| 成人漫画全彩无遮挡| av.在线天堂| 亚洲精品456在线播放app| 国产亚洲91精品色在线| 少妇熟女aⅴ在线视频| 亚洲国产精品sss在线观看| 亚洲在线自拍视频| 亚洲av熟女| 大型黄色视频在线免费观看| 国产av一区在线观看免费| 欧美性猛交╳xxx乱大交人| 中文亚洲av片在线观看爽| 夜夜夜夜夜久久久久| 日本精品一区二区三区蜜桃| 黄色一级大片看看| 久久精品夜色国产| 国产精品,欧美在线| 亚洲性久久影院| 好男人在线观看高清免费视频| 18禁黄网站禁片免费观看直播| 可以在线观看毛片的网站| 美女内射精品一级片tv| 国产三级在线视频| 欧美一级a爱片免费观看看| 69av精品久久久久久| 日韩精品青青久久久久久| 国产综合懂色| 日日干狠狠操夜夜爽| 国产乱人视频| 国产综合懂色| av在线播放精品| 日韩一本色道免费dvd| 亚洲人成网站在线播| 精品久久久久久久久亚洲| 国产日本99.免费观看| 亚洲国产高清在线一区二区三| 精品人妻熟女av久视频| 色在线成人网| 亚洲精品乱码久久久v下载方式| 日日撸夜夜添| 中文字幕av在线有码专区| 天美传媒精品一区二区| 欧美日韩综合久久久久久| 99在线视频只有这里精品首页| 婷婷精品国产亚洲av在线| 特大巨黑吊av在线直播| 免费电影在线观看免费观看| 国产欧美日韩精品一区二区| 亚洲av一区综合| 中文亚洲av片在线观看爽| 三级毛片av免费| 久久精品久久久久久噜噜老黄 | 精品乱码久久久久久99久播| 少妇的逼水好多| 女人十人毛片免费观看3o分钟| 亚洲精品影视一区二区三区av| 久久精品国产亚洲网站| 亚洲,欧美,日韩| 在线免费十八禁| 91麻豆精品激情在线观看国产| 搞女人的毛片| 午夜爱爱视频在线播放| 99久久中文字幕三级久久日本|