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

    中國老年人口死亡率的建模分析

    2017-11-01 09:47:56黃佩佩
    關(guān)鍵詞:余弦死亡率聚類

    黃佩佩,鄭 靜

    (杭州電子科技大學(xué)理學(xué)院,浙江 杭州 310018)

    中國老年人口死亡率的建模分析

    黃佩佩,鄭 靜

    (杭州電子科技大學(xué)理學(xué)院,浙江 杭州 310018)

    基于1963—2014年60歲到110歲人口死亡率,采用k-means聚類方法,以余弦相似度為距離函數(shù),將其聚為3類,針對每類采用時間序列分析的方法進(jìn)行建模.每類中分別選取62歲,86歲和94歲的死亡率數(shù)據(jù),利用1963—2006年數(shù)據(jù)進(jìn)行建模,并用2007—2014年數(shù)據(jù)進(jìn)行驗證,其樣本內(nèi)擬合精度和檢驗樣本預(yù)測精度都比較好,從而驗證了所建模型的有效性和可行性.最后用該模型預(yù)測了2015—2020年62歲,86歲和94歲的人口死亡率.

    死亡率;K-means聚類;ARIMA模型

    0 引 言

    人口死亡率一直是社會關(guān)注的一個焦點(diǎn)問題.隨著生活水平的提高和醫(yī)療衛(wèi)生條件的改善,人類死亡率呈下降趨勢,壽命不斷延長.1981—2010年,根據(jù)《中國統(tǒng)計年鑒》數(shù)據(jù)顯示,我國人均壽命從66.77歲上升至74.83歲,如果按同樣的速度增長,到2100年,我國的人均壽命將超過95歲.老年人口死亡率的準(zhǔn)確預(yù)測為我國社會養(yǎng)老及醫(yī)療衛(wèi)生預(yù)算提供了一定的理論依據(jù).

    目前,人類死亡率預(yù)測模型主要有2類,分別為Lee-Carter模型及其改進(jìn)和時間序列模型.Lee-Carter模型給出了不同年齡的對數(shù)死亡率關(guān)于時間的函數(shù)關(guān)系,其得到了廣泛的應(yīng)用,是當(dāng)今世界上最流行的死亡率預(yù)測模型.但其也有一定缺陷,用一個模型來模擬所有年齡人口死亡率,只能有局部代表意義,一些年齡模型的誤差較大,其人口死亡率預(yù)測不準(zhǔn)確.文獻(xiàn)[1]以中國人口生命表為基礎(chǔ),對2010年的人口按齡死亡模式進(jìn)行修正,研究發(fā)現(xiàn)2010年中國人口死亡率下降至5.58‰.文獻(xiàn)[2]介紹了預(yù)測死亡率常用的Lee-Carter模型的應(yīng)用,估計其方法及預(yù)測原理,并選用該模型對我國未來人口死亡率進(jìn)行了預(yù)測,結(jié)果表明未來人口死亡率將持續(xù)降低.文獻(xiàn)[3]對Lee-Carter模型進(jìn)行了完整的理論研究,并給出了完整的Lee-Carter模型理論分布和區(qū)間預(yù)測表達(dá)式,相比于傳統(tǒng)Lee-Carter預(yù)測方法,得到的預(yù)測區(qū)間較寬,長壽風(fēng)險評估更為準(zhǔn)確.基于以上討論并結(jié)合以往國內(nèi)外學(xué)者對人口死亡率的研究,本文采用時間序列分析的方法對60~110歲老年人口死亡率進(jìn)行了研究.首先對原序列采用k-means聚類的方法,利用余弦相似度對不同年齡死亡率進(jìn)行分類,進(jìn)而對每一類分別進(jìn)行分析建立模型,彌補(bǔ)了直接建模帶來的不足.相對以往研究,此方法建立的模型誤差較小,預(yù)測結(jié)果更為精確.

    1 原時間序列的聚類分析

    1.1 k-means聚類的基本思想

    聚類分析是直接比較各事物之間的性質(zhì),將性質(zhì)相近的歸為一類,性質(zhì)差別較大的歸為不同的類,以使簇內(nèi)有較高相似度,而簇間相似度較低[4].本文聚類采用k-means聚類,其基本思想是:從n個數(shù)據(jù)中任意選取c個對象用來作初始聚類中心,每個對象初始的代表一個簇的平均值或者中心,剩余n-c個對象,根據(jù)其到各個簇的距離,將它們歸類到距離最小的簇中心,然后重新計算每個簇的平均值,重復(fù)此過程,直至聚類準(zhǔn)則函數(shù)收斂;準(zhǔn)則函數(shù)通常為平方誤差準(zhǔn)則:

    (1)

    1.2 余弦相似度

    本文k-means聚類中距離采用余弦相似度度量.余弦相似度是用向量空間中2個向量的夾角余弦作為衡量2個向量之間的貼合程度,其相似度度量為:

    (2)

    其中,x1k,x2k為2個向量第k維的值,n為向量的個數(shù),且cosθ的取值范圍為[0,1].余弦值越大,對象間的相似度越大;余弦值越小,對象間的相似度越?。?/p>

    2 自回歸移動平均模型

    本文選用自回歸移動平均模型即ARMA模型[5],它是由自回歸模型AR和移動平均模型MA組成,其基本公式為:

    yt=φ1yt-1+φ2yt-2+…+φpyt-p+et-θ1et-1-…-θqet-q

    (3)

    其中,p為自回歸模型的階數(shù),φi為自回歸系數(shù),yt-i為時間序列在t-i期的觀察值,et為不能用模型說明的隨機(jī)擾動,q為移動平均模型的階數(shù),θi為移動平均系數(shù),et-i為第t-i個時期進(jìn)入系統(tǒng)的隨機(jī)擾動.建立ARMA模型要求時間序列具有平穩(wěn)性.原序列是不平穩(wěn)的,通過差分使序列達(dá)到平穩(wěn);經(jīng)過差分后建立的ARMA(p,q)模型稱為求和自回歸移動平均模型即ARIMA(p,d,q).

    3 人口死亡率建模分析

    3.1 數(shù)據(jù)分析與處理

    本文研究的老年人口死亡率數(shù)據(jù)從http://www.mortality.org/下載,數(shù)據(jù)為1963—2014年60~110歲人口死亡率,文中所有算法實現(xiàn)都是由R語言來實現(xiàn).?dāng)?shù)據(jù)較多,首先對數(shù)據(jù)采用k-means聚類方法.由于考慮到人口死亡率的走勢圖,將1963—2014年死亡率走勢相近的年齡聚為一類.這樣,每類建立一個模型,并用該模型預(yù)測該類中其他年齡死亡率.基于以上特點(diǎn),聚類中的距離采用余弦相似度,將數(shù)據(jù)聚為3類,分別為60~80歲、81~90歲、91~110歲.

    分別對3類數(shù)據(jù)采用時間序列方法進(jìn)行分析.分別選取62歲、86歲和94歲數(shù)據(jù)作為每一類的代表序列進(jìn)行建模,選取1963—2006年數(shù)據(jù)作為建模樣本,2007—2014年數(shù)據(jù)作為測試樣本來評估模型的優(yōu)劣,并由模型預(yù)測2015—2020年老年人口死亡率.由R軟件的單位根檢驗函數(shù)adf.test知,62歲、86歲和94歲死亡率數(shù)據(jù)的置信水平P值依次為0.838 8,0.380 7和0.539 5,從而數(shù)據(jù)為非平穩(wěn)序列.接著對序列進(jìn)行一階差分得到的新序列Yt的P值依次為0.021 2,0.202 5和0.156 0,從而得到62歲的一階差分序列為平穩(wěn)序列,但86歲和94歲的一階差分序列為非平穩(wěn)序列.然后再對86歲和94歲數(shù)據(jù)進(jìn)行二階差分,得到序列的P值均為0.010 0,則其為平穩(wěn)序列.

    3.2 模型建立與檢驗

    3.2.1 模型診斷

    模型診斷主要包括顯著性檢驗和殘差分析等[6].參數(shù)顯著性檢驗就是檢驗參數(shù)是否顯著異于零.如果參數(shù)不顯著,說明影響不明顯,將其從模型中刪除;如果參數(shù)顯著異于零,則將其保留在最終模型中.對于ARMA(p,q)模型,標(biāo)準(zhǔn)化殘差為:

    (4)

    H0:ARMA(p,q)是充分的.

    H1:ARMA(p,q)是不充分的.

    3.2.2 模型建立與檢驗

    經(jīng)過差分并由ADF檢驗得到62歲為一階平穩(wěn)序列,86歲和94歲為二階平穩(wěn)序列,即62歲、86歲和94歲序列的d分別為1,2,2,且新序列的ACF和PACF如圖1所示.

    圖1 ACF圖

    圖2 PACF圖

    由圖1、圖2可以看出,新序列的自相關(guān)函數(shù)和偏自相關(guān)函數(shù)都是拖尾的[7],因此可建立ARMA模型,計算所有模型的AIC值,記下AIC最小時所對應(yīng)的p,q.對62歲序列,當(dāng)p=1,q=2時,AIC值最小,故此模型為ARIMA(1,1,2);而86歲和94歲序列都是當(dāng)p=1,q=1時,AIC值最小,故這2類模型均為ARIMA(1,2,1).對模型的參數(shù)進(jìn)行最小二乘估計,得到估計參數(shù)結(jié)果如表1所示.

    表1 3個模型的參數(shù)估計

    在表1中,ARIMA(1,1,2)62,ARIMA(1,2,1)86和ARIMA(1,2,1)94分別為62歲、86歲和94歲老年人口死亡率模型的參數(shù)估計,從數(shù)值中可以看出,對于任一模型,參數(shù)都顯著異于零,則由參數(shù)顯著性檢驗知參數(shù)對模型影響顯著.

    然后用這3類模型分別預(yù)測2007—2014年62歲,86歲和94歲老年人口死亡率,預(yù)測結(jié)果如圖3所示.

    圖3 老年人口死亡率預(yù)測值、真實值及擬合值

    從圖3中可以看出,3個模型對于2007—2014年數(shù)據(jù)擬合效果較好,預(yù)測值與真實值在誤差允許的范圍內(nèi)幾乎接近一致.進(jìn)一步通過sarima函數(shù)檢驗3類模型的殘差自相關(guān)性,3類模型的白噪聲方差估計值分別為0.001 8,0.078 2,0.460 4.Ljung-Box統(tǒng)計量結(jié)果表明,3類模型均在自由度為6的情況下,殘差已經(jīng)不存在自相關(guān)性,且P值分別為0.216 8,0.140 6和0.151 0,所有值遠(yuǎn)遠(yuǎn)大于5%,則殘差有很好的獨(dú)立性.由模型診斷的理論可知,這3類模型都是正確設(shè)定的.其中94歲模型的死亡率出現(xiàn)突然上升的情況,并且上升的時間正好處于2003年附近,這與此時間段內(nèi)發(fā)生的非典疫情相吻合,而該年齡的人口基數(shù)小導(dǎo)致上升的最為明顯.

    3.2.3 類模型的誤差分析

    對上述所建3類模型進(jìn)行誤差分析,分別用R軟件計算出模型誤差與檢驗誤差,來檢驗?zāi)P偷膬?yōu)劣.3類模型的模型誤差與檢驗誤差如表2所示.

    表2 3個模型的模型誤差與檢驗誤差

    從表2可以看出,3類模型中,模型誤差和檢驗誤差都比較小,則表明樣本內(nèi)擬合精度和檢驗樣本預(yù)測精度都比較高,進(jìn)一步說明所建模型是正確的.

    3.3 老年人口死亡率預(yù)測值

    由以上分析,3類模型樣本內(nèi)擬合精度和檢驗樣本預(yù)測精度都比較好,實驗結(jié)果表明該模型構(gòu)造是充分的,所以其可用于預(yù)測.分別用此模型預(yù)測2015—2020年老年人口死亡率,如表3所示.

    表3 2015—2020年老年人口死亡率預(yù)測值 %

    由于聚類所選距離設(shè)定為余弦相似度,因此,其他年齡人口死亡率可代入相關(guān)參數(shù)來預(yù)測.由于模型樣本內(nèi)擬合精度和檢驗樣本預(yù)測精度都比較高,所以此預(yù)測值具有一定的參考意義.

    4 結(jié)束語

    本文提出了用ARIMA模型預(yù)測老年人口死亡率.通過研究表明,所建立的模型能較為準(zhǔn)確地描述老年人口死亡率所呈現(xiàn)的特征,為保險公司企業(yè)退休金及商業(yè)保險提供了理論依據(jù),對其制定決策有一定的參考價值;其次,預(yù)測值對國家社會養(yǎng)老保險及醫(yī)療衛(wèi)生預(yù)算也有重大意義.

    [1] 張文娟,魏蒙.中國人口的死亡水平及預(yù)期壽命評估——基于第六次人口普查數(shù)據(jù)的分析[J].人口學(xué)刊,2016,38(3):18-28.

    [2] 蘇華,趙文.我國死亡率預(yù)測與長壽風(fēng)險[J].經(jīng)營管理者,2016,31(6):3-4.

    [3] 王志剛,王曉軍,張學(xué)斌.Lee—Carter模型的理論分布和區(qū)間預(yù)測[J].?dāng)?shù)理統(tǒng)計與管理,2016,35(3):484-493.

    [4] TAN P N,STEINBACH M,KUMAR V.?dāng)?shù)據(jù)挖掘?qū)д?完整版)[M].北京:人民郵電出版社,2011:310-314.

    [5] Cryer J D,Chan K S.時間序列分析及應(yīng)用[M].潘紅宇,譯.北京:機(jī)械工業(yè)出版社,2011:63-133.

    [6] 趙華.時間序列數(shù)據(jù)分析:R軟件應(yīng)用[M].北京:清華大學(xué)出版社,2016:58-88.

    [7] 于寧莉,易東云,涂先勤.時間序列中自相關(guān)與偏相關(guān)函數(shù)分析[J].?dāng)?shù)學(xué)理論與應(yīng)用,2007,27(1):54-57.

    ModelingandAnalysisoftheMortalityoftheElderlyPopulation

    HUANG Peipei, ZHENG Jing

    (SchoolofScience,HangzhouDianziUniversity,HangzhouZhejiang310018,China)

    From 1963 to 2014 at the age of 60 to 110 year old population mortality based on the K-means clustering method, the cosine similarity as the distance function, which will be grouped into 3 categories, each kind of modeling method based on time series analysis. In each class was selected at the age of 62, 86 and 94 year old mortality data, the data from 1963 to 2006 were used to model and test the model with the data from 2007 to 2014. The fitting precision and prediction accuracy of test samples are good, which verifies the effectiveness and feasibility of the model. This model is also used to predict the class of 2015 to 2020 the elderly population mortality.

    mortality; K-means clustering; ARIMA model

    O213

    A

    1001-9146(2017)05-0092-05

    2016-12-14

    黃佩佩(1992-),女,山東濟(jì)寧人,碩士研究生,統(tǒng)計學(xué).通信作者:鄭靜副教授,E-mail:zhengjing@hdu.edu.cn.

    10.13954/j.cnki.hdu.2017.05.017

    猜你喜歡
    余弦死亡率聚類
    走路可以降低死亡率
    中老年保健(2021年4期)2021-08-22 07:07:02
    春季養(yǎng)雞這樣降低死亡率
    新冠肺炎的死亡率為何難確定?
    急性爛鰓、套腸、敗血癥…一旦治療不及時,死亡率或高達(dá)90%,叉尾鮰真的值得養(yǎng)嗎?
    基于DBSACN聚類算法的XML文檔聚類
    電子測試(2017年15期)2017-12-18 07:19:27
    兩個含余弦函數(shù)的三角母不等式及其推論
    分?jǐn)?shù)階余弦變換的卷積定理
    圖像壓縮感知在分?jǐn)?shù)階Fourier域、分?jǐn)?shù)階余弦域的性能比較
    基于改進(jìn)的遺傳算法的模糊聚類算法
    一種層次初始的聚類個數(shù)自適應(yīng)的聚類方法研究
    高清在线视频一区二区三区| 免费观看性生交大片5| 午夜日本视频在线| 精品一品国产午夜福利视频| 久久久久网色| 免费av中文字幕在线| 国产一区二区在线观看av| 男的添女的下面高潮视频| 免费大片黄手机在线观看| 亚洲成人一二三区av| 午夜老司机福利剧场| 黑人猛操日本美女一级片| 久久久久久伊人网av| 国产极品粉嫩免费观看在线 | 久久精品久久久久久噜噜老黄| 少妇人妻 视频| 最新中文字幕久久久久| 国产一区二区在线观看av| 最近最新中文字幕免费大全7| 国产成人精品婷婷| a级毛片免费高清观看在线播放| 在线观看国产h片| 精品少妇久久久久久888优播| 一级毛片我不卡| 亚洲欧美日韩另类电影网站| 日韩一区二区视频免费看| 亚洲av二区三区四区| 边亲边吃奶的免费视频| 成人亚洲欧美一区二区av| 久久鲁丝午夜福利片| 亚洲精品中文字幕在线视频| 精品少妇久久久久久888优播| 在线观看三级黄色| 欧美日韩在线观看h| 夜夜看夜夜爽夜夜摸| 岛国毛片在线播放| 国产淫语在线视频| 韩国高清视频一区二区三区| 国产69精品久久久久777片| 纯流量卡能插随身wifi吗| 老女人水多毛片| h视频一区二区三区| 国产又色又爽无遮挡免| 一个人看视频在线观看www免费| 五月开心婷婷网| 午夜福利在线观看免费完整高清在| 91久久精品电影网| 日本av免费视频播放| 日本vs欧美在线观看视频| 久久午夜综合久久蜜桃| 免费看av在线观看网站| 18+在线观看网站| 夫妻性生交免费视频一级片| 国产成人免费无遮挡视频| 又粗又硬又长又爽又黄的视频| 日本vs欧美在线观看视频| 亚州av有码| 精品久久蜜臀av无| 多毛熟女@视频| 91在线精品国自产拍蜜月| 美女cb高潮喷水在线观看| 免费黄色在线免费观看| 亚洲国产成人一精品久久久| 精品人妻一区二区三区麻豆| 精品99又大又爽又粗少妇毛片| 999精品在线视频| 丁香六月天网| 高清视频免费观看一区二区| 精品国产一区二区三区久久久樱花| 国产精品偷伦视频观看了| 亚洲精品乱久久久久久| 下体分泌物呈黄色| 色5月婷婷丁香| 黄片无遮挡物在线观看| 黄色毛片三级朝国网站| tube8黄色片| 你懂的网址亚洲精品在线观看| 中国三级夫妇交换| 国产黄色视频一区二区在线观看| 午夜福利,免费看| 久久99蜜桃精品久久| 人人妻人人添人人爽欧美一区卜| 亚洲国产成人一精品久久久| 夜夜看夜夜爽夜夜摸| 自拍欧美九色日韩亚洲蝌蚪91| 老司机亚洲免费影院| 日韩视频在线欧美| 天天躁夜夜躁狠狠久久av| av一本久久久久| 在线免费观看不下载黄p国产| 午夜免费男女啪啪视频观看| 少妇被粗大猛烈的视频| 一本大道久久a久久精品| 欧美一级a爱片免费观看看| 亚洲第一区二区三区不卡| 你懂的网址亚洲精品在线观看| 曰老女人黄片| 婷婷成人精品国产| 91成人精品电影| 精品少妇黑人巨大在线播放| 汤姆久久久久久久影院中文字幕| 亚洲一级一片aⅴ在线观看| 国产高清国产精品国产三级| 久久ye,这里只有精品| 亚洲高清免费不卡视频| 久久精品国产亚洲网站| 热re99久久精品国产66热6| 亚洲国产色片| 欧美日韩亚洲高清精品| 少妇的逼水好多| 插逼视频在线观看| 国产精品无大码| 国产国拍精品亚洲av在线观看| 久久97久久精品| 亚洲人与动物交配视频| 老女人水多毛片| 女人久久www免费人成看片| 肉色欧美久久久久久久蜜桃| 午夜av观看不卡| 最黄视频免费看| 乱人伦中国视频| 中文字幕亚洲精品专区| 国产欧美日韩一区二区三区在线 | 免费看光身美女| 日韩av不卡免费在线播放| 国产免费福利视频在线观看| av在线观看视频网站免费| 久久午夜综合久久蜜桃| 老熟女久久久| 永久免费av网站大全| av一本久久久久| 亚洲成人一二三区av| videosex国产| 日韩三级伦理在线观看| 欧美亚洲日本最大视频资源| 亚洲人成77777在线视频| 日韩大片免费观看网站| 亚洲精品成人av观看孕妇| 99re6热这里在线精品视频| 搡女人真爽免费视频火全软件| 日韩免费高清中文字幕av| 国产免费又黄又爽又色| 大香蕉久久成人网| 午夜福利影视在线免费观看| 蜜桃在线观看..| 五月伊人婷婷丁香| 亚洲国产精品成人久久小说| 精品酒店卫生间| 日本wwww免费看| av在线播放精品| 嫩草影院入口| 人妻少妇偷人精品九色| 亚洲欧美色中文字幕在线| 九九在线视频观看精品| 欧美日韩视频精品一区| av一本久久久久| 高清欧美精品videossex| 极品少妇高潮喷水抽搐| 99精国产麻豆久久婷婷| 亚洲国产精品专区欧美| 亚洲,一卡二卡三卡| 蜜桃国产av成人99| 18+在线观看网站| 精品熟女少妇av免费看| 国产熟女午夜一区二区三区 | 亚洲人与动物交配视频| 成年av动漫网址| 欧美一级a爱片免费观看看| 亚洲精品久久成人aⅴ小说 | 亚洲欧美一区二区三区黑人 | 天天躁夜夜躁狠狠久久av| 日韩成人av中文字幕在线观看| av在线app专区| 成人毛片a级毛片在线播放| 亚洲av电影在线观看一区二区三区| 国产爽快片一区二区三区| 午夜福利在线观看免费完整高清在| 久久ye,这里只有精品| 内地一区二区视频在线| freevideosex欧美| 国产精品久久久久成人av| 国产又色又爽无遮挡免| 一区二区av电影网| 亚洲精品第二区| 国产在线免费精品| 亚洲精品国产色婷婷电影| 亚洲国产成人一精品久久久| 亚洲精品自拍成人| 嫩草影院入口| 欧美人与善性xxx| 中国三级夫妇交换| 亚洲精品国产av蜜桃| 成人亚洲欧美一区二区av| 亚洲怡红院男人天堂| 国产欧美另类精品又又久久亚洲欧美| 亚洲内射少妇av| 午夜免费鲁丝| 97精品久久久久久久久久精品| 国产精品久久久久久久久免| 精品人妻在线不人妻| 天天躁夜夜躁狠狠久久av| 一区二区av电影网| 亚洲精品国产色婷婷电影| 人妻人人澡人人爽人人| 国产69精品久久久久777片| 午夜视频国产福利| 国产成人精品久久久久久| 五月玫瑰六月丁香| 精品一区二区三区视频在线| 免费高清在线观看日韩| 国产综合精华液| 国产一区有黄有色的免费视频| 免费高清在线观看视频在线观看| 性色avwww在线观看| 中国美白少妇内射xxxbb| 色吧在线观看| 另类精品久久| 亚洲精品自拍成人| 女性被躁到高潮视频| 不卡视频在线观看欧美| 欧美少妇被猛烈插入视频| 老熟女久久久| 美女主播在线视频| 国产精品国产三级国产专区5o| 欧美最新免费一区二区三区| 一级毛片 在线播放| 精品一区在线观看国产| 99久久综合免费| 建设人人有责人人尽责人人享有的| 亚洲av在线观看美女高潮| 精品一品国产午夜福利视频| 亚洲欧美精品自产自拍| 亚州av有码| 高清欧美精品videossex| 又黄又爽又刺激的免费视频.| 成年人免费黄色播放视频| 日本与韩国留学比较| 国语对白做爰xxxⅹ性视频网站| 日日摸夜夜添夜夜添av毛片| 国产一区亚洲一区在线观看| 国产爽快片一区二区三区| 亚洲精品亚洲一区二区| 在线 av 中文字幕| 成人亚洲精品一区在线观看| 一级a做视频免费观看| 99热这里只有精品一区| 午夜av观看不卡| a 毛片基地| 亚洲人成77777在线视频| 成人18禁高潮啪啪吃奶动态图 | 亚洲av国产av综合av卡| 亚洲精品国产av成人精品| 精品人妻熟女毛片av久久网站| 免费久久久久久久精品成人欧美视频 | 欧美日韩精品成人综合77777| a 毛片基地| 精品国产一区二区三区久久久樱花| 久久人人爽av亚洲精品天堂| 成人综合一区亚洲| 桃花免费在线播放| 亚洲不卡免费看| 亚洲精品,欧美精品| 国产精品久久久久成人av| 久久久久国产网址| 国产日韩一区二区三区精品不卡 | 少妇丰满av| 最近2019中文字幕mv第一页| 寂寞人妻少妇视频99o| 日日摸夜夜添夜夜爱| 国产不卡av网站在线观看| 97超视频在线观看视频| 亚洲国产精品国产精品| 久久久久精品久久久久真实原创| 亚洲五月色婷婷综合| 天堂中文最新版在线下载| 一区二区三区四区激情视频| 91精品三级在线观看| 欧美国产精品一级二级三级| 国产一区有黄有色的免费视频| 女性生殖器流出的白浆| 亚洲精品成人av观看孕妇| 乱码一卡2卡4卡精品| 婷婷色综合www| 女的被弄到高潮叫床怎么办| 午夜视频国产福利| 蜜臀久久99精品久久宅男| 国产国拍精品亚洲av在线观看| 黄色视频在线播放观看不卡| 欧美97在线视频| 亚洲精品,欧美精品| 春色校园在线视频观看| 国产又色又爽无遮挡免| 丝袜在线中文字幕| 国产一区二区三区综合在线观看 | 国产爽快片一区二区三区| 日日爽夜夜爽网站| 久久久久久久精品精品| 午夜福利视频在线观看免费| 日韩av不卡免费在线播放| 国产熟女午夜一区二区三区 | 欧美97在线视频| 久久人妻熟女aⅴ| 久久久久久伊人网av| 久久精品久久久久久久性| 成年美女黄网站色视频大全免费 | 国产成人午夜福利电影在线观看| 91久久精品电影网| 亚洲婷婷狠狠爱综合网| 亚洲国产毛片av蜜桃av| 久久热精品热| 18禁在线播放成人免费| 极品少妇高潮喷水抽搐| 日本av免费视频播放| 精品国产国语对白av| 国产高清三级在线| 边亲边吃奶的免费视频| 男女高潮啪啪啪动态图| 狠狠婷婷综合久久久久久88av| 日本免费在线观看一区| 亚洲av男天堂| 色婷婷久久久亚洲欧美| 国产淫语在线视频| 最黄视频免费看| 色婷婷久久久亚洲欧美| 久久久国产欧美日韩av| 2022亚洲国产成人精品| 在线观看免费高清a一片| 欧美97在线视频| 成人黄色视频免费在线看| 亚洲精品国产av成人精品| 久热这里只有精品99| 欧美日韩精品成人综合77777| 91精品伊人久久大香线蕉| 日韩视频在线欧美| 亚洲av成人精品一二三区| 一本一本综合久久| 日韩欧美一区视频在线观看| 午夜av观看不卡| 啦啦啦视频在线资源免费观看| 嘟嘟电影网在线观看| 性高湖久久久久久久久免费观看| 中文字幕最新亚洲高清| 国产午夜精品一二区理论片| 2018国产大陆天天弄谢| 日韩电影二区| 欧美少妇被猛烈插入视频| h视频一区二区三区| 国产国拍精品亚洲av在线观看| 亚洲第一av免费看| 亚洲内射少妇av| 日韩视频在线欧美| 亚洲欧洲日产国产| 久久毛片免费看一区二区三区| 男女边吃奶边做爰视频| 午夜免费男女啪啪视频观看| 亚洲国产精品一区三区| 久久久久久久久久成人| 久久97久久精品| 免费大片黄手机在线观看| 免费观看性生交大片5| 亚洲久久久国产精品| 亚洲欧洲精品一区二区精品久久久 | 天天影视国产精品| 一级片'在线观看视频| 一级爰片在线观看| 国产国拍精品亚洲av在线观看| 国产精品麻豆人妻色哟哟久久| 久久久久久伊人网av| 中文字幕人妻熟人妻熟丝袜美| 日日撸夜夜添| av在线老鸭窝| 成人国语在线视频| tube8黄色片| 各种免费的搞黄视频| 中文字幕人妻熟人妻熟丝袜美| 99热网站在线观看| 亚洲欧美色中文字幕在线| av又黄又爽大尺度在线免费看| 免费黄频网站在线观看国产| 欧美精品国产亚洲| 男人添女人高潮全过程视频| 亚洲内射少妇av| 秋霞伦理黄片| 菩萨蛮人人尽说江南好唐韦庄| 亚洲四区av| 国产午夜精品一二区理论片| 国产 一区精品| 大香蕉97超碰在线| 久久精品国产亚洲av涩爱| 国产国语露脸激情在线看| 老司机亚洲免费影院| 人人妻人人爽人人添夜夜欢视频| 少妇丰满av| 51国产日韩欧美| 女人久久www免费人成看片| 有码 亚洲区| 国产精品蜜桃在线观看| 免费高清在线观看日韩| 亚洲精品久久午夜乱码| 日韩一区二区视频免费看| 最黄视频免费看| 久久午夜福利片| 男人添女人高潮全过程视频| 永久免费av网站大全| 午夜免费鲁丝| 免费黄频网站在线观看国产| 国产极品天堂在线| 王馨瑶露胸无遮挡在线观看| 午夜视频国产福利| 3wmmmm亚洲av在线观看| 亚洲精品第二区| 熟女电影av网| 国产片特级美女逼逼视频| 伊人亚洲综合成人网| 亚洲五月色婷婷综合| 国产精品成人在线| 亚洲美女搞黄在线观看| 18禁裸乳无遮挡动漫免费视频| 天天影视国产精品| 激情五月婷婷亚洲| 一个人看视频在线观看www免费| 免费观看无遮挡的男女| 夜夜看夜夜爽夜夜摸| 男人爽女人下面视频在线观看| 国产精品 国内视频| 免费黄色在线免费观看| 一区在线观看完整版| 高清毛片免费看| 一本—道久久a久久精品蜜桃钙片| 久久久久久久亚洲中文字幕| 男女国产视频网站| 看免费成人av毛片| 亚洲精品国产av成人精品| 我的老师免费观看完整版| 黄色视频在线播放观看不卡| 国产精品无大码| 国产欧美日韩一区二区三区在线 | 亚洲国产精品国产精品| 高清av免费在线| 考比视频在线观看| 18禁在线播放成人免费| 色94色欧美一区二区| 99久久精品国产国产毛片| 人成视频在线观看免费观看| 欧美3d第一页| www.色视频.com| 成人综合一区亚洲| 9色porny在线观看| 看十八女毛片水多多多| 亚洲av免费高清在线观看| 日韩人妻高清精品专区| 男女边摸边吃奶| 国产在视频线精品| 亚洲精品乱码久久久久久按摩| 99热6这里只有精品| 日本午夜av视频| 欧美日韩亚洲高清精品| 在线观看免费日韩欧美大片 | 免费人妻精品一区二区三区视频| 麻豆乱淫一区二区| 一级,二级,三级黄色视频| 九色亚洲精品在线播放| 亚洲熟女精品中文字幕| av在线app专区| 亚洲内射少妇av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日产精品乱码卡一卡2卡三| 亚洲精品中文字幕在线视频| www.色视频.com| 一级毛片aaaaaa免费看小| 青春草亚洲视频在线观看| 久久久久久久久久成人| av电影中文网址| 日韩人妻高清精品专区| 中文字幕制服av| 久久久久视频综合| 中国美白少妇内射xxxbb| 免费观看无遮挡的男女| 亚洲av日韩在线播放| 亚洲精品第二区| 日韩伦理黄色片| 少妇精品久久久久久久| 在线观看国产h片| 欧美激情极品国产一区二区三区 | 天天影视国产精品| 99九九线精品视频在线观看视频| 免费av不卡在线播放| 中文乱码字字幕精品一区二区三区| 大香蕉久久成人网| 涩涩av久久男人的天堂| 搡老乐熟女国产| 久久久久久久久久久免费av| 美女国产高潮福利片在线看| av免费观看日本| av在线老鸭窝| 18禁在线播放成人免费| 美女脱内裤让男人舔精品视频| 国产 精品1| 国产午夜精品久久久久久一区二区三区| 熟妇人妻不卡中文字幕| 高清av免费在线| kizo精华| 嫩草影院入口| 91午夜精品亚洲一区二区三区| 在线亚洲精品国产二区图片欧美 | 另类亚洲欧美激情| 一二三四中文在线观看免费高清| 亚洲欧美精品自产自拍| 国产亚洲一区二区精品| 女人精品久久久久毛片| 国产又色又爽无遮挡免| 国产免费一区二区三区四区乱码| 午夜福利网站1000一区二区三区| 超碰97精品在线观看| 色吧在线观看| 熟女av电影| 99久久综合免费| 日韩人妻高清精品专区| 欧美成人午夜免费资源| 欧美另类一区| 女人久久www免费人成看片| 曰老女人黄片| 国产极品粉嫩免费观看在线 | 亚洲av成人精品一二三区| 男女无遮挡免费网站观看| 日本av手机在线免费观看| 永久网站在线| 香蕉精品网在线| 精品午夜福利在线看| 久久精品国产亚洲av涩爱| 大香蕉97超碰在线| 亚洲天堂av无毛| 在线观看人妻少妇| 黄片播放在线免费| 午夜福利在线观看免费完整高清在| 妹子高潮喷水视频| 爱豆传媒免费全集在线观看| 亚洲在久久综合| 免费日韩欧美在线观看| 日韩中字成人| 一区二区日韩欧美中文字幕 | 日韩免费高清中文字幕av| 2022亚洲国产成人精品| 丁香六月天网| 天堂8中文在线网| 国产有黄有色有爽视频| 妹子高潮喷水视频| 亚洲av福利一区| 久久久久国产网址| 水蜜桃什么品种好| www.色视频.com| a级片在线免费高清观看视频| av国产久精品久网站免费入址| 亚洲精品亚洲一区二区| 精品人妻一区二区三区麻豆| 美女xxoo啪啪120秒动态图| 国产av精品麻豆| 丰满饥渴人妻一区二区三| 黑人高潮一二区| 国产一区二区三区综合在线观看 | 精品一品国产午夜福利视频| av卡一久久| 18在线观看网站| 一个人免费看片子| av国产精品久久久久影院| 久久99一区二区三区| 香蕉精品网在线| 日本猛色少妇xxxxx猛交久久| 成人免费观看视频高清| 韩国av在线不卡| 国产精品国产三级国产av玫瑰| 精品少妇黑人巨大在线播放| 国产片特级美女逼逼视频| 日韩不卡一区二区三区视频在线| 免费看不卡的av| 我要看黄色一级片免费的| 美女主播在线视频| av视频免费观看在线观看| av国产久精品久网站免费入址| 91国产中文字幕| 男人添女人高潮全过程视频| 亚洲成人手机| 亚洲精品久久久久久婷婷小说| 老女人水多毛片| 免费观看性生交大片5| 国产精品 国内视频| 国产日韩欧美亚洲二区| 久久久久久久精品精品| 久久久亚洲精品成人影院| 丁香六月天网| 久久精品久久精品一区二区三区| 日韩av不卡免费在线播放| 精品亚洲成国产av| 午夜免费鲁丝| a级毛片在线看网站| 精品亚洲成国产av| 天美传媒精品一区二区| 国产永久视频网站| 丁香六月天网| 亚洲精品久久午夜乱码| a级毛片在线看网站| 精品亚洲乱码少妇综合久久| 天堂8中文在线网| 母亲3免费完整高清在线观看 | 久久人人爽av亚洲精品天堂| 亚洲精品aⅴ在线观看| 黑人巨大精品欧美一区二区蜜桃 | 天堂俺去俺来也www色官网| xxx大片免费视频| 国产成人精品在线电影| 99热这里只有是精品在线观看| 一级片'在线观看视频| 色网站视频免费| 飞空精品影院首页| 国产高清国产精品国产三级|