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

    帶時(shí)依協(xié)變量的重復(fù)測(cè)量資料的混合線性模型分析及其MIXED過(guò)程實(shí)現(xiàn)*

    2012-12-04 02:59:22上海交通大學(xué)醫(yī)學(xué)院生物統(tǒng)計(jì)學(xué)教研室200025張莉娜
    關(guān)鍵詞:協(xié)方差線性效應(yīng)

    上海交通大學(xué)醫(yī)學(xué)院生物統(tǒng)計(jì)學(xué)教研室(200025) 張莉娜

    重復(fù)測(cè)量資料是對(duì)同一觀察對(duì)象的同一觀察指標(biāo)在不同時(shí)間點(diǎn)上進(jìn)行的多次測(cè)量。其特點(diǎn)是存在時(shí)間上的自相關(guān)性?;谒迫缓瘮?shù)法原理的混合線性模型弱化了常規(guī)設(shè)計(jì)方法所要求的獨(dú)立性假定,是一般線性模型的擴(kuò)展,允許資料存在某種相關(guān)性以及協(xié)方差矩陣的多樣性,從而更好地適應(yīng)重復(fù)測(cè)量資料的特點(diǎn)。

    在重復(fù)測(cè)量設(shè)計(jì)中,除了由設(shè)計(jì)者所施加的處理因素外,還可能考慮到影響反應(yīng)變量的非處理因素,即協(xié)變量。協(xié)變量有兩種,一種是受試者在試驗(yàn)開(kāi)始前就有的,稱固定協(xié)變量或受試者別協(xié)變量;另一種協(xié)變量隨時(shí)間變化而變化,稱為時(shí)依協(xié)變量。

    本文通過(guò)一個(gè)臨床試驗(yàn)實(shí)例來(lái)對(duì)帶時(shí)依協(xié)變量的重復(fù)測(cè)量資料進(jìn)行混合線性模型分析,并給出MIXED的過(guò)程實(shí)現(xiàn)。

    資 料

    為了研究某新藥治療輕、中度原發(fā)性高血壓的療效和安全性,以坐位舒張壓為主要療效指標(biāo),采用隨機(jī)、多中心、陽(yáng)性藥物平行對(duì)照、雙盲雙模擬臨床試驗(yàn),收集234例輕、中度原發(fā)性高血壓患者隨機(jī)分到1,2兩組,分別采用對(duì)照藥和試驗(yàn)藥。所有符合試驗(yàn)方案、依從性良好、試驗(yàn)期間未服用禁止用藥、完成CRF的病例納入PP(per protocol)分析集,由于無(wú)效而提前退出的病例也納入PP分析。本試驗(yàn)的療效分析用PP分析。

    給藥方案:療程總計(jì)10周,兩組均按導(dǎo)入期(2周)和治療期(8周)進(jìn)行。給藥途徑為口服。受試者每天于早餐后(8時(shí)左右)服藥。導(dǎo)入期每日1次、每次試驗(yàn)和對(duì)照藥的安慰劑各1片;治療期頭2周每日1次、每次服試驗(yàn)或?qū)φ账幖捌淠M片各1片。若2、4、6周其中一個(gè)周末血壓≥140/90 mmHg,增加至每次各2片(仍為每日1次),直至總療程8周結(jié)束。

    分別在治療前、治療2周、治療4周、治療6周、治療8周5個(gè)時(shí)間點(diǎn)記錄每個(gè)受試者的坐位舒張壓和用藥量,數(shù)據(jù)結(jié)構(gòu)見(jiàn)表1。

    表1 帶時(shí)依協(xié)變量的重復(fù)測(cè)量資料數(shù)據(jù)結(jié)構(gòu)

    方 法

    混合線性模型(mixed linear model)

    記第i例受試者在時(shí)間點(diǎn)j(j=1,…,p)的反應(yīng)變量觀察值為yij,協(xié)變量觀察值為zij,則數(shù)據(jù)集表示為[(yi1,zi1)(yi2,zi2),…,(yip,zip)],i=1,…,n

    給定符號(hào)Y為反應(yīng)變量矩陣,X為設(shè)計(jì)矩陣,Z為時(shí)依協(xié)變量矩陣。建立線性關(guān)系式

    式中β為固定效應(yīng)的參數(shù)矩陣,γ為時(shí)依協(xié)變量的參數(shù)矩陣,其主對(duì)角線元素為 γjj,(j=1,…,p)、非主對(duì)角線元素為0。下標(biāo)為矩陣維數(shù)。展開(kāi)式為:

    其中等式右邊的第一部分是設(shè)計(jì)效應(yīng),中間部分是時(shí)依協(xié)變量效應(yīng)。γ11是第1個(gè)時(shí)間點(diǎn)的時(shí)依協(xié)變量 Z1=(z1p,z21,…,zn1)'的參數(shù),γpp是第 p 個(gè)時(shí)間點(diǎn)的時(shí)依協(xié)變量 Zp=(z1p,z21,…,znp)'的參數(shù)。

    在混合線性模型中,應(yīng)用最大似然法或約束最大似然法和最小方差二次無(wú)偏估計(jì)對(duì)參數(shù)進(jìn)行估計(jì)。

    比起一般線性模型,混合線性模型的優(yōu)勢(shì)在于假定協(xié)方差具有某種形式的結(jié)構(gòu),既不會(huì)像單變量分析方法那樣嚴(yán)格,也不會(huì)像多變量方差分析那樣對(duì)協(xié)方差矩陣完全無(wú)約束,從而更好地適應(yīng)重復(fù)測(cè)量資料的特點(diǎn)。而選擇正確的協(xié)方差結(jié)構(gòu)將影響固定效應(yīng)參數(shù)估計(jì)值的準(zhǔn)確性,尤其是影響其標(biāo)準(zhǔn)誤的計(jì)算,所以是擬合該混合模型時(shí)最關(guān)鍵的一步,本文給出幾種常用的協(xié)方差矩陣以備選用。

    of order 1[AR(1)]

    軟件實(shí)現(xiàn)

    利用SAS 9.1中的MIXED過(guò)程來(lái)實(shí)現(xiàn)帶時(shí)依協(xié)變量的重復(fù)測(cè)量資料混合線性模型擬合與參數(shù)估計(jì)。no表示受試者編號(hào),group表示分組(1=對(duì)照組,2=試驗(yàn)組),time表示不同重測(cè)時(shí)間點(diǎn),以各時(shí)間點(diǎn)坐位舒張壓的下降值為反應(yīng)變量,各時(shí)間點(diǎn)用藥量的改變值為時(shí)依協(xié)變量擬合模型。SAS程序如下:

    結(jié)果討論

    協(xié)方差結(jié)構(gòu)的選擇可以通過(guò)似然比統(tǒng)計(jì)量(-2Log Likelihood)、Akaike'Information Criterion(AIC)、Schwart's Bayesian Criterion(BIC)來(lái)判斷,其中主要是AIC和BIC,這兩個(gè)統(tǒng)計(jì)量的值越小說(shuō)明模型擬合得越好,從而選取相應(yīng)的協(xié)方差結(jié)構(gòu)。在統(tǒng)計(jì)量的值相近時(shí),則選取含參數(shù)個(gè)數(shù)最少的一個(gè)。

    表2 不同協(xié)方差結(jié)構(gòu)下的各種擬合優(yōu)度檢驗(yàn)統(tǒng)計(jì)量

    表2列出了不同協(xié)方差結(jié)構(gòu)下的各種擬合優(yōu)度檢驗(yàn)統(tǒng)計(jì)量,經(jīng)過(guò)綜合比較,選取無(wú)結(jié)構(gòu)型協(xié)方差矩陣(UN)。

    模型的最大似然檢驗(yàn):χ2=355.88,P <0.0001。說(shuō)明無(wú)結(jié)構(gòu)型協(xié)方差矩陣要優(yōu)于普通常方差的最小二乘估計(jì),模型擬合效果顯著。

    表3 固定效應(yīng)的3型檢驗(yàn)

    表3結(jié)果顯示:兩組平均坐位舒張壓下降值的差別無(wú)統(tǒng)計(jì)學(xué)意義(P=0.4475),即兩種藥對(duì)治療輕、中度原發(fā)性高血壓病的療效無(wú)差別;各時(shí)間點(diǎn)之間差異有統(tǒng)計(jì)學(xué)意義(P<0.0001);分組效應(yīng)與時(shí)間效應(yīng)的交互作用無(wú)統(tǒng)計(jì)學(xué)意義(P=0.1523),說(shuō)明試驗(yàn)組和對(duì)照組的坐位舒張壓下降值隨時(shí)間變化趨勢(shì)的差別無(wú)統(tǒng)計(jì)學(xué)意義;協(xié)變量×?xí)r間交互效應(yīng)有統(tǒng)計(jì)學(xué)意義(P=0.0481),說(shuō)明各時(shí)間點(diǎn)的用藥方案對(duì)降壓療效的影響有顯著效應(yīng);協(xié)變量×分組交互效應(yīng)有統(tǒng)計(jì)學(xué)意義(P=0.0075),說(shuō)明兩種藥的用藥方案對(duì)降壓療效的影響有顯著效應(yīng);時(shí)依協(xié)變量有統(tǒng)計(jì)學(xué)意義(P=0.0061),即給藥方案對(duì)治療輕、中度原發(fā)性高血壓病有顯著效應(yīng)。

    表4 固定效應(yīng)的參數(shù)估計(jì)值

    表4結(jié)果顯示:試驗(yàn)組相比對(duì)照組對(duì)坐位舒張壓下降值的影響為正效應(yīng),P=0.091接近顯著性水平,第8周相比第2周、第4周對(duì)坐位舒張壓下降值的影響為正效應(yīng),P值均小于0.05,第8周相比第6周對(duì)坐位舒張壓下降值的影響無(wú)統(tǒng)計(jì)學(xué)意義,對(duì)照組相比試驗(yàn)組給藥方案對(duì)坐位舒張壓下降值的影響為正效應(yīng),P=0.0072。

    表5 各時(shí)間點(diǎn)的兩兩比較

    表5結(jié)果顯示:第6周和第8周無(wú)統(tǒng)計(jì)學(xué)差異,其余各時(shí)間點(diǎn)都有統(tǒng)計(jì)學(xué)差異。

    表6結(jié)果顯示各組在各時(shí)間點(diǎn)上的坐位舒張壓下降值修正均數(shù)、標(biāo)準(zhǔn)誤、95%的可信區(qū)間。

    圖1從直觀上提示試驗(yàn)組在各時(shí)間點(diǎn)的坐位舒張壓下降值的增長(zhǎng)幅度略大于對(duì)照組。

    總之,藥物臨床試驗(yàn)大多需要定期收集病例數(shù)據(jù),在統(tǒng)計(jì)策略上應(yīng)該首先考慮混合模型。

    表6 經(jīng)x校正后兩組在各時(shí)間點(diǎn)的平均坐位舒張壓下降值

    混合線性模型中提供的多種協(xié)方差結(jié)構(gòu),在數(shù)據(jù)之間相關(guān)性較為復(fù)雜時(shí),能夠更好地?cái)M合數(shù)據(jù)之間存在的相關(guān)性,更適合于帶時(shí)依協(xié)變量的重復(fù)測(cè)量資料。SAS中的MIXED過(guò)程能充分挖掘數(shù)據(jù)信息,揭示其內(nèi)在特點(diǎn)并給出合理的統(tǒng)計(jì)結(jié)論。

    圖1 兩組坐位舒張壓下降值的修正均數(shù)在各時(shí)間點(diǎn)的比較

    1.Helen Brown,Robin Prescott.Applied Mixed Model in Medicine.Chichester:John Wiley & SonsLtd,1999,199-259.

    2.Littell RC,Miliken GA,Stroup WW,et al.SAS for Mixed Models.Second dition.Cary,NC,USA:SAS Institute Inc,2006:174-177.

    3.Khattree R,Naik DN.Applied Multivariate Statistics with SAS Software.Second Edition.Cary,NC,USA:SAS Institute Inc,1999:278-287.

    4.施紅英,沈毅.混合模型在臨床試驗(yàn)重復(fù)測(cè)量資料中的應(yīng)用.中國(guó)衛(wèi)生統(tǒng)計(jì),2007,24(2):140-142.

    5.高萌,張強(qiáng).方差分析和混合線性模型在重復(fù)測(cè)量數(shù)據(jù)中的應(yīng)用探討.現(xiàn)代預(yù)防醫(yī)學(xué),2008,35(7):1221-1223.

    6.余松林,向惠云,編著.重復(fù)測(cè)量資料分析方法與SAS程序.北京:科學(xué)出版社,2004,3.

    7.高惠璇編譯.SAS系統(tǒng)SAS/STAT軟件使用手冊(cè).北京:中國(guó)統(tǒng)計(jì)出版社,1995,4.

    猜你喜歡
    協(xié)方差線性效應(yīng)
    漸近線性Klein-Gordon-Maxwell系統(tǒng)正解的存在性
    鈾對(duì)大型溞的急性毒性效應(yīng)
    線性回歸方程的求解與應(yīng)用
    懶馬效應(yīng)
    二階線性微分方程的解法
    應(yīng)變效應(yīng)及其應(yīng)用
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    縱向數(shù)據(jù)分析中使用滑動(dòng)平均Cholesky分解對(duì)回歸均值和協(xié)方差矩陣進(jìn)行同時(shí)半?yún)?shù)建模
    關(guān)于協(xié)方差的U統(tǒng)計(jì)量檢驗(yàn)法
    亚洲欧美精品自产自拍| 熟女电影av网| 精品国产国语对白av| 亚洲欧洲日产国产| 寂寞人妻少妇视频99o| 免费高清在线观看视频在线观看| 国产成人精品无人区| 亚洲成色77777| 亚洲精品一二三| 国产精品久久久久久久久免| 热re99久久国产66热| 国产精品三级大全| 欧美日韩一级在线毛片| 香蕉丝袜av| 尾随美女入室| 一个人免费看片子| 少妇人妻精品综合一区二区| 人人澡人人妻人| 可以免费在线观看a视频的电影网站 | 婷婷色麻豆天堂久久| 午夜福利在线观看免费完整高清在| 男女午夜视频在线观看| 久久精品国产亚洲av天美| 三级国产精品片| 午夜福利影视在线免费观看| 毛片一级片免费看久久久久| 综合色丁香网| 久久久久久久精品精品| 日本wwww免费看| 97在线视频观看| 男女免费视频国产| 热99久久久久精品小说推荐| 熟妇人妻不卡中文字幕| 一级a爱视频在线免费观看| 一级毛片我不卡| 久久青草综合色| 黑人猛操日本美女一级片| 在线观看美女被高潮喷水网站| 七月丁香在线播放| 午夜日本视频在线| 国产女主播在线喷水免费视频网站| av有码第一页| 久久久欧美国产精品| 久久国产精品男人的天堂亚洲| 中文字幕制服av| 老汉色∧v一级毛片| 国产精品免费视频内射| a 毛片基地| 久久精品夜色国产| 亚洲精品第二区| 久久人人爽人人片av| 另类精品久久| 亚洲美女搞黄在线观看| 亚洲激情五月婷婷啪啪| 最近中文字幕2019免费版| 亚洲成人av在线免费| 岛国毛片在线播放| 性高湖久久久久久久久免费观看| 黄色一级大片看看| 看十八女毛片水多多多| 十八禁网站网址无遮挡| 高清视频免费观看一区二区| 久久久久精品人妻al黑| 国产精品.久久久| 亚洲成人av在线免费| 国产av码专区亚洲av| 成人亚洲精品一区在线观看| 看免费成人av毛片| 国产精品.久久久| 久久精品国产亚洲av高清一级| 天天操日日干夜夜撸| 丝袜美腿诱惑在线| 99re6热这里在线精品视频| 黄网站色视频无遮挡免费观看| 大片免费播放器 马上看| 少妇人妻久久综合中文| 国产极品粉嫩免费观看在线| 亚洲一码二码三码区别大吗| 中文字幕人妻丝袜制服| 国产女主播在线喷水免费视频网站| 久久精品国产亚洲av高清一级| 久久99一区二区三区| 纵有疾风起免费观看全集完整版| 在线天堂最新版资源| 久久精品人人爽人人爽视色| 久久久久久久久久人人人人人人| 十八禁网站网址无遮挡| 欧美日本中文国产一区发布| 久久亚洲国产成人精品v| 亚洲国产精品999| 巨乳人妻的诱惑在线观看| 亚洲图色成人| 国产在线一区二区三区精| 久久人人97超碰香蕉20202| 午夜精品国产一区二区电影| 妹子高潮喷水视频| 日韩欧美精品免费久久| 青青草视频在线视频观看| 男女高潮啪啪啪动态图| 国产激情久久老熟女| 又粗又硬又长又爽又黄的视频| 国产极品天堂在线| 欧美精品人与动牲交sv欧美| 国产免费福利视频在线观看| 18在线观看网站| 欧美成人精品欧美一级黄| 精品人妻一区二区三区麻豆| 在线观看人妻少妇| 亚洲,欧美精品.| 90打野战视频偷拍视频| av网站免费在线观看视频| 午夜福利在线观看免费完整高清在| 边亲边吃奶的免费视频| 日韩av不卡免费在线播放| 大香蕉久久成人网| 色婷婷av一区二区三区视频| 夫妻午夜视频| 久久热在线av| 国产日韩一区二区三区精品不卡| 国产免费福利视频在线观看| 一区二区三区精品91| 女人精品久久久久毛片| 大片电影免费在线观看免费| 大片免费播放器 马上看| 久久久久网色| 999精品在线视频| 久久精品国产a三级三级三级| 日韩一区二区视频免费看| 久久这里有精品视频免费| 欧美变态另类bdsm刘玥| 熟女av电影| 人体艺术视频欧美日本| 中文天堂在线官网| 免费人妻精品一区二区三区视频| 久久这里只有精品19| 亚洲综合色网址| 国产日韩一区二区三区精品不卡| 少妇猛男粗大的猛烈进出视频| 老司机亚洲免费影院| 妹子高潮喷水视频| 亚洲国产成人一精品久久久| 极品少妇高潮喷水抽搐| 另类亚洲欧美激情| 成人免费观看视频高清| 亚洲欧美精品自产自拍| 国产精品成人在线| 又大又黄又爽视频免费| 一级毛片 在线播放| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品国产av成人精品| 在线 av 中文字幕| 三级国产精品片| 久久久久久人妻| 欧美97在线视频| 国产不卡av网站在线观看| 欧美亚洲 丝袜 人妻 在线| 国产成人精品久久久久久| 久久久久久久久久久久大奶| 欧美亚洲 丝袜 人妻 在线| 一区二区av电影网| 青春草亚洲视频在线观看| 日日爽夜夜爽网站| 亚洲精品日本国产第一区| 国产成人欧美| 国产精品嫩草影院av在线观看| 咕卡用的链子| 91久久精品国产一区二区三区| 成人影院久久| 在线观看免费高清a一片| 婷婷色综合www| 国产精品国产三级国产专区5o| 国产精品久久久久久久久免| 久久精品国产a三级三级三级| 国产色婷婷99| 国产欧美亚洲国产| 美女大奶头黄色视频| 日韩av在线免费看完整版不卡| 男女边摸边吃奶| 国产精品一国产av| 久久久久视频综合| www.熟女人妻精品国产| 美女福利国产在线| 水蜜桃什么品种好| 男女午夜视频在线观看| 赤兔流量卡办理| 亚洲经典国产精华液单| 亚洲激情五月婷婷啪啪| 搡女人真爽免费视频火全软件| 国产一区二区在线观看av| 黄片无遮挡物在线观看| 观看av在线不卡| 日韩伦理黄色片| av网站在线播放免费| 天堂8中文在线网| 久久97久久精品| 美国免费a级毛片| 老司机亚洲免费影院| 人妻系列 视频| 18禁观看日本| 性色av一级| 亚洲国产色片| 日韩中文字幕欧美一区二区 | 中文字幕人妻丝袜制服| 亚洲人成电影观看| 美女国产视频在线观看| 高清不卡的av网站| 9热在线视频观看99| 久久精品国产鲁丝片午夜精品| 热99久久久久精品小说推荐| 国产高清不卡午夜福利| 亚洲av欧美aⅴ国产| 久久久久久伊人网av| 最近中文字幕高清免费大全6| a 毛片基地| 国产精品麻豆人妻色哟哟久久| 久久精品国产亚洲av涩爱| 欧美97在线视频| 赤兔流量卡办理| 亚洲精品aⅴ在线观看| 一级片'在线观看视频| 亚洲欧美中文字幕日韩二区| 一区福利在线观看| 夫妻性生交免费视频一级片| 1024视频免费在线观看| 亚洲av.av天堂| 啦啦啦在线免费观看视频4| 超碰成人久久| 啦啦啦啦在线视频资源| 亚洲精品久久成人aⅴ小说| 叶爱在线成人免费视频播放| 久久99一区二区三区| 国产成人91sexporn| 日本av手机在线免费观看| 日本欧美国产在线视频| 精品99又大又爽又粗少妇毛片| 久久久欧美国产精品| 大话2 男鬼变身卡| 日本午夜av视频| 激情五月婷婷亚洲| 免费日韩欧美在线观看| 亚洲少妇的诱惑av| 中文字幕人妻丝袜一区二区 | 少妇的逼水好多| 免费大片黄手机在线观看| 亚洲av日韩在线播放| 亚洲熟女精品中文字幕| 亚洲国产av新网站| 久久久久久久久免费视频了| 少妇 在线观看| 另类亚洲欧美激情| 精品福利永久在线观看| 最黄视频免费看| 91在线精品国自产拍蜜月| 久久精品久久久久久久性| 国产在线一区二区三区精| 亚洲成色77777| 日韩精品有码人妻一区| 亚洲综合色惰| 亚洲国产欧美网| 精品99又大又爽又粗少妇毛片| 久久这里有精品视频免费| 亚洲经典国产精华液单| 亚洲熟女精品中文字幕| 亚洲av男天堂| 黑人欧美特级aaaaaa片| 老司机亚洲免费影院| 免费黄网站久久成人精品| 精品少妇黑人巨大在线播放| 老司机影院毛片| 亚洲美女黄色视频免费看| 一本色道久久久久久精品综合| 免费观看性生交大片5| 少妇熟女欧美另类| 精品亚洲成国产av| 久久久久久人妻| 午夜免费男女啪啪视频观看| 一级毛片黄色毛片免费观看视频| 欧美日韩精品成人综合77777| 一级片'在线观看视频| 大片免费播放器 马上看| 性色av一级| 久久这里有精品视频免费| 日本-黄色视频高清免费观看| 爱豆传媒免费全集在线观看| 国产免费一区二区三区四区乱码| 免费黄网站久久成人精品| 日韩,欧美,国产一区二区三区| 哪个播放器可以免费观看大片| 夫妻午夜视频| 9热在线视频观看99| 欧美变态另类bdsm刘玥| 国产成人精品无人区| 国产男女内射视频| 亚洲,一卡二卡三卡| 国产日韩欧美视频二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区综合在线观看| 交换朋友夫妻互换小说| 亚洲av福利一区| 色婷婷久久久亚洲欧美| 亚洲精品在线美女| 久久鲁丝午夜福利片| 亚洲av男天堂| 日韩欧美精品免费久久| 亚洲图色成人| 亚洲成人av在线免费| 日韩一卡2卡3卡4卡2021年| 精品卡一卡二卡四卡免费| 三上悠亚av全集在线观看| 伦精品一区二区三区| 如日韩欧美国产精品一区二区三区| 人人澡人人妻人| 18禁动态无遮挡网站| 岛国毛片在线播放| 如日韩欧美国产精品一区二区三区| 丰满乱子伦码专区| 在线 av 中文字幕| 精品亚洲成a人片在线观看| 国产在线免费精品| 26uuu在线亚洲综合色| 天堂8中文在线网| 韩国高清视频一区二区三区| 日韩欧美一区视频在线观看| 免费观看在线日韩| 搡女人真爽免费视频火全软件| 看免费av毛片| 日韩中文字幕欧美一区二区 | 一区二区av电影网| 卡戴珊不雅视频在线播放| 欧美日韩亚洲高清精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 18在线观看网站| 男女国产视频网站| 日韩 亚洲 欧美在线| 欧美国产精品一级二级三级| 亚洲精品国产av成人精品| 欧美成人午夜免费资源| 美女高潮到喷水免费观看| 日韩伦理黄色片| 卡戴珊不雅视频在线播放| 日韩制服丝袜自拍偷拍| 国产男女超爽视频在线观看| 777米奇影视久久| 9热在线视频观看99| 国产淫语在线视频| 久久免费观看电影| av一本久久久久| 久热久热在线精品观看| 日本免费在线观看一区| 精品久久久久久电影网| 亚洲国产欧美网| 老女人水多毛片| 狠狠婷婷综合久久久久久88av| 亚洲成国产人片在线观看| 99久国产av精品国产电影| 欧美av亚洲av综合av国产av | 欧美日本中文国产一区发布| 国产极品粉嫩免费观看在线| 嫩草影院入口| 男女午夜视频在线观看| 你懂的网址亚洲精品在线观看| 大码成人一级视频| 久久久久精品性色| 午夜福利一区二区在线看| 日本-黄色视频高清免费观看| 国产成人精品婷婷| 色婷婷久久久亚洲欧美| 一区福利在线观看| 婷婷色麻豆天堂久久| 精品第一国产精品| 十八禁网站网址无遮挡| 免费高清在线观看日韩| 亚洲av男天堂| 一区二区三区乱码不卡18| 日韩精品免费视频一区二区三区| 欧美日韩成人在线一区二区| 啦啦啦视频在线资源免费观看| 久久久国产欧美日韩av| 午夜福利网站1000一区二区三区| 国产日韩一区二区三区精品不卡| 狠狠精品人妻久久久久久综合| 色94色欧美一区二区| 美女国产视频在线观看| 99热全是精品| 90打野战视频偷拍视频| 18禁裸乳无遮挡动漫免费视频| 人妻系列 视频| 免费高清在线观看日韩| 熟女少妇亚洲综合色aaa.| 亚洲欧美日韩另类电影网站| 精品酒店卫生间| 久久久久视频综合| 亚洲五月色婷婷综合| 亚洲国产日韩一区二区| 精品午夜福利在线看| 天堂俺去俺来也www色官网| videosex国产| 国产精品人妻久久久影院| 欧美成人午夜精品| 成人亚洲精品一区在线观看| 777米奇影视久久| 老汉色∧v一级毛片| 欧美激情高清一区二区三区 | 波多野结衣一区麻豆| 又黄又粗又硬又大视频| 精品人妻在线不人妻| av免费观看日本| 五月天丁香电影| 久久久久久久久久久免费av| 高清av免费在线| 999久久久国产精品视频| 伊人亚洲综合成人网| 在线观看美女被高潮喷水网站| 国产亚洲一区二区精品| av有码第一页| 亚洲欧洲精品一区二区精品久久久 | 大陆偷拍与自拍| av一本久久久久| 久久久久久人人人人人| 国产熟女欧美一区二区| 有码 亚洲区| 乱人伦中国视频| 欧美精品人与动牲交sv欧美| av在线播放精品| 日韩一区二区视频免费看| 日本av免费视频播放| 国产在视频线精品| 日韩一区二区视频免费看| 欧美亚洲 丝袜 人妻 在线| 最近中文字幕高清免费大全6| 国产乱人偷精品视频| 亚洲一区二区三区欧美精品| 黄色毛片三级朝国网站| 精品国产一区二区三区四区第35| 成人毛片60女人毛片免费| 欧美日韩视频高清一区二区三区二| 91成人精品电影| 永久免费av网站大全| 高清不卡的av网站| 伊人久久国产一区二区| 熟女电影av网| 亚洲情色 制服丝袜| 久久这里有精品视频免费| 亚洲国产av影院在线观看| 亚洲欧美一区二区三区国产| 亚洲精品美女久久av网站| 日韩中文字幕视频在线看片| 久久人人爽人人片av| 日本-黄色视频高清免费观看| 91精品国产国语对白视频| 80岁老熟妇乱子伦牲交| 亚洲美女视频黄频| 国产精品免费视频内射| 亚洲经典国产精华液单| 国产片内射在线| 亚洲少妇的诱惑av| 国产黄频视频在线观看| 久久午夜福利片| 欧美黄色片欧美黄色片| 成人国产av品久久久| 日韩电影二区| 日本色播在线视频| 亚洲四区av| 久久精品国产亚洲av涩爱| 各种免费的搞黄视频| 性色av一级| av.在线天堂| 一本大道久久a久久精品| 日韩电影二区| 亚洲综合色惰| 狠狠精品人妻久久久久久综合| 国产精品偷伦视频观看了| a 毛片基地| 大话2 男鬼变身卡| 一级毛片我不卡| 性高湖久久久久久久久免费观看| 中文字幕最新亚洲高清| 美女午夜性视频免费| 最近最新中文字幕免费大全7| 黑人巨大精品欧美一区二区蜜桃| 国产黄频视频在线观看| 一本大道久久a久久精品| 不卡视频在线观看欧美| 97人妻天天添夜夜摸| 国产在线一区二区三区精| 观看av在线不卡| 免费黄色在线免费观看| 在线天堂中文资源库| 亚洲精品乱久久久久久| 欧美成人午夜精品| 人人妻人人澡人人爽人人夜夜| 又粗又硬又长又爽又黄的视频| 国产亚洲一区二区精品| 成人亚洲欧美一区二区av| 欧美精品av麻豆av| 丝袜喷水一区| 美女主播在线视频| 国产精品久久久久久精品古装| 欧美日韩视频精品一区| 久久久久久免费高清国产稀缺| 丝袜人妻中文字幕| 一二三四在线观看免费中文在| 老鸭窝网址在线观看| 亚洲国产色片| freevideosex欧美| 一级黄片播放器| 久久久久久久久久久久大奶| 夜夜骑夜夜射夜夜干| 国产av码专区亚洲av| 亚洲 欧美一区二区三区| 国产在线视频一区二区| 久久久精品94久久精品| 欧美av亚洲av综合av国产av | 最近最新中文字幕大全免费视频 | 尾随美女入室| 欧美少妇被猛烈插入视频| 一区二区日韩欧美中文字幕| 国产激情久久老熟女| 秋霞伦理黄片| 亚洲综合色网址| 成年av动漫网址| 欧美日韩视频精品一区| 一边摸一边做爽爽视频免费| 亚洲国产精品一区二区三区在线| 咕卡用的链子| 丝袜美腿诱惑在线| 国产精品熟女久久久久浪| 国产福利在线免费观看视频| 日本欧美国产在线视频| 永久免费av网站大全| 亚洲欧美一区二区三区国产| 午夜日本视频在线| 精品少妇一区二区三区视频日本电影 | 91aial.com中文字幕在线观看| 激情视频va一区二区三区| 亚洲熟女精品中文字幕| 成年av动漫网址| 久久久久精品久久久久真实原创| 黑人巨大精品欧美一区二区蜜桃| 色吧在线观看| 精品国产一区二区三区四区第35| 亚洲国产欧美在线一区| 十八禁高潮呻吟视频| 亚洲伊人久久精品综合| 伊人久久国产一区二区| 亚洲综合色惰| 老女人水多毛片| 日韩三级伦理在线观看| 美女中出高潮动态图| 中文精品一卡2卡3卡4更新| 日本91视频免费播放| 中国三级夫妇交换| tube8黄色片| 精品午夜福利在线看| 成年女人在线观看亚洲视频| 午夜激情久久久久久久| 男女无遮挡免费网站观看| 18禁裸乳无遮挡动漫免费视频| 蜜桃国产av成人99| 成人手机av| 少妇人妻精品综合一区二区| 欧美精品人与动牲交sv欧美| 国产精品嫩草影院av在线观看| 肉色欧美久久久久久久蜜桃| 精品久久久精品久久久| 永久网站在线| www日本在线高清视频| 国产免费福利视频在线观看| 欧美 亚洲 国产 日韩一| 国产免费现黄频在线看| 久久久久久久精品精品| 成年人免费黄色播放视频| 天美传媒精品一区二区| 老汉色av国产亚洲站长工具| av天堂久久9| 99国产精品免费福利视频| 亚洲精品久久午夜乱码| videosex国产| 国语对白做爰xxxⅹ性视频网站| 国产一区二区激情短视频 | 青青草视频在线视频观看| 久久av网站| 久久99热这里只频精品6学生| 精品人妻熟女毛片av久久网站| 我的亚洲天堂| 欧美+日韩+精品| 综合色丁香网| 九草在线视频观看| av在线播放精品| 久久久久精品人妻al黑| 亚洲色图综合在线观看| 大香蕉久久网| 亚洲国产欧美网| 亚洲成人av在线免费| 夜夜骑夜夜射夜夜干| 国产成人精品久久久久久| 在现免费观看毛片| 永久网站在线| 女性被躁到高潮视频| 久久99一区二区三区| 曰老女人黄片| 男女啪啪激烈高潮av片| 一级黄片播放器| 亚洲精品,欧美精品| 一区二区三区激情视频| 中文欧美无线码| 一边摸一边做爽爽视频免费| 亚洲av国产av综合av卡| 秋霞在线观看毛片| 亚洲精品,欧美精品| 超碰97精品在线观看| 亚洲天堂av无毛| 18禁观看日本| 日韩精品有码人妻一区| 18禁观看日本|