MATLAB空间面板数据模型操作介绍

MATLAB 空间面板数据模型操作简介空间面板数据模型操作简介 MATLAB 安装 安装 在民主湖资源站上下载 MATLAB 2009a 或者 2010a 按照其中的安装说明 安装 MATLAB MATLAB 较大 占用内存较大 安装的话可能也要花费一定的时间 一 数据布局一 数据布局 首先我们说一下 MATLAB 处理空间面板数据时 数据文件是怎么布局的 熟悉 eviews 的同学 可能知道 eviews 中面板数据布局是 一个省份所有年份的数据作为一个单元 纵截面 一个时间 序列 然后再排放另一个省份所有年份的数据 依次将所有省份的数据排放完 如下图 红框中 1 94 1 95 1 96 1 97 中 1 是省份的代号 94 95 96 97 表示年份 eviews 是将每个省 份的数据放在一起 再将所有省份堆放在一起 与 eviews 不同 MATLAB 处理空间面板数据时 面板数据的布局是 在 excel 中说明 先 排放一个横截面上的数据 即某年所有省份的数据 再将不同年份的横截面按时间顺序堆放在一 起 如图 这里需要说明的是 MATLAB 中省份的序号需要与空间权重矩阵中省份一一对应 我们一般 就采用 中国统计年鉴 分地区数据中省份的排列顺序 二阶空间权重矩阵我会在附件中给出 二 数据的输入 二 数据的输入 MATLAB 与 excel 链接 在 excel 中点击 工具 加载宏 浏览 找到 MATLAB 的安装目 录 一般来说 如果安装时没有修改安装路径 此安装目录为 C Programfiles MATLAB R2009a toolbox exlink 点击 excllink xla 即可完成 excel 与 MATLAB 的链 接 这样的话 excel 中的数据就可以直接导入 MATLAB 中形成 MATLAB 的数据文件 操作完成后 excel 的加载宏界面如图 选中 Spreadsheet Link EX3 0 3 for use with MATLAB 即表示我们希望 excel 与 MATLAB 实现 链接 点击确定后 excel 窗口如图 红色框内选项为 startMATLAB 表示打开 MATLAB 一般来说 如果 excel 实现了与 MATLAB 的链接 当我们 打开 excel 时 excel 会自动启动 MATLAB 所以 excel 的打开速度很慢 Putmatrix 表示将 excel 中的数据输入 MATLAB 中 我们经常用的就是此项 选中需要导入 到 MATLAB 中的数据 就是一个矩阵 点击 putmatrix 即可 如 下表是一个数据表格 那么红 色框内的数据就是我们需要导入 MATLAB 中的数据 选中数据矩阵 点击 putmatrix 出现如下窗口 我们需要在弹出的窗口中输入变量名称 即我们选 中并输入到 MATLAB 中的数据矩阵在 MATLAB 中是以一个变量的形式存在 我们之后在使用时 需要调用该矩阵的某一列 后面详解 如我们在此输入变量名为 pc 点击确定 三 三 MATLAB 窗口窗口 此时原来已经打开的 MATLAB 窗口如图 标题栏 菜单栏 工具栏 控制窗口 工作空间 历史记录 这里我们简单介绍一下 MATLAB 的窗口 与很多我们见过的很多软件 如 excel word spss 一样 MATLAB 一样有标题栏 菜单栏 工具栏 我们主要讲 MATLAB 的工作窗口 如图 左侧是 控制窗口 在此可输入很多命令 如打 开某一个文件 查找函数 建立模型等等 右侧上面是 工作空间 一般是存储我们需要的数据文件 包括导入的原始数据 模型模拟 过程中产生的数据文件等 如我们之前导入了变量名为 pc 的一个数据矩阵 就存放在这里 双击 pc 就能查看该数据文件 右侧下方是 历史记录 你在控制窗口进行的任何一项操作都会在这里有显示 所以当你在 此需要进行原来某项操作时 其实可以点击历史窗口中相应的记录就可 控制窗口左侧隐藏了 current directory 这个是 当前目录 即我们现在所在的位置 点击 当前目录 就会以浮动的形式出现当前目录 或方便我们查找我们需要的目录 这是我习惯的一种桌面布局 你也可以根据需要调整 在菜单栏有一项 Desktop 桌面 点击 desktop desktop layout 桌面布局 会有几种布局形式 大家自己去看吧 四 与四 与 MATLAB 相关的重要操作相关的重要操作 1 数据文件存储 之前我们已经将模型需要的数据作为一个变量导入并存储在 MATLAB 的 工作空间 中 由 于我们是处理空间面板数据 所以我们还需要导入空间权重矩阵 如命名为 w1 方法一样 这时我们需要将包含空间权重矩阵 w1 和所有数据的 pc 的数据文件进行存储 方便以后调用 数据存储方法 File Save Workspace As 找到存储的地址 输入文件名 就形成了后缀名为 mat 的 MATLAB 能识别的数据文件 2 设置路径 MATLAB 并没有自带空间计量软件包 这个需要在网上下载 我打包给你们 当然我们也需 要将操作 MATLAB 过程中需要的数据文件 程序等保存在某个文件夹中 MATLAB 通过调用软件 包中已有的程序命令 以及自己保存的数据文件 来实现模型的模拟 此时我们需要告诉 MATLAB 它的权限在哪里 即它可以在哪里去查找这些程序 数据等 默认状况下 MATLAB 只 承认安装目录下文件夹中已有的程序 所以我们需要为 MATLAB 设置一条路径 让它能够使用 这条路径下包含的程序 数据 设置路径的方法 点击 File Set Path 出现 set path 窗口 如图 红框下按钮 Add with Subfolders 表示在添加路径时 在文件夹下的子目录也被 MATLAB 承 认 否则该文件夹下的子文件夹中的程序 数据将不能被 MATLAB 识别 而不能使用 点击此按 钮 选择文件夹位置 然后点击该窗口的 Save 即可实现路径的设置 如图 我将 D 盘下的 MATLAB 文件夹及其子文件夹添加到了 MATLAB 能识别的路径中 那么此文件夹中的所有与 MATLAB 相关的程序 数据等都能被 MATLAB 识别并使用 五 五 MATLAB 程序程序 空间计量工具箱空间计量工具箱 一般来说 我们不会在 MATLAB 的控制窗口一个一个的输入命令 而是在 MATLAB 的 Editor 编辑器 中将一个完整的程序编辑完毕然后进行存储 方便以后使用 当然我们更 一般的是使用别人已经编辑好的程序 进行一定的组合和修改来进行自己的研究 这里先讲 MATLAB 的空间计量工具箱 在 MATLAB 官方网上下载的最新版的空间计量工具箱 的名称为 jplv7 如上图 我在设置路径时 该工具箱已经包含在 D 盘下的 MATLAB 文件 夹下 这时 我们可以从红框所标识的按钮找到该工具箱 如图 左侧 current directory 当前目 录 即显示该目录下所有的文件夹 其中我们要使用的是 spatial 该文件夹下包含我们需要 的关于空间面数据模型相关的程序 而我们主要使用的完整程序完整程序是该文件夹下的 panel 空间空间 面板面板 空间计量空间计量 工具工具箱箱 我们可以先看一下以 demo 开头的 M 文件 MATLAB 程序 demo 即演示 如双击打开 demoLMsarsem panel 就会弹出该程序文件的编辑器 editor 点击红圈下的运行按钮 整个 程序的运行结果就会显示在 MATLAB 的 控制窗口 中 六 空间面板数据模型解读六 空间面板数据模型解读 空间计空间计 量模型量模型 帮助帮助 一 空间面板数据模型建立的一般步骤是 一 空间面板数据模型建立的一般步骤是 1 空间相关性检验 Moran s I 检验 LM 检验 这是我目前知道的两种检验方法 对于其他 的一些检验如 LR Walds 等本来就是用于截面空间数据空间相关性的检验 如何推广到空间面板 数据上来 我也不是很清楚 我也没找到相关的程序 只是看到有些人做空间面板数据模型时也给 出了这几种检验的结果 打包文件中有一个 作业 文件夹 里面有单独的 LM 与 moran 检验的 M 文件 对自己的数 据进行空间相关性检验可以直接使用这个程序 当然需要修改某些东西 后面有详解 2 空间模型的选择 即选择用 SAR 模型 还是 SEM 模型 同时还涉及到固定效应和随机效 应选择 以及在固定或随机效应模型中选择 地区固定 时期固定 双固定 还是混合面板 LM 检验的结果可以为 SAR 与 SEM 的选择提供依据 查看相关文献 面板数据的固定效应 SAR 模型直接调用 sar panel FE 程序 而固定效应 SEM 直接调用 sem panel FE 程序 这些 程序已经作为程序文件存在于某个子目录下 对于固定效应 地区固定 时间固定 双向固定 混合 的选择 其实只需要在引用上述程序 时 对某些指标进行设置 如下图 是截取的一段面板 SAR 模型的程序 其实我们可以从百分号 后的文字去理解这段程序在做什么 红框中是设定模型是那种固定效应 info model 0 表示此模型为混合模型 即没有固定效应 info model 1 表示此模型为地区固定效应模型 info model 2 表示时间固定效应模型 info model 3 表示双向固定效应模型 关于如何使用已有的完整程序来做自己的研究将在后面详解 二 学会使用 二 学会使用 帮助帮助 自学空间面板数据自学空间面板数据 关于 panel 文件夹下带有 demo 字样的 M 文件 一般都是一个完整的演示文件 可以运行 大家可以读一下 一般等号后面的语句就是对等号左边的变量进行赋值之类的 对于任何命令 都 可以在 MATLAB 的控制窗口输入 which 命令名命令名 进行查找 或者直接使用 MATLAB 提供的帮 助 看该命令是什么意思 该怎么使用 上上图给出了帮助的位置 下图是 帮助 窗口 在红 框内输入如 wk1read 点回车即可获取 wk1read 的相关解释 用法 描述 实例等等 目前关 于 MATLAB 用于空间面板计量方面的中文操作书我没有在网上找到 所以也只能通过这种方法来 学习 panel 下的每一个 M 文件到底是干什么用的 三 常见的命令 三 常见的命令 1 关于 后面可以输入中文 对已写程序进行描述和解释 后的文字是不被运行的 2 wk1read 读取后缀名为 wk1 的数据文件 wk1 数据文件是 MATLAB 能识别的数据文件 但具体该文件怎么形成的 我也不清楚 我一般使用的都是通过 EXCEL 链接导入到 MATLAB 中 的数据 该数据存储后后缀名为 mat 该文件的读取 直接使用命令 load 如读取名为 如读取名为 pc 的的 mat 文件 直接用文件 直接用 load pc mat 3 在任何一个空间面板数据模型的程序中 首先要做的 当然是读取数据 第二步就是对空 间面板数据进行描述 也叫定义空间面板数据 同时指定被解释变量 y 和解释变量 x 如下图 每一行分别表示 该面板数据的时期数为 30 T 30 该面板数据有 30 个地区 N 30 将空间权重矩阵标准化 W normw w1 将名为 consume 以矩阵形式出现在 MATLABA 中 的变量的第 3 列数据定义为被 解释变量 y 将名为 consume 的变量的第 4 5 7 8 9 10 12 列数据定义为解释变量矩阵 x 定义一个有 N T 行 1 列的全 1 矩阵 该矩阵名为 xconstant ones 即为全 1 矩阵 说明解释变量矩阵 x 的大小 有 nobs 行 K 列 size 为描述矩阵的大小 这里说明一下矩阵的引用 a 4 6 表示的是矩阵的第四到第六列 逗号前是行号 逗号后是列号 行号 列号必 须加中括号 如果使用连续几行或几列 如第一行到第七行 表示成 1 7 即可 如果几行将是间 断的 用逗号隔开 如果是全部的行或列 可省略数值和中括号 直接打冒号 即可 4 一般来说使用 panel 下的完整程序 在读懂意思后 我们用于自己的模型时 需要修改 的也就是读取的文件名 面板数据的定义和解释变量 x 被解释变量 y 的指定 以及在输出时变量的 名称 前面的都讲了 现在说说变量的名称 下面是一段最小二乘法回归以及 LM 检验 moran 检验的程序 results ols y xconstant x 表示以y作为被解释变量 x作为解释变量的含有 常数项的最小二乘法回归 xconstant即表示含有常数项 该命令运行结果全部存储在 results中 vnames strvcat gy intercept logy0 logpc gpc emp edu inv fir 即是对最终显示在控制窗口中的变量进行命名 将此名字赋值给 vnames 变 量 prt reg res