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

    禽流感疫情變動(dòng)對(duì)畜禽產(chǎn)品價(jià)格的動(dòng)態(tài)影響研究
    ——基于時(shí)變參數(shù)向量自回歸(TVP-VAR)模型

    2018-10-11 06:36:32鄭燕馬驥
    關(guān)鍵詞:疫情模型

    鄭燕,馬驥

    (中國農(nóng)業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院,北京 100083)

    近年來,禽流感疫情對(duì)我國畜禽市場(chǎng)沖擊較大,其已與金融危機(jī)、恐怖主義等不可預(yù)測(cè)事件,統(tǒng)稱為現(xiàn)今最具挑戰(zhàn)的市場(chǎng)風(fēng)險(xiǎn)之一[1]。已有研究表明,疫病是畜禽產(chǎn)品價(jià)格波動(dòng)的主要原因之一[2-3]。近幾年,有關(guān)禽流感疫情的影響逐漸受到國內(nèi)外學(xué)者的關(guān)注,相關(guān)研究主要集中在禽流感疫情對(duì)價(jià)格的影響[4-5]、禽流感疫情對(duì)經(jīng)濟(jì)的影響[6-7]以及禽流感防治措施及補(bǔ)償問題[8-10]等方面。在禽流感疫情數(shù)據(jù)的量化方面,學(xué)者們并沒有統(tǒng)一的指標(biāo)。當(dāng)禽流感造成大量家禽死亡時(shí),主要衡量指標(biāo)為代表禽流感疫情是否發(fā)生的虛擬變量[4]、禽流感發(fā)生次數(shù)[6]、畜禽撲殺量[11]等,但這些指標(biāo)并不能很好地反映禽流感發(fā)生的影響強(qiáng)度,而且數(shù)據(jù)的有效性和及時(shí)性較差。針對(duì)這一方面,國外學(xué)者采用指數(shù)構(gòu)建的方式進(jìn)行研究,如Roux[12]通過建立食品恐慌指數(shù)評(píng)估了禽流感疫情對(duì)埃及家禽行業(yè)垂直價(jià)格傳導(dǎo)的影響;Serra[13]通過建立信息索引變量評(píng)估了禽流感疫情對(duì)土耳其家禽產(chǎn)業(yè)市場(chǎng)價(jià)格傳遞的影響。但是,當(dāng)前我國并沒有建立這方面的指數(shù),這為研究禽流感疫情與畜禽產(chǎn)品價(jià)格波動(dòng)的關(guān)系方面帶來了不便。但值得慶幸的是,互聯(lián)網(wǎng)大數(shù)據(jù)的快速發(fā)展為相關(guān)研究帶來了新的思路,尤其是網(wǎng)絡(luò)搜索引擎的發(fā)展為經(jīng)濟(jì)學(xué)、社會(huì)學(xué)的研究開辟了新的領(lǐng)域。由于信息搜索是公眾為滿足一定的目標(biāo)需求而進(jìn)行的信息查尋活動(dòng),而用戶需求又受到外部環(huán)境的影響,因此,用戶關(guān)鍵詞的搜索與外部環(huán)境必然存在對(duì)應(yīng)的關(guān)系[14],而且根據(jù)搜索引擎提供的詞頻查詢工具獲得的關(guān)鍵詞搜索趨勢(shì)與流行趨勢(shì)保持了高度的一致性[15],能夠準(zhǔn)確地反映公眾對(duì)于某一事件的關(guān)注程度變化[16]。除此之外,由于禽流感關(guān)注度和畜禽產(chǎn)品市場(chǎng)價(jià)格波動(dòng)起伏,可能會(huì)出現(xiàn)結(jié)構(gòu)突變的現(xiàn)象,在不同時(shí)期,禽流感對(duì)畜禽產(chǎn)品價(jià)格的影響亦有所不同。以往較多研究利用線性或固定系數(shù)模型進(jìn)行分析,但利用傳統(tǒng)的線性或固定系數(shù)模型,不僅容易造成系統(tǒng)性偏誤,而且會(huì)影響估計(jì)效果,不利于準(zhǔn)確把握隨時(shí)間變化禽流感疫情對(duì)畜禽產(chǎn)品價(jià)格的影響。

    綜上所述,現(xiàn)有研究存在禽流感疫情度量以及研究方法的局限性,并不能在特定的社會(huì)背景下,分時(shí)期、分時(shí)點(diǎn)、分品種地研究禽流感疫情變動(dòng)對(duì)畜禽產(chǎn)品價(jià)格的沖擊作用,而這一點(diǎn)對(duì)我國政府制定相關(guān)政策、穩(wěn)定畜禽產(chǎn)品價(jià)格具有至關(guān)重要的戰(zhàn)略意義?;诖耍疚囊环矫娼梃b國外學(xué)者相關(guān)研究思路,借助全球最大的中文搜索引擎——百度,利用其提供的百度指數(shù)功能,通過構(gòu)建禽流感輿情指數(shù)作為禽流感疫情變動(dòng)的代理變量,以此彌補(bǔ)禽流感疫情度量指標(biāo)的不足;另一方面借助時(shí)變參數(shù)向量自回歸模型(TVP-VAR),分時(shí)期地研究禽流感疫情變動(dòng)對(duì)畜禽產(chǎn)品價(jià)格的動(dòng)態(tài)沖擊作用,深入剖析禽流感疫情變動(dòng)與畜禽價(jià)格波動(dòng)之間的關(guān)系,探究禽流感疫情變動(dòng)在社會(huì)輿論的“風(fēng)險(xiǎn)放大效應(yīng)”下對(duì)畜禽產(chǎn)品價(jià)格波動(dòng)的動(dòng)態(tài)影響。

    1 影響機(jī)制

    根據(jù)供求理論,畜禽產(chǎn)品的供求關(guān)系決定了畜禽產(chǎn)品的價(jià)格,以雞肉及其替代品豬肉為例,結(jié)合圖1來看,雞肉供給曲線為S0,雞肉需求曲線為D0,兩條曲線相交,形成雞肉的均衡價(jià)格P0。當(dāng)其他條件不變,雞肉的供求關(guān)系受到禽流感疫情的影響時(shí),則會(huì)引起消費(fèi)者恐慌,從而導(dǎo)致雞肉購買量下降,造成雞肉需求曲線下降到D1,同時(shí)由于禽流感疫情對(duì)肉雞的健康養(yǎng)殖也會(huì)造成一定影響,養(yǎng)殖場(chǎng)(戶)也可以接觸到禽流感疫情爆發(fā)的信息,根據(jù)自身生產(chǎn)預(yù)期,調(diào)整生產(chǎn),造成雞肉供給曲線下降到S1,雞肉需求和供給均下降的情況下,雞肉的價(jià)格也會(huì)發(fā)生變化,從均衡價(jià)格P0下降到P1。

    與此同時(shí),由于肉類消費(fèi)結(jié)構(gòu)的替代性,肉類產(chǎn)品之間存在“替代效應(yīng)”,在減少雞肉消費(fèi)情況下,會(huì)使得替代品豬肉的需求隨之上升,而在豬肉供給不變的情況下,豬肉價(jià)格也會(huì)隨之提高;但隨著時(shí)間的推移,禽流感疫情逐漸被控制或消退,公眾有關(guān)禽流感疫情關(guān)注度逐漸衰退,消費(fèi)者對(duì)雞肉的購買信心回升,雞肉需求增加,養(yǎng)殖場(chǎng)(戶)抓住補(bǔ)欄時(shí)機(jī),從而帶動(dòng)了雞肉供給的增加,即供給曲線由S1上升到S2、需求曲線由D1上升到D2,雞肉價(jià)格回升到P2的水平,形成新的均衡價(jià)格,而對(duì)豬肉需求也會(huì)隨著禽流感疫情的控制或消退而恢復(fù)正常,價(jià)格也會(huì)達(dá)到新的均衡價(jià)格。

    可以看出,由于市場(chǎng)信息的不對(duì)稱,消費(fèi)者單純依靠媒體途徑所獲得的動(dòng)物疫情信息進(jìn)行消費(fèi)行為決策,具有一定的盲目性。當(dāng)消費(fèi)者獲取到疫情發(fā)生的信息后,由于周邊環(huán)境和輿論媒體等的“風(fēng)險(xiǎn)放大效應(yīng)”,容易造成恐慌,消費(fèi)者消費(fèi)信心下降,改變?cè)械南M(fèi)產(chǎn)品結(jié)構(gòu),選擇其認(rèn)為安全的和偏好的其他產(chǎn)品消費(fèi),從而導(dǎo)致某些產(chǎn)品價(jià)格發(fā)生變動(dòng),但這種改變會(huì)隨著政府的應(yīng)急措施以及信息的時(shí)效性減弱而逐漸恢復(fù)。

    已有文獻(xiàn)研究表明,學(xué)者們對(duì)疫病會(huì)導(dǎo)致相關(guān)農(nóng)產(chǎn)品價(jià)格波動(dòng)的觀點(diǎn)基本達(dá)成共識(shí),但具體的影響關(guān)系還有一定的爭(zhēng)議[5,17],造成這種差異的原因可能是不同國家的人們飲食結(jié)構(gòu)的不同以及研究對(duì)象的不同。關(guān)于我國禽流感疫情變動(dòng)與畜禽產(chǎn)品價(jià)格方面還未有學(xué)者進(jìn)行深入研究,因此還需進(jìn)一步的探討。由此,提出本文的研究問題,即禽流感疫情變動(dòng)對(duì)不同畜禽產(chǎn)品價(jià)格的沖擊大小、沖擊方向和沖擊持續(xù)時(shí)間上有何差異?隨著時(shí)間的推移有何變化?不同時(shí)期是否存在差異?

    圖1 禽流感疫情變動(dòng)影響下的雞肉價(jià)格變動(dòng)圖Fig. 1 Changes of poultry price under the in fluence of avian in fluenza

    2 研究方法與數(shù)據(jù)來源

    借助網(wǎng)絡(luò)搜索關(guān)鍵詞的搜索量分析相關(guān)事物的趨勢(shì),已被廣泛應(yīng)用在不同領(lǐng)域[18],其主要集中在疾病監(jiān)測(cè)和經(jīng)濟(jì)金融等領(lǐng)域,研究結(jié)果顯示人們的網(wǎng)絡(luò)搜索行為與其社會(huì)活動(dòng)的相關(guān)性較高,因此可以通過這種相關(guān)性代替某些難以量化的指標(biāo),并預(yù)測(cè)未來事件的發(fā)展趨勢(shì),故本文通過構(gòu)建公眾禽流感關(guān)注度指數(shù)來表示禽流感疫情的變動(dòng)及程度大小。

    2.1 公眾禽流感關(guān)注度指數(shù)的構(gòu)建

    本文借鑒輿情指數(shù)的構(gòu)建方法[19],構(gòu)建公眾禽流感關(guān)注度指數(shù),步驟如下:首先選擇若干個(gè)與社會(huì)經(jīng)濟(jì)事件密切相關(guān)的關(guān)鍵詞,然后通過互聯(lián)網(wǎng)獲取不同時(shí)段內(nèi)該關(guān)鍵詞的搜索量,并通過一定算法將搜索量進(jìn)行合成,進(jìn)而得到該指數(shù),以此來反映公眾對(duì)該關(guān)鍵詞所代表的社會(huì)經(jīng)濟(jì)事件的關(guān)注度,從而反應(yīng)此社會(huì)經(jīng)濟(jì)事件影響。

    2.1.1 選擇關(guān)鍵詞 以往研究中關(guān)鍵詞的選取可以歸納為主觀選詞法和模型選詞法兩種,其中,主觀選詞法工作量相對(duì)較小,獲取較為方便,在實(shí)際研究中也運(yùn)用較多,且取得了較好的研究效果[19]。為此,本文采用主觀選詞法結(jié)合搜索引擎相關(guān)詞推薦功能以及百度搜索需求圖譜獲取與禽流感密切相關(guān)的關(guān)鍵詞,其中,關(guān)于禽流感疫情的關(guān)鍵詞主要包括以下三類:一是“禽流感”、“禽流感病毒”和“禽流感疫情”等反映禽流感疫情名稱的關(guān)鍵詞;二是包括“H5N1”、“H5N2”、“H7N7” 和“H7N9” 等 在內(nèi)的4種已經(jīng)確認(rèn)能夠感染人且為公眾熟知的禽流感病毒型號(hào)關(guān)鍵詞,以上兩類關(guān)鍵詞用于反映禽流感疫情本身;三是“預(yù)防禽流感”、“禽流感癥狀”和“禽流感最新消息”等反映公眾態(tài)度和關(guān)注傾向的關(guān)鍵詞。因此,本文共包括以上三類、10個(gè)關(guān)鍵詞。2.1.2 信息獲取 百度作為全球最大的中文搜索引擎,其市場(chǎng)占有率和市場(chǎng)份額一直位居我國搜索引擎市場(chǎng)第一,其搜索數(shù)據(jù)具有廣泛和持續(xù)的代表性[19],因此,本文選擇百度作為數(shù)據(jù)來源。通過提取上文10個(gè)關(guān)鍵詞搜索數(shù)據(jù)并進(jìn)行加總,考慮到公眾關(guān)注度變化較快,且為獲得更多的數(shù)據(jù)樣本,本文選取周度數(shù)據(jù)進(jìn)行分析,為消除異方差,對(duì)獲取的數(shù)據(jù)取對(duì)數(shù)處理,計(jì)算公式如下:

    式中:i=1,2,…,10。qlg表示禽流感輿情指數(shù),Xi表示第i個(gè)關(guān)鍵詞搜索數(shù)據(jù)。鑒于樣本可獲得性,樣本區(qū)間選取為2013年1月至2017年3月,圖2展示了公眾禽流感關(guān)注度指數(shù)的走勢(shì)。

    通過圖2可以看出,自2013年以來我國出現(xiàn)了5次禽流感關(guān)注度高峰期(2013年4月、2014年2月、2015年1月、2016年1月和2017年3月),分別對(duì)應(yīng)于禽流感疫情的5個(gè)發(fā)病高峰期:1)2013年我國發(fā)現(xiàn)人感染H7N9禽流感病例事件,并造成約有40人喪生,此次事件短時(shí)間內(nèi)迅速引起公眾廣泛關(guān)注,禽流感關(guān)注度短時(shí)間內(nèi)大幅上升,為歷年最高水平;2)2014年年初發(fā)現(xiàn)的多起人感染禽流感事件,在兩個(gè)月時(shí)間內(nèi),人感染禽流感病毒的病例達(dá)180多例,引起社會(huì)廣泛關(guān)注;3)2015年年初禽流感關(guān)注度高峰期,此次禽流感關(guān)注高峰主要是由于當(dāng)時(shí)世界多國和地區(qū)(日本、英國等)爆發(fā)的禽流感疫情;4)2016年上半年禽流感關(guān)注高峰期,主要是由于上海、江蘇、福建和廣東等地發(fā)生多起人感染禽流感事件,而且出現(xiàn)多起患者死亡事件,引起上半年公眾關(guān)注高峰;5)2017年以來關(guān)注高峰,此次高峰主要是由于2017年以來人感染禽流感疫情病例數(shù)顯著增加所引起。

    圖2 公眾禽流感關(guān)注度指數(shù)走勢(shì)Fig. 2 Public avian in fluenza concern index

    2.2 價(jià)格數(shù)據(jù)獲取與處理

    為考察禽流感疫情變動(dòng)對(duì)畜禽產(chǎn)品價(jià)格的影響,本文選擇雞肉、豬肉、牛肉和羊肉4種我國主要畜禽消費(fèi)產(chǎn)品的價(jià)格。與公眾禽流感關(guān)注度指數(shù)一樣,本文選取畜禽產(chǎn)品價(jià)格均為周度數(shù)據(jù),數(shù)據(jù)來源于中國畜牧業(yè)信息網(wǎng)。樣本區(qū)間為2013年1月—2017年3月。在獲取基礎(chǔ)價(jià)格數(shù)據(jù)后,為消除物價(jià)因素的影響,畜禽產(chǎn)品價(jià)格數(shù)據(jù)均以2013年1月為基期的CPI指數(shù)進(jìn)行平減,再計(jì)算出畜產(chǎn)品價(jià)格收益率序列,其中雞肉價(jià)格、豬肉價(jià)格、牛肉價(jià)格和羊肉價(jià)格的收益率分別用jr、zr、nr和yr表示。2013年以來,四者價(jià)格走勢(shì)如圖3所示。

    圖3 畜禽產(chǎn)品價(jià)格走勢(shì)Fig. 3 Livestock and poultry product prices

    從圖3可以看出,2013年以來我國畜禽產(chǎn)品價(jià)格總體較為平穩(wěn),但存在短期的波動(dòng)起伏。其中,雞肉和豬肉價(jià)格波動(dòng)幅度和頻率高于牛羊肉價(jià)格波動(dòng)幅度和頻率;2013年禽流感的大范圍爆發(fā),雞肉價(jià)格一度下跌嚴(yán)重,2013年5月8號(hào)跌至15.23元/kg,相比于同年2月13日下降17.94%,此時(shí)期,豬肉、牛肉、羊肉價(jià)格均有一個(gè)下降的過程,其中豬肉下降幅度最大,牛羊肉價(jià)格下降相對(duì)較小。2017年1月份爆發(fā)的H7N9疫情,雞肉價(jià)格再次出現(xiàn)大幅下降的現(xiàn)象,相比于年初下降了10%,豬肉、牛肉也同樣伴隨著價(jià)格的波動(dòng)。

    2.3 模型構(gòu)建

    由于禽流感關(guān)注度和畜禽產(chǎn)品市場(chǎng)價(jià)格波動(dòng)起伏,可能會(huì)出現(xiàn)結(jié)構(gòu)突變的現(xiàn)象,利用傳統(tǒng)的線性或固定系數(shù)模型,不僅會(huì)影響估計(jì)效果,而且易造成系統(tǒng)性偏誤,而借助時(shí)變性和非線性的時(shí)變參數(shù)向量自回歸(TVP-VAR)模型進(jìn)行分析可以彌補(bǔ)上述不足。TVP-VAR模型由Primiceri[20]于2005年提出,并由Nakajima[21]于2011年進(jìn)行完善,不同于傳統(tǒng)的VAR模型,該模型的優(yōu)點(diǎn)是模型的系數(shù)伴隨著沖擊大小的變動(dòng)以及傳播機(jī)制的變動(dòng)而變動(dòng),可以通過時(shí)變參數(shù)捕捉到時(shí)間序列上每個(gè)時(shí)間點(diǎn)的脈沖響應(yīng),從而有效的從動(dòng)態(tài)角度刻畫變量之間的時(shí)變作用關(guān)系[22]。因此,本文借助該模型分析不同時(shí)期禽流感疫情變動(dòng)對(duì)畜禽產(chǎn)品價(jià)格造成的動(dòng)態(tài)影響。

    根據(jù)Nakajima[21]的模型構(gòu)建思路,首先,考慮如下模型:

    式中:yt是一個(gè)n×1的可觀察內(nèi)生向量;ct是一個(gè)n×1的時(shí)變常數(shù)項(xiàng)向量;Bi,t,i=1,…,k是n×n的時(shí)變系數(shù)向量;μt衡量擁有協(xié)方差矩陣Ωt的不可觀察的沖擊向量,Ωt由以下方程表示其中,矩陣At具有下三角的矩陣形式:

    Σt為對(duì)稱矩陣 :Σt=diag(σ1,t,…,σn,t),通過把Bi,t,i=1,…,k轉(zhuǎn)為向量 βt, (1)式可以轉(zhuǎn)化為 :

    式中:I3為三維單位矩陣,Σβ、Σα和Σh均為正定矩陣。本文中模型的估計(jì)通過貝葉斯方法來完成,參數(shù)的后驗(yàn)數(shù)值估計(jì)通過馬爾可夫鏈蒙特卡洛(MCMC)方法完成。

    3 實(shí)證分析

    3.1 模型估計(jì)

    3.1.1 平穩(wěn)性檢驗(yàn) 在模型估計(jì)前,首先要對(duì)各變量進(jìn)行平穩(wěn)性檢驗(yàn),以免出現(xiàn)偽回歸。通過ADF檢驗(yàn)對(duì)各變量進(jìn)行平穩(wěn)性檢驗(yàn)(見表1),結(jié)果表明,5個(gè)變量均是平穩(wěn)的,其中,雞肉、豬肉、牛肉和羊肉價(jià)格收益率均在1%水平下平穩(wěn),公眾禽流感關(guān)注度指數(shù)在5%水平下平穩(wěn)。說明各變量均不存在單位根,變量可用于構(gòu)建TVP-VAR模型,能夠保證估計(jì)的有效性。

    表1 變量單位根檢驗(yàn)Table 1 Unit root test resutls

    3.1.2 模型估計(jì) 為判定模型的滯后階數(shù),借助一般VAR模型最優(yōu)滯后階數(shù)的選擇辦法,根據(jù)FPE、AIC、SC和HQ最小值準(zhǔn)則,可知模型最優(yōu)滯后期數(shù)為2(見表2)。在利用MCMC法進(jìn)行模擬前需先對(duì)參數(shù)賦初值,利用MCMC算法抽取M=10 000次樣本,并舍棄初始的1 000次抽樣,從而得到模型的有效樣本。表3是待估參數(shù)后驗(yàn)分布的均值、標(biāo)準(zhǔn)差、95%的置信區(qū)間以及收斂統(tǒng)計(jì)量。其中,參數(shù)的后驗(yàn)均值都位于95%的置信區(qū)間,Geweke診斷概率均大于10%,說明后驗(yàn)分布收斂于零的假設(shè)無法拒絕,而且從無效因子來看,除了(Σh)2的無效因子較大外,其他參數(shù)的無效因子均處于合理水平,說明本文模型參數(shù)的估計(jì)結(jié)果較為穩(wěn)健(見表3)。

    表2 模型滯后階數(shù)判定Table 2 Determination of the model lag

    表3 模型估計(jì)結(jié)果Table 3 Estimation results of the model

    3.2 時(shí)變脈沖響應(yīng)結(jié)果

    3.2.1 不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的等間隔脈沖響應(yīng) 圖4展示了不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的動(dòng)態(tài)脈沖響應(yīng)。首先,從圖4上可以看出畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)具有明顯的時(shí)變性特征,說明不同時(shí)期隨著禽流感疫情的變化,其對(duì)畜禽產(chǎn)品價(jià)格造成的影響不同。其次,從不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)脈沖響應(yīng)程度上看,雞肉價(jià)格受到的影響最大,其次為豬肉價(jià)格,牛羊肉價(jià)格受到的影響相對(duì)較小,這與以往研究結(jié)果基本一致[17,23],可能原因是禽流感爆發(fā)時(shí)公眾禽流感關(guān)注度快速上升,導(dǎo)致雞肉消費(fèi)大幅減少引起價(jià)格出現(xiàn)較大下降,豬肉作為雞肉最主要的替代品同樣受到較大的影響,而牛羊肉消費(fèi)群體較為固定,兩者與雞肉替代性較小,因此受到的影響也最小。

    再次,從不同滯后期畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)上看,滯后2期時(shí)畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的響應(yīng)最大,滯后4期的脈沖響應(yīng)明顯減弱,滯后6期的脈沖響應(yīng)基本趨于0,說明隨著滯后期的增加禽流感疫情變動(dòng)對(duì)畜禽產(chǎn)品價(jià)格的影響快速減弱,持續(xù)期并不長。最后,從不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)脈沖響應(yīng)的方向上看,短期內(nèi)(滯后2期)主要以負(fù)向響應(yīng)為主,雞肉和牛羊肉在滯后4期和6期時(shí)對(duì)禽流感疫情變動(dòng)的響應(yīng)也以負(fù)向響應(yīng)為主,值得注意的是,豬肉在滯后4期和6期時(shí)對(duì)禽流感疫情變動(dòng)的響應(yīng)主要以正向響應(yīng)為主,這與周海文等[24]的研究類似,解釋這一現(xiàn)象的原因可能是:禽流感疫情發(fā)生時(shí),公眾關(guān)于禽流感關(guān)注迅速增加,短期內(nèi)造成市場(chǎng)恐慌,降低畜禽產(chǎn)品消費(fèi),引起畜禽產(chǎn)品價(jià)格下降;但之后作為雞肉的主要替代品,豬肉消費(fèi)開始增加引起價(jià)格的上升,而牛羊肉由于與雞肉替代性較小并未出現(xiàn)價(jià)格上升。

    圖4 不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的動(dòng)態(tài)脈沖響應(yīng)Fig. 4 Dynamic impulse response of different livestock product prices to avian in fluenza epidemic

    3.2.2 不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的動(dòng)態(tài)脈沖響應(yīng) 具體來看不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的動(dòng)態(tài)脈沖響應(yīng)。首先,從圖4(a)可以看出,雞肉價(jià)格對(duì)禽流感疫情變動(dòng)的響應(yīng)在滯后2期、4期和6期時(shí)均為負(fù)向;自2013年以來經(jīng)歷了5次響應(yīng)高峰期,分別對(duì)應(yīng)于2013年以來的禽流感疫情關(guān)注度高峰期;其中2013年3月和2017年3月兩次脈沖響應(yīng)程度較高,主要原因是2013年3月初次發(fā)現(xiàn)人感染H7N9禽流感事件,引起市場(chǎng)恐慌,導(dǎo)致雞肉價(jià)格大幅度下降,而2017年3月的高峰期是由于2017年年初人感染禽流感疫情病例數(shù)顯著增加所引起,導(dǎo)致雞肉消費(fèi)下降進(jìn)而引起價(jià)格下降。從圖4(b)豬肉對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)上看,由于與雞肉之間存在較強(qiáng)的替代性,其脈沖響應(yīng)走勢(shì)與雞肉較為相似,同樣經(jīng)歷了幾次響應(yīng)高峰期。

    從圖4(c)可以看出,自2013年下半年以來,牛肉對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)均為負(fù)向,但同樣經(jīng)歷了幾個(gè)波動(dòng)的周期,在2014年初以及2015年初存在兩個(gè)較大的響應(yīng)高峰期,出現(xiàn)這一現(xiàn)象可能有兩個(gè)原因,一是與禽流感關(guān)注度波動(dòng)有關(guān),隨著關(guān)注度的波動(dòng)起伏,對(duì)牛肉價(jià)格影響也波動(dòng)起伏,在禽流感關(guān)注度高峰期對(duì)牛肉價(jià)格造成較大影響;二是與牛肉的季節(jié)性消費(fèi)有關(guān),第一和第四季度為消費(fèi)旺季,當(dāng)此時(shí)期發(fā)生禽流感疫情時(shí),公眾關(guān)注度增加,市場(chǎng)恐慌程度增加引起消費(fèi)量下降程度較大,導(dǎo)致價(jià)格下降程度較大。羊肉價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)可以劃分為兩個(gè)階段(見圖4(d)),一是2013年—2015年小幅波動(dòng)期,此階段羊肉價(jià)格對(duì)禽流感疫情變動(dòng)脈沖響應(yīng)較??;二是2015年以來的大幅波動(dòng)期,此階段羊肉價(jià)格對(duì)禽流感疫情變動(dòng)脈沖響應(yīng)較大,同時(shí)還可以發(fā)現(xiàn),2015年以來其他畜禽產(chǎn)品價(jià)格脈沖響應(yīng)也處在較高水平,造成這一現(xiàn)象的原因可能是隨著互聯(lián)網(wǎng)、多媒體等信息技術(shù)的快速發(fā)展,尤其是最近兩年智能手機(jī)、移動(dòng)終端的快速普及,使得公眾信息接收速度不斷加快,禽流感疫情的發(fā)生極易引起公眾關(guān)注,造成市場(chǎng)恐慌進(jìn)而導(dǎo)致畜禽產(chǎn)品市場(chǎng)不穩(wěn)定程度增加。

    3.2.3 不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的時(shí)點(diǎn)脈沖響應(yīng) 圖5展示了不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的時(shí)點(diǎn)脈沖響應(yīng),本文主要選取2013年第15周、2014年第5周和2016年第3周做時(shí)點(diǎn)脈沖響應(yīng)分析。從雞肉價(jià)格對(duì)禽流感疫情變動(dòng)在三個(gè)時(shí)點(diǎn)的脈沖響應(yīng)上看(圖5(a)),三個(gè)時(shí)點(diǎn)上雞肉價(jià)格對(duì)禽流感疫情變動(dòng)均產(chǎn)生了較大的負(fù)向響應(yīng),這與前文分析一致,但產(chǎn)生最大響應(yīng)期數(shù)不同,2016年第3周時(shí)的脈沖響應(yīng)峰值為當(dāng)期,而2013年和2014年兩個(gè)時(shí)點(diǎn)上響應(yīng)峰值均在第一期,這主要是由于近兩年互聯(lián)網(wǎng)信息傳播速度加快,人們關(guān)于禽流感疫情的響應(yīng)速度同樣加快,導(dǎo)致價(jià)格響應(yīng)時(shí)間縮短。圖5(b)展示了豬肉價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng),可以看出,豬肉價(jià)格對(duì)禽流感疫情變動(dòng)在短期為負(fù)向響應(yīng),之后出現(xiàn)反轉(zhuǎn),轉(zhuǎn)為正向響應(yīng),這進(jìn)一步證實(shí)了前文分析結(jié)果。

    從圖5(c)牛肉價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)上看,2013年第15周以正向響應(yīng)為主,2014年和2016年兩個(gè)時(shí)點(diǎn)均為負(fù)向響應(yīng)。羊肉價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)上看(圖5(d)),2013年第15周羊肉價(jià)格對(duì)禽流感疫情變動(dòng)為微弱正向響應(yīng),2014年第5周為正負(fù)交替的響應(yīng),2016年第3周為負(fù)向響應(yīng),這與前文結(jié)果一致。同時(shí),從圖上可以看出,豬肉和牛羊肉價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)的持續(xù)時(shí)間基本在6周左右,而雞肉對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)的持續(xù)時(shí)間相對(duì)較長,在9周左右。

    圖5 不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的時(shí)點(diǎn)脈沖響應(yīng)Fig. 5 Time point impulse response of different livestock product prices to avian in fluenza epidemic

    3.3 穩(wěn)健性檢驗(yàn)

    為進(jìn)一步驗(yàn)證前文實(shí)證結(jié)果的穩(wěn)健性,本文通過兩種方式來進(jìn)行,首先,通過調(diào)整估計(jì)模型的方式進(jìn)行穩(wěn)健性檢驗(yàn),通過采用BVAR模型分析禽流感疫情變動(dòng)對(duì)不同畜禽產(chǎn)品價(jià)格的影響,得到不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)(圖6)。從圖6可以看出,雞肉、牛羊肉價(jià)格對(duì)禽流感疫情變動(dòng)均為負(fù)向響應(yīng),而豬肉價(jià)格對(duì)禽流感疫情變動(dòng)則存在明顯的反轉(zhuǎn)效應(yīng),即短期為負(fù)向響應(yīng),之后轉(zhuǎn)為正向響應(yīng),該結(jié)果進(jìn)一步驗(yàn)證了前文實(shí)證結(jié)果。

    圖6 不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)脈沖響應(yīng)Fig. 6 Impulse response of different livestock product prices to avian in fluenza epidemic

    其次,通過構(gòu)建禽流感疫情發(fā)生變量替換上文禽流感代理變量,再次對(duì)模型進(jìn)行估計(jì),進(jìn)一步檢驗(yàn)上文結(jié)論的穩(wěn)健性。禽流感疫情變量數(shù)據(jù)來源于農(nóng)業(yè)部公布的《獸醫(yī)公報(bào)》,計(jì)算各地禽流感發(fā)病數(shù)、死亡數(shù)和撲殺數(shù)之和,并對(duì)其取對(duì)數(shù),最終獲取禽流感疫情發(fā)生變量(由于《獸醫(yī)公報(bào)》均為月度數(shù)據(jù),因此將各月所有周度數(shù)據(jù)均與月度數(shù)據(jù)相對(duì)應(yīng)進(jìn)行模型估計(jì))。并再次采用BVAR模型分析禽流感疫情變動(dòng)對(duì)不同畜禽產(chǎn)品價(jià)格的影響,得到不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)的脈沖響應(yīng)(圖7)。可以看出,脈沖響應(yīng)結(jié)果與代理變量估計(jì)結(jié)果較為一致,由此可以證明本文實(shí)證結(jié)果的穩(wěn)健性。

    4 結(jié)論與建議

    4.1 結(jié)論

    為分析禽流感疫情變動(dòng)是否會(huì)影響畜禽產(chǎn)品價(jià)格以及影響的大小,本文通過構(gòu)建禽流感疫情變動(dòng)代理變量——公眾禽流感關(guān)注度指數(shù),選取雞肉、豬肉、牛肉和羊肉價(jià)格,運(yùn)用2013年1月至2017年3月的周度數(shù)據(jù),采用時(shí)變參數(shù)向量自回歸模型(TVP-VAR)分析了禽流感疫情變動(dòng)對(duì)不同畜禽產(chǎn)品價(jià)格波動(dòng)的動(dòng)態(tài)影響,結(jié)果表明,禽流感疫情變動(dòng)會(huì)對(duì)畜禽產(chǎn)品價(jià)格造成影響,且影響具有時(shí)變性,禽流感疫情變動(dòng)對(duì)不同畜禽產(chǎn)品價(jià)格的沖擊大小、沖擊方向和沖擊持續(xù)時(shí)間上存在差異。其中,禽流感疫情對(duì)雞肉價(jià)格影響最大,其次為豬肉價(jià)格,牛羊肉價(jià)格受到的影響相對(duì)較?。磺萘鞲幸咔榈淖儎?dòng)對(duì)雞肉和牛羊肉價(jià)格以負(fù)向影響為主,而對(duì)豬肉價(jià)格的影響存在反轉(zhuǎn)效應(yīng),即短期存在負(fù)向影響,中長期存在正向影響;由于消費(fèi)季節(jié)性等特點(diǎn),禽流感疫情對(duì)畜禽產(chǎn)品價(jià)格的影響存在差異。上述結(jié)果通過調(diào)整估計(jì)方法以及調(diào)換變量等穩(wěn)健性檢驗(yàn)均得到了驗(yàn)證,從而保證了本文實(shí)證結(jié)果的穩(wěn)健性。

    4.2 建議

    基于本研究結(jié)論,認(rèn)為政府有關(guān)部門應(yīng)進(jìn)一步加大對(duì)禽流感疫情的監(jiān)測(cè)力度以及提高疫情處理能力。首先,應(yīng)進(jìn)一步完善禽流感預(yù)警系統(tǒng),加強(qiáng)對(duì)禽流感的監(jiān)測(cè)力度,建立畜禽產(chǎn)品價(jià)格聯(lián)動(dòng)監(jiān)測(cè)機(jī)制,防患于未然。其次,除關(guān)注禽流感疫情本身發(fā)生情況之外,一方面應(yīng)通過多渠道密切關(guān)注公眾禽流感輿情動(dòng)態(tài)變動(dòng)情況,尤其是借助互聯(lián)網(wǎng)渠道加強(qiáng)輿情監(jiān)督,做好輿情引導(dǎo);另一方面,加強(qiáng)禽流感疫情科學(xué)知識(shí)的宣傳,增強(qiáng)公眾對(duì)禽流感的認(rèn)知,減少公眾禽流感恐慌。最后,在疫情發(fā)生時(shí),除應(yīng)重點(diǎn)關(guān)注肉雞產(chǎn)業(yè)市場(chǎng)變化外,還應(yīng)關(guān)注豬肉產(chǎn)業(yè)的市場(chǎng)變化,做到積極正確的消費(fèi)及生產(chǎn)引導(dǎo)。

    圖7 不同畜禽產(chǎn)品價(jià)格對(duì)禽流感疫情變動(dòng)代理變量脈沖響應(yīng)Fig. 7 Impulse response of different livestock product prices to proxy variable of avian in fluenza epidemic

    參考文獻(xiàn):

    [1] Suder G, Inthavong S. New health risks and sociocultural contexts: Bird flu impacts on consumers and poultry businesses in Lao PDR[J]. Risk Analysis, 2008, 28(1): 1-12.

    [2] 王秀清, Weldegebriel H T, Rayner A J. 縱向關(guān)聯(lián)市場(chǎng)間的價(jià)格傳遞[J]. 經(jīng)濟(jì)學(xué)(季刊), 2007, 6(3): 885-898.Wang X Q, Weldegebriel H T, Rayner A J. Price transmission in vertically related markets[J]. China Economic Quarterly, 2007,6(3): 885-898.

    [3] Beach R H, Kuchler F, Leibtag E, et al. The effects of avian in fluenza news on consumer purchasing behavior: A case study of Italian consumers’ retail purchases[R]. Washington, 2008.

    [4] 周力, 劉常瑜. 禽流感風(fēng)險(xiǎn)下肉雞產(chǎn)業(yè)價(jià)格縱橫傳導(dǎo)研究[J].統(tǒng)計(jì)與決策, 2016(17): 93-96.Zhou L, Liu C Y. Vertical and horizontal transmission of broiler industry price under avian in fluenza risk[J]. Statistics & Decision,2016(17): 93-96.

    [5] 夏云. H7N9疫情對(duì)家禽產(chǎn)業(yè)的影響[D]. 廣州: 華南農(nóng)業(yè)大學(xué),2016.Xia Y. An empirical study on the impact of H7N9 influenza on poultry industry[D]. Guangzhou: South China Agricultural University, 2016.

    [6] 黃澤穎, 王濟(jì)民. 高致病性禽流感對(duì)我國肉雞產(chǎn)業(yè)的影響[J].中國農(nóng)業(yè)科技導(dǎo)報(bào), 2016, 18(1): 189-199.Huang Z Y, Wang J M. Economic impact of highly pathogenic avian influenza on broiler industry in China[J]. Journal of Agricultural Science and Technology, 2016, 18(1): 189-199.

    [7] 于樂榮, 李小云, 汪力斌, 等. 禽流感發(fā)生對(duì)家禽養(yǎng)殖農(nóng)戶的經(jīng)濟(jì)影響評(píng)估——基于兩期面板數(shù)據(jù)的分析[J]. 中國農(nóng)村經(jīng)濟(jì), 2009(7): 12-19, 30.Yu L R, Li X Y, Wang L B , et al. An assessment of the economic impacts of occurrence of avian influenza on poultry raising famers’ households: An analysis based on two period panel data[J]. Chinese Rural Economy, 2009(7): 12-19, 30.

    [8] 梅付春, 張陸彪. 禽流感疫區(qū)散養(yǎng)戶對(duì)撲殺補(bǔ)償政策配合意愿的實(shí)證分析[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題, 2008(S1): 173-177.Mei F C, Zhang L B. Empirical analysis on the small’s willingness to corporate with the compensation policy for culling in the bird flu-infected areas[J]. Issues in Agricultural Economy, 2008(S1):173-177.

    [9] 張莉琴, 康小瑋, 林萬龍. 高致病性禽流感疫情防制措施造成的養(yǎng)殖戶損失及政府補(bǔ)償分析[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題, 2009(12):28-33.Zhang L Q, Kang X W, Lin W L. Loss of livestock farms and government compensatory in the control of HPAI[J]. Issues in Agricultural Economy, 2009(12): 28-33.

    [10] Serra T. Food scare crises and price volatility: The case of the BSE in Spain[J]. Food Policy, 2011, 36(2): 179-185.

    [11] 劉明月, 陸遷. 禽流感疫情沖擊下疫區(qū)養(yǎng)殖戶生產(chǎn)恢復(fù)行為研究——以寧夏中衛(wèi)沙坡區(qū)為例[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題, 2016(5): 40-50, 111.Liu M Y, Lu Q. Study of farmers’ productive recovery behavior in avian-infected area under the shocks of avian in fluenza: A case study of Zhongwei in Ningxia[J]. Issues in Agricultural Economy,2016(5): 40-50, 111.

    [12] Roux A. Food scare crises and developing countries: The impact of avian in fluenza on vertical price transmission in the Egyptian poultry sector[J]. Food Policy, 2012, 37(3): 264-274.

    [13] Serra T. Vertical price transmission in the Turkish poultry market:The avian in fluenza crisis[J]. Applied Economics, 2015, 47(11):1106-1117.

    [14] 孫毅, 呂本富. 網(wǎng)絡(luò)搜索與經(jīng)濟(jì)行為相關(guān)性研究綜述[J]. 管理評(píng)論, 2011, 23(7): 72-77.Sun Y, Lü B F. A review of the researches on the correlation between internet search and economic behavior[J]. Management Review, 2011, 23(7): 72-77.

    [15] Rech J. Discovering trends in software engineering with google trend[J]. Acm Sigsoft Software Engineering Notes, 2007, 32(2): 1-2.

    [16] 張繼德, 廖微, 張榮武. 普通投資者關(guān)注對(duì)股市交易的量價(jià)影響——基于百度指數(shù)的實(shí)證研究[J]. 會(huì)計(jì)研究, 2014(8): 52-59, 97.Zhang J D, Liao W, Zhang R W. The effect of ordinary investors’attention on volume an price of stock market: empirical evidence based on baidu index [J]. Accounting Research, 2014(8): 52-59, 97.

    [17] 黃德林, 董蕾, 王濟(jì)民. 禽流感對(duì)養(yǎng)禽業(yè)和農(nóng)民收入的影響[J].農(nóng)業(yè)經(jīng)濟(jì)問題, 2004(6): 21-25, 79.Huang D L, Dong L, Wang J M. Influence of bird flu on the poultry industry and farmers’ income[J]. Issues in Agricultural Economy, 2004(6): 21-25, 79.

    [18] 盧洪濤, 李綱. 網(wǎng)絡(luò)搜索關(guān)鍵詞時(shí)序變化特征研究——以H7N9禽流感關(guān)鍵詞實(shí)驗(yàn)為例[J]. 情報(bào)雜志, 2014, 33(11): 175-180.Lu H T, Li G. Research on timing variation characteristics of keywords in web search keywords based on a case study about H7N9 avian influenza keywords experiment[J]. Journal of Intelligence, 2014, 33(11): 175-180.

    [19] 徐映梅, 高一銘. 基于互聯(lián)網(wǎng)大數(shù)據(jù)的CPI輿情指數(shù)構(gòu)建與應(yīng)用——以百度指數(shù)為例[J]. 數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究, 2017,34(1): 94-112.Xu Y M, Gao Y M. Construction of the public opinion index of CPI based on the internet big data[J]. The Journal of Quantitative& Technical Economics, 2017, 34(1): 94-112.

    [20] Primiceri. Time varying structural vector autoregressions and monetary policy[J]. Review of Economic Studies, 2005, 72(3):821-852.

    [21] Nakajima J. Time-varying parameter VAR model with stochastic volatility: An overview of methodology and empirical applications[J]. Jouchi Nakajima, 2011, 29: 107-142.

    [22] 錢燕, 萬解秋. 貨幣供應(yīng)、通貨膨脹與經(jīng)濟(jì)增長的互動(dòng)關(guān)系研究——基于時(shí)變參數(shù)VAR模型的實(shí)證檢驗(yàn)[J]. 西安財(cái)經(jīng)學(xué)院學(xué)報(bào), 2014, 27(1): 5-10.Qian Y, Wan J Q. Research on the interaction between money supply, in flation and economic growth: Based on the time-varying parameter VAR model[J]. Journal of Xi’an University of Finance and Economics, 2014, 27(1): 5-10.

    [23] Ishida T, Ishikawa N, Fukushige M. Impact of BSE and bird flu on consumers’ meat demand in Japan[J]. Applied Economics , 2010,42(1): 49-56.

    [24] 周海文, 王勁松, 王銳, 等. 外部沖擊對(duì)肉產(chǎn)品價(jià)格的影響——以疫病為例[J]. 世界農(nóng)業(yè), 2014(11): 76-82.Zhou H W, Wang J S, Wang R, et al. The impact of external shocks on meat products prices: A case study of epidemic diseases[J].World Agriculture, 2014(11): 76-82.

    猜你喜歡
    疫情模型
    一半模型
    戰(zhàn)疫情
    重要模型『一線三等角』
    抗疫情 顯擔(dān)當(dāng)
    疫情中的我
    疫情當(dāng)前 警察不退
    北極光(2020年1期)2020-07-24 09:04:04
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    待疫情散去 春暖花開
    文苑(2020年4期)2020-05-30 12:35:48
    37°女人(2020年5期)2020-05-11 05:58:52
    3D打印中的模型分割與打包
    国产成人a区在线观看| 欧美三级亚洲精品| 一本久久中文字幕| АⅤ资源中文在线天堂| 中文字幕av成人在线电影| 怎么达到女性高潮| 亚洲人成伊人成综合网2020| 国产精品久久久久久人妻精品电影| 午夜福利18| 国产毛片a区久久久久| 国产亚洲av嫩草精品影院| 国产爱豆传媒在线观看| 级片在线观看| 免费搜索国产男女视频| 99在线人妻在线中文字幕| 午夜免费成人在线视频| 欧美一级a爱片免费观看看| 十八禁人妻一区二区| 99国产精品一区二区蜜桃av| 国产精品三级大全| 又黄又爽又刺激的免费视频.| 极品教师在线免费播放| 亚洲美女视频黄频| 国产极品精品免费视频能看的| 国产精品电影一区二区三区| 国产精品美女特级片免费视频播放器| av在线天堂中文字幕| 亚洲精品色激情综合| 亚洲欧美日韩无卡精品| 免费看a级黄色片| 两个人的视频大全免费| 国产亚洲精品综合一区在线观看| 亚洲一区高清亚洲精品| 男女下面进入的视频免费午夜| 高潮久久久久久久久久久不卡| 免费看美女性在线毛片视频| 小说图片视频综合网站| 国产日本99.免费观看| 88av欧美| 免费高清视频大片| 深夜a级毛片| 一本综合久久免费| 真人做人爱边吃奶动态| 亚洲不卡免费看| 99精品久久久久人妻精品| 欧美激情在线99| 麻豆国产97在线/欧美| 国产精品久久久久久亚洲av鲁大| 国产欧美日韩精品一区二区| 欧美激情在线99| 麻豆久久精品国产亚洲av| 国产av一区在线观看免费| bbb黄色大片| av视频在线观看入口| 日本精品一区二区三区蜜桃| 无遮挡黄片免费观看| 我的老师免费观看完整版| 日日夜夜操网爽| 久久久久精品国产欧美久久久| 97超级碰碰碰精品色视频在线观看| 在线观看午夜福利视频| 90打野战视频偷拍视频| 国产精品久久久久久精品电影| 真人做人爱边吃奶动态| 无人区码免费观看不卡| 高清毛片免费观看视频网站| 亚洲最大成人av| 亚洲,欧美精品.| 久久精品国产自在天天线| 国产精品一区二区性色av| 久久久久久久久久成人| 99精品久久久久人妻精品| 性插视频无遮挡在线免费观看| 久久人人爽人人爽人人片va | 看免费av毛片| 男人狂女人下面高潮的视频| 在线观看午夜福利视频| 久久午夜福利片| 狠狠狠狠99中文字幕| 超碰av人人做人人爽久久| АⅤ资源中文在线天堂| 深夜精品福利| 黄片小视频在线播放| 精品久久久久久,| 首页视频小说图片口味搜索| 亚洲熟妇熟女久久| 动漫黄色视频在线观看| 男女那种视频在线观看| 亚洲av中文字字幕乱码综合| 69av精品久久久久久| 久久国产乱子伦精品免费另类| 亚洲专区中文字幕在线| 亚洲天堂国产精品一区在线| 三级毛片av免费| 精品一区二区三区视频在线| 国产三级黄色录像| 97超视频在线观看视频| 日本与韩国留学比较| 日韩精品青青久久久久久| 全区人妻精品视频| 国产av麻豆久久久久久久| 亚洲美女搞黄在线观看 | 免费人成视频x8x8入口观看| 动漫黄色视频在线观看| 国产毛片a区久久久久| 又爽又黄a免费视频| 国产麻豆成人av免费视频| 国产真实乱freesex| 亚洲五月婷婷丁香| 国模一区二区三区四区视频| 午夜亚洲福利在线播放| 亚洲狠狠婷婷综合久久图片| 亚洲国产欧美人成| 少妇被粗大猛烈的视频| 婷婷亚洲欧美| 特级一级黄色大片| 99热这里只有是精品50| 亚洲一区高清亚洲精品| 欧美成人免费av一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 一a级毛片在线观看| 婷婷精品国产亚洲av在线| 动漫黄色视频在线观看| 成年版毛片免费区| 黄色日韩在线| 欧美乱色亚洲激情| 国产成+人综合+亚洲专区| 亚洲中文日韩欧美视频| 窝窝影院91人妻| a级一级毛片免费在线观看| 亚洲在线观看片| 国产精品野战在线观看| 老司机深夜福利视频在线观看| 国产亚洲精品av在线| 伦理电影大哥的女人| 欧美激情久久久久久爽电影| 国产亚洲精品久久久久久毛片| 人妻夜夜爽99麻豆av| 亚洲第一区二区三区不卡| 久久精品国产99精品国产亚洲性色| 国产人妻一区二区三区在| 亚洲av熟女| 男女床上黄色一级片免费看| 最近视频中文字幕2019在线8| av天堂在线播放| 中文亚洲av片在线观看爽| 午夜两性在线视频| 婷婷亚洲欧美| 小蜜桃在线观看免费完整版高清| 欧美潮喷喷水| 丁香六月欧美| 特级一级黄色大片| 内射极品少妇av片p| 在线十欧美十亚洲十日本专区| 亚洲 欧美 日韩 在线 免费| 久久热精品热| 精品日产1卡2卡| 精品久久久久久成人av| 国产亚洲av嫩草精品影院| 午夜亚洲福利在线播放| 无遮挡黄片免费观看| 91午夜精品亚洲一区二区三区 | 欧美+日韩+精品| 精品国内亚洲2022精品成人| 日韩欧美在线乱码| 免费av观看视频| 国产亚洲精品综合一区在线观看| 午夜日韩欧美国产| 国产午夜精品论理片| 成人一区二区视频在线观看| 首页视频小说图片口味搜索| 久99久视频精品免费| 国产精品美女特级片免费视频播放器| 欧美xxxx性猛交bbbb| 男人狂女人下面高潮的视频| 别揉我奶头 嗯啊视频| 国产欧美日韩精品一区二区| 国产精品99久久久久久久久| 在线观看免费视频日本深夜| 在线十欧美十亚洲十日本专区| 国产毛片a区久久久久| 午夜日韩欧美国产| 成人永久免费在线观看视频| 精品不卡国产一区二区三区| 男人舔奶头视频| 亚洲国产精品sss在线观看| 国产真实伦视频高清在线观看 | 亚洲精品在线美女| 亚洲真实伦在线观看| 亚洲综合色惰| 男女床上黄色一级片免费看| 午夜福利在线观看吧| 日韩亚洲欧美综合| 中国美女看黄片| 亚洲最大成人手机在线| 能在线免费观看的黄片| 最新在线观看一区二区三区| 亚洲在线自拍视频| xxxwww97欧美| 久久精品国产99精品国产亚洲性色| 精品熟女少妇八av免费久了| 免费黄网站久久成人精品 | 免费看a级黄色片| av在线蜜桃| 日韩免费av在线播放| 精品福利观看| 久久欧美精品欧美久久欧美| 欧美三级亚洲精品| 国产野战对白在线观看| 欧美黑人欧美精品刺激| 一个人免费在线观看的高清视频| 十八禁人妻一区二区| 国产精品伦人一区二区| 欧洲精品卡2卡3卡4卡5卡区| 国产高清三级在线| 国产午夜精品论理片| 在线观看66精品国产| 天堂√8在线中文| 99久久成人亚洲精品观看| 亚洲av电影在线进入| 69人妻影院| 69av精品久久久久久| 免费观看的影片在线观看| 俺也久久电影网| 在线免费观看不下载黄p国产 | 亚洲最大成人手机在线| 亚洲成人久久爱视频| 国产单亲对白刺激| 久久久久久久亚洲中文字幕 | 99热只有精品国产| 国产一区二区亚洲精品在线观看| 日本免费a在线| 欧美一区二区亚洲| 99国产综合亚洲精品| 一级毛片久久久久久久久女| 亚洲精华国产精华精| 国产亚洲欧美98| 日本一本二区三区精品| 亚洲av二区三区四区| 一进一出抽搐gif免费好疼| 日韩欧美三级三区| 岛国在线免费视频观看| 一个人免费在线观看的高清视频| 精品免费久久久久久久清纯| 99riav亚洲国产免费| 99久久精品国产亚洲精品| 高潮久久久久久久久久久不卡| .国产精品久久| 国语自产精品视频在线第100页| 久久热精品热| 亚洲美女视频黄频| 午夜福利免费观看在线| 午夜免费激情av| 精品久久久久久成人av| 91在线观看av| 午夜福利高清视频| 亚洲av.av天堂| 综合色av麻豆| 18美女黄网站色大片免费观看| 99久国产av精品| 久久久色成人| 国产成人影院久久av| 欧美xxxx黑人xx丫x性爽| 一夜夜www| 午夜福利在线观看吧| 久久九九热精品免费| 老司机午夜福利在线观看视频| 日本撒尿小便嘘嘘汇集6| 欧美最黄视频在线播放免费| 国内精品久久久久精免费| 91狼人影院| 1000部很黄的大片| 黄色女人牲交| 国产精品久久久久久人妻精品电影| 成年女人毛片免费观看观看9| 在线播放无遮挡| 九九在线视频观看精品| 成人三级黄色视频| 国产精品日韩av在线免费观看| 夜夜爽天天搞| 国产精品嫩草影院av在线观看 | 免费一级毛片在线播放高清视频| 乱码一卡2卡4卡精品| 国内精品久久久久精免费| 欧美激情国产日韩精品一区| 91九色精品人成在线观看| 国产精品乱码一区二三区的特点| 欧美日韩黄片免| 国产伦精品一区二区三区视频9| 亚洲精品在线观看二区| 国产蜜桃级精品一区二区三区| 午夜福利在线在线| 最新中文字幕久久久久| 天堂√8在线中文| av视频在线观看入口| 国产69精品久久久久777片| 国产高清三级在线| 国语自产精品视频在线第100页| 精品国产亚洲在线| 一区二区三区高清视频在线| 欧美在线黄色| 综合色av麻豆| 亚洲国产精品成人综合色| 99热这里只有是精品在线观看 | 深夜a级毛片| 欧美一区二区国产精品久久精品| 欧美成人一区二区免费高清观看| 久久草成人影院| 搞女人的毛片| 国产男靠女视频免费网站| 精华霜和精华液先用哪个| 久久精品影院6| av视频在线观看入口| 亚洲成a人片在线一区二区| 天堂影院成人在线观看| 变态另类成人亚洲欧美熟女| 婷婷六月久久综合丁香| 露出奶头的视频| 成年人黄色毛片网站| 国产 一区 欧美 日韩| 欧美日本视频| 少妇人妻一区二区三区视频| 亚洲最大成人av| 精品人妻熟女av久视频| 国产黄a三级三级三级人| 99riav亚洲国产免费| 色综合欧美亚洲国产小说| 一个人看视频在线观看www免费| 午夜两性在线视频| 97人妻精品一区二区三区麻豆| 村上凉子中文字幕在线| 亚洲av美国av| 制服丝袜大香蕉在线| 精品欧美国产一区二区三| 女人被狂操c到高潮| 又黄又爽又刺激的免费视频.| 熟妇人妻久久中文字幕3abv| 在线a可以看的网站| 国产精品久久久久久亚洲av鲁大| 亚洲av日韩精品久久久久久密| 精品熟女少妇八av免费久了| 国产高清三级在线| 欧美精品啪啪一区二区三区| 日韩中字成人| 国产高清有码在线观看视频| 日本黄色片子视频| 成人精品一区二区免费| 2021天堂中文幕一二区在线观| 亚洲精品粉嫩美女一区| 真人一进一出gif抽搐免费| 欧美潮喷喷水| 欧美+日韩+精品| 欧美日韩国产亚洲二区| 别揉我奶头 嗯啊视频| 国产精品亚洲av一区麻豆| 成人美女网站在线观看视频| av在线老鸭窝| av天堂中文字幕网| 午夜福利欧美成人| 亚洲av第一区精品v没综合| 久久久国产成人精品二区| 欧美乱色亚洲激情| 少妇人妻一区二区三区视频| 亚洲狠狠婷婷综合久久图片| 亚洲人成伊人成综合网2020| 亚洲av中文字字幕乱码综合| 熟妇人妻久久中文字幕3abv| 精品久久久久久,| 欧美性猛交╳xxx乱大交人| av女优亚洲男人天堂| 免费看日本二区| 国产成+人综合+亚洲专区| 91字幕亚洲| 国产黄a三级三级三级人| 欧美xxxx黑人xx丫x性爽| 免费观看人在逋| 成人亚洲精品av一区二区| 国产精品亚洲一级av第二区| 男女床上黄色一级片免费看| 熟女人妻精品中文字幕| 又黄又爽又刺激的免费视频.| 色综合婷婷激情| 久久热精品热| 精品一区二区免费观看| 久久久久久久精品吃奶| 尤物成人国产欧美一区二区三区| 亚洲成av人片在线播放无| 久久欧美精品欧美久久欧美| 狠狠狠狠99中文字幕| 亚洲久久久久久中文字幕| 欧美激情国产日韩精品一区| 亚洲欧美日韩无卡精品| 久久天躁狠狠躁夜夜2o2o| 国产色婷婷99| 午夜福利高清视频| 亚洲av一区综合| 最新在线观看一区二区三区| 国产伦一二天堂av在线观看| 久久精品国产自在天天线| 国产一区二区亚洲精品在线观看| 一个人观看的视频www高清免费观看| 精品日产1卡2卡| 永久网站在线| 色吧在线观看| 日日干狠狠操夜夜爽| 真实男女啪啪啪动态图| 在线十欧美十亚洲十日本专区| 少妇被粗大猛烈的视频| 美女xxoo啪啪120秒动态图 | 国产精品综合久久久久久久免费| 在线免费观看的www视频| bbb黄色大片| av专区在线播放| 国产视频一区二区在线看| 亚洲内射少妇av| 成人性生交大片免费视频hd| 亚洲国产日韩欧美精品在线观看| 亚洲成人中文字幕在线播放| 97热精品久久久久久| 久久国产精品人妻蜜桃| 久久人人精品亚洲av| 色播亚洲综合网| 看免费av毛片| 五月伊人婷婷丁香| 日韩高清综合在线| 亚洲午夜理论影院| 欧美午夜高清在线| 在线看三级毛片| 亚洲经典国产精华液单 | 别揉我奶头 嗯啊视频| 全区人妻精品视频| 久久久久久国产a免费观看| 亚洲人成伊人成综合网2020| 亚洲精品一卡2卡三卡4卡5卡| 少妇的逼好多水| 欧美xxxx性猛交bbbb| 国产大屁股一区二区在线视频| 婷婷亚洲欧美| 国产高清视频在线观看网站| 国产淫片久久久久久久久 | 一个人看的www免费观看视频| 免费av毛片视频| 国产成人aa在线观看| 欧美日本视频| 九九久久精品国产亚洲av麻豆| 男人的好看免费观看在线视频| 欧美潮喷喷水| 午夜影院日韩av| 日本 欧美在线| 免费人成在线观看视频色| 日韩精品青青久久久久久| 亚洲av一区综合| 免费黄网站久久成人精品 | 国产探花极品一区二区| 国产乱人视频| 99在线人妻在线中文字幕| 欧美xxxx黑人xx丫x性爽| 又爽又黄a免费视频| 精品一区二区三区视频在线观看免费| 欧美精品国产亚洲| 日本撒尿小便嘘嘘汇集6| 欧美日韩综合久久久久久 | 99视频精品全部免费 在线| 国产精品日韩av在线免费观看| 黄色丝袜av网址大全| 精品人妻1区二区| 天天躁日日操中文字幕| 亚洲国产精品久久男人天堂| 三级国产精品欧美在线观看| 久久精品国产清高在天天线| 在线看三级毛片| 久久6这里有精品| 日韩欧美精品免费久久 | 日韩成人在线观看一区二区三区| 一区二区三区四区激情视频 | 亚洲第一欧美日韩一区二区三区| 国语自产精品视频在线第100页| 亚洲av免费在线观看| 日本 欧美在线| 国产午夜精品论理片| 亚洲综合色惰| 日本撒尿小便嘘嘘汇集6| 亚洲中文字幕一区二区三区有码在线看| 亚洲性夜色夜夜综合| 淫妇啪啪啪对白视频| 最新中文字幕久久久久| 国产淫片久久久久久久久 | 久久国产乱子免费精品| 桃红色精品国产亚洲av| 99国产精品一区二区三区| 国产精品三级大全| 亚洲av不卡在线观看| 亚洲成人久久性| 一个人观看的视频www高清免费观看| 最近最新中文字幕大全电影3| 女同久久另类99精品国产91| 深爱激情五月婷婷| 成人av在线播放网站| 亚洲av中文字字幕乱码综合| 少妇被粗大猛烈的视频| 精品欧美国产一区二区三| 人妻制服诱惑在线中文字幕| 国产国拍精品亚洲av在线观看| 国产精品久久久久久久久免 | 久久精品国产亚洲av香蕉五月| 精品无人区乱码1区二区| 在线天堂最新版资源| 女生性感内裤真人,穿戴方法视频| 国产精品永久免费网站| 精品人妻视频免费看| 国产精品嫩草影院av在线观看 | 亚洲男人的天堂狠狠| 国产三级在线视频| 最近在线观看免费完整版| h日本视频在线播放| 国产黄a三级三级三级人| 欧美高清性xxxxhd video| 国产乱人视频| 一边摸一边抽搐一进一小说| 成人高潮视频无遮挡免费网站| 亚洲av成人av| 中文字幕免费在线视频6| 一区二区三区高清视频在线| 久久天躁狠狠躁夜夜2o2o| 欧美+日韩+精品| 久久人人爽人人爽人人片va | 在线播放无遮挡| 深夜精品福利| 黄色丝袜av网址大全| 国产高清激情床上av| 国产精品久久久久久久久免 | 久久久久久大精品| 真实男女啪啪啪动态图| 亚洲精品亚洲一区二区| 97热精品久久久久久| 又黄又爽又免费观看的视频| 色播亚洲综合网| 在线a可以看的网站| 精品免费久久久久久久清纯| 午夜日韩欧美国产| 美女高潮喷水抽搐中文字幕| 天堂影院成人在线观看| 黄色女人牲交| 内地一区二区视频在线| 国产野战对白在线观看| 在线看三级毛片| 在线观看av片永久免费下载| 亚洲精品影视一区二区三区av| 好看av亚洲va欧美ⅴa在| 女生性感内裤真人,穿戴方法视频| 日本一本二区三区精品| 亚洲激情在线av| 在线播放国产精品三级| 午夜a级毛片| 老司机福利观看| av国产免费在线观看| 国产成人av教育| 国产精品综合久久久久久久免费| 欧美成狂野欧美在线观看| 亚洲人成网站在线播放欧美日韩| 欧美在线一区亚洲| 久久精品国产亚洲av涩爱 | 国产一区二区三区在线臀色熟女| 国产色爽女视频免费观看| 午夜影院日韩av| 欧美3d第一页| 不卡一级毛片| 久久天躁狠狠躁夜夜2o2o| 一区二区三区激情视频| 毛片一级片免费看久久久久 | 午夜激情欧美在线| 成年版毛片免费区| 成人无遮挡网站| 国产伦在线观看视频一区| 国产精品三级大全| 男女那种视频在线观看| 亚洲精品色激情综合| 国产 一区 欧美 日韩| av视频在线观看入口| 日日摸夜夜添夜夜添av毛片 | 国产一区二区在线观看日韩| 欧美成人性av电影在线观看| 欧美高清成人免费视频www| 国产私拍福利视频在线观看| x7x7x7水蜜桃| 1000部很黄的大片| 日本一二三区视频观看| 亚洲色图av天堂| 日本 欧美在线| 波野结衣二区三区在线| 99国产精品一区二区蜜桃av| 亚洲欧美日韩东京热| 成人午夜高清在线视频| 99国产精品一区二区三区| 天天躁日日操中文字幕| 老熟妇乱子伦视频在线观看| 69av精品久久久久久| 国产白丝娇喘喷水9色精品| 制服丝袜大香蕉在线| 天天一区二区日本电影三级| 久久久久久久精品吃奶| 成人国产综合亚洲| 成人无遮挡网站| 麻豆成人av在线观看| 18美女黄网站色大片免费观看| 亚洲专区国产一区二区| 性色avwww在线观看| 一个人看的www免费观看视频| 亚洲av第一区精品v没综合| 婷婷亚洲欧美| 午夜老司机福利剧场| av天堂中文字幕网| 性插视频无遮挡在线免费观看| 白带黄色成豆腐渣| 午夜老司机福利剧场|