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

    四種坐標(biāo)函數(shù)對流式細(xì)胞術(shù)數(shù)據(jù)可視化的影響

    2015-04-18 02:38:26杜慶華李慶山許艷麗
    海南醫(yī)學(xué) 2015年15期
    關(guān)鍵詞:對數(shù)線性運(yùn)算

    杜慶華,李慶山,許艷麗

    (廣州醫(yī)科大學(xué)附屬廣州市第一人民醫(yī)院血液內(nèi)科,廣東 廣州 510180)

    四種坐標(biāo)函數(shù)對流式細(xì)胞術(shù)數(shù)據(jù)可視化的影響

    杜慶華,李慶山,許艷麗

    (廣州醫(yī)科大學(xué)附屬廣州市第一人民醫(yī)院血液內(nèi)科,廣東 廣州 510180)

    目的 研究流式細(xì)胞術(shù)數(shù)據(jù)可視化中四種常見坐標(biāo)函數(shù)的特點(diǎn),探討這四種坐標(biāo)函數(shù)在數(shù)據(jù)可視化中的應(yīng)用。方法對比分析線性、對數(shù)、HyperLog以及Logical四個坐標(biāo)函數(shù)曲線的特點(diǎn),使用不同坐標(biāo)函數(shù)對同一個數(shù)據(jù)使用散點(diǎn)圖進(jìn)行顯示,對比圖形顯示的差異。結(jié)果線性函數(shù)對數(shù)據(jù)呈等比例顯示,但動態(tài)范圍不足。對數(shù)函數(shù)能拉伸1附近的數(shù)值,壓縮顯示數(shù)值大的數(shù)值,但不能顯示小于1的值。HyperLog與Logical函數(shù)數(shù)值小的時候近似線性顯示,數(shù)值大時近似對數(shù)顯示。結(jié)論不同的坐標(biāo)函數(shù)有不同的特點(diǎn),在進(jìn)行流式細(xì)胞術(shù)數(shù)據(jù)分析時,要根據(jù)需要選擇相應(yīng)的坐標(biāo)函數(shù)進(jìn)行顯示。

    流式細(xì)胞術(shù);坐標(biāo)函數(shù);數(shù)據(jù)可視化

    數(shù)據(jù)可視化是流式細(xì)胞術(shù)數(shù)據(jù)分析中的一個重要環(huán)節(jié),數(shù)據(jù)通過圖形顯示,然后才能通過設(shè)門操作進(jìn)行進(jìn)一步分析。流式數(shù)據(jù)的顯示樣式有直方圖、散點(diǎn)圖、等高圖或者密度圖等。每個細(xì)胞的熒光強(qiáng)度數(shù)據(jù)要在圖中顯示均需要通過坐標(biāo)函數(shù)計算其位置,因此坐標(biāo)函數(shù)的選擇會直接影響圖形的形狀。流式細(xì)胞術(shù)中常用的坐標(biāo)有線性坐標(biāo)、對數(shù)坐標(biāo)以及其他對數(shù)衍生的坐標(biāo)。由于細(xì)胞表面及內(nèi)部分子表達(dá)的差異很大,這些數(shù)據(jù)很難通過線性坐標(biāo)顯示[1-2],因此一直以來使用流式細(xì)胞術(shù)進(jìn)行免疫表型分析基本都使用對數(shù)坐標(biāo)以顯示數(shù)據(jù)。但因?yàn)閷?shù)無法恰當(dāng)?shù)仫@示小于1的數(shù)據(jù),故Bagwell及Parks建議使用對數(shù)衍生函數(shù)進(jìn)行坐標(biāo)轉(zhuǎn)換[3-4]。而這些坐標(biāo)函數(shù)間特性的比較暫時未見詳細(xì)報道。本文將深入探討不同坐標(biāo)函數(shù)的特點(diǎn),為坐標(biāo)函數(shù)的選擇提供指引。

    1 材料與方法

    1.1 儀器及軟件 流式細(xì)胞儀為FACS Canto (美國BD公司),數(shù)據(jù)為常規(guī)淋巴細(xì)胞亞群臨床檢測的流式數(shù)據(jù),軟件使用本實(shí)驗(yàn)室自行研發(fā)的流式數(shù)據(jù)分析軟件CFCS(軟件著作權(quán)登記號:2010SR064983)。

    1.2 方法

    1.2.1 函數(shù)曲線的繪制 對數(shù)函數(shù)及其反函數(shù)為:

    其中r為數(shù)據(jù)的分辨率,d為動態(tài)范圍的數(shù)量級。HypherLog為隱函數(shù)[3],故僅能寫出其反函數(shù):

    因此在進(jìn)行HypherLog函數(shù)運(yùn)算時,必須通過其反函數(shù)進(jìn)行求根運(yùn)算。函數(shù)中b為相關(guān)系數(shù),是影響線性范圍與對數(shù)范圍的參數(shù)。當(dāng)b=0時曲線最接近對數(shù)[3]。

    Logical與HypherLog同樣是隱函數(shù)[4-5],其反函數(shù)為:

    其中M為圖形顯示寬寬的數(shù)量級,而W=2p log (p)/(p+1),為線性范圍的寬度,其值可以根據(jù)以下公式求出:

    其中n為負(fù)值范圍參考點(diǎn)的值。

    1.2.2 坐標(biāo)函數(shù)曲線的繪制 使用MATLAB7.0繪制線性、對數(shù)函數(shù)、HyperLog函數(shù)以及Logical函數(shù)的函數(shù)曲線。

    1.2.3 不同坐標(biāo)函數(shù)對圖形的影響 同一個數(shù)據(jù)在CFCS上分別使用線性坐標(biāo)、對數(shù)坐標(biāo)、HyperLog以及Logical顯示,并對比四者圖形的區(qū)別。

    2 結(jié)果

    2.1 不同的坐標(biāo)函數(shù)曲線的特點(diǎn) 從函數(shù)曲線發(fā)現(xiàn)線性坐標(biāo)能顯示整個實(shí)數(shù)域的數(shù)據(jù),且具有對稱性。因其等比例的特性,當(dāng)數(shù)據(jù)動態(tài)范圍較大時會使數(shù)值小的數(shù)據(jù)被壓縮,造成顯示效果不佳。對數(shù)坐標(biāo)能保證數(shù)值小數(shù)據(jù)的顯示空間,但其缺點(diǎn)是僅僅能顯示≥1的數(shù)據(jù),對于<1的數(shù)值全按0來處理。HyperLog與Logical能顯示實(shí)數(shù)域的數(shù)據(jù),且二者在數(shù)值較小的時候函數(shù)曲線接近線性,在數(shù)值較大的時候接近對數(shù)曲線(見圖1)。

    圖1 線性、對數(shù)、HyperLog以及Logical的函數(shù)曲線

    2.2 不同的坐標(biāo)函數(shù)對流式數(shù)據(jù)顯示的影響 通過使用不同坐標(biāo)顯示相同的免疫熒光數(shù)據(jù),發(fā)現(xiàn)線性坐標(biāo)顯示時細(xì)胞群體往往聚集在數(shù)軸的一端,難以對群體進(jìn)行區(qū)分。對數(shù)坐標(biāo)顯示免疫熒光實(shí)驗(yàn)的數(shù)據(jù)時,細(xì)胞群體分布比較清晰,但是進(jìn)行多色熒光分析進(jìn)行補(bǔ)償后,有8%的細(xì)胞堆疊在基線上,造成視覺誤差,容易對分析結(jié)果造成影響。Logical與HypherLog兩者顯示效果無顯著差異,均能很好地控制堆疊在基線上細(xì)胞的數(shù)量,細(xì)胞群體分布清晰,補(bǔ)償后群體離散程度較對數(shù)坐標(biāo)明顯(見圖2)。

    圖2 使用4種坐標(biāo)函數(shù)對同一個數(shù)據(jù)生成散點(diǎn)圖的比較

    3 討論

    線性坐標(biāo)就是把數(shù)據(jù)直接按比例顯示在圖形上,其優(yōu)點(diǎn)是簡單直接,能顯示實(shí)數(shù)域的所有數(shù)據(jù),一般用于DNA含量的的測定,如細(xì)胞周期分析。但細(xì)胞表面分子表達(dá)情況差異相當(dāng)大,需要極寬的動態(tài)范圍才能顯示這類型的數(shù)據(jù),而且這類型的數(shù)據(jù)低熒光強(qiáng)度群體離散程度小并近似正態(tài)分布,高熒光強(qiáng)度群體離散程度大接近對數(shù)正態(tài)分布。線性坐標(biāo)顯示動態(tài)范圍寬的數(shù)據(jù)時,會壓縮數(shù)值小數(shù)據(jù)的顯示范圍,造成細(xì)胞群體往往聚集在數(shù)軸的一端。而對數(shù)坐標(biāo)能拉伸數(shù)值小數(shù)據(jù)的顯示范圍,壓縮數(shù)值大的數(shù)據(jù),擁有較寬的動態(tài)范圍,因此對數(shù)坐標(biāo)比線性坐標(biāo)更適合顯示免疫熒光的數(shù)據(jù)。但是對數(shù)坐標(biāo)的缺陷是無法對<1的值進(jìn)行變換,因此一般對于<1的值轉(zhuǎn)換為0。

    在實(shí)際應(yīng)用中為了消除通道間的熒光滲漏,往往需要對數(shù)據(jù)進(jìn)行補(bǔ)償運(yùn)算,該運(yùn)算是多個通道間的減法運(yùn)算。運(yùn)算過程中會把某一通道的檢測誤差引入到其他通道中,因此會造成數(shù)據(jù)離散增大,如果細(xì)胞群體被影響通道的熒光強(qiáng)度較低時,數(shù)據(jù)的離散將造成<1甚者負(fù)值數(shù)據(jù)的產(chǎn)生,這些數(shù)據(jù)在對數(shù)坐標(biāo)中堆疊在基線下,使該處細(xì)胞密度異常增高,造成假群體的產(chǎn)生[4,6]。從圖2中我們可看到8%的細(xì)胞堆在x坐標(biāo)上。如果設(shè)門分析時忽略了這些細(xì)胞,將對結(jié)果造成一定的影響。為了克服對數(shù)坐標(biāo)的這些缺陷,Parks提出理想的坐標(biāo)函數(shù)應(yīng)具備以下特點(diǎn):(1)函數(shù)可對不同的數(shù)據(jù)進(jìn)行顯示的優(yōu)化;(2)為了提供足夠?qū)挼膭討B(tài)范圍,更好地顯示高熒光強(qiáng)度的對數(shù)正態(tài)分布數(shù)據(jù),函數(shù)應(yīng)在隨數(shù)值增大逼近對數(shù)曲線;(3)函數(shù)應(yīng)在接近0的區(qū)域?qū)ΨQ并近似線性分布,這樣更適合顯示低熒光強(qiáng)度正態(tài)分布的數(shù)據(jù);(4)線性區(qū)域應(yīng)盡可能平滑地向?qū)?shù)區(qū)域過渡,以避免扭曲地顯示數(shù)據(jù);(5)線性化強(qiáng)度應(yīng)隨線性范圍增大而增大[4]。

    Bagwell及Parks分別提出了HypherLog與Logical函數(shù)。這兩個函數(shù)的特點(diǎn)是能顯示負(fù)值的數(shù)據(jù),且小數(shù)值區(qū)域接近線性分布,高數(shù)值區(qū)接近對數(shù)分布,且線性區(qū)與對數(shù)區(qū)能平滑過渡。在實(shí)際顯示中通過HypherLog及Logical兩個坐標(biāo)與對數(shù)坐標(biāo)比較,發(fā)現(xiàn)對于熒光強(qiáng)度大的細(xì)胞群體差異不大。HypherLog及Logical兩個坐標(biāo)下,熒光強(qiáng)度低的群體比對數(shù)坐標(biāo)更趨于聚集,更重要的是這兩個坐標(biāo)下能顯示對數(shù)坐標(biāo)不能顯示的數(shù)值<1的值。此外,HypherLog及Logical能識別一些因過補(bǔ)償而造成細(xì)胞群體落在坐標(biāo)軸上的數(shù)據(jù),并將之顯示出來。不同數(shù)據(jù)<0的值往往是不一樣的,故HypherLog與Logical需要適當(dāng)調(diào)整參數(shù)使得數(shù)據(jù)能恰當(dāng)?shù)仫@示。HypherLog通過相關(guān)系數(shù)b來控制線性區(qū)的斜率,b越大線性區(qū)域越寬,b為0時曲線最接近對數(shù)曲線。但Bagwell發(fā)現(xiàn)b太小會使陰性區(qū)域的原來為一群的群體在視覺上分裂為兩個,因此使用時要注意b的取值。Logical使用負(fù)值大小作為參數(shù)對曲線的先行范圍進(jìn)行控制,通過公式5計算W的值,-W與W之間即為線性區(qū)域。

    HypherLog的二階導(dǎo)數(shù)不為0而Logical為0,因此在0附近Logical比HypherLog更接近線性,且更快靠近對數(shù)曲線。從圖2實(shí)際顯示中看,這差異并不會對數(shù)據(jù)顯示分析造成很大的影響。

    由于不同的參數(shù)會對顯示效果造成一定影響,因此HypherLog與Logical兩個函數(shù)的穩(wěn)定性不如對數(shù)。當(dāng)進(jìn)行表達(dá)模式分析(如白血病免疫學(xué)分型)的時候,細(xì)胞群體的位置及走向是極其重要的信息。當(dāng)數(shù)據(jù)中含有個別極端數(shù)值時,要將其恰當(dāng)顯示則需要調(diào)整HypherLog及Logical函數(shù)的參數(shù),細(xì)胞群體的形狀及位置可能因此而發(fā)生扭曲改變(特別是熒光強(qiáng)度低的群體),容易導(dǎo)致分析錯誤。

    對數(shù)函數(shù)能直接運(yùn)算得到結(jié)果,而HypherLog與Logical無法顯式寫出其函數(shù)形式而無法直接進(jìn)行運(yùn)算,因此必須對方程進(jìn)行求根運(yùn)算。在分析軟件中,一般使用牛頓迭代法進(jìn)行求根運(yùn)算,運(yùn)算結(jié)果的精度隨迭代次數(shù)增加而增加。HypherLog與Logical二者運(yùn)算效率遠(yuǎn)低于對數(shù)運(yùn)算。Logical運(yùn)算的復(fù)雜程度大于HypherLog,故其運(yùn)算效率最低。在進(jìn)行大數(shù)據(jù)量的分析時我們可以預(yù)先建立一個函數(shù)表,進(jìn)行運(yùn)算時可通過二分查找法進(jìn)行檢索,從而極大地提高運(yùn)算的效率。

    變異范圍小和正態(tài)分布的數(shù)據(jù)(如細(xì)胞周期分析)宜使用線性坐標(biāo)顯示。免疫熒光這類動態(tài)范圍大的數(shù)據(jù)不宜使用線性顯示。如果這類數(shù)據(jù)無需考慮群體形狀走勢等信息,僅僅進(jìn)行群體劃分,則直接使用HypherLog或Logical進(jìn)行顯示設(shè)門即可。如果需要進(jìn)行表達(dá)模式分析,建議先使用對數(shù)顯示以了解群體分布情況,再使用HypherLog或Logical進(jìn)行設(shè)門,以避免負(fù)值數(shù)據(jù)的丟失。總之,不同的坐標(biāo)函數(shù)有其不同的特性,對數(shù)據(jù)的顯示及分析有著重要的影響。因此我們必須對不同坐標(biāo)函數(shù)的特性有所了解,分析時使用不同的坐標(biāo)函數(shù)顯示數(shù)據(jù),才能對數(shù)據(jù)有更全面的認(rèn)識,以避免坐標(biāo)函數(shù)選擇不當(dāng)而造成的視覺誤差,進(jìn)而影響分析結(jié)果。

    [1]Muirhead KA,Schmitt TC,Muirhead AR.Muirhead,determination of linear fluorescence intensities from flow cytometric data accumulated with logarithmic amplifiers[J].Cytometry,1983,3(4):251-256.

    [2]Watson JV,Chambers SH,Smith PJ.A pragmatic approach to the analysis of DNA histograms with a definable G1 peak[J].Cytometry,1987,8(1):1-8.

    [3]Bagwell CB.Hyperlog-a flexible log-like transform for negative,zero,and positive valued data[J].CytometryA,2005,64(1):34-42.

    [4]Parks DR,Roederer M,Moore WA.A new"Logicle"display method avoids deceptive effects of logarithmic scaling for low signals and compensated data[J].CytometryA,2006,69(6):541-551.

    [5]Moore WA,Parks DR.Update for the logicle data scale including operational code implementations[J].Cytometry A,2012,81(4): 273-277.

    [6]Finak G,Perez JM,Weng A,et al.Optimizing transformations for automated,high throughput analysis of flow cytometry data[J]. BMC Bioinformatics,2010,11:546.

    Effects of four scaling functions on flow cytometry data visualization.

    DU Qing-hua,LI Qing-shan,XU Yan-li.
    Department of Hematology,Guangzhou First People's Hospital,Guangzhou Medical University,Guangzhou 510180, Guangdong,CHINA

    ObjectiveTo explore the features of four scaling functions,and investigate the selection of 4 scaling functions in data visualization for flow cytometry.MethodsWe compared the features of 4 function curves, and dot plots from a same data file in 4 different scaling functions.ResultsLinear function displays data in a manner of proportion,but it is insufficient in dynamic range.Logarithmic function can stretch displaying the data value near zero,and can compress displaying the large data value,but it can not display data value less than zero.HyperLog and Logical function can display small data value similar to linear distribution,and display large data value similar to near logarithm.ConclusionThere are different features in different scaling functions.We should choose correct scaling function to display in data analyzing of flow cytometry.

    Flow cytometry;Scaling function;Data visualization

    R446

    A

    1003—6350(2015)15—2259—03

    10.3969/j.issn.1003-6350.2015.15.0814

    2015-03-02)

    2013年廣東省自然科學(xué)基金(編號:S2013010016726)

    李慶山。E-mail:qingshanli@hotmail.com

    猜你喜歡
    對數(shù)線性運(yùn)算
    漸近線性Klein-Gordon-Maxwell系統(tǒng)正解的存在性
    重視運(yùn)算與推理,解決數(shù)列求和題
    含有對數(shù)非線性項Kirchhoff方程多解的存在性
    指數(shù)與對數(shù)
    線性回歸方程的求解與應(yīng)用
    指數(shù)與對數(shù)
    有趣的運(yùn)算
    二階線性微分方程的解法
    對數(shù)簡史
    “整式的乘法與因式分解”知識歸納
    国产男女超爽视频在线观看| 午夜精品国产一区二区电影| 老司机深夜福利视频在线观看 | 亚洲精品国产色婷婷电影| √禁漫天堂资源中文www| 亚洲精品国产色婷婷电影| 亚洲av电影在线观看一区二区三区| 搡老乐熟女国产| 亚洲色图综合在线观看| 一本久久精品| 侵犯人妻中文字幕一二三四区| 中文字幕制服av| 亚洲成人免费av在线播放| 久久毛片免费看一区二区三区| 99九九在线精品视频| 精品福利观看| 日韩欧美国产一区二区入口| 黑人巨大精品欧美一区二区蜜桃| 在线永久观看黄色视频| 美女脱内裤让男人舔精品视频| av网站免费在线观看视频| 香蕉国产在线看| e午夜精品久久久久久久| 亚洲精品国产一区二区精华液| 欧美大码av| av网站在线播放免费| 狂野欧美激情性xxxx| 久久ye,这里只有精品| 高潮久久久久久久久久久不卡| 美国免费a级毛片| 国产成人av激情在线播放| 国产成人精品久久二区二区免费| 99热全是精品| 久久国产精品男人的天堂亚洲| 亚洲久久久国产精品| 黄频高清免费视频| 国产欧美日韩一区二区精品| 黑人猛操日本美女一级片| 美女高潮到喷水免费观看| 性色av乱码一区二区三区2| 这个男人来自地球电影免费观看| 老司机午夜福利在线观看视频 | 久9热在线精品视频| 狠狠狠狠99中文字幕| 91国产中文字幕| 亚洲免费av在线视频| 欧美大码av| 天天躁日日躁夜夜躁夜夜| 黑人巨大精品欧美一区二区mp4| 韩国精品一区二区三区| 欧美黄色片欧美黄色片| 纵有疾风起免费观看全集完整版| 精品久久久久久电影网| 久久久久久久精品精品| 侵犯人妻中文字幕一二三四区| 亚洲少妇的诱惑av| 在线观看免费视频网站a站| 亚洲第一av免费看| 99精品欧美一区二区三区四区| 九色亚洲精品在线播放| 99精品久久久久人妻精品| av不卡在线播放| 日本wwww免费看| 满18在线观看网站| 免费女性裸体啪啪无遮挡网站| 成在线人永久免费视频| 国产精品自产拍在线观看55亚洲 | 成年人午夜在线观看视频| 中文字幕精品免费在线观看视频| 欧美 亚洲 国产 日韩一| 美女视频免费永久观看网站| 亚洲第一青青草原| 高清黄色对白视频在线免费看| 一进一出抽搐动态| 亚洲 欧美一区二区三区| 狂野欧美激情性xxxx| 免费一级毛片在线播放高清视频 | 别揉我奶头~嗯~啊~动态视频 | 欧美av亚洲av综合av国产av| 国产男女超爽视频在线观看| 各种免费的搞黄视频| 久久精品国产综合久久久| 亚洲专区中文字幕在线| 免费高清在线观看日韩| 亚洲成国产人片在线观看| 色播在线永久视频| 在线观看人妻少妇| 大型av网站在线播放| 两人在一起打扑克的视频| 啦啦啦中文免费视频观看日本| 国产精品99久久99久久久不卡| 少妇精品久久久久久久| 中文字幕人妻丝袜制服| 最近最新免费中文字幕在线| 法律面前人人平等表现在哪些方面 | 啪啪无遮挡十八禁网站| 少妇的丰满在线观看| 精品欧美一区二区三区在线| 91精品国产国语对白视频| 97精品久久久久久久久久精品| 亚洲成人免费av在线播放| 亚洲国产日韩一区二区| 后天国语完整版免费观看| 黑人猛操日本美女一级片| 欧美人与性动交α欧美软件| 在线观看人妻少妇| 最近中文字幕2019免费版| 亚洲精品自拍成人| 日本猛色少妇xxxxx猛交久久| 美国免费a级毛片| 中文字幕制服av| 成年人免费黄色播放视频| 人成视频在线观看免费观看| av网站在线播放免费| 国产成+人综合+亚洲专区| 亚洲专区中文字幕在线| www.av在线官网国产| 中文欧美无线码| 国产精品欧美亚洲77777| av电影中文网址| 亚洲精品国产av蜜桃| 欧美激情 高清一区二区三区| 纵有疾风起免费观看全集完整版| 日本撒尿小便嘘嘘汇集6| 亚洲精品av麻豆狂野| 国产精品国产三级国产专区5o| 欧美+亚洲+日韩+国产| 人妻一区二区av| 母亲3免费完整高清在线观看| 夜夜夜夜夜久久久久| 黑人猛操日本美女一级片| 亚洲av片天天在线观看| 亚洲精品av麻豆狂野| 久久午夜综合久久蜜桃| 伦理电影免费视频| 成人三级做爰电影| av网站在线播放免费| 日本vs欧美在线观看视频| 久久这里只有精品19| 亚洲成国产人片在线观看| 国产成人啪精品午夜网站| 伦理电影免费视频| xxxhd国产人妻xxx| 夫妻午夜视频| 精品一区二区三区av网在线观看 | 亚洲精品一区蜜桃| 大陆偷拍与自拍| 精品国产乱码久久久久久小说| 亚洲第一欧美日韩一区二区三区 | 亚洲精品国产色婷婷电影| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品国产av蜜桃| 啦啦啦 在线观看视频| 又大又爽又粗| 日本av手机在线免费观看| 欧美日韩av久久| 国产日韩欧美亚洲二区| 国产伦人伦偷精品视频| 午夜激情久久久久久久| 黄色片一级片一级黄色片| 亚洲情色 制服丝袜| 国产老妇伦熟女老妇高清| 黑人巨大精品欧美一区二区蜜桃| 久久人妻熟女aⅴ| cao死你这个sao货| 国产片内射在线| 9191精品国产免费久久| 12—13女人毛片做爰片一| 十八禁网站网址无遮挡| 国产国语露脸激情在线看| 亚洲av国产av综合av卡| 免费av中文字幕在线| 狂野欧美激情性bbbbbb| 男女高潮啪啪啪动态图| 99久久国产精品久久久| 国产欧美日韩综合在线一区二区| 91老司机精品| 三上悠亚av全集在线观看| 亚洲成av片中文字幕在线观看| 最近最新免费中文字幕在线| 人人澡人人妻人| 亚洲精品中文字幕在线视频| 日韩,欧美,国产一区二区三区| 制服诱惑二区| 亚洲avbb在线观看| 国产亚洲精品一区二区www | 久久人妻福利社区极品人妻图片| 操美女的视频在线观看| 免费一级毛片在线播放高清视频 | 俄罗斯特黄特色一大片| av天堂在线播放| 999久久久精品免费观看国产| 建设人人有责人人尽责人人享有的| 在线永久观看黄色视频| 国产色视频综合| 久久综合国产亚洲精品| 久热爱精品视频在线9| 热re99久久精品国产66热6| 午夜福利,免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 可以免费在线观看a视频的电影网站| 夜夜骑夜夜射夜夜干| 丝袜在线中文字幕| 1024香蕉在线观看| 99国产精品免费福利视频| 波多野结衣一区麻豆| 黄色怎么调成土黄色| 18禁观看日本| 欧美黑人欧美精品刺激| 美女主播在线视频| 亚洲欧美日韩另类电影网站| 国产精品久久久久久人妻精品电影 | 精品乱码久久久久久99久播| 男女下面插进去视频免费观看| 中文精品一卡2卡3卡4更新| 欧美国产精品va在线观看不卡| 精品亚洲成国产av| 99国产综合亚洲精品| 视频在线观看一区二区三区| 热99国产精品久久久久久7| 在线观看人妻少妇| 满18在线观看网站| svipshipincom国产片| 美女国产高潮福利片在线看| 亚洲,欧美精品.| 老司机影院毛片| av一本久久久久| 日韩 亚洲 欧美在线| 五月开心婷婷网| 欧美中文综合在线视频| 国产精品久久久久久精品古装| 精品国产乱码久久久久久男人| 欧美日韩亚洲国产一区二区在线观看 | 在线精品无人区一区二区三| 人妻久久中文字幕网| 免费日韩欧美在线观看| 国产精品一区二区在线观看99| 国产一区二区在线观看av| 精品国产乱子伦一区二区三区 | 中国美女看黄片| 精品国产国语对白av| 国产亚洲av片在线观看秒播厂| 香蕉丝袜av| 黄片播放在线免费| 人人妻,人人澡人人爽秒播| 国产91精品成人一区二区三区 | 亚洲av男天堂| 曰老女人黄片| 久久女婷五月综合色啪小说| 两人在一起打扑克的视频| 高清欧美精品videossex| 午夜福利一区二区在线看| 亚洲人成电影观看| 亚洲欧美精品综合一区二区三区| 国产又色又爽无遮挡免| 操美女的视频在线观看| 久久久久久久国产电影| 欧美日韩中文字幕国产精品一区二区三区 | 国产麻豆69| 夜夜骑夜夜射夜夜干| 亚洲,欧美精品.| 一进一出抽搐动态| 亚洲中文字幕日韩| 亚洲av电影在线进入| 久久久国产成人免费| 下体分泌物呈黄色| 十八禁网站免费在线| √禁漫天堂资源中文www| 丰满饥渴人妻一区二区三| 高清视频免费观看一区二区| 国产男人的电影天堂91| 久久久久久免费高清国产稀缺| 欧美黑人精品巨大| 夜夜夜夜夜久久久久| 亚洲午夜精品一区,二区,三区| 久久狼人影院| 亚洲天堂av无毛| 国产1区2区3区精品| a级片在线免费高清观看视频| 亚洲熟女精品中文字幕| 亚洲自偷自拍图片 自拍| 日韩,欧美,国产一区二区三区| 国产成人啪精品午夜网站| 久久久久久免费高清国产稀缺| 啪啪无遮挡十八禁网站| 亚洲精品国产av蜜桃| 亚洲精品国产av蜜桃| 欧美在线黄色| 国产区一区二久久| 亚洲午夜精品一区,二区,三区| 日韩视频一区二区在线观看| 色婷婷av一区二区三区视频| 狠狠精品人妻久久久久久综合| 免费观看人在逋| 国产亚洲午夜精品一区二区久久| 日韩免费高清中文字幕av| 亚洲精品美女久久av网站| 欧美日韩视频精品一区| 日韩中文字幕欧美一区二区| 亚洲av日韩精品久久久久久密| svipshipincom国产片| 777久久人妻少妇嫩草av网站| 五月开心婷婷网| 成人国语在线视频| 丝瓜视频免费看黄片| 精品免费久久久久久久清纯 | 成人av一区二区三区在线看 | tocl精华| 搡老熟女国产l中国老女人| 纯流量卡能插随身wifi吗| 亚洲欧美精品自产自拍| 黄色a级毛片大全视频| 在线观看人妻少妇| 免费看十八禁软件| 美女国产高潮福利片在线看| 国产精品.久久久| 天堂8中文在线网| 久久精品国产a三级三级三级| 午夜精品久久久久久毛片777| www日本在线高清视频| 国产精品二区激情视频| 免费av中文字幕在线| 日韩制服丝袜自拍偷拍| 亚洲色图 男人天堂 中文字幕| 午夜日韩欧美国产| 男女床上黄色一级片免费看| 一级片'在线观看视频| 18禁观看日本| 亚洲一区二区三区欧美精品| 国产片内射在线| 久久亚洲精品不卡| 国产色视频综合| 国产在视频线精品| 一个人免费在线观看的高清视频 | 亚洲全国av大片| 最近中文字幕2019免费版| 一级片'在线观看视频| 在线观看一区二区三区激情| 亚洲色图 男人天堂 中文字幕| 丰满迷人的少妇在线观看| 亚洲欧美精品综合一区二区三区| 99精国产麻豆久久婷婷| 精品第一国产精品| tube8黄色片| 菩萨蛮人人尽说江南好唐韦庄| 久久久久国产精品人妻一区二区| 久久中文字幕一级| 中文字幕色久视频| 老司机在亚洲福利影院| 国产日韩欧美亚洲二区| 国产精品久久久人人做人人爽| 精品久久久久久久毛片微露脸 | 亚洲久久久国产精品| 久久精品国产亚洲av高清一级| 亚洲伊人久久精品综合| 精品国产乱子伦一区二区三区 | 黄片小视频在线播放| 天堂中文最新版在线下载| 国产麻豆69| 一区二区三区四区激情视频| 国产视频一区二区在线看| 免费少妇av软件| 俄罗斯特黄特色一大片| 久久香蕉激情| 午夜久久久在线观看| 国产成人免费无遮挡视频| 在线观看www视频免费| 老司机影院毛片| 伦理电影免费视频| 国产无遮挡羞羞视频在线观看| 交换朋友夫妻互换小说| 国产成人免费无遮挡视频| 久久亚洲精品不卡| 成人国语在线视频| 国产精品香港三级国产av潘金莲| 午夜老司机福利片| 国产亚洲一区二区精品| 肉色欧美久久久久久久蜜桃| 欧美 亚洲 国产 日韩一| 一二三四社区在线视频社区8| 欧美黄色淫秽网站| 在线观看www视频免费| 欧美成狂野欧美在线观看| 亚洲精品国产av成人精品| 国产精品二区激情视频| 老司机午夜十八禁免费视频| 老熟女久久久| 久久香蕉激情| 99久久综合免费| 亚洲性夜色夜夜综合| 色视频在线一区二区三区| 国产xxxxx性猛交| 黄片大片在线免费观看| 建设人人有责人人尽责人人享有的| 91av网站免费观看| 精品少妇内射三级| 一区二区日韩欧美中文字幕| 人人妻人人澡人人看| 精品人妻1区二区| 国产精品成人在线| 一级a爱视频在线免费观看| 亚洲精品粉嫩美女一区| 制服人妻中文乱码| 久久狼人影院| 在线av久久热| 国产主播在线观看一区二区| 性少妇av在线| 12—13女人毛片做爰片一| 人人妻人人澡人人看| 99热全是精品| 天天操日日干夜夜撸| av天堂久久9| 91精品三级在线观看| 欧美另类一区| 男女无遮挡免费网站观看| 日本猛色少妇xxxxx猛交久久| 男女下面插进去视频免费观看| 日本av手机在线免费观看| 国产欧美日韩一区二区三 | 国产免费视频播放在线视频| 亚洲男人天堂网一区| 久久久久久久精品精品| 亚洲,欧美精品.| 色综合欧美亚洲国产小说| 啦啦啦啦在线视频资源| 18禁黄网站禁片午夜丰满| 一级a爱视频在线免费观看| www.999成人在线观看| 日日夜夜操网爽| 热99久久久久精品小说推荐| 午夜久久久在线观看| 午夜两性在线视频| 99香蕉大伊视频| 久久久久国产精品人妻一区二区| 男女下面插进去视频免费观看| 亚洲专区国产一区二区| 99国产综合亚洲精品| 午夜免费鲁丝| 91麻豆精品激情在线观看国产 | 男男h啪啪无遮挡| 欧美午夜高清在线| 欧美另类亚洲清纯唯美| 丝袜喷水一区| 操出白浆在线播放| 精品亚洲乱码少妇综合久久| 亚洲熟女精品中文字幕| 麻豆乱淫一区二区| 亚洲七黄色美女视频| 亚洲午夜精品一区,二区,三区| 老司机影院毛片| 国产av国产精品国产| 亚洲人成77777在线视频| 色视频在线一区二区三区| 久久精品aⅴ一区二区三区四区| av福利片在线| 国产无遮挡羞羞视频在线观看| 亚洲一码二码三码区别大吗| 热re99久久国产66热| 欧美黄色片欧美黄色片| 丰满人妻熟妇乱又伦精品不卡| 黄片小视频在线播放| 黑人猛操日本美女一级片| 亚洲欧美色中文字幕在线| 国产黄色免费在线视频| 中文字幕制服av| 亚洲熟女精品中文字幕| 久久精品国产a三级三级三级| 国产成人av教育| 国产精品麻豆人妻色哟哟久久| 99香蕉大伊视频| 国产日韩欧美亚洲二区| 国产精品欧美亚洲77777| 美女中出高潮动态图| 欧美大码av| 久久毛片免费看一区二区三区| av国产精品久久久久影院| 亚洲av日韩精品久久久久久密| 精品第一国产精品| 90打野战视频偷拍视频| 色视频在线一区二区三区| 99re6热这里在线精品视频| 欧美日韩国产mv在线观看视频| 久久精品人人爽人人爽视色| 夜夜骑夜夜射夜夜干| 中文字幕人妻丝袜制服| 一级毛片精品| av有码第一页| 国产精品亚洲av一区麻豆| 久久久精品免费免费高清| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品一二三| 免费久久久久久久精品成人欧美视频| 亚洲,欧美精品.| 亚洲av电影在线观看一区二区三区| 大香蕉久久网| 一区二区三区四区激情视频| 51午夜福利影视在线观看| 欧美 亚洲 国产 日韩一| 亚洲av国产av综合av卡| 一边摸一边做爽爽视频免费| 视频区图区小说| 亚洲专区中文字幕在线| 日本猛色少妇xxxxx猛交久久| 国产淫语在线视频| 91老司机精品| 日本精品一区二区三区蜜桃| 一本综合久久免费| 久久中文字幕一级| 久久国产精品影院| 18在线观看网站| 久久久久久人人人人人| 久久精品久久久久久噜噜老黄| 亚洲第一av免费看| 欧美另类亚洲清纯唯美| 19禁男女啪啪无遮挡网站| 久久久久久久精品精品| 亚洲欧美精品自产自拍| 成年人午夜在线观看视频| 欧美成狂野欧美在线观看| 啦啦啦啦在线视频资源| 免费不卡黄色视频| 在线观看人妻少妇| 搡老乐熟女国产| 他把我摸到了高潮在线观看 | 亚洲欧美一区二区三区久久| 一级片免费观看大全| 搡老熟女国产l中国老女人| 中亚洲国语对白在线视频| 成人av一区二区三区在线看 | 成年av动漫网址| 欧美日韩一级在线毛片| 成人国产一区最新在线观看| 99精品久久久久人妻精品| 日韩制服丝袜自拍偷拍| a级毛片在线看网站| av不卡在线播放| 亚洲色图 男人天堂 中文字幕| 久久久国产欧美日韩av| 人妻久久中文字幕网| 成年人免费黄色播放视频| 国产成人免费无遮挡视频| 母亲3免费完整高清在线观看| 欧美日韩精品网址| 久久这里只有精品19| 亚洲中文字幕日韩| 女人精品久久久久毛片| 午夜福利视频在线观看免费| 搡老乐熟女国产| 午夜福利,免费看| 亚洲欧美清纯卡通| 在线看a的网站| 下体分泌物呈黄色| 丝袜人妻中文字幕| 午夜两性在线视频| 最新在线观看一区二区三区| 国产亚洲精品第一综合不卡| 精品第一国产精品| 日韩三级视频一区二区三区| a级毛片黄视频| 岛国毛片在线播放| 色婷婷久久久亚洲欧美| 久久免费观看电影| 午夜视频精品福利| 亚洲欧美成人综合另类久久久| 久久久欧美国产精品| 青青草视频在线视频观看| 狠狠婷婷综合久久久久久88av| 欧美国产精品va在线观看不卡| 亚洲熟女精品中文字幕| 每晚都被弄得嗷嗷叫到高潮| 大码成人一级视频| 无遮挡黄片免费观看| 亚洲av电影在线观看一区二区三区| 精品少妇黑人巨大在线播放| 欧美人与性动交α欧美软件| 日韩视频一区二区在线观看| 水蜜桃什么品种好| 男女之事视频高清在线观看| 国产成人欧美| 久久久国产成人免费| 日韩电影二区| 一区二区三区四区激情视频| 亚洲九九香蕉| 欧美日韩黄片免| 国产精品美女特级片免费视频播放器 | 特级一级黄色大片| 国产精品久久久久久人妻精品电影| 国内少妇人妻偷人精品xxx网站 | 亚洲一区高清亚洲精品| 五月玫瑰六月丁香| 老司机在亚洲福利影院| 97超级碰碰碰精品色视频在线观看| 午夜视频精品福利| 久久伊人香网站| 亚洲片人在线观看| 国产人伦9x9x在线观看| 免费看美女性在线毛片视频| 国语自产精品视频在线第100页| 免费在线观看亚洲国产| 巨乳人妻的诱惑在线观看| 免费在线观看日本一区| 免费av毛片视频| 久久热在线av| 国产三级黄色录像| 伦理电影免费视频| 久久亚洲精品不卡| 精品久久久久久久久久免费视频| 中文字幕人妻丝袜一区二区| 中文字幕高清在线视频| 亚洲av第一区精品v没综合| 正在播放国产对白刺激| 国产高清视频在线观看网站| 好看av亚洲va欧美ⅴa在| 久久午夜亚洲精品久久| 亚洲国产精品sss在线观看|