第二讲 Stata快速入门 高级计量经济学及Stata应用课件.pptx

第二讲 Stata快速入门 高级计量经济学及Stata应用课件.pptx

《第二讲 Stata快速入门 高级计量经济学及Stata应用课件.pptx》由会员分享，可在线阅读，更多相关《第二讲 Stata快速入门 高级计量经济学及Stata应用课件.pptx（77页珍藏版）》请在冰豆网上搜索。

22002200-0-50-525-25,陈强计量及Stata应用（c）,2014,1,1,高级计量经济学及Stata应用第二讲Stata快速入门,陈强山东大学经济学院,2020-05-25,陈强计量及Stata应用（c）2014,2,为什么使用Stata？

Stata是目前在欧美最为流行的计量软件，具有操作简单、功能强大的特点。

由于使用Stata的用户很多，对于最新的计量方法，常常可以下载由用户写的Stata命令程序（user-writtenStatacommand），十分方便。

而官方的Stata版本也经常更新，以适应计量经济学迅猛发展的需要。

本课主要介绍Stata12（Stata13刚推出）。

2020-05-25,陈强计量及Stata应用（c）2014,3,Stata的窗口,在最上方有一排菜单，即“FileEditDataGraphicsStatisticsUserWindowHelp”。

左上“Review”（历史窗口）：

此窗口记录着自启动Stata以来执行过的命令。

右上“Variables”（变量窗口）：

此窗口记录着目前Stata内存中的所有变量。

正上方“Results”（结果窗口）：

此窗口显示执行Stata命令后的输出结果。

正下方“Command”（命令窗口）：

在此窗口输入想要执行的Stata命令。

2020-05-25,陈强计量及Stata应用（c）2014,4,2020-05-25,陈强计量及Stata应用（c）2014,5,Stata操作实例,学习Stata的最方便方法大概是通过实例来学习。

因此，这里选取Nerlove（1963）对电力行业规模报酬的经典研究来介绍Stata的实际操作。

从本书的网站可以下载Nerlove（1963）论文原文，以及数据集“nerlove.xls”（Excel文件）。

该数据集包括了1955年美国145家电力企业的横截面数据。

2020-05-25,陈强计量及Stata应用（c）2014,6,将数据导入Stata,打开Stata软件后，点击DataEditor（Edit）图标（也可以点击菜单“Window”“DataEditor”），即可打开一个类似Excel的空白表格。

用Excel打开文件“nerlove.xls”，复制文件中的所有数据，并粘贴到DataEditor中。

导入数据的另一方法是，点击菜单“File”“Import”，然后导入各种格式的数据。

但这种方法有时不如直接从Excel表中粘贴数据来得方便直观。

2020-05-25,陈强计量及Stata应用（c）2014,7,变量窗口,关闭DataEditor后，即会看到右上方的“Variables”窗口出现了5个变量：

分别为tc（totalcost，总成本），q（totaloutput,总产量），pl（priceoflabor，小时工资率），pf（priceoffuel，燃料价格），与pk（usercostofcapital，资本的租赁价格。

2020-05-25,陈强计量及Stata应用（c）2014,8,存为dta数据文件,此时，可以点击Save图标（也可以点击菜单“File”“Save”），将数据存为Stata格式的文件（扩展名为dta），比如nerlove.dta。

以后就可以用Stata直接打开这个数据集了（不需要再从Excel表中粘贴过来）。

2020-05-25,陈强计量及Stata应用（c）2014,9,打开dta数据文件,打开的方式有三种：

点击Open图标（也可以点击菜单“File”“Open”），然后寻找要打开的dta文件的位置。

直接双击想要打开的dta文件在命令窗口输入以下命令（假设文件在E盘的根目录）并回车（按Enter键）,2020-05-25,陈强计量及Stata应用（c）2014,10,打开数据文件的命令,useE:

nerlove.dta,clear其中，选择项“clear”表示可以替代内存中的已有数据。

如果要关闭一个数据集（如果对数据集进行了改动，别忘了先存盘Save），以便使用另外一个数据集，可以在命令窗口输入clear这样，内存中所有的当前数据都被清空，然后可以再打开另外一个数据集。

2020-05-25,陈强计量及Stata应用（c）2014,11,变量的标签,在变量窗口，每个变量的“名字”（Name）旁边显示了其“标签”（label）。

但目前的标签过于简略，缺乏变量的解释信息。

如果想将变量“tc”的标签改为“totalcost”，可进行如下操作。

点击进入“变量管理器”（VariablesManager）图标（在DataEditor右侧）。

Stata中字母的大小写是严格区分的（casesensitive），因此Stata建议对于变量名一律使用小写字母。

2020-05-25,陈强计量及Stata应用（c）2014,12,审视数据,一个数据集可能很大，而我们常希望看到数据的概貌。

想看数据集中的变量名单、标签等，可以在命令窗口输入，describe其中，“describe”中的下划线表示，可以将该命令简写为“d”而得到同样的效果。

2020-05-25,陈强计量及Stata应用（c）2014,13,2020-05-25,陈强计量及Stata应用（c）2014,14,罗列数据,如果想看变量tc与q的具体数据，可使用命令，listtcq由于样本容量为145，要花一阵子才能显示完毕。

如果想中途停止该命令的执行，则可以点击Break图标，或直接在键盘上同时按“Ctrl+Break”。

