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

    隨機(jī)回歸模型在畜禽育種中的應(yīng)用

    2021-12-03 08:20:09邢杰
    中國畜禽種業(yè) 2021年8期
    關(guān)鍵詞:協(xié)方差方差遺傳

    邢杰

    (山西省右玉縣畜牧獸醫(yī)中心 037200)

    畜牧業(yè)生產(chǎn)中,遺傳改良是畜禽育種工作的重要環(huán)節(jié),利用統(tǒng)計模型進(jìn)行遺傳性狀評估是通過估計方差組分和遺傳參數(shù),采用合適的分析模型,在眾多動物群體中借助混合線性技術(shù),在較短時間內(nèi)篩選出優(yōu)良的目標(biāo)性狀,并加以固定。作為一種數(shù)據(jù)分析模型,隨機(jī)回歸模型(RRM) 有著廣泛的用途,應(yīng)用于社會經(jīng)濟(jì)諸多領(lǐng)域。早期RRM 主要應(yīng)用于奶牛泌乳、生長等性狀的預(yù)測評估,Meyer 等證明,協(xié)方差函數(shù)和RRM是等價的,并將其用于分析縱向數(shù)據(jù)或重復(fù)記錄,用來解釋由動態(tài)性狀在不同個體間動態(tài)變化而導(dǎo)致的遺傳差異[1]。隨著此模型的普及與優(yōu)化,人們逐漸將其用于肉牛、羊、豬及禽類(以雞為主)等其他畜禽相關(guān)性狀的遺傳評估,從對單一或少數(shù)性狀轉(zhuǎn)變到多個不同性狀的評估分析,這也是RRM 的不斷完善與推廣的漫長過程。

    1 隨機(jī)回歸模型的理論基礎(chǔ)

    1.1 回歸分析思想的建立與應(yīng)用

    回歸分析是一種預(yù)測性的建模技術(shù),它研究的是因變量和自變量之間的關(guān)系,通常用于預(yù)測分析時間序列模型及發(fā)現(xiàn)變量之間的因果關(guān)系。如通過對奶牛泌乳期各個階段產(chǎn)奶量與泌乳時間、環(huán)境等相關(guān)性的分析來確定產(chǎn)奶性狀的優(yōu)劣就是回歸分析方法的一種運(yùn)用?;貧w分析是建模與分析數(shù)據(jù)的重要手段,利用它能找出自變量與因變量之間的顯著關(guān)系,并顯示多個自變量對一個因變量的影響強(qiáng)度。在7 種回歸分析技術(shù)中,線性回歸和多項(xiàng)式回歸是常用于畜禽遺傳評估的兩種回歸分析技術(shù)。線性回歸的因變量是連續(xù)的,自變量可以是連續(xù)的,也可以是離散的,線性回歸使用最佳擬合直線建立因變量與自變量之間的關(guān)系,一元與多元線性回歸的區(qū)別在于,前者只有一個自變量,而后者有多個。多項(xiàng)式回歸則是一種運(yùn)用自變量指數(shù)大于1 的回歸方程,其最佳擬合線不是直線,而是用于擬合數(shù)據(jù)點(diǎn)的曲線。

    1.2 隨機(jī)回歸模型的提出

    隨機(jī)回歸模型(RRM)分析法最早由Henderson 等(1982年)提出[2],隨后Schaeffer 和Dekkers(1994 年)提議將此方法用于動物育種研究領(lǐng)域[3],Meyer 和Hill (1997 年) 證明,RRM 和協(xié)方差函數(shù)是等價的,并將其應(yīng)用于分析縱向數(shù)據(jù)或重復(fù)記錄,即用于分析動物生命過程中的一定時期收集的有關(guān)性狀記錄[1]。因此,Meyer 和Hill 的研究是對RRM 應(yīng)用的早期完善。

    1.3 隨機(jī)回歸模型的引入

    沿軌跡進(jìn)行的測定通??山槎x該軌跡參數(shù)的函數(shù)。重復(fù)測定經(jīng)常被人們用于描述性狀隨時間的變化,而處理重復(fù)測定的模型常常被用于動物生產(chǎn),如對奶牛多個泌乳記錄的分析就是經(jīng)常使用“重復(fù)模型”,從遺傳學(xué)的角度看,這種模型的典型特征是重復(fù)記錄被認(rèn)為是對同一性狀的表達(dá),即重復(fù)泌乳記錄之間的相關(guān)性被認(rèn)為是統(tǒng)一的。這類模型能專門描述按時間回歸的變化(如泌乳曲線和生長曲線),對這種變化的分析可能揭示其影響變量產(chǎn)生的原因,而相應(yīng)參數(shù)在品種內(nèi)也可能存在加性遺傳變異。分析途徑是估計每個動物的曲線參數(shù)和確定此類參數(shù)的變量組成,也可以使用模型來分析回歸系數(shù)在動物間改變,這樣的回歸系數(shù)不固定,但允許按指定分布而改變,從而表現(xiàn)為隨機(jī)回歸系數(shù)[4]。

    1.4 隨機(jī)回歸模型的分類

    1.4.1 基于重復(fù)率模型的隨機(jī)回歸模型

    這種模型有多種形式,Schaeffer 提出的RRM 是較為典型的一種,其表達(dá)式如下:

    其中y 為測定日產(chǎn)奶性狀觀察值,x 為對應(yīng)的協(xié)變量,HTD 為牛群的測定日固定環(huán)境效應(yīng);a,p 分別是加性遺傳效應(yīng)和隨機(jī)永久環(huán)境效應(yīng)的隨機(jī)回歸系數(shù),∑ajxj和∑pkxk分別為隨機(jī)回歸部分;b 為回歸系數(shù),∑bixi為固定回歸部分,用來擬合群的泌乳曲線;e 為隨機(jī)殘差。

    基于重復(fù)率模型的常用于分析奶牛的產(chǎn)奶記錄(泌乳曲線),另外Schaeffer 等還提出使用此模型對動物生存性狀進(jìn)行分析[5,6]。

    1.4.2 基于母體效應(yīng)

    此模型將基于重復(fù)率模型的RRM 的加性遺傳效應(yīng)分為直接遺傳效應(yīng)和母體遺傳效應(yīng),基于母體效應(yīng)的是基于重復(fù)率模型的一種,對應(yīng)母畜,還有一種基于公畜模型的RRM,其性質(zhì)與此模型類似[6]。

    1.4.3 二維隨機(jī)回歸模型

    將兩個隨時間改變的性狀對應(yīng)的時間尺度同時作為分析性狀測定的時間尺度,此模型即可稱為二維RRM。如使用此模型分析泌乳期記錄時,就需要將泌乳天數(shù)和胎次同時作為時間尺度。二維RRM 在早期應(yīng)用過程中,有時人們并不考慮不同性狀的時間尺度之間的相互作用,但實(shí)際上它們之間的確存在相關(guān)性,從而導(dǎo)致該模型的估計值不夠準(zhǔn)確[6]。

    1.5 隨機(jī)回歸模型經(jīng)常使用的統(tǒng)計理論

    1.5.1 協(xié)方差函數(shù)的應(yīng)用

    協(xié)方差函數(shù)可定義為給出軌跡不同點(diǎn)測定性狀的方差和協(xié)方差的連續(xù)函數(shù),此函數(shù)可用于描述表型協(xié)方差結(jié)構(gòu)。在RRM 中,使用協(xié)方差函數(shù)的優(yōu)點(diǎn)類似于使用回歸,因?yàn)閰f(xié)方差函數(shù)允許:(1)協(xié)方差隨時間逐漸變化;(2)協(xié)方差和時間差異之間存在關(guān)聯(lián)性;(3)預(yù)測軌跡點(diǎn)的方差和協(xié)方差函數(shù),或即使觀測點(diǎn)很少,但能使用來自所有其他測定的信息。協(xié)方差函數(shù)給出我們?nèi)绾斡米兓膮?shù)恰當(dāng)描述隨時間不斷改變的性狀,這一點(diǎn)對含有遺傳參數(shù)的評估模型尤為重要,協(xié)方差函數(shù)給出了我們處理此類變化的信息,即我們到底如何處理和分析生長曲線和泌乳曲線,協(xié)方差函數(shù)為分析相應(yīng)模型中的獨(dú)立方差分量提供了方法,這種方差分量揭示了性狀隨時間改變的特有模式[4]。

    1.5.2 方差分量(或遺傳參數(shù))的估計

    1.5.2.1 基于似然函數(shù)的估計

    當(dāng)似然函數(shù)取得最大值時,說明采樣數(shù)據(jù)發(fā)生的概率最大,我們稱其為最優(yōu)(或最大)似然函數(shù)。最大似然法(ML)、約束最大似然法(REML),ML 法和REML 法都是以使觀察值或觀察值的線性組合的似然函數(shù)達(dá)到極值的原則來估計方差組合。根據(jù)迭代法的不同,REML 又可分為最大期望約束最大似然法(EMREML)、非求導(dǎo)約束最大似然法(DFREML)和平均信息約束最大似然法(AIREML)[7]。

    1.5.2.2 基于線性模型的估計

    按照最佳線性無偏估計的原則去估計線性模型的固定效應(yīng)和隨機(jī)效應(yīng)的方差分量估計法稱為最佳線性無偏估計法(BLUP),該法使用的數(shù)學(xué)模型是混合線性效應(yīng)模型,因此,它將選擇指數(shù)法和最小二乘法有機(jī)結(jié)合起來。除了BLUP 法外,最小范數(shù)二次無偏估計法(MINQUE)和最小方差二次無偏估計法(MIVQUE)也是兩種方差分量的無偏估計法,但不如BLUP 法常用,這兩種方法均適用于所有混合模型。另外,最小二乘法也是一種模型曲線的擬合方法,最早由法國勒讓德提出,常用于估計方差分量[8]。

    1.5.2.3 貝葉斯法(Bayesian approach)估計

    貝葉斯法是一種重要的方差分量估計方法,此法將當(dāng)前信息和先驗(yàn)信息結(jié)合起來,先驗(yàn)密度與似然函數(shù)的乘積等于后驗(yàn)密度,此法能對數(shù)量性狀和閾性性狀進(jìn)行方差組分估計[8]。

    1.5.3 兩種常用的回歸方法

    1.5.3.1 勒讓德多項(xiàng)式(Legendre polynomials)回歸

    Legendre 多項(xiàng)式回歸是描述兩個變量在[-1,1]上正交的序列函數(shù)(或多項(xiàng)式集合),因變量通常是隨機(jī)的且服從某種分布,自變量通常由觀測者決定的一些時間點(diǎn)來校正,模型可以選擇與時間點(diǎn)個數(shù)相同的多項(xiàng)式階數(shù),這樣可以準(zhǔn)確描述數(shù)據(jù),雖然缺少一般性,但此類函數(shù)的優(yōu)點(diǎn)在于具有正交性,由于存在正交條件高階項(xiàng)系數(shù)趨于零,并且增加和刪除一個項(xiàng)對其他項(xiàng)沒有影響。在生物學(xué)中,絕大部分曲線都可使用三階、四階或五階多項(xiàng)式來表示。

    1.5.3.2 樣條回歸

    樣條回歸是一種特殊的函數(shù),由多項(xiàng)式分段定義,一般的樣條是自然的三次樣條。樣條回歸可以看作一個逐段回歸或分段回歸,即把數(shù)據(jù)分為幾段,在一段內(nèi)分別擬合模型,每一段內(nèi)可以擬合直線、二次項(xiàng)、三次項(xiàng)等,根據(jù)實(shí)際情況而定,與多項(xiàng)式回歸相比,樣條回歸擬合效果更優(yōu),而擬合模型的最終目的都是尋找相對較優(yōu)的模型。在特定回歸模型(如給定的RRM)中,經(jīng)常將樣條回歸用于擬合協(xié)方差結(jié)構(gòu)的分析中。

    2 隨機(jī)回歸模型的優(yōu)點(diǎn)與缺陷

    2.1 優(yōu)點(diǎn)方面

    當(dāng)RRM 側(cè)重于分析測定日記錄或縱向數(shù)據(jù)時,其優(yōu)點(diǎn)在于能以較少參數(shù)準(zhǔn)確擬合性狀的協(xié)方差,能靈活地與生產(chǎn)測定計劃相結(jié)合,能較準(zhǔn)確地消除測定日記錄中的環(huán)境效應(yīng),能用于擬合個體泌乳曲線與群體一般性狀產(chǎn)生的偏差,提高奶牛305d 產(chǎn)奶量的遺傳評估準(zhǔn)確度[6],RRM 的優(yōu)點(diǎn)還在于可利用特征值系數(shù)有效改變性狀曲線,能對不完整記錄進(jìn)行處理[9]。

    2.2 缺陷方面

    2.2.1 表現(xiàn)在等價協(xié)方差函數(shù)的使用上

    給定時間內(nèi)的多次測定值可視為多重性狀,但在時間推移中,測定值可能時間分布,理論上,多性狀方法包含無數(shù)個性狀,而實(shí)際上很多性狀的數(shù)據(jù)都可能會丟失(除非我們進(jìn)行連續(xù)測定),采取簡化的方式就是將特定時間的表達(dá)定義為獨(dú)立性狀,從而達(dá)到限制性狀數(shù)量的目的,但卻存在這些缺陷:(1)我們將協(xié)方差結(jié)構(gòu)擬合為不連續(xù)的結(jié)構(gòu),但實(shí)際上它是連續(xù)的;(2)在正確解釋同一時間段的更多測定值時,這種做法更顯單調(diào)煩瑣;(3)針對許多不同年齡性狀的多性狀模型,其缺陷在于相關(guān)性矩陣不具有結(jié)構(gòu)性[4]。

    2.2.2 表現(xiàn)在數(shù)據(jù)分析的便利程度上

    RRM 本質(zhì)上是一種帶有隨機(jī)效應(yīng)的回歸模型,一般來說,隨機(jī)效應(yīng)模型得出的結(jié)論偏向保守,置信區(qū)間較大,更難發(fā)現(xiàn)差異,帶給我們的信息是,如果各個試驗(yàn)的結(jié)果差異很大時,是否需要把各個試驗(yàn)合并往往需要我們慎重考慮,因而難以很快獲得較為客觀的結(jié)論。

    2.2.3 表現(xiàn)在使用隨機(jī)回歸模型需要的計算量上

    在進(jìn)行多元隨機(jī)回歸分析時,常出現(xiàn)計算量大幅上升現(xiàn)象,而擬合種群或亞種群平均增長曲線的參數(shù)也可以用于RRM,這同樣適合于擬合回歸模型所需的參數(shù)數(shù)量,但這種情況下,對于待分析(大型)動物育種數(shù)據(jù),更多的隨機(jī)回歸系數(shù)會使計算量迅速增加[4],這都是我們不愿面對的。另外,在使用基于RRM 的單步最佳線性無偏估計法(ssBLUP)時,其主要限制因素仍然是分析大數(shù)據(jù)集所需要的計算時間[10]。

    3 隨機(jī)回歸模型的優(yōu)化

    3.1 優(yōu)化手段

    從Legendre 多項(xiàng)式階此組合的優(yōu)選分析入手,尋找該模型的隨機(jī)效應(yīng)對性狀曲線(如生長曲線和泌乳曲線)的最佳擬合點(diǎn)。

    3.2 優(yōu)化目標(biāo)

    將更多的隨機(jī)回歸效應(yīng)納入至隨機(jī)回歸模型中,以期達(dá)到更加準(zhǔn)確的評估分析結(jié)果。

    4 隨機(jī)回歸模型在畜禽育種中應(yīng)用舉例

    4.1 隨機(jī)回歸模型在奶牛育種中的應(yīng)用

    4.1.1 隨機(jī)回歸—測定日模型的應(yīng)用

    奶牛育種的主要目的在于選擇不斷提高群體產(chǎn)奶性狀的遺傳水平[5],不過這也是奶牛遺傳育種的難點(diǎn)之一,產(chǎn)奶性狀極易受其他數(shù)量性狀和環(huán)境的共同影響,且環(huán)境效應(yīng)影響很大。性狀選擇遺傳評估主要是通過估計育種值(EBV)實(shí)現(xiàn)的,可以使用多性狀回歸測定日模型(MMRRM-TDM)進(jìn)行EBV 評估,MMRRM-TDM 利用來自測定日記錄的信息對個體進(jìn)行EBV 估計,進(jìn)而利用不同定義階段的EBV 值對種公畜進(jìn)行遺傳評定和性狀分析,通過持續(xù)力EBV 對泌乳曲線進(jìn)行遺傳改進(jìn)[10]。

    4.1.2 基于線性樣條的多性狀隨機(jī)回歸模型的應(yīng)用

    在不同遺傳組成的種群中使用不同的選擇強(qiáng)度和特定的選擇過程,純種和雜交奶牛之間的差異和協(xié)方差可能發(fā)生在多個品種群體中,制定和實(shí)施考慮基因型手段、親屬間遺傳力協(xié)方差的遺傳評估及每個不同品種群體存在的特定遺傳方差,有助于多品種群體的基因改良。線性樣條RRM(RRMLS)適用于獲得多品種群體的遺傳參數(shù),不同品種組成的后代群與種公牛的BV 會發(fā)生交互作用[11]。

    4.2 隨機(jī)回歸模型在肉牛育種中的應(yīng)用

    國外專家曾嘗試對該國的Hays Converter(HC)牛進(jìn)行遺傳參數(shù)和遺傳趨勢的評估,他們聯(lián)合應(yīng)用了MMT 和RRM 這兩種評估模型。研究顯示,外觀檢查和使用作為候選肉牛的個體表型可能不會對肉牛的性狀選擇產(chǎn)生最大限度的影響,使用估計育種值(EBV)作為選擇基礎(chǔ),預(yù)計將會更快地產(chǎn)生改良基因。該遺傳評估引入了當(dāng)代群體(CG)這個效應(yīng)值,該效應(yīng)在RRM 中被認(rèn)為是隨機(jī)的,則預(yù)測估計方差會減小,因?yàn)槟切└嘁咽占降男畔⒈挥糜陬A(yù)測EBV[12]。

    4.3 隨機(jī)回歸模型在蛋雞育種中的應(yīng)用

    體重也是雞體健康的標(biāo)志之一,與肉雞及其他家禽相比,從目前已取得的研究成果看,國內(nèi)專家曾以白來航雞與東鄉(xiāng)綠殼蛋雞F2資源群為研究對象,根據(jù)已收集的數(shù)據(jù),對這兩種蛋雞的遺傳參數(shù)作出評估。研究顯示,RRM 可用于蛋雞早期體重選育,但RRM 的應(yīng)用和優(yōu)化應(yīng)緊扣研究目的,且高階多項(xiàng)式適用于估計遺傳參數(shù),而低階多項(xiàng)式則更適用于EBV[9]。

    4.4 隨機(jī)回歸模型在綿羊育種中的應(yīng)用

    在家畜生長性狀中,體重是一個重要且容易測量的經(jīng)濟(jì)性狀,家畜生長發(fā)育過程中各階段的體重(包括初生重、斷奶重和周歲重等)是生長性狀遺傳力評估模型中的常用方差分量。國內(nèi)專家曾以考力代綿羊作為試驗(yàn)對象,對其體重生長的遺傳參數(shù)和協(xié)方差函數(shù)進(jìn)行分析和評估,他們根據(jù)似然比檢驗(yàn)(LRT)篩選出用于考力代綿羊多性狀遺傳分析和協(xié)方差函數(shù)估計的最佳模型[13]。

    5 隨機(jī)回歸模型應(yīng)用前景展望

    RRM 是處理包含隨機(jī)效應(yīng)的回歸分析方法,與具有等價作用的測定日模型(TDM)、協(xié)方差模型一樣,都是分析縱向數(shù)據(jù)的常用模型,但對特定時間點(diǎn)的縱向數(shù)據(jù)(針對多性狀模型)分析而言,很明顯只需要考慮性狀間隨機(jī)效應(yīng)(即加性遺傳效應(yīng)、永久環(huán)境效應(yīng)等)的相關(guān)性,因而,隨機(jī)參數(shù)數(shù)量有限,但在實(shí)際應(yīng)用中,這種情況幾乎不存在,因?yàn)椴豢疾煨誀钛貢r間改變的趨勢和走向,對性狀的預(yù)測和估計便無從談起,因而,重復(fù)測定是人們獲得可靠評估分析的必需,但卻是最簡單直接的手段,但一定限度內(nèi)的測定次數(shù),無論從計算量還是分析難度方面講,RRM 都能發(fā)揮其處理大樣本且結(jié)果準(zhǔn)確的優(yōu)點(diǎn),同時從建模計算的發(fā)展前景來看,借助計算機(jī)軟件處理復(fù)雜計算似乎早已成為趨勢,至少計算量大的問題一定限度內(nèi)是可以克服的,但隨著重復(fù)測定次數(shù)的增加,帶來的不僅是計算量的增加,而是性狀測定頻率難以把握和數(shù)據(jù)結(jié)果分析的復(fù)雜程度,因?yàn)橹貜?fù)測定的目的之一就是考察同一性狀在不同時間點(diǎn)的相關(guān)性高低,這是影響遺傳性狀評估準(zhǔn)確性的重要因素,如果時間間隔過大,相關(guān)性可能為零,過小一則測定工作量太大,再則最終產(chǎn)生的隨機(jī)效應(yīng)參數(shù)太多,從計算量上和數(shù)據(jù)分析上(每一協(xié)方差參數(shù)都要作出估計,都要與目標(biāo)曲線進(jìn)行擬合)都是巨大的挑戰(zhàn)。另外,重復(fù)測定的均勻程度也間接關(guān)系到最終結(jié)果的準(zhǔn)確性。另有研究顯示,RRM 對稀疏數(shù)據(jù)很敏感[12],這又間接說明,隨著重復(fù)測定次數(shù)的增多和時間間隔的縮小,RRM 對數(shù)量激增的隨機(jī)效應(yīng)參數(shù)的敏感度變小,這種情況同樣會影響最終分析結(jié)果,因而使用此模型可能難以獲得十分可靠的評估結(jié)果。對于RRM 的應(yīng)用前景展望,很多人都把此模型能處理更多的隨機(jī)效應(yīng)作為其準(zhǔn)確分析試驗(yàn)效應(yīng)的重要條件,但很少有人探討其具體的實(shí)現(xiàn)途徑,從上述分析可以看出,不斷優(yōu)化模型結(jié)構(gòu),尋求更為便捷、快速的計算處理手段,努力探索增加隨機(jī)效應(yīng),是提高評估準(zhǔn)確性與參數(shù)增多帶來分析不便之間的最佳解決方案(如使用樣條函數(shù)擬合協(xié)方差結(jié)構(gòu)就是一種有益嘗試),這幾點(diǎn)應(yīng)該是RRM 在未來發(fā)揮其應(yīng)有作用的關(guān)鍵所在。

    6 結(jié)語

    將數(shù)理統(tǒng)計學(xué)模型應(yīng)用于畜禽育種研究工作中,對現(xiàn)代畜牧業(yè)發(fā)展具有里程碑意義。無論哪個畜種,在具體育種工作中使用何種模型,應(yīng)視其研究目的和方向的不同靈活應(yīng)用,揚(yáng)長避短,只有這樣才能發(fā)揮其應(yīng)有作用,人們對隨機(jī)回歸模型的不斷優(yōu)化及其應(yīng)用畜種的不斷增多充分體現(xiàn)了這一點(diǎn),而隨著計算處理?xiàng)l件和分析手段的進(jìn)一步提高,相信該模型會受到業(yè)內(nèi)人士更多的關(guān)注。

    猜你喜歡
    協(xié)方差方差遺傳
    方差怎么算
    非遺傳承
    概率與統(tǒng)計(2)——離散型隨機(jī)變量的期望與方差
    還有什么會遺傳?
    還有什么會遺傳
    還有什么會遺傳?
    計算方差用哪個公式
    方差生活秀
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識方法
    av.在线天堂| av福利片在线观看| 男人狂女人下面高潮的视频| 亚洲欧美精品综合久久99| 嫩草影院新地址| 亚洲av男天堂| 天天躁夜夜躁狠狠久久av| 22中文网久久字幕| 亚洲欧美日韩东京热| 亚洲色图av天堂| 国产探花在线观看一区二区| 99热6这里只有精品| 国产精品女同一区二区软件| 女人十人毛片免费观看3o分钟| av天堂中文字幕网| 身体一侧抽搐| 波野结衣二区三区在线| 伦理电影大哥的女人| 亚洲欧美成人精品一区二区| 精品人妻一区二区三区麻豆| 亚洲精品国产av成人精品| 国模一区二区三区四区视频| 亚洲精品国产成人久久av| 亚洲国产欧洲综合997久久,| 亚洲精品国产成人久久av| 国产精品日韩av在线免费观看| www.色视频.com| 久久久久久久久久久免费av| 联通29元200g的流量卡| 国产午夜精品论理片| 超碰av人人做人人爽久久| 亚洲精品色激情综合| 久久精品影院6| 国产精品一区二区在线观看99 | 亚洲人成网站在线观看播放| 男人和女人高潮做爰伦理| 在线观看美女被高潮喷水网站| 国产亚洲5aaaaa淫片| 亚洲av一区综合| 99国产极品粉嫩在线观看| 成人亚洲欧美一区二区av| 美女脱内裤让男人舔精品视频 | 欧美潮喷喷水| 一夜夜www| 亚洲av熟女| 美女国产视频在线观看| 久久精品综合一区二区三区| 国产精品一区二区三区四区免费观看| 国产午夜精品论理片| 美女cb高潮喷水在线观看| 熟女人妻精品中文字幕| 嘟嘟电影网在线观看| 成人亚洲精品av一区二区| 99久国产av精品| 非洲黑人性xxxx精品又粗又长| 亚洲精品456在线播放app| 免费观看人在逋| 精品久久久久久久末码| av女优亚洲男人天堂| 2021天堂中文幕一二区在线观| 国产 一区精品| 成人永久免费在线观看视频| 欧美日韩一区二区视频在线观看视频在线 | 99久久九九国产精品国产免费| 国产精品国产高清国产av| 国产 一区 欧美 日韩| 精品久久久久久久人妻蜜臀av| 在线播放无遮挡| 亚洲精品乱码久久久v下载方式| 国产色婷婷99| av女优亚洲男人天堂| 国产精品美女特级片免费视频播放器| 中文字幕精品亚洲无线码一区| 韩国av在线不卡| 99久国产av精品| 成年av动漫网址| 亚洲精品乱码久久久v下载方式| 欧美xxxx黑人xx丫x性爽| 性插视频无遮挡在线免费观看| 色吧在线观看| 精品午夜福利在线看| 亚洲国产欧美在线一区| avwww免费| 中文精品一卡2卡3卡4更新| 久久精品久久久久久噜噜老黄 | 1000部很黄的大片| 18禁黄网站禁片免费观看直播| 日韩精品青青久久久久久| 美女大奶头视频| 亚洲国产精品sss在线观看| 日韩成人伦理影院| 精品少妇黑人巨大在线播放 | 又爽又黄a免费视频| 蜜桃亚洲精品一区二区三区| 又黄又爽又刺激的免费视频.| 国产精品乱码一区二三区的特点| 欧美xxxx性猛交bbbb| 人妻久久中文字幕网| 韩国av在线不卡| 日本一本二区三区精品| 精品久久久久久成人av| 啦啦啦观看免费观看视频高清| 国产高清三级在线| 观看免费一级毛片| 少妇人妻精品综合一区二区 | 精品久久久久久久人妻蜜臀av| 欧美最新免费一区二区三区| 又粗又爽又猛毛片免费看| 亚洲欧美成人综合另类久久久 | 狠狠狠狠99中文字幕| 18禁裸乳无遮挡免费网站照片| 最近2019中文字幕mv第一页| 国内久久婷婷六月综合欲色啪| 国产成人一区二区在线| 国产高清不卡午夜福利| 最好的美女福利视频网| 国产成年人精品一区二区| 亚洲欧美日韩卡通动漫| 在线免费十八禁| 欧美成人免费av一区二区三区| 午夜视频国产福利| 国产精品国产高清国产av| 中文欧美无线码| 国产黄a三级三级三级人| av视频在线观看入口| 麻豆国产av国片精品| 美女高潮的动态| 成人综合一区亚洲| 亚洲欧美日韩东京热| 午夜福利高清视频| 波多野结衣高清作品| 久久久久久国产a免费观看| 性色avwww在线观看| 91精品国产九色| 国产成人a区在线观看| 国产毛片a区久久久久| 国产高潮美女av| avwww免费| 亚洲av第一区精品v没综合| 国产片特级美女逼逼视频| 内地一区二区视频在线| a级一级毛片免费在线观看| 免费看光身美女| 日韩人妻高清精品专区| 国产蜜桃级精品一区二区三区| 国产精品精品国产色婷婷| 亚洲国产精品成人久久小说 | 美女国产视频在线观看| 精品免费久久久久久久清纯| 乱系列少妇在线播放| 老女人水多毛片| 成人亚洲欧美一区二区av| 国产真实伦视频高清在线观看| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩高清专用| 日本黄色片子视频| 日本一本二区三区精品| 人人妻人人澡人人爽人人夜夜 | 欧美xxxx性猛交bbbb| 老女人水多毛片| 人妻制服诱惑在线中文字幕| 欧美精品一区二区大全| 我的女老师完整版在线观看| 免费搜索国产男女视频| 成人漫画全彩无遮挡| 男插女下体视频免费在线播放| 国产亚洲欧美98| 看黄色毛片网站| 亚洲一区二区三区色噜噜| 国产伦精品一区二区三区视频9| 日韩成人av中文字幕在线观看| 少妇熟女欧美另类| 国产成人a区在线观看| 国产高清不卡午夜福利| 国内揄拍国产精品人妻在线| 成人午夜高清在线视频| 久久精品国产亚洲网站| 国产精品电影一区二区三区| av.在线天堂| 少妇丰满av| 免费人成视频x8x8入口观看| 久久草成人影院| 国产伦理片在线播放av一区 | av在线天堂中文字幕| 一进一出抽搐动态| 青春草视频在线免费观看| 在线观看免费视频日本深夜| 亚洲成人中文字幕在线播放| 成人特级av手机在线观看| 亚洲欧美日韩高清专用| 丝袜喷水一区| 亚洲人成网站在线播放欧美日韩| 国产高清激情床上av| 搡老妇女老女人老熟妇| 51国产日韩欧美| 午夜精品一区二区三区免费看| 亚洲国产高清在线一区二区三| 免费观看a级毛片全部| 国产精品日韩av在线免费观看| 麻豆国产av国片精品| 可以在线观看毛片的网站| 男的添女的下面高潮视频| 亚州av有码| 禁无遮挡网站| 国产91av在线免费观看| 亚洲国产精品合色在线| 人人妻人人澡欧美一区二区| 床上黄色一级片| 99久国产av精品国产电影| 天堂网av新在线| 黄色一级大片看看| 亚洲精品成人久久久久久| 日本爱情动作片www.在线观看| 99热这里只有是精品在线观看| 黄色欧美视频在线观看| 永久网站在线| 男女做爰动态图高潮gif福利片| 狠狠狠狠99中文字幕| 免费在线观看成人毛片| 亚洲国产日韩欧美精品在线观看| 91午夜精品亚洲一区二区三区| 久久久久久大精品| 久久精品国产亚洲网站| 久久亚洲国产成人精品v| 中文字幕av成人在线电影| 床上黄色一级片| 日韩一区二区三区影片| 国产精品人妻久久久久久| 欧美又色又爽又黄视频| 麻豆国产97在线/欧美| 久99久视频精品免费| 18+在线观看网站| 国产三级中文精品| 亚洲欧洲日产国产| 99久久精品国产国产毛片| 国产亚洲av片在线观看秒播厂 | av视频在线观看入口| 亚洲久久久久久中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 久久99精品国语久久久| 美女xxoo啪啪120秒动态图| 非洲黑人性xxxx精品又粗又长| 日韩欧美在线乱码| 久久99精品国语久久久| 青春草亚洲视频在线观看| 夜夜夜夜夜久久久久| 亚洲成av人片在线播放无| 日产精品乱码卡一卡2卡三| 九九热线精品视视频播放| 国产成人aa在线观看| 六月丁香七月| av黄色大香蕉| 女同久久另类99精品国产91| 国内揄拍国产精品人妻在线| 91久久精品国产一区二区三区| 国产精品久久久久久精品电影小说 | 久久精品影院6| 热99re8久久精品国产| 国产精品人妻久久久久久| 自拍偷自拍亚洲精品老妇| 久久人人爽人人爽人人片va| 国产亚洲av嫩草精品影院| 午夜精品一区二区三区免费看| 亚洲一区二区三区色噜噜| 丰满乱子伦码专区| 国产日本99.免费观看| 搞女人的毛片| 国产一区二区三区av在线 | 九草在线视频观看| 亚洲欧洲国产日韩| 性色avwww在线观看| 国产探花在线观看一区二区| 1000部很黄的大片| 天堂√8在线中文| 1024手机看黄色片| 中文字幕久久专区| 日韩av不卡免费在线播放| 少妇的逼好多水| 亚洲精品国产av成人精品| 国产成人a区在线观看| 国产精品一区二区性色av| 国产精品爽爽va在线观看网站| 久久久国产成人精品二区| 少妇熟女aⅴ在线视频| 亚洲自拍偷在线| 国产视频首页在线观看| 欧美成人免费av一区二区三区| 国产精品.久久久| 热99在线观看视频| 亚洲欧美日韩卡通动漫| 国产乱人偷精品视频| 亚洲第一区二区三区不卡| 深夜精品福利| 最近的中文字幕免费完整| 日日摸夜夜添夜夜添av毛片| 久久99蜜桃精品久久| 男女视频在线观看网站免费| 最近最新中文字幕大全电影3| 国产极品精品免费视频能看的| 在线免费观看的www视频| 99riav亚洲国产免费| 亚洲精品乱码久久久v下载方式| 国产精品一区二区在线观看99 | 成熟少妇高潮喷水视频| av福利片在线观看| 国产一区二区在线av高清观看| 少妇人妻精品综合一区二区 | 国内少妇人妻偷人精品xxx网站| 欧美性猛交╳xxx乱大交人| 中文字幕制服av| 不卡一级毛片| 我要搜黄色片| 一本精品99久久精品77| 色哟哟哟哟哟哟| 日韩强制内射视频| 欧美极品一区二区三区四区| 欧美三级亚洲精品| 亚洲最大成人av| 欧美高清性xxxxhd video| 精品一区二区三区人妻视频| or卡值多少钱| 中文字幕制服av| 精品熟女少妇av免费看| 一边亲一边摸免费视频| 欧美丝袜亚洲另类| 久久久久国产网址| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产人妻一区二区三区在| 亚洲国产日韩欧美精品在线观看| 看非洲黑人一级黄片| 最近手机中文字幕大全| 波多野结衣巨乳人妻| 26uuu在线亚洲综合色| 亚洲国产精品成人久久小说 | 欧美zozozo另类| 精品久久久噜噜| 日产精品乱码卡一卡2卡三| 精品久久久久久久末码| 国产欧美日韩精品一区二区| 91aial.com中文字幕在线观看| 插阴视频在线观看视频| 美女被艹到高潮喷水动态| 97人妻精品一区二区三区麻豆| videossex国产| 国产激情偷乱视频一区二区| 97人妻精品一区二区三区麻豆| 久久午夜福利片| 少妇的逼水好多| 久久亚洲精品不卡| 可以在线观看的亚洲视频| 少妇熟女aⅴ在线视频| 欧美精品国产亚洲| 婷婷色av中文字幕| 小说图片视频综合网站| 欧美一区二区国产精品久久精品| 成人鲁丝片一二三区免费| 亚洲中文字幕一区二区三区有码在线看| 黄色视频,在线免费观看| 老师上课跳d突然被开到最大视频| 久久久久久九九精品二区国产| 欧美成人a在线观看| 看片在线看免费视频| 成人国产麻豆网| 欧美xxxx性猛交bbbb| 亚洲成人av在线免费| 九九热线精品视视频播放| 久久精品国产亚洲av天美| 91久久精品电影网| 卡戴珊不雅视频在线播放| 国产精品电影一区二区三区| 欧美日韩精品成人综合77777| 亚洲欧美日韩高清专用| 国产毛片a区久久久久| 99热只有精品国产| 性欧美人与动物交配| 插逼视频在线观看| 舔av片在线| 天堂网av新在线| 久久精品国产自在天天线| 久99久视频精品免费| 免费无遮挡裸体视频| 综合色av麻豆| 亚洲不卡免费看| 日韩av不卡免费在线播放| 一区二区三区免费毛片| 亚洲国产精品久久男人天堂| 一区二区三区高清视频在线| 91久久精品国产一区二区成人| 99久久精品一区二区三区| 男女做爰动态图高潮gif福利片| 在线免费观看不下载黄p国产| 国产老妇女一区| 九九热线精品视视频播放| 久久久久久久久久久丰满| 蜜臀久久99精品久久宅男| 日韩高清综合在线| 日韩国内少妇激情av| 国产精品免费一区二区三区在线| 国产日本99.免费观看| 欧美精品一区二区大全| 男女视频在线观看网站免费| 欧美潮喷喷水| 草草在线视频免费看| 国产片特级美女逼逼视频| 国产色婷婷99| 国产午夜精品久久久久久一区二区三区| 欧美激情在线99| 中文字幕制服av| ponron亚洲| 成人无遮挡网站| 中文字幕av在线有码专区| 国产精品免费一区二区三区在线| 熟女人妻精品中文字幕| 国产av麻豆久久久久久久| 国产蜜桃级精品一区二区三区| 夜夜夜夜夜久久久久| 色哟哟哟哟哟哟| 亚洲五月天丁香| 午夜精品一区二区三区免费看| 一区福利在线观看| 亚洲成a人片在线一区二区| 国产一级毛片七仙女欲春2| 免费搜索国产男女视频| 免费av观看视频| 波多野结衣高清作品| 成年av动漫网址| 国产精品综合久久久久久久免费| 欧美激情在线99| 日韩成人伦理影院| 国产精品免费一区二区三区在线| 极品教师在线视频| av天堂在线播放| 麻豆成人午夜福利视频| 成年版毛片免费区| 亚洲最大成人av| 欧美一级a爱片免费观看看| 少妇丰满av| 国产69精品久久久久777片| 深夜精品福利| 国产淫片久久久久久久久| 夫妻性生交免费视频一级片| 一级毛片久久久久久久久女| 国产av麻豆久久久久久久| 色吧在线观看| 国产精品久久久久久久久免| 1000部很黄的大片| 美女cb高潮喷水在线观看| 婷婷亚洲欧美| 在线观看美女被高潮喷水网站| 久久这里有精品视频免费| 在线免费观看的www视频| 亚洲美女搞黄在线观看| 日日啪夜夜撸| 亚洲欧美日韩高清专用| 女同久久另类99精品国产91| 狠狠狠狠99中文字幕| 久久久久免费精品人妻一区二区| 精品午夜福利在线看| av黄色大香蕉| 久久精品国产自在天天线| 成人鲁丝片一二三区免费| 熟女电影av网| 嫩草影院入口| 久久精品国产亚洲av涩爱 | 精品久久久久久久人妻蜜臀av| 亚洲自偷自拍三级| 亚洲不卡免费看| 精品国产三级普通话版| 亚洲欧美精品综合久久99| 2021天堂中文幕一二区在线观| 国产真实乱freesex| 插逼视频在线观看| 欧美成人一区二区免费高清观看| 欧美人与善性xxx| 亚洲精品粉嫩美女一区| 国产在线精品亚洲第一网站| 亚洲欧美日韩无卡精品| 最近的中文字幕免费完整| 身体一侧抽搐| 久久久久九九精品影院| 久久鲁丝午夜福利片| 日本免费a在线| 18+在线观看网站| 成年女人永久免费观看视频| 日韩视频在线欧美| 亚洲综合色惰| 国产精品久久久久久久久免| 国产精品福利在线免费观看| 国产黄色视频一区二区在线观看 | 伦理电影大哥的女人| 超碰av人人做人人爽久久| 亚洲三级黄色毛片| 丰满人妻一区二区三区视频av| 尤物成人国产欧美一区二区三区| 亚洲欧洲国产日韩| 欧美另类亚洲清纯唯美| 久久久久久九九精品二区国产| 久久九九热精品免费| 夜夜看夜夜爽夜夜摸| 午夜免费激情av| 亚洲中文字幕日韩| 亚洲精品影视一区二区三区av| 白带黄色成豆腐渣| 卡戴珊不雅视频在线播放| 国产精品国产高清国产av| 91在线精品国自产拍蜜月| 国产伦在线观看视频一区| 深爱激情五月婷婷| 日韩在线高清观看一区二区三区| 综合色av麻豆| 97超碰精品成人国产| 亚洲精品粉嫩美女一区| 可以在线观看毛片的网站| 国产一区二区在线av高清观看| 免费一级毛片在线播放高清视频| 五月玫瑰六月丁香| 国内精品久久久久精免费| 爱豆传媒免费全集在线观看| 国产又黄又爽又无遮挡在线| 久久中文看片网| 国产色婷婷99| 少妇的逼好多水| 少妇的逼水好多| 国产一区二区亚洲精品在线观看| 国产探花在线观看一区二区| 国产av不卡久久| 天堂av国产一区二区熟女人妻| 一级毛片电影观看 | 最近最新中文字幕大全电影3| 国产精品,欧美在线| 国产 一区 欧美 日韩| 99热只有精品国产| 亚洲欧美日韩高清在线视频| 一边亲一边摸免费视频| 日韩一区二区视频免费看| 亚洲内射少妇av| 啦啦啦韩国在线观看视频| 麻豆国产97在线/欧美| 成人午夜精彩视频在线观看| 97在线视频观看| 国产伦理片在线播放av一区 | 亚洲国产精品国产精品| 亚洲美女视频黄频| 国产精品野战在线观看| 九色成人免费人妻av| 欧美变态另类bdsm刘玥| 中国国产av一级| 国产免费男女视频| av在线播放精品| 白带黄色成豆腐渣| 国产精品久久久久久av不卡| 在线a可以看的网站| 午夜福利视频1000在线观看| 欧美精品一区二区大全| 欧美色欧美亚洲另类二区| 哪个播放器可以免费观看大片| 高清毛片免费看| 色吧在线观看| 久久这里有精品视频免费| 18禁在线播放成人免费| 少妇高潮的动态图| 亚洲va在线va天堂va国产| 日韩欧美 国产精品| 麻豆精品久久久久久蜜桃| 少妇被粗大猛烈的视频| 亚洲第一电影网av| 国产成人精品婷婷| 国产成人午夜福利电影在线观看| av在线天堂中文字幕| 日本熟妇午夜| 好男人视频免费观看在线| 成年女人看的毛片在线观看| 日韩 亚洲 欧美在线| 日本与韩国留学比较| 极品教师在线视频| 久久精品国产清高在天天线| 一本精品99久久精品77| 亚洲18禁久久av| 18+在线观看网站| 国产黄色小视频在线观看| h日本视频在线播放| 给我免费播放毛片高清在线观看| 成人毛片60女人毛片免费| 一级黄片播放器| 亚洲欧美精品专区久久| 高清毛片免费看| 长腿黑丝高跟| 18+在线观看网站| 少妇猛男粗大的猛烈进出视频 | 99热只有精品国产| 国产伦精品一区二区三区四那| 我要搜黄色片| 国产精品一区二区在线观看99 | 日本黄色视频三级网站网址| 国产精华一区二区三区| 亚洲激情五月婷婷啪啪| 免费看日本二区| 国产一区二区在线观看日韩| 国产又黄又爽又无遮挡在线| 亚洲精品粉嫩美女一区| av专区在线播放| 国产黄色小视频在线观看| 伦理电影大哥的女人| 国产淫片久久久久久久久| 久久99蜜桃精品久久| 免费黄网站久久成人精品| 国产精品国产三级国产av玫瑰| 免费一级毛片在线播放高清视频| 免费大片18禁| 免费av不卡在线播放| 亚洲国产色片| 成人三级黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 美女黄网站色视频| 国产精品人妻久久久久久|