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

    SPSS進(jìn)行面板數(shù)據(jù)分析的基本方法

    2021-12-06 03:25:14李國柱李子寧李從欣
    內(nèi)蒙古統(tǒng)計(jì) 2021年4期
    關(guān)鍵詞:估計(jì)量組內(nèi)對話框

    ○ 文/李國柱 李子寧 李從欣

    文章梳理了面板數(shù)據(jù)尤其是短面板數(shù)據(jù)的估計(jì)方法,提出了用SPSS軟件進(jìn)行了混合回歸、組間估計(jì)、組內(nèi)估計(jì)、一階差分估計(jì)、最小二乘虛擬變量估計(jì)的操作流程,并與STATA軟件的估計(jì)結(jié)果進(jìn)行了對比。結(jié)果證實(shí),用SPSS估計(jì)面板數(shù)據(jù)可以得到與STATA完全相同的結(jié)果,拓展了SPSS軟件的使用范圍。

    SPSS是世界上應(yīng)用最廣泛的統(tǒng)計(jì)軟件之一,其特點(diǎn)是操作簡單、統(tǒng)計(jì)方法比較齊全,在視窗版軟件中得到廣大用戶的認(rèn)可。但SPSS軟件也有一個(gè)大的缺點(diǎn),其包括的統(tǒng)計(jì)方法主要是基于截面數(shù)據(jù)的統(tǒng)計(jì)方法,在時(shí)間序列數(shù)據(jù)方面僅包括自相關(guān)圖、互相關(guān)圖、專家建模器等,沒有單位根檢驗(yàn)、協(xié)整分析、向量自回歸模型等估計(jì)方法,而對于面板數(shù)據(jù)更是絲毫沒有涉及。面板數(shù)據(jù)由于可以解決遺漏變量問題、提供更多動(dòng)態(tài)行為信息以及使樣本容量增大,在實(shí)證研究中得到越來越廣泛的應(yīng)用。那么能否用SPSS進(jìn)行面板數(shù)據(jù)分析呢?雖然不能直接估計(jì),但只要熟知面板數(shù)據(jù)的估計(jì)方法,通過適當(dāng)?shù)牟僮?,SPSS也可以估計(jì)面板數(shù)據(jù)。

    一、SPSS輸入面板數(shù)據(jù)的方法

    SPSS默認(rèn)輸入的是截面數(shù)據(jù),當(dāng)輸入時(shí)間序列數(shù)據(jù)時(shí),必須首先通過“數(shù)據(jù)/定義日期”過程定義時(shí)間變量,然后通過“轉(zhuǎn)換/創(chuàng)建時(shí)間序列”功能創(chuàng)建一個(gè)新的時(shí)間序列。而對于面板數(shù)據(jù),可采用類似STATA軟件的輸入方式,除包括建模的變量外,還需要?jiǎng)?chuàng)建一個(gè)日期變量和截面變量。本文以2010-2017年我國31個(gè)?。ㄗ灾螀^(qū)、直轄市)的數(shù)據(jù)為例,說明面板數(shù)據(jù)的輸入方法以及各種估計(jì)方法。

    SPSS面板數(shù)據(jù)輸入格式

    假定要分析投資對國內(nèi)生產(chǎn)總值的影響,需要在SPSS中新建四個(gè)變量,其中province代表不同的省份(即截面單元),year代表年份(即時(shí)間單元),gdp代表國內(nèi)生產(chǎn)總值,invest代表投資。上圖是與STATA軟件相同的輸入格式。

    當(dāng)然,也可以把上圖中的數(shù)據(jù)格式按變量year進(jìn)行排序,此時(shí)數(shù)據(jù)中前31行為2010年各省份的數(shù)據(jù),32-62行是2011年各省份的數(shù)據(jù),以此類推。面板數(shù)據(jù)的兩種排列方式并不影響分析結(jié)果。

    二、混合回歸模型的參數(shù)估計(jì)

    如果從時(shí)間上看,不同個(gè)體之間不存在顯著差異,從截面上看,不同個(gè)體之間也不存在顯著差異,這種模型與一般回歸模型無本質(zhì)差異,只要滿足經(jīng)典回歸基本假定條件,就可以把所有數(shù)據(jù)放在一起進(jìn)行混合最小二乘估計(jì),估計(jì)量具有線性、無偏和有效性。

    SPSS進(jìn)行混合回歸很簡單,在菜單中依次選擇“分析/回歸/線性回歸”,在彈出的對話框中,因變量選擇gdp,解釋變量選擇invest,其他選擇默認(rèn)設(shè)置,點(diǎn)“確定”按鈕即可得到混合回歸結(jié)果,匯總?cè)绫?所示。

    表1 SPSS與STATA混合回歸結(jié)果

    由表1可以看出,SPSS混合回歸結(jié)果和STATA混合回歸結(jié)果相同,解釋變量invest的系數(shù)均為1.323,t統(tǒng)計(jì)量值均為25.535,p值均為0。

    三、固定效應(yīng)模型的參數(shù)估計(jì)

    面板數(shù)據(jù)模型可表示為:

    其中,yit是被解釋變量,xit是隨個(gè)體和時(shí)間而變化的解釋變量,zi是不隨時(shí)間而變的個(gè)體特征;不可觀測的隨機(jī)變量αi代表個(gè)體異質(zhì)性,隨個(gè)體變化,但不隨時(shí)間變化,即個(gè)體效應(yīng);β是回歸系數(shù),對于不同個(gè)體回歸系數(shù)β是相同的;uit是隨機(jī)誤差項(xiàng)。在以后的實(shí)證分析中,將用具體變量名代替yit和xit。

    若αi與某個(gè)解釋變量相關(guān),則稱為固定效應(yīng)模型;若αi與xit和zi均不相關(guān),則稱為隨機(jī)效應(yīng)模型。因此固定效應(yīng)模型和隨機(jī)效應(yīng)模型在表達(dá)式上是相同的,都是模型(1),只是根據(jù)個(gè)體效應(yīng)與解釋變量是否相關(guān)區(qū)分為固定效應(yīng)模型和隨機(jī)效應(yīng)模型。

    (一)組內(nèi)估計(jì)量

    對于固定效應(yīng)模型,由于αi與某個(gè)解釋變量相關(guān),故OLS不是一致的估計(jì)。解決方法之一是對面板數(shù)據(jù)中每個(gè)個(gè)體的觀測值針對時(shí)間求其平均值,利用離差變換消掉個(gè)體效應(yīng)αi,采用OLS方法估計(jì)模型回歸系數(shù)β。

    給定個(gè)i體,將模型(1)兩邊對時(shí)間求平均可得:

    將模型(1)減模型(2)可得:

    在采用SPSS計(jì)算組內(nèi)估計(jì)量時(shí),共需要三個(gè)步驟。第一步要計(jì)算每個(gè)地區(qū)的均值,在菜單中依次選擇“數(shù)據(jù)/分類匯總”,在彈出的對話框中,將provicne選入“分組變量”復(fù)選框,將gdp和invest選入“匯總變量”復(fù)選框。點(diǎn)擊“確定”按鈕后,SPSS數(shù)據(jù)窗口增加了兩個(gè)均值變量,均值變量名SPSS默認(rèn)為gdp_mean和invest_mean。第二步在菜單中依次選擇“轉(zhuǎn)換/計(jì)算變量”,在彈出的對話框中將目標(biāo)變量定義為gdpstar,數(shù)學(xué)表達(dá)式定義為“gdp-gdp_mean”,點(diǎn)擊“確定”后可得到中心化后的因變量值;按類似的方法可得到invest中心化后的變量investstar。第三步在菜單中依次選擇“分析/回歸/線性回歸”,在彈出的對話框中,因變量選擇gdpstar,解釋變量選擇investstar,其他選擇默認(rèn)設(shè)置,點(diǎn)“確定”按鈕即可得到組內(nèi)估計(jì)量結(jié)果,匯總?cè)绫?所示。

    表2 SPSS與STATA組內(nèi)估計(jì)量

    由表2可知,SPSS計(jì)算的組內(nèi)估計(jì)量中,解釋變量的系數(shù)為0.845,t統(tǒng)計(jì)量值為26.743,p值為0;STATA計(jì)算的組內(nèi)估計(jì)量中,解釋變量的系數(shù)為0.845,t統(tǒng)計(jì)量值為25.06,p值為0;除了t統(tǒng)計(jì)量值稍有差異外,系數(shù)估計(jì)值與p值完全相同。需要注意的是,兩種軟件對常數(shù)項(xiàng)的估計(jì)并不相同,原因在于兩種軟件對常數(shù)項(xiàng)的定義不同。在用SPSS進(jìn)行組內(nèi)估計(jì)時(shí),由于采用按時(shí)間平均然后中心化,個(gè)體效應(yīng)已經(jīng)被消掉了,也就是組內(nèi)回歸沒有常數(shù)項(xiàng),這也是SPSS組內(nèi)估計(jì)結(jié)果常數(shù)項(xiàng)的p值為1的原因。而STATA結(jié)果中的常數(shù)項(xiàng)實(shí)際上是個(gè)體效應(yīng)的均值,即,各個(gè)地區(qū)的個(gè)體效應(yīng)為,根據(jù)離差和恒等于零的性質(zhì),個(gè)體效應(yīng)之和等于零,因此STATA結(jié)果和SPSS結(jié)果并不矛盾。

    (二)一階差分估計(jì)量

    削掉個(gè)體效應(yīng)的另外一種方法是一階差分法。對于模型(1),取其滯后一期的關(guān)系式,得

    將模型(1)減模型(4)可得:

    SPSS只有定義時(shí)間序列變量后才能產(chǎn)生滯后變量和差分變量,而對于面板數(shù)據(jù),SPSS很難直接產(chǎn)生滯后變量和差分變量,可以借助STATA軟件產(chǎn)生差分變量,然后復(fù)制到SPSS軟件中,當(dāng)然也可以采用EXCEL,只不過稍微繁瑣一點(diǎn)。令gdp差分后的變量為dgdp,invest差分后的變量為dinvest。在菜單中依次選擇“分析/回歸/線性回歸”,在彈出的對話框中,因變量選擇dgdp,解釋變量選擇dinvest,然后點(diǎn)擊“選項(xiàng)”按鈕,將“包含常數(shù)項(xiàng)”前面的對鉤去掉,表示模型不包含常數(shù)項(xiàng);其他選擇默認(rèn)設(shè)置,點(diǎn)“確定”按鈕即可得到一階差分估計(jì)量結(jié)果,匯總?cè)绫?所示。

    表3 SPSS與STATA一階差分估計(jì)量

    由表3可以看出,SPSS與STATA一階差分估計(jì)結(jié)果完全相同。應(yīng)當(dāng)注意的是此處STATA軟件估計(jì)采用的是因變量差分對解釋變量差分直接回歸的方法,如果采用xtserial命令,系數(shù)估計(jì)值和SPSS方法完全相同,但由于該命令采用了穩(wěn)健標(biāo)準(zhǔn)誤,計(jì)算出來的t統(tǒng)計(jì)量值會(huì)有所不同。

    (三)最小二乘虛擬變量回歸

    對于固定效應(yīng)模型(1)中的個(gè)體固定效應(yīng)αi,將其視為個(gè)體i的截矩項(xiàng),即個(gè)體i的待估參數(shù)。對于n位個(gè)體的n個(gè)不同的截矩項(xiàng),可以通過在模型 (1)中引入n-1個(gè)虛擬變量來估計(jì)截矩項(xiàng),估計(jì)模型如下:

    常數(shù)項(xiàng)α1表示被遺漏的虛擬變量D1所對應(yīng)個(gè)體1的截矩項(xiàng),而個(gè)體i(i>1)的截矩項(xiàng)是α1+αi。

    對模型(6)采用OLS方法估計(jì),稱為“最小二乘虛擬變量法”。如果模型是正確設(shè)定的,且符合模型全部假定條件。則回歸系數(shù)估計(jì)量是無偏的、有效的、一致的估計(jì)。

    在SPSS主菜單中選擇“轉(zhuǎn)換/創(chuàng)建虛擬變量”,在彈出的“創(chuàng)建虛擬變量”對話框中,在“針對下列變量創(chuàng)建虛擬變量”框下選入變量province,在“主效應(yīng)虛擬變量-根名稱”下輸入“D”,表示虛擬變量,變量名以D開頭,點(diǎn)擊“確定”按鈕后,在SPSS中變量窗口即會(huì)出現(xiàn)31個(gè)虛擬變量。虛擬變量名和數(shù)據(jù)中省份的輸入順序有關(guān),上海對應(yīng)的虛擬變量為D_1,云南對應(yīng)的虛擬變量為D_2,以此類推。在采用虛擬變量回歸時(shí),如果回歸方程包括常數(shù)項(xiàng),為了避免“虛擬變量陷阱”,31個(gè)省份只能加入30個(gè)虛擬變量。STATA默認(rèn)將數(shù)據(jù)集中的第一個(gè)省份(本例為上海)作為比較的基礎(chǔ),因此為了和STATA估計(jì)結(jié)果進(jìn)行比較,在使用SPSS進(jìn)行最小二乘虛擬變量回歸時(shí),不包括虛擬變量D_1。在菜單中依次選擇“分析/回歸/線性回歸”,在彈出的對話框中,因變量選擇gdp,解釋變量依次選擇invest以及D_2-D_31,其他選擇默認(rèn)設(shè)置,點(diǎn)“確定”按鈕即可得到最小二乘虛擬變量回歸結(jié)果(如表4)。

    表4 SPSS與STATA最小二乘虛擬變量回歸

    由表4可以看出,除個(gè)別系數(shù)因計(jì)算精度和四舍五入原因稍有差異外,SPSS最小二乘虛擬變量回歸結(jié)果和STATA結(jié)果完全相同。表中的常數(shù)項(xiàng)18243.87代表上海市的截距,其他省份對應(yīng)的系數(shù)代表與上海的差距,如云南對應(yīng)的系數(shù)-1578.25代表云南比上海低1578.25。由表中p值可知只有浙江和上海沒有顯著性差異,江蘇、廣東、山東顯著高于上海。另外,最小二乘虛擬變量回歸在解釋變量invest的估計(jì)上和組內(nèi)估計(jì)量相同,均為0.845。

    四、隨機(jī)效應(yīng)模型的參數(shù)估計(jì)

    隨機(jī)效應(yīng)模型的參數(shù)估計(jì)方法包括廣義最小二乘法和組間估計(jì)法。雖然廣義最小二乘法最小效,但該方法涉及復(fù)雜的矩陣運(yùn)算和較多參數(shù)估計(jì),使用SPSS比較繁瑣,此處僅介紹組間估計(jì)方法。

    對于模型(1),如果對每個(gè)個(gè)體取時(shí)間平均值,得如下模型:

    對模型(7)使用OLS進(jìn)行參數(shù)估計(jì),得到的估計(jì)量稱為組間估計(jì)量,記為。由于包含了xit的信息,如果αi與解釋變量{xit,zi}相關(guān),則不一致。因此不能在固定效應(yīng)模型下使用組間估計(jì)法。組間估計(jì)法相當(dāng)于面板數(shù)據(jù)被壓縮為截面數(shù)據(jù)。

    在采用組內(nèi)估計(jì)量時(shí),已經(jīng)通過“分類匯總”功能得到各省在gdp和invest兩個(gè)變量上的均值gdp_mean和invest_mean,但我們并不能采用這兩個(gè)變量直接回歸,原因在于每個(gè)省的均值都有8個(gè)相同的值(從2010年至2018年),可以采用等距抽樣的方法使每個(gè)省只有一個(gè)均值。為了避免觀測順序打亂,首先在數(shù)據(jù)集增加一個(gè)編號(hào)變量(ID),打開需要添加ID號(hào)的數(shù)據(jù)集,選擇“文件/新建/語法”,打開彈出的“語法”窗口,并輸入以下語句:

    DATASET ACTIVATE 數(shù)據(jù)集名.

    點(diǎn)執(zhí)行按鈕。返回到原先的數(shù)據(jù)集窗口,則可以看到id號(hào)變量了。點(diǎn)擊“轉(zhuǎn)換/計(jì)算變量”,彈出計(jì)算變量對話框,在“目標(biāo)變量”的空白框中輸入新變量名newid,在“函數(shù)組”選項(xiàng)框中點(diǎn)擊all,在“函數(shù)和特殊變量”選項(xiàng)框中點(diǎn)擊mod,把mod送入“數(shù)字表達(dá)式”空白框中,把mod(?,?)的第一個(gè)問號(hào)改為id-4(由于前8個(gè)數(shù)值相同,用id減去1-8中的任意一個(gè)數(shù)均不會(huì)影響最終結(jié)果),第二個(gè)問號(hào)改為8。點(diǎn)擊“確定”按鈕,數(shù)據(jù)集中newid為0的即是要抽選的樣本單元。

    選擇“數(shù)據(jù)/選擇個(gè)案”,打開選擇個(gè)案對話框,選擇第二項(xiàng)“如果條件滿足”,單擊“如果”按鈕后彈出If對話框,將其中的條件設(shè)置為“newid=0”,即可得到等距抽樣的所有樣本,該樣本即為進(jìn)行組間估計(jì)的樣本。在該樣本范圍內(nèi)依次選擇“分析/回歸/線性回歸”,在彈出的對話框中,因變量選擇gdp_mean,解釋變量選擇invest_mean,其他選擇默認(rèn)設(shè)置,點(diǎn)“確定”按鈕即可得到組間估計(jì)回歸結(jié)果,匯總?cè)绫?所示。

    表5 SPSS與STATA組間估計(jì)量

    從表5可以看出,SPSS與STATA估計(jì)組間估計(jì)量的結(jié)果相同。

    五、結(jié)論與啟示

    從實(shí)證研究結(jié)果來看,在混合回歸、組內(nèi)估計(jì)、組間估計(jì)、一階差分估計(jì)、最小二乘虛擬變量回歸方面,SPSS可以得到與STATA軟件相同的結(jié)果。只不過STATA等軟件有內(nèi)嵌的命令,操作起來很簡單,而SPSS由于沒有內(nèi)嵌選項(xiàng),操作起來相對復(fù)雜而已。實(shí)際上只要理解各種統(tǒng)計(jì)方法的基本原理,以及熟悉SPSS的基本操作,就可以使用SPSS解決大多數(shù)沒有內(nèi)嵌的統(tǒng)計(jì)方法。

    猜你喜歡
    估計(jì)量組內(nèi)對話框
    用心說題 提高效率 培養(yǎng)能力
    正常恢復(fù)虛擬機(jī)
    Bootlace Worms’Secret etc.
    淺談估計(jì)量的優(yōu)良性標(biāo)準(zhǔn)
    合作學(xué)習(xí)組內(nèi)交流討論時(shí)間的遵循原則
    合作學(xué)習(xí)“組內(nèi)交流討論時(shí)間”注意問題
    基于配網(wǎng)先驗(yàn)信息的諧波狀態(tài)估計(jì)量測點(diǎn)最優(yōu)配置
    電測與儀表(2015年6期)2015-04-09 12:00:50
    合作學(xué)習(xí)組內(nèi)交流討論時(shí)間探究
    負(fù)極值指標(biāo)估計(jì)量的漸近性質(zhì)
    使用變異系數(shù)和Kurtosis系數(shù)的雙輔助變量的比估計(jì)
    日日啪夜夜撸| 男女下面进入的视频免费午夜| 熟女电影av网| 在线观看66精品国产| av在线亚洲专区| 麻豆成人午夜福利视频| 一区二区三区四区激情视频 | 亚洲乱码一区二区免费版| 亚洲精品日韩av片在线观看| 欧美bdsm另类| 国产成人福利小说| 日韩成人av中文字幕在线观看 | 特级一级黄色大片| 亚洲欧美日韩东京热| 免费人成在线观看视频色| 在线国产一区二区在线| 亚洲最大成人中文| 在线天堂最新版资源| 日韩精品青青久久久久久| 97碰自拍视频| 此物有八面人人有两片| 国产精品综合久久久久久久免费| 亚洲第一区二区三区不卡| 日本黄色视频三级网站网址| .国产精品久久| 欧美人与善性xxx| 中文资源天堂在线| 精品久久久久久久人妻蜜臀av| 亚洲成av人片在线播放无| 一个人看视频在线观看www免费| 欧美不卡视频在线免费观看| 精品福利观看| 变态另类成人亚洲欧美熟女| 国产老妇女一区| 麻豆久久精品国产亚洲av| 久久人人精品亚洲av| 精品久久久久久久久av| 三级毛片av免费| 啦啦啦观看免费观看视频高清| 欧美日韩综合久久久久久| 日本与韩国留学比较| 国产国拍精品亚洲av在线观看| 一边摸一边抽搐一进一小说| 国产亚洲精品av在线| 国产亚洲精品久久久久久毛片| 亚洲精品成人久久久久久| 亚洲国产欧洲综合997久久,| 床上黄色一级片| 国产精品人妻久久久影院| 最好的美女福利视频网| 午夜福利在线观看免费完整高清在 | 一边摸一边抽搐一进一小说| 久久久精品欧美日韩精品| 色综合站精品国产| 免费大片18禁| 少妇猛男粗大的猛烈进出视频 | 高清日韩中文字幕在线| 亚洲18禁久久av| 国产成人福利小说| 我要看日韩黄色一级片| 国产成人a∨麻豆精品| 美女cb高潮喷水在线观看| 国产麻豆成人av免费视频| 十八禁国产超污无遮挡网站| 日本成人三级电影网站| 男女下面进入的视频免费午夜| 两个人视频免费观看高清| 亚洲一区高清亚洲精品| 成人二区视频| 夜夜夜夜夜久久久久| 国产精品永久免费网站| 又黄又爽又刺激的免费视频.| 国产aⅴ精品一区二区三区波| 级片在线观看| 久久久久久久亚洲中文字幕| 欧美潮喷喷水| 成人三级黄色视频| 成人国产麻豆网| 精品午夜福利在线看| 久久精品影院6| 中文字幕熟女人妻在线| 99九九线精品视频在线观看视频| 一区二区三区免费毛片| 国产久久久一区二区三区| 久久久久久久久久成人| 亚洲成人精品中文字幕电影| 成人二区视频| 亚洲色图av天堂| 久久精品国产亚洲av香蕉五月| 岛国在线免费视频观看| 男人狂女人下面高潮的视频| 成人永久免费在线观看视频| 久久热精品热| 一区二区三区高清视频在线| 97超级碰碰碰精品色视频在线观看| 久99久视频精品免费| 尾随美女入室| 色综合色国产| 精品不卡国产一区二区三区| 午夜福利在线观看免费完整高清在 | 麻豆乱淫一区二区| 色av中文字幕| 国产色爽女视频免费观看| 国产激情偷乱视频一区二区| 看免费成人av毛片| 国产高潮美女av| 亚洲精品一区av在线观看| 婷婷六月久久综合丁香| 午夜免费观看性视频| 夜夜爽夜夜爽视频| h视频一区二区三区| 人妻制服诱惑在线中文字幕| 国产av精品麻豆| freevideosex欧美| 久久久久视频综合| 亚洲精品国产av成人精品| 精品亚洲成国产av| 日日啪夜夜撸| 综合色丁香网| 一个人免费看片子| 男女免费视频国产| 噜噜噜噜噜久久久久久91| a级毛色黄片| 国产伦在线观看视频一区| 黄片无遮挡物在线观看| 精品卡一卡二卡四卡免费| 国产深夜福利视频在线观看| 精品人妻偷拍中文字幕| 91在线精品国自产拍蜜月| 精品99又大又爽又粗少妇毛片| 亚洲欧美日韩东京热| 一区二区av电影网| 亚洲四区av| 国产色婷婷99| 午夜福利网站1000一区二区三区| 欧美另类一区| 亚洲人成网站在线播| 大陆偷拍与自拍| 中国三级夫妇交换| 国产亚洲一区二区精品| 最近中文字幕2019免费版| 亚洲av.av天堂| 午夜福利影视在线免费观看| 色5月婷婷丁香| 欧美bdsm另类| 欧美日韩av久久| 国产精品人妻久久久久久| 国产精品偷伦视频观看了| 亚洲av电影在线观看一区二区三区| 人妻系列 视频| 亚洲国产精品999| 久久精品久久久久久噜噜老黄| 亚洲婷婷狠狠爱综合网| 精品一品国产午夜福利视频| 国产一区二区三区av在线| 在线观看免费高清a一片| 观看免费一级毛片| 日日啪夜夜撸| 国模一区二区三区四区视频| 午夜免费观看性视频| 只有这里有精品99| 欧美日韩国产mv在线观看视频| 青春草视频在线免费观看| a级毛片在线看网站| 国产日韩一区二区三区精品不卡 | 三级国产精品片| 亚洲久久久国产精品| 国产伦精品一区二区三区视频9| 亚洲婷婷狠狠爱综合网| 少妇被粗大的猛进出69影院 | 久久精品国产亚洲av涩爱| 日本午夜av视频| 91在线精品国自产拍蜜月| 狠狠精品人妻久久久久久综合| 成年女人在线观看亚洲视频| av播播在线观看一区| 欧美少妇被猛烈插入视频| 大香蕉97超碰在线| 人人妻人人爽人人添夜夜欢视频 | 黄色视频在线播放观看不卡| 少妇高潮的动态图| 亚洲欧洲日产国产| 成年av动漫网址| 欧美区成人在线视频| 国产亚洲欧美精品永久| 嘟嘟电影网在线观看| 日韩熟女老妇一区二区性免费视频| 成人国产av品久久久| 80岁老熟妇乱子伦牲交| 久久av网站| 亚洲国产精品一区二区三区在线| 午夜av观看不卡| 欧美精品国产亚洲| 日本av手机在线免费观看| 成年人午夜在线观看视频| 成年av动漫网址| 99re6热这里在线精品视频| 天天操日日干夜夜撸| 国产成人免费观看mmmm| 久久青草综合色| 久久久久国产精品人妻一区二区| 日日啪夜夜爽| 91久久精品电影网| 美女福利国产在线| 少妇被粗大的猛进出69影院 | 九九久久精品国产亚洲av麻豆| 亚洲av男天堂| 欧美少妇被猛烈插入视频| 丰满少妇做爰视频| 欧美成人精品欧美一级黄| 夜夜爽夜夜爽视频| 大香蕉久久网| a级毛片在线看网站| 亚洲av男天堂| 视频中文字幕在线观看| 国产精品久久久久久精品电影小说| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线观看播放| 熟女电影av网| 少妇丰满av| 一级毛片久久久久久久久女| 色视频www国产| 久久精品熟女亚洲av麻豆精品| 99精国产麻豆久久婷婷| 日本wwww免费看| 丰满人妻一区二区三区视频av| 亚洲成色77777| 国产免费一级a男人的天堂| 成人国产麻豆网| 中文字幕av电影在线播放| 超碰97精品在线观看| 尾随美女入室| 国内少妇人妻偷人精品xxx网站| 男女国产视频网站| 欧美成人午夜免费资源| 久久99精品国语久久久| 免费观看的影片在线观看| 久久久欧美国产精品| 亚洲欧美成人综合另类久久久| 亚洲国产精品成人久久小说| 夫妻性生交免费视频一级片| 一级二级三级毛片免费看| 99久久人妻综合| 成人毛片a级毛片在线播放| 最近中文字幕高清免费大全6| 黄色欧美视频在线观看| 色婷婷av一区二区三区视频| 18禁动态无遮挡网站| 91午夜精品亚洲一区二区三区| 卡戴珊不雅视频在线播放| 中文字幕久久专区| 日韩伦理黄色片| 水蜜桃什么品种好| 国产精品久久久久久久电影| 视频中文字幕在线观看| 久久久欧美国产精品| 七月丁香在线播放| 纯流量卡能插随身wifi吗| 欧美区成人在线视频| 99久久精品热视频| 性色avwww在线观看| 亚洲第一av免费看| 国产成人精品久久久久久| 国产亚洲最大av| 国产精品嫩草影院av在线观看| 欧美+日韩+精品| 街头女战士在线观看网站| 9色porny在线观看| 亚洲国产欧美日韩在线播放 | 午夜福利视频精品| 亚洲精品一二三| 亚洲欧洲国产日韩| 免费在线观看成人毛片| 永久网站在线| 国精品久久久久久国模美| 精品一品国产午夜福利视频| 人人妻人人澡人人看| 丝袜脚勾引网站| 亚洲经典国产精华液单| 国产精品女同一区二区软件| 男人添女人高潮全过程视频| 亚洲综合色惰| 狠狠精品人妻久久久久久综合| 久久热精品热| 国内精品宾馆在线| 女性生殖器流出的白浆| 成人18禁高潮啪啪吃奶动态图 | 午夜视频国产福利| 精品午夜福利在线看| 国产精品蜜桃在线观看| 亚洲精品日韩在线中文字幕| 国产精品嫩草影院av在线观看| 欧美激情极品国产一区二区三区 | 最近的中文字幕免费完整| 女性生殖器流出的白浆| 久久青草综合色| 两个人免费观看高清视频 | 久久久精品免费免费高清| 日韩精品有码人妻一区| 另类精品久久| 日本与韩国留学比较| 国产黄色视频一区二区在线观看| 免费观看av网站的网址| 大片免费播放器 马上看| 亚洲欧美日韩东京热| 成人影院久久| a级一级毛片免费在线观看| 亚洲国产精品专区欧美| 18禁动态无遮挡网站| 免费高清在线观看视频在线观看| 国产精品久久久久成人av| 国产精品无大码| 一二三四中文在线观看免费高清| 观看免费一级毛片| 大话2 男鬼变身卡| 国产精品99久久久久久久久| 国产精品.久久久| 日韩av不卡免费在线播放| 亚洲内射少妇av| 啦啦啦啦在线视频资源| 成人免费观看视频高清| 五月玫瑰六月丁香| 免费在线观看成人毛片| 五月开心婷婷网| 亚洲av福利一区| 秋霞在线观看毛片| 国产精品一区二区三区四区免费观看| 亚洲国产精品专区欧美| 国产淫片久久久久久久久| 午夜免费男女啪啪视频观看| 3wmmmm亚洲av在线观看| 国产亚洲精品久久久com| 熟女av电影| 极品人妻少妇av视频| 免费看av在线观看网站| 在线观看免费日韩欧美大片 | 韩国高清视频一区二区三区| 久久青草综合色| 日本欧美视频一区| 亚洲av免费高清在线观看| 18+在线观看网站| 亚洲欧美精品专区久久| 91精品伊人久久大香线蕉| 久久综合国产亚洲精品| 亚洲人与动物交配视频| 三级经典国产精品| 黑人高潮一二区| 久久99热6这里只有精品| av视频免费观看在线观看| 91aial.com中文字幕在线观看| 国产日韩一区二区三区精品不卡 | 18禁裸乳无遮挡动漫免费视频| 99热全是精品| 国产精品.久久久| 最新中文字幕久久久久| 一级毛片我不卡| 在线播放无遮挡| 最近最新中文字幕免费大全7| 久久99热6这里只有精品| 国产高清不卡午夜福利| 有码 亚洲区| 狠狠精品人妻久久久久久综合| 亚洲av国产av综合av卡| 啦啦啦中文免费视频观看日本| av黄色大香蕉| 国产永久视频网站| 在线观看免费高清a一片| 丝袜喷水一区| 国产精品成人在线| 国产精品久久久久久av不卡| 人妻 亚洲 视频| 精品国产一区二区久久| 一区二区三区乱码不卡18| 亚洲无线观看免费| 午夜久久久在线观看| 国产精品久久久久久久久免| 美女脱内裤让男人舔精品视频| 国产亚洲午夜精品一区二区久久| 亚洲不卡免费看| 亚洲真实伦在线观看| av卡一久久| 国产又色又爽无遮挡免| 精品久久久久久久久亚洲| 美女福利国产在线| 赤兔流量卡办理| 只有这里有精品99| 王馨瑶露胸无遮挡在线观看| av专区在线播放| 免费观看a级毛片全部| 国模一区二区三区四区视频| 国产淫语在线视频| av播播在线观看一区| 亚洲va在线va天堂va国产| 一级毛片aaaaaa免费看小| 91aial.com中文字幕在线观看| 一区二区三区免费毛片| h视频一区二区三区| 久久精品国产自在天天线| 日本欧美国产在线视频| www.av在线官网国产| 亚洲综合色惰| 性色av一级| 在线观看av片永久免费下载| 97精品久久久久久久久久精品| 国产极品天堂在线| 男人狂女人下面高潮的视频| 精品国产一区二区久久| 亚洲av成人精品一区久久| 久久久国产一区二区| 日韩欧美一区视频在线观看 | 久久鲁丝午夜福利片| 久久国产亚洲av麻豆专区| 欧美最新免费一区二区三区| 在现免费观看毛片| 日日啪夜夜撸| 中文字幕亚洲精品专区| 日韩大片免费观看网站| 日韩精品免费视频一区二区三区 | 国产黄频视频在线观看| 如何舔出高潮| 一级毛片 在线播放| 亚洲第一av免费看| 狂野欧美激情性xxxx在线观看| 伦精品一区二区三区| 久久精品久久精品一区二区三区| av天堂中文字幕网| 亚洲国产色片| 中文欧美无线码| 中文字幕免费在线视频6| 亚洲av男天堂| 国产欧美日韩精品一区二区| 欧美日韩视频高清一区二区三区二| 精品人妻熟女毛片av久久网站| 亚洲色图综合在线观看| 国产av码专区亚洲av| 精品一区二区三区视频在线| 午夜日本视频在线| 亚洲国产成人一精品久久久| 亚洲一级一片aⅴ在线观看| 成年女人在线观看亚洲视频| 人人澡人人妻人| 精品酒店卫生间| 一区二区三区乱码不卡18| 午夜福利网站1000一区二区三区| 亚洲精品成人av观看孕妇| 哪个播放器可以免费观看大片| 熟妇人妻不卡中文字幕| 亚洲av.av天堂| 美女中出高潮动态图| 国产精品伦人一区二区| 亚洲国产色片| 亚洲精品一二三| 婷婷色综合大香蕉| 精品久久久久久电影网| 欧美精品国产亚洲| 成年人午夜在线观看视频| 少妇被粗大的猛进出69影院 | 日日摸夜夜添夜夜爱| 亚洲精品乱码久久久v下载方式| 蜜桃在线观看..| 亚洲精品色激情综合| 国产美女午夜福利| 亚洲精品一二三| 3wmmmm亚洲av在线观看| 日本黄色片子视频| 伊人久久精品亚洲午夜| 伦理电影免费视频| 晚上一个人看的免费电影| 日本与韩国留学比较| 日日撸夜夜添| 中文字幕av电影在线播放| 中文字幕亚洲精品专区| 中文天堂在线官网| 一区二区三区四区激情视频| 赤兔流量卡办理| 国产黄频视频在线观看| 亚洲成人手机| 中国美白少妇内射xxxbb| 日韩不卡一区二区三区视频在线| 国产无遮挡羞羞视频在线观看| 久久久久久久久久人人人人人人| 美女主播在线视频| 国产精品国产三级国产专区5o| 亚洲经典国产精华液单| 天美传媒精品一区二区| 高清毛片免费看| 一级,二级,三级黄色视频| 色婷婷久久久亚洲欧美| 99国产精品免费福利视频| 国内揄拍国产精品人妻在线| 国产一级毛片在线| 亚洲电影在线观看av| 毛片一级片免费看久久久久| 日韩电影二区| 大码成人一级视频| 91精品国产九色| 久久99一区二区三区| 免费观看的影片在线观看| 一级毛片黄色毛片免费观看视频| 九色成人免费人妻av| 久久国产精品男人的天堂亚洲 | 亚洲不卡免费看| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲av男天堂| 成年人免费黄色播放视频 | 午夜激情福利司机影院| 3wmmmm亚洲av在线观看| 热re99久久精品国产66热6| 久久久久久久久大av| 亚洲国产精品国产精品| 3wmmmm亚洲av在线观看| 男女无遮挡免费网站观看| 日韩中文字幕视频在线看片| 久久久久久久久久成人| 亚洲,一卡二卡三卡| 免费观看在线日韩| 在现免费观看毛片| av免费观看日本| 如何舔出高潮| 777米奇影视久久| 国产日韩欧美亚洲二区| 成人18禁高潮啪啪吃奶动态图 | 亚洲av综合色区一区| 99久久精品热视频| 中文字幕人妻熟人妻熟丝袜美| 日韩强制内射视频| 国产又色又爽无遮挡免| 王馨瑶露胸无遮挡在线观看| 国产av精品麻豆| 麻豆精品久久久久久蜜桃| 亚洲av中文av极速乱| 免费黄网站久久成人精品| 久久精品国产自在天天线| 国产淫语在线视频| 亚洲内射少妇av| 伦理电影大哥的女人| 最近中文字幕高清免费大全6| 亚洲国产精品一区二区三区在线| 一本一本综合久久| 美女主播在线视频| 国产欧美日韩一区二区三区在线 | 欧美另类一区| 少妇的逼水好多| 热re99久久国产66热| av有码第一页| 久久精品夜色国产| 久久99精品国语久久久| 少妇精品久久久久久久| 大片电影免费在线观看免费| 亚洲中文av在线| 国产黄片视频在线免费观看| 纵有疾风起免费观看全集完整版| 国内精品宾馆在线| 3wmmmm亚洲av在线观看| 亚洲av不卡在线观看| 国产精品三级大全| 精品一品国产午夜福利视频| 99久久综合免费| 岛国毛片在线播放| 18禁在线无遮挡免费观看视频| 日韩欧美一区视频在线观看 | 夜夜爽夜夜爽视频| 久久人人爽av亚洲精品天堂| 国产午夜精品一二区理论片| 亚洲va在线va天堂va国产| 老司机影院毛片| 天堂俺去俺来也www色官网| 久久久午夜欧美精品| 日本91视频免费播放| 亚洲精品乱久久久久久| 久久鲁丝午夜福利片| 国产av国产精品国产| 久久热精品热| 我要看黄色一级片免费的| 精华霜和精华液先用哪个| 女的被弄到高潮叫床怎么办| 黑人猛操日本美女一级片| 在线观看一区二区三区激情| 久久6这里有精品| 国产精品一区二区三区四区免费观看| 51国产日韩欧美| 日韩人妻高清精品专区| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩视频精品一区| 秋霞伦理黄片| 十分钟在线观看高清视频www | 波野结衣二区三区在线| 人妻系列 视频| 免费人妻精品一区二区三区视频| 欧美一级a爱片免费观看看| 国产精品免费大片| a级一级毛片免费在线观看| 能在线免费看毛片的网站| 国产精品一区二区在线观看99| 一本大道久久a久久精品| av福利片在线| 乱人伦中国视频| 观看美女的网站| 色94色欧美一区二区| 免费黄色在线免费观看| 国产 一区精品| 成年女人在线观看亚洲视频| 精品久久久久久久久av| 久久女婷五月综合色啪小说| 王馨瑶露胸无遮挡在线观看| 欧美日韩视频精品一区| 午夜视频国产福利| 9色porny在线观看| 成人毛片a级毛片在线播放| 建设人人有责人人尽责人人享有的| 免费不卡的大黄色大毛片视频在线观看| 精品亚洲乱码少妇综合久久| 亚洲精品久久午夜乱码| 国产黄片美女视频| 国产视频内射|