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

    利用SAS軟件快速方差分析多組資料

    2018-06-13 06:58:12張白鴿
    生物學(xué)雜志 2018年3期
    關(guān)鍵詞:語句葉綠素辣椒

    宋 釗, 張白鴿, 曹 健

    (廣東省農(nóng)業(yè)科學(xué)院蔬菜研究所 廣東省蔬菜新技術(shù)研究重點實驗室,廣州 510640)

    SAS的英文全名是Statistical Analysis System,即統(tǒng)計分析系統(tǒng),從20世紀(jì)60年代末期由美國北卡羅來納州立大學(xué)兩位教授開始發(fā)展的,是當(dāng)今國際上最著名的數(shù)據(jù)分析軟件之一,目前已經(jīng)與SPSS(Statistical product and service solutions)和BMDP(Bio Medical Data Processing)并駕齊驅(qū)[1],在國際上被譽為數(shù)據(jù)處理和統(tǒng)計分析的標(biāo)準(zhǔn)軟件系統(tǒng),在國內(nèi)廣泛用于農(nóng)業(yè)試驗統(tǒng)計分析[2]。

    對于SAS,SPSS和EXCEL等幾個數(shù)據(jù)統(tǒng)計軟件的使用已經(jīng)有多篇報道從各方面進行了較多研究[3-5],比較各個軟件在實際使用過程的優(yōu)缺點。由于SAS系統(tǒng)為全英文界面,對于英語基礎(chǔ)比較薄弱的使用者存在一定困難,并且程序模塊化,對于沒有計算機程序語言編程基礎(chǔ)的人員來說不易掌握[5],所以很長一段時間內(nèi)國內(nèi)研究者在進行單因素或者多因素方差分析時大多使用在農(nóng)業(yè)生物統(tǒng)計中得到廣泛應(yīng)用的Excel[6-7]。但是Excel只能給出方差分析表,不能進行平均數(shù)的多重比較和字母標(biāo)記差異顯著性結(jié)果[5],不能同時對多組觀察值一次性全部進行方差分析,如果樣本數(shù)較多將費時費力并且容易出錯。本文介紹利用SAS統(tǒng)計軟件同時對20組辣椒品種在3個不同脅迫條件下辣椒葉片葉綠素含量數(shù)據(jù)進行方差分析,檢測同一個品種在不同澇漬脅迫條件葉片葉綠素含量差異顯著性,同時對20個品種間的葉綠素含量差異進行比較。

    1 材料與方法

    1.1 實驗材料

    總計20份辣椒材料,為辣椒屬一年生種(CapsicumannuumL.),按照果實形狀可以分為尖椒、甜椒、美人椒和線椒4個類型,均為在華南地區(qū)栽培較廣,具有一定代表性的栽培種。

    1.2 實驗方法

    試驗在廣州市天河區(qū)五山路廣東省農(nóng)業(yè)科學(xué)院蔬菜研究所旁邊科研基地大棚中進行,2015年4月15號播種,5月12間苗,5月15號開始脅迫處理。試驗設(shè)置1個對照和澇漬脅迫3、6和9 d等3個處理,然后恢復(fù)3 d,每個處理3個重復(fù),每個重復(fù)每個品種10株,穴盤1孔放1株,品種隨機排序。對照進行正常田間管理,缺水后即澆水。澇漬脅迫方法參照尹冬梅等[8]的方法加以改進,水層高于辣椒植株根部2~3 cm,用帶體積刻度的量杯倒入等體積水量開始脅迫處理。

    1.3 葉綠素測定

    使用日本柯尼卡美能達(Konica Minolta)SPAD-502 PLUS葉綠素?zé)晒鈨x對葉片的葉綠素含量進行測定,該儀器通過測量葉子對兩個波長段里的吸收率,來評估當(dāng)前葉子中的葉綠素的相對含量,用SPAD值表示目前葉子中葉綠素含量相對應(yīng)的參數(shù)。

    1.4 數(shù)據(jù)分析軟件和方法

    The SAS System for Windows Version 8.01 TS Level 01M0;SAS Institute,Cary,NC,USA(美國北卡羅來納州SAS研究所出品)。本試驗為平衡試驗設(shè)計,采用過程步PROC的ANOVA進行方差分析,同時進行Duncan′s multiple-range test多重比較。

    2 結(jié)果與分析

    2.1 辣椒葉片葉綠素測量結(jié)果

    表1為澇漬脅迫處理后對照和3個處理的植株葉片葉綠素SPAD值,每個品種每個處理測量15個數(shù)據(jù),20個品種4組數(shù)據(jù)總計1200個數(shù)據(jù)。限于篇幅表1每個品種只列出了5個重復(fù)數(shù)據(jù),省去了小數(shù)點后面數(shù)值。

    表1 20份辣椒材料澇漬脅迫下葉片葉綠素?zé)晒釹PAD值Table 1 Chlorophyll SPAD values of leaves of 20 pepper varieties under waterlogging stress

    2.2 比較每個品種處理間的差異顯著性

    2.2.1 建立SAS數(shù)據(jù)集

    如表1數(shù)據(jù)在SAS系統(tǒng)中是不能直接使用的,需要轉(zhuǎn)換成SAS的格式,也就是建立數(shù)據(jù)集,我們通過SAS語言的數(shù)據(jù)步DATA中的INPUT函數(shù)和CARDS選項建立數(shù)據(jù)集。

    核心代碼及說明如下:

    DATA Work.yelvsu; /*輸入20個辣椒品種葉綠素數(shù)據(jù)*/

    DO trt=1 TO 4; /*4個處理數(shù)序號,保留*/

    DO t=1 TO 15;

    INPUT brd1-brd20@@;

    OUTPUT;

    END;

    END;

    DROP t; /*15個重復(fù)數(shù)序號,丟棄*/

    CARDS;

    38 32 32 42 37 34 29 32 32 38 37 37 35 44 40 35 28 30 37 30

    40 35 33 38 40 36 37 29 39 38 36 36 38 41 42 39 31 33 37 34

    40 34 35 39 32 36 27 27 43 34 39 35 37 38 35 39 38 33 33 33

    39 39 29 35 37 39 31 29 36 33 42 36 34 36 38 39 39 37 34 31

    41 36 38 37 32 33 33 29 31 25 38 37 39 37 37 40 28 28 36 34

    ……(省略部分類似數(shù)據(jù))

    ;

    RUN;

    “DATA Work.yelvsu;”語句建立臨時庫Work的數(shù)據(jù)集yelvsu;“DO trt=1 TO 4;”語句表示在最外層進行4次循環(huán),因為有4組數(shù)據(jù);“DO t=1 TO 15;”語句表示連續(xù)寫入15個重復(fù)值;“DROP t;”語句則表示重復(fù)次數(shù)1到15不需要寫入數(shù)據(jù)集,只寫入相對應(yīng)的SPAD值;“INPUT brd1-brd20@@;”語句表示一行寫入20個品種的SPAD值;“OUTPUT;”語句表示將過程步中的數(shù)據(jù)輸出到數(shù)據(jù)集中保存起來;“CARDS;”語句則列出INPUT所有變量對應(yīng)的數(shù)據(jù),以“;”表示數(shù)據(jù)的結(jié)束。執(zhí)行之后自動建立trt和brd1,brd2,brd3…,brd20等21列數(shù)據(jù)。只有當(dāng)原始數(shù)據(jù)排列格式如表1時,使用上述SAS程序語句就可以直接建立適合SAS進行統(tǒng)計分析的數(shù)據(jù)集。

    2.2.2 對20個辣椒品種同時進行方差分析

    在Excel中只能對20組辣椒澇漬脅迫實驗數(shù)據(jù)逐個進行方差分析,在SAS中可以一次性對20組數(shù)據(jù)同時進行分析。使用SAS的過程步PROC對數(shù)據(jù)進行分析,同時進行DUNCAN多重比較,顯著性水平0.05。

    核心代碼及注釋如下:

    PROC ANOVA data=WORK.Yelvsu; /*調(diào)用上一步建立的數(shù)據(jù)集Yelvsu*/

    class TRT;

    model BRD1 BRD2 BRD3 BRD4 BRD5 BRD6 BRD7 BRD8 BRD9 BRD10 BRD11 BRD12 BRD13

    BRD14 BRD15 BRD16 BRD17 BRD18 BRD19 BRD20=TRT; /*對20個品種進行方差分析*/

    means TRT / HOVTEST=LEVENE; /*進行方差齊性檢驗*/

    means trt / DUNCAN alpha=0.05; /*進行DUNCAN多重比較*/

    RUN;

    對本試驗而言,辣椒葉片葉綠素SPAD值是因變量,進行澇漬脅迫使得辣椒葉片葉綠素含量在處理間出現(xiàn)差異,那么脅迫處理為“因素”,加上對照,此“因素”具有4個“水平”,在單因素方差分析中該“因素”為自變量,自變量是唯一的。在SAS語句中,通過“class TRT;”語句將TRT設(shè)置為水平,對4個水平進行方差分析和比較。按照蓋鈞鎰[9]對實驗統(tǒng)計方法的描述,本試驗中對1個辣椒品種進行4個水平方差分析為“組內(nèi)觀察值數(shù)目相等的單向分組資料的方差分析”,那么上述SAS語句“model BRD1 BRD2 BRD3 BRD4 BRD5 BRD6 BRD7 BRD8 BRD9 BRD10 BRD11 BRD12 BRD13 BRD14 BRD15 BRD16 BRD17 BRD18 BRD19 BRD20=TRT; ”可以一次性完成20組單向分組資料的方差分析,這樣就一次性對20個品種在澇漬脅迫3、6和9 d后葉片的葉綠素含量差異進行了顯著性分析,而不用重復(fù)20次相同的工作?!癿eans TRT / HOVTEST=LEVENE;”表示在方差分析之前先進行方差的齊性檢驗。“means trt / DUNCAN alpha=0.05;”表示執(zhí)行唐肯氏多范圍檢定(Duncan’s Multiple-Range Test),當(dāng)進行顯著性水平為0.05的多重比較時,“alpha=0.05”參數(shù)可以省略。

    2.3 比較辣椒品種間的差異顯著性

    2.3.1 建立數(shù)據(jù)集

    比較20個辣椒品種在3個處理與對照之間的葉片葉綠素SPAD含量差異顯著性,需要重新建立數(shù)據(jù)集。

    核心語句及注釋如下:

    DATA Work.yelvsu1; /*輸入所有品種葉綠素數(shù)據(jù)*/

    DO trt=1 TO 4; /*4個處理數(shù),保留*/

    DO t=1 TO 15; /*15個重復(fù)數(shù),不寫入數(shù)據(jù)庫*/

    DO brd=1 TO 20; /*20個品種數(shù),保留*/

    INPUT spad@@;

    OUTPUT;

    END;

    END;

    END;

    DROP t;

    CARDS;

    38 32 32 42 37 34 29 32 32 38 37 37 35 44 40 35 28 30 37 30

    40 35 33 38 40 36 37 29 39 38 36 36 38 41 42 39 31 33 37 34

    40 34 35 39 32 36 27 27 43 34 39 35 37 38 35 39 38 33 33 33

    39 39 29 35 37 39 31 29 36 33 42 36 34 36 38 39 39 37 34 31

    41 36 38 37 32 33 33 29 31 25 38 37 39 37 37 40 28 28 36 34

    ……(省略部分類似數(shù)據(jù))

    ;

    RUN;

    由于是比較品種間在不同澇漬脅迫條件下葉片葉綠素含量的差異,那么這里可以把20個品種看成一個因素的20個水平,所以需要將所有20個品種葉片葉綠素SPAD值放到SAS數(shù)據(jù)集單獨一列中,以便進行處理間的方差比較,在SAS語句中通過添加一個循環(huán)來實現(xiàn),具體語句為“DO brd=1 TO 20;”。執(zhí)行之后數(shù)據(jù)集yelvsu1會產(chǎn)生trt,brd和spad共3列數(shù)據(jù)。

    2.3.2進行方差分析

    比較澇漬脅迫下辣椒品種間葉片葉綠素含量差異,并進行DUNCAN多重比較。

    核心代碼和注釋如下:

    PROC anova data=Work.yelvsu1; /*調(diào)用上一步建立的數(shù)據(jù)集Yelvsu1*/

    class BRD;

    model spad=BRD;

    means BRD / HOVTEST=LEVENE;

    means brd / DUNCAN; /*默認(rèn)顯著性水平為0.05*/

    RUN;

    “class BRD;”語句是把所有20個辣椒品種看作一個因素的20個水平,比較20個辣椒品種間在澇漬脅迫下葉綠素含量是否存在差異顯著性。如果把所有的辣椒品種看成一個整體,比較4個處理之間在澇漬脅迫下葉片葉綠素含量是否有差異,在程序中把class的值修改為TRT即可。在進行多重比較時,不寫alpha參數(shù)值則默認(rèn)顯著性水平為0.05。如果想同時比較不同品種間和不同處理間辣椒葉片SPAD值差異,按照蓋鈞鎰[9]對實驗統(tǒng)計方法的描述,此類分析屬于“組內(nèi)有重復(fù)觀察值的兩向分組資料的方差分析”[9],即多因素方差分析。與楊孔雀等[10]使用的方法不同,SAS系統(tǒng)PROC步的anova分析主要研究單個因素對因變量的影響[11],多因素方差分析建議使用PROC步的GLM模型進行參數(shù)估計。

    3 討論

    在SAS數(shù)據(jù)集的建立方法上,除了本文所述采用DATA步的input函數(shù)方法外,還可以使用從外部導(dǎo)入數(shù)據(jù)文件的方法,但是在導(dǎo)入數(shù)據(jù)之前需要將數(shù)據(jù)整理成SAS所需要的格式。在農(nóng)業(yè)生物統(tǒng)計中常常需要對大量的數(shù)據(jù)進行統(tǒng)計分析,在本篇文章中有1500個葉綠素SPAD值數(shù)據(jù),在將該數(shù)據(jù)在Excel中整理成SAS所需格式時,由于數(shù)據(jù)超長,1500個數(shù)據(jù)就需要占用1500行,超過一個屏幕的高度,這樣在整理數(shù)據(jù)容易造成數(shù)據(jù)的丟失、重疊等問題。當(dāng)使用SAS分析較大數(shù)量的數(shù)據(jù)時建議使用DATA步來建立數(shù)據(jù)集,盡量不要采用人工整理數(shù)據(jù)然后導(dǎo)入SAS的方式,避免人為錯誤導(dǎo)致分析結(jié)果的失真。

    4 小結(jié)

    使用SAS同時對20個辣椒品種進行差異顯著性分析,明顯比用Excel逐個進行方差分析更為方便快捷,并且能用字母標(biāo)注法顯示各個品種在各處理間的顯著性差異。在SAS系統(tǒng)中Work數(shù)據(jù)庫是一個臨時庫,當(dāng)SAS系統(tǒng)關(guān)閉時,臨時數(shù)據(jù)庫中的內(nèi)容會被自動清空,所以在實際使用過程中建議單獨建立一個永久庫,這樣避免每次進行數(shù)據(jù)分析時都要重新建立數(shù)據(jù)集的麻煩。DATA步INPUT語句中的“@@”符號表示順序依次連續(xù)讀入數(shù)據(jù),如果沒有此符號則按照行讀入數(shù)據(jù)。在實際SAS分析數(shù)據(jù)時我們往往遇到組內(nèi)觀察值數(shù)據(jù)不完整的情況,這時我們可以將空格數(shù)據(jù)用英文狀態(tài)的“.”來代替,但是過程步PROC這里不能使用ANOVA,ANOVA主要對平衡實驗設(shè)計的數(shù)據(jù)執(zhí)行變異數(shù)分析,對于不平衡試驗設(shè)計數(shù)據(jù)則建議用一般線性模型GLM進行變異數(shù)分析。

    [1]彭紹英.世界統(tǒng)計與分析全才 SAS 系統(tǒng)應(yīng)用指南[M].北京: 希望電子出版社, 2000: 1-868.

    [2]顏亭玉, 杜曉林.SAS聚類分析過程在農(nóng)業(yè)試驗統(tǒng)計分析中的應(yīng)用[J].北京農(nóng)學(xué)院學(xué)報, 2009,24(1): 73-75.

    [3]周 倩, 張晉昕.含缺失值的重復(fù)測量資料分析在SPSS和SAS中的實現(xiàn)[J].循證醫(yī)學(xué), 2013,13(2): 120-123.

    [4]林 潔, 孫志明.SAS、PASS、Stata三種常用軟件樣本量計算方法及結(jié)果差異的比較[J].中國醫(yī)藥導(dǎo)報, 2015,12(18): 133-137.

    [5]詹秋文.Excel和SAS在生物統(tǒng)計學(xué)的應(yīng)用比較[J].生物學(xué)雜志, 2009,26(1): 74-76.

    [6]霍志軍, 李菊艷, 潘曉琳.Excel在農(nóng)業(yè)生物統(tǒng)計分析中的應(yīng)用[J].現(xiàn)代化農(nóng)業(yè), 2003(9): 28-30.

    [7]范 平, 崔黨群, 詹克慧, 等.Excel軟件在生物統(tǒng)計實驗教學(xué)中的綜合開發(fā)應(yīng)用[J].實驗技術(shù)與管理, 2003,20(2): 65-69.

    [8]尹冬梅, 管志勇, 陳素梅, 等.菊花及其近緣種屬植物耐澇評價體系建立及耐澇性鑒定[J].植物遺傳資源學(xué)報, 2009, 10(3): 399-404.

    [9]蓋鈞鎰.試驗統(tǒng)計方法[M].北京: 中國農(nóng)業(yè)出版社, 2000: 100-127.

    [10]楊孔雀, 孫占育, 蒙惠軍.用SAS軟件進行組內(nèi)有重復(fù)的兩向分組資料的方差分析[J].河北農(nóng)業(yè)科學(xué), 2010,14(10): 165-166.

    [11]阮 敬.SAS統(tǒng)計分析從入門到精通[M].北京: 人民郵電出版社, 2009: 99-122.

    猜你喜歡
    語句葉綠素辣椒
    原來,你還是這樣的辣椒
    你的辣椒結(jié)出果實了嗎?
    辣椒也瘋狂
    提取葉綠素
    重點:語句銜接
    揀辣椒
    中外文摘(2020年9期)2020-06-01 13:47:56
    桃樹葉綠素含量與SPAD值呈極顯著正相關(guān)
    精彩語句
    葉綠素家族概述
    由松針制取三種葉綠素鈉鹽及其穩(wěn)定性的研究
    欧洲精品卡2卡3卡4卡5卡区| 亚洲av.av天堂| 国产成年人精品一区二区| av天堂中文字幕网| 亚洲第一电影网av| 久久久久九九精品影院| 国产精品爽爽va在线观看网站| 国产一区二区在线观看日韩| 欧美日韩瑟瑟在线播放| 成年女人毛片免费观看观看9| 一级毛片久久久久久久久女| 在线观看一区二区三区| 国产高清不卡午夜福利| 午夜久久久久精精品| 精品一区二区三区视频在线| 色哟哟·www| 亚洲精品在线观看二区| 国产精品一区二区三区四区免费观看 | 欧美成人免费av一区二区三区| 听说在线观看完整版免费高清| 精华霜和精华液先用哪个| 免费看a级黄色片| 亚洲第一电影网av| 99精品久久久久人妻精品| 乱人视频在线观看| 尾随美女入室| 在线看三级毛片| 国产精品爽爽va在线观看网站| 午夜a级毛片| av在线观看视频网站免费| 久久欧美精品欧美久久欧美| 国产精品亚洲一级av第二区| av专区在线播放| 国产综合懂色| 老女人水多毛片| 久久精品国产99精品国产亚洲性色| 精品人妻熟女av久视频| 日韩欧美精品免费久久| 国内毛片毛片毛片毛片毛片| 有码 亚洲区| 欧美bdsm另类| 男人舔奶头视频| 欧美日韩综合久久久久久 | 波多野结衣巨乳人妻| 观看免费一级毛片| 两人在一起打扑克的视频| 中国美女看黄片| 69av精品久久久久久| 久久久精品欧美日韩精品| 国产一区二区在线av高清观看| 99热只有精品国产| 国产不卡一卡二| 神马国产精品三级电影在线观看| 婷婷亚洲欧美| 天堂av国产一区二区熟女人妻| 麻豆av噜噜一区二区三区| 男插女下体视频免费在线播放| 久久6这里有精品| 国产男人的电影天堂91| 国产伦精品一区二区三区视频9| 嫩草影视91久久| 全区人妻精品视频| 国产 一区 欧美 日韩| 女人被狂操c到高潮| 别揉我奶头~嗯~啊~动态视频| 99久久无色码亚洲精品果冻| 听说在线观看完整版免费高清| 国产男人的电影天堂91| 亚洲国产欧洲综合997久久,| 五月玫瑰六月丁香| 日韩精品有码人妻一区| 欧美+亚洲+日韩+国产| 91久久精品电影网| 在线观看午夜福利视频| 噜噜噜噜噜久久久久久91| 长腿黑丝高跟| 长腿黑丝高跟| 99在线人妻在线中文字幕| 少妇人妻一区二区三区视频| 12—13女人毛片做爰片一| 一级黄色大片毛片| 欧美绝顶高潮抽搐喷水| 91在线观看av| 色综合婷婷激情| 午夜日韩欧美国产| av在线老鸭窝| 精品人妻一区二区三区麻豆 | 国产私拍福利视频在线观看| 日本 av在线| www日本黄色视频网| 精品人妻视频免费看| 亚洲av中文av极速乱 | 亚洲欧美日韩卡通动漫| 69人妻影院| 成人特级av手机在线观看| АⅤ资源中文在线天堂| 国产精品一及| 自拍偷自拍亚洲精品老妇| 久久热精品热| 成年女人看的毛片在线观看| 九九在线视频观看精品| 成人高潮视频无遮挡免费网站| 国产探花在线观看一区二区| 亚洲在线自拍视频| 中文字幕人妻熟人妻熟丝袜美| 欧洲精品卡2卡3卡4卡5卡区| 老司机午夜福利在线观看视频| 美女 人体艺术 gogo| 久9热在线精品视频| 99国产精品一区二区蜜桃av| 乱人视频在线观看| 国内毛片毛片毛片毛片毛片| 成人特级av手机在线观看| 亚洲性夜色夜夜综合| 观看美女的网站| 男女之事视频高清在线观看| 哪里可以看免费的av片| 日日撸夜夜添| 亚洲av.av天堂| 欧美黑人欧美精品刺激| 免费不卡的大黄色大毛片视频在线观看 | 日日啪夜夜撸| 日本黄大片高清| 中文字幕久久专区| 亚洲中文日韩欧美视频| a级一级毛片免费在线观看| 精品国内亚洲2022精品成人| 欧美日韩亚洲国产一区二区在线观看| 欧美3d第一页| 夜夜夜夜夜久久久久| 欧美xxxx性猛交bbbb| 真人一进一出gif抽搐免费| 日本 欧美在线| 色综合站精品国产| 97超视频在线观看视频| 国产91精品成人一区二区三区| 亚洲乱码一区二区免费版| 婷婷六月久久综合丁香| 2021天堂中文幕一二区在线观| 1000部很黄的大片| 国内精品美女久久久久久| av在线蜜桃| 久久欧美精品欧美久久欧美| 成人国产一区最新在线观看| 亚洲成人久久性| 欧美三级亚洲精品| 久久精品91蜜桃| 欧美+亚洲+日韩+国产| 亚洲精品一区av在线观看| 午夜免费成人在线视频| 人人妻,人人澡人人爽秒播| 成年女人看的毛片在线观看| 日本免费一区二区三区高清不卡| av在线观看视频网站免费| 国产精品av视频在线免费观看| 最新中文字幕久久久久| 黄片wwwwww| 日韩中字成人| 97超视频在线观看视频| 99久久精品热视频| 国产蜜桃级精品一区二区三区| 99久久精品热视频| 五月玫瑰六月丁香| 亚洲精品国产成人久久av| 国产精品精品国产色婷婷| 国产熟女欧美一区二区| 看免费成人av毛片| 黄色一级大片看看| 伦理电影大哥的女人| 亚洲色图av天堂| 动漫黄色视频在线观看| 亚洲国产高清在线一区二区三| 女同久久另类99精品国产91| 一个人看的www免费观看视频| 国产视频一区二区在线看| 精品一区二区三区视频在线观看免费| 男女啪啪激烈高潮av片| 亚洲精品一区av在线观看| 欧美日本亚洲视频在线播放| 久久精品久久久久久噜噜老黄 | av.在线天堂| 12—13女人毛片做爰片一| 亚洲真实伦在线观看| 日韩精品有码人妻一区| 美女高潮喷水抽搐中文字幕| 国产一区二区三区在线臀色熟女| 亚洲精华国产精华精| 99在线视频只有这里精品首页| 成年免费大片在线观看| 欧美极品一区二区三区四区| 中文字幕人妻熟人妻熟丝袜美| 啦啦啦观看免费观看视频高清| 窝窝影院91人妻| 18+在线观看网站| 国产亚洲精品综合一区在线观看| 真人做人爱边吃奶动态| 搡老熟女国产l中国老女人| 午夜福利在线观看免费完整高清在 | 最好的美女福利视频网| 精品久久久久久,| 日韩精品有码人妻一区| 亚洲18禁久久av| 深爱激情五月婷婷| 国产视频一区二区在线看| 国产精品女同一区二区软件 | 国产高清三级在线| 中文字幕av在线有码专区| av专区在线播放| 成人高潮视频无遮挡免费网站| 国产高清不卡午夜福利| 在线观看av片永久免费下载| 国产午夜精品久久久久久一区二区三区 | 两个人视频免费观看高清| www.www免费av| 可以在线观看的亚洲视频| 亚洲欧美精品综合久久99| av在线老鸭窝| 日日摸夜夜添夜夜添小说| 黄色配什么色好看| 亚洲中文字幕一区二区三区有码在线看| 午夜福利在线在线| 久久这里只有精品中国| 舔av片在线| 欧美黑人巨大hd| 成人性生交大片免费视频hd| 亚洲美女搞黄在线观看 | 日韩在线高清观看一区二区三区 | 精品久久久久久久久av| 国产一区二区三区在线臀色熟女| 在线观看午夜福利视频| 国模一区二区三区四区视频| 国产美女午夜福利| 啦啦啦韩国在线观看视频| 久久精品国产亚洲av涩爱 | 91在线精品国自产拍蜜月| 69人妻影院| 我要搜黄色片| 国产真实乱freesex| 蜜桃亚洲精品一区二区三区| 亚洲三级黄色毛片| 伦理电影大哥的女人| av天堂在线播放| 日韩欧美一区二区三区在线观看| 成人特级av手机在线观看| 97热精品久久久久久| 欧美xxxx性猛交bbbb| 亚洲天堂国产精品一区在线| 久久久久久久精品吃奶| 蜜桃亚洲精品一区二区三区| 色综合婷婷激情| 亚洲av.av天堂| 国产69精品久久久久777片| 国产午夜福利久久久久久| 美女xxoo啪啪120秒动态图| 国产成人a区在线观看| 3wmmmm亚洲av在线观看| 久久久久性生活片| 日韩欧美 国产精品| 天天躁日日操中文字幕| 黄色配什么色好看| 制服丝袜大香蕉在线| 国产精品电影一区二区三区| 国产蜜桃级精品一区二区三区| 国产视频内射| 欧美一区二区精品小视频在线| 欧美xxxx黑人xx丫x性爽| 国产精品不卡视频一区二区| 欧美日韩中文字幕国产精品一区二区三区| 级片在线观看| 黄色女人牲交| 亚洲av电影不卡..在线观看| 成人三级黄色视频| 少妇人妻精品综合一区二区 | 成人亚洲精品av一区二区| 此物有八面人人有两片| 真人一进一出gif抽搐免费| 人妻制服诱惑在线中文字幕| 极品教师在线免费播放| 啦啦啦韩国在线观看视频| 久久欧美精品欧美久久欧美| 99久久无色码亚洲精品果冻| 日韩精品中文字幕看吧| 91在线精品国自产拍蜜月| 赤兔流量卡办理| 精品午夜福利在线看| 精品99又大又爽又粗少妇毛片 | 精品久久久久久久久久久久久| 国产三级在线视频| 在线看三级毛片| 99久久久亚洲精品蜜臀av| 久久久久国内视频| 国产精品国产高清国产av| 听说在线观看完整版免费高清| 欧美一级a爱片免费观看看| 日本 av在线| 国产精品伦人一区二区| 国产真实伦视频高清在线观看 | 久久亚洲精品不卡| 亚洲avbb在线观看| 在线观看一区二区三区| 欧美极品一区二区三区四区| 午夜福利欧美成人| 成年女人毛片免费观看观看9| 成人国产一区最新在线观看| bbb黄色大片| 日本-黄色视频高清免费观看| 久久天躁狠狠躁夜夜2o2o| 亚洲最大成人中文| 91在线精品国自产拍蜜月| 在线观看舔阴道视频| 国产精品日韩av在线免费观看| 久久久精品大字幕| 99热只有精品国产| 国产成人a区在线观看| www.色视频.com| 欧美另类亚洲清纯唯美| 最新中文字幕久久久久| 亚洲男人的天堂狠狠| 精品人妻视频免费看| 亚洲欧美激情综合另类| 自拍偷自拍亚洲精品老妇| 国产成人av教育| 99久久久亚洲精品蜜臀av| 国内精品宾馆在线| av在线蜜桃| 嫁个100分男人电影在线观看| 国产毛片a区久久久久| 国产精品精品国产色婷婷| 少妇的逼水好多| 国产精品三级大全| 成人精品一区二区免费| 狂野欧美激情性xxxx在线观看| 国产探花在线观看一区二区| 有码 亚洲区| 中文字幕久久专区| 国产黄a三级三级三级人| 成人美女网站在线观看视频| 十八禁网站免费在线| 国产精品一区二区免费欧美| 99热只有精品国产| 黄色欧美视频在线观看| 亚洲内射少妇av| 日韩大尺度精品在线看网址| 有码 亚洲区| 亚洲三级黄色毛片| 亚洲国产欧美人成| 麻豆一二三区av精品| 亚洲av成人av| 国产欧美日韩精品一区二区| 精品久久国产蜜桃| 午夜日韩欧美国产| 我要看日韩黄色一级片| 特级一级黄色大片| 久久久久久久精品吃奶| 波野结衣二区三区在线| 一本久久中文字幕| 动漫黄色视频在线观看| 亚洲国产色片| 国内少妇人妻偷人精品xxx网站| 午夜精品久久久久久毛片777| 国产午夜福利久久久久久| 国产亚洲av嫩草精品影院| 日韩中文字幕欧美一区二区| 久久国内精品自在自线图片| 国产一区二区三区视频了| 桃红色精品国产亚洲av| 午夜福利在线观看免费完整高清在 | 国产不卡一卡二| 午夜福利成人在线免费观看| 久久精品久久久久久噜噜老黄 | 午夜福利在线观看免费完整高清在 | 在线国产一区二区在线| 偷拍熟女少妇极品色| 日日摸夜夜添夜夜添av毛片 | 看十八女毛片水多多多| 日韩欧美精品免费久久| 精品一区二区三区视频在线观看免费| av女优亚洲男人天堂| 三级国产精品欧美在线观看| 简卡轻食公司| 婷婷精品国产亚洲av在线| 一级黄片播放器| 人人妻,人人澡人人爽秒播| 国产视频内射| eeuss影院久久| 精品久久久久久成人av| 日韩欧美一区二区三区在线观看| 日本一本二区三区精品| 精品久久国产蜜桃| 中文资源天堂在线| 女人十人毛片免费观看3o分钟| 色综合色国产| 日韩欧美精品免费久久| 国产精品一区二区性色av| av女优亚洲男人天堂| 日韩人妻高清精品专区| 日本撒尿小便嘘嘘汇集6| 亚洲乱码一区二区免费版| 亚洲精品456在线播放app | 日本一本二区三区精品| 麻豆国产97在线/欧美| 少妇的逼水好多| 国产精品精品国产色婷婷| 欧美人与善性xxx| 少妇裸体淫交视频免费看高清| 国产精品乱码一区二三区的特点| 精品久久久久久久末码| 不卡一级毛片| 国产成人aa在线观看| 人人妻,人人澡人人爽秒播| 干丝袜人妻中文字幕| 精品久久久久久久久亚洲 | 变态另类丝袜制服| 欧美日本视频| 亚洲七黄色美女视频| 亚洲色图av天堂| 午夜老司机福利剧场| 国产极品精品免费视频能看的| av女优亚洲男人天堂| 成人国产一区最新在线观看| 一级av片app| 午夜福利欧美成人| 国产aⅴ精品一区二区三区波| 亚洲国产欧洲综合997久久,| 国产精品人妻久久久久久| 国产精品精品国产色婷婷| 国产成人一区二区在线| 免费人成视频x8x8入口观看| av.在线天堂| 在线免费十八禁| 久久久午夜欧美精品| 国产亚洲欧美98| 3wmmmm亚洲av在线观看| 久久人妻av系列| 少妇的逼好多水| 成人精品一区二区免费| а√天堂www在线а√下载| 亚洲经典国产精华液单| 国产精品嫩草影院av在线观看 | 午夜福利视频1000在线观看| 91午夜精品亚洲一区二区三区 | www日本黄色视频网| 直男gayav资源| 精品午夜福利在线看| 少妇的逼水好多| 男女之事视频高清在线观看| 国产老妇女一区| 麻豆成人午夜福利视频| 国产av在哪里看| 人妻久久中文字幕网| 国产蜜桃级精品一区二区三区| 又黄又爽又免费观看的视频| 嫩草影院新地址| .国产精品久久| 精品欧美国产一区二区三| 国产三级中文精品| 久久精品国产亚洲网站| 成人国产麻豆网| 日本一本二区三区精品| 亚洲狠狠婷婷综合久久图片| 白带黄色成豆腐渣| 1024手机看黄色片| 色在线成人网| 国产在视频线在精品| 国产国拍精品亚洲av在线观看| 中文字幕高清在线视频| 91久久精品国产一区二区成人| 成人特级黄色片久久久久久久| 可以在线观看的亚洲视频| 午夜福利18| 国产精品永久免费网站| 女生性感内裤真人,穿戴方法视频| 国产伦在线观看视频一区| 欧美zozozo另类| 欧美三级亚洲精品| 精品无人区乱码1区二区| 欧美bdsm另类| 中文字幕久久专区| 超碰av人人做人人爽久久| 免费无遮挡裸体视频| 久久午夜亚洲精品久久| 色在线成人网| 国产aⅴ精品一区二区三区波| 色综合站精品国产| 欧美不卡视频在线免费观看| 亚洲最大成人手机在线| 99国产极品粉嫩在线观看| 亚洲久久久久久中文字幕| 亚洲国产欧洲综合997久久,| 国产精品99久久久久久久久| 一本久久中文字幕| 日本五十路高清| 久久精品国产亚洲av天美| 国产精品一区二区三区四区久久| 国产亚洲精品av在线| 亚洲五月天丁香| 欧美一区二区精品小视频在线| 免费看日本二区| 丰满人妻一区二区三区视频av| 亚洲中文日韩欧美视频| 性插视频无遮挡在线免费观看| 亚洲成人久久性| 亚洲精品456在线播放app | 亚洲国产日韩欧美精品在线观看| 亚洲国产精品成人综合色| 我的老师免费观看完整版| 女人被狂操c到高潮| 丰满的人妻完整版| 看片在线看免费视频| 丝袜美腿在线中文| 成人毛片a级毛片在线播放| 欧美3d第一页| 久久久久久久久久黄片| 国产亚洲精品av在线| 欧美xxxx黑人xx丫x性爽| 国产乱人视频| 亚洲精品乱码久久久v下载方式| 九色成人免费人妻av| 亚洲最大成人中文| 欧美色视频一区免费| 亚洲精品日韩av片在线观看| 久久人妻av系列| 俄罗斯特黄特色一大片| 看免费成人av毛片| 最近在线观看免费完整版| 搡老岳熟女国产| 国产高清三级在线| 日本免费a在线| 夜夜看夜夜爽夜夜摸| 午夜激情福利司机影院| 久久99热这里只有精品18| 观看免费一级毛片| 亚洲av一区综合| 亚洲电影在线观看av| 啦啦啦啦在线视频资源| 美女大奶头视频| 婷婷丁香在线五月| 99热网站在线观看| 欧美日本视频| 嫩草影院入口| 久久精品91蜜桃| 亚洲av免费高清在线观看| 国产一区二区三区av在线 | 无遮挡黄片免费观看| 中文字幕免费在线视频6| 国产精品人妻久久久影院| 日韩国内少妇激情av| 最近在线观看免费完整版| 黄片wwwwww| 99久久精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 老熟妇仑乱视频hdxx| 色av中文字幕| 国产黄a三级三级三级人| 亚洲色图av天堂| 午夜激情福利司机影院| 99久久中文字幕三级久久日本| 国产成人一区二区在线| 日韩亚洲欧美综合| 最新中文字幕久久久久| 午夜免费男女啪啪视频观看 | 一本精品99久久精品77| 超碰av人人做人人爽久久| 日本黄大片高清| 午夜福利成人在线免费观看| 久久精品综合一区二区三区| 欧美日韩综合久久久久久 | 精品人妻偷拍中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 噜噜噜噜噜久久久久久91| 在线天堂最新版资源| 99热这里只有是精品在线观看| 美女高潮喷水抽搐中文字幕| 欧美精品啪啪一区二区三区| 国产成人aa在线观看| www.色视频.com| 国产午夜精品论理片| 日韩亚洲欧美综合| 日本 av在线| 国产精品不卡视频一区二区| 久久九九热精品免费| 国产不卡一卡二| 精品国内亚洲2022精品成人| 亚洲精品影视一区二区三区av| 精品人妻1区二区| 成人美女网站在线观看视频| 亚洲国产精品成人综合色| 日日啪夜夜撸| 九色国产91popny在线| 一个人免费在线观看电影| 国内精品美女久久久久久| 99热这里只有精品一区| or卡值多少钱| 亚洲中文字幕一区二区三区有码在线看| 老司机福利观看| 狂野欧美激情性xxxx在线观看| 免费观看在线日韩| 婷婷丁香在线五月| 91久久精品国产一区二区三区| 一区二区三区激情视频| 欧美日韩亚洲国产一区二区在线观看| a级毛片a级免费在线| 精品一区二区三区视频在线观看免费| 亚洲av电影不卡..在线观看| 久久久久久久午夜电影| 一区二区三区激情视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲最大成人av| 欧美三级亚洲精品| 国产精品久久电影中文字幕| 18禁在线播放成人免费| 三级国产精品欧美在线观看| 欧美国产日韩亚洲一区| 一个人免费在线观看电影| 日韩人妻高清精品专区|