运行结果如下。

中断命令,2020-05-25,陈强计量及Stata应用（c）2014,15,2020-05-25,陈强计量及Stata应用（c）2014,16,旧命令的调用,如果改变主意，仍然希望显示变量tc与q的全部数据，那么只要把光标放在命令窗口，并按键盘上的“PageUp”键即可调用上一个命令（反之，使用“PageDown”键可调用下一个命令）。

另一种简便的方法是，在左上角的历史窗口点击任何曾用过的命令：

如果用鼠标单击旧命令，则会把旧命令重新调入命令窗口，按回车后即执行，或将旧命令进行编辑后再执行；如果用鼠标双击旧命令，则将马上自动执行。

定义子集,有时我们想对数据集的一部分执行命令，比如只想看变量tc与q的前5个数据，则可输入命令：

listtcqin1/5,2020-05-25,陈强计量及Stata应用（c）2014,17,逻辑关系,也可以通过逻辑关系来定义数据集的子集。

如果要列出所有满足条件“”的变量tc与q的数据，则可以使用以下命令，listtcqifq=10000,2020-05-25,陈强计量及Stata应用（c）2014,18,2020-05-25,陈强计量及Stata应用（c）2014,19,只对子集进行运算,如果想删除满足“q=10000”条件的观测值，则可使用命令，dropifq=10000反之，如果只想保留满足“q=10000”条件的观测值，而删去所有其他观测值，keepifq=10000,考察变量的统计特征,如果想看变量q的统计特征，可输入命令，summarizeq,2020-05-25,陈强计量及Stata应用（c）2014,20,考察变量的统计特征（续）,如果不指明变量，则将显示数据集中所有变量的统计指标。

su,2020-05-25,陈强计量及Stata应用（c）2014,21,考察变量的统计特征（续2）,如想看更多的统计指标，可使用命令suq,detail,2020-05-25,陈强计量及Stata应用（c）2014,22,经验累积分布函数,如果要显示变量pl的经验累积分布函数（empiricalcumulativedistributionfunction），可使用命令，tabulatepl,2020-05-25,陈强计量及Stata应用（c）2014,23,相关系数,如果要显示内存中5个变量之间的相关系数，可输入命令，correlatetcqplpfpk,2020-05-25,陈强计量及Stata应用（c）2014,24,2020-05-25,陈强计量及Stata应用（c）2014,25,相关系数（续）,pwcorrplpfpk,sigstar（.05）“pw”表示pairwise。

选择项“sig”表示显示相关系数的显著性水平（即p值，列在相关系数的下方），选择项“star（.05）”表示给所有显著性水平小于或等于5%的相关系数打上星号。

如pwcorr之后没有指定变量，显示所有变量的相关系数,2020-05-25,陈强计量及Stata应用（c）2014,26,直方图,如果想看变量q的直方图（假定组宽为1000），可输入以下命令（也可通过菜单来输入此命令）：

histogramq,width（1000）frequency,2020-05-25,陈强计量及Stata应用（c）2014,27,散点图,如果要画tc与q之间的散点图，则可输入以下命令：

scattertcq,2020-05-25,陈强计量及Stata应用（c）2014,28,2020-05-25,陈强计量及Stata应用（c）2014,29,在散点图上标注观测值,在上页的散点图中，无法知道每个点分别对应哪个观测值为此，首先定义一个新变量“n”来表示第n个观测值。

genn=_n其中，“_n”即表示第n个观测值。

scattertcq,mlabel（n）mlabpos（6）选择项“mlabel（n）”表示以变量“n”作为“marklabel”（标签）；“mlabpos（6）”（marklabelposition）表示将此标签放在散点正下方（6点钟的位置），默认位置为散点的右边（3点钟）。

2020-05-25,陈强计量及Stata应用（c）2014,30,在散点图上画回归直线,twoway（scattertcq）（lfittcq）其中，“lfit”表示“linearfit”（线性拟合）graphsavescatter1,2020-05-25,陈强计量及Stata应用（c）2014,31,在散点图上画二次回归曲线,twoway（scattertcq）（qfittcq）其中，“qfit”表示“quadraticfit”（二次拟合）。

graphsavescatter2,2020-05-25,陈强计量及Stata应用（c）2014,32,将两图合并,graphcombinescatter1.gphscatter2.gph,2020-05-25,陈强计量及Stata应用（c）2014,33,Nerlove（1963）对电力企业的研究,假设生产函数为Cobb-Douglas函数：

可以证明成本函数也为Cobb-Douglas函数：

为规模效应。

取对数可得：

2020-05-25,陈强计量及Stata应用（c）2014,34,2020-05-25,陈强计量及Stata应用（c）2014,35,生成新变量,可使用命令“generate”来生成新变量。

glntc=log（tc）glnq=log（q）glnpl=log（pl）glnpf=log（pf）glnpk=log（pk）,2020-05-25,陈强计量及Stata应用（c）2014,36,生成新变量（续）,如果需要q的非线性平方项，可使用命令：

gq2=q2如果要生成lnpl与lnpk的互动项（interactionterm），则可以使用命令，glnplpk=lnpl*lnpk,2020-05-25,陈强计量及Stata应用（c）2014,37,生成虚拟变量,假设希望定义“q=10000”为大企业，并使用“虚拟变量”（dummyvariable）