在空間OLS回歸分析中如果得到LM檢驗(yàn)并且判斷得到應(yīng)該使用空間誤差SEM模型時(shí),接著本文檔介紹空間誤差SEM模型。首先空間誤差SEM模型的數(shù)學(xué)模式公式如下:
y = βk * x + u , u = λ * Wu + µ(µ為擾動(dòng)項(xiàng)),Wu為誤差(擾動(dòng)項(xiàng))空間滯后變量,λ為其回歸系數(shù)值
空間誤差模型時(shí),其將誤差項(xiàng)納入模型中,其實(shí)質(zhì)是將誤差項(xiàng)自相關(guān)納入考慮中(空間滯后模型是將因變量空間滯后變量納入模型)。
空間誤差模型SEM案例
- 1、背景
當(dāng)前有一份空間數(shù)據(jù),其為美國(guó)哥倫布市49個(gè)社區(qū)的相關(guān)數(shù)據(jù),包括犯罪率(crime)、房?jī)r(jià)(hoval)和家庭收入(income),當(dāng)前希望研究房?jī)r(jià)和家庭收入對(duì)于犯罪率的影響關(guān)系,并且在研究這一影響關(guān)系時(shí),考慮空間性,并且使用空間誤差SEM模型進(jìn)行分析(具體是否應(yīng)該使用空間誤差模型,通常以空間OLS回歸的LM檢驗(yàn)進(jìn)行分析判斷決定)。部分?jǐn)?shù)據(jù)如下圖所示:
上面展示的是‘分析數(shù)據(jù)’,共有49個(gè)社區(qū),該49個(gè)社區(qū)對(duì)應(yīng)的‘空間權(quán)重矩陣’如下圖所示:
圖中數(shù)字1表示兩個(gè)空間點(diǎn)(社區(qū))之間相鄰,數(shù)字0表示兩個(gè)社區(qū)不相鄰?臻g權(quán)重矩陣數(shù)據(jù)可點(diǎn)擊此處下載。
- 2、理論
空間誤差模型SEM的自變量包括誤差項(xiàng),其意義為X對(duì)于Y無(wú)法解決的部分,其具有空間效應(yīng)關(guān)系。其數(shù)學(xué)模式公式如下:
y = βk * x + u , u = λ * Wu + µ(µ為擾動(dòng)項(xiàng)),Wu為誤差(擾動(dòng)項(xiàng))空間滯后變量,λ為其回歸系數(shù)值
- 3、操作
本例子操作如下:
下拉選擇‘空間權(quán)重矩陣’文檔即spatialweight這份數(shù)據(jù),默認(rèn)對(duì)空間權(quán)重矩陣行標(biāo)準(zhǔn)化處理,需要注意的是,空間權(quán)重矩陣通常需要進(jìn)行行標(biāo)準(zhǔn)化處理。
另需要提示的是,在使用空間計(jì)量相關(guān)的方法時(shí),其均需要‘空間權(quán)重矩陣’和‘分析數(shù)據(jù)’兩份數(shù)據(jù),并且均需要單獨(dú)上傳到SPSSAU中,并且對(duì)‘分析數(shù)據(jù)’進(jìn)行分析時(shí),下拉選擇對(duì)應(yīng)的‘空間權(quán)重矩陣’,操作上分為以下3個(gè)步驟。
第1:上傳‘空間權(quán)重矩陣’文檔
此處需要注意:上傳的數(shù)據(jù)需要為n*n階格式,而且第1行為空間點(diǎn)的名稱(比如31省市的名稱)。類似下圖格式:
第2:上傳‘分析數(shù)據(jù)’文檔
此處需要注意:比如31省市數(shù)據(jù),‘空間權(quán)重矩陣’有著該31個(gè)空間點(diǎn)的順序比如北京-》天津-》河北-》山西-》…,那么‘分析數(shù)據(jù)’的31行數(shù)據(jù)也需要按此順序才可以。
第3:針對(duì)‘分析數(shù)據(jù)’進(jìn)行分析,并且選擇‘空間權(quán)重矩陣’文檔
此處需要注意:進(jìn)行某空間研究方法時(shí)需要下拉選擇‘空間權(quán)重矩陣’,選擇后,SPSSAU會(huì)自動(dòng)判斷其是否為‘空間權(quán)重矩陣’格式,包括是否為n*n階結(jié)構(gòu),是否具有對(duì)稱性等。如果不是則會(huì)進(jìn)行信息提示,請(qǐng)勿必注意空間權(quán)重矩陣數(shù)據(jù)格式。
- 4、SPSSAU輸出結(jié)果
SPSSAU共輸出6個(gè)表格,分別是模型基本參數(shù)等、空間誤差SEM模型分析結(jié)果、空間誤差SEM模型相關(guān)檢驗(yàn)匯總、信息準(zhǔn)則指標(biāo)結(jié)果、空間效應(yīng)分析和空間誤差SEM模型分析結(jié)果-簡(jiǎn)化格式表格,如下所述。
表格
說(shuō)明
模型基本參數(shù)等
輸出模型的基礎(chǔ)參數(shù)值信息等
空間誤差SEM模型分析結(jié)果
輸出模型的分析結(jié)果,包括回歸系數(shù)和顯著性檢驗(yàn)結(jié)果等
空間誤差SEM模型相關(guān)檢驗(yàn)匯總
輸出相關(guān)的檢驗(yàn)比如異方差檢驗(yàn)等
信息準(zhǔn)則指標(biāo)結(jié)果
如果是極大似然ML法時(shí)則會(huì)輸出信息準(zhǔn)則指標(biāo)等
空間效應(yīng)分析
輸出空間效應(yīng)分析表格
空間誤差SEM模型分析結(jié)果-簡(jiǎn)化格式
輸出模型結(jié)果的簡(jiǎn)化表格格式
- 5、文字分析
上表格模型的基本參數(shù)信息,包括具體的空間計(jì)量模型名稱,是否使用穩(wěn)健標(biāo)準(zhǔn)誤差,空間權(quán)重矩陣名稱及是否對(duì)其進(jìn)行標(biāo)準(zhǔn)化處理等,模型估計(jì)方法等,表格中僅展示模型的參數(shù)信息等無(wú)特別分析意義。需要注意的是,當(dāng)前默認(rèn)使用ML極大似然法進(jìn)行估計(jì),但當(dāng)選中Robust穩(wěn)健標(biāo)準(zhǔn)誤法時(shí),則使用GMM估計(jì),GMM估計(jì)法時(shí)不會(huì)輸出llf指標(biāo)等,即其會(huì)影響到后續(xù)輸出信息準(zhǔn)則指標(biāo)表格。
上表格展示空間誤差SEM模型回歸結(jié)果,其數(shù)學(xué)模型為y = β * x + u, u = λ * Wu + μ (其中β表示X的回歸系數(shù),Wu表示u的空間滯后變量,λ表示W(wǎng)u的回歸系數(shù),u和μ為誤差擾動(dòng)項(xiàng)),結(jié)合當(dāng)前數(shù)據(jù),其公式為:crime = 66.756-0.223*hoval-1.616*income+0.361*殘差空間滯后變量。
具體針對(duì)各項(xiàng)的影響關(guān)系來(lái)看:hoval的回歸系數(shù)值為-0.223,并且呈現(xiàn)出0.05水平顯著性(p =0.017<0.05),意味著hoval會(huì)對(duì)crime產(chǎn)生顯著的負(fù)向影響關(guān)系,即說(shuō)明房?jī)r(jià)會(huì)負(fù)向影響犯罪率,房?jī)r(jià)越高犯罪率越低。income的回歸系數(shù)值為-1.616,并且呈現(xiàn)出0.01水平顯著性(p =0.000<0.01),意味著income會(huì)對(duì)crime產(chǎn)生顯著的負(fù)向影響關(guān)系,家庭收入越高犯罪率越低。誤差項(xiàng)空間滯后變量的回歸系數(shù)lambda值為0.361,并且呈現(xiàn)出0.05水平顯著性(p =0.034<0.05),意味著誤差項(xiàng)有著空間項(xiàng),使用空間誤差模型較為適合。
上表格展示異方差White檢驗(yàn)和JB檢驗(yàn)等,分別用于異方差和正態(tài)性檢驗(yàn),空間計(jì)量模型時(shí)對(duì)于空間作用的關(guān)注力度明顯最高,對(duì)于異方差和正態(tài)性關(guān)注度相對(duì)較低,從上表格可以看到,有著一定的異方差問(wèn)題,因而進(jìn)一步分析時(shí)可考慮選擇穩(wěn)健標(biāo)準(zhǔn)誤法,最終使用GMM估計(jì)法得到更科學(xué)的分析結(jié)果。
上表格展示信息準(zhǔn)則結(jié)果表格,包括llf值和另外兩個(gè)值即AIC值和Schwarz準(zhǔn)則值,llf值通常越大越好,但是AIC值和Schwarz準(zhǔn)則值均是越小越好,如果希望對(duì)比模型優(yōu)劣,可考慮使用上述三個(gè)指標(biāo),但需要注意的是,極大似然法估計(jì)ML法時(shí)才會(huì)輸出上述指標(biāo),如果是比如GMM估計(jì)則沒(méi)有輸出上述指標(biāo)。
上表格展示空間效應(yīng)分析結(jié)果,直接效應(yīng)ADI反映自變量X對(duì)于自身區(qū)域Y的平均影響效應(yīng)情況,間接(溢出)效應(yīng)AII反應(yīng)自變量X對(duì)其它區(qū)域Y的平均影響效應(yīng)情況,總效應(yīng)ATI=直接效應(yīng)ADI+間接(溢出)效應(yīng)AII。在空間誤差SEM時(shí),其僅考慮誤差項(xiàng)空間滯后變量放入模型中,單獨(dú)對(duì)于自變量其不會(huì)有間接溢出效應(yīng),因而AII值全部均為0,總效應(yīng)ATI值即直接效應(yīng)ADI值。其空間效應(yīng)的計(jì)算公式如下:
上表格展示模型的簡(jiǎn)化表格格式,不再重復(fù)分析。
- 6、剖析
- 涉及以下幾個(gè)關(guān)鍵點(diǎn),分別如下:
通常是使用空間OLS回歸分析得到LM檢驗(yàn),并且結(jié)合LM檢驗(yàn)判斷后決定是否使用空間滯后模型。
疑難解惑
- 空間誤差SEM模型時(shí)誤差項(xiàng)空間滯后變量的意義?
誤差項(xiàng)空間滯后變量,是指X對(duì)于Y無(wú)法解決部分的空間滯后項(xiàng),即Y無(wú)法解決的其它部分,它具有空間效應(yīng)關(guān)系,通常只需要考慮是否需要使用空間誤差SEM模型,并且觀察誤差項(xiàng)空間滯后項(xiàng)是否具有顯著性即可,如果呈現(xiàn)出顯著性則意味著當(dāng)前模型應(yīng)該是適合的,因?yàn)榭紤]該項(xiàng)時(shí)該項(xiàng)呈現(xiàn)出顯著性。