• <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在线视频| 日韩大尺度精品在线看网址 | 国产欧美日韩一区二区三| 俄罗斯特黄特色一大片| 99精品久久久久人妻精品| 国产精品 欧美亚洲| 无限看片的www在线观看| 国产高清视频在线播放一区| 欧美黑人欧美精品刺激| 天天添夜夜摸| 99香蕉大伊视频| 亚洲成av人片免费观看| 久久香蕉激情| 久久久久久亚洲精品国产蜜桃av| 国产激情久久老熟女| 亚洲欧美精品综合久久99| 99久久精品国产亚洲精品| 精品无人区乱码1区二区| 18禁黄网站禁片午夜丰满| 中文字幕人妻熟女乱码| 亚洲一区高清亚洲精品| 国产高清videossex| 国产成人啪精品午夜网站| 精品欧美一区二区三区在线| 校园春色视频在线观看| 欧美国产日韩亚洲一区| 久9热在线精品视频| 国内毛片毛片毛片毛片毛片| 非洲黑人性xxxx精品又粗又长| 日韩大尺度精品在线看网址 | 天天添夜夜摸| 电影成人av| 亚洲成国产人片在线观看| 高清在线国产一区| 少妇熟女aⅴ在线视频| 亚洲中文字幕日韩| 欧美午夜高清在线| 99riav亚洲国产免费| 国产欧美日韩综合在线一区二区| 中文字幕高清在线视频| 亚洲精品中文字幕在线视频| 少妇被粗大的猛进出69影院| 一区二区三区精品91| 成人国语在线视频| 亚洲一区高清亚洲精品| 日日爽夜夜爽网站| 女同久久另类99精品国产91| 国产精品久久久久久精品电影 | 99在线人妻在线中文字幕| 欧美成人性av电影在线观看| 法律面前人人平等表现在哪些方面| 国产精品野战在线观看| 欧美成人午夜精品| 男人舔女人下体高潮全视频| 视频在线观看一区二区三区| 人妻久久中文字幕网| 不卡av一区二区三区| 黄色成人免费大全| 精品国产亚洲在线| bbb黄色大片| 丝袜美腿诱惑在线| 欧美激情高清一区二区三区| а√天堂www在线а√下载| 亚洲人成77777在线视频| 99国产精品免费福利视频| 国产精品,欧美在线| 国产一区二区三区综合在线观看| 男女做爰动态图高潮gif福利片 | 99香蕉大伊视频| 国产精品永久免费网站| 亚洲一区中文字幕在线| 国产精品久久久久久精品电影 | 欧美黄色片欧美黄色片| 精品乱码久久久久久99久播| 这个男人来自地球电影免费观看| 午夜成年电影在线免费观看| 国产亚洲精品一区二区www| 午夜亚洲福利在线播放| 1024视频免费在线观看| 50天的宝宝边吃奶边哭怎么回事| 变态另类成人亚洲欧美熟女 | 国产成人精品无人区| 禁无遮挡网站| 国产精品 欧美亚洲| 在线观看免费视频日本深夜| 日韩欧美一区视频在线观看| а√天堂www在线а√下载| 国产精品秋霞免费鲁丝片| 午夜久久久在线观看| 丝袜美足系列| 婷婷精品国产亚洲av在线| 日韩欧美三级三区| 99riav亚洲国产免费| 校园春色视频在线观看| 国产亚洲精品一区二区www| 亚洲一区二区三区色噜噜| 午夜a级毛片| or卡值多少钱| 可以在线观看的亚洲视频| 欧美大码av| 欧美日韩瑟瑟在线播放| 神马国产精品三级电影在线观看 | 欧美日韩瑟瑟在线播放| 亚洲精品在线观看二区| or卡值多少钱| 最新在线观看一区二区三区| 国产精品免费视频内射| 免费少妇av软件| 国产伦人伦偷精品视频| 欧美中文日本在线观看视频| 国产精品99久久99久久久不卡| 久久国产乱子伦精品免费另类| 黄色视频不卡| 在线观看日韩欧美| 伦理电影免费视频| 又黄又爽又免费观看的视频| 亚洲成人精品中文字幕电影| 久久中文字幕一级| 亚洲人成伊人成综合网2020| 久久 成人 亚洲| 免费无遮挡裸体视频| 日韩欧美三级三区| 午夜免费激情av| 日韩大尺度精品在线看网址 | tocl精华| 成人18禁高潮啪啪吃奶动态图| 国产99白浆流出| 国产私拍福利视频在线观看| 黄色片一级片一级黄色片| 久久午夜综合久久蜜桃| 婷婷精品国产亚洲av在线| 亚洲精品国产精品久久久不卡| 女警被强在线播放| 9191精品国产免费久久| 激情视频va一区二区三区| АⅤ资源中文在线天堂| 欧美乱妇无乱码| 丝袜美足系列| 久久国产乱子伦精品免费另类| 国产一区二区三区综合在线观看| 国产亚洲欧美98| 精品久久久久久久人妻蜜臀av | 亚洲少妇的诱惑av| 老汉色av国产亚洲站长工具| 午夜福利视频1000在线观看 | 亚洲人成网站在线播放欧美日韩| 亚洲欧美激情在线| 在线免费观看的www视频| 国产一区二区三区在线臀色熟女| 成人国语在线视频| 亚洲一区二区三区色噜噜| 免费在线观看完整版高清| 国产精品美女特级片免费视频播放器 | 中出人妻视频一区二区| 亚洲专区字幕在线| 美女午夜性视频免费| 久久精品国产亚洲av高清一级| 热99re8久久精品国产| 亚洲在线自拍视频| 成人三级黄色视频| 欧美在线一区亚洲| 免费在线观看黄色视频的| av电影中文网址| 亚洲专区中文字幕在线| 亚洲av成人av| 99久久国产精品久久久| 亚洲国产日韩欧美精品在线观看 | 国产伦一二天堂av在线观看| 最好的美女福利视频网| 18禁美女被吸乳视频| 午夜免费激情av| 午夜免费激情av| 亚洲人成电影观看| 亚洲五月婷婷丁香| 真人一进一出gif抽搐免费| 成人国语在线视频| 欧美另类亚洲清纯唯美| 老司机福利观看| 91大片在线观看| www.熟女人妻精品国产| bbb黄色大片| 欧美日韩福利视频一区二区| 免费在线观看亚洲国产| 久久伊人香网站| 老司机在亚洲福利影院| 午夜免费观看网址| 日本a在线网址| 亚洲精品国产一区二区精华液| 亚洲自偷自拍图片 自拍| 美女免费视频网站| 嫁个100分男人电影在线观看| 精品国内亚洲2022精品成人| 免费看a级黄色片| 日本 欧美在线| 久久性视频一级片| 欧美 亚洲 国产 日韩一| 久久天堂一区二区三区四区| 在线av久久热| av免费在线观看网站| 99国产精品免费福利视频| 大香蕉久久成人网| 久久精品国产综合久久久| 巨乳人妻的诱惑在线观看| 成人手机av| 91大片在线观看| 国产亚洲精品第一综合不卡| 国产精品久久久久久人妻精品电影| 国产日韩一区二区三区精品不卡| 嫩草影视91久久| 女生性感内裤真人,穿戴方法视频| 精品欧美一区二区三区在线| 又紧又爽又黄一区二区| 国产成人精品在线电影| 久久精品国产99精品国产亚洲性色 | 国产日韩一区二区三区精品不卡| 亚洲美女黄片视频| 精品国产乱码久久久久久男人| 久久草成人影院| 亚洲国产日韩欧美精品在线观看 | 99久久99久久久精品蜜桃| 亚洲美女黄片视频| 中国美女看黄片| 午夜久久久久精精品| 亚洲电影在线观看av| 欧美黄色片欧美黄色片| 国产欧美日韩综合在线一区二区| 欧美激情久久久久久爽电影 | 久久国产精品人妻蜜桃| 一区二区三区激情视频| 日本 av在线| bbb黄色大片| 亚洲五月色婷婷综合| 免费在线观看完整版高清| 女人精品久久久久毛片| 日本黄色视频三级网站网址| 欧美日韩乱码在线| 精品第一国产精品| 国产xxxxx性猛交| 电影成人av| 亚洲av成人一区二区三| 在线永久观看黄色视频| 嫁个100分男人电影在线观看| 久久精品国产亚洲av高清一级| 少妇粗大呻吟视频| 精品电影一区二区在线| 国产成人欧美| 亚洲精品在线观看二区| 色哟哟哟哟哟哟| 女同久久另类99精品国产91| 国产精品永久免费网站| 欧美另类亚洲清纯唯美| 99久久综合精品五月天人人| 日韩欧美在线二视频| 成人精品一区二区免费| 亚洲精品av麻豆狂野| 手机成人av网站| 亚洲 欧美 日韩 在线 免费| 亚洲一区二区三区色噜噜| 亚洲av成人不卡在线观看播放网| 在线观看免费视频网站a站| 国产精品爽爽va在线观看网站 | 亚洲熟女毛片儿| 一级a爱视频在线免费观看| 日本黄色视频三级网站网址| 欧美+亚洲+日韩+国产| av天堂在线播放| 亚洲欧美日韩无卡精品| 成人三级黄色视频| 午夜免费观看网址| 国产亚洲精品久久久久5区| 一本久久中文字幕| 午夜老司机福利片| 亚洲第一电影网av| 免费一级毛片在线播放高清视频 | aaaaa片日本免费| 一卡2卡三卡四卡精品乱码亚洲| 国产精品野战在线观看| 啦啦啦观看免费观看视频高清 | 亚洲国产高清在线一区二区三 | 一进一出好大好爽视频| 最近最新中文字幕大全免费视频| 女人精品久久久久毛片| 亚洲免费av在线视频| √禁漫天堂资源中文www| 热re99久久国产66热| 黑人巨大精品欧美一区二区mp4| 亚洲国产精品合色在线| 精品国产国语对白av| 9色porny在线观看| 成人国语在线视频| 欧美一级毛片孕妇| 日本三级黄在线观看| 国产在线精品亚洲第一网站| 波多野结衣巨乳人妻| 亚洲熟女毛片儿| 日韩精品免费视频一区二区三区| 色综合站精品国产| 日韩欧美免费精品| 好男人在线观看高清免费视频 | 国产一卡二卡三卡精品| 韩国精品一区二区三区| 国产精品九九99| av有码第一页| 成人永久免费在线观看视频| 欧美乱码精品一区二区三区| 国产亚洲av高清不卡| 日本免费一区二区三区高清不卡 | 午夜免费观看网址| 在线免费观看的www视频| 国产又色又爽无遮挡免费看| 美女午夜性视频免费| 男人舔女人下体高潮全视频| 无遮挡黄片免费观看| 国产精品九九99| 最近最新免费中文字幕在线| 黄色丝袜av网址大全| 国内精品久久久久精免费| 欧美乱色亚洲激情| 亚洲少妇的诱惑av| 一区二区日韩欧美中文字幕| 色综合站精品国产| 精品一区二区三区四区五区乱码| 一级作爱视频免费观看| 国产亚洲欧美精品永久| 成人精品一区二区免费| 夜夜看夜夜爽夜夜摸| 久久狼人影院| 久久久久久久午夜电影| 给我免费播放毛片高清在线观看| 欧美中文日本在线观看视频| 好看av亚洲va欧美ⅴa在| 亚洲全国av大片| 午夜福利一区二区在线看| 久久中文字幕一级| 欧美黑人精品巨大| 成人av一区二区三区在线看| 亚洲国产高清在线一区二区三 | 午夜激情av网站| 欧美一区二区精品小视频在线| 午夜福利欧美成人| 久久久精品欧美日韩精品| 51午夜福利影视在线观看| 成年人黄色毛片网站| 欧美中文综合在线视频| 免费在线观看日本一区| 亚洲 欧美一区二区三区| www日本在线高清视频| 午夜影院日韩av| 亚洲久久久国产精品| 电影成人av| 久久精品91无色码中文字幕| 精品久久久久久久久久免费视频| 亚洲成人免费电影在线观看| 国产欧美日韩一区二区三| 一级片免费观看大全| 多毛熟女@视频| 久久久久久久久免费视频了| av电影中文网址| 每晚都被弄得嗷嗷叫到高潮| 亚洲男人天堂网一区| 亚洲av第一区精品v没综合| 免费在线观看亚洲国产| 日日夜夜操网爽| 色综合婷婷激情| 日韩中文字幕欧美一区二区| 久久中文看片网| 高清黄色对白视频在线免费看| 人人妻人人澡人人看| 啦啦啦韩国在线观看视频| 亚洲成人免费电影在线观看| 88av欧美| 午夜精品久久久久久毛片777| 人人妻人人澡人人看| 国产成人免费无遮挡视频| 亚洲欧美精品综合一区二区三区| 神马国产精品三级电影在线观看 | 亚洲人成网站在线播放欧美日韩| 美女大奶头视频| 自线自在国产av| av欧美777| 午夜精品久久久久久毛片777| 久久伊人香网站| 老司机深夜福利视频在线观看| 午夜福利18| 黑丝袜美女国产一区| 好男人在线观看高清免费视频 | 亚洲男人天堂网一区| 午夜福利在线观看吧| 热99re8久久精品国产| 中出人妻视频一区二区| 欧美最黄视频在线播放免费| 777久久人妻少妇嫩草av网站| 精品久久久久久,| 大型av网站在线播放| 欧美色欧美亚洲另类二区 | 免费在线观看亚洲国产| 视频在线观看一区二区三区| 日日干狠狠操夜夜爽| 免费在线观看黄色视频的| 国产精品久久电影中文字幕| 91九色精品人成在线观看| 少妇 在线观看| 波多野结衣av一区二区av| 亚洲精品国产精品久久久不卡| 亚洲av熟女| 亚洲一码二码三码区别大吗| 久久久久九九精品影院| 免费在线观看黄色视频的| 高潮久久久久久久久久久不卡| 精品电影一区二区在线| 欧美黑人欧美精品刺激| 99香蕉大伊视频| 欧美一级a爱片免费观看看 | 18禁国产床啪视频网站| 久久精品国产亚洲av高清一级| 欧美色视频一区免费| 美女扒开内裤让男人捅视频| 亚洲 欧美一区二区三区| 午夜福利高清视频| 亚洲熟妇中文字幕五十中出| bbb黄色大片| 亚洲三区欧美一区| 丁香欧美五月| 人人妻人人澡欧美一区二区 | 日韩av在线大香蕉| 国产亚洲欧美98| 国产精品乱码一区二三区的特点 | 宅男免费午夜| 久久久久久久精品吃奶| 亚洲欧美精品综合久久99| 99国产精品99久久久久| 欧美绝顶高潮抽搐喷水| 97人妻精品一区二区三区麻豆 | 亚洲av电影不卡..在线观看| 搡老熟女国产l中国老女人| 777久久人妻少妇嫩草av网站| 好男人电影高清在线观看| 亚洲最大成人中文| 成人亚洲精品av一区二区| 无限看片的www在线观看| 他把我摸到了高潮在线观看| 热99re8久久精品国产| 我的亚洲天堂| 国产一卡二卡三卡精品| av网站免费在线观看视频| 精品高清国产在线一区| 美女午夜性视频免费| 后天国语完整版免费观看| 精品第一国产精品| 国产高清视频在线播放一区| 亚洲熟妇熟女久久| 正在播放国产对白刺激| 国产伦一二天堂av在线观看| 99久久久亚洲精品蜜臀av| 久久中文字幕人妻熟女| 日韩欧美在线二视频| 国产午夜福利久久久久久| 国产黄a三级三级三级人| 很黄的视频免费| 成人永久免费在线观看视频| 一级a爱片免费观看的视频| 99国产精品一区二区蜜桃av| 正在播放国产对白刺激| 999久久久国产精品视频| 国产黄a三级三级三级人| 男女下面插进去视频免费观看| 成年女人毛片免费观看观看9| 久久精品国产综合久久久| 男人舔女人的私密视频| 制服丝袜大香蕉在线| 久久性视频一级片| 国内久久婷婷六月综合欲色啪| 亚洲三区欧美一区| av有码第一页| 欧美成狂野欧美在线观看| 成人精品一区二区免费| 波多野结衣高清无吗| 少妇裸体淫交视频免费看高清 | 激情视频va一区二区三区| 国产欧美日韩精品亚洲av| 国产麻豆69| 少妇熟女aⅴ在线视频| 国产成年人精品一区二区| 黄色片一级片一级黄色片| 美国免费a级毛片| 国产一区在线观看成人免费| 性色av乱码一区二区三区2| 免费在线观看日本一区| 亚洲精品中文字幕在线视频| 免费女性裸体啪啪无遮挡网站| 黑人操中国人逼视频| 女生性感内裤真人,穿戴方法视频| 麻豆国产av国片精品| 亚洲成国产人片在线观看| 成人特级黄色片久久久久久久| 色精品久久人妻99蜜桃| 最新美女视频免费是黄的| 成人三级黄色视频| 免费人成视频x8x8入口观看| 亚洲欧美日韩高清在线视频| 韩国精品一区二区三区| 亚洲精华国产精华精| 两个人看的免费小视频| 一区福利在线观看| 精品熟女少妇八av免费久了| 欧美日本视频| 色在线成人网| 久久性视频一级片| 国产成人免费无遮挡视频| 夜夜看夜夜爽夜夜摸| 久久精品国产清高在天天线| 久久精品人人爽人人爽视色| 老鸭窝网址在线观看| 久久青草综合色| 久久欧美精品欧美久久欧美| 亚洲aⅴ乱码一区二区在线播放 | 亚洲欧美日韩高清在线视频| 在线观看一区二区三区| 亚洲五月婷婷丁香| 精品国产国语对白av| 宅男免费午夜| 一级片免费观看大全| 亚洲中文字幕日韩| 欧美激情 高清一区二区三区| 校园春色视频在线观看| 在线观看66精品国产| 国产精品精品国产色婷婷| 午夜久久久久精精品| 亚洲美女黄片视频| 波多野结衣一区麻豆| 欧美日韩瑟瑟在线播放| 亚洲精品粉嫩美女一区| 狂野欧美激情性xxxx| 亚洲av五月六月丁香网| 宅男免费午夜| 美女 人体艺术 gogo| 国产精品98久久久久久宅男小说| 久久久久久大精品| 免费在线观看日本一区| 99国产精品免费福利视频| 丰满人妻熟妇乱又伦精品不卡| 丰满的人妻完整版| 级片在线观看| 色综合亚洲欧美另类图片| 中文字幕最新亚洲高清| 国产精品二区激情视频| 亚洲情色 制服丝袜| 人人妻人人爽人人添夜夜欢视频| 久久影院123| 在线观看免费视频网站a站| 国产亚洲欧美精品永久| 亚洲自拍偷在线| 国产97色在线日韩免费| 久久久久亚洲av毛片大全| 在线av久久热| 国产野战对白在线观看| 亚洲成av人片免费观看| 美女国产高潮福利片在线看| 黄色视频不卡| 国产伦人伦偷精品视频| 久久国产亚洲av麻豆专区| 丁香欧美五月| 精品无人区乱码1区二区| 色老头精品视频在线观看| 欧美黑人欧美精品刺激| 首页视频小说图片口味搜索| 纯流量卡能插随身wifi吗| 久久人妻熟女aⅴ| 亚洲欧美日韩高清在线视频| 久久久久久人人人人人| 极品教师在线免费播放| 99精品在免费线老司机午夜| 午夜免费成人在线视频| 久久亚洲真实| 真人做人爱边吃奶动态| 91麻豆精品激情在线观看国产| 亚洲色图 男人天堂 中文字幕| 国产免费男女视频| 91精品三级在线观看| 国产人伦9x9x在线观看| av在线播放免费不卡| 日本撒尿小便嘘嘘汇集6| 少妇 在线观看| 国产午夜精品久久久久久| 免费看美女性在线毛片视频| 搞女人的毛片| 美女高潮喷水抽搐中文字幕| 1024视频免费在线观看| 国产1区2区3区精品| 99国产极品粉嫩在线观看| 国产不卡一卡二| 欧美日韩福利视频一区二区| 日本黄色视频三级网站网址| 久久精品91无色码中文字幕| 少妇熟女aⅴ在线视频| 国产亚洲欧美在线一区二区| 天堂影院成人在线观看| 亚洲成人免费电影在线观看| 国产色视频综合| 麻豆久久精品国产亚洲av| 亚洲欧美激情综合另类| 涩涩av久久男人的天堂| 女人被狂操c到高潮| 国产成人欧美| 热99re8久久精品国产| 亚洲中文字幕日韩|