MATLAB空間面板數(shù)據(jù)模型操作簡(jiǎn)介

word格式-可編輯-感謝下載支持MATLAB空間面板數(shù)據(jù)模型操作簡(jiǎn)介MATLAB安裝：在民主湖資源站上下載2009a2010a（MATLAB較大，占用內(nèi)存較大，安裝的話(huà)可能也要花費(fèi)一定的時(shí)間）一、數(shù)據(jù)布局 ：首先我們說(shuō)一下處理空間面板數(shù)據(jù)時(shí)，數(shù)據(jù)文件是怎么布局的，熟悉eviews的同學(xué)可能eviews（縱截面：一個(gè)時(shí)間序列1-91-9”“1-961-97”中，1是省份的代號(hào)，94,95,96,97表示年份，eviews是將每個(gè)省份的數(shù)據(jù)放在一起，。與eviewsMATLAB處理空間面板數(shù)據(jù)時(shí)面板數(shù)據(jù)的布局是（在excel中說(shuō)明個(gè)橫截面上的數(shù)據(jù)。如圖：word格式-可編輯-感謝下載支持這里需要說(shuō)明的是，MATLAB中省份的序號(hào)需要與空間權(quán)重矩陣中省份一一對(duì)應(yīng)用《中國(guó)統(tǒng)計(jì)年鑒》分地區(qū)數(shù)據(jù)中省份的排列順序（二階空間權(quán)重矩陣我會(huì)在附件中給出。二、數(shù)據(jù)的輸入：MATLAB與excel鏈接：在excel中點(diǎn)擊“工具→加載宏→瀏覽MATLAB的安裝目錄，一點(diǎn)擊excllink.xla即可完成excel與MATLAB的鏈接excel中的數(shù)據(jù)就可以直接導(dǎo)入中形成的數(shù)據(jù)文件。操作完成后excel的加載宏界面如圖：SpreadsheetLinkEX3.0.3forusewithMATLA即表示我們希望excel與MATLAB點(diǎn)擊確定后，excel窗口如圖：word格式-可編輯-感謝下載支持紅色框內(nèi)選項(xiàng)為：excel實(shí)現(xiàn)了與excel時(shí)，excel會(huì)自動(dòng)啟動(dòng)MATLAB，所以excel的打開(kāi)速度很慢。PutmatrixexcelMATLAB中，我們經(jīng)常用的就是此項(xiàng)，選中需要導(dǎo)入到MATLAB中的數(shù)據(jù)就是一個(gè)矩陣點(diǎn)擊putmatrix即可的數(shù)據(jù)就是我們需要導(dǎo)入MATLAB中的數(shù)據(jù)。word格式-可編輯-感謝下載支持選中數(shù)據(jù)矩陣，點(diǎn)擊putmatrix，即我們選中并輸入到中的數(shù)據(jù)矩陣在該矩陣的某一列（后面詳解pword格式-可編輯-感謝下載支持三、MATLAB窗口此時(shí)原來(lái)已經(jīng)打開(kāi)的MATLAB窗口如圖：

標(biāo)題欄菜單欄工具欄工作空間控制窗口歷史記錄這里我們簡(jiǎn)單介紹一下MATLAB的窗口，與很多我們見(jiàn)過(guò)的很多軟件，如excel，word，spss一樣，MATLAB一樣有標(biāo)題欄、菜單欄、工具欄。我們主要講MATLAB的工作窗口，如圖，左側(cè)是“控制窗口個(gè)文件，查找函數(shù)，建立模型等等。

word格式-可編輯-感謝下載支持右側(cè)上面是“工作空間（產(chǎn)生的數(shù)據(jù)文件等p”的一個(gè)數(shù)據(jù)矩陣，就存放在這里，雙擊“p進(jìn)行原來(lái)某項(xiàng)操作時(shí)，其實(shí)可以點(diǎn)擊歷史窗口中相應(yīng)的記錄就可。currentdirector前目錄”就會(huì)以浮動(dòng)的形式出現(xiàn)當(dāng)前目錄，或方便我們查找我們需要的目錄。這是我習(xí)慣的一種桌面布局Deskto（桌面desktop→desktoplayou（桌面布局，會(huì)有幾種布局形式，大家自己去看吧。四、與MATLAB相關(guān)的重要操作1、數(shù)據(jù)文件存儲(chǔ)之前我們已經(jīng)將模型需要的數(shù)據(jù)作為一個(gè)變量導(dǎo)入并存儲(chǔ)在MATLAB的“工作空間”中，由于我們是（w這時(shí)我們需要將包含空間權(quán)重矩陣“w1”和所有數(shù)據(jù)的“pc”的數(shù)據(jù)文件進(jìn)行存儲(chǔ)，方便以后調(diào)用“Fil→SaveWorkspaceA輸入文件名.mat的MATLAB能識(shí)別的數(shù)據(jù)文件。2、設(shè)置路徑MATLAB（我打包給你們當(dāng)然我們也需要將MATLABMATLAB通過(guò)調(diào)用軟件包中已有我們需要告訴MATLAB它的權(quán)限在哪里，即它可以在哪里去查找這些程序、數(shù)據(jù)等默認(rèn)狀況下，MATLAB只承認(rèn)安裝目錄下文件夾中已有的程序我們需要為MATLAB設(shè)置一條路徑Fil→SetPatsetpat”窗口，如圖：空間計(jì)量工具箱紅框下按鈕“AddwithSubfolders”表示在添加路徑時(shí)，在文件夾下的子目錄也被MATLAB承認(rèn)，D盤(pán)下的文件夾及其子文件夾添加到了MATLAB等都能被MATLAB識(shí)別并使用。

