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

    回歸校準和模擬外推對測量誤差的校正效果研究*

    2014-03-10 07:03:55霞張橋陳永杰李慧婷張秋菊劉美娜
    中國衛(wèi)生統(tǒng)計 2014年5期
    關(guān)鍵詞:外推法測量誤差回歸系數(shù)

    陳 霞張 橋陳永杰李慧婷張秋菊劉美娜△

    回歸校準和模擬外推對測量誤差的校正效果研究*

    陳 霞1張 橋2陳永杰1李慧婷1張秋菊1劉美娜1△

    目的探討回歸校準法(RC)和模擬外推法(SIMEX)對logistic回歸中測量誤差的校正效果。方法通過SAS軟件產(chǎn)生有測量誤差的模擬數(shù)據(jù),用RC和SIMEX對測量誤差進行校正,對比設(shè)定的真實β值和校正后的β*值之間的差別,以評價校正效果。結(jié)果當X可精確測量時,在設(shè)定的條件下,RC1的校正效果較好;P-SIMEX僅當很小時校正效果較好。當X不可精確測量時,隨著測量誤差的增大,E-SIMEX的校正效果降低,而RC2的校正效果相對較穩(wěn)定。結(jié)論無論X是否可測,在經(jīng)典測量誤差模型前提下,RC對logistic回歸模型中測量誤差的校正效果優(yōu)于SIMEX,建議應(yīng)用RC校正測量誤差。

    測量誤差 回歸校準法 模擬外推法 logistic回歸回歸稀釋

    統(tǒng)計分析過程中建立線性回歸或非線性回歸模型時,都假定自變量和應(yīng)變量無測量誤差,所得數(shù)據(jù)是變量的真實值。但在科研領(lǐng)域中,通過某種方式或者途徑所獲得的變量值未必是此變量的真實值,尤其是環(huán)境流行病學(xué)、膳食問卷調(diào)查等觀察性研究中,測量誤差是廣泛存在的問題[1,2]。數(shù)據(jù)分析時若忽略測量誤差的影響,會導(dǎo)致參數(shù)估計有偏,統(tǒng)計推斷不準確,因此需要對測量誤差進行校正[3]。關(guān)于測量誤差校正方法的研究,最經(jīng)典的是由Rosner提出的回歸校準法[4]和由Cook提出的模擬外推法[5],可通過SAS、Stata等統(tǒng)計軟件實現(xiàn)分析過程。本文通過模擬實驗,探討不同條件下回歸校準法和模擬外推法對logistic回歸中自變量測量誤差的校正效果。

    基本原理與方法

    1.測量誤差

    測量誤差是指科研過程中由于信息收集、固有的生物學(xué)變異和一些變量的主觀性等造成的誤差[6]。傳統(tǒng)回歸模型,如式(1):

    此模型中,變量Yi,Xi是可測的,εi是隨機誤差,β是待估參數(shù),fXi,()β是關(guān)于X的已知函數(shù)。通??筛鶕?jù)回歸函數(shù)fXi,()β的不同形式,估計出相應(yīng)的參數(shù)β。但在實際建模實踐時,自變量可能不可測,或者其測量帶有不可忽略的測量誤差。在一些觀察性研究中,如膳食頻率調(diào)查(FFQ),自變量X為真實值,是被調(diào)查者過去一年中的日平均膳食攝入量,通過問卷調(diào)查獲取的值為真實值X的估計值W,W和X之間的差別就是測量誤差[7]。

    2.回歸校準(regression calibration,RC)

    回歸校準的基本原理[8]是:先建立自變量真實值X關(guān)于觀測值W和協(xié)變量Z的回歸函數(shù),以此函數(shù)代替真實值X,再進行常規(guī)的回歸分析。RC的前提條件是W和X為線性關(guān)系,X關(guān)于W和Z的回歸函數(shù)叫做回歸校準函數(shù)。根據(jù)已知數(shù)據(jù)的特點,有RC1和RC2兩種方法建立回歸校準函數(shù)。

    RC1:當X可測,且測量值W與真實值X滿足非差分測量誤差模型[9](non-differentialmeasurement errormodel)W=γ0+γX+U,其中且U與X、Y相互獨立??蓮闹鲾?shù)據(jù)集中按一定比例隨機抽取一個同時有W和X值的數(shù)據(jù)集,稱效度數(shù)據(jù)集(validation data)。其較準過程分三步:首先利用效度數(shù)據(jù)集,建立X與W的線性模型X=λ0+λ1W,估算出系數(shù)然后根據(jù)校準模型利用主數(shù)據(jù)集,估算出每個個體校準后的XRC1;最后建立Y與XRC1的logistic回歸模型Pr(Y=1|X)=exp(β0+β1XRC1),估算出系數(shù)即為校準后的回歸系數(shù)。此種回歸校準簡稱RC1。

    3.模擬外推(simulation extrapolation,SIMEX)

    模擬外推是一種基于模擬的參數(shù)估計方法,其基本原理是:對有測量誤差的自變量W額外加一個逐級增大的擬測量誤差,基于擬數(shù)據(jù)集進行參數(shù)估計,再擬合這些參數(shù)的變化規(guī)律,最后外推出正確的參數(shù)值。此方法前提條件是W與真實值X滿足經(jīng)典測量誤差模型W=X+U,其中且U與X、Y相互獨立。根據(jù)是否已知,SIMEX可分為參數(shù)型和經(jīng)驗型[10]。參數(shù)型SIMEX(P-SIMEX)的基本步驟如下:

    模擬:在自變量觀測值W原有測量誤差U的基礎(chǔ)上,加上逐級增大的測量誤差,增加的幅度為θ,形成一個新模擬的其中測量誤差U的大小用其方差度量,即的測量誤差為

    擬估計:設(shè)θ為誤差增加幅度,b為擬估計模擬的次數(shù),定義

    數(shù)據(jù)模擬

    1.模擬設(shè)計

    真實數(shù)據(jù)集參數(shù)設(shè)置:樣本量N=1000,自變量X~N(0,1),協(xié)變量Z~N(0,1),通過logistic回歸模型Pr(Y=1|X,Z)=exp(β0+βxX+βzZ),產(chǎn)生因變量Y。β0=-2.5,βX=(0.5,1),βZ=1,分別對應(yīng)有90%、85%的Y=0。

    非差分測量誤差模型參數(shù)設(shè)置:γ0=0,γ=(0.5,根據(jù)模型W=γ0+γX+u產(chǎn)生自變量測量值W。從真實數(shù)據(jù)集中簡單隨機抽取100個觀測值作為效度數(shù)據(jù)集。此時應(yīng)用的校正方法為RC1、P-SIMEX(僅γ0=0,γ=1,且假定已知)。

    為確保模擬結(jié)果的穩(wěn)定性,對上述隨機變量的模擬均采用蒙特卡洛方法[11],模擬次數(shù)為100次,模擬結(jié)果取均值。

    2.評價指標

    評價校正效果的指標為校正后的回歸系數(shù)β*與設(shè)定真實值β真之差的絕對值和β*的標準差。β*與β真之差的絕對值越小,校正效果越好;β*的標準差越小,校正效果越穩(wěn)定。本文擬規(guī)定:0.05]時,校正效果理想;時,校正效果較好。

    結(jié) 果

    1.非差分測量誤差模型校正效果

    不同測量誤差條件下RC1校正效果如表1、圖1所示。將含測量誤差的測量值W作為自變量納入模型時,回歸系數(shù)βnaive都小于β真;隨著測量誤差的增大,βnaive越偏離β真。當非差分測量誤差模型中γ=0.5時,無論測量誤差的大小,經(jīng)RC1校正的雖然較βnaive好,但校正效果都不夠理想;當γ=0.75時,只有測量誤差較小1),RC1才能達到理想的校正效果。僅當W與X的關(guān)系滿足經(jīng)典測量誤差模型(即γ=1),RC1都能達到理想的校正效果。

    基于P-SIMEX的前提條件,僅當γ0=0,γ=1,且假定已知時,應(yīng)用此種方法進行校正。當非差分測量誤差模型中γ=1時,RC1和P-SIMEX如表1、圖2所示。僅當時,P-SIMEX能達到較好的校正效果。隨著測量誤差的增大,P-SIMEX校正的雖然較βnaive大,但校正效果都不好。在相同條件下,隨著的增大,經(jīng)RC1校正的的標準誤逐漸增大,而經(jīng)P-SIMEX校正的的標準誤逐漸減小。

    表1 RC1和P-SIMEX對測量誤差的校正效果

    圖1 非差分測量誤差條件下RC1的校正效果

    圖2 RC1(γ=1)和P-SIMEX的校正效果

    2.經(jīng)典測量誤差模型校正效果

    不同測量誤差條件下RC2和E-SIMEX的校正效果如表2、圖3所示。無論測量誤差多大,經(jīng)RC2校正后校正效果均理想。僅當時,E-SIMEX能達到較好的校正效果。在相同條件下,隨著的增大,經(jīng)RC2校正的的標準誤逐漸增大,而經(jīng)E-SIMEX校正的的標準誤逐漸減小。

    表2 RC2和E-SIMEX對測量誤差的校正效果

    圖3 經(jīng)典測量誤差條件下RC2及E-SIMEX的校正效果

    討 論

    測量誤差在環(huán)境暴露和營養(yǎng)膳食等流行病學(xué)調(diào)查中是普遍存在的,但在大多數(shù)研究數(shù)據(jù)處理分析時,往往被忽視,對其造成的影響很少進行調(diào)整或糾正。直接用含測量誤差的數(shù)據(jù)進行回歸分析時會產(chǎn)生一定的偏倚,導(dǎo)致參數(shù)估計有偏,或者推斷不準確,其效應(yīng)稱為回歸稀釋[12]。根據(jù)測量誤差的來源,可產(chǎn)生不同的稀釋效應(yīng)[13]:自變量X的測量誤差會造成回歸系數(shù)的減小;而應(yīng)變量Y的測量誤差不影響回歸系數(shù),但會造成其標準誤增大,使估計的精度下降。本文通過模擬試驗比較回歸校準法與模擬外推法兩種方法對logistic回歸中自變量測量誤差的校正效果結(jié)果顯示,當自變量有測量誤差時,未校正的βnaive都小于設(shè)定的β真,符合回歸稀釋的效應(yīng)。

    根據(jù)自變量X是否可精確測量,測量誤差的校正可分為RC1、P-SIMEX及RC2、E-SIMEX兩大類。當X可精確測量時,在設(shè)定的條件下,RC1的校正效果較好,P-SIMEX僅當很小時校正效果較好。當X不可精確測量時,隨著測量誤差的增大,E-SIMEX的校正效果降低,而RC2的校正效果相對較穩(wěn)定。無論X是否可測,在經(jīng)典測量誤差模型前提下,RC對logistic回歸模型中測量誤差的校正效果優(yōu)于SIMEX。但SIMEX也有自身的優(yōu)勢,其不僅可獲得校正后的回歸系數(shù),還可通過外推圖示直觀地體現(xiàn)出測量誤差對回歸系數(shù)的效應(yīng)。與RC相比,SIMEX更適用于對非線性回歸模型中的測量誤差進行校正[8]。

    基于本文模擬結(jié)果,對logistic回歸模型中測量誤差進行校正時,在經(jīng)典測量誤差模型前提下,建議應(yīng)用回歸校正法。本研究僅對模型中單個自變量(正態(tài)分布)含加法測量誤差的情況進行討論。但在實際應(yīng)用中,模型中含測量誤差自變量的個數(shù)、自變量的分布以及測量誤差模型的結(jié)構(gòu)更加復(fù)雜,回歸校準與模擬外推的校正效果如何,還需進一步研究。

    1.Pollack AZ,Perkins NJ,Mum ford SL,etal.Correlated biomarkermeasurement error:an important threat to inference in environmental epidem iology.American journal of epidem iology,2013,177(1):84-92.

    2.Day NE,Wong MY,Bingham S,et al.Correlated measurement error implications for nutritional epidem iology.Int J Epidem iol,2004;33(6):1373-1381.

    3.張橋,張秋菊,陳霞,等.回歸校準法對線性回歸中信息偏倚校正的模擬研究.中國衛(wèi)生統(tǒng)計,2013,30(3):327-329.

    4.Rosner B,Spiegelman D,W illettWC.Correction of logistic regression relative risk estimates and confidence intervals for measurement error:the case ofmultiple covariatesmeasured w ith error.American Journal of Epidem iology,1990,132(4):734-745.

    5.Cook JR,Stefanski LA.Simulation-extrapolation estimation in parametric measurement errormodels.Journal of the American Statistical Association,1994,89(428):1314-1328.

    6.Buonaccorsi JP.Measurement error:models,methods and applications. Chapman and Hall/CRC,2010.

    7.Jaceldo-Siegl K,Knutsen SF,Sabate′J,et al.Validation of nutrient intake using an FFQ and repeated 24 h recalls in black and white subjects of the Adventist Health Study-2(AHS-2).Public Health Nutr,2010,13(6):812-819.

    8.Carroll RJ,Ruppert D,Stefanski LA,etal.Measurementerror in nonlinearmodels:amodern perspective.Second Edition,Chapman and Hall/CRC,2006.

    9.Armstrong BG.Effect ofmeasurement error on epidem iological studies of environmental and occupational exposures.Occupational and environmentalmedicine,1998,55(10):651-656.

    10.Devanarayan V,Stefanski LA.Empirical simulation extrapolation for measurement error models with replicate measurements.Statistics&Probability Letters,2002,59(3):219-225.

    11.Rubinstein RY.Simulation and the Monte Carlo method.Wiley-interscience,2009.

    12.MacMahon S,Peto R,Cutler J,et al.Blood pressure,stroke,and coronary heart disease.Part1,prolonged differences in blood pressure:prospective observational studies corrected for the regression dilution bias. Lancet,1990,335:765-774.

    13.Hutcheon JA,Chiolero A,Hanley JA.Random measurement error and regression dilution bias.BMJ,2010,340(7761):1402-1406.

    (責(zé)任編輯:丁海龍)

    The Correction Effect Analysis of Regression Calibration and Simulation Extrapolation for M easurement Error

    Chen Xia,Zhang Qiao,Chen Yongjie,et al(BiostatisticsDepartment,HarbinMedicalUniversity(150081),Harbin)

    ObjectiveTo estimate the correction effect of regression calibration(RC)and simulation extrapolation(SIMEX)formeasurement error in logistic regression model.MethodsWe simulated datasets including given measurement error by SAS software,and applied RC and SIMEX to correctmeasurementerror.Then evaluated the correction effectby comparing the difference between trueβvalue and correctedβ*value.ResultsIfXcould be accuratelymeasured,RC1 could achieve preferable correction effect under the given,while P-SIMEX could only do wellwhen the valuewere ratherm inimal.IfXwas unobservable,the effect of E-SIMEX appeared worse w ith the increase ofmeasurement error,while RC2 performed relatively stable w ith differentvalues.ConclusionWhether or notXcould be accurately measured,RC worked better than SIMEX in the classicalmeasurement errormodel condition.We recommend RC on the basis of our simulation results.

    Measurementerror;Regression calibration;Simulation extrapolation;Logistic regression;Regression dilution

    國家科技支撐計劃(2011BAIO9B02)

    1.哈爾濱醫(yī)科大學(xué)公共衛(wèi)生學(xué)院衛(wèi)生統(tǒng)計學(xué)教研室(150081)

    2.哈爾濱醫(yī)科大學(xué)醫(yī)務(wù)處

    △通信作者:劉美娜,E-mail:liumeina369@163.com

    猜你喜歡
    外推法測量誤差回歸系數(shù)
    一種基于機械臂的外推天線增益標定系統(tǒng)
    中國新通信(2023年3期)2023-06-24 03:00:06
    密度測量誤差分析
    縱向數(shù)據(jù)下變系數(shù)測量誤差模型的漸近估計
    濱海相河口大型沉井基礎(chǔ)設(shè)計荷載效應(yīng)估計
    人民黃河(2018年8期)2018-09-10 16:05:55
    多元線性回歸的估值漂移及其判定方法
    電導(dǎo)法協(xié)同Logistic方程進行6種蘋果砧木抗寒性的比較
    多元線性模型中回歸系數(shù)矩陣的可估函數(shù)和協(xié)方差陣的同時Bayes估計及優(yōu)良性
    用熱線法測定氣體熱導(dǎo)率的儀器改進
    牽引變壓器功率測量誤差分析
    稀溶液與固體反應(yīng)焓變的測量和校準
    18禁动态无遮挡网站| 久久女婷五月综合色啪小说| 久久97久久精品| 亚洲精品色激情综合| 亚洲真实伦在线观看| 夜夜看夜夜爽夜夜摸| 亚洲国产日韩一区二区| 午夜视频国产福利| 美女视频免费永久观看网站| 热99国产精品久久久久久7| 日日摸夜夜添夜夜爱| 爱豆传媒免费全集在线观看| 精品国产乱码久久久久久小说| 这个男人来自地球电影免费观看 | 精品久久久久久久久亚洲| 最近中文字幕2019免费版| av有码第一页| kizo精华| 成人午夜精彩视频在线观看| 在线观看一区二区三区激情| 丝瓜视频免费看黄片| 天堂俺去俺来也www色官网| 亚洲人与动物交配视频| 我要看黄色一级片免费的| 好男人视频免费观看在线| 王馨瑶露胸无遮挡在线观看| 91精品伊人久久大香线蕉| 人人妻人人添人人爽欧美一区卜| 欧美国产精品一级二级三级 | 久热久热在线精品观看| 久久精品国产a三级三级三级| 国产深夜福利视频在线观看| 99热这里只有精品一区| 新久久久久国产一级毛片| 国产日韩欧美亚洲二区| 一个人免费看片子| 成年女人在线观看亚洲视频| 97超碰精品成人国产| 伊人久久精品亚洲午夜| 久久精品国产亚洲av涩爱| 久久亚洲国产成人精品v| 国产精品人妻久久久久久| 中文资源天堂在线| 久久毛片免费看一区二区三区| av专区在线播放| 少妇的逼水好多| 亚洲国产精品999| 欧美97在线视频| a级一级毛片免费在线观看| 日本欧美视频一区| 秋霞伦理黄片| 亚洲精品一二三| 日韩,欧美,国产一区二区三区| 伦精品一区二区三区| 日韩,欧美,国产一区二区三区| 国产白丝娇喘喷水9色精品| 亚洲性久久影院| 国产在线视频一区二区| 大香蕉久久网| 国产高清不卡午夜福利| 国产成人精品福利久久| 欧美最新免费一区二区三区| 精品熟女少妇av免费看| 国产精品一区二区性色av| 一级,二级,三级黄色视频| 国产免费福利视频在线观看| 一级毛片黄色毛片免费观看视频| 高清午夜精品一区二区三区| 精品熟女少妇av免费看| 国内少妇人妻偷人精品xxx网站| 亚洲精品自拍成人| 中文欧美无线码| 精品久久久久久久久亚洲| 亚洲国产av新网站| 一级毛片 在线播放| 久久久久视频综合| 一本色道久久久久久精品综合| 亚州av有码| 在线观看三级黄色| 国产欧美另类精品又又久久亚洲欧美| 在线观看三级黄色| 日日撸夜夜添| 精品午夜福利在线看| av在线观看视频网站免费| 在线看a的网站| 免费看不卡的av| 中国美白少妇内射xxxbb| 国产成人一区二区在线| 男女无遮挡免费网站观看| 亚洲人成网站在线观看播放| 国产伦理片在线播放av一区| 亚洲在久久综合| 国产成人freesex在线| 亚洲激情五月婷婷啪啪| 久久热精品热| 成人国产麻豆网| 又粗又硬又长又爽又黄的视频| 九色成人免费人妻av| 国国产精品蜜臀av免费| 中国国产av一级| 综合色丁香网| 久久国产亚洲av麻豆专区| 卡戴珊不雅视频在线播放| 国产日韩一区二区三区精品不卡 | 少妇高潮的动态图| 看十八女毛片水多多多| 老司机影院毛片| 国产熟女欧美一区二区| 男的添女的下面高潮视频| 亚洲精品日韩在线中文字幕| 只有这里有精品99| 国产精品成人在线| 观看av在线不卡| 久久久久久久久大av| 精品一区二区免费观看| 日本vs欧美在线观看视频 | 亚洲精品一区蜜桃| 中文字幕久久专区| 日韩成人伦理影院| www.av在线官网国产| 色94色欧美一区二区| 三级国产精品欧美在线观看| freevideosex欧美| 日韩伦理黄色片| 免费不卡的大黄色大毛片视频在线观看| 免费久久久久久久精品成人欧美视频 | 伦精品一区二区三区| .国产精品久久| 欧美精品人与动牲交sv欧美| 国产有黄有色有爽视频| www.色视频.com| 亚州av有码| 亚洲色图综合在线观看| 免费看日本二区| 国产毛片在线视频| 亚洲久久久国产精品| 老司机影院毛片| 久久久久久久久久久久大奶| av在线观看视频网站免费| 日本色播在线视频| 日韩大片免费观看网站| 亚洲欧洲精品一区二区精品久久久 | 久久ye,这里只有精品| 18禁裸乳无遮挡动漫免费视频| 国产精品欧美亚洲77777| 黄色毛片三级朝国网站 | 欧美+日韩+精品| h日本视频在线播放| 日本-黄色视频高清免费观看| 欧美亚洲 丝袜 人妻 在线| 亚洲伊人久久精品综合| 有码 亚洲区| 六月丁香七月| 最近的中文字幕免费完整| 丝瓜视频免费看黄片| 蜜桃在线观看..| 国产av精品麻豆| 街头女战士在线观看网站| 99久久精品热视频| 亚洲精品国产成人久久av| 一本久久精品| 欧美日韩亚洲高清精品| 亚洲av不卡在线观看| 夜夜爽夜夜爽视频| 久久久久久久久久成人| 国产av精品麻豆| 91aial.com中文字幕在线观看| 亚洲一级一片aⅴ在线观看| 国产精品女同一区二区软件| av不卡在线播放| 久久精品久久久久久久性| 日韩视频在线欧美| 成人午夜精彩视频在线观看| 国产高清不卡午夜福利| 日韩人妻高清精品专区| 成人国产av品久久久| 成人国产av品久久久| 热99国产精品久久久久久7| 国产日韩欧美视频二区| 国产欧美日韩精品一区二区| 亚洲欧美成人精品一区二区| 黑丝袜美女国产一区| 国产极品天堂在线| 一本—道久久a久久精品蜜桃钙片| 国产综合精华液| 一区二区三区乱码不卡18| 热re99久久精品国产66热6| 国产一区有黄有色的免费视频| 亚洲国产日韩一区二区| 天天操日日干夜夜撸| 在线看a的网站| 欧美日韩在线观看h| 日韩大片免费观看网站| 久久6这里有精品| 国产精品一区www在线观看| 亚洲国产精品国产精品| 精品少妇黑人巨大在线播放| 亚洲经典国产精华液单| 日本与韩国留学比较| 最近最新中文字幕免费大全7| 欧美另类一区| av国产久精品久网站免费入址| 亚洲精品中文字幕在线视频 | 大话2 男鬼变身卡| 婷婷色综合大香蕉| 亚洲欧洲日产国产| 精品卡一卡二卡四卡免费| 久久精品国产亚洲网站| 国产日韩欧美在线精品| 黄色一级大片看看| 亚洲一区二区三区欧美精品| 国内揄拍国产精品人妻在线| av一本久久久久| 中文字幕精品免费在线观看视频 | 少妇人妻 视频| 国产日韩一区二区三区精品不卡 | 日韩强制内射视频| 丰满少妇做爰视频| 国产精品久久久久久久久免| 免费人成在线观看视频色| 老司机影院毛片| 看非洲黑人一级黄片| 久久ye,这里只有精品| 国内精品宾馆在线| 成年av动漫网址| 中文资源天堂在线| 91午夜精品亚洲一区二区三区| 一级片'在线观看视频| 中文字幕精品免费在线观看视频 | 欧美亚洲 丝袜 人妻 在线| 中文精品一卡2卡3卡4更新| 久久 成人 亚洲| 国产精品久久久久久精品古装| 欧美变态另类bdsm刘玥| 国产女主播在线喷水免费视频网站| 亚洲综合色惰| 一二三四中文在线观看免费高清| 国产亚洲一区二区精品| 久热这里只有精品99| 日本-黄色视频高清免费观看| 草草在线视频免费看| 色视频在线一区二区三区| 午夜激情福利司机影院| 中文字幕久久专区| 精品国产一区二区三区久久久樱花| 视频区图区小说| 91久久精品国产一区二区三区| 久久国内精品自在自线图片| 精品人妻一区二区三区麻豆| 看十八女毛片水多多多| 国产精品无大码| 18禁在线播放成人免费| 欧美xxxx性猛交bbbb| h视频一区二区三区| 亚洲国产精品一区二区三区在线| 26uuu在线亚洲综合色| 亚洲精品久久午夜乱码| 午夜激情福利司机影院| 在线观看免费视频网站a站| 国产免费福利视频在线观看| 亚洲av男天堂| 国产精品久久久久久av不卡| 亚洲人成网站在线播| 精品人妻熟女毛片av久久网站| 日日啪夜夜撸| 国产亚洲午夜精品一区二区久久| 热re99久久国产66热| 极品教师在线视频| 麻豆成人av视频| 国产黄片美女视频| 男的添女的下面高潮视频| 大片电影免费在线观看免费| av卡一久久| 亚洲精品乱码久久久v下载方式| 国产毛片在线视频| 亚洲av国产av综合av卡| 精品午夜福利在线看| 亚洲综合精品二区| 亚洲精品国产av蜜桃| 看免费成人av毛片| 精品国产乱码久久久久久小说| 啦啦啦中文免费视频观看日本| 欧美成人精品欧美一级黄| 99国产精品免费福利视频| 少妇人妻一区二区三区视频| 亚洲欧洲精品一区二区精品久久久 | 如日韩欧美国产精品一区二区三区 | 成年女人在线观看亚洲视频| 亚洲成人av在线免费| 成人国产麻豆网| av卡一久久| 老熟女久久久| 99精国产麻豆久久婷婷| 美女国产视频在线观看| 午夜老司机福利剧场| 极品少妇高潮喷水抽搐| 国产av码专区亚洲av| 十分钟在线观看高清视频www | 又黄又爽又刺激的免费视频.| 久久国产精品大桥未久av | 久热久热在线精品观看| 一个人看视频在线观看www免费| 日日啪夜夜撸| 丰满迷人的少妇在线观看| 美女视频免费永久观看网站| 人妻制服诱惑在线中文字幕| 日本色播在线视频| 91久久精品国产一区二区成人| 国产深夜福利视频在线观看| 色网站视频免费| 色婷婷av一区二区三区视频| 国语对白做爰xxxⅹ性视频网站| 亚洲第一av免费看| 免费人妻精品一区二区三区视频| 91久久精品国产一区二区成人| 久久久久久久亚洲中文字幕| 亚洲国产精品一区二区三区在线| 成人18禁高潮啪啪吃奶动态图 | 美女福利国产在线| 中文字幕制服av| 国产男女超爽视频在线观看| 99久久中文字幕三级久久日本| 久久国产亚洲av麻豆专区| 国产精品伦人一区二区| 精品酒店卫生间| 精品国产一区二区久久| 午夜91福利影院| 国产高清三级在线| 国产成人精品久久久久久| 丰满人妻一区二区三区视频av| 男人和女人高潮做爰伦理| 久久婷婷青草| 国产黄色免费在线视频| 午夜免费鲁丝| 久久精品久久久久久久性| 99re6热这里在线精品视频| 少妇人妻 视频| 夜夜骑夜夜射夜夜干| 一级爰片在线观看| 亚洲精品乱码久久久久久按摩| 中国国产av一级| 中文字幕亚洲精品专区| 国产白丝娇喘喷水9色精品| 国产伦精品一区二区三区视频9| 热re99久久国产66热| 亚洲图色成人| 免费观看av网站的网址| 国产精品免费大片| 国产亚洲一区二区精品| 伦理电影大哥的女人| xxx大片免费视频| 99热6这里只有精品| 欧美激情国产日韩精品一区| 视频区图区小说| 99久国产av精品国产电影| 亚洲精品一二三| 久久久久精品久久久久真实原创| 国模一区二区三区四区视频| 内地一区二区视频在线| av有码第一页| 日韩大片免费观看网站| 国产精品嫩草影院av在线观看| 久热这里只有精品99| 国产精品偷伦视频观看了| 欧美性感艳星| 中文字幕亚洲精品专区| 欧美区成人在线视频| 啦啦啦在线观看免费高清www| 亚洲精品成人av观看孕妇| 日韩,欧美,国产一区二区三区| 又大又黄又爽视频免费| videossex国产| 成年女人在线观看亚洲视频| 中文欧美无线码| 在线天堂最新版资源| 国产一区二区在线观看日韩| 亚洲va在线va天堂va国产| 亚洲国产毛片av蜜桃av| 久久久久久久久久久久大奶| 美女xxoo啪啪120秒动态图| 亚洲精品国产av蜜桃| 人体艺术视频欧美日本| 国产精品欧美亚洲77777| 国产精品蜜桃在线观看| 一级毛片黄色毛片免费观看视频| 国产精品免费大片| 色视频www国产| 色网站视频免费| 97超视频在线观看视频| 亚洲国产最新在线播放| h视频一区二区三区| 日本与韩国留学比较| 一区二区三区免费毛片| a级一级毛片免费在线观看| 亚洲精品日韩av片在线观看| 亚洲婷婷狠狠爱综合网| 欧美97在线视频| 欧美亚洲 丝袜 人妻 在线| av国产精品久久久久影院| 欧美97在线视频| 久久国产乱子免费精品| 欧美 亚洲 国产 日韩一| 久久国产精品男人的天堂亚洲 | 最近的中文字幕免费完整| 伦理电影大哥的女人| 有码 亚洲区| 欧美日韩亚洲高清精品| 久久99一区二区三区| 少妇裸体淫交视频免费看高清| 97在线人人人人妻| 两个人的视频大全免费| 视频中文字幕在线观看| 一个人免费看片子| 如日韩欧美国产精品一区二区三区 | 成人黄色视频免费在线看| 观看美女的网站| 久久这里有精品视频免费| 女人精品久久久久毛片| 亚洲av不卡在线观看| 2022亚洲国产成人精品| 在现免费观看毛片| 色网站视频免费| 国产成人免费无遮挡视频| 亚洲性久久影院| 亚洲欧洲精品一区二区精品久久久 | 精品国产乱码久久久久久小说| 黄色毛片三级朝国网站 | 99热网站在线观看| 国产中年淑女户外野战色| 久久久久久久久久久免费av| 免费久久久久久久精品成人欧美视频 | 大码成人一级视频| 亚洲丝袜综合中文字幕| 成人毛片a级毛片在线播放| 最近中文字幕高清免费大全6| 老女人水多毛片| 亚洲精品亚洲一区二区| 在线 av 中文字幕| 国产av精品麻豆| 亚洲精品久久午夜乱码| 国产日韩欧美在线精品| 蜜臀久久99精品久久宅男| 亚洲成人手机| 中文字幕久久专区| 国产精品偷伦视频观看了| 欧美日韩国产mv在线观看视频| 赤兔流量卡办理| av福利片在线| 五月天丁香电影| 亚洲性久久影院| 国产亚洲欧美精品永久| 日韩 亚洲 欧美在线| 亚洲欧美精品自产自拍| 一级毛片久久久久久久久女| 男女无遮挡免费网站观看| 亚洲美女黄色视频免费看| 欧美日韩av久久| 高清午夜精品一区二区三区| 卡戴珊不雅视频在线播放| 又黄又爽又刺激的免费视频.| 亚洲精品日韩av片在线观看| 美女内射精品一级片tv| 成年女人在线观看亚洲视频| 97超视频在线观看视频| 中文字幕免费在线视频6| 欧美日韩一区二区视频在线观看视频在线| 18禁在线播放成人免费| 国产精品麻豆人妻色哟哟久久| 一区二区av电影网| 久久99热6这里只有精品| 伊人亚洲综合成人网| 亚洲成人手机| 国产成人午夜福利电影在线观看| 九色成人免费人妻av| 涩涩av久久男人的天堂| 国产极品粉嫩免费观看在线 | 黄色怎么调成土黄色| 亚洲人成网站在线播| 只有这里有精品99| 亚洲欧美一区二区三区国产| 91久久精品电影网| 色婷婷久久久亚洲欧美| 精品国产露脸久久av麻豆| 国产成人免费无遮挡视频| 精品少妇黑人巨大在线播放| 国产精品欧美亚洲77777| 国产免费又黄又爽又色| 亚州av有码| 99久久人妻综合| 久久精品久久久久久噜噜老黄| 亚洲欧美中文字幕日韩二区| 精品酒店卫生间| 久久久精品94久久精品| 青春草视频在线免费观看| 成人国产av品久久久| 久久久午夜欧美精品| 色婷婷久久久亚洲欧美| 国产白丝娇喘喷水9色精品| 久久 成人 亚洲| 久久鲁丝午夜福利片| 午夜免费鲁丝| 亚洲,欧美,日韩| 久久午夜综合久久蜜桃| 久久久精品免费免费高清| 久久人人爽人人爽人人片va| 国产免费又黄又爽又色| 国产极品粉嫩免费观看在线 | 一级毛片我不卡| 成人亚洲精品一区在线观看| 久久影院123| 另类亚洲欧美激情| 亚洲内射少妇av| 精品视频人人做人人爽| 精品一区二区免费观看| 久久久久久久国产电影| 欧美日韩精品成人综合77777| 啦啦啦啦在线视频资源| 日日啪夜夜撸| 少妇熟女欧美另类| 蜜臀久久99精品久久宅男| 韩国高清视频一区二区三区| 只有这里有精品99| 久久99热6这里只有精品| 22中文网久久字幕| 亚洲人与动物交配视频| 国产精品久久久久久久电影| 能在线免费看毛片的网站| 久久国产精品大桥未久av | 波野结衣二区三区在线| 久久 成人 亚洲| 国产精品.久久久| 午夜老司机福利剧场| 国产成人freesex在线| 18禁在线播放成人免费| 不卡视频在线观看欧美| 卡戴珊不雅视频在线播放| 欧美最新免费一区二区三区| 亚洲高清免费不卡视频| av天堂中文字幕网| 插阴视频在线观看视频| 亚洲欧美一区二区三区黑人 | 一个人免费看片子| 美女国产视频在线观看| videos熟女内射| 日韩中文字幕视频在线看片| 新久久久久国产一级毛片| 欧美精品国产亚洲| 大又大粗又爽又黄少妇毛片口| 大码成人一级视频| 国产欧美日韩精品一区二区| a级毛片免费高清观看在线播放| 97超碰精品成人国产| 偷拍熟女少妇极品色| 少妇被粗大的猛进出69影院 | 秋霞在线观看毛片| 永久免费av网站大全| 免费黄色在线免费观看| 成年人免费黄色播放视频 | 免费在线观看成人毛片| 国产精品成人在线| 中文字幕亚洲精品专区| 极品人妻少妇av视频| 老女人水多毛片| 国产男人的电影天堂91| a级一级毛片免费在线观看| 91久久精品国产一区二区三区| 91精品一卡2卡3卡4卡| 亚洲情色 制服丝袜| 日韩,欧美,国产一区二区三区| 91精品国产国语对白视频| 亚洲av二区三区四区| 午夜福利影视在线免费观看| 永久免费av网站大全| 日韩一区二区三区影片| 美女主播在线视频| 中文字幕免费在线视频6| 成人国产av品久久久| 蜜臀久久99精品久久宅男| 国产精品一二三区在线看| 欧美日韩精品成人综合77777| 国产爽快片一区二区三区| 国产一区二区在线观看av| 91成人精品电影| 国产精品免费大片| 少妇的逼水好多| 国产综合精华液| h日本视频在线播放| 日韩欧美精品免费久久| 亚洲成人av在线免费| 亚洲四区av| 欧美精品人与动牲交sv欧美| a级毛色黄片| 国产 一区精品| 亚洲人成网站在线播| 国产极品粉嫩免费观看在线 | 青春草国产在线视频| 一区二区三区精品91| 欧美丝袜亚洲另类| 欧美日韩精品成人综合77777| 中文字幕制服av| 亚洲欧美成人综合另类久久久| 三上悠亚av全集在线观看 | 极品少妇高潮喷水抽搐| 有码 亚洲区| 不卡视频在线观看欧美| 一本久久精品| 欧美精品一区二区大全| 日韩一区二区三区影片| 看非洲黑人一级黄片| 免费看av在线观看网站| 国产欧美日韩综合在线一区二区 | a级一级毛片免费在线观看| 国产成人91sexporn| 高清av免费在线| av视频免费观看在线观看|