• <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
    喹硫平治療酒精致焦慮、抑郁的臨床效果觀察
    99视频精品全部免费 在线| 伊人久久国产一区二区| 五月天丁香电影| 天堂8中文在线网| 丝袜美足系列| 狂野欧美白嫩少妇大欣赏| 卡戴珊不雅视频在线播放| 国产精品人妻久久久影院| 汤姆久久久久久久影院中文字幕| 久久99蜜桃精品久久| 伊人亚洲综合成人网| 日产精品乱码卡一卡2卡三| 男女边摸边吃奶| 久久99蜜桃精品久久| 成人手机av| 国产淫语在线视频| 日韩制服骚丝袜av| 秋霞伦理黄片| 大又大粗又爽又黄少妇毛片口| 在线观看免费视频网站a站| 人人妻人人爽人人添夜夜欢视频| 国产精品一区二区在线观看99| av有码第一页| 国产黄频视频在线观看| 人人妻人人添人人爽欧美一区卜| 免费看av在线观看网站| 大香蕉久久网| 日韩制服骚丝袜av| 2021少妇久久久久久久久久久| 新久久久久国产一级毛片| 五月开心婷婷网| 欧美人与性动交α欧美精品济南到 | 久久久久久久久久久免费av| 国产欧美日韩综合在线一区二区| 久久精品久久精品一区二区三区| 在线观看www视频免费| 欧美老熟妇乱子伦牲交| 婷婷色综合大香蕉| 午夜福利视频精品| 在线免费观看不下载黄p国产| 亚洲av福利一区| 最近的中文字幕免费完整| 在线观看国产h片| 久久久久久久久大av| 一本色道久久久久久精品综合| 伊人亚洲综合成人网| 成人免费观看视频高清| 国产探花极品一区二区| 99精国产麻豆久久婷婷| 18在线观看网站| 少妇被粗大猛烈的视频| 伊人亚洲综合成人网| 一边摸一边做爽爽视频免费| 国产免费视频播放在线视频| 国产精品久久久久久av不卡| 久久久久网色| 少妇人妻精品综合一区二区| 亚洲欧美一区二区三区国产| 日本欧美国产在线视频| 国产 精品1| 日本黄色片子视频| 91成人精品电影| 亚洲精华国产精华液的使用体验| 亚洲人成77777在线视频| 视频在线观看一区二区三区| 草草在线视频免费看| 欧美国产精品一级二级三级| 亚洲成人手机| 人妻系列 视频| 老女人水多毛片| 久久久久久久久久成人| 亚洲国产精品一区二区三区在线| 欧美人与性动交α欧美精品济南到 | 亚洲国产欧美在线一区| 最近中文字幕高清免费大全6| 老司机亚洲免费影院| 欧美日韩亚洲高清精品| 久久精品熟女亚洲av麻豆精品| 女人久久www免费人成看片| 国产精品 国内视频| 亚洲国产精品专区欧美| 丰满迷人的少妇在线观看| av国产久精品久网站免费入址| 99热全是精品| 又黄又爽又刺激的免费视频.| 一个人看视频在线观看www免费| 三级国产精品欧美在线观看| 七月丁香在线播放| 中文精品一卡2卡3卡4更新| 黑丝袜美女国产一区| 亚洲四区av| 国产精品人妻久久久影院| 久久久久久久久久人人人人人人| 综合色丁香网| 欧美少妇被猛烈插入视频| 久久青草综合色| 老司机亚洲免费影院| 国产精品蜜桃在线观看| 亚洲国产欧美日韩在线播放| 另类精品久久| 大香蕉97超碰在线| 国产成人精品婷婷| 午夜日本视频在线| 精品酒店卫生间| 十八禁网站网址无遮挡| 久久精品熟女亚洲av麻豆精品| 国产无遮挡羞羞视频在线观看| 亚洲人成网站在线播| 99久久精品一区二区三区| 老司机影院毛片| 五月伊人婷婷丁香| 狂野欧美激情性xxxx在线观看| 国精品久久久久久国模美| 久久精品国产亚洲av天美| 高清午夜精品一区二区三区| 久久午夜综合久久蜜桃| 精品亚洲乱码少妇综合久久| 欧美xxⅹ黑人| 天堂8中文在线网| 国产欧美日韩一区二区三区在线 | 午夜免费鲁丝| 欧美 日韩 精品 国产| 亚洲经典国产精华液单| 有码 亚洲区| 人妻系列 视频| 人妻制服诱惑在线中文字幕| 中文精品一卡2卡3卡4更新| 久久久久国产精品人妻一区二区| 国产成人精品福利久久| 一本色道久久久久久精品综合| 国产精品久久久久久精品电影小说| 国产成人午夜福利电影在线观看| 亚洲怡红院男人天堂| 欧美 亚洲 国产 日韩一| 日韩一区二区视频免费看| 日产精品乱码卡一卡2卡三| 久久青草综合色| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品aⅴ在线观看| 国产亚洲欧美精品永久| 亚洲精品成人av观看孕妇| 精品人妻熟女毛片av久久网站| 国产一区二区在线观看av| 美女cb高潮喷水在线观看| 久久这里有精品视频免费| 国产午夜精品久久久久久一区二区三区| 91久久精品国产一区二区成人| 精品亚洲乱码少妇综合久久| 亚洲国产精品一区二区三区在线| 成人漫画全彩无遮挡| 欧美日韩亚洲高清精品| 男的添女的下面高潮视频| 中国国产av一级| 久久久久久久国产电影| 午夜福利在线观看免费完整高清在| av一本久久久久| 韩国高清视频一区二区三区| 新久久久久国产一级毛片| 国产片特级美女逼逼视频| 22中文网久久字幕| 这个男人来自地球电影免费观看 | 亚洲精品,欧美精品| 国产亚洲欧美精品永久| 丰满少妇做爰视频| 免费黄网站久久成人精品| 91精品三级在线观看| 亚洲第一区二区三区不卡| 国产亚洲最大av| 在线亚洲精品国产二区图片欧美 | 少妇人妻 视频| 久久久精品免费免费高清| 人人妻人人爽人人添夜夜欢视频| 老司机影院毛片| 久久久久久久久久久免费av| 日韩av在线免费看完整版不卡| 七月丁香在线播放| 黄色视频在线播放观看不卡| 日韩中文字幕视频在线看片| 色网站视频免费| 少妇丰满av| av电影中文网址| 乱人伦中国视频| 精品国产露脸久久av麻豆| 欧美日韩在线观看h| 国产亚洲午夜精品一区二区久久| 欧美 亚洲 国产 日韩一| 国产爽快片一区二区三区| 久久久精品免费免费高清| 欧美最新免费一区二区三区| 国产精品国产三级国产av玫瑰| 人妻人人澡人人爽人人| 亚洲五月色婷婷综合| 中文字幕免费在线视频6| 亚洲精品456在线播放app| 亚洲国产精品国产精品| 一个人看视频在线观看www免费| 2022亚洲国产成人精品| 激情五月婷婷亚洲| 少妇人妻精品综合一区二区| 久久狼人影院| 亚洲无线观看免费| 国产男女超爽视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 在线观看免费高清a一片| 欧美成人精品欧美一级黄| 国产探花极品一区二区| 欧美亚洲 丝袜 人妻 在线| 国模一区二区三区四区视频| 国产毛片在线视频| 高清av免费在线| 欧美3d第一页| 欧美xxⅹ黑人| 免费观看性生交大片5| 久久久久久久亚洲中文字幕| 欧美日韩精品成人综合77777| 国产老妇伦熟女老妇高清| 欧美+日韩+精品| 亚洲欧洲国产日韩| 亚洲欧美日韩卡通动漫| 久久99热6这里只有精品| 999精品在线视频| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人精品一区二区| 亚洲国产精品999| 天堂中文最新版在线下载| av免费在线看不卡| 久久综合国产亚洲精品| 欧美精品亚洲一区二区| 亚洲av不卡在线观看| 狠狠精品人妻久久久久久综合| 青春草国产在线视频| 91精品一卡2卡3卡4卡| 最新的欧美精品一区二区| 狂野欧美激情性xxxx在线观看| 插阴视频在线观看视频| 欧美丝袜亚洲另类| 久热这里只有精品99| 日韩一本色道免费dvd| 各种免费的搞黄视频| 国产69精品久久久久777片| 国产免费福利视频在线观看| 亚洲成人手机| 久久久久久久亚洲中文字幕| 国产伦理片在线播放av一区| 精品熟女少妇av免费看| 日韩伦理黄色片| 久久青草综合色| 国产69精品久久久久777片| 国产日韩欧美亚洲二区| 日韩电影二区| 精品国产一区二区三区久久久樱花| 精品一区二区免费观看| 春色校园在线视频观看| 黄色欧美视频在线观看| 黑丝袜美女国产一区| 日韩三级伦理在线观看| 91久久精品国产一区二区成人| 亚洲av综合色区一区| 少妇高潮的动态图| 18禁在线无遮挡免费观看视频| 韩国高清视频一区二区三区| 国产又色又爽无遮挡免| 男女免费视频国产| 欧美精品国产亚洲| 亚洲精品日本国产第一区| 秋霞在线观看毛片| 国产精品99久久99久久久不卡 | 免费观看无遮挡的男女| 内地一区二区视频在线| 成年美女黄网站色视频大全免费 | 极品少妇高潮喷水抽搐| 国产有黄有色有爽视频| 99热网站在线观看| 99久久综合免费| 自线自在国产av| 亚洲经典国产精华液单| 69精品国产乱码久久久| 热re99久久国产66热| 少妇精品久久久久久久| 男的添女的下面高潮视频| 视频中文字幕在线观看| 亚洲av在线观看美女高潮| 国产av精品麻豆| 国产亚洲欧美精品永久| 国产精品熟女久久久久浪| 一级爰片在线观看| 母亲3免费完整高清在线观看 | 久久韩国三级中文字幕| 亚洲五月色婷婷综合| 一区二区av电影网| 国产成人免费无遮挡视频| 精品亚洲乱码少妇综合久久| 美女脱内裤让男人舔精品视频| 777米奇影视久久| 丝瓜视频免费看黄片| 欧美日韩在线观看h| 亚洲人与动物交配视频| 我的老师免费观看完整版| 能在线免费看毛片的网站| 欧美激情极品国产一区二区三区 | 亚洲高清免费不卡视频| 久久久国产精品麻豆| 欧美日韩视频高清一区二区三区二| 成人午夜精彩视频在线观看| √禁漫天堂资源中文www| 欧美精品高潮呻吟av久久| 亚洲美女视频黄频| 九九爱精品视频在线观看| 亚洲熟女精品中文字幕| 美女视频免费永久观看网站| 日本av免费视频播放| 水蜜桃什么品种好| a级毛片在线看网站| 国产精品嫩草影院av在线观看| 色网站视频免费| 91精品国产国语对白视频| 99re6热这里在线精品视频| 亚洲欧洲精品一区二区精品久久久 | 久久久久久久久久久免费av| 卡戴珊不雅视频在线播放| 一个人看视频在线观看www免费| 搡老乐熟女国产| 午夜免费观看性视频| 亚洲av福利一区| 午夜福利视频在线观看免费| 亚洲一级一片aⅴ在线观看| 亚洲精品乱码久久久v下载方式| 欧美人与性动交α欧美精品济南到 | 哪个播放器可以免费观看大片| 99热全是精品| 一区二区三区精品91| 久久精品夜色国产| 最新的欧美精品一区二区| 国产黄色视频一区二区在线观看| 一级,二级,三级黄色视频| 精品一区二区三区视频在线| 中文天堂在线官网| 51国产日韩欧美| 亚洲av欧美aⅴ国产| 老司机影院成人| 精品熟女少妇av免费看| 日韩亚洲欧美综合| 天天操日日干夜夜撸| 国产精品嫩草影院av在线观看| 欧美日韩一区二区视频在线观看视频在线| 高清欧美精品videossex| 国产av码专区亚洲av| 丰满少妇做爰视频| 狠狠婷婷综合久久久久久88av| 日韩欧美精品免费久久| 一级毛片黄色毛片免费观看视频| 国产精品久久久久成人av| 一级毛片黄色毛片免费观看视频| 欧美 亚洲 国产 日韩一| 全区人妻精品视频| 国产精品嫩草影院av在线观看| 精品熟女少妇av免费看| 中文乱码字字幕精品一区二区三区| 日本-黄色视频高清免费观看| 亚洲欧美清纯卡通| 三级国产精品欧美在线观看| 国产 精品1| 欧美精品人与动牲交sv欧美| 少妇 在线观看| 男女边摸边吃奶| 久久精品夜色国产| 国产精品偷伦视频观看了| 韩国av在线不卡| 日产精品乱码卡一卡2卡三| 韩国av在线不卡| 最新的欧美精品一区二区| 秋霞在线观看毛片| 91aial.com中文字幕在线观看| 热re99久久精品国产66热6| av播播在线观看一区| 欧美丝袜亚洲另类| 看非洲黑人一级黄片| 三上悠亚av全集在线观看| 久久久久精品久久久久真实原创| 伊人久久国产一区二区| 日韩中文字幕视频在线看片| 精品一品国产午夜福利视频| 搡老乐熟女国产| 中国美白少妇内射xxxbb| 9色porny在线观看| 国产探花极品一区二区| 久久韩国三级中文字幕| 国产有黄有色有爽视频| 成年人午夜在线观看视频| 成人毛片a级毛片在线播放| 水蜜桃什么品种好| 国产精品偷伦视频观看了| 啦啦啦啦在线视频资源| 少妇丰满av| 国产深夜福利视频在线观看| 97超视频在线观看视频| 亚洲精品久久成人aⅴ小说 | 日韩视频在线欧美| av一本久久久久| 熟女av电影| 午夜免费鲁丝| 国产精品一区二区在线观看99| 成人无遮挡网站| av.在线天堂| 久久久久网色| 国产精品国产三级专区第一集| 日本黄色片子视频| 人人妻人人澡人人看| 新久久久久国产一级毛片| 97超视频在线观看视频| h视频一区二区三区| a级毛片免费高清观看在线播放| 97超碰精品成人国产| 天堂8中文在线网| 国产欧美亚洲国产| 一级毛片我不卡| 亚洲欧美精品自产自拍| 少妇人妻精品综合一区二区| 亚洲精品日韩av片在线观看| 精品亚洲乱码少妇综合久久| 午夜激情久久久久久久| 日韩成人伦理影院| 久久国产亚洲av麻豆专区| 免费高清在线观看日韩| 五月天丁香电影| 九色成人免费人妻av| 久久久精品免费免费高清| 在线观看国产h片| 一个人看视频在线观看www免费| 中文字幕av电影在线播放| 免费高清在线观看日韩| 中文乱码字字幕精品一区二区三区| 日韩一本色道免费dvd| 久久99热6这里只有精品| 日韩中字成人| 九九爱精品视频在线观看| 日本wwww免费看| 国产一区有黄有色的免费视频| 免费久久久久久久精品成人欧美视频 | 国产免费视频播放在线视频| 久久鲁丝午夜福利片| 久久久久人妻精品一区果冻| 伦精品一区二区三区| 91精品国产国语对白视频| 国产精品99久久久久久久久| 人人澡人人妻人| 五月开心婷婷网| 日韩免费高清中文字幕av| 国产成人精品无人区| 韩国av在线不卡| 欧美日韩精品成人综合77777| 久久久久久久久久久丰满| 大码成人一级视频| 日韩成人伦理影院| 亚洲成人手机| 天天操日日干夜夜撸| 建设人人有责人人尽责人人享有的| 亚洲熟女精品中文字幕| 乱码一卡2卡4卡精品| 午夜久久久在线观看| 国产精品不卡视频一区二区| av卡一久久| 亚洲av日韩在线播放| 韩国av在线不卡| 日韩伦理黄色片| 国产色爽女视频免费观看| 夫妻午夜视频| 九草在线视频观看| 大香蕉97超碰在线| 免费人妻精品一区二区三区视频| 国产黄色视频一区二区在线观看| 亚洲av欧美aⅴ国产| 亚洲精品日本国产第一区| 夜夜看夜夜爽夜夜摸| 毛片一级片免费看久久久久| 国产成人freesex在线| 亚洲av电影在线观看一区二区三区| 美女中出高潮动态图| 中国三级夫妇交换| 国产乱人偷精品视频| 亚洲美女视频黄频| 亚洲国产精品国产精品| 日韩不卡一区二区三区视频在线| 91精品国产九色| 久久99蜜桃精品久久| 亚洲精品成人av观看孕妇| 精品人妻偷拍中文字幕| 夜夜骑夜夜射夜夜干| 免费观看av网站的网址| 老司机亚洲免费影院| 国产成人精品久久久久久| 一本大道久久a久久精品| 久久免费观看电影| 18禁观看日本| 高清av免费在线| 亚洲在久久综合| 国产在视频线精品| 久久99精品国语久久久| 国产乱人偷精品视频| 欧美国产精品一级二级三级| 大香蕉久久网| av女优亚洲男人天堂| 欧美日韩视频精品一区| 久久久久久久亚洲中文字幕| 久久人人爽av亚洲精品天堂| 国产有黄有色有爽视频| 人妻系列 视频| 亚洲国产av影院在线观看| 99久久精品一区二区三区| 一区二区三区四区激情视频| 久久99一区二区三区| 国产精品久久久久久精品古装| 成人国产av品久久久| 中国三级夫妇交换| a级片在线免费高清观看视频| 大片免费播放器 马上看| 91久久精品电影网| 国产欧美日韩一区二区三区在线 | 亚洲性久久影院| 最后的刺客免费高清国语| 日韩人妻高清精品专区| 九九爱精品视频在线观看| 午夜久久久在线观看| www.av在线官网国产| 国产片特级美女逼逼视频| 哪个播放器可以免费观看大片| 欧美精品高潮呻吟av久久| 久久人人爽av亚洲精品天堂| av.在线天堂| 国产69精品久久久久777片| 十分钟在线观看高清视频www| 亚洲欧美成人精品一区二区| 欧美日韩国产mv在线观看视频| 亚洲一区二区三区欧美精品| 精品亚洲乱码少妇综合久久| 国产高清国产精品国产三级| 99国产综合亚洲精品| 在线免费观看不下载黄p国产| 搡女人真爽免费视频火全软件| 精品人妻在线不人妻| av卡一久久| 亚洲精品乱码久久久久久按摩| 亚洲精品乱码久久久v下载方式| 黑人猛操日本美女一级片| 亚洲av二区三区四区| 黄色视频在线播放观看不卡| 一区二区av电影网| 国产不卡av网站在线观看| 一区二区三区免费毛片| 超色免费av| 久久久国产欧美日韩av| 天堂8中文在线网| 制服丝袜香蕉在线| 亚洲av不卡在线观看| 精品久久国产蜜桃| 免费av不卡在线播放| 麻豆乱淫一区二区| 97超视频在线观看视频| 伊人久久精品亚洲午夜| 亚洲欧美中文字幕日韩二区| 搡女人真爽免费视频火全软件| 欧美变态另类bdsm刘玥| 午夜老司机福利剧场| 欧美日韩成人在线一区二区| 少妇精品久久久久久久| 最近中文字幕高清免费大全6| 高清黄色对白视频在线免费看| 女的被弄到高潮叫床怎么办| 多毛熟女@视频| 国产高清不卡午夜福利| 热re99久久精品国产66热6| 在线观看三级黄色| 五月开心婷婷网| 亚洲五月色婷婷综合| 日韩大片免费观看网站| 免费看光身美女| 国产精品一区www在线观看| 久久久久精品久久久久真实原创| 国产午夜精品一二区理论片| 国语对白做爰xxxⅹ性视频网站| 超色免费av| 国产成人精品一,二区| 亚洲精品av麻豆狂野| 日韩av不卡免费在线播放| 性色av一级| 一区二区av电影网| 久久午夜综合久久蜜桃| 成人午夜精彩视频在线观看| 中国国产av一级| 黄色配什么色好看| 91精品伊人久久大香线蕉| 美女脱内裤让男人舔精品视频| 国产精品熟女久久久久浪| 国产永久视频网站| 国产视频内射| 亚洲天堂av无毛| 国产乱人偷精品视频| 国产欧美另类精品又又久久亚洲欧美| 美女中出高潮动态图| 美女xxoo啪啪120秒动态图| 免费看光身美女| 亚洲婷婷狠狠爱综合网| 久久人妻熟女aⅴ| 欧美精品亚洲一区二区| 欧美精品一区二区免费开放| 婷婷成人精品国产| 欧美bdsm另类| videos熟女内射| 亚洲国产欧美在线一区| 免费人成在线观看视频色| 男女国产视频网站| 伊人久久精品亚洲午夜| videossex国产| 两个人的视频大全免费| 一边摸一边做爽爽视频免费|