word格式-可編輯-感謝下載支持五、MATLAB程序——空間計(jì)量工具箱一般來(lái)說(shuō)我們不會(huì)在的控制窗口一個(gè)一個(gè)的輸入命令而是在MATLAB的 “Editor”（編輯器）已經(jīng)編輯好的程序，進(jìn)行一定的組合和修改來(lái)進(jìn)行自己的研究。官方網(wǎng)上下載的最新版的空間計(jì)量工具箱的名稱(chēng)jplvD盤(pán)下的MATLABcurrentdirector（當(dāng)前目錄）即顯示該目錄下spatia完整程序pane。幫助量模型我們可以先看一下以“demo”開(kāi)頭的M文件（MATLAB程序即演示，如雙擊打開(kāi)“demoLMsarsem_paneedito運(yùn)行結(jié)果就會(huì)顯示在的“控制窗口”中。word格式-可編輯-感謝下載支持六、空間面板數(shù)據(jù)模型解讀（一）空間面板數(shù)據(jù)模型建立的一般步驟是：1MorasILM檢驗(yàn)（這是我目前知道的兩種檢驗(yàn)方法，對(duì)于其他的一些檢驗(yàn)如也不是很清楚，我也沒(méi)找到相關(guān)的程序，只是看到有些人做空間面板數(shù)據(jù)模型時(shí)也給出了這幾種檢驗(yàn)的結(jié)果。LMmoran檢驗(yàn)的M空間相關(guān)性檢驗(yàn)可以直接使用這個(gè)程序，當(dāng)然需要修改某些東西（后面有詳解。2SARSEM以及在固定或隨機(jī)效應(yīng)模型中選擇（地區(qū)固定、時(shí)期固定、雙固定、還是混合面板。LM檢驗(yàn)的結(jié)果可以為SAR與SEM（查看相關(guān)文獻(xiàn)SAR模型直接調(diào)用“sar_panel_FE”程序，而固定效應(yīng)SEM直接調(diào)用“sem_panel_FE”程序（為程序文件存在于某個(gè)子目錄下。對(duì)于固定效應(yīng)（地區(qū)固定，時(shí)間固定，雙向固定，混合）的選擇，其實(shí)只需要在引用上述程序時(shí)，對(duì)某些指標(biāo)進(jìn)行設(shè)置。如下圖，是截取的一段面板SAR文字去理解這段程序在做什么。紅框中是設(shè)定模型是那種固定效應(yīng)：info.model=0：表示此模型為混合模型，即沒(méi)有固定效應(yīng)；info.model=1表示雙向固定效應(yīng)模型。關(guān)于如何使用已有的完整程序來(lái)做自己的研究將在后面詳解。（二）學(xué)會(huì)使用“幫助”自學(xué)空間面板數(shù)據(jù)關(guān)于“paneldemoMword格式-可編輯-感謝下載支持MATLAB命令名”進(jìn)行查找，或者直接使用MATLAB提供的幫助，看該命令（上上圖給出了幫助的位置MATLAB用于空間面板計(jì)量”下的每一個(gè)M文件到底是干什么用的。（三）常見(jiàn)的命令：1、關(guān)于%：%后面可以輸入中文，對(duì)已寫(xiě)程序進(jìn)行描述和解釋?zhuān)?后的文字是不被運(yùn)行的。2、wk1read.wk1（wk1數(shù)據(jù)文件是MATLAB能識(shí)別的數(shù)據(jù)文件，但具EXCEL鏈接導(dǎo)入到MATLAB.mat，該文件的讀取，直接使用命令loadmat文件，直接用loadpc.mat）3數(shù)據(jù)進(jìn)行描述（也叫定義空間面板數(shù)據(jù)，同時(shí)指定被解釋變量y和解釋變量x：如下圖：每一行分別表示：該面板數(shù)據(jù)的時(shí)期數(shù)為該面板數(shù)據(jù)有30個(gè)地區(qū)N=3，將空間權(quán)重矩陣標(biāo)準(zhǔn)化W=normw(w1，將名為consume（以矩陣形式出現(xiàn)在MATLABA中）的變量的第3列數(shù)據(jù)定義為被解釋變量y,

word格式-可編輯-感謝下載支持將名為consume45、、、91012列數(shù)據(jù)定義為解釋變量矩陣x，定義一個(gè)有N*T1列的全1xconstanones即為全1矩陣說(shuō)明解釋變量矩陣x的大?。河衝obs行，K（size為描述矩陣的大小。這里說(shuō)明一下矩陣的引用：（[6）括號(hào)[][1:7]即可；如果幾行將是間斷的，用逗號(hào)隔開(kāi)；如果是全部的行或列，可省略數(shù)值和中括號(hào)，直接打冒號(hào)“4是讀取的文件名，面板數(shù)據(jù)的定義和解釋變量x被解釋變量y的指定，以及在輸出時(shí)變量的名稱(chēng)，前面的都講了，現(xiàn)在說(shuō)說(shuō)變量的名稱(chēng)。下面是一段最小二乘法回歸以及LM檢驗(yàn)，moran檢驗(yàn)的程序：“results=ols(y,[xconstantx]);”表示以y作為被解釋變量，x作為解釋變量的含有常數(shù)項(xiàng)的最小二乘法回歸（xconstant即表示含有常數(shù)項(xiàng)），該命令運(yùn)行結(jié)果全部存儲(chǔ)在results中;“vnames=strvcat('gy','intercept','logy0','logpc','gpc','emp','edu','inv','fir');”即是對(duì)最終顯示在控制窗口中的變量進(jìn)行命名，將此名字賦值