第七講空間計(jì)量經(jīng)濟(jì)學(xué)模型的matlab估計(jì).

1、空間計(jì)量經(jīng)濟(jì)學(xué)基本模型的matlab估計(jì)空間滯后模型sar（）函數(shù)功能估計(jì)空間滯后模型（空間自回歸-回歸模型）y=WyxN（0廣2In）中的未知參數(shù)p、B和d2。使用方法res=sar（y,x,W,info）*res:存儲結(jié)果的變量；y:被解釋變量；x:解釋變量；w：空間權(quán)重矩陣；info：結(jié)構(gòu)化參數(shù)，具體可使用helpsar語句查看ii注意事項(xiàng)1）WW為權(quán)重矩陣，因?yàn)槭窍∈杈仃嚕紨?shù)據(jù)通常以nX3的數(shù)組形式存儲，需要用sparse函數(shù)轉(zhuǎn)換為矩陣形式。*2）ydev（不再需要）sar函數(shù)求解的標(biāo)準(zhǔn)模型可以包含常數(shù)項(xiàng)，被解釋變量（因變量）y，不再需要轉(zhuǎn)換為離差形式（ydev）。*3）x需要注意

2、x的生成方式，應(yīng)將常數(shù)項(xiàng)包括在內(nèi)。*4）infoinfo為結(jié)構(gòu)化參數(shù)，事前賦值；通常調(diào)整info.1flag（標(biāo)準(zhǔn)n?1000）、nfo.rmin和info.rmax。*5）vnames在輸出結(jié)果中說明被解釋變量。使用方法：vnames=strvcat（variablenamel,variablename2）;*6）Asymptotict-stat（漸進(jìn)t統(tǒng)計(jì)量）rho的檢驗(yàn)：漸進(jìn)t分布，估計(jì)值的顯著性使用相應(yīng)的Z概率表示。應(yīng)用實(shí)例估計(jì)地區(qū)投票率受周邊地區(qū)投票率的影響程度案例素材1997年，Pace等人研究了美國3107個(gè)縣的選舉投票率影響因素，運(yùn)用的是美國1980年大選的公開投票數(shù)據(jù)，形成了

3、一個(gè)包含3107個(gè)樣本數(shù)據(jù)的截面數(shù)據(jù)集elect.dat可以通過matlab軟件打開elect.dat查看，并打開elect.txt查看各列數(shù)據(jù)的含義。計(jì)量模型認(rèn)為各縣的投票率受到相鄰地區(qū)投票率的影響，同時(shí)，還受到選民教育水平、選民住房情況、選民收入水平的影響，據(jù)此得到如下計(jì)量模型：y=（3o+pWy+x（3+eLN（0,異In）轉(zhuǎn)換為：y=pWy+1x3o3eN（0,異In）程序語句1）近似估計(jì)缺省設(shè)置：info.lflag=1注意取對數(shù)值，得到y(tǒng),x。2）精確估計(jì)infoflag=0運(yùn)行結(jié)果xy2cont（）函數(shù)功能：使用地區(qū)x坐標(biāo)和y坐標(biāo)，生成空間鄰接矩陣。使用方法：W1W2W3=xy2

4、cont（x,y）其中，W2是行標(biāo)準(zhǔn)化后的空間鄰接矩陣。一個(gè)例子：使用anselin數(shù)據(jù)，生成w，并與wmat比較其差異。141、空間誤差模型sem()函數(shù)功能估計(jì)空間誤差模型廠Xuu=WuNg2In)中的未知參數(shù)伏入和*。使用方法res=sem，x，W，info)*res:存儲結(jié)果的變量；y:被解釋變量；x:解釋變量；w:空間權(quán)重矩陣；info:結(jié)構(gòu)化參數(shù)，具體可使用helpsem語句查看注意事項(xiàng)1) xx應(yīng)將常數(shù)項(xiàng)包括在內(nèi)。*2) infoinfo為結(jié)構(gòu)化參數(shù)，事前賦值；通常調(diào)整info.1flag(標(biāo)準(zhǔn)n?1000)、nfo.rmin和info.rmax。*3) vnames在輸出結(jié)果中

5、說明被解釋變量。使用方法：vnames=strvcat(variablenamel,variablename2);*應(yīng)用實(shí)例估計(jì)地區(qū)投票率受周邊地區(qū)投票率的影響程度案例素材1997年，Pace等人研究了美國3107個(gè)縣的選舉投票率影響因素，運(yùn)用的是美國1980年大選的公開投票數(shù)據(jù)，形成了一個(gè)包含3107個(gè)樣本數(shù)據(jù)的截面數(shù)據(jù)集elect.dat可以通過matlab軟件打開elect.dat查看，并打開elect.txt查看各列數(shù)據(jù)的含義。計(jì)量模型認(rèn)為各縣的投票率受到相鄰地區(qū)投票率的影響，同時(shí)，還受到選民教育水平、選民住房情況、選民收入水平的影響，據(jù)此得到如下計(jì)量模型：y=po+xB+uu=AWu

6、+&N(0,(72In)轉(zhuǎn)換為：y=1xBoBuu=AWu+&N(0,72In)程序語句1) 近似估計(jì)缺省設(shè)置：info.lflag=1注意取對數(shù)值，得到y(tǒng),x。2) 精確估計(jì)infoflag=0運(yùn)行結(jié)果誤差項(xiàng)空間依賴性的檢驗(yàn)*1)MoranI統(tǒng)計(jì)量檢驗(yàn)使用方法:res=moran(y,x,W);prt(res);*2) 似然比檢驗(yàn)(Iratios)使用方法：res=lratios(y,x,W);prt(res);*3) 拉格朗日乘子(LM)檢驗(yàn)使用方法：res=lmerror(y,x,W);prt(res);*4) 沃德(Walds)檢驗(yàn)使用方法：res=walds(y,x,W);prt(r

7、es);*5) 基于sar殘差的檢驗(yàn)使用方法：res=lmsar(y,x,W1,W2);prt(res);三、空間杜賓模型sdm()函數(shù)功能估計(jì)空間杜賓模型yWyX1000，只能使用近似估計(jì)。運(yùn)行結(jié)果四、廣義空間模型sac()函數(shù)功能估計(jì)廣義空間模型y=pWiy+xB+uU=2W2U+N（0,異In）中的未知參數(shù)p、B、入和o2使用方法res=sac,X,Wi,W2,info)高階鄰接矩陣的生成slag()Wp=slag(W,p)注意：不包括低階鄰接。應(yīng)用實(shí)例1) 估計(jì)地區(qū)犯罪率受周邊地區(qū)犯罪率的影響程度2) 估計(jì)地區(qū)投票率受周邊地區(qū)投票率的影響程度練習(xí)作業(yè)1）估計(jì)地區(qū)犯罪率受周邊地區(qū)犯罪率以及其他因素的影響程度數(shù)據(jù)：Anselin（1980）。方法：采用精確估計(jì)，分別采用空間滯后、空間誤差、空間杜賓和廣義空間模型估計(jì)。具體任務(wù)：分別按照不同的模型，（1）編寫運(yùn)算語句；（2）得到運(yùn)算結(jié)果，并保存在作業(yè)里；（3）分析不同方法得到的結(jié)果差異。2）估計(jì)地區(qū)投票率受周邊地區(qū)投票率