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

    基于PCA挖掘的SLE疾病的數(shù)據(jù)相關(guān)性研究

    2018-01-02 08:44:34魏汝哲王劍平歐陽(yáng)鑫楊曉洪車國(guó)霖
    軟件 2017年12期
    關(guān)鍵詞:原始數(shù)據(jù)降維理工大學(xué)

    魏汝哲,王劍平,付 萍,張 果,歐陽(yáng)鑫,楊曉洪,車國(guó)霖

    (1. 昆明理工大學(xué) 信息工程與自動(dòng)化學(xué)院,云南 昆明 650500;2. 昆明醫(yī)科大第二附屬醫(yī)院 風(fēng)濕免疫科,云南 昆明 650500)

    基于PCA挖掘的SLE疾病的數(shù)據(jù)相關(guān)性研究

    魏汝哲1,王劍平1,付 萍2,張 果1,歐陽(yáng)鑫1,楊曉洪1,車國(guó)霖1

    (1. 昆明理工大學(xué) 信息工程與自動(dòng)化學(xué)院,云南 昆明 650500;2. 昆明醫(yī)科大第二附屬醫(yī)院 風(fēng)濕免疫科,云南 昆明 650500)

    本文針對(duì) SLE患者并發(fā)繼發(fā)性腎炎發(fā)病原因復(fù)雜多樣、影響因素眾多的問題,提出了一種主成分分析降維算法,用于分析這些數(shù)據(jù),找出最相關(guān)指標(biāo)及相關(guān)規(guī)律。運(yùn)用 python軟件運(yùn)行主成分分析算法的程序,降維得出影響SLE患者并發(fā)繼發(fā)性腎炎主要指標(biāo)。結(jié)果表明,利用主成分分析算法可以找出影響SLE患者并發(fā)繼發(fā)性腎炎主要指標(biāo),為診治提供科學(xué)依據(jù)。

    SLE患者并發(fā)繼發(fā)性腎炎,主成分分析(pca),降維

    0 引言

    隨著信息時(shí)代的到來(lái), 數(shù)據(jù)挖掘被越來(lái)越多地應(yīng)用于臨床實(shí)踐。利用信息技術(shù),醫(yī)療記錄和隨訪數(shù)據(jù)可以更有效地被存儲(chǔ)和提取。同時(shí),從醫(yī)學(xué)數(shù)據(jù)中尋找潛在的關(guān)系或規(guī)律,從而獲得有效的對(duì)病人進(jìn)行診斷、治療的知識(shí);增加對(duì)疾病的預(yù)測(cè)準(zhǔn)確性,在早期發(fā)現(xiàn)疾病,提高治愈率。

    主成分分析在醫(yī)療方面主要應(yīng)用于中藥的成分吸收的藥效,任衛(wèi)高等在主成分分析用于清肺消痤凝膠體外經(jīng)皮滲透的研究[1]中,主成分分析法可以較全面的反應(yīng)接受液中藥物的信息,以此建立的評(píng)價(jià)方法可準(zhǔn)確的優(yōu)選出透皮吸收促進(jìn)劑,并為后續(xù)中藥復(fù)方經(jīng)皮給藥提供了可借鑒的研究思路。張衛(wèi)國(guó)等在基于PCA多特征融合的肝臟B超臨床醫(yī)學(xué)圖像識(shí)別研究與實(shí)現(xiàn)研究[2]中,針對(duì)任何單一特征都不能完整地表示醫(yī)學(xué)圖像的信息,提取醫(yī)學(xué)圖像的顏色特征、紋理特征以及區(qū)域形狀特征,更多地保留圖像的各種信息。并對(duì)提取的特征利用主成分分析(PCA)方法進(jìn)行特征級(jí)的數(shù)據(jù)融合。

    系統(tǒng)性紅斑狼瘡(Systemic Lupus Erythematosus,SLE)是一種臨床表現(xiàn)涉及多種器官、多個(gè)系統(tǒng)受損的慢性系統(tǒng)性自身性免疫性疾病,其具有臨床表現(xiàn)復(fù)雜多樣,輕重緩急不一,免疫調(diào)節(jié)功能紊亂,血清中含有以抗核抗體為主的多種致病性自身抗體,病程多呈現(xiàn)出病情緩解和急性發(fā)作交替等特點(diǎn)[3-4]。SLE是一種育齡期女性好發(fā)、多臟器累及的自身免疫性疾病。狼瘡性腎炎(LN)是指在SLE的基礎(chǔ)上有腎臟疾病臨床表現(xiàn)和(或)伴有腎功能異常,或僅在腎活檢時(shí)發(fā)現(xiàn)有腎小球腎炎病變的系統(tǒng)性紅斑狼瘡患者。LN實(shí)際上是一種腎病,即病人的腎臟器官發(fā)生的病變,但是當(dāng)這種腎病和SLE結(jié)合在一起的時(shí)候就變成了一種非常難以診斷同時(shí)又難以治療的繼發(fā)性疾病[5-7]。LN是SLE繼發(fā)性疾病中發(fā)病率最高的,同時(shí)由于其發(fā)病原因復(fù)雜多樣、影響因素眾多而變得非常難以確診。因此,尋找影響 LN的主要指標(biāo),對(duì) LN疾病的診斷具有非常重要的醫(yī)學(xué)意義。

    PCA是一種掌握主要矛盾的統(tǒng)計(jì)分析方法 ,能夠通過簡(jiǎn)化數(shù)據(jù)(即用較少的綜合指標(biāo)代替原來(lái)具有一定相關(guān)性的較多的指標(biāo))來(lái)反映原來(lái)多變量的大部分信息。本文從昆明醫(yī)學(xué)院第二附屬醫(yī)院風(fēng)濕免疫科近10年關(guān)于系統(tǒng)性紅斑狼瘡的病例中,提取出了758個(gè)已確診患有SLE并發(fā)腎炎的病例。由于該病是SLE繼發(fā)性疾病中發(fā)病率最高的,同時(shí)由于其發(fā)病原因復(fù)雜多樣、影響因素眾多而變得非常難以醫(yī)治,而且每個(gè)病人的整理病歷有280個(gè)指標(biāo),醫(yī)生很難從中挑選出比較重要的指標(biāo)對(duì)病人進(jìn)行醫(yī)治,對(duì)病人的治療造成很大的困難。

    本文通過運(yùn)用 pca降維方法,對(duì)病人的指標(biāo)數(shù)據(jù)進(jìn)行降維,通過降維來(lái)尋找影響LN的主要指標(biāo),找出的指標(biāo)來(lái)對(duì)醫(yī)生提供幫助制定治療方案。

    1 主成分分析算法(pca)概述

    主成分分析(Principal Component Analysis, PCA)是一種常用的數(shù)據(jù)分析方法。PCA通過線性變換將原始數(shù)據(jù)變換為一組各維度線性無(wú)關(guān)的表示,可用于提取數(shù)據(jù)的主要特征分量,常用于高維數(shù)據(jù)的降維。在PCA中,數(shù)據(jù)從原來(lái)的坐標(biāo)系轉(zhuǎn)換到了新的坐標(biāo)系,新坐標(biāo)系的選擇是由數(shù)據(jù)本身決定的。第一個(gè)新坐標(biāo)軸選擇的是原始數(shù)據(jù)中方差最大的方向,第二個(gè)新坐標(biāo)軸的選擇和第一個(gè)新坐標(biāo)軸正交而且具有最大方差的方向。該過程一直重復(fù),重復(fù)次數(shù)為原始數(shù)據(jù)中特征的數(shù)目。我們會(huì)發(fā)現(xiàn),大部分方差都包含在最前面的幾個(gè)新坐標(biāo)軸中。因此我們可以忽略余下的坐標(biāo)軸,即對(duì)數(shù)據(jù)進(jìn)行了降維處理。

    設(shè)F1表示原變量的第一個(gè)線性組合所形成的主成分指標(biāo),即由數(shù)學(xué)知識(shí)可知,每一個(gè)主成分所提取的信息量可用其方差來(lái)度量,其方差 V ar(F1) 越大,表示 F1包含的信息越多。常常希望第一主成分 F1所含的信息量最大,因此在所有的線性組合中選取 F1應(yīng)該是的所有線性組合中方差最大的,故稱 F1為第一主成分。如果第一主成分不足以代表原來(lái)p個(gè)指標(biāo)的信息,再考慮選取第二個(gè)主成分指標(biāo) F2,為有效地反應(yīng)原信息, F1已有的信息就不需要再出現(xiàn)在 F2中,即 F2是與 F1不相關(guān)的 X1, X2,… ,Xp的所有線性組合中方差最大的,故稱 F2為第二主成分,依此類推構(gòu)造出的為原變量指標(biāo)第一、第二、…、第m個(gè)主成分。

    根據(jù)以上分析得知:

    (1) Fi與 Fj互不相關(guān),即 C ov( Fi, Fj) = 0

    2 基于PCA的sle的數(shù)據(jù)挖掘

    2.1 實(shí)驗(yàn)環(huán)境

    本文實(shí)驗(yàn)中所用的計(jì)算機(jī)是 Intel (R) Pentium(R) CPU G2020 @ 2.90 GHz,8 GB內(nèi)存,500 GB硬盤,Windows操作系統(tǒng)。在python環(huán)境下實(shí)現(xiàn)實(shí)驗(yàn)。

    2.2 數(shù)據(jù)來(lái)源及預(yù)處理

    本文從昆明醫(yī)學(xué)院第二附屬醫(yī)院風(fēng)濕免疫科近10年關(guān)于系統(tǒng)性紅斑狼瘡的病例中,提取出758個(gè)已確診患有SLE并發(fā)腎炎的病例。每個(gè)病人的整理病歷有280個(gè)指標(biāo),但是該數(shù)據(jù)包含有很多的缺失值。處理數(shù)據(jù)中的缺失值,我們有一些可選的做法:1.使用可用特征的均值來(lái)填補(bǔ)缺失值;2.使用特殊值來(lái)填補(bǔ)缺失值,如-1;3.忽略有缺失值的樣本。在280個(gè)指標(biāo)里幾乎所有的樣本都有缺失值,因此去除不完整的樣本不太現(xiàn)實(shí)。在后面的算法中要進(jìn)行去均值化,所以也沒有采取特征的均值來(lái)填補(bǔ)缺失值。這里選擇實(shí)數(shù)0來(lái)替換所有的缺失值。

    本文在實(shí)驗(yàn)仿真中使用了 python第三方庫(kù)numpy,由于每組數(shù)據(jù)數(shù)值差異較大,有可能影響預(yù)測(cè)效果,因此需要對(duì)數(shù)據(jù)進(jìn)行歸一化處理,使每組數(shù)據(jù)的參數(shù)都在 0~1 之間。

    2.3 pca降維

    用代碼實(shí)現(xiàn)pca降維的流程圖如下:

    圖1 pca降維流程圖Fig.1 Pca dimension reduction flow chart

    Python代碼如下:

    def loadDataSet(fileName, delim=' '):

    fr = open(fileName)

    stringArr = [line.strip().split(delim) for line in fr.readlines()]

    datArr = [map(float,line) for line in stringArr]

    return mat(datArr)

    def pca(dataMat, topNfeat=9999999):

    meanVals = mean(dataMat, axis=0)

    meanRemoved = dataMat - meanVals

    covMat = cov(meanRemoved, rowvar=0)

    eigVals,eigVects = linalg.eig(mat(covMat))

    eigValInd = argsort(eigVals)

    eigValInd = eigValInd[:-(topNfeat+1):-1]

    redEigVects = eigVects[:,eigValInd1]

    lowDDataMat = meanRemoved * redEigVects

    reconMat = (lowDDataMat * redEigVects.T)+ meanVals

    return lowDDataMat, reconMat,eigVals,eigValInd

    loadDataSet()函數(shù)是用來(lái)將數(shù)據(jù)轉(zhuǎn)換為 pca()函數(shù)可以操作的數(shù)據(jù)集。

    pca()函數(shù)有兩個(gè)參數(shù):第一個(gè)參數(shù)是調(diào)用loadDataSet()后得到的數(shù)據(jù)集,第二個(gè)參數(shù)topNfeat則是一個(gè)可選參數(shù),即應(yīng)用的N個(gè)特征值。如果不指定 topNfeat的值,那么函數(shù)將會(huì)返回前 9999999個(gè)特征值,或者原始數(shù)據(jù)中的全部的特征。這里我們根據(jù)和醫(yī)院的交流選取前5個(gè)特征。

    首先計(jì)算并減去原始數(shù)據(jù)集的平均值,mean()函數(shù)求平均值。然后通過cov()函數(shù)來(lái)計(jì)算協(xié)方差矩陣,Linalg.eig()函數(shù)計(jì)算其特征值和特征向量,接著利用啊如果 argsort()函數(shù)對(duì)特征值進(jìn)行從小到大的排序。根據(jù)特征值排序結(jié)果的逆序就可以得到topNfeat個(gè)最大的特征向量。這些特征向量將構(gòu)成后面對(duì)數(shù)據(jù)進(jìn)行轉(zhuǎn)換的矩陣,該矩陣則利用N個(gè)特征將原始數(shù)據(jù)轉(zhuǎn)換到新空間中。最后,原始數(shù)據(jù)被重構(gòu)后返回用于調(diào)試,同時(shí)降維之后的數(shù)據(jù)集也被返回了。

    3 挖掘的結(jié)果及分析

    本文是通過降維來(lái)尋找影響 LN的主要指標(biāo),通過 pca降維之后數(shù)據(jù)轉(zhuǎn)換到了新的空間,所以用來(lái)構(gòu)建新的空間的topNfeat個(gè)最大的特征向量是影響LN的主要指標(biāo),topNfeat個(gè)最大的特征向量也就是topNfeat個(gè)最大的特征值。

    pca()函數(shù)最后輸出的 eigValInd就是 topNfeat個(gè)最大的特征值的下標(biāo),通過數(shù)據(jù)可視化,圖1就是最相關(guān)的 5個(gè)理化指標(biāo),這 5個(gè)指標(biāo)分別是dsDNA 抗體、clq 抗體[8-9]、ANCA[10]、SM 抗體[11-12]、ALB。

    圖2 數(shù)據(jù)挖掘結(jié)果Fig.2 Data mining results

    通過實(shí)驗(yàn)得到dsDNA抗體、SM抗體、clq抗體、ANCA抗體、ALB這5個(gè)影響LN的主要指標(biāo),根據(jù)實(shí)驗(yàn)結(jié)果可以幫助醫(yī)生得到最佳的治療方案。

    4 結(jié)語(yǔ)

    本文通過運(yùn)用 pca降維方法,對(duì)病人的指標(biāo)數(shù)據(jù)進(jìn)行降維,通過降維來(lái)尋找影響LN的主要指標(biāo),找出的指標(biāo)來(lái)對(duì)醫(yī)生提供幫助制定治療方案?,F(xiàn)將本文的主要內(nèi)容總結(jié)如下:

    (1)本文對(duì)來(lái)自醫(yī)院的病例資料進(jìn)行數(shù)據(jù)預(yù)處理,包括數(shù)據(jù)清洗、數(shù)據(jù)插補(bǔ)、數(shù)據(jù)變換,為下一步的研究做準(zhǔn)備工作;

    (2)基于預(yù)處理后的數(shù)據(jù),通過python運(yùn)算pca算法,對(duì)原始數(shù)據(jù)降維,獲得影響LN的主要 指標(biāo)。

    (3)將運(yùn)行算法獲得的結(jié)果與醫(yī)學(xué)文獻(xiàn)所得結(jié)果相對(duì)比,驗(yàn)證了算法的正確,也為治療 LN提供更加科學(xué)的依據(jù)。

    [1] 任衛(wèi)高, 彭琳秀, 雷飛飛, 等. 主成分分析用于清肺消痤凝膠體外經(jīng)皮滲透的研究[J]. 中國(guó)中藥雜志. 2015, 40(2):231-235.

    [2] 張衛(wèi)國(guó), 王桂花. 基于PCA多特征融合的肝臟B超臨床醫(yī)學(xué)圖像識(shí)別研究與實(shí)現(xiàn)[J]. 計(jì)算機(jī)應(yīng)用與軟件. 2014, (4):239-243.

    [3] Hochberg MC.Updating the American College of Rheum atology revised criteria for the classification of systemic lupus erythem atosus[J]. Arthritis Rheum, 1997, 40(9): 1725.

    [4] Kozyrev SV, Abelson AK, Wojcik J, et al. Functional variants in the B-cell gene BANK1 are associated with systemic lupus erythematosus[J]. Nat Genet, 2008, 40(2): 211-216.

    [5] Biesecker G, Katz S, Koffler D. Renal localization of the membrane attack complex in systemic lupus erythematosus nephritis[J]. Journal of Experimental Medicine, 1982, 154(6):1779-1794.

    [6] HOUMAN MH, SMITI- KHANFIR M, BEN GHORBELL I,et al.Sysmetic lupus erythematosus in Tunisia: Demographic and clinical analysis of 100 patients[J]. Lupus, 2004, 13(3):204-211.

    [7] RESENDE AL, TITAN SM, BARROS RT, et al. Worse renal outcome of lupus nephritis in male patients: a case- control study[J]. Lupus, 2011, 20(6): 561- 567.

    [8] Botto M, Walport MJ. C1q, autoimmunity and apotosis[J].Immunobiol, 2002, 205: 395-406.

    [9] Braun A, Sis J, Max R, et al. Anti-chromatin and anti-C1q antibodies in systemic lupus erythematosus compared to other systemic autoimmune diseases[J]. Scand J Rheumatol,2007, 36: 291-298.

    [10] 劉軍生, 肖湘成, 陳立平, 等. 狼瘡腎炎抗中性粒細(xì)胞胞漿抗體與臨床表現(xiàn)的相關(guān)性研究[J]. 醫(yī)學(xué)臨床研究, 2006,23(10): 1606-1610.

    [11] MANSON JJ, MA A, ROGERS P, MASON LJ, et al.Relationship between anti-dsDNA, anti- nucleosome and anti- alpha-actinin antibodies and markers of renal disease in patients with lupus nephritis: a prospective longitudinal study[J]. Arthritis Res Ther, 2009, 11(5): 154.

    [12] SUH-LAILAM BB, CHIARO TR, DAVIS K W, et al. Evaluation of a high avidity anti- dsDNA IgG enzyme- linked immunosorbent assay for the diagnosis of systemic lupus erythematosus[J]. Int J Clin Exp Pathol, 2011, 4(8): 748-754.

    Correlation of PCA-based Data Mining SLE Disease

    WEI Ru-zhe, WANG Jian-ping, FU Ping, ZHANG Guo
    (Kunming University of Science and Technology, Kunming 650500, China; Intelligent Measurement and Control Institute, Kunming 650500)

    In this paper, we propose a principal component analysis reduction algorithm to solve the problem of complex and diverse influencing factors of SLE patients with secondary nephritis. The program of running the principal component analysis algorithm using python software was used to reduce the influence of SLE Patients with secondary primary symptoms of secondary nephritis. The results show that the principal component analysis algorithm can find out the main indexes that affect the secondary nephritis in SLE patients, and provide scientific basis for diagnosis and treatment.

    SLE patients with secondary nephritis; Principal component analysis; Dimensionality reduction

    TP181

    A

    10.3969/j.issn.1003-6970.2017.12.018

    本文著錄格式:魏汝哲,王劍平,付萍,等. 基于PCA挖掘的SLE疾病的數(shù)據(jù)相關(guān)性研究[J]. 軟件,2017,38(12):95-97

    國(guó)家自然科學(xué)基金(61364008);云南省應(yīng)用基礎(chǔ)研究重點(diǎn)項(xiàng)目(2014FA029);云南省教育廳重點(diǎn)基金項(xiàng)目(2013Z127);昆明理工大學(xué)人才培養(yǎng)項(xiàng)目(14118596);昆明理工大學(xué)復(fù)雜工業(yè)控制學(xué)科方向團(tuán)隊(duì)建設(shè)計(jì)劃

    魏汝哲(1992-),男,山東聊城人,碩士,主要研究方向?yàn)獒t(yī)療數(shù)據(jù)挖掘,檢測(cè)技術(shù)。單位:昆明理工大學(xué);張果,男,研究生導(dǎo)師,昆明理工大學(xué);歐陽(yáng)鑫,男,研究生導(dǎo)師,昆明理工大學(xué);楊曉洪,女,研究生導(dǎo)師,昆明理工大學(xué);車國(guó)霖,男,研究生導(dǎo)師,昆明理工大學(xué)。

    王劍平,男,研究生導(dǎo)師,昆明理工大學(xué);付萍,女,研究生導(dǎo)師,昆明醫(yī)科大學(xué)。

    猜你喜歡
    原始數(shù)據(jù)降維理工大學(xué)
    Three-Body’s epic scale and fiercely guarded fanbase present challenges to adaptations
    昆明理工大學(xué)
    GOLDEN OPPORTUNITY FOR CHINA-INDONESIA COOPERATION
    受特定變化趨勢(shì)限制的傳感器數(shù)據(jù)處理方法研究
    降維打擊
    海峽姐妹(2019年12期)2020-01-14 03:24:40
    昆明理工大學(xué)
    昆明理工大學(xué)
    浙江理工大學(xué)
    全新Mentor DRS360 平臺(tái)借助集中式原始數(shù)據(jù)融合及直接實(shí)時(shí)傳感技術(shù)實(shí)現(xiàn)5 級(jí)自動(dòng)駕駛
    汽車零部件(2017年4期)2017-07-12 17:05:53
    拋物化Navier-Stokes方程的降維仿真模型
    免费高清视频大片| 久久久久性生活片| 欧美日韩综合久久久久久 | 欧洲精品卡2卡3卡4卡5卡区| 午夜影院日韩av| 97超级碰碰碰精品色视频在线观看| 麻豆久久精品国产亚洲av| 成人无遮挡网站| 国语自产精品视频在线第100页| 淫秽高清视频在线观看| 免费一级毛片在线播放高清视频| 长腿黑丝高跟| 久久久久久久亚洲中文字幕 | 国产日本99.免费观看| 99国产极品粉嫩在线观看| 简卡轻食公司| 两个人的视频大全免费| 能在线免费观看的黄片| 又粗又爽又猛毛片免费看| 中文字幕精品亚洲无线码一区| 亚洲av成人不卡在线观看播放网| 简卡轻食公司| 久久久久亚洲av毛片大全| 最新在线观看一区二区三区| 亚洲国产精品sss在线观看| 老司机福利观看| 搡女人真爽免费视频火全软件 | 级片在线观看| 亚洲国产欧洲综合997久久,| 欧美潮喷喷水| 日韩欧美免费精品| 国产精品一区二区三区四区免费观看 | 综合色av麻豆| 精品一区二区三区视频在线观看免费| 日本黄色视频三级网站网址| 午夜日韩欧美国产| 亚洲美女搞黄在线观看 | 婷婷精品国产亚洲av| 午夜精品在线福利| 人人妻,人人澡人人爽秒播| 亚洲av电影不卡..在线观看| 国产精品一及| 小蜜桃在线观看免费完整版高清| 中文字幕熟女人妻在线| 久久久久国产精品人妻aⅴ院| 一边摸一边抽搐一进一小说| 亚洲中文字幕一区二区三区有码在线看| 成人av一区二区三区在线看| 精品乱码久久久久久99久播| 一级毛片久久久久久久久女| 成年人黄色毛片网站| 成年女人毛片免费观看观看9| 精品久久久久久久久久久久久| 日本一本二区三区精品| 国产精品99久久久久久久久| 国内精品美女久久久久久| 热99在线观看视频| 亚洲av日韩精品久久久久久密| 日韩欧美精品免费久久 | 夜夜看夜夜爽夜夜摸| 亚洲国产精品sss在线观看| 搡老岳熟女国产| 色综合亚洲欧美另类图片| 噜噜噜噜噜久久久久久91| 男女之事视频高清在线观看| 好男人在线观看高清免费视频| 国产亚洲欧美98| 精品人妻一区二区三区麻豆 | 欧美日韩中文字幕国产精品一区二区三区| 性色avwww在线观看| aaaaa片日本免费| 久久精品国产亚洲av天美| 脱女人内裤的视频| 亚洲最大成人中文| 99久久久亚洲精品蜜臀av| 午夜福利欧美成人| 亚洲狠狠婷婷综合久久图片| 日韩 亚洲 欧美在线| 亚洲精品乱码久久久v下载方式| 精品国内亚洲2022精品成人| 成人国产综合亚洲| 亚洲美女黄片视频| 久久久精品欧美日韩精品| 亚洲最大成人手机在线| 久久草成人影院| 日日摸夜夜添夜夜添av毛片 | 亚洲av成人不卡在线观看播放网| 男女床上黄色一级片免费看| 日韩欧美免费精品| bbb黄色大片| 亚洲在线自拍视频| 熟妇人妻久久中文字幕3abv| 此物有八面人人有两片| 一进一出好大好爽视频| 成年女人永久免费观看视频| 女人十人毛片免费观看3o分钟| bbb黄色大片| 在线播放国产精品三级| 国产精品亚洲美女久久久| 女同久久另类99精品国产91| 精品人妻1区二区| 丰满的人妻完整版| 舔av片在线| 在现免费观看毛片| 欧美成人免费av一区二区三区| 亚洲人成网站高清观看| 免费av观看视频| 一个人看的www免费观看视频| 精品人妻偷拍中文字幕| 国内揄拍国产精品人妻在线| 亚洲熟妇中文字幕五十中出| 90打野战视频偷拍视频| www.www免费av| 女人被狂操c到高潮| 国产色婷婷99| 国产探花极品一区二区| 国产午夜福利久久久久久| 亚洲自拍偷在线| 国产精品,欧美在线| 少妇裸体淫交视频免费看高清| 18禁裸乳无遮挡免费网站照片| 国产精品伦人一区二区| 99久久99久久久精品蜜桃| 无遮挡黄片免费观看| 夜夜爽天天搞| 可以在线观看毛片的网站| 日本五十路高清| 欧美日韩乱码在线| 国产aⅴ精品一区二区三区波| 国产高潮美女av| 国产精品爽爽va在线观看网站| 免费看日本二区| 亚洲av五月六月丁香网| 亚洲精品粉嫩美女一区| 91麻豆av在线| 色哟哟哟哟哟哟| 我要看日韩黄色一级片| 欧美不卡视频在线免费观看| 国产伦精品一区二区三区视频9| 久久天躁狠狠躁夜夜2o2o| 国内精品久久久久精免费| 亚洲精品亚洲一区二区| 免费搜索国产男女视频| 国产欧美日韩精品一区二区| 一a级毛片在线观看| 99riav亚洲国产免费| 丁香六月欧美| 成人国产一区最新在线观看| 90打野战视频偷拍视频| 精品久久久久久久久久久久久| 国产午夜精品论理片| 久久久久亚洲av毛片大全| 99国产精品一区二区蜜桃av| 成人欧美大片| 亚洲美女视频黄频| 能在线免费观看的黄片| 黄色日韩在线| 亚洲性夜色夜夜综合| 国产成人啪精品午夜网站| 国产一区二区三区视频了| 亚洲最大成人av| 人妻久久中文字幕网| 亚洲第一区二区三区不卡| 90打野战视频偷拍视频| 在线播放国产精品三级| 午夜影院日韩av| 久久午夜福利片| 国产av麻豆久久久久久久| 久久国产乱子伦精品免费另类| 69av精品久久久久久| 久久久久久久亚洲中文字幕 | 99久国产av精品| 国产亚洲精品久久久久久毛片| 久久久久久久午夜电影| 内地一区二区视频在线| 日韩免费av在线播放| 天美传媒精品一区二区| av在线老鸭窝| 成人三级黄色视频| 亚洲国产色片| 亚洲国产精品合色在线| 淫秽高清视频在线观看| 国产精品一区二区免费欧美| 欧美bdsm另类| aaaaa片日本免费| h日本视频在线播放| 亚洲内射少妇av| 又爽又黄无遮挡网站| 亚洲人成网站在线播放欧美日韩| 69人妻影院| 亚洲av成人精品一区久久| 午夜免费成人在线视频| 美女cb高潮喷水在线观看| 在线播放无遮挡| 久久精品夜夜夜夜夜久久蜜豆| 国产欧美日韩一区二区精品| 亚洲专区国产一区二区| 舔av片在线| 日韩欧美国产在线观看| 亚洲国产精品sss在线观看| 国产亚洲av嫩草精品影院| 国产亚洲精品av在线| www.色视频.com| 国产成人av教育| 精品福利观看| 在线观看美女被高潮喷水网站 | 热99在线观看视频| 成人鲁丝片一二三区免费| 99视频精品全部免费 在线| 欧美一区二区国产精品久久精品| 中文在线观看免费www的网站| 午夜福利在线在线| 国产精品嫩草影院av在线观看 | 永久网站在线| 亚洲黑人精品在线| 国产精品不卡视频一区二区 | 午夜日韩欧美国产| 欧美精品国产亚洲| av黄色大香蕉| 色噜噜av男人的天堂激情| 一级作爱视频免费观看| 色5月婷婷丁香| 午夜福利18| 精品免费久久久久久久清纯| 午夜精品久久久久久毛片777| 亚洲人成网站高清观看| 999久久久精品免费观看国产| 国产av不卡久久| 99在线视频只有这里精品首页| 麻豆一二三区av精品| 在线观看av片永久免费下载| 国产精品三级大全| 国产日本99.免费观看| 国产毛片a区久久久久| 亚洲精品色激情综合| 国产欧美日韩一区二区精品| 国产色婷婷99| 日本黄色视频三级网站网址| 深夜精品福利| 非洲黑人性xxxx精品又粗又长| 99久久99久久久精品蜜桃| 亚洲美女黄片视频| 制服丝袜大香蕉在线| 男人狂女人下面高潮的视频| 搡老熟女国产l中国老女人| 在线观看66精品国产| 久久99热这里只有精品18| 国产精品自产拍在线观看55亚洲| 久久久久性生活片| av专区在线播放| 欧美三级亚洲精品| 男女做爰动态图高潮gif福利片| 免费看美女性在线毛片视频| 一级黄片播放器| 搡老熟女国产l中国老女人| 人人妻人人澡欧美一区二区| 偷拍熟女少妇极品色| 中文字幕精品亚洲无线码一区| 真人做人爱边吃奶动态| 欧美一区二区精品小视频在线| 国产一区二区亚洲精品在线观看| x7x7x7水蜜桃| 欧美日韩国产亚洲二区| 亚洲国产精品999在线| 精品一区二区三区视频在线| 欧美区成人在线视频| 免费看美女性在线毛片视频| 国产精品人妻久久久久久| 美女黄网站色视频| 精品久久久久久成人av| 亚洲欧美日韩无卡精品| 国产真实乱freesex| 日本成人三级电影网站| 欧美成人a在线观看| 国产伦一二天堂av在线观看| 有码 亚洲区| 在线播放国产精品三级| 久久亚洲精品不卡| 校园春色视频在线观看| 欧美又色又爽又黄视频| 欧美日韩福利视频一区二区| 免费在线观看影片大全网站| 在现免费观看毛片| 国产午夜福利久久久久久| 国产精品久久久久久人妻精品电影| 韩国av一区二区三区四区| 精品午夜福利在线看| 国产久久久一区二区三区| 精品久久久久久成人av| 日韩中字成人| 国产成人欧美在线观看| 国产三级在线视频| 免费在线观看亚洲国产| 国产精品久久久久久亚洲av鲁大| 亚洲人成电影免费在线| 国产精品综合久久久久久久免费| 国产淫片久久久久久久久 | 十八禁网站免费在线| 在线播放国产精品三级| www.色视频.com| 午夜福利成人在线免费观看| 成人三级黄色视频| 国产成人欧美在线观看| 日韩成人在线观看一区二区三区| 首页视频小说图片口味搜索| 国产男靠女视频免费网站| 简卡轻食公司| 午夜福利高清视频| 久久久国产成人免费| 夜夜看夜夜爽夜夜摸| 亚洲五月天丁香| 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 在线播放无遮挡| 国产亚洲欧美在线一区二区| 亚洲三级黄色毛片| 亚洲18禁久久av| 最近最新中文字幕大全电影3| 亚洲美女视频黄频| 久久精品国产亚洲av天美| 欧美在线黄色| 精品99又大又爽又粗少妇毛片 | 制服丝袜大香蕉在线| 亚洲成av人片在线播放无| 亚洲精品粉嫩美女一区| 真人做人爱边吃奶动态| 欧美日韩中文字幕国产精品一区二区三区| 麻豆一二三区av精品| 亚洲综合色惰| 在现免费观看毛片| 久久久久久久久中文| 免费看a级黄色片| 夜夜爽天天搞| 99热精品在线国产| 午夜a级毛片| 国产成人aa在线观看| 在线十欧美十亚洲十日本专区| 波多野结衣高清作品| 最新在线观看一区二区三区| 中文字幕高清在线视频| 国产高清视频在线播放一区| 天堂动漫精品| 性色avwww在线观看| 别揉我奶头 嗯啊视频| 亚洲一区二区三区色噜噜| 俄罗斯特黄特色一大片| 色av中文字幕| 日韩欧美 国产精品| 亚洲av成人av| 好男人在线观看高清免费视频| 日韩欧美国产在线观看| 免费av不卡在线播放| 日本 av在线| 亚洲国产精品成人综合色| 亚洲在线观看片| 亚洲美女黄片视频| 他把我摸到了高潮在线观看| 天堂√8在线中文| 听说在线观看完整版免费高清| 91在线精品国自产拍蜜月| www.www免费av| 欧美xxxx性猛交bbbb| 如何舔出高潮| 国产aⅴ精品一区二区三区波| 国产精品野战在线观看| 国产免费男女视频| 三级男女做爰猛烈吃奶摸视频| 国产免费一级a男人的天堂| 午夜影院日韩av| 欧美午夜高清在线| 亚洲内射少妇av| 亚洲七黄色美女视频| 成年人黄色毛片网站| 亚洲成a人片在线一区二区| 欧洲精品卡2卡3卡4卡5卡区| 国产色婷婷99| 国产免费男女视频| 欧美三级亚洲精品| 久久久久亚洲av毛片大全| 在线观看一区二区三区| 国产高清视频在线观看网站| av在线天堂中文字幕| 99久久无色码亚洲精品果冻| 亚洲国产色片| 亚洲av电影不卡..在线观看| 午夜精品一区二区三区免费看| 国产熟女xx| 亚洲中文字幕日韩| 2021天堂中文幕一二区在线观| 精品一区二区免费观看| 久久久久久大精品| 91久久精品国产一区二区成人| 亚洲国产色片| 中文字幕免费在线视频6| 成人一区二区视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区激情视频| av在线观看视频网站免费| 国产91精品成人一区二区三区| 亚洲狠狠婷婷综合久久图片| 18美女黄网站色大片免费观看| 99热6这里只有精品| 国产aⅴ精品一区二区三区波| 韩国av一区二区三区四区| 国产精品久久久久久久电影| 神马国产精品三级电影在线观看| 久久亚洲精品不卡| or卡值多少钱| 禁无遮挡网站| 欧美极品一区二区三区四区| 真人做人爱边吃奶动态| av福利片在线观看| 亚洲美女黄片视频| 亚洲av成人av| 亚洲精品影视一区二区三区av| 51国产日韩欧美| 国产高清三级在线| 国产精品一区二区三区四区免费观看 | 亚洲,欧美,日韩| av在线蜜桃| 亚洲国产欧美人成| 午夜福利成人在线免费观看| 欧美一区二区精品小视频在线| 欧美黑人巨大hd| x7x7x7水蜜桃| 人妻丰满熟妇av一区二区三区| 婷婷亚洲欧美| 亚洲成人精品中文字幕电影| АⅤ资源中文在线天堂| 一区二区三区免费毛片| 欧美中文日本在线观看视频| 免费观看精品视频网站| 亚洲专区国产一区二区| 无遮挡黄片免费观看| 国产一区二区激情短视频| 嫁个100分男人电影在线观看| 午夜福利成人在线免费观看| 在线播放无遮挡| 日韩免费av在线播放| 88av欧美| 最后的刺客免费高清国语| 欧美极品一区二区三区四区| 搡老熟女国产l中国老女人| 国产精品综合久久久久久久免费| 国产91精品成人一区二区三区| 99国产精品一区二区蜜桃av| 国产69精品久久久久777片| 99久久精品一区二区三区| 永久网站在线| 黄片小视频在线播放| 亚洲精品影视一区二区三区av| 免费大片18禁| 一本一本综合久久| 九色国产91popny在线| 在线观看66精品国产| 男女下面进入的视频免费午夜| 日韩欧美免费精品| 国产日本99.免费观看| 中文字幕av在线有码专区| 久久久久久九九精品二区国产| 日韩欧美三级三区| 成人无遮挡网站| 丰满人妻熟妇乱又伦精品不卡| 欧美激情久久久久久爽电影| 亚洲美女黄片视频| 亚洲黑人精品在线| 亚洲,欧美,日韩| 亚洲精品色激情综合| 国产单亲对白刺激| 永久网站在线| а√天堂www在线а√下载| 十八禁国产超污无遮挡网站| 婷婷精品国产亚洲av在线| 啪啪无遮挡十八禁网站| 在线观看舔阴道视频| 日本五十路高清| 欧美日韩乱码在线| 成熟少妇高潮喷水视频| 免费无遮挡裸体视频| av欧美777| 最近最新免费中文字幕在线| 日韩欧美在线二视频| 麻豆成人av在线观看| 免费电影在线观看免费观看| 国产精华一区二区三区| 夜夜躁狠狠躁天天躁| 国产精品综合久久久久久久免费| 琪琪午夜伦伦电影理论片6080| 三级毛片av免费| 中文资源天堂在线| 搡老岳熟女国产| 免费搜索国产男女视频| av黄色大香蕉| 国产三级在线视频| 日韩 亚洲 欧美在线| 国产亚洲精品久久久久久毛片| 欧美一区二区亚洲| 天堂网av新在线| 精品国产三级普通话版| 国内精品久久久久久久电影| 熟女电影av网| 嫩草影院新地址| 麻豆一二三区av精品| 高清日韩中文字幕在线| 久久久久精品国产欧美久久久| 精品人妻偷拍中文字幕| 欧美日韩乱码在线| 国产精品永久免费网站| 国产精品一区二区性色av| 国产亚洲欧美在线一区二区| 国产亚洲欧美98| 一个人看的www免费观看视频| 亚州av有码| 精品无人区乱码1区二区| 日本熟妇午夜| 我的女老师完整版在线观看| 麻豆成人午夜福利视频| 国产亚洲精品av在线| 最新中文字幕久久久久| 久久九九热精品免费| 在线播放无遮挡| 精品国产三级普通话版| 真人一进一出gif抽搐免费| 韩国av一区二区三区四区| 搡老妇女老女人老熟妇| 亚洲欧美日韩东京热| 日本a在线网址| 欧美精品啪啪一区二区三区| 欧美在线一区亚洲| 久久久久久久亚洲中文字幕 | 国产视频内射| 久久人人精品亚洲av| 日韩亚洲欧美综合| 久久久久免费精品人妻一区二区| 男女下面进入的视频免费午夜| 毛片女人毛片| 9191精品国产免费久久| 久久婷婷人人爽人人干人人爱| 别揉我奶头~嗯~啊~动态视频| 久久久久久久久久成人| 国产精品,欧美在线| 久久精品久久久久久噜噜老黄 | 尤物成人国产欧美一区二区三区| 国内精品美女久久久久久| 观看美女的网站| 国产精品一区二区免费欧美| 97碰自拍视频| 18美女黄网站色大片免费观看| 校园春色视频在线观看| 少妇熟女aⅴ在线视频| 在线免费观看的www视频| 色精品久久人妻99蜜桃| 久久久精品大字幕| 一级a爱片免费观看的视频| 久久午夜亚洲精品久久| 伊人久久精品亚洲午夜| 亚洲熟妇中文字幕五十中出| 亚洲精品456在线播放app | 搡老熟女国产l中国老女人| 91麻豆精品激情在线观看国产| 在线看三级毛片| 亚洲电影在线观看av| 国内毛片毛片毛片毛片毛片| 日本一本二区三区精品| 搡老岳熟女国产| 简卡轻食公司| 欧美日韩综合久久久久久 | 婷婷亚洲欧美| 淫秽高清视频在线观看| 亚洲天堂国产精品一区在线| 国产真实伦视频高清在线观看 | 我要看日韩黄色一级片| 又黄又爽又刺激的免费视频.| 深夜a级毛片| 欧美日韩福利视频一区二区| 波野结衣二区三区在线| av在线天堂中文字幕| 国产色婷婷99| 99久久精品热视频| 欧美黄色淫秽网站| 成人av一区二区三区在线看| 亚洲精品影视一区二区三区av| 国产熟女xx| 哪里可以看免费的av片| 黄色丝袜av网址大全| 亚洲精品日韩av片在线观看| 一级作爱视频免费观看| 又紧又爽又黄一区二区| 天美传媒精品一区二区| 看黄色毛片网站| 长腿黑丝高跟| 中文字幕人妻熟人妻熟丝袜美| 美女免费视频网站| 无人区码免费观看不卡| 精品一区二区三区人妻视频| 18美女黄网站色大片免费观看| 搞女人的毛片| a级毛片a级免费在线| 亚洲三级黄色毛片| 色吧在线观看| 别揉我奶头 嗯啊视频| 亚洲成人中文字幕在线播放| 很黄的视频免费| ponron亚洲| 真人一进一出gif抽搐免费| 国产精品亚洲一级av第二区| 丝袜美腿在线中文| 一个人看视频在线观看www免费| 久久久久久大精品| 国产精品电影一区二区三区| 亚州av有码| 国产aⅴ精品一区二区三区波| 国产主播在线观看一区二区| 婷婷丁香在线五月| 少妇人妻一区二区三区视频|