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

    Pareto分布中尺度參數(shù)的區(qū)間估計

    2015-03-10 01:43:55徐付霞
    關鍵詞:樞軸樣本容量正態(tài)

    王 娟,徐付霞

    (天津工業(yè)大學 理學院,天津 300387)

    Pareto分布中尺度參數(shù)的區(qū)間估計

    王 娟,徐付霞

    (天津工業(yè)大學 理學院,天津 300387)

    研究了樞軸量法、極大似然估計的漸近正態(tài)性法和輪廓似然函數(shù)法,求解了Pareto分布中尺度參數(shù)的置信區(qū)間.通過生成隨機數(shù)進行數(shù)值模擬分析,得出區(qū)間估計和區(qū)間長度,對這幾種區(qū)間估計方法進行了比較.

    Pareto分布;樞軸量;極大似然;輪廓似然

    近年來,一些文獻已對比較常見分布的未知參數(shù)的區(qū)間估計問題進行了研究,解決了很多與之相關的實際問題,如正態(tài)分布中均值的區(qū)間估計和學生成績之間的關系;指數(shù)分布中,參數(shù)的區(qū)間估計和組件壽命之間的關系等[1].Pareto分布作為一種收入分布由意大利經濟學家Pareto提出,它不僅在經濟收入模型中得到應用,在其他領域中也得到了廣泛的應用.可以用它描述各種社會經濟、物理以及生物現(xiàn)象,例如城市人口、股票價格、商業(yè)失效、江河流量和某種藥理過后病人的存活時間等,而且在軍事、天文領域也得到了廣泛應用[2].在此之前,很多文章只研究了Pareto分布參數(shù)的多種點估計方法,而研究Pareto分布中參數(shù)的區(qū)間估計更有意義,便于人們對各種實際問題進行更加全面深入地探討.因此,本文給出Pareto分布中尺度參數(shù)的幾種區(qū)間估計方法,重點研究之前不常見的極大似然估計的漸近正態(tài)性法和輪廓似然函數(shù)法,并說明這幾種方法的適應范圍及優(yōu)缺點.

    1 預備知識

    1.1 樞軸量[3]

    一般來說,欲構造θ的置信水平為1-α的置信區(qū)間,我們首先考慮θ的極大似然估計MLE,或θ的充分統(tǒng)計量T=T(X),并基于T(X)尋找樞軸量,然后構造θ的置信區(qū)間.在統(tǒng)計量T(X)是連續(xù)型隨機變量時,設T(X)的分布函數(shù)為G(t;θ)=Pθ{T(X)≤t},那么G(T(x)t;θ)~U(0,1),U(0,1)為(0,1)上的均勻分布,此時G(T(X);θ)可被取為樞軸量.

    選取兩個常數(shù)c和d(0

    Pθ{c≤G(T(X);θ)≤d}=1-α

    (1)

    通過求解由式(1)確定的關于θ的不等式c≤G(T(X);θ)≤d就可得到θ的置信水平為1-α的置信區(qū)間.

    1.2 極大似然估計的漸近正態(tài)性[4]

    極大似然估計是使用極為廣泛的一種參數(shù)估計方法,使用該估計的最大優(yōu)點是能得到抽樣分布的漸近分布,從而得到標準誤和置信區(qū)間的近似值.設隨機變量X的分布函數(shù)為F∈{F(x;θ)∶θ∈Rd},(d≥1).X1,…,Xn是X的一個樣本,(·)和分別表示對數(shù)似然函數(shù)和參數(shù)向量θ=(θ1,…,θd)的極大似然估計,則在一定的正則條件下,當n充分大時,有

    其中“A”表示漸近分布,Nd表示d維正態(tài)分布,矩陣

    IE(θ)度量了對數(shù)似然函數(shù)的期望曲率,稱為Fisher信息陣.

    其中:Z1-α/2表示標準正態(tài)分布的1-α/2分位數(shù),以下同.

    (2)

    通常這個區(qū)間比用Fisher信息陣得到的區(qū)間更精確.

    1.3 輪廓似然函數(shù)[5]

    當總體分布含有兩個以上的參數(shù),而我們只對其中的某個或某幾個感興趣時,可以用輪廓似然函數(shù)方法求解其置信區(qū)間.這種方法可以避免未研究的參數(shù)對所研究參數(shù)的區(qū)間估計造成影響.

    設θ=(θ(1),θ(2)),θ(1)表示θ中感興趣的k維分量,θ(2)表示θ的其他d-k維分量,也稱為討厭參數(shù).θ(1)的輪廓似然函數(shù)定義為

    即對每個θ(1)的值,輪廓似然函數(shù)就是除θ(1)之外,似然函數(shù)對所有其他分量求最大值.在適當?shù)恼齽t條件下,當n充分大時,有

    Dp(θ(1))=2{

    Dp(θ(1))=2{p(θ(1))}稱為偏差度函數(shù)表示自由度為k的卡方分布.對于給定的顯著性水平α,有

    成立,由此可得θ(1)的置信水平為1-α的置信區(qū)間為

    (3)

    2 Pareto分布中尺度參數(shù)的區(qū)間估計

    設總體X服從參數(shù)為β,θ的Pareto分布,其概率密度為

    (4)

    其中:β為形狀參數(shù),θ為尺度參數(shù),或稱門限參數(shù).本文對尺度參數(shù)θ作區(qū)間估計,對于形狀參數(shù)β分兩種情況研究.

    2.1 β已知時,θ的區(qū)間估計

    我們可以用兩種方法得到θ的區(qū)間估計:

    2.1.1樞軸量法

    (5)

    其中I(θ,+∞)(x(1)是示性函數(shù),X(1)=min{X1,…,Xn},x(1)是X(1)的觀測值.

    使式(5)達到最大,就得到參數(shù)θ的極大似然估計

    . (6)

    下面求X(1)的分布函數(shù),由式(4)可得Pareto分布的分布函數(shù)為

    于是X(1)的密度函數(shù)和分布函數(shù)分布為

    任選取兩個常數(shù)c和d(0

    從而θ的置信水平為1-α的置信區(qū)間為

    區(qū)間長度為

    可見區(qū)間平均長度關于c遞增,取最小的c值,c=0,就得到平均長度最短的置信度為1-α的置信區(qū)間為

    (7)

    2.1.2極大似然估計的漸近正態(tài)性法

    相應的1階觀測信息矩陣為

    由1.2節(jié)的式(2)得θ的近似置信度為1-α的置信區(qū)間為

    (8)

    多數(shù)情況下,參數(shù)β事先未給定,這時β不能直接用于對參數(shù)θ的估計,我們可以采用以下方法求解θ的估計.

    2.2 當未β知時,θ的區(qū)間估計

    2.2.1極大似然估計的漸近正態(tài)性法

    (9)

    使式(9)取得最大值,就得到的極大似然估計

    (10)

    (11)

    由1.2節(jié)的式(2)得θ,β的置信度為1-α的置信區(qū)間分別為

    (12)

    2.2.2輪廓似然函數(shù)法

    由式(9)得θ的輪廓似然函數(shù)為

    對任意θ,欲使(θ,β)取最大值,需對(θ,β)關于β求導,并令其為0,即得

    (13)

    偏差度函數(shù)為

    Dp(θ)=2{

    則對給定的α,有

    成立,即

    (14)

    3 數(shù)值模擬分析

    3.1 β已知,對θ的區(qū)間估計的數(shù)值模擬

    假設β=0.5已知(為0.5),通過Matlab軟件生成n=16個服從參數(shù)θ=10的Pareto分布的隨機數(shù)如下[6]:

    32.425,145.863,1448.172,23.907,67.620,14.705,28.806,286.610,230.913,28.327,883.347,10.201,1635.349,2403.728,11.355,13.507

    進一步,選取不同的樣本容量n,生成服從參數(shù)β=0.5,θ=10的Pareto分布的隨機數(shù),得到兩種估計方法對應的置信區(qū)間,其長度分別為L1,L2,見表1.

    表1 β已知時,兩種估計方法的區(qū)間長度

    為了更直觀的看出L隨樣本容量n的變化趨勢,畫出n與兩種方法得到的區(qū)間長度之間的變化關系圖,如圖1(A).

    圖1 區(qū)間長度隨n的變化關系圖

    由圖1可見,樞軸量法得到的區(qū)間估計長度L1比極大似然估計的漸近正態(tài)性法得到的區(qū)間估計長度L2短.當樣本容量n較小時,樞軸量法和極大似然估計的漸近正態(tài)性法所得的區(qū)間估計長度都較大,隨著樣本容量n的逐步增大,兩種估計方法的區(qū)間長度逐漸減小,精度提高,兩者越來越接近且趨于平穩(wěn)(隨著的增大,圖中曲線斜率變化幅度越來越小).為了保證估計的精度,我們多數(shù)采用樞軸量法,但當樞軸量過于復雜難于求解,并且樣本容量足夠大時,極大似然估計的漸近正態(tài)性法也不失為一種好的區(qū)間估計方法.

    3.2 β未知,θ對的區(qū)間估計的數(shù)值模擬

    利用Matlab軟件,生成了n=16個服從參數(shù)β=0.1,θ=0.2的Pareto分布的隨機數(shù)為: 1.274×104,155.352,19.291,18.349,393.295,3.866×107,18.327,2.333, 1.076×1010, 2.209, 2.098×104, 1.847×106, 3.631, 8.923,47.963,722.575.

    再應用輪廓似然函數(shù)法得到偏差度函數(shù)

    對于給定的顯著性水平α有

    可見,當β未知時,通過上述兩種方法得到的θ的區(qū)間估計,輪廓似然函數(shù)法的區(qū)間長度較極大似然估計的漸近正態(tài)性法的區(qū)間長度h2短.為了進一步對兩種方法作出比較,現(xiàn)取不同的樣本容量n,得到兩種估計方法對應的區(qū)間長度,分別為h1,h2.如表2.

    表2 未知時,兩種估計方法的區(qū)間長度

    相應地,區(qū)間長度h隨樣本容量的變化關系如圖1(B) .

    [1] 蔣福坤, 劉正春. 指數(shù)分布參數(shù)的最短區(qū)間估計[J]. 數(shù)理統(tǒng)計與管理, 2004, 23(3): 43-45.

    [2] 李海芬.Pareto分布的統(tǒng)計分析. [D]. 上海: 華東師范大學, 2004.

    [3] 茆詩松, 王靜龍. 高等數(shù)理統(tǒng)計[M]. 北京: 高等教育出版社, 2004.

    [4]STUARTC.AnintroductiontoStatisticalModelingofExtremeValue[M].NewYork:Springer, 1990.

    [5] 史道濟. 實用極值統(tǒng)計方法[M]. 天津:天津科學技術出版社, 2005.

    [6]NELSENRB.Anintroductiontocopulas[M].NewYork:Springer, 1998.

    [7] 茆詩松 ,程依明. 概率論與數(shù)理統(tǒng)計教程[M]. 北京: 高等教育出版社, 2004.

    Interval estimates on shape parameter of Pareto distribution

    WANG Juan, XU Fu-xia

    (School of Science, Tianjin Polytechnic University, Tianjin 300387, China)

    In this paper, the methods of pivot quantity, the asymptotic normality of maximum likelihood estimate and the profile likelihood function were used to solve the confidence interval of the shape parameter of Pareto distribution. Calculated the interval estimation and interval length by generating more groups of random numbers and carry on the comparison to the interval estimation methods.

    Pareto distribution; pivot quantity; max likelihood; profile likelihood

    2015-02-27.

    王 娟(1990-),女,碩士,研究方向:極值統(tǒng)計.

    O212

    A

    1672-0946(2015)05-0629-05

    猜你喜歡
    樞軸樣本容量正態(tài)
    WK-35 電鏟中央樞軸液氮冷裝工藝研究
    面向神經機器翻譯的樞軸方法研究綜述
    探討參數(shù)區(qū)間估計中樞軸量的選取——以單個正態(tài)總體均值為例
    采用無核密度儀檢測壓實度的樣本容量確定方法
    雙冪變換下正態(tài)線性回歸模型參數(shù)的假設檢驗
    基于泛正態(tài)阻抗云的諧波發(fā)射水平估計
    半?yún)?shù)EV模型二階段估計的漸近正態(tài)性
    抽水蓄能電站球閥樞軸軸套故障分析及改造
    基于正態(tài)變換的貸款組合定價模型構建及實證
    技術經濟(2014年8期)2014-02-28 01:29:42
    廣義高斯分布參數(shù)估值與樣本容量關系
    測繪通報(2013年2期)2013-12-11 07:27:44
    国产真实伦视频高清在线观看| 中文欧美无线码| 美女黄网站色视频| 国产伦理片在线播放av一区| 欧美一区二区亚洲| 高清视频免费观看一区二区 | 欧美变态另类bdsm刘玥| 欧美另类一区| 九草在线视频观看| 欧美性猛交╳xxx乱大交人| 毛片女人毛片| 日韩不卡一区二区三区视频在线| 久久久精品免费免费高清| 51国产日韩欧美| 成人毛片60女人毛片免费| 校园人妻丝袜中文字幕| 国产老妇女一区| 亚洲不卡免费看| 精品久久久久久久久亚洲| 亚洲av中文字字幕乱码综合| 午夜免费激情av| 天美传媒精品一区二区| 18+在线观看网站| 成人无遮挡网站| 国产黄片视频在线免费观看| 国产高潮美女av| 亚洲av不卡在线观看| 日日撸夜夜添| 成人欧美大片| 搡女人真爽免费视频火全软件| 久久久久久久久中文| 99久久精品国产国产毛片| 午夜日本视频在线| 日本黄色片子视频| 精品人妻一区二区三区麻豆| 亚洲av一区综合| 国产老妇女一区| 乱系列少妇在线播放| 免费av不卡在线播放| 久99久视频精品免费| 精品熟女少妇av免费看| av女优亚洲男人天堂| 亚洲在久久综合| 好男人视频免费观看在线| 中文精品一卡2卡3卡4更新| 天天躁夜夜躁狠狠久久av| 国产成人午夜福利电影在线观看| 国产视频内射| 一区二区三区免费毛片| 3wmmmm亚洲av在线观看| 欧美最新免费一区二区三区| 啦啦啦啦在线视频资源| 2021天堂中文幕一二区在线观| 精品欧美国产一区二区三| 高清午夜精品一区二区三区| 少妇人妻精品综合一区二区| 真实男女啪啪啪动态图| 高清av免费在线| 亚洲成色77777| 最近的中文字幕免费完整| 美女主播在线视频| 狂野欧美白嫩少妇大欣赏| 久久久亚洲精品成人影院| 好男人视频免费观看在线| 午夜福利视频精品| 久久久久九九精品影院| av国产免费在线观看| 欧美高清性xxxxhd video| 一级a做视频免费观看| 国产一级毛片在线| 美女黄网站色视频| 啦啦啦啦在线视频资源| 赤兔流量卡办理| 大陆偷拍与自拍| 水蜜桃什么品种好| 嫩草影院新地址| 国产一级毛片七仙女欲春2| 亚洲国产成人一精品久久久| or卡值多少钱| 最近最新中文字幕大全电影3| 国产午夜福利久久久久久| 久久久久精品性色| 亚洲欧美成人精品一区二区| 亚洲欧美日韩卡通动漫| 草草在线视频免费看| www.av在线官网国产| 精品久久久久久成人av| 国产视频首页在线观看| 真实男女啪啪啪动态图| 欧美日韩精品成人综合77777| 亚洲精品日韩在线中文字幕| 成人二区视频| 精品久久久久久电影网| 天堂中文最新版在线下载 | 极品少妇高潮喷水抽搐| 高清av免费在线| 麻豆久久精品国产亚洲av| 小蜜桃在线观看免费完整版高清| 欧美成人精品欧美一级黄| or卡值多少钱| 国产精品福利在线免费观看| 精品国产露脸久久av麻豆 | 亚洲av免费高清在线观看| 男女边吃奶边做爰视频| 久久精品国产自在天天线| 国产午夜精品一二区理论片| 婷婷色av中文字幕| 日韩亚洲欧美综合| 一级毛片久久久久久久久女| 搡女人真爽免费视频火全软件| 亚洲不卡免费看| 听说在线观看完整版免费高清| 深夜a级毛片| 亚洲欧美精品自产自拍| av福利片在线观看| 色综合色国产| 中文资源天堂在线| 国产一级毛片在线| 简卡轻食公司| 国产免费视频播放在线视频 | 美女内射精品一级片tv| 国产 亚洲一区二区三区 | av专区在线播放| 欧美3d第一页| 成人特级av手机在线观看| 国产不卡一卡二| 国产毛片a区久久久久| 欧美成人精品欧美一级黄| 成人午夜高清在线视频| 色吧在线观看| 久久99精品国语久久久| 三级经典国产精品| 日本三级黄在线观看| 在现免费观看毛片| 午夜福利视频精品| 亚洲av电影在线观看一区二区三区 | 蜜桃亚洲精品一区二区三区| 国产成人a∨麻豆精品| 亚洲精品影视一区二区三区av| 国产三级在线视频| 亚洲av日韩在线播放| 国产白丝娇喘喷水9色精品| 国产又色又爽无遮挡免| 国产亚洲av嫩草精品影院| 男女国产视频网站| 99久久中文字幕三级久久日本| 婷婷色综合大香蕉| 亚洲精品乱码久久久久久按摩| 国产高清不卡午夜福利| 日韩一区二区三区影片| 午夜免费激情av| 欧美激情在线99| 嫩草影院精品99| 亚洲欧美日韩无卡精品| 亚洲第一区二区三区不卡| 日日摸夜夜添夜夜添av毛片| 成年av动漫网址| 麻豆成人午夜福利视频| 毛片一级片免费看久久久久| av福利片在线观看| 亚洲丝袜综合中文字幕| 成年女人在线观看亚洲视频 | 久久久久久九九精品二区国产| 看免费成人av毛片| 精品久久久久久久人妻蜜臀av| 九九爱精品视频在线观看| 三级国产精品片| 如何舔出高潮| 久久久久久久大尺度免费视频| 18禁在线播放成人免费| 黑人高潮一二区| 亚洲一级一片aⅴ在线观看| 欧美日韩精品成人综合77777| 亚洲精品日韩在线中文字幕| 边亲边吃奶的免费视频| 欧美三级亚洲精品| 国产精品.久久久| 高清毛片免费看| 亚洲精品第二区| 免费观看的影片在线观看| 在线免费十八禁| 国产精品久久久久久久电影| 欧美激情国产日韩精品一区| 国产探花极品一区二区| 小蜜桃在线观看免费完整版高清| 丝袜喷水一区| 欧美日韩亚洲高清精品| 国产中年淑女户外野战色| 天堂俺去俺来也www色官网 | 亚洲18禁久久av| 久久久成人免费电影| 亚洲精品乱码久久久久久按摩| 看十八女毛片水多多多| 国产成人一区二区在线| 黄色欧美视频在线观看| 丝瓜视频免费看黄片| 九九爱精品视频在线观看| av卡一久久| 少妇的逼好多水| 日韩大片免费观看网站| 九九久久精品国产亚洲av麻豆| 国产精品av视频在线免费观看| 亚洲欧美日韩无卡精品| 亚洲精品aⅴ在线观看| 国产乱人视频| 别揉我奶头 嗯啊视频| 日韩亚洲欧美综合| 亚洲精品第二区| 狠狠精品人妻久久久久久综合| 免费看日本二区| 亚洲最大成人手机在线| 成年人午夜在线观看视频 | 免费黄网站久久成人精品| 成人美女网站在线观看视频| 成人午夜精彩视频在线观看| 国产伦理片在线播放av一区| 91精品一卡2卡3卡4卡| 高清视频免费观看一区二区 | 岛国毛片在线播放| 小蜜桃在线观看免费完整版高清| 大又大粗又爽又黄少妇毛片口| 特大巨黑吊av在线直播| 建设人人有责人人尽责人人享有的 | 97精品久久久久久久久久精品| 久久久久性生活片| 美女内射精品一级片tv| videos熟女内射| 亚洲精品亚洲一区二区| 2021天堂中文幕一二区在线观| 亚洲av福利一区| 国语对白做爰xxxⅹ性视频网站| 国产亚洲一区二区精品| 丝袜喷水一区| 成人亚洲精品一区在线观看 | 97精品久久久久久久久久精品| 国产黄色视频一区二区在线观看| 人妻少妇偷人精品九色| 久久久久性生活片| av在线老鸭窝| 国产亚洲精品久久久com| 国产精品久久久久久精品电影小说 | 在线观看一区二区三区| 国产极品天堂在线| 亚洲av国产av综合av卡| 日韩av不卡免费在线播放| 老司机影院毛片| 精品一区二区三区人妻视频| 亚洲精品第二区| 午夜免费激情av| 国产成人aa在线观看| 中文字幕制服av| 91精品国产九色| 非洲黑人性xxxx精品又粗又长| 非洲黑人性xxxx精品又粗又长| 一本久久精品| 亚洲精品久久午夜乱码| 国产不卡一卡二| 国产精品av视频在线免费观看| 99热网站在线观看| 国产一区二区亚洲精品在线观看| 床上黄色一级片| 狂野欧美激情性xxxx在线观看| 久久久精品94久久精品| 夫妻午夜视频| 我的女老师完整版在线观看| 精品久久久久久电影网| 国产精品国产三级专区第一集| 免费观看在线日韩| 22中文网久久字幕| 纵有疾风起免费观看全集完整版 | 久久国内精品自在自线图片| 六月丁香七月| 在线观看美女被高潮喷水网站| 一夜夜www| 高清欧美精品videossex| 两个人视频免费观看高清| 国产黄a三级三级三级人| 九色成人免费人妻av| 一级毛片电影观看| 嫩草影院精品99| 国产综合精华液| 日本免费a在线| 91久久精品国产一区二区成人| 九色成人免费人妻av| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品专区欧美| 日韩精品有码人妻一区| 婷婷六月久久综合丁香| 精华霜和精华液先用哪个| 一级毛片黄色毛片免费观看视频| 国产成人精品久久久久久| 成人亚洲欧美一区二区av| 成年女人看的毛片在线观看| 最近中文字幕2019免费版| 午夜福利在线观看免费完整高清在| 国产一级毛片七仙女欲春2| 精品人妻偷拍中文字幕| 亚洲自拍偷在线| 在线观看一区二区三区| 欧美最新免费一区二区三区| 国产爱豆传媒在线观看| 国产一级毛片七仙女欲春2| 亚洲乱码一区二区免费版| 精品国产一区二区三区久久久樱花 | 国产一区二区亚洲精品在线观看| 草草在线视频免费看| 中国国产av一级| 亚洲四区av| 狂野欧美激情性xxxx在线观看| 亚洲av免费在线观看| 我要看日韩黄色一级片| 一本久久精品| 亚洲不卡免费看| 国产伦一二天堂av在线观看| 在线天堂最新版资源| 婷婷色综合大香蕉| 91久久精品国产一区二区三区| 18禁裸乳无遮挡免费网站照片| 成人鲁丝片一二三区免费| 国产一区有黄有色的免费视频 | 免费少妇av软件| 久久久久免费精品人妻一区二区| 天堂中文最新版在线下载 | 国产成人freesex在线| 午夜日本视频在线| 亚洲国产欧美在线一区| 欧美zozozo另类| 秋霞在线观看毛片| 国产成人a区在线观看| 大香蕉久久网| 久久久午夜欧美精品| 亚洲国产欧美在线一区| 久久午夜福利片| 婷婷色av中文字幕| 免费看不卡的av| 亚洲欧美日韩东京热| 我的老师免费观看完整版| 日韩制服骚丝袜av| 女人十人毛片免费观看3o分钟| 午夜激情福利司机影院| 中文字幕人妻熟人妻熟丝袜美| 免费av观看视频| 99久久九九国产精品国产免费| 国内揄拍国产精品人妻在线| 国产黄频视频在线观看| 最近中文字幕高清免费大全6| 亚洲av不卡在线观看| 欧美潮喷喷水| 人妻制服诱惑在线中文字幕| www.av在线官网国产| 人人妻人人澡欧美一区二区| 日韩av免费高清视频| 午夜免费激情av| av线在线观看网站| 国产成人freesex在线| 麻豆国产97在线/欧美| 久久久久精品久久久久真实原创| 国产精品三级大全| 午夜福利视频1000在线观看| 人妻少妇偷人精品九色| 一级av片app| 国产成人精品一,二区| 小蜜桃在线观看免费完整版高清| 身体一侧抽搐| 色综合亚洲欧美另类图片| 久久精品久久精品一区二区三区| 久久99热6这里只有精品| 亚洲人成网站在线观看播放| 免费观看在线日韩| 久久精品国产自在天天线| 亚洲精品久久午夜乱码| 国产乱来视频区| 在线播放无遮挡| 午夜福利网站1000一区二区三区| 97超视频在线观看视频| 能在线免费看毛片的网站| 国产91av在线免费观看| 人妻制服诱惑在线中文字幕| 久久99精品国语久久久| a级毛片免费高清观看在线播放| 国产精品久久久久久精品电影小说 | 亚洲国产日韩欧美精品在线观看| 国产一级毛片七仙女欲春2| 久久久久久久久久黄片| 久久午夜福利片| 精品亚洲乱码少妇综合久久| 岛国毛片在线播放| 成人国产麻豆网| 久久99热这里只频精品6学生| 麻豆成人午夜福利视频| 一级爰片在线观看| 天堂中文最新版在线下载 | 嘟嘟电影网在线观看| 欧美变态另类bdsm刘玥| 亚洲精品亚洲一区二区| 男女边吃奶边做爰视频| 久99久视频精品免费| 午夜视频国产福利| 日韩av不卡免费在线播放| av一本久久久久| 欧美xxxx性猛交bbbb| 亚洲丝袜综合中文字幕| 日韩精品青青久久久久久| xxx大片免费视频| 男人爽女人下面视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 在线播放无遮挡| 一区二区三区高清视频在线| 国产综合精华液| 国产综合懂色| 人人妻人人澡欧美一区二区| 一级毛片我不卡| 国产精品不卡视频一区二区| 国内揄拍国产精品人妻在线| 在线 av 中文字幕| 欧美高清成人免费视频www| 3wmmmm亚洲av在线观看| 毛片一级片免费看久久久久| 国产日韩欧美在线精品| 日韩精品青青久久久久久| 亚洲高清免费不卡视频| 欧美一级a爱片免费观看看| 亚洲精品日韩在线中文字幕| 嫩草影院精品99| 干丝袜人妻中文字幕| 欧美高清成人免费视频www| 蜜臀久久99精品久久宅男| 舔av片在线| 日本黄色片子视频| 亚洲av二区三区四区| 亚洲自偷自拍三级| 非洲黑人性xxxx精品又粗又长| 日本熟妇午夜| 国产精品一及| 亚洲av免费高清在线观看| 日韩av不卡免费在线播放| 午夜日本视频在线| 亚洲av电影在线观看一区二区三区 | 国产精品99久久久久久久久| av又黄又爽大尺度在线免费看| av在线蜜桃| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久久人妻蜜臀av| 亚洲av男天堂| 国产91av在线免费观看| 国产永久视频网站| 成人无遮挡网站| 色5月婷婷丁香| 久久久精品欧美日韩精品| 天堂影院成人在线观看| 国产精品无大码| 亚洲人与动物交配视频| 欧美日韩一区二区视频在线观看视频在线 | 精品一区二区三卡| 欧美日韩精品成人综合77777| 韩国高清视频一区二区三区| 国产有黄有色有爽视频| 中国国产av一级| 天堂影院成人在线观看| 亚洲乱码一区二区免费版| 亚洲精品国产av成人精品| 一个人看视频在线观看www免费| 亚洲av男天堂| 日本三级黄在线观看| 国产高潮美女av| 久久99热这里只频精品6学生| 国产乱人偷精品视频| 国产成人精品久久久久久| 综合色av麻豆| 亚洲精品国产av蜜桃| 亚洲国产日韩欧美精品在线观看| 草草在线视频免费看| 久久人人爽人人片av| 99视频精品全部免费 在线| 嫩草影院精品99| 一级毛片aaaaaa免费看小| 久久久成人免费电影| 中文字幕制服av| 街头女战士在线观看网站| 国产老妇伦熟女老妇高清| 亚洲av日韩在线播放| 成人av在线播放网站| 国产精品一区二区性色av| 最后的刺客免费高清国语| 成年版毛片免费区| 青春草亚洲视频在线观看| 天堂网av新在线| 国产一区二区在线观看日韩| 欧美xxxx性猛交bbbb| av.在线天堂| 国产av国产精品国产| 秋霞伦理黄片| 日本猛色少妇xxxxx猛交久久| 日本一本二区三区精品| 精华霜和精华液先用哪个| 一级黄片播放器| 日韩精品青青久久久久久| 能在线免费观看的黄片| 夜夜看夜夜爽夜夜摸| 亚洲国产最新在线播放| 伦理电影大哥的女人| 中国国产av一级| 国产激情偷乱视频一区二区| 亚洲内射少妇av| 老女人水多毛片| 看非洲黑人一级黄片| 亚洲精品成人av观看孕妇| 久久精品国产自在天天线| 国内精品宾馆在线| 白带黄色成豆腐渣| 一级片'在线观看视频| 黄色欧美视频在线观看| 最近最新中文字幕免费大全7| 蜜桃亚洲精品一区二区三区| 亚洲激情五月婷婷啪啪| 国产精品熟女久久久久浪| 国产成人精品福利久久| 亚洲av一区综合| 国产色婷婷99| 亚洲欧美成人综合另类久久久| 建设人人有责人人尽责人人享有的 | 亚洲国产精品国产精品| 亚洲精品视频女| 简卡轻食公司| 插逼视频在线观看| 一夜夜www| 天堂av国产一区二区熟女人妻| 国产 亚洲一区二区三区 | 一级黄片播放器| 永久免费av网站大全| 亚洲婷婷狠狠爱综合网| 成年版毛片免费区| 91精品伊人久久大香线蕉| 免费看光身美女| 激情 狠狠 欧美| 精品欧美国产一区二区三| 日韩国内少妇激情av| 国产成人福利小说| 免费黄频网站在线观看国产| 又大又黄又爽视频免费| 99热这里只有是精品在线观看| 91久久精品电影网| 精品久久久精品久久久| 人妻制服诱惑在线中文字幕| 91久久精品国产一区二区三区| 青春草视频在线免费观看| 国产成人免费观看mmmm| 亚洲久久久久久中文字幕| 亚洲精品日韩av片在线观看| 麻豆乱淫一区二区| 亚洲精品aⅴ在线观看| 国产在视频线精品| 亚洲成人av在线免费| 一区二区三区四区激情视频| 欧美人与善性xxx| 最近视频中文字幕2019在线8| 免费观看精品视频网站| xxx大片免费视频| 十八禁国产超污无遮挡网站| 久久久久久久久久黄片| 亚洲天堂国产精品一区在线| 五月伊人婷婷丁香| 亚洲高清免费不卡视频| 亚洲无线观看免费| 噜噜噜噜噜久久久久久91| 欧美潮喷喷水| 又粗又硬又长又爽又黄的视频| 国产亚洲av片在线观看秒播厂 | 欧美精品国产亚洲| 乱码一卡2卡4卡精品| 亚洲国产av新网站| 亚洲美女搞黄在线观看| 看非洲黑人一级黄片| 欧美日韩亚洲高清精品| 深夜a级毛片| 免费看av在线观看网站| 日日啪夜夜撸| 国产精品国产三级国产专区5o| av网站免费在线观看视频 | 国产高清三级在线| 国产精品一及| 老女人水多毛片| 在线免费观看的www视频| 老女人水多毛片| www.色视频.com| 精品人妻一区二区三区麻豆| 日本色播在线视频| 晚上一个人看的免费电影| 精品人妻偷拍中文字幕| 国产黄色视频一区二区在线观看| 成人毛片a级毛片在线播放| 亚洲国产精品成人久久小说| 精品熟女少妇av免费看| 日韩亚洲欧美综合| 精品一区二区三区人妻视频| 麻豆成人午夜福利视频| 久久久亚洲精品成人影院| 能在线免费观看的黄片| 中文天堂在线官网| 国产成人aa在线观看| 男女啪啪激烈高潮av片| 午夜激情福利司机影院| 1000部很黄的大片| 在线免费十八禁| 久久久久久九九精品二区国产| av女优亚洲男人天堂| 男的添女的下面高潮视频| 97精品久久久久久久久久精品| 神马国产精品三级电影在线观看| 亚洲人与动物交配视频| 汤姆久久久久久久影院中文字幕 | 欧美一级a爱片免费观看看| 欧美精品国产亚洲| 99久久精品一区二区三区| 成人漫画全彩无遮挡|