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

    線性回歸模型中非正態(tài)數據的處理

    2017-07-18 11:52:49丘甜華偉平李寶銀江希鈿
    武夷學院學報 2017年6期
    關鍵詞:正態(tài)因變量正態(tài)分布

    丘甜,華偉平,李寶銀,江希鈿

    (1.武夷學院商學院,福建武夷山354300;2.武夷學院生態(tài)與資源工程學院,福建武夷山354300;3.福建江夏學院,福建福州350108;4.福建農林大學林學院,福建福州350002)

    線性回歸模型中非正態(tài)數據的處理

    丘甜1,華偉平2,李寶銀3,江希鈿4

    (1.武夷學院商學院,福建武夷山354300;2.武夷學院生態(tài)與資源工程學院,福建武夷山354300;3.福建江夏學院,福建福州350108;4.福建農林大學林學院,福建福州350002)

    為了對非正態(tài)數據進行線性回歸分析,需要對非正態(tài)數據的處理方法進行研究。在Box-Cox變換的基礎上改進的雙冪變換是一種有效的處理方法。結合Matlab軟件給出了雙冪變換下線性回歸模型中參數的極大似然估計與最小二乘估計方法,并通過實例研究顯示:雙冪變換使非正態(tài)數據服從正態(tài)分布,對于異常數據的處理有一定的效果,是數據正態(tài)變換的理想工具。

    非正態(tài)數據;雙冪變換;線性回歸模型;極大似然估計;最小二乘估計

    線性回歸模型的因變量假定來源于正態(tài)分布的總體。在這一假定前提下,通常的做法是采用極大似然法或最小二乘法給出參數的估計。如果通過隨機抽樣的數據非正態(tài)分布,則會使顯著性檢驗程序是無效的,實際上估計出來的參數沒有意義,直接影響回歸分析結果。這意味著,當數據不滿足正態(tài)性假定時,可能拒絕實際上好的模型,大大地增加了統計推斷中所犯第一類錯誤的概率。因此,非正態(tài)數據的處理是應用線性回歸模型時需要解決的問題。

    1 處理非正態(tài)數據的常用方法

    1.1 Box-Cox變換

    在考察可觀測隨機因變量Y和自變量X之間的關系時,經常采用如下正態(tài)線性回歸模型[1-2]:

    其中β∈Rp為回歸系數,ε為不可觀測隨機誤差向量。易知,模型(1)實際上需要滿足Gauss-Markov條件:Y~Nn(Xβ,σ2In)。針對所獲得的數據,對其進行回歸診斷,若不滿足Gauss-Markov條件,研究者在統計推斷時有可能會增加犯錯誤的概率。常用的處理方式是對數據采取某種“治療”措施,其中,數據變換就是一種常用的處理非正態(tài)數據的辦法。至于采取何種變換更為有效,這取決于數據本身的特點,也成為了如今國內外統計學界研究的熱點問題之一。迄今為止,如下由觀測值Y到Y(λ)的Box-Cox變換[3](依賴于未知參數λ):

    式(2)是研究得最為透徹的變換。通過參數的選擇,達到對非正態(tài)數據的綜合治理。雖然Box-Cox變換具有許多優(yōu)點,但注意到

    這表明Box-Cox變換存在截斷問題,即當λ<0時,Y(λ)在-1/λ處右截斷,而λ>0時,Y(λ)是在-1/λ處左截斷的,只有在λ=0時,Y(λ)取值范圍是(-∞,+∞),從而認為Y經過變換(2)得到的f(Y,λ)是一組正態(tài)樣本是不正確的(除了λ=0這種情況)。

    1.2 雙冪變換

    為了克服Box-Cox變換中的截斷問題,諸多學者進行了研究。目前,Yang[4-5]引入了如下變換。

    變換(3)是在Box-Cox變換的基礎上進行修正而得到的一種新變換,被稱為雙冪變換。非正態(tài)數據通過雙冪變換后服從以下正態(tài)線性回歸模型:

    式中:X是已知n×p列滿秩設計陣;ε是不可觀測的n維隨機誤差向量;λ∈(-∞,+∞),β,σ2是未知參數。

    對于參數λ∈(-∞,+∞)均有g(Y,-λ)=g(Y,λ)成立,這意味著根據樣本Y=(Y1,…,Yn)T不能唯一確定模型(4)中的參數λ。因此,為了消除這種不確定性,將參數λ限制在區(qū)間[0,+∞)內。

    2 非正態(tài)數據處理下參數的估計

    2.1 極大似然估計

    由模型(4)可得未知參數λ,β,σ2基于Y=(Y1,…,Yn)T的對數似然函數為:

    其中J(Y,λ)表示變量Y(λ)=(Y1(λ),…,Yn(λ))T和Y=(Y1,…,Yn)T之間的變換Jacob行列式,其表達式為:

    當λ≥0已知時,由(5)利用熟知的正態(tài)線性回歸模型極大似然估計的結果可得知(β,σ2)的極大似然估計分別為:

    其中PX=(XTX)-1XT為一正交投影陣。將(7)式代入(5)中可得:

    這說明了參數λ∈[0,+∞)的極大似然估計是以下極值問題的解∈[0,+∞),而參數β∈Rp,σ2>0的極大似然估計分別為用Matlab中內嵌的函數fminbnd可方便地算出[6-7]。

    2.2 最小二乘估計

    最小二乘估計作為另一種參數估計方法,也經常被人們所運用。作為比較,本文同時考慮模型(4)中參數λ,β,σ2的最小二乘估計法。此時,相應的誤差平方和為:

    由此及熟知的線性回歸模型最小二乘估計的結果可得已知λ≥0時,β的最小二乘估計為:

    3 數值實驗

    在實際分析時,常常考慮多個變量,并且這些變量中有一個變量是特別關心的,稱為因變量Y,其他變量作為影響因變量的自變量(考慮三個自變量X1、X2和X3)。為了分析自變量對因變量的影響,按照隨機原則抽取了36個樣本數據,如表1所示。

    表1 樣本數據Table 1 Sample data

    3.1 非正態(tài)數據診斷

    數據非正態(tài)性診斷[8-10]的常用方法是K-S檢驗。該檢驗的原假設認為總體符合正態(tài)分布。在顯著性水平取0.05時,若檢驗的P值小于0.05,則否定原假設,認為總體呈現非正態(tài)分布。

    通過對因變量Y進行基本統計分析,統計結果(表2)顯示K-S檢驗的P值為0.031,小于0.05,說明數據總體不符合正態(tài)分布,并且偏態(tài)系數為0.912,表現出一定程度的右偏,這在圖1有更直觀的表現。通過統計軟件SPSS操作,并輸出結果??梢耘袛嘣撘蜃兞渴欠钦龖B(tài)數據。因此,在進行回歸分析前必須進行數據正態(tài)性處理。

    表2 因變量Y的基本統計量Table 2 Basic statistics of dependent variable Y

    3.2 數據變換及正態(tài)性檢驗

    采用雙冪變換對因變量Y進行修正得到Y(λ),用極大似然法估計出變換參數λ^ML為1.128 6,用最小二乘法估計λ^LS為0.883 4。對變換的效果進行正態(tài)性檢驗。通過軟件SPSS進行基本統計分析,得到變換后的偏度系數及K-S檢驗結果(表3),發(fā)現右偏程度有所減少,且變換值能順利通過K-S檢驗。

    表3 原數據的雙冪變換及正態(tài)分布檢驗結果Table 3 The results of the dual power transformation and normal distribution of the original data

    為更直觀的反映雙冪變換下兩種估計的正態(tài)效果,給出變換后的與的正態(tài)Q-Q圖,同變換前進行比較,見圖2。

    3.3 雙冪變換后線性回歸分析

    由于線性回歸分析[11]的思路是一致的,目的在于考察變量之間的數量關系,并通過一定的數學表達式即回歸方程將關系描述出來,進而確定自變量對因變量的影響程度。以通過極大似然估計下雙冪變換后線性回歸分析為例。以Y()為因變量,X1、X2和X3為自變量進行多元線性回歸分析,具體結果見表4、表5。

    圖2 原數據與雙冪變換后數據Q-Q圖比較Figure 2 Comparison of the original data and the dual power transformation data Q-Q

    表4 方差分析Table 4 Variance analysis

    從表4中得到,模型復相關系數達到0.998,而決定系數為0.996,取得了較好的擬合優(yōu)度。檢驗回歸方程的P值為0,小于顯著性水平0.05,通過了回歸方程的顯著性檢驗,即所得到的回歸方程有統計學意義。

    標準回歸系數的絕對值反映了影響因變量的程度,絕對值越大,則有越大的控制。由表5可看出影響程度從大到小依次為X3、X2、X1。回歸方程為:

    表5 參數估計Table 5 Parameter estimation

    回歸系數的顯著性檢驗是要檢驗自變量對因變量的影響是否顯著。由表5可知,在顯著性水平為0.05時,三個自變量X1、X2、X3的回歸系數顯著性水平t檢驗的P值都小于0.05,所以拒絕原假設,即認為這些回歸系數和0有顯著差異。

    4 小結

    當搜集回來的樣本數據不滿足正態(tài)分布時,用線性回歸模型分析所得到的結論是不準確的。因此,需要進行非正態(tài)數據的正態(tài)變換。目前,Box-Cox變換是一種研究的最為透徹的方法之一,但存在截斷問題。而雙冪變換能夠克服截斷問題。結合Matlab軟件給出了雙冪變換下線性回歸模型中參數的極大似然估計與最小二乘估計的求法。在實例分析中,采用統計軟件SPSS對數據診斷其正態(tài)性后,用雙冪變換對該數據進行正態(tài)變換,結果表明雙冪變換具有較強的正態(tài)變換能力,是非正態(tài)數據正態(tài)變換的理想工具。

    [1]王桂松,史建紅.線性模型引論[M].北京:科學出版社,2004:175-178.

    [2]王松桂.線性統計模型:線性回歸與方差分析[M].北京:高等教育出版社,1999:1-20.

    [3]BOX G E P,COX D R.An analysis of transformation[J].Journal of the Royal Statistical Society B,2012(26):211-252.

    [4]YANG Z L.A modified family of power transformations[J].Economics Letters,2006,92(1):14-19.

    [5]YANG Z L,Anthony F.Inference for general parametric functions in Box-Cox-type transformation models[J].Canadian Journal of Statistics,2008,36(2):301-319.

    [6]張學敏.Matlab基礎及應用[M].北京:中國電力出版社,2009:201-205.

    [7]趙芳芳,賈翔宇,許作良.CIR模型參數校準的極大似然法[J].統計與信息論壇,2015(9):3-7.

    [8]李曉暉,袁峰,白曉宇,等.典型礦區(qū)非正態(tài)分布土壤元素數據的正態(tài)變換方法對比研究[J].地理與地理信息科學,2010(6):102-105.

    [9]莊泓剛.基于非正態(tài)分布的動態(tài)金融波動性模型研究[D].天津:天津大學,2009.

    [10]焦璨,張敏強,黃慶均,等.非正態(tài)分布測量數據對克隆巴赫信度α系數的影響[J].應用心理學,2008(3):276-281.

    [11]劉兆君.伴隨置信度的線性回歸模型[J].統計與信息論壇,2015(7):3-7.

    (責任編輯:葉麗娜)

    Processing of Non-normal Data in Linear Regression M odel

    QIU Tian1,HUAWeiping2,LIBaoyin3,JIANG Xidian4
    (1.School of Business,Wuyi University,Wuyishan,Fujian 354300;2.School of Ecology Resource Engineering,Wuyi University,Wuyishan,Fujian 354300;3.Fujian Jiangxia University,Fuzhou,Fujian 350108;4.School of Forestry,Fujian Agriculture and Forestry University,Fuzhou,Fujian 350002)

    In order tomake a linear regression analysis on the non-normal data,it is necessary to study the non-normal data processingmethod.Based on the Box-Cox transform,the improved dual power transformation is an effectivemethod.Themaximum likelihood estimation and least square estimation of the parameters in the linear regression model are given by the Matlab software.The case studies show that the non-normal data is subject to normal distribution with the dual power transformation,which has a certain effect on the processing of abnormal data and is an ideal tool for the normal transformation of data.

    non-normal data;dual power transformation;linear regressionmodel;maximum likelihood estimate;least squares estimation

    O212

    A

    :1674-2109(2017)06-0053-05

    2017-02-16

    南平市科技計劃項目(N2014Z01);??蒲许椖浚╔LZ201401,XL201512S);福建省生態(tài)產業(yè)綠色技術重點實驗資助項目(WYKF2017-8)。

    丘甜(1988-),女,漢族,助教,主要從事空間計量模型的統計推斷研究。

    江希鈿(1958-),男,漢族,教授,主要從事生物數學建模研究。

    猜你喜歡
    正態(tài)因變量正態(tài)分布
    調整有限因變量混合模型在藥物經濟學健康效用量表映射中的運用
    中國藥房(2022年7期)2022-04-14 00:34:30
    適應性回歸分析(Ⅳ)
    ——與非適應性回歸分析的比較
    偏最小二乘回歸方法
    文理導航(2017年20期)2017-07-10 23:21:03
    基于對數正態(tài)分布的出行時長可靠性計算
    雙冪變換下正態(tài)線性回歸模型參數的假設檢驗
    統計與決策(2017年2期)2017-03-20 15:25:23
    正態(tài)分布及其應用
    基于泛正態(tài)阻抗云的諧波發(fā)射水平估計
    正態(tài)分布題型剖析
    半參數EV模型二階段估計的漸近正態(tài)性
    χ2分布、t 分布、F 分布與正態(tài)分布間的關系
    十八禁高潮呻吟视频| 亚洲国产精品国产精品| 国产日韩欧美在线精品| 成年人免费黄色播放视频| 日韩av在线免费看完整版不卡| 中文字幕色久视频| 在现免费观看毛片| tube8黄色片| 性高湖久久久久久久久免费观看| 大香蕉久久成人网| 一级片'在线观看视频| 黄色毛片三级朝国网站| 韩国精品一区二区三区| 女人久久www免费人成看片| 久久九九热精品免费| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区久久| 老熟女久久久| 最近中文字幕2019免费版| 免费在线观看黄色视频的| 91老司机精品| 亚洲国产成人一精品久久久| 成人免费观看视频高清| 老司机午夜十八禁免费视频| 男男h啪啪无遮挡| 悠悠久久av| 深夜精品福利| 1024视频免费在线观看| 中国国产av一级| 亚洲精品日本国产第一区| 乱人伦中国视频| 亚洲国产精品国产精品| 欧美日韩综合久久久久久| 国产亚洲精品久久久久5区| 国产精品一国产av| 精品欧美一区二区三区在线| 中文字幕色久视频| 日本欧美视频一区| 欧美国产精品va在线观看不卡| 成年人午夜在线观看视频| 黄网站色视频无遮挡免费观看| 亚洲国产成人一精品久久久| 美女中出高潮动态图| 久久久久精品人妻al黑| 中文字幕另类日韩欧美亚洲嫩草| 免费看十八禁软件| 欧美日韩福利视频一区二区| 国产成人欧美在线观看 | 黄色视频不卡| 男女床上黄色一级片免费看| 日韩精品免费视频一区二区三区| 热re99久久精品国产66热6| 亚洲国产精品成人久久小说| 在现免费观看毛片| 精品高清国产在线一区| 午夜免费鲁丝| 91字幕亚洲| 十八禁网站网址无遮挡| 18禁黄网站禁片午夜丰满| 国精品久久久久久国模美| 婷婷成人精品国产| 一二三四社区在线视频社区8| 一区二区三区精品91| 欧美精品av麻豆av| 老汉色∧v一级毛片| 高清视频免费观看一区二区| 久久99精品国语久久久| 久久午夜综合久久蜜桃| 少妇粗大呻吟视频| 久久女婷五月综合色啪小说| 一级片'在线观看视频| 国产精品久久久久久精品电影小说| 国产日韩欧美在线精品| 亚洲精品久久久久久婷婷小说| 国产片内射在线| 久久久久国产精品人妻一区二区| 国产又爽黄色视频| 久久久亚洲精品成人影院| 国产亚洲午夜精品一区二区久久| 激情视频va一区二区三区| 午夜福利一区二区在线看| 天天躁夜夜躁狠狠久久av| 男女边吃奶边做爰视频| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲,欧美精品.| a级毛片在线看网站| 免费av中文字幕在线| 国产精品免费视频内射| 蜜桃在线观看..| 免费高清在线观看日韩| 色网站视频免费| 国语对白做爰xxxⅹ性视频网站| 亚洲综合色网址| 成年女人毛片免费观看观看9 | 亚洲av男天堂| 亚洲国产精品成人久久小说| 精品人妻熟女毛片av久久网站| 午夜av观看不卡| 久久青草综合色| 久久久久久人人人人人| 国产一区二区激情短视频 | 精品久久久精品久久久| 国产主播在线观看一区二区 | 久久久久久久国产电影| 黄频高清免费视频| 中文字幕最新亚洲高清| 国产爽快片一区二区三区| 精品久久久久久久毛片微露脸 | 国产精品二区激情视频| 亚洲五月色婷婷综合| 热99国产精品久久久久久7| 精品高清国产在线一区| 黄片播放在线免费| 日韩免费高清中文字幕av| 国产爽快片一区二区三区| 亚洲成色77777| 亚洲精品乱久久久久久| 一级a爱视频在线免费观看| 亚洲精品一二三| 久久99精品国语久久久| 中文精品一卡2卡3卡4更新| 亚洲专区中文字幕在线| 悠悠久久av| 日本欧美视频一区| 九草在线视频观看| 亚洲精品日韩在线中文字幕| 两人在一起打扑克的视频| av一本久久久久| 欧美乱码精品一区二区三区| 中文字幕精品免费在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 一二三四在线观看免费中文在| 欧美激情 高清一区二区三区| 叶爱在线成人免费视频播放| 亚洲av电影在线进入| 9色porny在线观看| 国产成人a∨麻豆精品| 麻豆乱淫一区二区| 另类亚洲欧美激情| 狠狠婷婷综合久久久久久88av| 波多野结衣一区麻豆| 国产高清国产精品国产三级| 国产视频一区二区在线看| 亚洲第一av免费看| 日韩,欧美,国产一区二区三区| 久久免费观看电影| 一本—道久久a久久精品蜜桃钙片| 精品国产国语对白av| 国产野战对白在线观看| 午夜福利免费观看在线| 制服诱惑二区| 亚洲国产欧美在线一区| 国产免费一区二区三区四区乱码| 人妻 亚洲 视频| 久久精品亚洲av国产电影网| av电影中文网址| 成人午夜精彩视频在线观看| 久久这里只有精品19| 亚洲一区中文字幕在线| 欧美黄色片欧美黄色片| 亚洲综合色网址| 韩国高清视频一区二区三区| 精品国产一区二区三区久久久樱花| 一级毛片 在线播放| 老司机亚洲免费影院| 亚洲视频免费观看视频| 亚洲伊人久久精品综合| a 毛片基地| 99国产精品免费福利视频| 欧美人与性动交α欧美精品济南到| 亚洲av片天天在线观看| 午夜免费观看性视频| 亚洲黑人精品在线| 欧美中文综合在线视频| 好男人电影高清在线观看| 性色av乱码一区二区三区2| 男女高潮啪啪啪动态图| videos熟女内射| 久久精品国产亚洲av高清一级| 亚洲av电影在线观看一区二区三区| 成人国产一区最新在线观看 | 91精品伊人久久大香线蕉| 欧美精品高潮呻吟av久久| 国产精品 国内视频| 亚洲精品久久午夜乱码| 亚洲成人免费av在线播放| 色精品久久人妻99蜜桃| 亚洲国产欧美网| 80岁老熟妇乱子伦牲交| 菩萨蛮人人尽说江南好唐韦庄| 九色亚洲精品在线播放| 叶爱在线成人免费视频播放| av国产精品久久久久影院| av在线老鸭窝| 亚洲色图 男人天堂 中文字幕| 日日爽夜夜爽网站| 亚洲视频免费观看视频| 成人国产av品久久久| 9色porny在线观看| 男人添女人高潮全过程视频| 菩萨蛮人人尽说江南好唐韦庄| 五月天丁香电影| 亚洲成av片中文字幕在线观看| 性少妇av在线| 亚洲精品第二区| 免费少妇av软件| 亚洲九九香蕉| 欧美人与性动交α欧美精品济南到| 大片电影免费在线观看免费| 2021少妇久久久久久久久久久| 中文字幕人妻丝袜一区二区| 国产激情久久老熟女| 三上悠亚av全集在线观看| 大型av网站在线播放| 亚洲 欧美一区二区三区| 久久狼人影院| 最近手机中文字幕大全| 欧美在线黄色| 亚洲人成电影观看| 中文乱码字字幕精品一区二区三区| 亚洲成av片中文字幕在线观看| 成年女人毛片免费观看观看9 | 少妇裸体淫交视频免费看高清 | 中国美女看黄片| 欧美黄色淫秽网站| 国产高清视频在线播放一区 | 国产亚洲午夜精品一区二区久久| 丝袜喷水一区| 18禁黄网站禁片午夜丰满| 日韩 亚洲 欧美在线| 性高湖久久久久久久久免费观看| 18禁黄网站禁片午夜丰满| 国产一区二区在线观看av| 看十八女毛片水多多多| 亚洲欧美精品自产自拍| 狂野欧美激情性xxxx| 免费一级毛片在线播放高清视频 | 波多野结衣av一区二区av| 亚洲,欧美,日韩| 国产精品 国内视频| 国产激情久久老熟女| 一二三四在线观看免费中文在| 午夜免费观看性视频| 日韩一本色道免费dvd| 精品第一国产精品| 黄色 视频免费看| 18禁国产床啪视频网站| 精品福利永久在线观看| 日韩av免费高清视频| 午夜福利影视在线免费观看| 国产亚洲一区二区精品| 久久毛片免费看一区二区三区| 丝袜美腿诱惑在线| 日韩制服骚丝袜av| 亚洲精品av麻豆狂野| 欧美日韩亚洲高清精品| 看免费成人av毛片| 欧美人与性动交α欧美精品济南到| 中文字幕高清在线视频| 涩涩av久久男人的天堂| 欧美xxⅹ黑人| 亚洲综合色网址| 一边摸一边抽搐一进一出视频| 日韩精品免费视频一区二区三区| 麻豆av在线久日| 亚洲av电影在线观看一区二区三区| 超碰97精品在线观看| 在线看a的网站| 丝袜美足系列| 19禁男女啪啪无遮挡网站| www.熟女人妻精品国产| 精品一区二区三卡| 精品国产一区二区久久| 欧美 日韩 精品 国产| 亚洲国产欧美一区二区综合| 欧美人与性动交α欧美软件| 丝袜美腿诱惑在线| 午夜91福利影院| 亚洲,欧美,日韩| 国产精品 欧美亚洲| 精品人妻一区二区三区麻豆| 中文字幕人妻熟女乱码| 美女视频免费永久观看网站| 只有这里有精品99| 亚洲欧美激情在线| 男女下面插进去视频免费观看| 中文字幕高清在线视频| 亚洲中文av在线| 久久午夜综合久久蜜桃| 在线观看免费视频网站a站| 最新在线观看一区二区三区 | 纵有疾风起免费观看全集完整版| 一本综合久久免费| 19禁男女啪啪无遮挡网站| 久久午夜综合久久蜜桃| 国产一卡二卡三卡精品| 亚洲国产中文字幕在线视频| 国产精品香港三级国产av潘金莲 | 久久久欧美国产精品| 久久久久视频综合| 国产野战对白在线观看| 女人精品久久久久毛片| 精品亚洲成a人片在线观看| 一区在线观看完整版| 国产视频一区二区在线看| 国产精品久久久久久人妻精品电影 | 国产成人精品久久久久久| 丝袜美足系列| 精品久久久久久久毛片微露脸 | 欧美 亚洲 国产 日韩一| av有码第一页| 午夜福利视频精品| 国产免费现黄频在线看| 亚洲中文日韩欧美视频| 看免费成人av毛片| 日日摸夜夜添夜夜爱| 国产伦理片在线播放av一区| 中文字幕最新亚洲高清| 久久亚洲国产成人精品v| √禁漫天堂资源中文www| 一个人免费看片子| 日韩中文字幕视频在线看片| 免费少妇av软件| 午夜老司机福利片| 乱人伦中国视频| 操出白浆在线播放| 自线自在国产av| 中文字幕人妻熟女乱码| 亚洲精品国产av蜜桃| 成人黄色视频免费在线看| 亚洲五月色婷婷综合| 老鸭窝网址在线观看| 国产片特级美女逼逼视频| 天堂中文最新版在线下载| 王馨瑶露胸无遮挡在线观看| 国产精品一区二区在线观看99| 午夜免费鲁丝| 九草在线视频观看| 中文字幕高清在线视频| 国产高清不卡午夜福利| 久久精品国产综合久久久| 大香蕉久久网| 看免费av毛片| 五月开心婷婷网| 人人妻人人添人人爽欧美一区卜| 亚洲,欧美,日韩| 久久综合国产亚洲精品| 国产又色又爽无遮挡免| 久久ye,这里只有精品| 成年动漫av网址| 久久久久国产一级毛片高清牌| 午夜福利视频在线观看免费| 亚洲精品乱久久久久久| 国产真人三级小视频在线观看| 久久久久视频综合| 性色av一级| 两个人免费观看高清视频| 十八禁高潮呻吟视频| 51午夜福利影视在线观看| 操美女的视频在线观看| 国产1区2区3区精品| 七月丁香在线播放| 午夜福利影视在线免费观看| 午夜福利一区二区在线看| 亚洲色图 男人天堂 中文字幕| 日韩免费高清中文字幕av| 久久亚洲国产成人精品v| 桃花免费在线播放| 性色av一级| 大型av网站在线播放| 亚洲色图综合在线观看| 国产淫语在线视频| 精品熟女少妇八av免费久了| 国产精品久久久久久精品电影小说| 岛国毛片在线播放| 国产成人啪精品午夜网站| 亚洲av在线观看美女高潮| 婷婷色麻豆天堂久久| 一区二区三区乱码不卡18| 熟女av电影| 又粗又硬又长又爽又黄的视频| 国产无遮挡羞羞视频在线观看| 国产女主播在线喷水免费视频网站| 日韩视频在线欧美| 后天国语完整版免费观看| 久久久久精品国产欧美久久久 | 男人操女人黄网站| 亚洲精品日本国产第一区| 午夜福利免费观看在线| 一级毛片女人18水好多 | 久久久精品94久久精品| 蜜桃国产av成人99| 欧美在线一区亚洲| 咕卡用的链子| 久久精品国产亚洲av涩爱| 国产高清videossex| 精品国产一区二区三区四区第35| 日本av免费视频播放| www日本在线高清视频| 欧美激情极品国产一区二区三区| 多毛熟女@视频| 久久久精品区二区三区| 一二三四在线观看免费中文在| 男女床上黄色一级片免费看| 一本—道久久a久久精品蜜桃钙片| 好男人电影高清在线观看| 国产熟女欧美一区二区| 免费久久久久久久精品成人欧美视频| 国产视频一区二区在线看| 亚洲精品一二三| 女人精品久久久久毛片| 黄色 视频免费看| 亚洲国产欧美一区二区综合| 性色av一级| 99精品久久久久人妻精品| 丝瓜视频免费看黄片| 国产一区二区在线观看av| av天堂在线播放| 七月丁香在线播放| 午夜福利乱码中文字幕| 制服人妻中文乱码| 国产在线视频一区二区| 久久精品国产综合久久久| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久精品电影小说| 亚洲激情五月婷婷啪啪| 国产成人欧美| 精品一区二区三区四区五区乱码 | 久久中文字幕一级| 欧美国产精品va在线观看不卡| 满18在线观看网站| 国产精品99久久99久久久不卡| 超色免费av| 欧美日韩黄片免| 人人妻人人澡人人看| 涩涩av久久男人的天堂| 男女边吃奶边做爰视频| 久久av网站| 国产成人免费无遮挡视频| 99国产综合亚洲精品| av不卡在线播放| 久久天躁狠狠躁夜夜2o2o | 菩萨蛮人人尽说江南好唐韦庄| 91成人精品电影| 国产成人91sexporn| 久久国产精品男人的天堂亚洲| 高清av免费在线| 国产一卡二卡三卡精品| 91精品国产国语对白视频| 欧美日韩亚洲综合一区二区三区_| 久久久久网色| 成年人黄色毛片网站| 黑人巨大精品欧美一区二区蜜桃| 少妇的丰满在线观看| 80岁老熟妇乱子伦牲交| 亚洲,一卡二卡三卡| 免费高清在线观看日韩| 欧美日韩av久久| 国产精品成人在线| 精品人妻1区二区| 午夜激情久久久久久久| 蜜桃在线观看..| 高清不卡的av网站| 国产免费福利视频在线观看| 亚洲精品第二区| av福利片在线| 日本欧美视频一区| 精品一区二区三区av网在线观看 | 午夜免费成人在线视频| 国产一区亚洲一区在线观看| 欧美激情 高清一区二区三区| 国产成人免费无遮挡视频| 少妇裸体淫交视频免费看高清 | 精品亚洲成国产av| 欧美精品av麻豆av| 成人亚洲精品一区在线观看| www日本在线高清视频| 超色免费av| www.自偷自拍.com| 亚洲精品乱久久久久久| 欧美av亚洲av综合av国产av| 七月丁香在线播放| 亚洲成人免费av在线播放| 十分钟在线观看高清视频www| 999久久久国产精品视频| 午夜av观看不卡| 婷婷色av中文字幕| 亚洲一区二区三区欧美精品| 下体分泌物呈黄色| 成年美女黄网站色视频大全免费| 91精品国产国语对白视频| 人人妻人人爽人人添夜夜欢视频| 夫妻午夜视频| 午夜福利乱码中文字幕| 多毛熟女@视频| 午夜免费观看性视频| 国产男女内射视频| 亚洲国产看品久久| xxx大片免费视频| 精品亚洲成a人片在线观看| 国产国语露脸激情在线看| 999精品在线视频| 成人黄色视频免费在线看| 婷婷成人精品国产| 国产一区二区三区av在线| 亚洲中文字幕日韩| 亚洲精品日韩在线中文字幕| 亚洲五月色婷婷综合| 欧美激情高清一区二区三区| 精品高清国产在线一区| 捣出白浆h1v1| 丰满人妻熟妇乱又伦精品不卡| 波多野结衣av一区二区av| 精品久久久久久电影网| 国产无遮挡羞羞视频在线观看| 中国美女看黄片| 欧美精品啪啪一区二区三区 | 一级毛片我不卡| 黄色视频不卡| 久久综合国产亚洲精品| a级片在线免费高清观看视频| 久久久久精品国产欧美久久久 | 亚洲一区二区三区欧美精品| 欧美日韩av久久| 一级片免费观看大全| 女人爽到高潮嗷嗷叫在线视频| 亚洲人成77777在线视频| 精品国产国语对白av| 日韩,欧美,国产一区二区三区| 久久影院123| 久久精品aⅴ一区二区三区四区| 亚洲成国产人片在线观看| 18禁观看日本| 观看av在线不卡| 国产一级毛片在线| 国产一区二区三区av在线| 午夜两性在线视频| 国产精品 国内视频| 国产精品久久久久久精品电影小说| 日本a在线网址| 无限看片的www在线观看| 波野结衣二区三区在线| 一级毛片女人18水好多 | bbb黄色大片| 妹子高潮喷水视频| 91国产中文字幕| 丰满饥渴人妻一区二区三| 国产精品国产三级专区第一集| 97精品久久久久久久久久精品| 天天操日日干夜夜撸| 中文字幕精品免费在线观看视频| 视频区欧美日本亚洲| 真人做人爱边吃奶动态| 日韩熟女老妇一区二区性免费视频| 国精品久久久久久国模美| 看免费成人av毛片| av在线播放精品| 19禁男女啪啪无遮挡网站| 亚洲色图综合在线观看| 高清视频免费观看一区二区| 午夜免费鲁丝| 尾随美女入室| 91麻豆av在线| 新久久久久国产一级毛片| 日韩一卡2卡3卡4卡2021年| 国产在线一区二区三区精| 国产麻豆69| 极品少妇高潮喷水抽搐| 国产成人av教育| 国产爽快片一区二区三区| 一区在线观看完整版| 性色av一级| 欧美日韩成人在线一区二区| 建设人人有责人人尽责人人享有的| 欧美亚洲 丝袜 人妻 在线| 国产日韩欧美视频二区| 亚洲激情五月婷婷啪啪| 亚洲欧美日韩另类电影网站| 纵有疾风起免费观看全集完整版| 我要看黄色一级片免费的| 99国产精品一区二区蜜桃av | 精品久久久精品久久久| 老司机影院毛片| 我要看黄色一级片免费的| 99国产精品一区二区蜜桃av | 岛国毛片在线播放| 精品少妇久久久久久888优播| 日韩熟女老妇一区二区性免费视频| 中文乱码字字幕精品一区二区三区| 制服诱惑二区| 男人添女人高潮全过程视频| 亚洲成人手机| 国产精品成人在线| 国产精品香港三级国产av潘金莲 | 欧美精品人与动牲交sv欧美| 韩国精品一区二区三区| 久久热在线av| 99国产精品一区二区蜜桃av | 欧美国产精品一级二级三级| 美女脱内裤让男人舔精品视频| 国产一级毛片在线| 久久久精品94久久精品| 中文字幕人妻丝袜制服| 纯流量卡能插随身wifi吗| 国产高清国产精品国产三级| 久久精品国产亚洲av高清一级| 伊人久久大香线蕉亚洲五| 国产日韩一区二区三区精品不卡| 国产伦人伦偷精品视频| 男人舔女人的私密视频| 欧美大码av| 国产女主播在线喷水免费视频网站| 人人妻,人人澡人人爽秒播 | 热99久久久久精品小说推荐| 国产精品熟女久久久久浪|