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

    傾向性評分匹配法在不良反應(yīng)信號檢測中的應(yīng)用*

    2012-12-04 02:59:50許金芳葉小飛杜文民
    中國衛(wèi)生統(tǒng)計 2012年6期
    關(guān)鍵詞:傾向性喹硫平數(shù)據(jù)挖掘

    王 超 吳 騁 許金芳 錢 維 葉小飛 杜文民 賀 佳

    藥品不良反應(yīng)(adverse drug reaction,ADR)監(jiān)測最主要的數(shù)據(jù)來源是自發(fā)呈報系統(tǒng)(spontaneous reporting system,SRS),由各個地區(qū)的衛(wèi)生部門根據(jù)各自監(jiān)測到的不良反應(yīng)報告及時自發(fā)地上報得來,截止2010年該數(shù)據(jù)已達到692 904。要從如此海量而又復雜的數(shù)據(jù)中尋求信息并發(fā)現(xiàn)藥品與不良反應(yīng)之間的關(guān)聯(lián),需要借助有效的數(shù)據(jù)挖掘和統(tǒng)計分析方法才能實現(xiàn)。自1998年以來,對于不良反應(yīng)信號的數(shù)據(jù)挖掘工作已成為WHO烏普薩拉監(jiān)測中心(UMC)信號檢測過程中的一項常規(guī)工作〔1〕。如今,越來越多的國家已經(jīng)在用數(shù)據(jù)挖掘工具來提取信息。比較成熟的挖掘方法包括荷蘭藥物警戒中心的報告比數(shù)比法(reporting odds ratio,ROR);英國藥品和保健產(chǎn)品管理局的綜合標準法(Medicines and Healthcare Products Regulatory Agency MHRA,或 MCA);WHO的貝葉斯可信神經(jīng)網(wǎng)絡(luò)法(Bayesian confidence propagation neuralnetwork,BCPNN)和美國FDA的多項伽馬泊松收縮估計法(Multi-item gamma poisson shrinker,MGPS)等〔2-5〕。

    應(yīng)用數(shù)據(jù)挖掘方法得出的不成比例報告僅僅是知識發(fā)現(xiàn)過程的第一步。對分析觀察性研究的數(shù)據(jù),一個最大的阻礙就是可能存在可測量的或不可測量的混雜,進而扭曲或遮蔽了真實的內(nèi)部關(guān)聯(lián)〔1〕。在ADR監(jiān)測工作中,目前大多情況仍是在基線不均衡的狀態(tài)下分析SRS數(shù)據(jù)。如老年人服用某種藥品發(fā)生的某不良反應(yīng)概率更高,但在SRS數(shù)據(jù)中某段時間內(nèi)上報的老年人很多,因此在沒有將年齡均衡的情況下得出的該藥和不良反應(yīng)的關(guān)聯(lián)是不可信的。如果能有效控制其中可能存在的混雜因素,將會使不同特征的數(shù)據(jù)具有可比性,從而更好地挖掘出數(shù)據(jù)間隱藏的信息。對于分層、匹配或多元回歸等傳統(tǒng)的去除混雜的方法,當面對像SRS這樣混雜因素較多的數(shù)據(jù)時,若同時分析這些混雜因素會出現(xiàn)過度分層或過匹配的現(xiàn)象,logistic回歸在自變量較多的情況下也會出現(xiàn)共線性等問題〔6〕。傾向性評分法(propensity score,PS)作為一種均衡基線的新方法〔7〕經(jīng)Rosenbaum和Rubin于上世紀80年代提出后引起了廣泛關(guān)注,主要應(yīng)用于觀察性和臨床非隨機化數(shù)據(jù)的研究。本文將介紹應(yīng)用PS匹配法對SRS數(shù)據(jù)進行混雜偏倚矯正并探討其適用性與效果。

    原 理

    1.PS的基本原理

    Rosenbaum 和 Rubin 把傾向性評分定義為〔7-8〕:在可觀察的協(xié)變量(Xi)條件下,研究對象i(i=1,2,…,N)接受某種處理(或暴露)因素(Zi=1)而非對照因素(Zi=0)的條件概率,計算公式為:

    可將公式進一步轉(zhuǎn)化為:

    這里的e(x)就是傾向性評分??梢钥闯?,PS概括了所有協(xié)變量的作用并整合為一個綜合的分值。因此如果處理和對照有近似的PS值,我們就可以認為他們在這個具體的PS水平上是均衡可比的。

    2.PS的應(yīng)用

    通過PS去除混雜主要包括以下幾個步驟〔9-10〕:第一,建立一個模型來估計傾向性評分值;第二,評價模型的適用性;第三,通過合適的方法應(yīng)用傾向性評分值估計處理(或暴露)因素的作用。通常應(yīng)用PS的方法主要為匹配,分層和回歸矯正。這三種方法計算PS的過程是相同的,待PS估算后就可應(yīng)用不同的方法來處理。傾向性評分分層法和回歸矯正法簡單易行,其中分層法應(yīng)用較廣,它是將所得的PS值分成若干層,一般分五層就能降低90%以上的偏倚〔11-12〕,再在每一層里進行處理效果的估計,之后加權(quán)這五層的估計值就能得到矯正混雜后的總結(jié)果。許多研究表明〔13-14〕,PS匹配法近年來受到更多的關(guān)注,是一種最能均衡各組之間分布與構(gòu)成的方法。并且由于SRS特點,對照組的人群數(shù)量遠遠大于處理組,此時更適合采用匹配的方法分析PS〔7〕。因此本研究將采用傾向性評分匹配法來處理SRS的數(shù)據(jù)。

    資料和方法

    1.模擬試驗

    為了驗證PS應(yīng)用在SRS中的可行性,本文首先在模擬數(shù)據(jù)集中進行了探討。模擬數(shù)據(jù)集根據(jù)Brookhart〔15〕的模擬多變量數(shù)據(jù)方法以及SRS自身的特點產(chǎn)生:

    首先,利用正態(tài)分布離差隨機數(shù)RANNOR,模擬SRS年齡的變量生成連續(xù)協(xié)變量X1,利用二項分布隨機數(shù)RANBIN模擬性別生成二分類協(xié)變量X2,并且X1,X2相互獨立。

    其次,在X1,X2變量下,應(yīng)用logistic回歸產(chǎn)生暴露變量 Drug1、Drug2、Drug3、Drug4。暴露模型如下:

    最后,以 X1、X2、Drug1、Drug2、Drug3、Drug4為自變量,應(yīng)用logistic回歸生成二分類結(jié)果變量ADR,模型如下:

    基于模擬數(shù)據(jù),通過 ROR、PRR、BCPNN〔2-5〕三種方法分別挖掘數(shù)據(jù),并根據(jù)其各自的標準閾值進行過濾篩選。選擇所有挖掘方法均能檢測出的可疑組合進行下一步的PS矯正。將PS匹配法應(yīng)用于SRS主要分以下幾步:

    第一,應(yīng)用logistic回歸納入候選協(xié)變量建立模型,得出每個個體的傾向性評分值。在選擇協(xié)變量時,我們應(yīng)該盡量選取既與某藥因素有關(guān)又與某不良反應(yīng)變量有關(guān)的協(xié)變量納入模型〔15〕,并且根據(jù)C統(tǒng)計值判斷模型的適用性。

    第二,將包含有傾向評分的全部觀察對象按是否服用某藥劃分為兩個數(shù)據(jù)集,再依次從用藥組選出一個個體,并從對照組找出和該個體的PS值最為接近的全部個體(小于所規(guī)定的選擇標準),再隨機從這些選定的對象中抽取一個作為對照。就這樣直至符合選擇標準的觀察對象全部抽取。

    第三,按照抽取好的樣本再次通過統(tǒng)計算法(如ROR、PRR)等分析某藥和某不良反應(yīng)的關(guān)聯(lián)〔16〕。具體步驟如圖1所示。匹配方法通過SAS9.1.3編程實現(xiàn),本文參考 Parsons的貪婪法進行1∶1匹配〔17〕。

    圖1 應(yīng)用PS均衡SRS數(shù)據(jù)的具體步驟

    2.實例應(yīng)用

    實際數(shù)據(jù)來自于上海市2009年藥品不良反應(yīng)監(jiān)測中心SRS的報告,共24 297例。數(shù)據(jù)集的初步清理包括,根據(jù)國家衛(wèi)生部發(fā)布的常用處方藥品通用名錄規(guī)范和編碼藥品的通用名稱,并采用MedDRA14.0中文版的優(yōu)先級(PT)規(guī)范了ADR名稱。應(yīng)用PS均衡混雜協(xié)變量的方法與模擬試驗類似,此處不再贅述。

    3.統(tǒng)計軟件

    應(yīng)用SAS9.1.3軟件實現(xiàn)SRS數(shù)據(jù)的模擬、PS匹配、logistic回歸模型的建立、χ2檢驗和t檢驗等。

    結(jié) 果

    1.模擬試驗結(jié)果

    模擬試驗產(chǎn)生1000例數(shù)據(jù),其中結(jié)果變量模型參數(shù):β0ADR= -12.458,β1ADR=1.125,β2ADR=2.979,β3ADR=6.362,β4ADR=6.467。Drug1 的 OR=193.268,Drug2的 OR=216.613;Drug3、Drug4的 OR值約為1,沒有進入模型。具體結(jié)果如表1所示,可以看出Drug1-ADR、Drug2-ADR組合為真陽性信號,PS匹配前后的信號值變化不大,而假陽性的 Drug3-ADR、Drug4-ADR在PS前的信號值較強,但在考慮了基線均衡后便無信號產(chǎn)生。

    表1 模擬的四個Drug-ADR組合在PS均衡協(xié)變量前后挖掘算法值的變化

    2.實際應(yīng)用結(jié)果應(yīng)用數(shù)據(jù)挖掘算法對09年上海市的SRS數(shù)據(jù)分析,將初步得到的可疑組合按信號值大小排列,其中部分組合如表2所示。

    表2 四種數(shù)據(jù)挖掘算法對上海市09年SRS初步挖掘結(jié)果舉例

    在表2中列出的可疑組合除“喹硫平-閉經(jīng)”外,其他均有說明書及文獻支持,為加強此組合的可信性,我們應(yīng)用SRS現(xiàn)存的基線協(xié)變量信息,通過傾向性評分法來對其進一步分析。通過專業(yè)知識和變量選擇策略〔15〕,我們從喹硫平組和非喹硫平組(對照組)的原始數(shù)據(jù)中選擇了年齡、體重、性別、ADR發(fā)生時間和是否有合并用藥等協(xié)變量納入模型。對年齡以國際通用的標準(WONCA)〔18〕劃分成八個等級:<1 歲,1 ~4 歲,5~14歲,15~24歲,25~44歲,45~64歲,65~74歲,>75歲。ADR發(fā)生時間按照季節(jié)性劃分為四段。此外,根據(jù) Rubin的理論,當 C-statistics在 0.65到0.85時為判定模型較好的標準〔19〕。該模型的C-statistics=0.750,因此擬合效果較好,待估算PS之后進行匹配。表3為比較對PS值進行匹配前后協(xié)變量在組間的分布情況。如表3所示,除性別外其他四個協(xié)變量在PS匹配前兩組中的分布都是不均衡的P<0.05,差異均有統(tǒng)計學意義;當進行PS匹配后,協(xié)變量在兩組間的分布達到均衡,P>0.1。再將匹配后的數(shù)據(jù)用三種挖掘算法分析喹硫平與閉經(jīng)的關(guān)系。表4則展示了三種信號挖掘算法在PS匹配之后,其可疑信號值均顯著降低,信號消失。

    表3 PS匹配前后喹硫平組與其他藥品組間協(xié)變量分布的均衡性

    表4 三種數(shù)據(jù)挖掘方法在PS匹配前后計算“喹硫平-閉經(jīng)”組合的信號值

    討 論

    數(shù)據(jù)挖掘方法提供的自動信號產(chǎn)生過程對藥品安全性的監(jiān)測起到了關(guān)鍵的作用。但由于SRS數(shù)據(jù)量巨大和藥品上市后的一些限制,導致這些方法因不能有效的體現(xiàn)隱藏在數(shù)據(jù)中的所有信息而產(chǎn)生諸多假陽性。傾向性評分法的優(yōu)點在于〔17〕:首先,它能將多個協(xié)變量綜合為一個值來分析,既能避免了過度分層和過分匹配等問題,同樣也避免了自變量間的共線性;其次能使估計PS前后的不均衡程度更明顯,因為當我們應(yīng)用傳統(tǒng)的回歸模型矯正混雜時,這種不均衡程度經(jīng)常被掩蓋。再者,PS簡化了多重結(jié)果需要考慮的事項,因為傾向性評分僅僅估計一次就可以分別的由每個結(jié)果所應(yīng)用。此外,用標準的統(tǒng)計軟件來實施會相對簡單易懂。本研究在進行混雜協(xié)變量矯正之前,三種數(shù)據(jù)挖掘方法均對“喹硫平-閉經(jīng)”組合提示為可疑信號,其協(xié)變量除了性別外,年齡、體重、是否合并用藥和不良反應(yīng)發(fā)生時間在兩組間分布都不均衡,待矯正之后便無信號產(chǎn)生。經(jīng)分析原始數(shù)據(jù)及詳細查閱咨詢后得知,由于喹硫平組中有較多合并用藥,其多數(shù)為奮乃靜等傳統(tǒng)的抗精神病類藥品,他們一定程度上會引起高催乳素血癥〔20〕,繼而可能出現(xiàn)閉經(jīng)。而喹硫平乃新一代抗精神病類藥品,許多實驗證明其導致高催乳素血癥及閉經(jīng)的可能性很小。因此,不能忽略合并用藥因素來單獨分析喹硫平與閉經(jīng)反應(yīng)之間的關(guān)聯(lián)。

    然而一些不足也難以避免:一方面,如同傳統(tǒng)的矯正方法不能去除非觀察性的偏倚,此方法也不能去除由這種混雜導致的偏倚;另一方面,如果協(xié)變量選擇不恰當會降低效能,例如只選擇與處理因素有關(guān)而與結(jié)果無關(guān)的協(xié)變量有時會降低估計的效能,甚至會增加偏倚,但如果樣本量足夠大,即使協(xié)變量的選擇不當也不會受到較大影響〔17,19,21〕。

    PS應(yīng)用于SRS中還存在一定的阻礙。由于SRS屬于觀察性研究的數(shù)據(jù),不能保證每種藥品的用藥人數(shù)都足夠,樣本量太少會導致PS效能降低甚至無法分析;其次是自發(fā)上報的數(shù)據(jù)都會存在漏報問題,如果納入模型的協(xié)變量有缺失數(shù)據(jù),勢必會影響模型擬合效果,此時必須將缺失數(shù)據(jù)進行填補后再納入分析。目前,SRS中的現(xiàn)有信息并沒有完全的應(yīng)用于藥品安全性研究中,在不改變挖掘算法截斷值的條件下要減少假陽性信號的產(chǎn)生,就必須添加更多的混雜因素進行分析〔6〕。若能將現(xiàn)有的信息全部運用使基線均衡,再得出的信號值其可信性就會顯著提高,進而為專家評價其信號關(guān)聯(lián)提供更多有力可靠的證據(jù)。相信隨著PS的不斷改進和SRS數(shù)據(jù)質(zhì)量的逐步提高,傾向性評分法將會成為均衡SRS數(shù)據(jù)的有力工具。

    1.Hopstadius J,Niklas NG,Bate A,et al.Impact of stratification on adverse drug reaction surveillance drug safety,2008,31(11):1035-1048.

    2.van Puijenbreek EP,Bate A,Leufkens HG,et al.A comparison of measures of disproportionality for signal detection in spontaneous reporting systems for adverse drug reactions.Pharmacoepidemiology and drug safety,2002,11:3-10.

    3.Bate A,Lindquist M,Edwards IR,et a1.A Bayesian neural network for adverse drug reaction signal generation European Journal of Clinical Pharmacology,1998,54:315-321.

    4.DuMouchel W,Pregibon D.Empirical bayes screening for multi-item association.Conference on knowledge discovery in data,proceedings of the seventh ACM SIGKDD international conference on knowledge discovery and data mining.

    5.Hauben M,Madigan D,Gerrits CM,et al.The role of data mining in pharmacovigilance.Drug safety,2005,4(5):929-948.

    6.Stephen JW.Stratification for spontaneous report database.Drug safety,2008,31(11):1049-1052.

    7.Rosenbaum PR,Rubin DB.The central role of the propensity score in observational studies for causal effects.Biometrika,1983,70:41-55.

    8.D'Agostino RB.Propensity score methods for bias reduction in the comparison of a treatment to a non-randomized control group.Stat Med,1998,17:2265-2281.

    9.Sekula Dipl Math P,Caputo A,Dunant A,et al.An application of propensity score method to estimate the treatment effect of corticosteroids in patients with severe cutaneous adverse reactions.Pharmacoepidemiology and drug safety,2010,19:10-18.

    10.Perkins SM,Tu W,Underhill MG,et al.The use of propensity scores in pharmacoepidemiology research.Pharmacoepidemiology and drug safety,2000,9:93-101.

    11.Rosenbaum PR,Rubin DB.Reducing bias in observational studies using subclassification on the propensity score.American Statistical Association,1984,79(387):516-524.

    12.Cochran WG.The effectiveness of adjustment by subclassification in removing bias in observational studies.Biometrics,1968,24:295-313.

    13.Austin PC.A critical appraisal of propensity score matching in the medical literature between 1996 and 2003.Stat Med,2008,27:2037-2049.

    14.Austin PC.Propensity-score matching in the cardiovascular surgery literature from 2004 to 2006:a systematic review and suggestions for improvement.Thorac Cardiovasc Surg,2007,134:1128-1135.

    15.Brookhart MA,Schneeweiss S,Rothman KJ,et al.Variable selection for propensity score models.Am J Epidemiol,2006,163:1149-1156.

    16.Rubin DB,Thomas N.Matching using estimated propensity scores:relating theory to practice.Biometrics,1996,52:249-264.

    17.Parsons LS,Ovation Research Group,Seattle WA.Reducing bias in a propensity score match-pair sample using greedy matching technique.http://www2.sas.com/proceedings/sugi26/p214-26.pdf.2010-9-22.

    18.Bentzen N.An international glossary for general/family practice.Fam Pract,1995,12(3):341-369.

    19.Rubin DB.Estimating causal effects from large data sets using propensity scores.Ann Intern Med,1997,127:757-763.

    20.吳莉珍,安彬.抗精神病藥物對女性精神病患者血清催乳素的影響.臨床心身疾病雜志.2007,13(5):247.

    21.Wang J,Donnan PT.Propensity score methods in drug safety studies:practice,strengths and limitations.Pharmacoepidemiology and drug safety,2001,10:341-344.

    猜你喜歡
    傾向性喹硫平數(shù)據(jù)挖掘
    小劑量喹硫平對文拉法辛治療抑郁癥增效作用的比較研究
    基于模糊數(shù)學法的阿舍勒銅礦深部巖體巖爆傾向性預(yù)測
    探討人工智能與數(shù)據(jù)挖掘發(fā)展趨勢
    阿立哌唑與喹硫平用于晚發(fā)型精神分裂癥的療效對比
    基于并行計算的大數(shù)據(jù)挖掘在電網(wǎng)中的應(yīng)用
    電力與能源(2017年6期)2017-05-14 06:19:37
    探討舍曲林合并喹硫平治療強迫癥的療效和安全性
    一種基于Hadoop的大數(shù)據(jù)挖掘云服務(wù)及應(yīng)用
    關(guān)于醫(yī)患沖突報道的傾向性分析——以“湘潭產(chǎn)婦死亡案”為例
    “沒準兒”“不一定”“不見得”和“說不定”的語義傾向性和主觀性差異
    語言與翻譯(2015年4期)2015-07-18 11:07:43
    喹硫平治療酒精致焦慮、抑郁的臨床效果觀察
    欧美黑人欧美精品刺激| 美女午夜性视频免费| 欧美黑人精品巨大| 成人三级做爰电影| 视频区欧美日本亚洲| 欧美亚洲日本最大视频资源| 亚洲专区字幕在线| 黄色视频在线播放观看不卡| 亚洲国产欧美网| 精品少妇黑人巨大在线播放| 精品亚洲成a人片在线观看| 免费黄频网站在线观看国产| 国产又色又爽无遮挡免费看| 亚洲精品国产精品久久久不卡| 久久久精品区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 中文字幕最新亚洲高清| 亚洲色图 男人天堂 中文字幕| 久久人妻熟女aⅴ| 美女午夜性视频免费| 黑人巨大精品欧美一区二区mp4| 最近最新中文字幕大全电影3 | 亚洲 欧美一区二区三区| 亚洲欧美一区二区三区久久| 午夜福利影视在线免费观看| 亚洲中文字幕日韩| 黄网站色视频无遮挡免费观看| 97在线人人人人妻| h视频一区二区三区| 亚洲熟妇熟女久久| 亚洲欧美一区二区三区黑人| 亚洲天堂av无毛| 久久精品国产99精品国产亚洲性色 | 国产精品久久久久成人av| 国产精品1区2区在线观看. | 日韩免费高清中文字幕av| 建设人人有责人人尽责人人享有的| 亚洲精品国产区一区二| av在线播放免费不卡| 女人久久www免费人成看片| 亚洲精品美女久久av网站| 老汉色av国产亚洲站长工具| 国产成人精品久久二区二区91| 啪啪无遮挡十八禁网站| 免费在线观看黄色视频的| 国产成人啪精品午夜网站| 啦啦啦在线免费观看视频4| 国产精品久久久久久精品电影小说| 久久婷婷成人综合色麻豆| 国产淫语在线视频| 国产午夜精品久久久久久| 精品国产一区二区三区久久久樱花| 一边摸一边抽搐一进一小说 | 女人久久www免费人成看片| 亚洲一码二码三码区别大吗| 视频区图区小说| 精品久久蜜臀av无| 久久99热这里只频精品6学生| 国产一区二区 视频在线| 亚洲国产欧美网| 亚洲五月色婷婷综合| 久久久久久人人人人人| 国产免费av片在线观看野外av| 啦啦啦 在线观看视频| 精品乱码久久久久久99久播| 久久久国产精品麻豆| 99国产综合亚洲精品| 亚洲av欧美aⅴ国产| 免费观看a级毛片全部| 高清毛片免费观看视频网站 | 九色亚洲精品在线播放| 婷婷成人精品国产| 亚洲第一av免费看| 视频区欧美日本亚洲| 在线观看免费日韩欧美大片| 亚洲国产av新网站| 欧美日韩亚洲国产一区二区在线观看 | 日韩欧美国产一区二区入口| 国产精品电影一区二区三区 | 欧美黄色片欧美黄色片| 天堂动漫精品| 久久精品人人爽人人爽视色| 国产一区二区三区在线臀色熟女 | 国产高清videossex| 成人国语在线视频| 久久精品成人免费网站| 亚洲国产欧美一区二区综合| 在线观看免费日韩欧美大片| 成人国产av品久久久| 国产麻豆69| 热99re8久久精品国产| 国产精品免费大片| 久久精品国产综合久久久| 999精品在线视频| 午夜老司机福利片| e午夜精品久久久久久久| 日日摸夜夜添夜夜添小说| 丰满饥渴人妻一区二区三| 国产又色又爽无遮挡免费看| h视频一区二区三区| 男女床上黄色一级片免费看| 大型av网站在线播放| 一本色道久久久久久精品综合| 丰满迷人的少妇在线观看| 热99国产精品久久久久久7| 在线观看www视频免费| 18禁观看日本| 亚洲成人国产一区在线观看| 老鸭窝网址在线观看| 老司机福利观看| 啦啦啦中文免费视频观看日本| 欧美日韩av久久| 好男人电影高清在线观看| 91精品三级在线观看| 又大又爽又粗| 大码成人一级视频| 搡老熟女国产l中国老女人| 自线自在国产av| 亚洲男人天堂网一区| a在线观看视频网站| 国产精品免费大片| 亚洲午夜精品一区,二区,三区| 18禁观看日本| 十八禁网站免费在线| 一夜夜www| 日本一区二区免费在线视频| 欧美国产精品一级二级三级| 女性被躁到高潮视频| 999久久久国产精品视频| www.自偷自拍.com| 亚洲专区中文字幕在线| 伊人久久大香线蕉亚洲五| 久久婷婷成人综合色麻豆| 欧美午夜高清在线| 免费观看av网站的网址| 制服诱惑二区| a级片在线免费高清观看视频| 久久久久久人人人人人| 少妇精品久久久久久久| 久久久精品区二区三区| 国产老妇伦熟女老妇高清| 中文字幕av电影在线播放| 国产精品亚洲av一区麻豆| 欧美激情极品国产一区二区三区| 日韩免费av在线播放| 欧美日韩av久久| 久久精品亚洲精品国产色婷小说| 国产在线精品亚洲第一网站| 777米奇影视久久| 免费观看a级毛片全部| av国产精品久久久久影院| 久久久久国产一级毛片高清牌| 国内毛片毛片毛片毛片毛片| 脱女人内裤的视频| 午夜久久久在线观看| 亚洲国产中文字幕在线视频| 黑人猛操日本美女一级片| 黄色丝袜av网址大全| 国产99久久九九免费精品| 手机成人av网站| 欧美亚洲日本最大视频资源| av欧美777| 成在线人永久免费视频| 日韩欧美一区视频在线观看| 欧美精品高潮呻吟av久久| 国产一区二区三区综合在线观看| 电影成人av| 99久久精品国产亚洲精品| 大片电影免费在线观看免费| 亚洲色图 男人天堂 中文字幕| 欧美精品一区二区免费开放| 黄片大片在线免费观看| 精品国产一区二区久久| 777米奇影视久久| 亚洲精品国产色婷婷电影| 中文字幕最新亚洲高清| 色视频在线一区二区三区| 制服人妻中文乱码| 搡老乐熟女国产| 国产色视频综合| 欧美精品一区二区大全| 视频在线观看一区二区三区| 亚洲国产av新网站| 欧美国产精品一级二级三级| 脱女人内裤的视频| 又黄又粗又硬又大视频| 在线观看免费高清a一片| 中文字幕精品免费在线观看视频| 国产成人欧美| 欧美 日韩 精品 国产| 女人被躁到高潮嗷嗷叫费观| 日韩中文字幕欧美一区二区| 亚洲精品国产色婷婷电影| 午夜老司机福利片| 国产精品秋霞免费鲁丝片| 亚洲一卡2卡3卡4卡5卡精品中文| 91麻豆av在线| 亚洲精品国产一区二区精华液| 亚洲精品国产区一区二| 自线自在国产av| 国产日韩欧美在线精品| 久久精品国产亚洲av香蕉五月 | 乱人伦中国视频| 欧美日韩亚洲高清精品| 欧美成人午夜精品| 色视频在线一区二区三区| 在线十欧美十亚洲十日本专区| 国产精品亚洲一级av第二区| 亚洲成人手机| 激情在线观看视频在线高清 | 亚洲欧美日韩另类电影网站| 国产又色又爽无遮挡免费看| 男女免费视频国产| 桃花免费在线播放| 中文字幕人妻熟女乱码| 国产精品久久久久久人妻精品电影 | 国产精品二区激情视频| 黄频高清免费视频| 免费观看a级毛片全部| 久久久精品94久久精品| 日韩免费高清中文字幕av| 女人被躁到高潮嗷嗷叫费观| 一区福利在线观看| 久久中文字幕人妻熟女| 69精品国产乱码久久久| 国产精品欧美亚洲77777| 亚洲国产欧美网| 悠悠久久av| 国产精品久久久av美女十八| 一边摸一边抽搐一进一出视频| 亚洲国产成人一精品久久久| 久久天堂一区二区三区四区| 日本精品一区二区三区蜜桃| 黄频高清免费视频| 精品一区二区三区四区五区乱码| 黄色a级毛片大全视频| 精品午夜福利视频在线观看一区 | 亚洲国产欧美一区二区综合| 免费观看a级毛片全部| 最新在线观看一区二区三区| 精品国产乱子伦一区二区三区| 成年人午夜在线观看视频| 亚洲av成人不卡在线观看播放网| 久久中文看片网| 亚洲av日韩在线播放| av片东京热男人的天堂| 日本av免费视频播放| 一边摸一边抽搐一进一小说 | 久久国产精品大桥未久av| 亚洲国产欧美一区二区综合| 人人妻人人澡人人看| 最近最新中文字幕大全电影3 | 亚洲精品国产一区二区精华液| 亚洲人成电影免费在线| 成人18禁在线播放| 精品卡一卡二卡四卡免费| 婷婷成人精品国产| kizo精华| 国产成人精品无人区| 国产单亲对白刺激| 中文字幕精品免费在线观看视频| 国产深夜福利视频在线观看| 亚洲久久久国产精品| 最新的欧美精品一区二区| 黄片大片在线免费观看| 亚洲av电影在线进入| 俄罗斯特黄特色一大片| 久9热在线精品视频| 成人18禁高潮啪啪吃奶动态图| 九色亚洲精品在线播放| 一区二区三区激情视频| 黄片播放在线免费| 老司机深夜福利视频在线观看| 国产精品免费视频内射| 新久久久久国产一级毛片| 国产午夜精品久久久久久| 久久免费观看电影| 欧美一级毛片孕妇| 在线观看舔阴道视频| 美女视频免费永久观看网站| 久久精品熟女亚洲av麻豆精品| 日韩视频一区二区在线观看| 成在线人永久免费视频| 国产精品免费视频内射| 国产精品久久久av美女十八| 亚洲精品中文字幕一二三四区 | 色94色欧美一区二区| 69精品国产乱码久久久| 久久久久视频综合| 国产精品久久久久成人av| 免费看a级黄色片| 国产亚洲欧美在线一区二区| 黄色片一级片一级黄色片| 色视频在线一区二区三区| 麻豆乱淫一区二区| 午夜福利乱码中文字幕| 99国产精品一区二区蜜桃av | 黑人欧美特级aaaaaa片| 肉色欧美久久久久久久蜜桃| 久久九九热精品免费| 国产精品久久电影中文字幕 | 欧美精品人与动牲交sv欧美| 成年女人毛片免费观看观看9 | 国产亚洲午夜精品一区二区久久| 亚洲精品久久午夜乱码| 一本综合久久免费| 99国产精品一区二区三区| 日韩成人在线观看一区二区三区| 国产一区二区三区在线臀色熟女 | 国产精品98久久久久久宅男小说| 自线自在国产av| 久久国产精品影院| 另类亚洲欧美激情| 又大又爽又粗| 丰满人妻熟妇乱又伦精品不卡| 黄网站色视频无遮挡免费观看| 免费黄频网站在线观看国产| 妹子高潮喷水视频| 1024香蕉在线观看| 国产黄频视频在线观看| 午夜免费鲁丝| 国产无遮挡羞羞视频在线观看| 亚洲精品美女久久久久99蜜臀| 嫁个100分男人电影在线观看| 最新美女视频免费是黄的| 欧美日韩av久久| 涩涩av久久男人的天堂| 久久精品国产亚洲av高清一级| 亚洲欧洲精品一区二区精品久久久| 黄色 视频免费看| 国产精品1区2区在线观看. | 电影成人av| 国产精品 欧美亚洲| 久久国产精品影院| 天堂动漫精品| 久久这里只有精品19| netflix在线观看网站| 国产在线一区二区三区精| 日韩欧美免费精品| 亚洲第一欧美日韩一区二区三区 | 老司机靠b影院| 国产欧美日韩精品亚洲av| 国产极品粉嫩免费观看在线| 免费在线观看日本一区| av网站在线播放免费| 亚洲欧美一区二区三区久久| 欧美精品亚洲一区二区| 老司机亚洲免费影院| 国产精品影院久久| 女性生殖器流出的白浆| 别揉我奶头~嗯~啊~动态视频| 国产精品一区二区免费欧美| 欧美日韩黄片免| 一本一本久久a久久精品综合妖精| 如日韩欧美国产精品一区二区三区| 亚洲 欧美一区二区三区| 18禁观看日本| 香蕉丝袜av| 久久ye,这里只有精品| 日本五十路高清| 一本一本久久a久久精品综合妖精| 国产91精品成人一区二区三区 | 亚洲精品美女久久av网站| 麻豆av在线久日| 国产单亲对白刺激| 国产精品久久久久成人av| 亚洲第一av免费看| 久久久久久久精品吃奶| 久久影院123| 日韩免费av在线播放| 精品免费久久久久久久清纯 | 国产成人一区二区三区免费视频网站| 欧美人与性动交α欧美软件| 中文字幕制服av| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区激情短视频| 国产精品秋霞免费鲁丝片| 在线 av 中文字幕| 他把我摸到了高潮在线观看 | 久久中文字幕一级| 欧美日韩一级在线毛片| 国产91精品成人一区二区三区 | 精品欧美一区二区三区在线| 国产精品国产高清国产av | 精品人妻熟女毛片av久久网站| 久热这里只有精品99| 视频在线观看一区二区三区| 国产1区2区3区精品| 天天影视国产精品| 久久久久久久久久久久大奶| 欧美日韩亚洲综合一区二区三区_| 色综合欧美亚洲国产小说| 日本黄色视频三级网站网址 | 91精品国产国语对白视频| 国产又爽黄色视频| 欧美日韩一级在线毛片| 一边摸一边抽搐一进一出视频| 19禁男女啪啪无遮挡网站| 人妻久久中文字幕网| 老熟女久久久| 两性夫妻黄色片| 久久久久网色| 免费在线观看完整版高清| 亚洲精品乱久久久久久| 另类精品久久| 啪啪无遮挡十八禁网站| 亚洲成人手机| kizo精华| 精品欧美一区二区三区在线| 91字幕亚洲| 久久久久国内视频| 午夜福利视频精品| 人人妻,人人澡人人爽秒播| 肉色欧美久久久久久久蜜桃| 免费av中文字幕在线| 午夜免费成人在线视频| 国产成人影院久久av| 久久久久久久久免费视频了| 亚洲avbb在线观看| 动漫黄色视频在线观看| 国内毛片毛片毛片毛片毛片| 丰满迷人的少妇在线观看| 婷婷成人精品国产| 亚洲国产精品一区二区三区在线| 高潮久久久久久久久久久不卡| 韩国精品一区二区三区| 电影成人av| 午夜福利视频精品| 国产又色又爽无遮挡免费看| 欧美日韩视频精品一区| 中文字幕精品免费在线观看视频| 中文字幕人妻丝袜一区二区| 中文字幕人妻熟女乱码| 99精品久久久久人妻精品| 在线十欧美十亚洲十日本专区| 免费看a级黄色片| 日韩成人在线观看一区二区三区| 亚洲美女黄片视频| 亚洲综合色网址| 国产欧美日韩一区二区三区在线| 51午夜福利影视在线观看| 青青草视频在线视频观看| 精品一品国产午夜福利视频| 99久久精品国产亚洲精品| 999久久久精品免费观看国产| 狠狠婷婷综合久久久久久88av| 人人澡人人妻人| 久久久精品免费免费高清| av片东京热男人的天堂| 大片免费播放器 马上看| 一二三四在线观看免费中文在| 在线观看免费视频网站a站| tocl精华| 超碰成人久久| 99精国产麻豆久久婷婷| 国产精品亚洲av一区麻豆| 一级毛片女人18水好多| 久久热在线av| 亚洲五月色婷婷综合| 黑丝袜美女国产一区| 丰满迷人的少妇在线观看| 日本欧美视频一区| 午夜精品久久久久久毛片777| 曰老女人黄片| 人人妻人人澡人人爽人人夜夜| 久久久久久久久免费视频了| 午夜福利视频在线观看免费| 麻豆国产av国片精品| 色婷婷久久久亚洲欧美| 久久国产精品大桥未久av| 啦啦啦中文免费视频观看日本| a级片在线免费高清观看视频| 欧美成人免费av一区二区三区 | 欧美激情久久久久久爽电影 | 久久久水蜜桃国产精品网| 久久精品国产99精品国产亚洲性色 | 久久久久久久久久久久大奶| 午夜福利视频精品| 日本精品一区二区三区蜜桃| 午夜福利一区二区在线看| 男女之事视频高清在线观看| 成人手机av| 在线十欧美十亚洲十日本专区| 性高湖久久久久久久久免费观看| 色综合欧美亚洲国产小说| 自线自在国产av| 丝袜人妻中文字幕| 国产成+人综合+亚洲专区| 岛国毛片在线播放| 在线播放国产精品三级| 丁香六月天网| 十八禁人妻一区二区| 曰老女人黄片| 俄罗斯特黄特色一大片| 国产免费视频播放在线视频| 久久人妻福利社区极品人妻图片| 久久国产精品男人的天堂亚洲| 大型av网站在线播放| 在线观看舔阴道视频| 女人久久www免费人成看片| 久久人人爽av亚洲精品天堂| 国产主播在线观看一区二区| 18禁裸乳无遮挡动漫免费视频| 亚洲国产成人一精品久久久| 成人手机av| 国产精品免费视频内射| 国产免费福利视频在线观看| 久久影院123| 精品人妻在线不人妻| 婷婷丁香在线五月| 精品人妻1区二区| 久9热在线精品视频| 久久久久久免费高清国产稀缺| 欧美精品人与动牲交sv欧美| 天天躁日日躁夜夜躁夜夜| 99国产精品免费福利视频| 亚洲精品国产色婷婷电影| 国产精品香港三级国产av潘金莲| 在线观看人妻少妇| www日本在线高清视频| 少妇粗大呻吟视频| 黄频高清免费视频| 久久国产精品人妻蜜桃| 国产亚洲精品第一综合不卡| 夜夜爽天天搞| 啦啦啦中文免费视频观看日本| 日韩 欧美 亚洲 中文字幕| 久9热在线精品视频| 亚洲va日本ⅴa欧美va伊人久久| 国产精品国产av在线观看| e午夜精品久久久久久久| 国产亚洲av高清不卡| videosex国产| 亚洲国产中文字幕在线视频| 飞空精品影院首页| 午夜福利在线免费观看网站| av一本久久久久| 自线自在国产av| 69av精品久久久久久 | 欧美黄色片欧美黄色片| 久久亚洲真实| 久久人妻av系列| 免费看a级黄色片| 国产成人精品久久二区二区91| 成人影院久久| 亚洲国产成人一精品久久久| 叶爱在线成人免费视频播放| 国产高清国产精品国产三级| 国产欧美日韩一区二区三| 午夜成年电影在线免费观看| 午夜福利视频在线观看免费| aaaaa片日本免费| 男人舔女人的私密视频| 欧美日韩一级在线毛片| 午夜视频精品福利| 国产国语露脸激情在线看| 欧美+亚洲+日韩+国产| 天天躁狠狠躁夜夜躁狠狠躁| av福利片在线| 999久久久国产精品视频| 日本精品一区二区三区蜜桃| 欧美日韩一级在线毛片| 精品国产亚洲在线| 大香蕉久久成人网| 中文亚洲av片在线观看爽 | 男女无遮挡免费网站观看| 欧美黄色淫秽网站| 精品亚洲成国产av| 国产欧美日韩综合在线一区二区| 蜜桃国产av成人99| 亚洲色图 男人天堂 中文字幕| 丰满饥渴人妻一区二区三| 欧美乱妇无乱码| 国产无遮挡羞羞视频在线观看| 久9热在线精品视频| 三上悠亚av全集在线观看| 十八禁网站网址无遮挡| 国产精品秋霞免费鲁丝片| www.自偷自拍.com| 亚洲第一欧美日韩一区二区三区 | 一区二区日韩欧美中文字幕| 制服人妻中文乱码| 欧美黑人欧美精品刺激| 最近最新中文字幕大全电影3 | 狠狠婷婷综合久久久久久88av| 热99国产精品久久久久久7| 亚洲人成77777在线视频| 青青草视频在线视频观看| 80岁老熟妇乱子伦牲交| 18禁裸乳无遮挡动漫免费视频| 日韩免费av在线播放| 久久性视频一级片| 精品一区二区三区视频在线观看免费 | 99re6热这里在线精品视频| 成人精品一区二区免费| 中文字幕人妻丝袜制服| 777久久人妻少妇嫩草av网站| 免费人妻精品一区二区三区视频| 动漫黄色视频在线观看| 中文字幕人妻熟女乱码| 伊人久久大香线蕉亚洲五| 老司机午夜福利在线观看视频 | 免费看十八禁软件| 亚洲精品国产区一区二| 菩萨蛮人人尽说江南好唐韦庄| 日韩视频在线欧美| 欧美大码av| 国产免费现黄频在线看| 老司机亚洲免费影院| 18在线观看网站| 亚洲国产欧美日韩在线播放| 国产av又大| 亚洲av成人不卡在线观看播放网| 色综合欧美亚洲国产小说| 国产高清激情床上av|