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

    不同統(tǒng)計(jì)方法在手足口病季節(jié)性分析中的比較研究*

    2015-01-27 10:31:04唐咸艷歐陽(yáng)云仇小強(qiáng)黃高明周紅霞
    關(guān)鍵詞:季節(jié)性南寧市圓形

    唐咸艷 周 穎 李 嶠 歐陽(yáng)云 仇小強(qiáng) 黃高明 周紅霞

    廣西醫(yī)科大學(xué)公共衛(wèi)生學(xué)院流行病與衛(wèi)生統(tǒng)計(jì)學(xué)教研室(530021)

    不同統(tǒng)計(jì)方法在手足口病季節(jié)性分析中的比較研究*

    唐咸艷 周 穎 李 嶠 歐陽(yáng)云 仇小強(qiáng) 黃高明 周紅霞△

    廣西醫(yī)科大學(xué)公共衛(wèi)生學(xué)院流行病與衛(wèi)生統(tǒng)計(jì)學(xué)教研室(530021)

    目的 比較不同統(tǒng)計(jì)方法在手足口病季節(jié)性研究中的應(yīng)用價(jià)值,為手足口病季節(jié)性防控提供科學(xué)依據(jù)。方法 以南寧市手足口病為例,收集2006-2009年南寧市各鄉(xiāng)鎮(zhèn)、城區(qū)手足口病的新發(fā)病例數(shù),利用流行曲線、季節(jié)指數(shù)、變異系數(shù)、圓形分布法以及掃描統(tǒng)計(jì)量等方法分析手足口病的季節(jié)性趨勢(shì)。結(jié)果 流行曲線直觀描述南寧市手足口病高發(fā)月份集中在4-7月,12月和1月呈現(xiàn)小高峰。季節(jié)指數(shù)判斷1月、4-7月為南寧市手足口病的高發(fā)季節(jié),以5-6月份發(fā)生強(qiáng)度最大。圓形分布顯示南寧市手足口病的高發(fā)日期為5月10日,高發(fā)時(shí)段跨度大。變異系數(shù)顯示5-6月和9月的變異程度較低,為潛在的季節(jié)性月份。掃描統(tǒng)計(jì)量顯示南寧市手足口病高發(fā)時(shí)段集中在4月下旬到7月上旬之間,以5-6月份最為高發(fā);高發(fā)時(shí)段的相對(duì)風(fēng)險(xiǎn)為非高發(fā)時(shí)段的1.33~1.85倍。結(jié)論 不同統(tǒng)計(jì)方法在檢測(cè)手足口病季節(jié)性規(guī)律時(shí)的側(cè)重面、精確性和實(shí)際意義有所不同,探測(cè)到高發(fā)時(shí)段的價(jià)值和意義亦有所不同。變異系數(shù)描述疾病在某時(shí)段內(nèi)流行強(qiáng)度的變異度,變異度高說(shuō)明疾病流行強(qiáng)度的強(qiáng)弱交替,但變異度低可能意味著普遍高發(fā)或普遍低發(fā)。流行曲線、季節(jié)指數(shù)雖可直觀描述手足口病的季節(jié)趨勢(shì),但缺乏統(tǒng)計(jì)推斷,可能受隨機(jī)誤差的影響。圓形分布可定量探測(cè)出具有統(tǒng)計(jì)學(xué)意義的手足口病高發(fā)日期和高發(fā)時(shí)段,但可能因高發(fā)時(shí)段跨度較大,而失去實(shí)際意義。掃描統(tǒng)計(jì)量兼顧手足口病的時(shí)空信息,在校正單純性時(shí)間變化趨勢(shì)和周內(nèi)效應(yīng)后,可探測(cè)到更精確、穩(wěn)定和更具實(shí)際意義的高發(fā)時(shí)段,在探討手足口病的季節(jié)趨勢(shì)中具有優(yōu)勢(shì)。

    手足口病 圓形分布法 掃描統(tǒng)計(jì)量 季節(jié)分析

    手足口病呈現(xiàn)一定的季節(jié)性分布,但各地區(qū)間的高發(fā)期不盡相同[1-2]。手足口病在地區(qū)間的季節(jié)性差異可能是其傳播機(jī)制與當(dāng)?shù)厣鷳B(tài)因素相互作用的結(jié)果,也可能是探討手足口病季節(jié)性所用的方法不同而引起的差異。因此,運(yùn)用合理的統(tǒng)計(jì)方法,及時(shí)查明當(dāng)?shù)厥肿憧诓〉募竟?jié)性特征,是科學(xué)防控手足口病的重要前提。本研究以南寧市手足口病為例,旨在探討傳統(tǒng)統(tǒng)計(jì)方法(流行曲線、季節(jié)指數(shù)、變異系數(shù)、圓形分布法),以及掃描統(tǒng)計(jì)量在手足口病季節(jié)性分析中的價(jià)值,為手足口病的季節(jié)性防控提供依據(jù)。

    資料與方法

    1.資料來(lái)源

    2006年1月1日-2009年12月31日逐日南寧市各鄉(xiāng)鎮(zhèn)(城區(qū))手足口病發(fā)病資料來(lái)源于國(guó)家法定傳染病網(wǎng)絡(luò)直報(bào)系統(tǒng),人口數(shù)據(jù)來(lái)源于廣西統(tǒng)計(jì)年鑒。通過(guò)發(fā)病數(shù)和人口數(shù)計(jì)算出南寧市各鄉(xiāng)鎮(zhèn)、城區(qū)手足口病的發(fā)病率(1/10萬(wàn))。

    2.統(tǒng)計(jì)方法

    (1)流行曲線:直觀描述疾病的發(fā)生水平隨時(shí)間的變化趨勢(shì)。

    (2)季節(jié)指數(shù)

    季節(jié)指數(shù)是時(shí)間數(shù)列中將一定時(shí)期內(nèi)因季節(jié)變動(dòng)而出現(xiàn)的有規(guī)律的周期性變動(dòng)指標(biāo)進(jìn)行分析[3]。季節(jié)指數(shù)如果大于100%表示旺季,小于100%表示淡季,在100%左右表示平季[4]。

    (3)圓形分布法

    圓形分布法首先將周期性時(shí)間資料轉(zhuǎn)換為圓的角度,通過(guò)計(jì)算平均角及平均角的標(biāo)準(zhǔn)差,對(duì)手足口病發(fā)生的季節(jié)性進(jìn)行統(tǒng)計(jì)描述。此外,通過(guò)對(duì)平均角進(jìn)行假設(shè)檢驗(yàn),判斷手足口病的時(shí)間分布是否具有季節(jié)性[5]。

    (4)變異系數(shù)

    變異系數(shù)反映數(shù)據(jù)變異程度,其大小不僅受變量值離散程度的影響,還受變量值平均水平大小的影響。變異系數(shù)是標(biāo)準(zhǔn)差和均數(shù)之比,用百分?jǐn)?shù)表示,即CV=(標(biāo)準(zhǔn)差/平均值)×100%。各年間某月份手足口病的變異系數(shù)越小,說(shuō)明該月份具有季節(jié)性,可能是季節(jié)性高發(fā)月份,也可能是季節(jié)性低發(fā)月份。若正確判定,需要結(jié)合流行曲線。

    (5)掃描統(tǒng)計(jì)量

    建立活動(dòng)性掃描窗口,窗口的大小和位置處于動(dòng)態(tài)的改變之中。窗口的圓心根據(jù)地圖中網(wǎng)格線或地理單位中心變動(dòng);窗口的半徑在0至總?cè)丝谝欢ū壤脑O(shè)定值之間變動(dòng)[6]。掃描窗口半徑每變動(dòng)一次,通過(guò)統(tǒng)計(jì)量LLR(log likelihood ratio,LLR)計(jì)算窗口內(nèi)外區(qū)域間發(fā)病率的差異大小,最終選取LLR最大的掃描窗口作為時(shí)空聚集區(qū)。

    3.軟件實(shí)現(xiàn)

    手足口病的流行曲線、季節(jié)指數(shù)、變異系數(shù)和圓形分布法通過(guò)excel 2007實(shí)現(xiàn);掃描統(tǒng)計(jì)量在SatScan 9.2中實(shí)現(xiàn)。

    結(jié) 果

    1.南寧市手足口病的流行曲線

    2006-2009年南寧市手足口病流行曲線(圖1)呈“W”型,結(jié)果直觀顯示高發(fā)月份主要集中在4-7月;1月和12月又出現(xiàn)一個(gè)小高峰。

    2.南寧市手足口病的季節(jié)指數(shù)

    南寧市手足口病1-12月份的季節(jié)指數(shù)見(jiàn)表1。結(jié)果說(shuō)明,1月、4-7月為南寧市手足口病的高發(fā)月份,其中5、6月份發(fā)生強(qiáng)度較大。

    3.南寧市手足口病的圓形分布法

    南寧市手足口病的圓形分布法分析結(jié)果見(jiàn)表2:sinα=0.7901,cosα=-0.6128,因此平均角α=127.7968°,落在第二象限內(nèi);α對(duì)應(yīng)的日期為5月10日,故南寧市手足口病的發(fā)病高峰時(shí)點(diǎn)為5月10日。γ值較小,s值較大,說(shuō)明季節(jié)性高峰平緩,高發(fā)時(shí)段分散太長(zhǎng)(α±s)對(duì)應(yīng)于1月11日-9月5日。Z=218.4525,P<0.05,說(shuō)明南寧市手足口病的季節(jié)性差異具有統(tǒng)計(jì)學(xué)意義。

    4.南寧市手足口病的變異系數(shù)

    2006-2009年南寧市手足口病發(fā)病的變異系數(shù)見(jiàn)表3。結(jié)果顯示:5-6月,以及9月的CV分別為11.84%、21.83%和16.17%,這說(shuō)明5-6月和9月的離散程度較小,具有較明顯的季節(jié)性高發(fā)或低發(fā)。結(jié)合流行曲線圖,可知5-6月份為手足口病的高發(fā)月份。

    5.南寧市手足口病的掃描統(tǒng)計(jì)量

    2006-2009年南寧市手足口病的單純性時(shí)間掃描統(tǒng)計(jì)量見(jiàn)表4。盡管各年的高發(fā)期不盡相同,但較相似。高發(fā)時(shí)段大都集中在4月下旬到7月上旬之間,發(fā)生風(fēng)險(xiǎn)位于1.33~1.85之間。

    討 論

    目前,探討手足口病季節(jié)性的方法不僅有傳統(tǒng)的統(tǒng)計(jì)描述、統(tǒng)計(jì)推斷法,還有時(shí)空分析技術(shù)。季節(jié)性分析中,季節(jié)性研究常用的統(tǒng)計(jì)描述法有流行曲線、變異系數(shù)、季節(jié)指數(shù)、集中度(M值)法[2-4]。這些方法主要用圖表和統(tǒng)計(jì)指標(biāo)揭示有無(wú)季節(jié)性趨勢(shì),但并不進(jìn)行統(tǒng)計(jì)推斷,因此研究結(jié)果可能受到隨機(jī)誤差的影響。圓形分布法和游程檢驗(yàn)進(jìn)行季節(jié)性檢驗(yàn)[5,7],可消除隨機(jī)誤差的干擾。

    本研究中,流行曲線直觀揭示南寧市手足口病的流行高峰期為4-7月,其次12月和1月也有小的升高。但是,這種高峰期的鑒定缺乏統(tǒng)計(jì)量化標(biāo)準(zhǔn),可能存在“假象”。本研究中流行曲線的比例為相對(duì)構(gòu)成比,其數(shù)值大小不僅受本月手足口病發(fā)生例數(shù)大小的影響,還可能受到其他月份手足口病發(fā)生例數(shù)變化的影響。因此,這種高峰不穩(wěn)定,可能會(huì)呈現(xiàn)虛假高發(fā)的現(xiàn)象。季節(jié)指數(shù)用比值大小量化手足口病各月份發(fā)生的強(qiáng)度大小,但也存在一定的局限性。一方面,季節(jié)指數(shù)僅僅是描述,缺乏統(tǒng)計(jì)推斷,其結(jié)果可能受到隨機(jī)誤差的干擾;另一方面,季節(jié)指數(shù)以月為單位,時(shí)間尺度并不精確,可能會(huì)帶來(lái)高估或低估的問(wèn)題。對(duì)某些月份而言,其高發(fā)日期僅從下旬開始,若季節(jié)指數(shù)把該月視為高發(fā)月份的話,那么該月上旬就被錯(cuò)誤地視為高發(fā)日期了。變異系數(shù)雖可判斷出變異程度較低的月份,但變異程度低可能意味著普遍高發(fā)或普遍低發(fā)。因此,變異系數(shù)自身無(wú)法定性判定是高發(fā)月份還是低發(fā)月份,故常需與流行曲線或季節(jié)指數(shù)聯(lián)合使用。本研究中,變異系數(shù)探測(cè)出變異程度較低的月份有5月、6月和9月,但這并不代表這三個(gè)月都是手足口病的高發(fā)月份。實(shí)際上,9月的發(fā)生強(qiáng)度并不高,約只占全年的5%左右,因此9月份應(yīng)當(dāng)屬于低發(fā)月份。與前面的季節(jié)性方法相比,圓形分布法不僅可探索出高發(fā)日期和高發(fā)時(shí)段,并對(duì)其進(jìn)行統(tǒng)計(jì)推斷,排除了隨機(jī)誤差的影響。但必須注意圓形分布法界定的高發(fā)時(shí)段的實(shí)際意義大小。當(dāng)圓形分布中的角標(biāo)準(zhǔn)差(s)非常大時(shí),這必將導(dǎo)致高發(fā)時(shí)段(α±s)變得非常寬廣;此時(shí),圓形分布界定的高發(fā)時(shí)段可能就失去了實(shí)際意義。

    掃描統(tǒng)計(jì)量是一種較為新穎的時(shí)空技術(shù),不僅可用來(lái)探測(cè)疾病的時(shí)空聚集性,還可用于疾病的早期預(yù)警研究。與傳統(tǒng)的季節(jié)性研究方法相比,掃描統(tǒng)計(jì)量兼顧手足口病的時(shí)空信息,在校正單純性的時(shí)間變化趨勢(shì)和周內(nèi)效應(yīng)后,掃描統(tǒng)計(jì)量通過(guò)對(duì)數(shù)似然比法LLR評(píng)價(jià)高發(fā)時(shí)段是否具有統(tǒng)計(jì)學(xué)意義,并給出相對(duì)危險(xiǎn)度(relative risk,RR)。此外,掃描統(tǒng)計(jì)通過(guò)蒙特卡洛法(Monte Carlo)進(jìn)行模擬,探索到的高發(fā)時(shí)段較圓形分布法可能更精確、穩(wěn)定和更具實(shí)際意義。本研究中,在校正時(shí)間趨勢(shì)和周內(nèi)變化趨勢(shì)后,掃描統(tǒng)計(jì)量探測(cè)到的手足口病高發(fā)時(shí)段穩(wěn)定在4月下旬至7月上旬之間(較圓形分布和季節(jié)指數(shù)的高發(fā)時(shí)段精確),并以5月和6月發(fā)生強(qiáng)度最大;高發(fā)時(shí)段的相對(duì)風(fēng)險(xiǎn)為非高發(fā)時(shí)段的1.33~1.85倍之間。

    手足口病的季節(jié)性在地區(qū)間有所不同,這可能與當(dāng)?shù)氐臍庀髼l件有一定的關(guān)系。北方地區(qū),手足口病的高發(fā)月份集中在6-8月份[1];而南方地區(qū),其高發(fā)期始于4月,止于9月,并可能在1月和12月出現(xiàn)小的高峰[2]。概括而言,因腸道病毒受氣溫和相對(duì)濕度的影響,手足口病的高發(fā)月份在南方地區(qū)出現(xiàn)早、終止晚,流行期相對(duì)較長(zhǎng);而北方地區(qū),其高發(fā)期來(lái)的相對(duì)遲、終止的也早,且不易在冬春季節(jié)出現(xiàn)小的高峰。本研究中,南寧地區(qū)4-5月正值梅雨季節(jié),氣溫高、濕度大,非常適合腸道病毒的繁殖。這也導(dǎo)致了手足口病的高發(fā)期在南寧地區(qū)出現(xiàn)相對(duì)較早(4月中下旬)。值得注意的是2008年,其流行情況與歷年不同:發(fā)生高峰期始于5月初,較別的年份有所遲緩;而且在1月份,沒(méi)出現(xiàn)小的高峰。這可能與南寧地區(qū)在2008年1月遭遇罕見(jiàn)冰凍災(zāi)害有關(guān),氣溫驟降,腸道病毒不宜孳生。罕見(jiàn)冰凍氣象災(zāi)害的出現(xiàn),也在一定程度上延遲了4月氣溫的回升。

    當(dāng)然,本研究也存在一定的局限性。首先,文中所探討的方法僅是季節(jié)性常用的一些方法,集中度、時(shí)間序列等其他方法并沒(méi)涉及,因此研究結(jié)果可能不夠全面。其次,季節(jié)性是手足口病的重要特性,在短時(shí)期內(nèi),季節(jié)性出現(xiàn)較大的短期波動(dòng)的可能性不大;若從長(zhǎng)期趨勢(shì)來(lái)看,全球氣候的變暖,以及兒童健康行為的改變,可能會(huì)導(dǎo)致手足口病的季節(jié)性出現(xiàn)一定的變動(dòng)。本研究探討2006-2009年的手足口病資料,沒(méi)及時(shí)納入2011-2012年的病例資料,這可能會(huì)限制我們探究手足口病季節(jié)性的短期波動(dòng)。最后,本研究推測(cè)氣象條件可能與手足口病季節(jié)性的形成及其波動(dòng)有關(guān),但缺乏數(shù)據(jù)支持。因此后續(xù)研究中,我們擬采用空間統(tǒng)計(jì)學(xué),從時(shí)空維度和鄉(xiāng)鎮(zhèn)級(jí)小尺度上探討氣象條件與手足口病的關(guān)系。

    [1]白婧,王菊光,劉方遙,等.2011年北京市海淀區(qū)手足口病實(shí)驗(yàn)室檢測(cè)結(jié)果分析.中國(guó)衛(wèi)生統(tǒng)計(jì),2013,30(1):95-97.

    [2]曾舸,高立冬,黃威,等.湖南省2008-2010年手足口病的流行病學(xué)分析.現(xiàn)代預(yù)防醫(yī)學(xué),2013,40(6):1137-1139.

    [3]解華山,劉凌云.荊門市常見(jiàn)傳染病季節(jié)性分布特征動(dòng)態(tài)分析.公共衛(wèi)生與預(yù)防醫(yī)學(xué),2011,22(2):45-48.

    [4]曹世義,尹平,盧祖洵.應(yīng)用改良的集中度法分析我國(guó)結(jié)核病發(fā)病的季節(jié)性特征.中華疾病控制雜志,2011,15(12):1081-1083.

    [5]李艷琳,周雨.應(yīng)用圓形分布分析我區(qū)近年甲型肝炎患者發(fā)病的季節(jié)性趨勢(shì).中國(guó)衛(wèi)生統(tǒng)計(jì),2010,27(3):303-304.

    [6]方立群,曹春香,陳國(guó)勝,等.地理信息系統(tǒng)應(yīng)用于中國(guó)大陸高致病性禽流感的空間分布及環(huán)境因素分析.中華流行病學(xué)雜志,2005,26(11):839-842.

    [7]許陽(yáng)婷,洪鐳.南京市2001-2011年16種傳染病時(shí)間聚集性分布特征動(dòng)態(tài)分析.現(xiàn)代預(yù)防醫(yī)學(xué),2013,40(18):3346-3349.

    (責(zé)任編輯:丁海龍)

    Comparison of Statistical Methods for Detecting Seasonality of HFMD

    Tang Xianyan,Zhou Ying,Li Qiao,et al.

    (Department of Epidemiology and Biostatistics,School of Public Health,Guangxi Medical University(530021),Nanning)

    Objective To compare the utility of statistical methods for detecting seasonality of HFMD.Methods Daily data on notified cases of HFMD over the period of 2006-2009 was collected from Nanning prefecture.Epidemic curve,seasonality index,coefficient of variance,circular distribution and scan statistics were employed to detect the seasonality of HFMD.Results Period of April to July was found to be the months with higher HFMD incidence by epidemic curve,as well as December & January with secondary incidence peak.Seasonality index found that January and April to July were the months of higher HFMD incidence,May and June were the incidence peak of HFMD as well.Circular distribution showed that 10th May was the peak onset of HFMD,and the duration of higher HFMD incidence was over 11th January to 5th September.May and June were the potential months with higher HFMD incidence,as well as lowest coefficient of variance.And period of late April to early July was detected as higher HFMD incidence cluster by scan statistics,with relative risk of 1.33~1.85.Conclusion As far as various statistical methods were concerned,the capacity of detecting seasonality was inconsistent.Coefficient of variance just showed the magnitude of variance of HFMD incidence,rather than identifying months with higher incidence of HFMD.Considering the shortage of statistic inference,the seasonality of HFMD detected by seasonality index would be disturbed by random error,although it could statistically describe the seasonal trend of HFMD incidence.To some extent,the capacity of circular distribution was determined by parameter,and the outcome of detecting seasonality might be impractical.In contrast,scan statistics might be more powerful,practical and robust to detect seasonality of HFMD,by integrating spatiotemporal information of HFMD and adjusting for effect of day-of-week,as well as effect of purely temporal trend.

    HFMD;Circular distribution;Scan statistics;Seasonality analysis

    廣西科學(xué)研究與技術(shù)開發(fā)計(jì)劃項(xiàng)目(桂科攻10124001A-65);廣西自然科學(xué)基金(2013GXNSFBA019125);廣西研究生教育創(chuàng)新計(jì)劃資助項(xiàng)目(YCSZ2013029)

    △通信作者:周紅霞,E-mail:gmuies@163.com

    猜你喜歡
    季節(jié)性南寧市圓形
    粕類季節(jié)性規(guī)律:豆粕篇
    湖南飼料(2021年3期)2021-07-28 07:05:58
    Sam’s Family
    季節(jié)性需求放緩 鉀肥價(jià)格下行
    關(guān) 公
    寶藏(2019年3期)2019-03-28 05:24:16
    蔬菜價(jià)格呈季節(jié)性回落
    遠(yuǎn)離季節(jié)性過(guò)敏
    Coco薇(2017年12期)2018-01-03 21:34:42
    為什么窨井蓋大多都是圓形的
    肥皂泡為什么是圓形?
    圓形題
    圓形變身喵星人
    老鸭窝网址在线观看| 国产精品蜜桃在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费播放大片免费观看视频在线观看| 久久久久久久精品精品| 97精品久久久久久久久久精品| 午夜福利视频在线观看免费| 免费女性裸体啪啪无遮挡网站| 日韩中字成人| 午夜免费鲁丝| 中文字幕亚洲精品专区| 欧美老熟妇乱子伦牲交| 国产精品香港三级国产av潘金莲 | 丰满少妇做爰视频| videos熟女内射| 99re6热这里在线精品视频| 国产黄色免费在线视频| 9色porny在线观看| 如何舔出高潮| 岛国毛片在线播放| 成人国产麻豆网| 久久精品久久精品一区二区三区| 国产精品国产三级国产专区5o| 久久久久久久亚洲中文字幕| 最近最新中文字幕大全免费视频 | 亚洲熟女精品中文字幕| 91在线精品国自产拍蜜月| 麻豆精品久久久久久蜜桃| 人妻 亚洲 视频| 最近最新中文字幕大全免费视频 | 亚洲av成人精品一二三区| 国产av国产精品国产| 亚洲成人一二三区av| 欧美老熟妇乱子伦牲交| 婷婷色综合大香蕉| 丝袜美足系列| 国产日韩欧美在线精品| 午夜福利网站1000一区二区三区| www.熟女人妻精品国产| 欧美成人午夜精品| 欧美日本中文国产一区发布| 国产麻豆69| 91精品国产国语对白视频| 黄色 视频免费看| 伦精品一区二区三区| 91国产中文字幕| 成人漫画全彩无遮挡| 99热全是精品| 日韩,欧美,国产一区二区三区| 国产一区二区三区av在线| 精品第一国产精品| 狠狠精品人妻久久久久久综合| 免费大片黄手机在线观看| 国产 一区精品| 在线观看美女被高潮喷水网站| 天堂俺去俺来也www色官网| 久久亚洲国产成人精品v| 视频在线观看一区二区三区| 综合色丁香网| av线在线观看网站| 亚洲av电影在线进入| 国产麻豆69| 一级毛片我不卡| 亚洲伊人色综图| 国产高清不卡午夜福利| 欧美日韩亚洲国产一区二区在线观看 | 满18在线观看网站| 欧美精品国产亚洲| 五月天丁香电影| 久久精品熟女亚洲av麻豆精品| 亚洲av日韩在线播放| 高清av免费在线| 国产精品99久久99久久久不卡 | 亚洲av中文av极速乱| 丰满乱子伦码专区| 国产精品久久久久久久久免| 亚洲成人av在线免费| 性少妇av在线| 国产 精品1| 少妇被粗大的猛进出69影院| 69精品国产乱码久久久| av国产久精品久网站免费入址| 亚洲精品一区蜜桃| 精品国产露脸久久av麻豆| 日韩大片免费观看网站| 狂野欧美激情性bbbbbb| 人妻一区二区av| 国产熟女午夜一区二区三区| 婷婷色综合www| 国产又爽黄色视频| 亚洲国产精品一区二区三区在线| 多毛熟女@视频| 少妇熟女欧美另类| 日韩,欧美,国产一区二区三区| 国产国语露脸激情在线看| 久久午夜福利片| 精品亚洲成国产av| 日韩一本色道免费dvd| 91久久精品国产一区二区三区| 欧美成人午夜精品| 性高湖久久久久久久久免费观看| 亚洲精品乱久久久久久| 亚洲精品第二区| 中文字幕人妻丝袜一区二区 | 卡戴珊不雅视频在线播放| 精品亚洲乱码少妇综合久久| 男女下面插进去视频免费观看| 一级片免费观看大全| 亚洲内射少妇av| 一级黄片播放器| 欧美+日韩+精品| 国产极品粉嫩免费观看在线| 欧美 亚洲 国产 日韩一| 欧美av亚洲av综合av国产av | √禁漫天堂资源中文www| 日韩 亚洲 欧美在线| kizo精华| 亚洲国产av影院在线观看| 高清欧美精品videossex| 亚洲国产看品久久| 日本av免费视频播放| 国产1区2区3区精品| 亚洲国产欧美日韩在线播放| 免费观看在线日韩| 丰满乱子伦码专区| 亚洲美女搞黄在线观看| 寂寞人妻少妇视频99o| 亚洲欧洲日产国产| 久久久久久久大尺度免费视频| 久久精品人人爽人人爽视色| 精品午夜福利在线看| 少妇 在线观看| 亚洲男人天堂网一区| av福利片在线| 女性被躁到高潮视频| 赤兔流量卡办理| 中文欧美无线码| 高清不卡的av网站| 国产精品香港三级国产av潘金莲 | 边亲边吃奶的免费视频| 日本爱情动作片www.在线观看| 久久人人97超碰香蕉20202| 丝袜美腿诱惑在线| 精品第一国产精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产黄频视频在线观看| 欧美人与性动交α欧美软件| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 97人妻天天添夜夜摸| 人妻一区二区av| 搡女人真爽免费视频火全软件| av不卡在线播放| 欧美变态另类bdsm刘玥| 亚洲国产av影院在线观看| 亚洲精品日本国产第一区| 18+在线观看网站| 观看av在线不卡| 久久精品国产综合久久久| 亚洲三级黄色毛片| 免费播放大片免费观看视频在线观看| 一二三四在线观看免费中文在| 亚洲,欧美,日韩| 国产成人91sexporn| 国产片特级美女逼逼视频| 欧美激情 高清一区二区三区| av福利片在线| 亚洲欧美精品综合一区二区三区 | 欧美精品av麻豆av| 国产深夜福利视频在线观看| 亚洲人成电影观看| 亚洲欧美一区二区三区久久| 伦精品一区二区三区| 美女国产视频在线观看| 久久久久精品性色| 五月开心婷婷网| 中文字幕色久视频| 三上悠亚av全集在线观看| av.在线天堂| 国产高清国产精品国产三级| 久久久国产欧美日韩av| 如日韩欧美国产精品一区二区三区| 精品卡一卡二卡四卡免费| 久久婷婷青草| 少妇人妻久久综合中文| 在线观看www视频免费| 久久女婷五月综合色啪小说| 久久影院123| 自线自在国产av| 天天操日日干夜夜撸| 青青草视频在线视频观看| 日日啪夜夜爽| 免费观看a级毛片全部| 亚洲精品国产av成人精品| 国产成人午夜福利电影在线观看| 日韩精品有码人妻一区| 日韩一区二区三区影片| 国产熟女欧美一区二区| 国精品久久久久久国模美| 日本欧美视频一区| 日韩精品有码人妻一区| 天堂俺去俺来也www色官网| 99国产综合亚洲精品| 国产亚洲最大av| 久久久久久人妻| 久久免费观看电影| 日韩欧美精品免费久久| 久久亚洲国产成人精品v| 久久久久久伊人网av| 最近的中文字幕免费完整| 女人久久www免费人成看片| 国产成人精品婷婷| 韩国精品一区二区三区| 国产精品偷伦视频观看了| 一本久久精品| 女人久久www免费人成看片| 国产精品欧美亚洲77777| 伊人久久大香线蕉亚洲五| 电影成人av| 日日爽夜夜爽网站| 在线观看免费视频网站a站| 久久精品久久久久久噜噜老黄| 久久久久精品性色| 成人亚洲欧美一区二区av| 有码 亚洲区| 免费黄网站久久成人精品| 91午夜精品亚洲一区二区三区| 午夜福利视频精品| 欧美+日韩+精品| 亚洲国产欧美在线一区| 亚洲国产成人一精品久久久| 欧美亚洲 丝袜 人妻 在线| 色播在线永久视频| 亚洲,欧美精品.| 最近最新中文字幕免费大全7| 国产精品国产三级国产专区5o| 香蕉丝袜av| 日本爱情动作片www.在线观看| 日韩精品免费视频一区二区三区| 国产成人免费观看mmmm| 飞空精品影院首页| 欧美最新免费一区二区三区| 日本vs欧美在线观看视频| 成年人午夜在线观看视频| 美女高潮到喷水免费观看| 国产激情久久老熟女| 亚洲av在线观看美女高潮| 丝袜人妻中文字幕| 两个人看的免费小视频| 亚洲图色成人| 只有这里有精品99| 在线 av 中文字幕| av.在线天堂| 啦啦啦在线免费观看视频4| 午夜福利一区二区在线看| 日日撸夜夜添| 伊人久久大香线蕉亚洲五| 亚洲人成电影观看| 久久久久国产精品人妻一区二区| 最近最新中文字幕大全免费视频 | 9色porny在线观看| 老汉色∧v一级毛片| 亚洲av欧美aⅴ国产| 亚洲婷婷狠狠爱综合网| av有码第一页| 成人18禁高潮啪啪吃奶动态图| 一个人免费看片子| 七月丁香在线播放| av线在线观看网站| 日本色播在线视频| 国产精品一二三区在线看| 精品国产乱码久久久久久男人| 久久久久久久久免费视频了| 18禁国产床啪视频网站| 亚洲欧洲日产国产| 日韩精品免费视频一区二区三区| 在线观看三级黄色| 国产一区有黄有色的免费视频| 最新中文字幕久久久久| 国产一级毛片在线| 天天躁夜夜躁狠狠久久av| 丰满少妇做爰视频| 久久国产精品大桥未久av| 亚洲精品久久午夜乱码| 亚洲国产精品999| 亚洲av福利一区| 国产爽快片一区二区三区| 国产精品 国内视频| 最近手机中文字幕大全| 99国产综合亚洲精品| 一区二区三区四区激情视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产爽快片一区二区三区| 久久99精品国语久久久| 日韩一区二区三区影片| 18在线观看网站| 美女国产高潮福利片在线看| 成年美女黄网站色视频大全免费| 丰满少妇做爰视频| 最近中文字幕2019免费版| 一区二区三区四区激情视频| 亚洲国产日韩一区二区| 成人18禁高潮啪啪吃奶动态图| 在线观看免费视频网站a站| 亚洲精品乱久久久久久| 久久久国产欧美日韩av| 色婷婷av一区二区三区视频| 高清在线视频一区二区三区| 一区二区av电影网| 久久久久久久国产电影| 国产极品粉嫩免费观看在线| 亚洲精品国产av蜜桃| 亚洲美女黄色视频免费看| 国产精品免费视频内射| 国产精品 国内视频| 老司机亚洲免费影院| 国产av精品麻豆| 国产1区2区3区精品| 亚洲久久久国产精品| 国产片内射在线| 日韩视频在线欧美| 熟女少妇亚洲综合色aaa.| 国产野战对白在线观看| 18禁裸乳无遮挡动漫免费视频| 久久久久久人妻| 国产熟女午夜一区二区三区| 欧美日韩亚洲高清精品| 五月伊人婷婷丁香| 国产精品蜜桃在线观看| 26uuu在线亚洲综合色| 美女视频免费永久观看网站| 免费高清在线观看视频在线观看| 建设人人有责人人尽责人人享有的| 成人毛片a级毛片在线播放| 日韩免费高清中文字幕av| 午夜福利在线免费观看网站| 精品国产国语对白av| 国产成人精品福利久久| 久久韩国三级中文字幕| www.av在线官网国产| 日本爱情动作片www.在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 女人高潮潮喷娇喘18禁视频| av片东京热男人的天堂| 亚洲少妇的诱惑av| 少妇的丰满在线观看| 国产亚洲av片在线观看秒播厂| 亚洲欧洲国产日韩| 亚洲精品一区蜜桃| 最近最新中文字幕免费大全7| 久久精品国产鲁丝片午夜精品| 久久精品国产自在天天线| 亚洲av综合色区一区| 只有这里有精品99| 老汉色∧v一级毛片| 黄片无遮挡物在线观看| 国产精品久久久久久精品古装| 2022亚洲国产成人精品| 亚洲五月色婷婷综合| 老鸭窝网址在线观看| 人妻少妇偷人精品九色| 电影成人av| 欧美日韩精品网址| 丝袜人妻中文字幕| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 蜜桃国产av成人99| 国产极品粉嫩免费观看在线| 亚洲国产精品国产精品| 一本色道久久久久久精品综合| 波多野结衣av一区二区av| 亚洲精品乱久久久久久| 日韩电影二区| 国产人伦9x9x在线观看 | 毛片一级片免费看久久久久| 免费高清在线观看日韩| 一本大道久久a久久精品| av.在线天堂| av不卡在线播放| 午夜福利在线观看免费完整高清在| 美国免费a级毛片| 精品国产一区二区三区久久久樱花| 9191精品国产免费久久| 免费观看av网站的网址| 在线 av 中文字幕| 日韩av在线免费看完整版不卡| 午夜激情av网站| 一级,二级,三级黄色视频| 精品少妇一区二区三区视频日本电影 | 久久国产精品大桥未久av| 亚洲欧美一区二区三区久久| 国产免费福利视频在线观看| 曰老女人黄片| 久久久久久久久久久久大奶| 国产精品一二三区在线看| 国产一区二区三区av在线| 亚洲国产欧美日韩在线播放| 亚洲国产av影院在线观看| av在线app专区| 成人国产av品久久久| 久久精品国产亚洲av涩爱| 我的亚洲天堂| videosex国产| 免费在线观看视频国产中文字幕亚洲 | 中文字幕av电影在线播放| 精品酒店卫生间| 久久韩国三级中文字幕| 欧美日韩一区二区视频在线观看视频在线| 在现免费观看毛片| 亚洲av国产av综合av卡| 亚洲av综合色区一区| 欧美bdsm另类| 久久久国产欧美日韩av| 亚洲欧洲日产国产| 高清视频免费观看一区二区| 一二三四中文在线观看免费高清| 熟女电影av网| 9热在线视频观看99| a级毛片在线看网站| 亚洲三区欧美一区| 精品国产乱码久久久久久小说| 久久精品国产亚洲av天美| 日韩大片免费观看网站| 人妻 亚洲 视频| 久久热在线av| 伦理电影大哥的女人| 亚洲国产最新在线播放| 边亲边吃奶的免费视频| 久久久精品区二区三区| 亚洲美女视频黄频| 国产精品久久久av美女十八| 色哟哟·www| 1024香蕉在线观看| 国产成人免费无遮挡视频| 成人亚洲欧美一区二区av| 日本欧美国产在线视频| 亚洲av国产av综合av卡| 中文字幕亚洲精品专区| 一本色道久久久久久精品综合| 精品视频人人做人人爽| 人人澡人人妻人| 女性生殖器流出的白浆| 日本欧美国产在线视频| 日产精品乱码卡一卡2卡三| 老汉色∧v一级毛片| 亚洲国产av影院在线观看| 一级爰片在线观看| av天堂久久9| 国产一区二区三区综合在线观看| 咕卡用的链子| 天天操日日干夜夜撸| 免费在线观看视频国产中文字幕亚洲 | 精品一区在线观看国产| 亚洲精品美女久久av网站| 久久午夜福利片| 高清黄色对白视频在线免费看| 在线免费观看不下载黄p国产| 亚洲色图 男人天堂 中文字幕| 九九爱精品视频在线观看| 亚洲成国产人片在线观看| 久久这里有精品视频免费| 青春草亚洲视频在线观看| 精品人妻熟女毛片av久久网站| 久久午夜福利片| 男女下面插进去视频免费观看| 欧美成人午夜精品| 成人手机av| 久久精品国产亚洲av天美| 777米奇影视久久| 午夜福利乱码中文字幕| 91在线精品国自产拍蜜月| 国产av精品麻豆| 人妻少妇偷人精品九色| av免费观看日本| 老女人水多毛片| 黑人猛操日本美女一级片| 亚洲精品美女久久久久99蜜臀 | 啦啦啦中文免费视频观看日本| 国产精品一区二区在线观看99| 老司机影院成人| 在线观看免费日韩欧美大片| 久久精品国产a三级三级三级| 亚洲综合色网址| 18禁国产床啪视频网站| 亚洲国产欧美网| 国产女主播在线喷水免费视频网站| 丰满饥渴人妻一区二区三| 日韩不卡一区二区三区视频在线| 男女无遮挡免费网站观看| 黑人欧美特级aaaaaa片| 国产精品av久久久久免费| 狂野欧美激情性bbbbbb| 视频区图区小说| 男女无遮挡免费网站观看| 色视频在线一区二区三区| 美女脱内裤让男人舔精品视频| 老司机影院成人| 两性夫妻黄色片| 精品亚洲成a人片在线观看| 多毛熟女@视频| 久久久久国产网址| 亚洲四区av| 只有这里有精品99| 精品亚洲成a人片在线观看| 欧美国产精品一级二级三级| 国产一区二区激情短视频 | 91aial.com中文字幕在线观看| 亚洲国产精品国产精品| 伦理电影免费视频| 亚洲av电影在线进入| 欧美精品亚洲一区二区| 久久精品久久久久久噜噜老黄| 欧美97在线视频| 亚洲情色 制服丝袜| 久久国产精品大桥未久av| 国产精品免费视频内射| 久久久精品国产亚洲av高清涩受| 午夜免费观看性视频| av国产精品久久久久影院| 国产成人精品久久久久久| 啦啦啦啦在线视频资源| 日韩三级伦理在线观看| 久久av网站| 99九九在线精品视频| 一级片免费观看大全| 18禁观看日本| 国产精品欧美亚洲77777| 美女高潮到喷水免费观看| 久久久久精品性色| 国产极品粉嫩免费观看在线| 一区二区日韩欧美中文字幕| av卡一久久| 日韩av在线免费看完整版不卡| 免费黄网站久久成人精品| 国产亚洲午夜精品一区二区久久| 黄色配什么色好看| 国产 精品1| 亚洲欧洲国产日韩| 国产午夜精品一二区理论片| 精品国产国语对白av| 这个男人来自地球电影免费观看 | 国产熟女欧美一区二区| 丰满饥渴人妻一区二区三| 一本—道久久a久久精品蜜桃钙片| 免费看不卡的av| 一级毛片我不卡| 大香蕉久久成人网| 欧美精品高潮呻吟av久久| 观看美女的网站| 最近最新中文字幕大全免费视频 | 国产精品国产三级国产专区5o| 日日撸夜夜添| 欧美+日韩+精品| 老司机亚洲免费影院| 欧美日韩一级在线毛片| 精品视频人人做人人爽| 91精品伊人久久大香线蕉| 男女边摸边吃奶| 十八禁高潮呻吟视频| 久久热在线av| 男女免费视频国产| 999久久久国产精品视频| 一级毛片电影观看| 色网站视频免费| 激情视频va一区二区三区| 久久久久国产一级毛片高清牌| 亚洲精品美女久久av网站| 国产 精品1| 菩萨蛮人人尽说江南好唐韦庄| 熟女av电影| 亚洲四区av| 久久精品国产自在天天线| 两个人免费观看高清视频| 国产免费一区二区三区四区乱码| 亚洲精品,欧美精品| 欧美精品高潮呻吟av久久| 亚洲色图综合在线观看| 欧美精品人与动牲交sv欧美| 久久精品国产a三级三级三级| 久久影院123| 欧美xxⅹ黑人| 黄色 视频免费看| √禁漫天堂资源中文www| 成人亚洲欧美一区二区av| 热re99久久精品国产66热6| 又粗又硬又长又爽又黄的视频| 岛国毛片在线播放| 亚洲精品第二区| 亚洲成国产人片在线观看| 国产欧美亚洲国产| 九色亚洲精品在线播放| 99久久中文字幕三级久久日本| 高清视频免费观看一区二区| av女优亚洲男人天堂| 精品一区二区免费观看| 91精品三级在线观看| 水蜜桃什么品种好| 少妇人妻久久综合中文| 亚洲精品美女久久久久99蜜臀 | 高清欧美精品videossex| av福利片在线| 1024香蕉在线观看| 色婷婷久久久亚洲欧美| 成人国产麻豆网| 欧美黄色片欧美黄色片| 亚洲国产精品成人久久小说| 极品人妻少妇av视频| 久久精品国产亚洲av涩爱| 日韩 亚洲 欧美在线| 亚洲一级一片aⅴ在线观看| www.av在线官网国产| 欧美成人午夜精品| av网站在线播放免费| 午夜福利在线免费观看网站| 超碰97精品在线观看| 97在线视频观看|