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

    一種基于單形空間缺失數據的補全方法

    2022-12-08 17:03:08李瑋琦
    現代計算機 2022年19期
    關鍵詞:共線性對數向量

    劉 冰,李瑋琦

    (達州職業(yè)技術學院人工智能學院,達州 635001)

    0 引言

    成分數據是一種具有比例結構的多維數據,其數學形式定義為

    則稱向量X為D維成分數據[1],向量空間SD為單形空間,其中X中的每一個元素xi表示其在整體中所占的比重。相較于普通數據而言,成分數據除了用于分析整體中各部分數據間的相對關系外,還有利于揭示普通數據所隱藏的相對信息,在諸如社會學、經濟學、氣象學、地質學、醫(yī)學等領域都有十分廣泛的應用[2-4]。

    由公式(1)可以看出,單形空間中的數據受到兩個條件的限制,一個是有界約束,另一個是定和約束,而現有一般的統(tǒng)計分析方法對被分析數據是沒有約束要求的。顯然,這就導致現有一般的統(tǒng)計分析方法在單形空間中無效,一個最主要的原因在于單形空間中數據的協方差矩陣通常是奇異矩陣,其含義與普通的數據不同[5];另一個原因在于單形空間中的數據總體一般不滿足多元正態(tài)分布的假設,這就會導致建立模型十分困難。為了能運用現有一般的統(tǒng)計分析方法去分析單形空間中的數據,通常的做法是先進行預處理,即經過一定的變換使之成為無約束數據。文獻[1]首次提出通過對成分數據的對數比變換建立成分數據的邏輯正態(tài)分布模型,有效地解決了現有一般的統(tǒng)計分析方法對單形空間中數據有界定和約束的限制,但建立的模型對數據的解釋性較差;文獻[6]提出了一種新變換,即對稱對數比變換,該方法能很好地解釋數據,但在某些情況下容易造成變換后的數據間存在較高的冗余性,從而損失部分數據信息;文獻[7-8]在此基礎上又進行了改進,即對稱等距對數比變換(isometric logratio transformations,ILR),但該方法對數據的完整性要求較高,若存在數據缺失,當補全為0時會造成變換后的數據為無窮的情況,顯然失去了實際的意義。對于成分數據的缺失值補全,目前一般采用單形空間均值(SM)補全法、極大似然補全法、期望最大填補法、k近鄰補全法(KNN),等等。這些方法對于回答信息而言其實現較為容易,但穩(wěn)健性差,結果偏差較大,補全后的數據冗余度高,缺少解釋性。

    消除變換后的數據冗余性的最有效方法是實現單形空間到歐氏空間的正交變換,為此,本文首先給出單形空間的代數運算體系,在給出文獻[4,9-11]的相關變換過程的基礎上,著重研究對于成分數據存在缺失值時通過對數比變換后存在的多重共線性數據的填補方法,并結合主成分分析法提出了一種較為有效的參考解決路徑。

    1 單形空間的代數體系

    為了實現單形空間到歐氏空間的正交變換,本文給出單形空間對于向量的加法、數乘、內積以及距離的如下定義。

    (1)向量的加法:對于任意X,Y∈SD,向量X,Y的加法運算⊕定義為

    式中,A[.]為封閉運算,即

    (2)向量的數乘:對于任意X∈SD,任意實數a∈R,a與向量X的數乘運算?定義為

    (3)向量的內積:對于任意X,Y∈SD,向量X,Y的內積定義為

    (4)向量的距離:對于任意X,Y∈SD,向量X,Y的Aitchison距離[10]定義為

    2 成分數據的缺失

    在實際工作中,待分析處理數據集中某些數據或屬性值缺失的原因是多方面的,或者是調查者基于主觀的判斷認為不重要而丟棄某些數據;或者是由于客觀的問卷設計存在瑕疵、錄入失誤、受訪者拒絕回答而沒能采集到某些數據;或者是在原始數據的存儲過程中,由于設備的故障造成存儲數據的不全或失敗而丟失某些數據等,使得沒能滿足設計預期獲得詳細而全面的資料數據,顯然,如果缺失數據占比較大,對于后續(xù)的數據分析處理會造成難以估計的影響。

    從式(1)易知,若一個成分向量中只有一個元素值缺失,即可根據定和限制求出該缺失值,因此一般來說,成分數據的缺失值是指某個樣本或屬性值中至少有兩個或兩個以上的缺失值。其數學形式化定義如下:

    若某個向量Xk(k= 1,2,…,n)中至少有兩個元素值存在缺失,則稱Z為缺失數據矩陣。

    3 缺失數據的補全

    由于主客觀等因素的影響,經常會碰到待分析處理的數據集中某個數據或某些屬性值出現為零或缺失的情況。對于前者,通常的做法是將其處理為缺失值;而對于后者,一般先要考慮缺失數據的占比情況,若某行(列)缺失數據比超過90%,一般進行剔除處理,或重新進行該行(列)數據的采集。對于缺失數據比小于90%的情況,則對缺失數據進行某種策略的填補?!N經典的填補方法是基于k近鄰(KNN)方法[12],即用通過某缺失值的k個最近鄰樣本信息來估算該缺失值;另一種是把缺失值當作一類隨機變量或者隱變量,建立概率隱變量模型,然后通過EM、VI(Variational Inference)或者MCI(Monte Carlo Inference)來估計缺失值的分布,具體做法是:

    對于式(1),進行如下的處理:

    其中,xOi為非缺失值數據,xMj為缺失值數據,則缺失值xMj的分布估計為

    然后計算該分布的期望值,并將其置為缺失值的估計值。

    對于簡單的模型,其解析解可用EM 算法求解;若模型復雜,則可借助MCI 去進行逼近求解,但無法解決結構帶來的不實用的問題。

    此外,對于多元線性回歸模型,若變量之間線性無關,還可采用回歸估計法對缺失值進行填補。但變量之間完全線性無關僅僅具有理論上可能,在實際情況下,變量之間往往存在多重共線性,若直接采用回歸估計法,其估計結果會與實際情況相去甚遠。

    對于多重補全法,文獻[9-11,13]給出的方法較有代表性,下面作簡要介紹:

    3.1 非對稱對數比變換

    對于式(1),定義如下變換:

    其逆變換式為xD=,進而有:

    從式(9)中可以看到,該變換是一個從單形空間SD到歐氏空間RD-1上的線性變換,而非正交變換,變換后的yi與變換前的xi不具有一一對應的關系,即存在非對稱關系,這就會導致建立的模型不能合理準確地解釋數據。

    3.2 對稱對數比變換

    針對非對稱對數比變換存在的缺陷,張堯庭[6]在《成分數據統(tǒng)計分析引論》中提出了一種新變換,使得變換后的yi與變換前的xi具有了一一對應的關系,即存在對稱關系,這就使得建立的模型具有了一定的可解釋性。其具體變換式如下:

    對于式(1),定義如下變換:

    其逆變換式為xD=,進而有:

    從式(10)看到,該變換是正交變換,但當0<xi<1 時,?[α1,α2,…,αD]T∈SD,αiyi≠0,即變換得到的數據存在一定的相關性,導致了變量間協方差矩陣不滿秩,從而使得基于協方差結構的統(tǒng)計方法無效,在實際應用中,應當避免使用該變換對成分數據進行預處理。

    3.3 對稱等距對數比變換

    文獻[4,11,13]又在對稱對數比變換基礎上進行了改進,即對稱等距對數比變換(isometric log-ratio transformations,ILR),具體如下:

    對于式(1),定義如下變換:

    容易得出:式(11)的逆變換式為

    進而有:

    從式(11)可以看出,該變換實現了從單形空間SD到歐氏空間RD-1的正交變換,確保了在變換后的空間中運用傳統(tǒng)的統(tǒng)計分析方法進行合理的模型建立。但在xi= 0時,對應的yi的結果將為無窮,失去了實際的意義,對后續(xù)的進一步分析處理造成了障礙。

    4 成分數據的補全

    一般情況下,對于缺失數據不宜貿然進行刪除處理,通常需要采用某種方法進行補全操作。常用的方法有:均值補全法、極大似然估計法、多重補全法等,其中多重補全法是通過估計出待補全的值加上不同的噪聲來得到補全值。對于成分數據缺失值的補全,Hron 等[12]提出的k 近鄰法較有代表性,該方法是通過用Aitchison 距離來尋找到含缺失值樣本的k 個近鄰,并用該k 個近鄰的中位數來進行初始補全,然后用最小二乘法來進行迭代補全。本文在前面定義的單形空間的加法運算⊕以及數乘運算?的基礎上結合文獻[12]的方法,提出一種基于單形空間缺失成分數據的補全方法,同時運用主成分分析法,處理將成分數據變換為一般數據后可能存在的多重共線性的情況。

    4.1 多重共線性

    在進行多元回歸分析時,若某些解釋變量之間存在嚴格或近似的線性關系,其樣本點或屬性值的一個微小改變都會極大地擾動回歸系數的估計值,使得回歸系數極不穩(wěn)定[14]。因為某些解釋變量之間存在的強相關關系將極大地降低ZTZ的可逆性,大多數情況變得不可逆,即使通過某種計算使其變得可逆,其逆矩陣的特征值也往往會較大,導致標準誤差值也較大,進而降低了參數估計值的精度,無法得出穩(wěn)定的回歸模型,回歸系數及符號也與實際情況相去甚遠。

    檢測多重共線性的方法主要有:

    (1)通過計算自變量間的相關系數與顯著性來進行判斷,即若某些變量間的相關系數顯著,則認為它們之間可能存在多重共線性問題。

    (2)使用回歸分析中的方差膨脹系數(Variance inflation factor,VIF)值來進行判斷,VIF的計算公式為VIF=1/ (1 -)。其中,Ri為負相關系數。自變量之間共線性程度與VIF 值存在較強的正相關關系。根據Hair(1995)標準,當VIF≤10 時,模型的多重共線性較弱;當10 <VIF≤100 時,模型的多重共線性較為嚴重;當VIF>100時,模型的多重共線性很嚴重。

    (3)容忍值(Tolerance)法,也是較為常用的方法。其計算公式為Tol=1/VIF。顯然,其與方差膨脹系數法的判定標準相反,自變量之間共線性程度與Tol值存在較強的負相關關系。在實際中,通常為Tol指定一個閾值,確保小于閾值的變量間的相關系數矩陣可逆,使回歸系數的估計值具有較強的穩(wěn)定性。該方法的缺陷在于Tol閾值的確定存在隨意性,沒有一個統(tǒng)一的標準。

    (4)主成分回歸法,對于矩陣(6),設ZTZ的特征值為λ1≥λ1≥… ≥λn>0,稱h=λ1/λn為ZTZ的條件數,一般地,若h<100,則認為模型的多重共線性程度較??;若100 <λ1<1000則認為模型的多重共線性程度較強;若h>1000,則認為模型的多重共線性程度嚴重。

    需要說明的是,在現實工作中,獲得的數據集一般都存在多重共線性,只是程度不同而已,對于共線性程度較小或一般的問題可以不必采取措施。另外,如果學得模型的擬合度好,也可不需處理多重共線性問題。

    4.2 單形空間上的均值補全

    根據公式(5)易知,若兩個樣本各自成分數據子集相似,則它們之間的Aitchison 距離可以用其對應子集的Aitchison 距離大約表示。即dA(xi,xj)≈dA(XMi,XMj)≈dA(XOi,XOj),其中 :XMi,XMj和XOi,XOj分別是樣本xi,xj各自所對應的缺失值和確定值成分數據子集。

    下面根據第2節(jié)的相關定義及文獻[12]的方法給出xi的某一缺失成分xmi∈XMi,m∈M的補全步驟:

    (1)根據Aitchison 距離找到含缺失值xi相應子成分XMi的k(k<n)個最近鄰, 并記為其對應的k個全樣本依次為。

    (2)根據定義1 和定義2,計算出k個全樣本的均值:

    (3)求出xmi的補全值:

    其中:Ij=(0,…,1,…,0)T∈Rn的第j個元素為1,j= 1,2,…,n。

    4.3 基于主成分的補全

    在大多數情況下,在單形空間上由子成分的Aitchison 距離對缺失值進行均值補全后的成分數據存在多重共線性,基于4.1 節(jié)所述,下面采用主成分回歸分析法對上節(jié)初始補全后的成分數據再次進行修正補全,主要步驟如下:

    (1)將含有缺失成分的樣本xi和其k個最近鄰樣本xi[1],xi[2],…,xi[k]組成一個單形空間矩陣,并將缺失值xmi初始補全后的變換到第1 行第1列,記為:A(k+1)×D。

    (2)根據公式(11),將單形空間矩陣轉換為歐式空間矩陣,如下所示:

    其中,α= irl(x?mi),A為一k×(D- 1)階矩陣,令:

    (3)對矩陣A做主成分分析,其協方差矩陣記為Λ=,其中,zu,zv為A的行向量,為A的行向量均值。

    (4)計算Λ的前p個主成分,依次為λ1≥λ1≥…≥λp≥0,則響應變量Y與間的關系為

    其中:m= 1,2,…,k,ε為誤差項。

    (5)通過上式得到βj的估計值,計算到缺失值xmi的補全值為

    (6)運用公式(12)將數據xmi還原為成分數據,并通過第(1)步將其調回到原始位置。

    5 評價與比較

    5.1 評價指標

    為了有效地評價上述方法對數據集中缺失值的補全效果,本文采用正規(guī)化方均根差(the normalized root mean squares error,NRMSE)作為判別準則,即:

    其中,vg為補全值,vr為真實值,μ(.)為均值,σ(.)為方差。NRMSE值的大小反映了真實值與補全值之間差距,若NRMSE值較大,則說明補全值與真實值存在較大差距;若NRMSE接近于0,則說明補全值非常接近真實值。

    5.2 比較分析

    為了驗證前述方法的有效性,選用文獻[1]中Hongite 數據,該數據集包含25個樣本,每個樣本包含5 個特征:ablite,blandite,cornite,daubite,endite,根據4.1 節(jié)所述計算得到條件數h=2747.238>>1000,即認為該數據集存在嚴重的多重共線性。下面假定該數據集的ablite 和cornite 特征數據缺失,分別運用k近鄰補全法(KNN),單形空間均值(SM)和主成分補全法(PCA)對缺失值進行補全,得到比較結果見表1。

    表1 KNN,SM和PCA補全操作比較結果表

    從表1可以看出,當條件數h>>1000時,用PCA 方法進行補全的結果最好,KNN 的補全結果最差。

    根據文獻[10]的結論模擬100 個5 維的成分數據x~N5φ(μ,∑),其中μ=(0,0,0,0)T,∑是一個主對角線上全為1,其余全為p的4 階方陣。在假定p的取值分別為0.3、0.7、0.995,缺失率(MR)分別為10%、20%和30%情況下,分別運用KNN、SM 和PCA 方法進行缺失值補全,并用NRMSE進行評價比較,結果如圖1所示。

    圖1 MR與P分別取值時三種方法的補全比較結果

    其中,在圖1 中的圖a1~a3 是在MR不同p一定時三種方法的補全比較結果,圖b1~b3 是在p不同MR一定時三種方法的補全比較結果。

    圖1 MR與p分別取值時三種方法的補全比較結果(續(xù))

    從圖b1 與圖b2 中可以看出,PCA 比KNN,SM 的結果都要好。這三種補全法在MR一定時,p與NRMSE呈負相關關系,也就是說若數據間的多重共線性程度越大,無論哪種方法的補全效果都越好。而在p一定時,三種方法的MR與NRMSE呈正相關關系。作為初始的補全法,KNN 法明顯比SM 差,隨著MR的増大,結果會更差。但在MR變大時,PCA 法明顯比SM 效果好,隨著MR的増大,結果會更明顯。

    6 結語

    基于單行空間完備的代數體系提出的等距對數比變換是一個正交變換,該變換既克服了非對稱對數比變換改變內積及距離等幾何概念的缺陷,同時,又避免了對數比變換導致的多重共線性給多元分析方法帶來的不利影響。對于含有缺失值的多元數據來說,無論是基于模型還是基于距離,多變量補全法比單變量補全法結果更為準確:在單形空間上先進行均值補全,然后運用等距對數比對補全后的數據進行變換,最后再對變換后的數據運用主成分法進行第二次補全,實例分析表明,再次運用主成分法進行二次補全要比其他方法的效果更好。

    猜你喜歡
    共線性對數向量
    向量的分解
    含有對數非線性項Kirchhoff方程多解的存在性
    指數與對數
    聚焦“向量與三角”創(chuàng)新題
    銀行不良貸款額影響因素分析
    科學與財富(2021年3期)2021-03-08 10:56:02
    指數與對數
    文氏圖在計量統(tǒng)計類課程教學中的應用
    ——以多重共線性內容為例
    不完全多重共線性定義存在的問題及其修正建議
    對數簡史
    向量垂直在解析幾何中的應用
    男人添女人高潮全过程视频| 丁香六月天网| 久久人人97超碰香蕉20202| 乱人伦中国视频| 国产亚洲精品第一综合不卡| 一个人免费看片子| 在线观看一区二区三区激情| 国产精品二区激情视频| 交换朋友夫妻互换小说| 精品久久蜜臀av无| 色婷婷久久久亚洲欧美| 极品少妇高潮喷水抽搐| 夫妻性生交免费视频一级片| 精品一区二区三区四区五区乱码 | 国产99久久九九免费精品| 精品少妇一区二区三区视频日本电影 | 在线免费观看不下载黄p国产| 女人高潮潮喷娇喘18禁视频| 黑丝袜美女国产一区| 只有这里有精品99| 国产99久久九九免费精品| 妹子高潮喷水视频| 国产伦理片在线播放av一区| 国产 一区精品| 制服人妻中文乱码| 亚洲成人一二三区av| 熟女av电影| 视频在线观看一区二区三区| 操美女的视频在线观看| 亚洲欧美清纯卡通| avwww免费| 天天添夜夜摸| 亚洲av福利一区| 最近手机中文字幕大全| 天天影视国产精品| 电影成人av| 色吧在线观看| 国产精品久久久久成人av| 女人高潮潮喷娇喘18禁视频| 成人三级做爰电影| 亚洲色图综合在线观看| 一本久久精品| 黄色视频在线播放观看不卡| 国产亚洲av高清不卡| 国产亚洲av高清不卡| 午夜精品国产一区二区电影| 在线观看免费视频网站a站| 看免费成人av毛片| 亚洲精品aⅴ在线观看| 亚洲av综合色区一区| 精品第一国产精品| 精品一区二区三区av网在线观看 | 99国产综合亚洲精品| 亚洲国产欧美一区二区综合| 国产黄频视频在线观看| 精品国产一区二区久久| 午夜激情久久久久久久| 亚洲av电影在线观看一区二区三区| 啦啦啦视频在线资源免费观看| 永久免费av网站大全| 婷婷色麻豆天堂久久| 亚洲欧洲国产日韩| xxxhd国产人妻xxx| 婷婷成人精品国产| 十八禁网站网址无遮挡| 国产精品无大码| 日本av手机在线免费观看| 婷婷色av中文字幕| 精品国产一区二区三区久久久樱花| 国产一区二区在线观看av| 国产免费一区二区三区四区乱码| 国产深夜福利视频在线观看| 男女免费视频国产| 国产激情久久老熟女| 欧美日韩av久久| 考比视频在线观看| 最黄视频免费看| 丁香六月欧美| 丁香六月天网| 制服丝袜香蕉在线| 久久影院123| 国产精品国产三级国产专区5o| 国产一卡二卡三卡精品 | 久久精品亚洲熟妇少妇任你| 赤兔流量卡办理| 伊人亚洲综合成人网| 99热网站在线观看| 各种免费的搞黄视频| 七月丁香在线播放| 久久久久国产一级毛片高清牌| 欧美另类一区| 人人澡人人妻人| 啦啦啦视频在线资源免费观看| 日本欧美国产在线视频| 晚上一个人看的免费电影| 巨乳人妻的诱惑在线观看| 尾随美女入室| 亚洲欧美精品自产自拍| 99精品久久久久人妻精品| av又黄又爽大尺度在线免费看| 久久久久久久久免费视频了| 婷婷色综合www| 狂野欧美激情性bbbbbb| 嫩草影视91久久| avwww免费| 国产精品av久久久久免费| 性少妇av在线| 亚洲欧美一区二区三区国产| 亚洲国产精品999| 一区在线观看完整版| 91精品国产国语对白视频| 国产深夜福利视频在线观看| 一级片免费观看大全| 中文字幕最新亚洲高清| 久久人人97超碰香蕉20202| 丝袜美腿诱惑在线| 国产片内射在线| 免费看av在线观看网站| 免费观看av网站的网址| 成人漫画全彩无遮挡| 麻豆精品久久久久久蜜桃| 亚洲精品美女久久久久99蜜臀 | 啦啦啦在线免费观看视频4| 亚洲专区中文字幕在线 | 国产精品免费视频内射| 国产在线视频一区二区| 国产一区二区在线观看av| 日本欧美国产在线视频| 亚洲精品美女久久久久99蜜臀 | 9191精品国产免费久久| 日本欧美视频一区| a级毛片在线看网站| 国产一区二区 视频在线| 中文字幕亚洲精品专区| 最近中文字幕2019免费版| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美日韩另类电影网站| 90打野战视频偷拍视频| 男女边吃奶边做爰视频| 成年av动漫网址| 一级毛片黄色毛片免费观看视频| 成人18禁高潮啪啪吃奶动态图| 999精品在线视频| 成人免费观看视频高清| www.熟女人妻精品国产| 亚洲精品一区蜜桃| av在线观看视频网站免费| 观看美女的网站| 午夜精品国产一区二区电影| av网站免费在线观看视频| 久久久久久久久免费视频了| av在线app专区| 最近最新中文字幕大全免费视频 | 一边摸一边做爽爽视频免费| 国产精品二区激情视频| 亚洲综合色网址| 国产精品免费视频内射| 90打野战视频偷拍视频| 日本猛色少妇xxxxx猛交久久| 久久久精品区二区三区| 免费女性裸体啪啪无遮挡网站| 一级毛片电影观看| 精品亚洲乱码少妇综合久久| 亚洲精品久久久久久婷婷小说| av在线播放精品| 亚洲少妇的诱惑av| 十八禁人妻一区二区| 中文字幕av电影在线播放| 91成人精品电影| 1024视频免费在线观看| 丁香六月欧美| 亚洲人成网站在线观看播放| 菩萨蛮人人尽说江南好唐韦庄| 国产片特级美女逼逼视频| 黄片小视频在线播放| 亚洲av成人精品一二三区| 日本欧美国产在线视频| 国产成人精品福利久久| 欧美中文综合在线视频| 精品亚洲乱码少妇综合久久| 亚洲av国产av综合av卡| 中文字幕高清在线视频| 久久狼人影院| 一区在线观看完整版| 又黄又粗又硬又大视频| 18禁国产床啪视频网站| 性少妇av在线| 亚洲欧美日韩另类电影网站| 你懂的网址亚洲精品在线观看| 男女高潮啪啪啪动态图| 波多野结衣av一区二区av| 亚洲人成电影观看| 在现免费观看毛片| 欧美日本中文国产一区发布| 美女扒开内裤让男人捅视频| 自线自在国产av| 看免费成人av毛片| 精品久久蜜臀av无| 在线亚洲精品国产二区图片欧美| 18在线观看网站| 亚洲成人手机| 亚洲欧美成人精品一区二区| av又黄又爽大尺度在线免费看| 婷婷色av中文字幕| 大话2 男鬼变身卡| 少妇猛男粗大的猛烈进出视频| 日韩电影二区| 日本vs欧美在线观看视频| 韩国av在线不卡| 久久久久久人人人人人| 人人澡人人妻人| 97精品久久久久久久久久精品| 女人高潮潮喷娇喘18禁视频| 久久国产精品大桥未久av| 男人操女人黄网站| 中国国产av一级| 日韩 亚洲 欧美在线| 国产又爽黄色视频| 大陆偷拍与自拍| 看非洲黑人一级黄片| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区精品视频观看| 波多野结衣av一区二区av| 90打野战视频偷拍视频| 91国产中文字幕| 成人18禁高潮啪啪吃奶动态图| 激情视频va一区二区三区| 久久久久精品人妻al黑| 一区二区av电影网| 日本av手机在线免费观看| 日本欧美视频一区| 又大又黄又爽视频免费| 国产精品99久久99久久久不卡 | 性少妇av在线| 晚上一个人看的免费电影| 国产女主播在线喷水免费视频网站| 国产1区2区3区精品| avwww免费| 午夜老司机福利片| 免费看av在线观看网站| 又大又爽又粗| 久久久久精品性色| 一级毛片黄色毛片免费观看视频| 国产成人a∨麻豆精品| 国产又色又爽无遮挡免| av网站在线播放免费| 狠狠精品人妻久久久久久综合| 男的添女的下面高潮视频| 亚洲成人手机| 国产成人精品福利久久| 国产成人av激情在线播放| 欧美在线黄色| 日本wwww免费看| 久久人人爽人人片av| 狠狠婷婷综合久久久久久88av| 色视频在线一区二区三区| 久久精品熟女亚洲av麻豆精品| av福利片在线| 亚洲成人一二三区av| 最近最新中文字幕大全免费视频 | 国产深夜福利视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 天堂俺去俺来也www色官网| 男女免费视频国产| 最近手机中文字幕大全| 久久精品人人爽人人爽视色| 亚洲情色 制服丝袜| 日本欧美视频一区| 18禁裸乳无遮挡动漫免费视频| 欧美 日韩 精品 国产| 日韩av在线免费看完整版不卡| 亚洲综合精品二区| 黑丝袜美女国产一区| 97在线人人人人妻| 日韩精品有码人妻一区| 国产亚洲av高清不卡| 国产免费福利视频在线观看| 一本久久精品| 色网站视频免费| av国产精品久久久久影院| 人妻 亚洲 视频| 亚洲av男天堂| 久久精品国产a三级三级三级| 在线观看免费日韩欧美大片| 狂野欧美激情性bbbbbb| 欧美精品人与动牲交sv欧美| 久久精品亚洲熟妇少妇任你| 欧美激情高清一区二区三区 | 男女之事视频高清在线观看 | 亚洲国产精品成人久久小说| 成人亚洲欧美一区二区av| 久久青草综合色| 国产精品国产av在线观看| 少妇的丰满在线观看| 狠狠精品人妻久久久久久综合| 亚洲精品,欧美精品| 丝袜脚勾引网站| 国产有黄有色有爽视频| 女人爽到高潮嗷嗷叫在线视频| 黄色视频不卡| 日本爱情动作片www.在线观看| 国产午夜精品一二区理论片| 欧美少妇被猛烈插入视频| 欧美精品高潮呻吟av久久| 日韩精品有码人妻一区| 狠狠精品人妻久久久久久综合| 国产黄频视频在线观看| 建设人人有责人人尽责人人享有的| 国产福利在线免费观看视频| 日本欧美国产在线视频| 成年美女黄网站色视频大全免费| 亚洲综合精品二区| 久久久国产精品麻豆| 天美传媒精品一区二区| 在线精品无人区一区二区三| 国产亚洲午夜精品一区二区久久| 一区福利在线观看| 国产亚洲欧美精品永久| 国产色婷婷99| 大片电影免费在线观看免费| 日韩电影二区| 欧美中文综合在线视频| 久久精品aⅴ一区二区三区四区| 十分钟在线观看高清视频www| 成年av动漫网址| 日本色播在线视频| 老汉色∧v一级毛片| 大香蕉久久成人网| 亚洲三区欧美一区| 欧美精品一区二区免费开放| 热re99久久精品国产66热6| 亚洲精品自拍成人| 亚洲五月色婷婷综合| 母亲3免费完整高清在线观看| 无限看片的www在线观看| 欧美日韩综合久久久久久| 男男h啪啪无遮挡| 亚洲精品美女久久久久99蜜臀 | 亚洲成人一二三区av| 欧美亚洲 丝袜 人妻 在线| 97人妻天天添夜夜摸| 亚洲七黄色美女视频| 国产在视频线精品| 51午夜福利影视在线观看| 亚洲国产最新在线播放| 人人妻人人澡人人爽人人夜夜| 啦啦啦视频在线资源免费观看| 国产在线视频一区二区| 啦啦啦在线免费观看视频4| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| 亚洲成人国产一区在线观看 | 久久久精品区二区三区| 2021少妇久久久久久久久久久| 欧美黑人欧美精品刺激| 91老司机精品| 欧美人与善性xxx| 黑人猛操日本美女一级片| 男女无遮挡免费网站观看| 啦啦啦在线观看免费高清www| tube8黄色片| 女人精品久久久久毛片| 免费少妇av软件| 亚洲精品乱久久久久久| 精品一区二区三卡| 欧美少妇被猛烈插入视频| 午夜福利一区二区在线看| 欧美中文综合在线视频| 亚洲视频免费观看视频| 中国国产av一级| 18禁国产床啪视频网站| 日韩成人av中文字幕在线观看| 男女床上黄色一级片免费看| 久久久久久久精品精品| 精品酒店卫生间| 中文字幕av电影在线播放| 国产在线免费精品| 王馨瑶露胸无遮挡在线观看| 精品国产乱码久久久久久男人| 亚洲精品国产av成人精品| 国产熟女欧美一区二区| 不卡av一区二区三区| 天堂中文最新版在线下载| 免费日韩欧美在线观看| 纵有疾风起免费观看全集完整版| 777久久人妻少妇嫩草av网站| 天天躁夜夜躁狠狠久久av| 一级a爱视频在线免费观看| 久久久久久久久久久久大奶| 亚洲熟女毛片儿| 国产精品偷伦视频观看了| 国产精品一二三区在线看| 久久精品亚洲av国产电影网| 999精品在线视频| 免费观看a级毛片全部| 如何舔出高潮| 日韩不卡一区二区三区视频在线| 人人妻,人人澡人人爽秒播 | 国产精品女同一区二区软件| 精品人妻熟女毛片av久久网站| 一区二区三区四区激情视频| 老司机影院毛片| netflix在线观看网站| 美女脱内裤让男人舔精品视频| 夜夜骑夜夜射夜夜干| 欧美激情高清一区二区三区 | 街头女战士在线观看网站| 亚洲精品在线美女| 亚洲欧美成人综合另类久久久| 亚洲美女搞黄在线观看| 亚洲精品乱久久久久久| 久久久精品免费免费高清| 久久久久人妻精品一区果冻| 久久鲁丝午夜福利片| 中文字幕av电影在线播放| 亚洲人成77777在线视频| 最近手机中文字幕大全| 国产成人精品久久二区二区91 | 中文字幕另类日韩欧美亚洲嫩草| 曰老女人黄片| 只有这里有精品99| 男人舔女人的私密视频| 一本色道久久久久久精品综合| 精品国产一区二区三区四区第35| 中文字幕最新亚洲高清| 黄色视频不卡| 亚洲欧美清纯卡通| 色精品久久人妻99蜜桃| 国产精品久久久久久精品古装| 少妇被粗大猛烈的视频| 亚洲成人免费av在线播放| 国产成人欧美在线观看 | a级毛片在线看网站| 亚洲三区欧美一区| bbb黄色大片| 一个人免费看片子| 男女国产视频网站| 成人午夜精彩视频在线观看| 亚洲综合色网址| 丁香六月天网| 久久久精品国产亚洲av高清涩受| 建设人人有责人人尽责人人享有的| 亚洲av日韩在线播放| 精品卡一卡二卡四卡免费| 五月开心婷婷网| 久久影院123| 欧美在线黄色| 色吧在线观看| 久久久久久人人人人人| 大香蕉久久成人网| 欧美日韩精品网址| 中文字幕另类日韩欧美亚洲嫩草| 精品第一国产精品| 亚洲精品一二三| 美女扒开内裤让男人捅视频| 两个人看的免费小视频| 丝袜美足系列| 高清视频免费观看一区二区| 亚洲国产欧美在线一区| 人人妻,人人澡人人爽秒播 | 国产一卡二卡三卡精品 | 亚洲三区欧美一区| 免费人妻精品一区二区三区视频| 亚洲综合精品二区| 国产在线免费精品| 人人妻人人添人人爽欧美一区卜| 日韩一区二区视频免费看| 成年女人毛片免费观看观看9 | 国产精品免费大片| 国产xxxxx性猛交| 日韩中文字幕视频在线看片| 乱人伦中国视频| 亚洲欧洲国产日韩| 18禁观看日本| 观看av在线不卡| 老汉色∧v一级毛片| 亚洲美女黄色视频免费看| 激情视频va一区二区三区| bbb黄色大片| 国产又爽黄色视频| 天堂8中文在线网| 黑丝袜美女国产一区| 黄色视频在线播放观看不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美av亚洲av综合av国产av | 国产在线免费精品| 色婷婷久久久亚洲欧美| 香蕉国产在线看| 极品少妇高潮喷水抽搐| 男女国产视频网站| 久久久国产精品麻豆| 国产老妇伦熟女老妇高清| 只有这里有精品99| av在线app专区| 男女高潮啪啪啪动态图| 久久久国产精品麻豆| 国产片特级美女逼逼视频| 街头女战士在线观看网站| 伊人久久国产一区二区| 国产精品免费大片| 岛国毛片在线播放| 在线 av 中文字幕| 制服诱惑二区| 在线观看一区二区三区激情| 黄片播放在线免费| 国产av码专区亚洲av| 国产黄色视频一区二区在线观看| 免费在线观看黄色视频的| 欧美在线黄色| 国产99久久九九免费精品| 天堂中文最新版在线下载| 久久精品久久久久久久性| 桃花免费在线播放| 国产伦理片在线播放av一区| 人妻人人澡人人爽人人| 汤姆久久久久久久影院中文字幕| 十八禁高潮呻吟视频| 男女国产视频网站| 欧美最新免费一区二区三区| 少妇 在线观看| 午夜老司机福利片| 另类亚洲欧美激情| 热re99久久国产66热| 电影成人av| 国产爽快片一区二区三区| 又粗又硬又长又爽又黄的视频| 99re6热这里在线精品视频| 一区二区三区乱码不卡18| 亚洲人成电影观看| 91成人精品电影| 久久热在线av| av.在线天堂| 精品午夜福利在线看| 尾随美女入室| 一本大道久久a久久精品| 日韩大码丰满熟妇| 丁香六月天网| 黄片播放在线免费| 国产伦人伦偷精品视频| 国产精品一区二区精品视频观看| 亚洲av成人精品一二三区| 考比视频在线观看| av国产久精品久网站免费入址| 麻豆av在线久日| 高清在线视频一区二区三区| 国产精品欧美亚洲77777| 黄色视频不卡| 777久久人妻少妇嫩草av网站| 视频区图区小说| av视频免费观看在线观看| 精品卡一卡二卡四卡免费| 国产黄色免费在线视频| 99久久人妻综合| av天堂久久9| 国产有黄有色有爽视频| 亚洲av综合色区一区| 亚洲av福利一区| 九草在线视频观看| 国产黄色视频一区二区在线观看| 日韩制服骚丝袜av| 免费在线观看完整版高清| 超碰成人久久| 日本色播在线视频| avwww免费| 中文字幕高清在线视频| 这个男人来自地球电影免费观看 | 丰满乱子伦码专区| 国产男女超爽视频在线观看| 韩国高清视频一区二区三区| 精品久久蜜臀av无| 国产成人精品福利久久| 一级a爱视频在线免费观看| 性少妇av在线| 亚洲,欧美,日韩| 国产老妇伦熟女老妇高清| 久久久久精品久久久久真实原创| 国产午夜精品一二区理论片| 一区二区三区乱码不卡18| 制服诱惑二区| 在线观看免费视频网站a站| 大片免费播放器 马上看| 亚洲人成电影观看| 亚洲第一av免费看| 啦啦啦视频在线资源免费观看| 日韩av在线免费看完整版不卡| 91成人精品电影| 久久精品久久精品一区二区三区| 天天躁夜夜躁狠狠久久av| videos熟女内射| 久久青草综合色| 一级,二级,三级黄色视频| 国产亚洲av片在线观看秒播厂| 天堂8中文在线网| 欧美精品人与动牲交sv欧美| 日日啪夜夜爽| 青青草视频在线视频观看| 纵有疾风起免费观看全集完整版| av线在线观看网站| 亚洲情色 制服丝袜| 国产女主播在线喷水免费视频网站| 亚洲av成人不卡在线观看播放网 | 一边亲一边摸免费视频| 亚洲av电影在线进入| 青草久久国产| 国产免费又黄又爽又色| 卡戴珊不雅视频在线播放| 午夜日本视频在线| 国产精品麻豆人妻色哟哟久久| 校园人妻丝袜中文字幕| 亚洲国产欧美一区二区综合| 91成人精品电影| 中文字幕av电影在线播放| 免费高清在线观看视频在线观看| 国产色婷婷99| 99九九在线精品视频|