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

    地區(qū)生產(chǎn)總值與固定資產(chǎn)投資額的數(shù)據(jù)分析

    2017-09-12 03:09:04呂林靜
    關(guān)鍵詞:投資額離群線圖

    呂林靜

    (海南大學(xué)金融統(tǒng)計(jì)系, ???70228)

    地區(qū)生產(chǎn)總值與固定資產(chǎn)投資額的數(shù)據(jù)分析

    呂林靜

    (海南大學(xué)金融統(tǒng)計(jì)系, ???70228)

    運(yùn)用探索性數(shù)據(jù)分析的方法定量地分析各地全社會(huì)固定資產(chǎn)投資對(duì)于地區(qū)經(jīng)濟(jì)增長(zhǎng)的拉動(dòng)作用;分別運(yùn)用最小二乘法和三組耐抗線的方法進(jìn)行回歸擬合,并且通過(guò)將原始?xì)埐钸M(jìn)行不同尺度下的轉(zhuǎn)換、對(duì)每個(gè)樣本杠桿率進(jìn)行計(jì)算以及殘差箱線圖判斷異常值點(diǎn),綜合分析之后找出異常值,然后對(duì)去掉異常值點(diǎn)后的數(shù)據(jù)進(jìn)行最小二乘法的回歸,求出擬合直線,最后分別求出兩種方法下的均方誤差并對(duì)擬合模型的優(yōu)劣進(jìn)行評(píng)價(jià)。

    地區(qū)生產(chǎn)總值; OLS; 三組耐抗線; 異常值

    1 地區(qū)生產(chǎn)總值與固定資產(chǎn)投資概述

    國(guó)內(nèi)生產(chǎn)總值是指在一定時(shí)期內(nèi)一個(gè)國(guó)家或地區(qū)的經(jīng)濟(jì)中所生產(chǎn)出的全部最終產(chǎn)品和勞務(wù)的價(jià)值,是衡量國(guó)家經(jīng)濟(jì)狀況的最佳指標(biāo),不但可反映一個(gè)國(guó)家的經(jīng)濟(jì)表現(xiàn),更可以反映一國(guó)的國(guó)力與財(cái)富。實(shí)際統(tǒng)計(jì)中,一般以支出法為標(biāo)準(zhǔn)來(lái)計(jì)算國(guó)內(nèi)生產(chǎn)總值。在支出法中,投資支出直接構(gòu)成了國(guó)內(nèi)生產(chǎn)總值的一部分,而固定資產(chǎn)投資又是投資支出的重要組成部分。地區(qū)生產(chǎn)總值是指本地區(qū)所有常住單位在一定時(shí)期內(nèi)生產(chǎn)活動(dòng)的最終成果,是國(guó)內(nèi)生產(chǎn)總值在不同地區(qū)的具體化。因此,研究地區(qū)生產(chǎn)總值與固定資產(chǎn)投資的關(guān)系是有必要的,本文主要對(duì)地區(qū)生產(chǎn)總值和固定資產(chǎn)投資總額的關(guān)系進(jìn)行研究。

    2 數(shù)據(jù)來(lái)源

    為了定量地研究地區(qū)生產(chǎn)總值與固定資產(chǎn)投資額的關(guān)系,筆者選擇了地區(qū)生產(chǎn)總值以及各地區(qū)全社會(huì)固定資產(chǎn)投資額作為研究對(duì)象,選取的樣本為2014年中國(guó) 31 個(gè)省市的地區(qū)生產(chǎn)總值與全社會(huì)固定資產(chǎn)投資額。相關(guān)數(shù)據(jù)來(lái)源于統(tǒng)計(jì)局發(fā)行的 《中國(guó)統(tǒng)計(jì)年鑒——2015》。其中,31個(gè)省市的生產(chǎn)總值在《中國(guó)統(tǒng)計(jì)年鑒——2015》“國(guó)民經(jīng)濟(jì)核算”下的 “地區(qū)生產(chǎn)總值和指數(shù)”中,全社會(huì)固定資產(chǎn)投資在《中國(guó)固定資產(chǎn)投資統(tǒng)計(jì)年鑒——2015》“固定資產(chǎn)投資”下的“各地區(qū)全社會(huì)固定資產(chǎn)投資”中。

    3 地區(qū)生產(chǎn)總值與固定資產(chǎn)投資數(shù)據(jù)分析

    為了便于分析,把各地區(qū)的地區(qū)生產(chǎn)總值記作Y,各地區(qū)全社會(huì)固定資產(chǎn)投資記作X。為了初步分析地區(qū)生產(chǎn)總值(Y)與全社會(huì)固定資產(chǎn)投資額(X)的關(guān)系,可以作出以X為橫坐標(biāo),Y為縱坐標(biāo)的散點(diǎn)圖。通過(guò)對(duì)散點(diǎn)圖的觀察可以發(fā)現(xiàn)兩者間存在著某種線性關(guān)系,于是假設(shè)地區(qū)生產(chǎn)總值與全社會(huì)固定資產(chǎn)投資額有著密切聯(lián)系,下文進(jìn)一步探討與研究。

    3.1 圖形分析

    通過(guò)繪制地區(qū)生產(chǎn)總值和全社會(huì)固定資產(chǎn)投資額這兩個(gè)箱線圖,如圖1、圖2所示。對(duì)比可得知,地區(qū)生產(chǎn)總值和全社會(huì)固定資產(chǎn)投資額數(shù)據(jù)的中位數(shù)離下四分位數(shù)比離上四分位數(shù)近得多,數(shù)據(jù)均呈現(xiàn)出左偏分布,并且地區(qū)生產(chǎn)總值有三個(gè)明顯的離群值。三個(gè)離群值分別為廣東省,江蘇省和山東省的地區(qū)生產(chǎn)總值,筆者認(rèn)為,是由于上述三省為沿海省份,經(jīng)濟(jì)更為發(fā)達(dá)而導(dǎo)致的。

    圖1 各地區(qū)生產(chǎn)總值與全社會(huì)固定資產(chǎn)投資額的散點(diǎn)圖

    圖2 各地區(qū)生產(chǎn)總值與全社會(huì)固定資產(chǎn)投資額的箱線圖

    通過(guò)對(duì)各地區(qū)生產(chǎn)總值和全社會(huì)固定資產(chǎn)投資額莖葉圖的繪制,如圖3、圖4所示,可以直觀地看出,這兩個(gè)批的分布均是偏斜的,且各地生產(chǎn)總值還有三個(gè)邊遠(yuǎn)數(shù)據(jù)點(diǎn)。與從箱線圖中得出的結(jié)論一致。

    圖3 地區(qū)生產(chǎn)總值莖葉圖

    圖4 全社會(huì)固定資產(chǎn)投資莖葉圖

    3.2 線性分析——最小二乘法

    通過(guò)對(duì)地區(qū)生產(chǎn)總值與全社會(huì)固定資產(chǎn)投資額的簡(jiǎn)單線性分析可知,地區(qū)生產(chǎn)總值與全社會(huì)固定資產(chǎn)投資額之間的相關(guān)系數(shù)為0.777 854,調(diào)整后的相關(guān)系數(shù)為0.770 194,且通過(guò)了t檢驗(yàn)和F檢驗(yàn),符合假設(shè)檢驗(yàn)的條件。但由于調(diào)整后的相關(guān)系數(shù)只有0.770 194,故擬合效果不是很理想。

    表2 最小二乘法回歸結(jié)果

    根據(jù)相關(guān)數(shù)據(jù)進(jìn)行最小二乘回歸分析,可以得到的最小二乘法模型為:

    (1)

    (2726.172) (0.141 014)

    t=(-0.406 127) (10.076 95)

    R2=0.777 854F=101.545 0n=31

    圖5 最小二乘法下擬合直線

    由于原始數(shù)據(jù)具有一定的非對(duì)稱性,并且在地區(qū)生產(chǎn)總值數(shù)據(jù)右側(cè)的尾部有離群值存在,無(wú)法正確判斷兩者的數(shù)據(jù)分布情況,考慮到嚴(yán)謹(jǐn)性,我們嘗試進(jìn)行數(shù)據(jù)變換。

    4 數(shù)據(jù)變換

    霍格林指出:數(shù)據(jù)變換可以保持批中數(shù)據(jù)的次序;可以保持一個(gè)批的字母值,不計(jì)由于數(shù)據(jù)點(diǎn)之間作內(nèi)插造成的小差異;保證了原始批中非常靠近的點(diǎn),在重新表達(dá)后批中也非常靠近,至少相對(duì)于新尺度很靠近;保證函數(shù)沒有尖銳的拐角;并且由初等函數(shù)來(lái)規(guī)定,以便差不多用最簡(jiǎn)單的手持計(jì)算器就可以快而方便地做重新表達(dá)。因此數(shù)據(jù)變換可以達(dá)到令人滿意的效果,下面進(jìn)行一系列的數(shù)據(jù)變換,以期將不對(duì)稱的數(shù)據(jù)變?yōu)閷?duì)稱的數(shù)據(jù),并且消除異常值。

    4.1 以10為底的對(duì)數(shù)尺度數(shù)據(jù)變換

    將地區(qū)生產(chǎn)總值和全社會(huì)固定資產(chǎn)投資額進(jìn)行對(duì)數(shù)變換后,分別對(duì)數(shù)據(jù)做出箱線圖,如圖6所示。

    圖6 分別對(duì)地區(qū)生產(chǎn)總值和全社會(huì)固定 資產(chǎn)投資額進(jìn)行對(duì)數(shù)變換后的箱線圖

    通過(guò)箱線圖可以看出,經(jīng)過(guò)對(duì)數(shù)變換后,這兩個(gè)批的數(shù)據(jù)與變換前相比較為對(duì)稱,但是離群值變多。由于最小二乘回歸線不提供耐抗性,一個(gè)野的數(shù)據(jù)點(diǎn)就可以支配控制擬合線,若用這兩個(gè)新的數(shù)據(jù)批,估計(jì)出的參數(shù)就不穩(wěn)定,顯然,對(duì)數(shù)變換的效果不理想。

    4.2 平方根變換

    為方便表示,記A=Y^1/2,B=X^1/2,做出箱線圖,如圖7所示。

    圖7 分別對(duì)地區(qū)生產(chǎn)總值和全社會(huì)固定 資產(chǎn)投資額開平方根后的箱線圖

    通過(guò)箱線圖可以看出,經(jīng)過(guò)平方根變換后,這兩個(gè)批的數(shù)據(jù)均比較對(duì)稱,且對(duì)全社會(huì)固定資產(chǎn)投資額進(jìn)行平方根變換后,沒有離群值的存在,但地區(qū)生產(chǎn)總值經(jīng)平方根變換后存在三個(gè)異常值。因此,平方根變換的效果不夠理想。

    4.3 四次方根變換

    為方便表示,記M=Y^1/4,N=X^1/4。對(duì)數(shù)據(jù)做出箱線圖,如圖8所示。

    圖8 分別對(duì)地區(qū)生產(chǎn)總值和全社會(huì)固定 資產(chǎn)投資額開四次方根后的箱線圖

    觀察箱線圖,可以發(fā)現(xiàn)經(jīng)過(guò)四次方根變換后,這兩個(gè)批的數(shù)據(jù)均比較對(duì)稱,但地區(qū)生產(chǎn)總值經(jīng)四次方根變換后存在兩個(gè)異常值,全社會(huì)固定資產(chǎn)投資額經(jīng)四次方根變換后存在一個(gè)離群值。因此,經(jīng)過(guò)四次方根變換后,這兩個(gè)批的異常值并不能消失,由于野值使OLS法估計(jì)出的參數(shù)不穩(wěn)定,因此,四次方根變換同樣不能達(dá)到理想的效果。

    5 三組耐抗線分析

    經(jīng)過(guò)上述數(shù)據(jù)變換,這兩個(gè)批中的異常值并不能被消除,考慮到對(duì)經(jīng)濟(jì)意義的解釋,我們不再繼續(xù)進(jìn)行數(shù)據(jù)變換。由于最小二乘回歸線不提供耐抗性,一個(gè)野的數(shù)據(jù)點(diǎn)就可以支配控制擬合線,使它給出一個(gè)完全引入歧途的Y和X之間關(guān)系的概括公式,而三組耐抗線則避免了這個(gè)困難,因此我們采用三組耐抗線的方法對(duì)各地區(qū)生產(chǎn)總值(Y)和全社會(huì)固定資產(chǎn)投資額(X)的關(guān)系進(jìn)行探索。

    對(duì)于數(shù)據(jù)(x,y),首先把x的值按升序排序使得x1≤x2≤…≤xn,然后在這些排過(guò)序的值的基礎(chǔ)上,將(xi,yi)分割成三個(gè)組(左組,中組,右組),使組的大小盡可能的相等,見表3。

    表3 分組后數(shù)據(jù)

    在三個(gè)組中,我們分別求出x和y的中位數(shù),得到三個(gè)組的總括點(diǎn)(XL,XL),(XM,XM), 以及(XR,YR)。組內(nèi)總括點(diǎn)使得擬合線具有耐抗性,為了更好地?cái)M合數(shù)據(jù),我們以斜率b0和中心值XM來(lái)進(jìn)行擬合。求得的(XL,YL),(XM,YM),(XR,YR),分別為(6 470.33,8 051.605),(13 843.22,15 714.63),(25 512.365,32 179.695)。因此初始斜率為:

    (2)

    (3)

    數(shù)據(jù)點(diǎn)分成的三個(gè)組以及初始線的殘差,見表4。

    表4 地區(qū)生產(chǎn)總值和全社會(huì)固定資產(chǎn)投資額——三個(gè)組和初始?xì)埐?/p>

    (4)

    (5)

    因此,一次調(diào)整后的斜率為:

    (6)

    (7)

    可以看到,δ1按絕對(duì)值大小比b0小得多,但超過(guò)了b0絕對(duì)值的1 %,還不是可以忽略不計(jì)的。于是,做第二步迭代。分別求出X和R(1)的中位數(shù),得到三個(gè)組的總括點(diǎn):(6 470.33,133.475 271 3),(13 843.22,-534.939 879 3),(25 512.365,178.002 558 6),通過(guò)計(jì)算可以得出:

    (8)

    (9)

    因此,二次迭代后斜率為:

    (10)

    (11)

    可以看到,δ2的絕對(duì)值比b0的絕對(duì)值小得多,且小于b0絕對(duì)值的1 %(0.012 671),于是停止迭代。因此,結(jié)果的擬合線為:

    (12)

    圖9 三組耐抗線下擬合直線

    離這條線的殘差對(duì)全社會(huì)固定資產(chǎn)投資額(X)的散點(diǎn)圖如圖10所示。

    圖10 擬合耐抗線后地區(qū)生產(chǎn)總值殘差 對(duì)全社會(huì)固定資產(chǎn)投資額散點(diǎn)圖

    總體來(lái)看,除去在外值之后,散點(diǎn)圖較為平坦??梢钥闯?,稍微異常的數(shù)據(jù)點(diǎn)幾乎對(duì)于概括了數(shù)據(jù)主體部分的直線沒有影響。

    圖11 最小二乘回歸殘差對(duì)全社會(huì)固定資產(chǎn)投資額散點(diǎn)圖

    上圖為最小二乘殘差對(duì)全社會(huì)固定資產(chǎn)總額的散點(diǎn)圖,可以看出它與圖10很相似,但是這個(gè)圖提示一個(gè)略微向下的趨勢(shì)。這就是說(shuō),如果去掉一個(gè)稍微負(fù)斜率的直線后,最小二乘殘差看起來(lái)會(huì)更水平一些。

    6 最小二乘分析模型優(yōu)化

    6.1 異常點(diǎn)分析

    6.1.1 殘差圖

    (13)

    將各省的全社會(huì)固定資產(chǎn)投資(X)與其殘差(ri)畫圖分析,見圖12。

    圖12 最小二乘回歸原始?xì)埐顖D

    由圖12可看到有一個(gè)點(diǎn)明顯超出10 000的范圍,而另兩個(gè)點(diǎn)也略超出10 000,則可判定有異常值的存在。嚴(yán)重超出1 000的樣本為廣東,略超出10 000的有上海、北京。

    6.1.2 杠桿率圖

    (14)

    圖13 杠桿率圖

    在本實(shí)例中,由于此模型為一元線性回歸模型,只有截距與斜率兩個(gè)參數(shù),又因?yàn)橐还灿?1個(gè)樣本數(shù)據(jù),因此這31個(gè)殘差的平均杠桿率為2 /31,即0. 064 5。而從圖中可以看出有兩個(gè)殘差的杠桿率遠(yuǎn)遠(yuǎn)大于其平均杠桿率,即山東和江蘇,有一個(gè)殘差的杠桿率略微大于其平均杠桿率,為西藏。

    6.1.3 標(biāo)準(zhǔn)化殘差圖

    (15)

    (16)

    (17)

    (18)

    從圖14可以看出,由于標(biāo)準(zhǔn)化后的殘差消除了異方差,且服從正態(tài)分布,則其標(biāo)準(zhǔn)化殘差的絕對(duì)值超過(guò)Za/2(1.96)的概率為0.05,而由圖看出有一個(gè)值明顯超出了1.96,還有一個(gè)值略微超過(guò)1.96,可找到超過(guò)1.96的異常值分別為廣東、上海。

    圖14 標(biāo)準(zhǔn)化殘差圖

    6.1.4.殘差箱線圖

    從上四分?jǐn)?shù)Fu向上跨出dF的3/2倍,從下四分?jǐn)?shù)向下跨出dF的3/2倍,得到兩個(gè)截?cái)帱c(diǎn),把這兩個(gè)截?cái)帱c(diǎn)之外的觀測(cè)值當(dāng)作離群值。根據(jù)這一定義,直接畫出殘差的箱線圖即可一目了然地發(fā)現(xiàn)離群值,即廣東、上海、北京、浙江四地,如圖15所示。

    圖15 殘差箱線圖

    通過(guò)以上對(duì)于異常點(diǎn)的不同分析,可看出廣東、上海在4種方式的異常點(diǎn)判斷中均出現(xiàn)了三次,但廣東每次都是遠(yuǎn)遠(yuǎn)超出標(biāo)準(zhǔn),而上海有兩次都是超出標(biāo)準(zhǔn)不多;北京出現(xiàn)了兩次,且每次都是稍微超出標(biāo)準(zhǔn)而被選為異常點(diǎn);山東、江蘇、浙江、西藏僅僅出現(xiàn)了一次。綜合以上分析,筆者認(rèn)為應(yīng)將廣東省剔除,以避免異常值對(duì)于回歸擬合的影響。

    6.2 模型優(yōu)化

    將廣東省對(duì)應(yīng)的數(shù)據(jù)剔除之后,再次進(jìn)行最小二乘回歸,回歸結(jié)果如下:

    表5 改進(jìn)后的最小二乘法回歸結(jié)果

    (19)

    (185 9.479) (0.097 576)

    t=(-0.295 814) (13.531 30)

    R2=0.867 359 F=183.096 1 n=30

    圖16 改進(jìn)后的散點(diǎn)圖和擬合直線

    6.3 模型優(yōu)化效果評(píng)價(jià)

    用均方誤差來(lái)比較直接用最小二乘法擬合的直線與剔除異常點(diǎn)后用最小二乘法擬合的直線效果的優(yōu)劣。

    通過(guò)計(jì)算得出最初的最小二乘回歸的均方誤差為6.193 5×107,優(yōu)化后的均方誤差為2.873 3×107,由此明顯看出去掉一個(gè)離群點(diǎn)后的擬合偏誤要小得多,模型的優(yōu)化是有效的。

    7 三組耐抗線與最小二乘法的比較

    由三組耐抗線算得的均方誤差為7.332 9×107,大于用最初的最小二乘法算出的均方誤差6.193 5×107。

    圖17 原始數(shù)據(jù)下三組耐抗線與最小二乘法殘差的比較

    將原始樣本數(shù)據(jù)的最小二乘法計(jì)算于三組耐抗線計(jì)算所得到的殘差畫在一張圖中,得到圖17。其中藍(lán)色的圈代表運(yùn)用最小二乘法擬合得到的殘差,綠色的圈表示運(yùn)用三組耐抗線進(jìn)行擬合后得到的殘差。

    通過(guò)觀察殘差對(duì)X(全社會(huì)固定資產(chǎn)投資額)的散點(diǎn)圖可以發(fā)現(xiàn),運(yùn)用三組耐抗線方法得到的殘差中有距離均值非常遠(yuǎn)的點(diǎn),而最小二乘法得到的殘差到中心的距離則沒有三組耐抗線的殘差距中心遠(yuǎn)。由此可以看出,三組耐抗線本身具有耐抗性,它不因個(gè)別偏離較大的點(diǎn)而劇烈影響所擬合的參數(shù)值。因此即便有差異很大的點(diǎn),它也會(huì)將離群點(diǎn)對(duì)擬合參數(shù)的影響降到很小,因此在殘差圖上就可能出現(xiàn)某些離群點(diǎn)距離擬合直線上對(duì)應(yīng)點(diǎn)的距離很大的情況。而最小二乘法則不然,因?yàn)樗沟脷埐钇椒胶妥钚?,因而考慮到了包括異常值在內(nèi)的所有的數(shù)。對(duì)于越偏遠(yuǎn)的點(diǎn),最小二乘法的擬合效果越差,因?yàn)樗鼮榱耸箽埐畹钠椒胶妥钚。赡軙?huì)舍去其他擬合效果較好的模型。因此對(duì)于異常值較多或者偏離較遠(yuǎn)的情況下,三組耐抗線的方法相對(duì)于最小二乘法更加可靠。但在上述數(shù)據(jù)中,由于僅有一個(gè)野值,并且其離群的程度不是很大,因此由三組耐抗線算得的均方誤差要大于最小二乘法算出的均方誤差。

    8 結(jié)論

    [1] (美)(D.C.霍格林)David C.Hoaglin. 探索性數(shù)據(jù)分析[M].陳忠璉,郭德媛,譯.北京:中國(guó)統(tǒng)計(jì)出版社, 1998.

    [2] 江克寧.地區(qū)生產(chǎn)總值與固定資產(chǎn)投資關(guān)系的實(shí)證分析[J].統(tǒng)計(jì)與決策,2015(17):97-100.

    [3] 吳建群.探索性數(shù)據(jù)分析在能源消費(fèi)中的應(yīng)用探討——人均生活能源消費(fèi)和人均GDP的關(guān)系研究[J].現(xiàn)代商貿(mào)工業(yè),2011(17):46-47.

    2017-03-14

    呂林靜(1997-),女,海南??谌耍芯糠较?yàn)榻?jīng)濟(jì)學(xué),電話:17889982336。

    F127;P224

    A

    1671-4733(2017)04-0118-05

    猜你喜歡
    投資額離群線圖
    個(gè)體化預(yù)測(cè)結(jié)腸癌術(shù)后發(fā)生并發(fā)癥風(fēng)險(xiǎn)列線圖模型的建立
    基于箱線圖的出廠水和管網(wǎng)水水質(zhì)分析
    Китай вышел на второе место в мире по объемам прямых внешних инвестиций
    中亞信息(2016年9期)2017-01-18 05:50:04
    東山頭遺址采集石器線圖
    離群數(shù)據(jù)挖掘在發(fā)現(xiàn)房產(chǎn)銷售潛在客戶中的應(yīng)用
    離群的小雞
    2015 年1~9月機(jī)械工業(yè)固定資產(chǎn)投資增長(zhǎng) 8.7%
    2015年1~6月機(jī)械工業(yè)固定資產(chǎn)投資增長(zhǎng)9.55%
    應(yīng)用相似度測(cè)量的圖離群點(diǎn)檢測(cè)方法
    一種基于核空間局部離群因子的離群點(diǎn)挖掘方法
    日本五十路高清| 免费看十八禁软件| 成年免费大片在线观看| 亚洲av中文字字幕乱码综合| 精品国产亚洲在线| 久久中文看片网| 日本a在线网址| 操出白浆在线播放| 国产激情偷乱视频一区二区| 日韩欧美 国产精品| 丰满的人妻完整版| 国产精品爽爽va在线观看网站| www.www免费av| 色av中文字幕| 免费在线观看视频国产中文字幕亚洲| 男人舔奶头视频| www.999成人在线观看| 免费看十八禁软件| 国产高清三级在线| 欧美三级亚洲精品| 国产高清视频在线观看网站| 久久九九热精品免费| 日韩欧美在线二视频| 又黄又粗又硬又大视频| 国产精品久久久av美女十八| 亚洲九九香蕉| 亚洲精品456在线播放app | 90打野战视频偷拍视频| 在线免费观看不下载黄p国产 | 久久久精品欧美日韩精品| 午夜福利视频1000在线观看| 亚洲欧美日韩高清在线视频| 999久久久精品免费观看国产| 在线永久观看黄色视频| 十八禁人妻一区二区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕日韩| 五月玫瑰六月丁香| 精品熟女少妇八av免费久了| 色综合婷婷激情| 国产精品综合久久久久久久免费| 欧美日韩中文字幕国产精品一区二区三区| 午夜精品久久久久久毛片777| 黑人欧美特级aaaaaa片| av在线蜜桃| 亚洲av美国av| 日韩有码中文字幕| 一级毛片高清免费大全| 日本五十路高清| 亚洲狠狠婷婷综合久久图片| www国产在线视频色| 中文资源天堂在线| 高清在线国产一区| 国产精品久久久久久久电影 | 亚洲精品中文字幕一二三四区| 色视频www国产| 在线国产一区二区在线| 又粗又爽又猛毛片免费看| 99精品欧美一区二区三区四区| 十八禁人妻一区二区| svipshipincom国产片| 成人国产综合亚洲| 久久精品国产清高在天天线| 一本一本综合久久| 国产激情偷乱视频一区二区| 久久国产乱子伦精品免费另类| 亚洲精品粉嫩美女一区| 一个人免费在线观看电影 | 日韩人妻高清精品专区| 久久久久久大精品| 中出人妻视频一区二区| 变态另类成人亚洲欧美熟女| 国产免费av片在线观看野外av| 欧美激情在线99| 午夜福利成人在线免费观看| 在线a可以看的网站| 黄色女人牲交| 日本 欧美在线| 国产高清videossex| 在线a可以看的网站| 免费av毛片视频| 国产精品久久视频播放| 精品国内亚洲2022精品成人| 麻豆久久精品国产亚洲av| 亚洲熟妇中文字幕五十中出| 国产成人av教育| 少妇的逼水好多| 97碰自拍视频| 免费大片18禁| 老司机福利观看| 成人国产综合亚洲| 国产免费av片在线观看野外av| 99国产精品一区二区蜜桃av| 熟女少妇亚洲综合色aaa.| 成人国产综合亚洲| 亚洲国产欧洲综合997久久,| 国产日本99.免费观看| 亚洲av成人一区二区三| 国产一区二区在线观看日韩 | 精品一区二区三区视频在线 | av视频在线观看入口| 欧美极品一区二区三区四区| 中文亚洲av片在线观看爽| 免费av毛片视频| 久久性视频一级片| 9191精品国产免费久久| 99re在线观看精品视频| 青草久久国产| 91在线观看av| 久久亚洲精品不卡| 午夜精品在线福利| 欧美乱码精品一区二区三区| 欧美日本亚洲视频在线播放| 精品免费久久久久久久清纯| 真实男女啪啪啪动态图| 男人舔奶头视频| 免费大片18禁| 国产又黄又爽又无遮挡在线| 国产三级中文精品| 99国产精品一区二区蜜桃av| 日韩大尺度精品在线看网址| 美女大奶头视频| 亚洲欧美日韩卡通动漫| 久久久久精品国产欧美久久久| 欧美三级亚洲精品| 国产精品精品国产色婷婷| 又爽又黄无遮挡网站| 少妇熟女aⅴ在线视频| 亚洲片人在线观看| 99在线视频只有这里精品首页| 中文字幕av在线有码专区| 男女午夜视频在线观看| 久久亚洲真实| 日日夜夜操网爽| 人人妻,人人澡人人爽秒播| 国内毛片毛片毛片毛片毛片| 亚洲天堂国产精品一区在线| 嫩草影视91久久| 免费观看人在逋| 欧美另类亚洲清纯唯美| 日本黄色视频三级网站网址| 舔av片在线| 老司机福利观看| 国产视频一区二区在线看| 亚洲18禁久久av| 国产麻豆成人av免费视频| 国产亚洲精品一区二区www| 99re在线观看精品视频| 别揉我奶头~嗯~啊~动态视频| 中文字幕高清在线视频| 国产高清激情床上av| 免费观看的影片在线观看| 国产av麻豆久久久久久久| 久久人人精品亚洲av| 夜夜爽天天搞| 亚洲av中文字字幕乱码综合| 国产蜜桃级精品一区二区三区| 免费观看的影片在线观看| 99国产精品99久久久久| 男女做爰动态图高潮gif福利片| 国产69精品久久久久777片 | 成年人黄色毛片网站| 国产激情欧美一区二区| 99re在线观看精品视频| 操出白浆在线播放| 在线观看66精品国产| 欧美黑人欧美精品刺激| 国产v大片淫在线免费观看| 在线观看免费视频日本深夜| 国产一区在线观看成人免费| 欧美黑人巨大hd| 国产aⅴ精品一区二区三区波| 日韩国内少妇激情av| 少妇人妻一区二区三区视频| 嫁个100分男人电影在线观看| 一级毛片精品| 九九久久精品国产亚洲av麻豆 | 一个人看的www免费观看视频| 88av欧美| 久久久国产成人精品二区| 一个人看的www免费观看视频| 亚洲国产欧洲综合997久久,| 欧美色欧美亚洲另类二区| 国产探花在线观看一区二区| 操出白浆在线播放| 99在线视频只有这里精品首页| 美女扒开内裤让男人捅视频| 欧美精品啪啪一区二区三区| 日韩精品青青久久久久久| 巨乳人妻的诱惑在线观看| 国产午夜精品久久久久久| а√天堂www在线а√下载| 午夜福利视频1000在线观看| 99久久综合精品五月天人人| 国产激情欧美一区二区| 婷婷精品国产亚洲av| 波多野结衣高清无吗| 免费电影在线观看免费观看| 少妇裸体淫交视频免费看高清| 99在线视频只有这里精品首页| 欧美国产日韩亚洲一区| 日本熟妇午夜| 日韩人妻高清精品专区| 亚洲av成人av| 1024手机看黄色片| 精品一区二区三区av网在线观看| 国产亚洲精品久久久久久毛片| 久久精品国产清高在天天线| 俺也久久电影网| 高清在线国产一区| 又紧又爽又黄一区二区| 国产蜜桃级精品一区二区三区| 国产黄色小视频在线观看| x7x7x7水蜜桃| 免费在线观看亚洲国产| 香蕉国产在线看| 天堂影院成人在线观看| www日本在线高清视频| 18禁黄网站禁片免费观看直播| 国产精品自产拍在线观看55亚洲| 国产69精品久久久久777片 | 亚洲av美国av| 熟女人妻精品中文字幕| 中文字幕av在线有码专区| 色哟哟哟哟哟哟| 国产黄a三级三级三级人| 久久九九热精品免费| 精品久久久久久久人妻蜜臀av| 悠悠久久av| 琪琪午夜伦伦电影理论片6080| 三级毛片av免费| 男人舔奶头视频| 久久亚洲精品不卡| 黄色日韩在线| 好男人电影高清在线观看| 欧美性猛交黑人性爽| 色视频www国产| 免费无遮挡裸体视频| 无遮挡黄片免费观看| cao死你这个sao货| 蜜桃久久精品国产亚洲av| 亚洲人成网站在线播放欧美日韩| 宅男免费午夜| 三级毛片av免费| 成人永久免费在线观看视频| 精华霜和精华液先用哪个| 免费搜索国产男女视频| 欧美国产日韩亚洲一区| 亚洲熟女毛片儿| 黄色女人牲交| 欧美zozozo另类| 亚洲一区高清亚洲精品| 99国产精品99久久久久| 大型黄色视频在线免费观看| 亚洲av免费在线观看| 国产成人影院久久av| 久久久久精品国产欧美久久久| 日韩精品中文字幕看吧| 国产一区在线观看成人免费| 日韩欧美免费精品| 老司机午夜十八禁免费视频| 一级毛片高清免费大全| 好男人电影高清在线观看| 欧美日韩中文字幕国产精品一区二区三区| 在线观看日韩欧美| 波多野结衣巨乳人妻| 欧美av亚洲av综合av国产av| 99国产精品一区二区蜜桃av| 无限看片的www在线观看| 国产欧美日韩精品亚洲av| 亚洲成人久久性| 午夜亚洲福利在线播放| 看免费av毛片| 国产精品久久久久久久电影 | 香蕉av资源在线| 亚洲成人久久爱视频| 在线观看66精品国产| 最近视频中文字幕2019在线8| 午夜久久久久精精品| 在线视频色国产色| 一进一出抽搐动态| 不卡一级毛片| 国产精品 国内视频| 亚洲一区二区三区不卡视频| 90打野战视频偷拍视频| 桃色一区二区三区在线观看| 12—13女人毛片做爰片一| 欧美黑人巨大hd| 无限看片的www在线观看| 国产一区二区激情短视频| 成熟少妇高潮喷水视频| 精品人妻1区二区| 国产亚洲精品久久久com| www日本在线高清视频| 精品久久久久久久毛片微露脸| 国产在线精品亚洲第一网站| 在线永久观看黄色视频| 日本在线视频免费播放| 午夜a级毛片| 成人鲁丝片一二三区免费| 亚洲专区中文字幕在线| 免费观看的影片在线观看| 欧美日韩精品网址| 51午夜福利影视在线观看| 真人做人爱边吃奶动态| 变态另类丝袜制服| 国内精品久久久久精免费| 免费av毛片视频| 午夜福利成人在线免费观看| 国产一区在线观看成人免费| 老司机午夜福利在线观看视频| 国产三级中文精品| 国产成人aa在线观看| 黑人巨大精品欧美一区二区mp4| 两个人视频免费观看高清| 91九色精品人成在线观看| 国产黄a三级三级三级人| xxx96com| 日本 av在线| 欧美中文日本在线观看视频| 人人妻人人澡欧美一区二区| 99久久精品国产亚洲精品| 18禁观看日本| 99久久国产精品久久久| www.自偷自拍.com| 国产成人一区二区三区免费视频网站| 999久久久国产精品视频| 国产精华一区二区三区| 国产午夜福利久久久久久| 在线观看午夜福利视频| 国产视频一区二区在线看| 啪啪无遮挡十八禁网站| 婷婷精品国产亚洲av| 两个人视频免费观看高清| 欧美乱色亚洲激情| 在线观看美女被高潮喷水网站 | 99精品久久久久人妻精品| 成年女人看的毛片在线观看| 久久午夜亚洲精品久久| 午夜福利在线观看免费完整高清在 | 国产乱人伦免费视频| 在线观看美女被高潮喷水网站 | av片东京热男人的天堂| 婷婷精品国产亚洲av在线| 一级作爱视频免费观看| 亚洲中文字幕一区二区三区有码在线看 | 中文字幕av在线有码专区| 床上黄色一级片| 国产综合懂色| 国产成年人精品一区二区| 最近最新中文字幕大全电影3| 精品国产三级普通话版| 一卡2卡三卡四卡精品乱码亚洲| 亚洲美女黄片视频| 熟女少妇亚洲综合色aaa.| 亚洲精品在线观看二区| 国产日本99.免费观看| 日韩精品青青久久久久久| 可以在线观看的亚洲视频| 黄色 视频免费看| 欧美乱妇无乱码| 久久精品91蜜桃| 一本久久中文字幕| 精品国产亚洲在线| 久久久国产成人免费| 国产99白浆流出| 两性夫妻黄色片| 亚洲成人免费电影在线观看| 日本一二三区视频观看| 老汉色∧v一级毛片| tocl精华| 国产精品综合久久久久久久免费| 黄色片一级片一级黄色片| 日韩人妻高清精品专区| 免费在线观看成人毛片| 两性夫妻黄色片| 黑人操中国人逼视频| 丰满人妻一区二区三区视频av | 精品熟女少妇八av免费久了| 女同久久另类99精品国产91| a级毛片a级免费在线| 久久国产精品人妻蜜桃| 精华霜和精华液先用哪个| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久久电影 | tocl精华| 老司机深夜福利视频在线观看| 成人特级av手机在线观看| 99久久成人亚洲精品观看| 精品久久久久久,| 国产淫片久久久久久久久 | 哪里可以看免费的av片| 又爽又黄无遮挡网站| 日本黄大片高清| 99国产精品一区二区蜜桃av| 丁香六月欧美| 久久99热这里只有精品18| xxxwww97欧美| 999久久久精品免费观看国产| 精品免费久久久久久久清纯| 国产探花在线观看一区二区| 久久亚洲精品不卡| aaaaa片日本免费| 999久久久精品免费观看国产| 两个人视频免费观看高清| 成人欧美大片| 国产蜜桃级精品一区二区三区| 99热6这里只有精品| 国产视频一区二区在线看| 狠狠狠狠99中文字幕| 精品久久久久久久人妻蜜臀av| 亚洲激情在线av| 欧美一级毛片孕妇| www.自偷自拍.com| 两人在一起打扑克的视频| 国产野战对白在线观看| 欧美日韩亚洲国产一区二区在线观看| 九九久久精品国产亚洲av麻豆 | 亚洲熟妇中文字幕五十中出| 国产精品一区二区免费欧美| www.999成人在线观看| 成人欧美大片| 国内毛片毛片毛片毛片毛片| 亚洲av五月六月丁香网| 国产精品99久久99久久久不卡| 99久久无色码亚洲精品果冻| 成人特级黄色片久久久久久久| 后天国语完整版免费观看| 精品电影一区二区在线| av片东京热男人的天堂| 色噜噜av男人的天堂激情| 亚洲第一电影网av| 一二三四社区在线视频社区8| 首页视频小说图片口味搜索| 少妇丰满av| 日韩欧美国产在线观看| 九色成人免费人妻av| 日韩欧美在线乱码| 成人午夜高清在线视频| 在线观看免费视频日本深夜| 精品熟女少妇八av免费久了| 久久人人精品亚洲av| 人妻久久中文字幕网| 国产精品久久久久久精品电影| 精品久久久久久久久久免费视频| 国产一区二区激情短视频| 中文资源天堂在线| 夜夜爽天天搞| 老汉色∧v一级毛片| 国产成人一区二区三区免费视频网站| 成人三级黄色视频| 级片在线观看| 国产精品九九99| 国产午夜精品论理片| 亚洲aⅴ乱码一区二区在线播放| 一二三四在线观看免费中文在| 可以在线观看毛片的网站| 午夜福利欧美成人| 九色国产91popny在线| 麻豆国产97在线/欧美| 韩国av一区二区三区四区| 日日干狠狠操夜夜爽| 一二三四在线观看免费中文在| 香蕉国产在线看| 日韩精品青青久久久久久| 免费看日本二区| 久久婷婷人人爽人人干人人爱| 91老司机精品| 一个人免费在线观看的高清视频| 国产毛片a区久久久久| 中国美女看黄片| 色噜噜av男人的天堂激情| 中文在线观看免费www的网站| 老汉色av国产亚洲站长工具| 欧美乱色亚洲激情| 一边摸一边抽搐一进一小说| 国产精品综合久久久久久久免费| www国产在线视频色| 日韩欧美一区二区三区在线观看| 窝窝影院91人妻| 欧美日韩乱码在线| 亚洲人成网站在线播放欧美日韩| 99久久国产精品久久久| 国产又色又爽无遮挡免费看| 久久热在线av| a级毛片a级免费在线| 久久久久国产精品人妻aⅴ院| 母亲3免费完整高清在线观看| 身体一侧抽搐| www国产在线视频色| 久久天堂一区二区三区四区| 少妇的逼水好多| 欧美zozozo另类| 国产午夜福利久久久久久| x7x7x7水蜜桃| 99久久精品热视频| 国产激情偷乱视频一区二区| 最近最新中文字幕大全免费视频| 动漫黄色视频在线观看| 精品一区二区三区四区五区乱码| 国产高清激情床上av| 人妻丰满熟妇av一区二区三区| 黄色视频,在线免费观看| 国产精品98久久久久久宅男小说| 成人av在线播放网站| 国产成人系列免费观看| 国产又黄又爽又无遮挡在线| 天天添夜夜摸| 亚洲精品一卡2卡三卡4卡5卡| 成人av在线播放网站| 97人妻精品一区二区三区麻豆| 亚洲欧洲精品一区二区精品久久久| 国产午夜精品久久久久久| 日本熟妇午夜| 亚洲色图 男人天堂 中文字幕| 亚洲美女黄片视频| 国产 一区 欧美 日韩| 成人18禁在线播放| 国产三级黄色录像| 岛国在线观看网站| 午夜精品在线福利| 91av网一区二区| 精品一区二区三区av网在线观看| 在线观看美女被高潮喷水网站 | 制服人妻中文乱码| 久久久国产精品麻豆| 色老头精品视频在线观看| 国产精华一区二区三区| 午夜福利18| 日韩欧美国产在线观看| 国产精品自产拍在线观看55亚洲| 精品不卡国产一区二区三区| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩一区二区三| 日韩欧美在线乱码| 这个男人来自地球电影免费观看| 精品人妻1区二区| 亚洲欧洲精品一区二区精品久久久| a级毛片a级免费在线| 亚洲精品一卡2卡三卡4卡5卡| 一进一出抽搐动态| 亚洲美女视频黄频| 人人妻人人看人人澡| 午夜福利在线观看吧| 日韩欧美国产一区二区入口| 亚洲欧洲精品一区二区精品久久久| 噜噜噜噜噜久久久久久91| 小说图片视频综合网站| 97超级碰碰碰精品色视频在线观看| 国产欧美日韩一区二区三| netflix在线观看网站| 午夜福利在线观看吧| 国产精品久久久久久亚洲av鲁大| 亚洲欧洲精品一区二区精品久久久| 国产一区二区三区在线臀色熟女| 99视频精品全部免费 在线 | 香蕉av资源在线| 91在线观看av| 午夜福利视频1000在线观看| 国产成人系列免费观看| 国产探花在线观看一区二区| 免费搜索国产男女视频| 男人舔奶头视频| 中文字幕最新亚洲高清| 我要搜黄色片| 我的老师免费观看完整版| 欧美极品一区二区三区四区| 麻豆国产97在线/欧美| 国产精品98久久久久久宅男小说| 国内精品久久久久久久电影| 久久天堂一区二区三区四区| 国产伦在线观看视频一区| 久久午夜亚洲精品久久| 一a级毛片在线观看| 成人精品一区二区免费| 久久久久国内视频| 无遮挡黄片免费观看| 日韩大尺度精品在线看网址| 亚洲黑人精品在线| 久久精品人妻少妇| 中文字幕最新亚洲高清| 一个人免费在线观看电影 | 精品久久蜜臀av无| 国产午夜精品久久久久久| 国产精品1区2区在线观看.| 很黄的视频免费| 欧美黄色片欧美黄色片| 不卡av一区二区三区| 可以在线观看的亚洲视频| 一进一出好大好爽视频| 99久久综合精品五月天人人| 女生性感内裤真人,穿戴方法视频| 亚洲专区国产一区二区| 色视频www国产| 久久精品影院6| 国产精品98久久久久久宅男小说| 欧美色欧美亚洲另类二区| 亚洲18禁久久av| 老司机在亚洲福利影院| 两人在一起打扑克的视频| 亚洲av成人av| 免费av不卡在线播放| 首页视频小说图片口味搜索| 国产成人福利小说| 在线观看免费视频日本深夜| 亚洲精品456在线播放app | 又黄又爽又免费观看的视频| 午夜激情欧美在线| 宅男免费午夜| 国内精品一区二区在线观看| 久久这里只有精品19| 欧美黑人欧美精品刺激| 老司机在亚洲福利影院| 色视频www国产| 国产一区二区在线av高清观看|