• <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)計類課程教學中的應用
    ——以多重共線性內容為例
    不完全多重共線性定義存在的問題及其修正建議
    對數簡史
    向量垂直在解析幾何中的應用
    亚洲国产色片| 欧美日韩精品成人综合77777| 久久鲁丝午夜福利片| 国产精品久久久久久久久免| 精品人妻熟女av久视频| 给我免费播放毛片高清在线观看| 国产成人91sexporn| 精品熟女少妇av免费看| 亚洲av五月六月丁香网| 啦啦啦韩国在线观看视频| 91精品国产九色| 国产成人aa在线观看| 给我免费播放毛片高清在线观看| or卡值多少钱| 免费不卡的大黄色大毛片视频在线观看 | 男人的好看免费观看在线视频| 嫩草影院入口| 成年免费大片在线观看| 国内久久婷婷六月综合欲色啪| 内地一区二区视频在线| 中文字幕av在线有码专区| 老熟妇乱子伦视频在线观看| 日韩欧美 国产精品| 深夜a级毛片| 精品乱码久久久久久99久播| 亚洲av一区综合| 亚洲va在线va天堂va国产| 51国产日韩欧美| 99久国产av精品| 淫妇啪啪啪对白视频| 中文资源天堂在线| 欧美日韩综合久久久久久| 在线观看av片永久免费下载| 日韩中字成人| 成人特级av手机在线观看| 精品久久国产蜜桃| 一级毛片我不卡| 女人被狂操c到高潮| 久久久精品大字幕| 波野结衣二区三区在线| 日本与韩国留学比较| 日韩欧美三级三区| 毛片一级片免费看久久久久| 麻豆成人午夜福利视频| 国产乱人偷精品视频| 欧美日本视频| 亚洲欧美精品自产自拍| 精华霜和精华液先用哪个| 波多野结衣高清作品| 97超视频在线观看视频| 国产精品久久久久久av不卡| 桃色一区二区三区在线观看| 国产大屁股一区二区在线视频| 亚洲在线自拍视频| 精品不卡国产一区二区三区| 精品少妇黑人巨大在线播放 | 欧美最黄视频在线播放免费| 99国产极品粉嫩在线观看| a级毛片a级免费在线| 亚洲成人久久爱视频| 99热精品在线国产| 男女啪啪激烈高潮av片| 久久久久久久午夜电影| 午夜免费激情av| 午夜久久久久精精品| 日本成人三级电影网站| 国产高清激情床上av| 美女内射精品一级片tv| 天堂网av新在线| 身体一侧抽搐| 中文字幕熟女人妻在线| 精品日产1卡2卡| 亚洲欧美日韩东京热| 亚洲成人中文字幕在线播放| 日韩高清综合在线| av专区在线播放| 日韩av不卡免费在线播放| 久久久精品94久久精品| 亚洲av电影不卡..在线观看| 亚洲中文字幕日韩| 免费高清视频大片| 97超级碰碰碰精品色视频在线观看| 午夜精品一区二区三区免费看| 日本成人三级电影网站| 九色成人免费人妻av| 亚洲成人久久爱视频| 欧美高清成人免费视频www| 99国产极品粉嫩在线观看| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 国产伦精品一区二区三区视频9| 日本a在线网址| 成熟少妇高潮喷水视频| 高清午夜精品一区二区三区 | 禁无遮挡网站| 精品一区二区三区av网在线观看| 国产精品久久久久久久电影| 91久久精品国产一区二区三区| 久99久视频精品免费| 看片在线看免费视频| 亚洲精品一区av在线观看| 国产精品女同一区二区软件| 九九爱精品视频在线观看| 欧美最新免费一区二区三区| 精品久久久久久久久av| 国产v大片淫在线免费观看| 国产精品一区二区三区四区免费观看 | 麻豆国产97在线/欧美| 国产成人a∨麻豆精品| 国产欧美日韩一区二区精品| 亚洲成人久久爱视频| 免费黄网站久久成人精品| 国产精品野战在线观看| 波野结衣二区三区在线| 亚洲真实伦在线观看| 亚洲熟妇中文字幕五十中出| av视频在线观看入口| 乱系列少妇在线播放| 在线播放国产精品三级| 在线播放无遮挡| 亚洲中文字幕日韩| 国产不卡一卡二| 国产一区二区三区在线臀色熟女| 老司机影院成人| 亚洲人成网站在线观看播放| 在线观看一区二区三区| 久久久国产成人免费| 午夜日韩欧美国产| 一a级毛片在线观看| 18禁裸乳无遮挡免费网站照片| 久久久a久久爽久久v久久| 国产精品一区二区性色av| 淫秽高清视频在线观看| 亚洲图色成人| 中文字幕av成人在线电影| 国产精品女同一区二区软件| 内地一区二区视频在线| 91午夜精品亚洲一区二区三区| 噜噜噜噜噜久久久久久91| 99久国产av精品| 欧美区成人在线视频| 免费高清视频大片| 欧美成人精品欧美一级黄| 婷婷精品国产亚洲av| 黑人高潮一二区| 热99re8久久精品国产| 熟妇人妻久久中文字幕3abv| 精品午夜福利在线看| 少妇人妻精品综合一区二区 | 欧美最新免费一区二区三区| 日韩强制内射视频| 中文在线观看免费www的网站| 久久精品91蜜桃| 国产一区二区在线av高清观看| 国产中年淑女户外野战色| 日韩精品中文字幕看吧| 国产一区二区亚洲精品在线观看| 色视频www国产| 国产精品国产高清国产av| 国产精品免费一区二区三区在线| 国产精品久久视频播放| 91在线观看av| 色吧在线观看| 女人被狂操c到高潮| 成人亚洲欧美一区二区av| 国产亚洲91精品色在线| 欧美高清性xxxxhd video| 日本黄大片高清| 久久久精品欧美日韩精品| 高清毛片免费看| 夜夜看夜夜爽夜夜摸| 久久久久久大精品| 老师上课跳d突然被开到最大视频| 色哟哟哟哟哟哟| 嫩草影院精品99| 精品99又大又爽又粗少妇毛片| 在线观看免费视频日本深夜| 久久久精品欧美日韩精品| 内地一区二区视频在线| 天天躁夜夜躁狠狠久久av| 成人午夜高清在线视频| 大又大粗又爽又黄少妇毛片口| 麻豆久久精品国产亚洲av| 国产精品女同一区二区软件| 国产探花在线观看一区二区| 久久久成人免费电影| 国产真实乱freesex| 国产精品亚洲美女久久久| 一进一出好大好爽视频| 女人被狂操c到高潮| 夜夜看夜夜爽夜夜摸| 中国国产av一级| 欧美国产日韩亚洲一区| 22中文网久久字幕| 欧美成人一区二区免费高清观看| 最近最新中文字幕大全电影3| 精品一区二区三区人妻视频| 国产成人a∨麻豆精品| 免费一级毛片在线播放高清视频| 少妇被粗大猛烈的视频| 能在线免费观看的黄片| 黄色日韩在线| 日本爱情动作片www.在线观看 | 亚洲熟妇熟女久久| 亚洲美女搞黄在线观看 | 欧美丝袜亚洲另类| 亚洲专区国产一区二区| 舔av片在线| 国产大屁股一区二区在线视频| 日本免费一区二区三区高清不卡| 亚洲av一区综合| 国产高清三级在线| 男女下面进入的视频免费午夜| 国产亚洲av嫩草精品影院| 97超碰精品成人国产| 久久精品国产自在天天线| 91在线精品国自产拍蜜月| 女的被弄到高潮叫床怎么办| 不卡视频在线观看欧美| 精华霜和精华液先用哪个| 岛国在线免费视频观看| 日本在线视频免费播放| 精品不卡国产一区二区三区| 国产蜜桃级精品一区二区三区| 国产男人的电影天堂91| 日本一二三区视频观看| 国产色爽女视频免费观看| 夜夜看夜夜爽夜夜摸| 中文在线观看免费www的网站| 搞女人的毛片| 91av网一区二区| 欧美在线一区亚洲| 亚洲不卡免费看| 国产探花极品一区二区| 一区福利在线观看| 午夜福利在线观看吧| 国产高清三级在线| 99视频精品全部免费 在线| 国产亚洲91精品色在线| 成熟少妇高潮喷水视频| 欧美日本亚洲视频在线播放| 禁无遮挡网站| 在线a可以看的网站| 久久人人爽人人爽人人片va| 国产精品永久免费网站| 看非洲黑人一级黄片| 99在线视频只有这里精品首页| 伊人久久精品亚洲午夜| 狂野欧美激情性xxxx在线观看| 欧美xxxx性猛交bbbb| 成人欧美大片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产高清在线一区二区三| 一a级毛片在线观看| 嫩草影院新地址| 欧美丝袜亚洲另类| 日韩成人伦理影院| 亚洲欧美日韩无卡精品| 国产黄色小视频在线观看| 精品日产1卡2卡| 波野结衣二区三区在线| 亚洲av美国av| av女优亚洲男人天堂| 久久久久久久久中文| 国产三级在线视频| av在线播放精品| 欧美国产日韩亚洲一区| 三级国产精品欧美在线观看| 国产一区二区三区在线臀色熟女| 少妇人妻精品综合一区二区 | 久久午夜福利片| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品一卡2卡三卡4卡5卡| www日本黄色视频网| 男女视频在线观看网站免费| 精品人妻一区二区三区麻豆 | 国产一区亚洲一区在线观看| 日韩 亚洲 欧美在线| 伊人久久精品亚洲午夜| 99久久精品热视频| 亚洲最大成人手机在线| 小蜜桃在线观看免费完整版高清| 村上凉子中文字幕在线| 69人妻影院| 亚洲第一电影网av| 我的老师免费观看完整版| 久久久精品欧美日韩精品| 亚洲人成网站在线播| 免费观看在线日韩| 亚洲国产欧美人成| 成年免费大片在线观看| 成年版毛片免费区| 国产精品一区二区三区四区久久| 国产精品久久电影中文字幕| 亚洲成人中文字幕在线播放| 啦啦啦啦在线视频资源| 日本五十路高清| 久久久国产成人免费| 九九在线视频观看精品| 18禁在线无遮挡免费观看视频 | 美女内射精品一级片tv| 久久鲁丝午夜福利片| 有码 亚洲区| 丝袜喷水一区| 日日撸夜夜添| 国产在视频线在精品| 国产不卡一卡二| 日本黄色片子视频| 日产精品乱码卡一卡2卡三| 简卡轻食公司| 99久久久亚洲精品蜜臀av| 18禁裸乳无遮挡免费网站照片| 国产探花极品一区二区| 一夜夜www| 成人精品一区二区免费| 亚洲无线在线观看| 成人精品一区二区免费| 亚洲无线在线观看| 国产精品人妻久久久久久| 麻豆国产97在线/欧美| 成人美女网站在线观看视频| 亚洲18禁久久av| 又爽又黄a免费视频| 日韩精品有码人妻一区| 免费黄网站久久成人精品| 一级a爱片免费观看的视频| 夜夜夜夜夜久久久久| 日韩av不卡免费在线播放| 精品一区二区三区视频在线| 一进一出抽搐gif免费好疼| 男插女下体视频免费在线播放| 国内精品一区二区在线观看| 亚洲av熟女| 欧美激情国产日韩精品一区| 嫩草影院新地址| 特级一级黄色大片| 一夜夜www| 国产伦一二天堂av在线观看| 日本免费a在线| 亚洲成av人片在线播放无| 亚洲精品粉嫩美女一区| 成人精品一区二区免费| 久久久精品欧美日韩精品| 久久精品国产亚洲av天美| 1000部很黄的大片| 在线免费十八禁| 亚洲经典国产精华液单| 亚洲最大成人av| 国内精品久久久久精免费| 黑人高潮一二区| 男女那种视频在线观看| 深夜精品福利| 国产精品电影一区二区三区| 欧美性猛交黑人性爽| 真人做人爱边吃奶动态| 久久久久久大精品| av在线亚洲专区| 舔av片在线| 91在线观看av| 男女之事视频高清在线观看| 色5月婷婷丁香| 99九九线精品视频在线观看视频| 97超碰精品成人国产| 亚洲精品国产成人久久av| 国产在线男女| 亚洲最大成人av| 欧美高清成人免费视频www| 九九热线精品视视频播放| 99久久精品一区二区三区| 亚洲内射少妇av| 一级毛片aaaaaa免费看小| 国产成人福利小说| 波多野结衣巨乳人妻| 亚洲精品粉嫩美女一区| 中文字幕av在线有码专区| 九九热线精品视视频播放| 亚洲中文日韩欧美视频| 听说在线观看完整版免费高清| 一级毛片aaaaaa免费看小| 有码 亚洲区| 真人做人爱边吃奶动态| 俄罗斯特黄特色一大片| 国产伦在线观看视频一区| 国产精品女同一区二区软件| 久久鲁丝午夜福利片| 欧美激情在线99| 91精品国产九色| 中文亚洲av片在线观看爽| 一本久久中文字幕| 国产成人精品久久久久久| 一本一本综合久久| 欧美三级亚洲精品| 俺也久久电影网| 夜夜看夜夜爽夜夜摸| 亚洲最大成人手机在线| 久久热精品热| 久久人人爽人人片av| 人人妻人人澡欧美一区二区| 国产v大片淫在线免费观看| 在线观看一区二区三区| 国产大屁股一区二区在线视频| 男女做爰动态图高潮gif福利片| 成人亚洲精品av一区二区| 成人国产麻豆网| 99久久成人亚洲精品观看| 亚洲av.av天堂| 久久人妻av系列| 国产精品国产三级国产av玫瑰| 欧美国产日韩亚洲一区| 国产成年人精品一区二区| 青春草视频在线免费观看| 黑人高潮一二区| 成人美女网站在线观看视频| a级毛片免费高清观看在线播放| 亚洲丝袜综合中文字幕| 99九九线精品视频在线观看视频| 欧美+日韩+精品| 亚洲欧美日韩卡通动漫| 国产亚洲精品久久久久久毛片| 六月丁香七月| 国产午夜福利久久久久久| av女优亚洲男人天堂| 嫩草影视91久久| 国产美女午夜福利| 成人av一区二区三区在线看| 高清毛片免费观看视频网站| 国产黄片美女视频| 亚洲最大成人av| 日本一本二区三区精品| 午夜影院日韩av| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲av涩爱 | 色吧在线观看| 晚上一个人看的免费电影| 亚洲激情五月婷婷啪啪| 最后的刺客免费高清国语| 久久久久国产网址| 日韩强制内射视频| 联通29元200g的流量卡| 国产精品,欧美在线| 蜜桃久久精品国产亚洲av| 国产精品无大码| 美女 人体艺术 gogo| 国产亚洲精品久久久久久毛片| 亚洲av一区综合| 少妇的逼水好多| АⅤ资源中文在线天堂| 中文在线观看免费www的网站| 国产午夜精品久久久久久一区二区三区 | 久久午夜福利片| 少妇高潮的动态图| 美女内射精品一级片tv| 12—13女人毛片做爰片一| 秋霞在线观看毛片| 精华霜和精华液先用哪个| 久久精品国产鲁丝片午夜精品| 国产精品一及| 久久久久久久久中文| 天美传媒精品一区二区| 麻豆久久精品国产亚洲av| 最后的刺客免费高清国语| 12—13女人毛片做爰片一| 麻豆精品久久久久久蜜桃| 在线观看一区二区三区| 午夜精品一区二区三区免费看| 国产成人福利小说| 伊人久久精品亚洲午夜| 日本一二三区视频观看| 99国产极品粉嫩在线观看| 人妻久久中文字幕网| 免费av毛片视频| 精品国内亚洲2022精品成人| 欧美日韩乱码在线| 精品不卡国产一区二区三区| 国产麻豆成人av免费视频| 美女免费视频网站| 亚洲人成网站高清观看| 乱人视频在线观看| 亚洲色图av天堂| 久久九九热精品免费| 男女啪啪激烈高潮av片| 精品99又大又爽又粗少妇毛片| 久久久久久伊人网av| 国产精品野战在线观看| 亚洲欧美日韩卡通动漫| 三级国产精品欧美在线观看| 久久久国产成人精品二区| 色综合亚洲欧美另类图片| 成人一区二区视频在线观看| 亚洲七黄色美女视频| 你懂的网址亚洲精品在线观看 | 欧美成人一区二区免费高清观看| 国内精品美女久久久久久| 国产日本99.免费观看| 国产不卡一卡二| 精品一区二区三区av网在线观看| 综合色av麻豆| 国产黄片美女视频| 熟女人妻精品中文字幕| 久久精品久久久久久噜噜老黄 | 午夜激情欧美在线| 在线观看一区二区三区| 亚洲在线观看片| 免费大片18禁| 高清午夜精品一区二区三区 | 精品人妻一区二区三区麻豆 | 天堂影院成人在线观看| 成人特级av手机在线观看| 国产探花极品一区二区| 日本欧美国产在线视频| 午夜影院日韩av| 精品人妻一区二区三区麻豆 | 久久久久精品国产欧美久久久| 久久精品国产清高在天天线| 久久久欧美国产精品| 一个人观看的视频www高清免费观看| 国产老妇女一区| 国产美女午夜福利| 男女做爰动态图高潮gif福利片| 天堂网av新在线| 午夜亚洲福利在线播放| 精品一区二区三区视频在线观看免费| 国产黄色视频一区二区在线观看 | 成人毛片a级毛片在线播放| 91精品国产九色| 天堂动漫精品| 精品少妇黑人巨大在线播放 | 69av精品久久久久久| 国产一级毛片七仙女欲春2| 性欧美人与动物交配| 亚洲一级一片aⅴ在线观看| 亚洲无线观看免费| 免费观看的影片在线观看| 免费看a级黄色片| 久久人妻av系列| 免费看美女性在线毛片视频| 免费观看在线日韩| 一级av片app| 亚洲熟妇中文字幕五十中出| 看免费成人av毛片| 简卡轻食公司| 草草在线视频免费看| 99久久精品一区二区三区| 国产成人a∨麻豆精品| 日韩av不卡免费在线播放| 亚洲图色成人| 中国国产av一级| 夜夜看夜夜爽夜夜摸| 干丝袜人妻中文字幕| 亚洲18禁久久av| 久久鲁丝午夜福利片| 亚洲欧美中文字幕日韩二区| 精品国产三级普通话版| 久久久久久久久中文| 在线观看午夜福利视频| 国产av在哪里看| 又黄又爽又刺激的免费视频.| a级毛片免费高清观看在线播放| 国内精品一区二区在线观看| 国产精品美女特级片免费视频播放器| 日韩中字成人| 赤兔流量卡办理| 亚洲婷婷狠狠爱综合网| 亚洲真实伦在线观看| 国产在线精品亚洲第一网站| 一级毛片久久久久久久久女| 麻豆久久精品国产亚洲av| 国产91av在线免费观看| 九九热线精品视视频播放| 亚洲人成网站在线播放欧美日韩| 热99re8久久精品国产| 日韩欧美一区二区三区在线观看| 午夜视频国产福利| 午夜福利在线观看吧| 国产男靠女视频免费网站| 精品一区二区三区av网在线观看| 欧美日韩乱码在线| 悠悠久久av| 亚洲国产精品久久男人天堂| 最近视频中文字幕2019在线8| 一个人看的www免费观看视频| 精品人妻熟女av久视频| 国产精品久久久久久精品电影| 小蜜桃在线观看免费完整版高清| 久久韩国三级中文字幕| 亚洲精品在线观看二区| 国产单亲对白刺激| 免费高清视频大片| 国产精品1区2区在线观看.| 97超级碰碰碰精品色视频在线观看| 一进一出好大好爽视频| 国产精品99久久久久久久久| av女优亚洲男人天堂| 日韩欧美一区二区三区在线观看| 在线国产一区二区在线| 亚洲精品国产av成人精品 | av黄色大香蕉| 综合色丁香网| 在线观看av片永久免费下载| 老女人水多毛片| 亚洲欧美日韩卡通动漫| 日日啪夜夜撸| 欧美一级a爱片免费观看看| 欧美激情在线99| 蜜臀久久99精品久久宅男| 国内精品久久久久精免费| 嫩草影院入口| 色在线成人网| 免费观看精品视频网站| 91在线精品国自产拍蜜月| 波多野结衣高清无吗| 男女那种视频在线观看| 免费搜索国产男女视频| 日产精品乱码卡一卡2卡三| 中文字幕久久专区|