• <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)的聚類方法研究
    亚洲最大成人中文| 97超级碰碰碰精品色视频在线观看| 国产蜜桃级精品一区二区三区| 亚洲美女黄片视频| 日韩一本色道免费dvd| 久久精品国产鲁丝片午夜精品| 毛片一级片免费看久久久久| 偷拍熟女少妇极品色| 美女xxoo啪啪120秒动态图| 国产精品女同一区二区软件| 久久久精品94久久精品| 欧美又色又爽又黄视频| 中文资源天堂在线| 日韩欧美国产在线观看| 国产精品伦人一区二区| 成人性生交大片免费视频hd| 国产淫片久久久久久久久| 中国美女看黄片| 听说在线观看完整版免费高清| 久久久久久久午夜电影| 中文字幕熟女人妻在线| 国产成人福利小说| 精品午夜福利在线看| 亚洲av成人av| 日本精品一区二区三区蜜桃| 男女下面进入的视频免费午夜| 国产一区二区三区av在线 | 久久久久九九精品影院| 亚洲国产精品成人综合色| 99国产精品一区二区蜜桃av| 偷拍熟女少妇极品色| 午夜影院日韩av| av天堂中文字幕网| 亚洲专区国产一区二区| 亚洲av成人av| 亚洲av.av天堂| 一区福利在线观看| 深夜a级毛片| 国产高清视频在线播放一区| 高清午夜精品一区二区三区 | 日韩欧美在线乱码| 亚洲欧美清纯卡通| 国产一区二区亚洲精品在线观看| 嫩草影视91久久| 欧美高清成人免费视频www| 精品人妻一区二区三区麻豆 | 亚洲欧美日韩卡通动漫| 99久久久亚洲精品蜜臀av| 亚洲一区高清亚洲精品| 日日撸夜夜添| 久久精品国产亚洲av涩爱 | 激情 狠狠 欧美| 丰满乱子伦码专区| 草草在线视频免费看| 97人妻精品一区二区三区麻豆| 全区人妻精品视频| 麻豆久久精品国产亚洲av| 波多野结衣高清无吗| 国产av在哪里看| 一级a爱片免费观看的视频| 老熟妇仑乱视频hdxx| 伦精品一区二区三区| 特级一级黄色大片| 亚洲av电影不卡..在线观看| 插阴视频在线观看视频| 日韩三级伦理在线观看| 久久精品91蜜桃| 直男gayav资源| 国产真实乱freesex| 亚州av有码| av.在线天堂| 久久久精品大字幕| 男人的好看免费观看在线视频| 国产高清激情床上av| 一个人看视频在线观看www免费| 热99re8久久精品国产| 国产成年人精品一区二区| h日本视频在线播放| 99热这里只有精品一区| 免费人成视频x8x8入口观看| 午夜a级毛片| 白带黄色成豆腐渣| 级片在线观看| 97热精品久久久久久| 看非洲黑人一级黄片| 欧美激情国产日韩精品一区| 欧美丝袜亚洲另类| 亚洲综合色惰| 亚洲精品一区av在线观看| 国产一级毛片七仙女欲春2| 国产久久久一区二区三区| а√天堂www在线а√下载| 亚洲精品日韩在线中文字幕 | 三级毛片av免费| 午夜老司机福利剧场| 日韩亚洲欧美综合| 一本一本综合久久| 久久精品夜色国产| 国产片特级美女逼逼视频| 免费电影在线观看免费观看| 搡老熟女国产l中国老女人| 桃色一区二区三区在线观看| 我要看日韩黄色一级片| 色综合色国产| 国产单亲对白刺激| 少妇的逼好多水| 国产精品1区2区在线观看.| 国产综合懂色| 成人av在线播放网站| 菩萨蛮人人尽说江南好唐韦庄 | 99热这里只有是精品在线观看| 男人舔女人下体高潮全视频| 日韩在线高清观看一区二区三区| 校园人妻丝袜中文字幕| 一个人看的www免费观看视频| 日韩强制内射视频| 久久精品国产亚洲av天美| 久久精品综合一区二区三区| 亚洲精品亚洲一区二区| 在线观看66精品国产| 亚洲欧美成人精品一区二区| 99久国产av精品| 亚洲精品影视一区二区三区av| 久久久成人免费电影| 色5月婷婷丁香| 校园春色视频在线观看| 国产精品久久久久久久久免| 国产美女午夜福利| 一区二区三区四区激情视频 | 99精品在免费线老司机午夜| 欧美色欧美亚洲另类二区| av在线播放精品| 国产真实伦视频高清在线观看| 亚洲电影在线观看av| 看十八女毛片水多多多| 老师上课跳d突然被开到最大视频| 麻豆乱淫一区二区| 一个人免费在线观看电影| 久久精品国产自在天天线| 成人av一区二区三区在线看| 成人综合一区亚洲| 国产精品一区二区免费欧美| 精品久久久久久久久久久久久| 熟女人妻精品中文字幕| 国产视频内射| 亚洲第一电影网av| 十八禁国产超污无遮挡网站| 最近手机中文字幕大全| 男女视频在线观看网站免费| 亚洲婷婷狠狠爱综合网| 成年免费大片在线观看| 国产激情偷乱视频一区二区| 黄色配什么色好看| 国产精品综合久久久久久久免费| 国产极品精品免费视频能看的| 国产免费一级a男人的天堂| 精品人妻偷拍中文字幕| av在线蜜桃| 国产女主播在线喷水免费视频网站 | 日本免费一区二区三区高清不卡| 男女下面进入的视频免费午夜| 色噜噜av男人的天堂激情| 十八禁网站免费在线| 一级毛片久久久久久久久女| 欧美激情在线99| 国产精品久久视频播放| 亚洲最大成人手机在线| 毛片女人毛片| 非洲黑人性xxxx精品又粗又长| 极品教师在线视频| 欧美潮喷喷水| 我的老师免费观看完整版| 日韩中字成人| 五月玫瑰六月丁香| 卡戴珊不雅视频在线播放| 午夜福利视频1000在线观看| 男人狂女人下面高潮的视频| 淫妇啪啪啪对白视频| 人妻制服诱惑在线中文字幕| 婷婷亚洲欧美| 天天躁夜夜躁狠狠久久av| 最好的美女福利视频网| 美女免费视频网站| 有码 亚洲区| 亚洲最大成人手机在线| 99久久成人亚洲精品观看| 国产精品电影一区二区三区| 白带黄色成豆腐渣| 久久精品91蜜桃| 国内精品美女久久久久久| 久久热精品热| 日韩制服骚丝袜av| 深夜精品福利| 女生性感内裤真人,穿戴方法视频| www日本黄色视频网| 九九爱精品视频在线观看| 少妇熟女欧美另类| 亚洲成人久久爱视频| 欧美丝袜亚洲另类| 午夜老司机福利剧场| 日韩欧美三级三区| 99精品在免费线老司机午夜| 久久久色成人| 国产精华一区二区三区| 国产精品伦人一区二区| 啦啦啦韩国在线观看视频| 国产毛片a区久久久久| 免费av不卡在线播放| 插逼视频在线观看| av在线观看视频网站免费| 日韩av在线大香蕉| 国产日本99.免费观看| 成年av动漫网址| 免费观看精品视频网站| 久久久久免费精品人妻一区二区| 国产午夜福利久久久久久| 97人妻精品一区二区三区麻豆| 91久久精品国产一区二区三区| 国产高清有码在线观看视频| 国产精品一区www在线观看| 亚洲国产精品成人久久小说 | 亚洲av成人av| 日韩高清综合在线| 非洲黑人性xxxx精品又粗又长| 久久精品国产鲁丝片午夜精品| 国产精品99久久久久久久久| 国产午夜福利久久久久久| 色吧在线观看| 噜噜噜噜噜久久久久久91| 好男人在线观看高清免费视频| 直男gayav资源| 国产毛片a区久久久久| 蜜臀久久99精品久久宅男| 亚洲熟妇中文字幕五十中出| 午夜福利在线观看吧| 午夜精品国产一区二区电影 | 亚洲人成网站高清观看| h日本视频在线播放| 日本黄大片高清| 久久久久国产网址| 欧美激情国产日韩精品一区| 两个人的视频大全免费| 五月伊人婷婷丁香| 日韩精品青青久久久久久| 国产高潮美女av| 成年女人看的毛片在线观看| h日本视频在线播放| 内地一区二区视频在线| 亚洲av熟女| 在线天堂最新版资源| 成人欧美大片| videossex国产| 噜噜噜噜噜久久久久久91| 91av网一区二区| 欧美区成人在线视频| a级一级毛片免费在线观看| 日本黄色片子视频| 日本爱情动作片www.在线观看 | 久久久精品欧美日韩精品| 两个人的视频大全免费| 日本熟妇午夜| 大又大粗又爽又黄少妇毛片口| 精品一区二区三区视频在线| 成人漫画全彩无遮挡| 国内揄拍国产精品人妻在线| 成人精品一区二区免费| 国产久久久一区二区三区| 午夜爱爱视频在线播放| 卡戴珊不雅视频在线播放| 亚洲熟妇中文字幕五十中出| 日本五十路高清| 人人妻人人澡人人爽人人夜夜 | 国产视频一区二区在线看| 在线播放国产精品三级| 久久精品国产亚洲网站| 国产成人福利小说| 欧美日韩在线观看h| 中文亚洲av片在线观看爽| 乱码一卡2卡4卡精品| 三级毛片av免费| 亚洲三级黄色毛片| 日本成人三级电影网站| 国产不卡一卡二| 欧美一区二区亚洲| 亚洲人与动物交配视频| 精品午夜福利视频在线观看一区| 日日摸夜夜添夜夜添av毛片| 欧美日韩综合久久久久久| 波多野结衣高清无吗| 国产伦精品一区二区三区视频9| a级毛片a级免费在线| 简卡轻食公司| 国产激情偷乱视频一区二区| 久久久久久九九精品二区国产| 97在线视频观看| 国产美女午夜福利| 精品人妻一区二区三区麻豆 | 国产色婷婷99| 欧美不卡视频在线免费观看| 亚洲精品日韩av片在线观看| 亚洲性夜色夜夜综合| 国产v大片淫在线免费观看| 亚洲专区国产一区二区| 禁无遮挡网站| 亚洲av中文字字幕乱码综合| 人人妻人人澡欧美一区二区| 丰满乱子伦码专区| 91精品国产九色| 欧美成人免费av一区二区三区| 精华霜和精华液先用哪个| 亚洲av中文字字幕乱码综合| 午夜爱爱视频在线播放| 日本免费一区二区三区高清不卡| 可以在线观看毛片的网站| 国产伦精品一区二区三区四那| 日韩成人伦理影院| 最好的美女福利视频网| 亚洲av免费在线观看| 好男人在线观看高清免费视频| 成人特级黄色片久久久久久久| 18禁在线无遮挡免费观看视频 | 少妇裸体淫交视频免费看高清| 色噜噜av男人的天堂激情| 亚洲国产精品sss在线观看| 国产精品日韩av在线免费观看| 国产美女午夜福利| 十八禁国产超污无遮挡网站| 我的老师免费观看完整版| 免费观看人在逋| 国产久久久一区二区三区| 久久久久免费精品人妻一区二区| 亚洲中文日韩欧美视频| 成年av动漫网址| 国产精品一及| 高清毛片免费观看视频网站| 午夜激情福利司机影院| 精品人妻偷拍中文字幕| 菩萨蛮人人尽说江南好唐韦庄 | 大型黄色视频在线免费观看| 舔av片在线| 久久草成人影院| 网址你懂的国产日韩在线| 精品一区二区免费观看| 一级毛片久久久久久久久女| 国产精品伦人一区二区| 黄色日韩在线| 草草在线视频免费看| 国产私拍福利视频在线观看| 欧美成人a在线观看| 亚洲精品456在线播放app| 亚洲人成网站在线观看播放| 国内少妇人妻偷人精品xxx网站| 激情 狠狠 欧美| 国产真实乱freesex| 亚洲av中文字字幕乱码综合| 国产精品无大码| 热99re8久久精品国产| 精品午夜福利在线看| 熟女人妻精品中文字幕| 久久久久久久久久成人| 国产真实伦视频高清在线观看| 亚洲aⅴ乱码一区二区在线播放| 淫妇啪啪啪对白视频| 国产精品一及| 国产v大片淫在线免费观看| 国产亚洲精品久久久久久毛片| 日韩制服骚丝袜av| av卡一久久| 成年av动漫网址| 国产探花极品一区二区| 一进一出好大好爽视频| 国产淫片久久久久久久久| 91在线观看av| 成人三级黄色视频| 国产成人影院久久av| 中文字幕久久专区| 中文字幕av成人在线电影| 在线观看美女被高潮喷水网站| 国产精品一二三区在线看| 色综合色国产| 亚洲av免费在线观看| 最近的中文字幕免费完整| 成人高潮视频无遮挡免费网站| 免费电影在线观看免费观看| 别揉我奶头~嗯~啊~动态视频| 亚洲av不卡在线观看| 免费看av在线观看网站| 俄罗斯特黄特色一大片| 亚洲综合色惰| 日韩强制内射视频| 国产又黄又爽又无遮挡在线| 久久国产乱子免费精品| 在线免费十八禁| 色综合色国产| 俄罗斯特黄特色一大片| 欧美潮喷喷水| 丰满人妻一区二区三区视频av| 亚洲欧美日韩卡通动漫| 日本欧美国产在线视频| 最近手机中文字幕大全| 99久国产av精品国产电影| 两个人视频免费观看高清| 欧美色欧美亚洲另类二区| 日韩人妻高清精品专区| 一级黄片播放器| 神马国产精品三级电影在线观看| 99在线视频只有这里精品首页| ponron亚洲| 在线免费观看不下载黄p国产| 欧美+亚洲+日韩+国产| 韩国av在线不卡| 亚洲美女黄片视频| 乱系列少妇在线播放| 亚洲无线观看免费| 亚洲美女视频黄频| 国产精品三级大全| 性色avwww在线观看| av黄色大香蕉| 中文字幕久久专区| 两个人视频免费观看高清| 国产麻豆成人av免费视频| 欧美激情久久久久久爽电影| 国产精品一区二区三区四区免费观看 | 亚洲国产欧美人成| 久久天躁狠狠躁夜夜2o2o| 国产91av在线免费观看| 亚洲第一区二区三区不卡| 国产欧美日韩精品亚洲av| 波多野结衣高清无吗| 免费电影在线观看免费观看| 一区二区三区四区激情视频 | 九九热线精品视视频播放| 色播亚洲综合网| 精品日产1卡2卡| 搡老熟女国产l中国老女人| 国产一区二区在线观看日韩| 国产精品av视频在线免费观看| 人妻夜夜爽99麻豆av| 国产高清视频在线观看网站| 91久久精品国产一区二区成人| 欧美极品一区二区三区四区| 3wmmmm亚洲av在线观看| 国内久久婷婷六月综合欲色啪| av免费在线看不卡| 真人做人爱边吃奶动态| 最近视频中文字幕2019在线8| 国产成人一区二区在线| 黄色视频,在线免费观看| 人人妻人人看人人澡| 美女xxoo啪啪120秒动态图| 成人永久免费在线观看视频| 夜夜爽天天搞| 国产精品一区二区性色av| 精品无人区乱码1区二区| 欧美日韩在线观看h| 欧美高清性xxxxhd video| 日韩欧美三级三区| 中文字幕av在线有码专区| 校园人妻丝袜中文字幕| 成人亚洲欧美一区二区av| 99在线视频只有这里精品首页| 99国产极品粉嫩在线观看| 欧美成人精品欧美一级黄| 菩萨蛮人人尽说江南好唐韦庄 | 男女之事视频高清在线观看| 一个人看的www免费观看视频| 插阴视频在线观看视频| 又粗又爽又猛毛片免费看| 免费看日本二区| 99热网站在线观看| 国产探花在线观看一区二区| 日韩欧美免费精品| 床上黄色一级片| 国产亚洲精品综合一区在线观看| 亚洲国产精品成人综合色| 99视频精品全部免费 在线| 好男人在线观看高清免费视频| 人人妻,人人澡人人爽秒播| 亚洲一区高清亚洲精品| 3wmmmm亚洲av在线观看| 天堂网av新在线| 黄色视频,在线免费观看| 伦理电影大哥的女人| 午夜免费激情av| 国产高清视频在线播放一区| 色视频www国产| 久久久成人免费电影| 一级毛片aaaaaa免费看小| 在线观看美女被高潮喷水网站| 欧美最新免费一区二区三区| 国产精品美女特级片免费视频播放器| 精品人妻视频免费看| 中文在线观看免费www的网站| 精品久久久久久久久av| 精品少妇黑人巨大在线播放 | 99九九线精品视频在线观看视频| .国产精品久久| 人妻少妇偷人精品九色| 日韩欧美精品v在线| 午夜久久久久精精品| 国产精品日韩av在线免费观看| 少妇熟女aⅴ在线视频| 国产亚洲精品久久久com| 亚洲真实伦在线观看| 国产精品永久免费网站| 国产精品人妻久久久影院| 亚洲中文日韩欧美视频| 国产精品久久久久久av不卡| 久久久久久九九精品二区国产| 亚洲欧美精品综合久久99| 亚洲一级一片aⅴ在线观看| 真人做人爱边吃奶动态| 久久韩国三级中文字幕| 国内少妇人妻偷人精品xxx网站| 91精品国产九色| 亚洲精品色激情综合| 色5月婷婷丁香| 国产成人影院久久av| 人人妻,人人澡人人爽秒播| 亚洲图色成人| 色哟哟·www| 国产精品国产高清国产av| 国产片特级美女逼逼视频| 色在线成人网| 亚洲av成人精品一区久久| 一级毛片久久久久久久久女| 久久精品91蜜桃| 啦啦啦啦在线视频资源| 亚洲国产高清在线一区二区三| 亚洲精品成人久久久久久| 婷婷色综合大香蕉| 亚洲欧美清纯卡通| 亚洲成av人片在线播放无| 欧美又色又爽又黄视频| 深夜精品福利| 久久精品夜色国产| av黄色大香蕉| 亚洲乱码一区二区免费版| 亚洲人与动物交配视频| 我的老师免费观看完整版| 美女内射精品一级片tv| 久久久久国产网址| 岛国在线免费视频观看| 国产精品99久久久久久久久| 18禁在线播放成人免费| 成人特级黄色片久久久久久久| 91在线精品国自产拍蜜月| 在线播放无遮挡| 99久久精品热视频| 能在线免费观看的黄片| 久久亚洲国产成人精品v| 国产三级中文精品| 久久久久九九精品影院| 99久久成人亚洲精品观看| av在线观看视频网站免费| 99热这里只有是精品在线观看| 国产av在哪里看| 久久久久性生活片| 欧美日韩综合久久久久久| 国产亚洲精品久久久com| 22中文网久久字幕| 午夜激情欧美在线| 无遮挡黄片免费观看| 精品午夜福利视频在线观看一区| 久久精品国产亚洲av香蕉五月| 国产aⅴ精品一区二区三区波| 国产精品人妻久久久影院| 精品久久久噜噜| 久久久成人免费电影| 直男gayav资源| 久久精品久久久久久噜噜老黄 | 又粗又爽又猛毛片免费看| 色视频www国产| 国内揄拍国产精品人妻在线| 国产 一区 欧美 日韩| 少妇的逼好多水| 国产精品不卡视频一区二区| 嫩草影院精品99| av在线亚洲专区| 我要看日韩黄色一级片| 国内精品一区二区在线观看| 成人无遮挡网站| 免费在线观看影片大全网站| 久久久久性生活片| 国内揄拍国产精品人妻在线| 国产色爽女视频免费观看| 久久久精品94久久精品| 国产免费男女视频| 美女内射精品一级片tv| 九九在线视频观看精品| 美女大奶头视频| 国产精品美女特级片免费视频播放器| 精品一区二区三区视频在线观看免费| 一区二区三区免费毛片| 悠悠久久av| 给我免费播放毛片高清在线观看| 精品人妻偷拍中文字幕| 极品教师在线视频| 国产私拍福利视频在线观看| 成人精品一区二区免费| 午夜亚洲福利在线播放| 青春草视频在线免费观看| 三级男女做爰猛烈吃奶摸视频| 97热精品久久久久久| 中文亚洲av片在线观看爽| 99riav亚洲国产免费| 精品一区二区三区人妻视频| 欧美xxxx性猛交bbbb| 婷婷色综合大香蕉| 日韩欧美国产在线观看| 哪里可以看免费的av片| 亚洲性久久影院| 夜夜爽天天搞| 国产高清视频在线播放一区| 久久久国产成人精品二区|