Stata上机实验笔记
实验报告模板作业
实验报告
一、实验目的和要求
学习stata软件的应用要求自带电脑和老师要求准备的课件
二、实验原理
stata的上机操作
三、主要仪器设备、试剂或材料
电脑、课件以及上课要用到的关于stata的文件
四、实验方法与步骤
输入list
输入list make
输入scatter wage educ
输入line wage educ
输入line wage educ
输入twoway area wage educ,sort
输入histogram wage
输入histogram lwage
输入histogram lwage if married==1
输入histogram lwage if married==0
五、实验数据记录、处理及结果分析
在stata的应用过程中数据都记录在软件中
六、讨论、心得
通过对stata软件的上机练习，对其有了初步的认识，对计量经济学的研究学习有更进一步的帮助。
stata上机实验第五讲 工具变量(IV)
究竟该用OLS 还是IV
即解释变量是否真的存在内生性？ 假设能够找到方程外的工具变量。 1。如果所有解释变量都是外生变量，则OLS
比IV 更有效。在这种情况下使用IV，虽然估 计量仍然是一致的，会增大估计量的方差。2。 如果存在内生解释变量，则OLS 是不一致的， 而IV 是一致的。
豪斯曼检验（Hausman specification test）原假设： H0 ：所有解释变量均为外生变量。 H1：至少有一个解释变量为内生变量。
检验方法： estat firststage 1。初步判断可以用偏R2(partial R2) (剔除掉模型中原有外生变量的影响)。 2。 Minimum eigenvalue statistic(最小特征 值统计量)，经验上此数应该大于10。
ivregress 2sls lw80 expr80 tenure80 (s80 iq=med kww mrt age), first
使用grilic.dta估计教育投资的回报率。
变量说明：lw80（80年工资对数），s80 （80年时受教育年限），expr80（80年时工 龄），tenure80（80年时在现单位工作年 限）， iq（智商），med（母亲的教育年 限），kww（在‘knowledge of the World of Work’测试中的成绩），mrt（婚姻虚拟变量， 已婚=1），age（年龄）。
ivregress 2sls lw80 expr80 tenure80 (s80 iq=med kww mrt age), first estat overid ivregress gmm lw80 expr80 tenure80 (s80 iq=med kww mrt age) estat overid
(最新整理)stata上机实验第五讲——面板数据的处理..
• corr(u_i, Xb) 个体效应与解释变量的相关系数，相关 系数为0或者接近于0，可以使用随机效应模型；相 关系数不为0，需要使用固定效应模型。u-i不表示残 差，表示个体效应。
2021/7/26
19
模型选择
• 固定效应还是混合OLS？ 可以直接观测F值
• 随机效应还是混合OLS？ 先用随机效应回归，然后运行xttest0
• xtgls Panel-data models using GLS
• xtpcse OLS or Prais-Winsten models with panelcorrected standard errors
• xtrchh Hildreth-Houck random coefficients models
（完整版）Stata学习笔记和国贸理论总结
（完整版）Stata学习笔记和国贸理论总结Stata学习笔记⼀、认识数据（⼀）向stata中导⼊txt、csv格式的数据1．这两种数据可以⽤⽂本⽂档打开，新建记事本，然后将相应⽂档拖⼊记事本即可打开数据，复制2．按下stata中的edit按钮，右键选择paste special3．*.xls/*.xlsx数据仅能⽤Excel打开，不可⽤记事本打开，打开后会出现乱码，也不要保存，否则就恢复不了。
逗号分隔的数据常为csv数据。
（⼆）⽹页数据⽹页上的表格只要能选中的，都能复制到excel中；⽹页数据的下载可以通过百度“国家数据”进⾏搜索、下载⼆、Do-file 和log⽂件打开stata后，第⼀步就要do-file，记录步骤和历史记录，⽅便⽇后查看。
Stata处理中保留的三种⽂件：原始数据(*.dta)，记录处理步骤(*.do)，以及处理的历史记录(*.smcl)。
三、导⼊StataStata不识别带有中⽂的变量，如果导⼊的数据第⼀⾏有中⽂就没法导⼊。
但是对于列来说不会出现这个问题，不分析即可（Stata不分析字符串，红⾊⽂本显⽰；被分析的数据，⿊⾊显⽰）；第⼀⾏是英⽂变量名，选择“Treat first row as variable names”在导⼊新数据的时候，需要清空原有数据，clear命令。
导⼊空格分隔数据：复制——Stata中选择edit按钮或输⼊相应命令——右键选择paste special——并选择，确定；导⼊Excel 中数据，复制粘贴即可；逗号分隔数据，选择paste special后点击comma，然后确定。
Stata数据格式为*.dta，导⼊后统⼀使⽤此格式。
四、基本操作（⼏个命令）（⼀）use auto，clear 。
在清空原有数据的同时，导⼊新的auto数据。
（⼆）browse 。
浏览数据。
（三）describe 和list。
查看数据，describe 和list 使⽤list命令能使我们根据⾃⼰的需要选择数据（例如其与in/if语句的结合使⽤）。
stata笔记
1.一般检验假设系数为0，t比较大则拒绝假设，认为系数不为0.假设系数为0，P比较小则拒绝假设，认为系数不为0.假设方程不显着，F比较大则拒绝假设，认为方程显着。
2.小样本运用OLS进行估计的前提条件为：（1）线性假定。
即解释变量与被解释变量之间为线性关系。
这一前提可以通过将非线性转换为线性方程来解决。
（2）严格外生性。
即随机扰动项独立于所有解释变量：与解释变量之间所有时候都是正交关系，随机扰动项期望为0。
(工具变量法解决)（3）不存在严格的多重共线性。
一般在现实数据中不会出现，但是设置过多的虚拟变量时，可能会出现这种现象。
Stata可以自动剔除。
（4）扰动项为球型扰动项，即随即扰动项同方差，无自相关性。
3.大样本估计时，一般要求数据在30个以上就可以称为大样本了。
大样本的前提是（1）线性假定（2）渐进独立的平稳过程（3）前定解释变量，即解释变量与同期的扰动项正交。
（4）E（XiXit）为非退化矩阵。
（5）gt为鞅差分序列，且其协方差矩阵为非退化矩阵。
与小样本相比，其不需要严格的外生性和正太随机扰动项的要求。
4.命令稳健标准差回归：reg y x1 x2 x3, robust 回归系数与OLS一样，但标准差存在差异。
如果认为存在异方差，则使用稳健标准差。
使用稳健标准差可以对大样本进行检验。
只要样本容量足够大，在模型出现异方差的情况下，使用稳健标准差时参数估计、假设检验等均可正常进行，即可以很大程度上消除异方差带来的副作用对单个系数进行检验：test lnq=1线性检验：testnl _b[lnpl]=_b[lnq]^25.如果回归模型为非线性，不方便使用OLS,则可以采取最大似然估计法（MLE）,或者非线性最小二乘法（NLS）6.违背经典假设，即存在异方差的情况。
截面数据通常会出现异方差。
因此检验异方差可以：（1）看残差图，但只是直观，可能并不准确。
rvfplot (residual-versus-fitted plot) 与拟合值的散点图rvpplot varname (residual-versus-predictor plot) 与解释变量的散点图扰动项的方差随观测值而变动，表示可能存在异方差。
stata学习笔记（stata学习笔记）
stata学习笔记（stata学习笔记）data managementCreate a new dataEdit / / variables in the data table and the creation of open dataInput x1 x2......Set OBS 10Gen x1=_nGen, x2=seq ()Egen, x3=seq (), B (5) t (5)Egen x4=fill (3434)Rename X1 pop / / variable VAR1 renamed popRename x2 placeMax C= (1,0.8\0.8,1)Drawnorm, x1, X2, means (1,10), SDS (0.3,2), corr (C), n (500)Gen x1=invnormal (uniform ())Gen roll=1+trunc (uniform () *6) randomly generates 1-6 randomnumbersGen x=exp (uniform ())Gen x=-3ln (uniform ())Gen x= (invnorm (uniform ())) ^2 chi square distributionGen, x=invttail (DF, uniform ()) t distributionGen, x=invFtail (DF1, df2, uniform ()) F distributionSample 10, countLabel variable pop population in 1000s, 1995 "/ / add tags for the variable popLabel define, sex_label 1, "male", 2 "female""Label values sex sex_label / / add value labels for the variable sexSave AAA / / keep the aaa.dta fileSave, replaceMerge dataUse a.datAppend using B.datUse a.datSort placeSave, replaceUse B.datSort placeMerge place using a.datReshape, long, grow, I (ID), J (year)Reshppe, wide, grow, I (ID), J (year)ClearCD f:\ statistics \stataUse AAASort pop / / as the pop variable orderingOrder place pop place pop / / variables were placed in the first, second positionDescrible / / description variable informationList / / show variable and variable valuesList, Sep (3) is shown separately in each of the 3 linesList, sepby (VaR) is shown as bounded by the VaR variableSummarize X / / display basic information variables, can add "d" to display detailed informationBy, VAR1, var2, sort:su, X (by can be used for Su, CI, centile, etc.)Tabstat, x, stats (mean, median, SD,, VaR, skewness, kurtosis, IQR, CV, semean, P2, etc)Collapse (sum), VAR1, var2 (SD), var3 (mean), newvar1=var4 (median), newvar2=var5A subset of variables (used by if and in)List, pop, place, sex, in, 1/50Sort popList pop place in -4/1 / / four observation shows that the value of pop maximumSummarize if pop<1000Summarize if place = = "China""Summarize, pop, place, sex, if, pop>100 & pop<1000Summarize place sex if pop<100 pop>1000 |Summarize place if pop<. / / the missing value is bigger than any numericalDrop, pop, if, place==, "China""KeepCreate and replace variables1, use, canada1, clearGenerate gap=flife-mlife"Label variable gap" "flife-mlife gap life""Format gap%4.1f / / fixed width of 4 decimal 1Other%4.1g (width 4, decimal part at least 1, can be displayed by decimal or scientific notation),%4.1eFormat only changes the display and does not affect the calculationUse, canada1, clearGenerate type=1Replace, type=2, if, place==, "Canada""Replace, type=3, if, place==, "Yukou""operator+ * / ^ mod (x, y)Use function(ABS)ACOS () //di ACOS (0.5) *180/_piSin, cos, asin, atan, atan2 () y/x's tangent functionSqrt, log (), ==ln (), log10, expThe smallest integer of ceil (x) >xThe maximum integer of floor (x) <xRound (x) four into fiveComb () lnfactorial ()distribution functionProbability of Ttail (DF, t) t>t0.05 (Dan Ce)Invttail (DF, P) calculates the T value based on the probability, and P is the right probabilityF (DF1, df2, f) left probability invF (N1, N2, P)Ftail (DF1, df2, f) the right probability invFtail (N1, N2, P)Chi2 (DF, x) left probabilityChi2tail (DF, x) right probabilityBinomial (n, x, P), n trials, x times and smaller probability1-binomial (n, X-1, P)Normal (z) standard normal distribution, left, cumulative probabilityDate function(1) assume that the numeric variable a is 20100312Gen str str_a=string (a,%10.0f) / / a conversion to character variableGene _ to date = DATE ("STR _, Ymd") / / 转换str _ a为日期变量, 返回值为当前日期 - 1960年1月1日的数值FORMAT DATE _% TD / / 转换date _ a的格式为日期12may2010假设有数值变量a格式为20100312101205STR str Gene _ = String ("% 16.0f")To _ = Clock Gene Double Date (STR _, "ymdhms")_% TC to date format假设有三个数值变量m、d、y分别表示月、日、年Gene _ date to mdy = (m, d)EgenEgen = seq (x t), B (3) (2) 111222111222Egen fill (x = 100,98) 100 98 94 96X = (0,2,7,0,2,7 egne fill)Rowmean egen x = (x1, X2, x3) 产生新变量, 其值为x1x2x3各行的均值Rowsum egen x = (x1, X2, x3) 产生新变量, 其值为x1x2x3各行的和Egen = STD X (a)Num 1: 15 for STD / egen xx = (AX)Xrank egen = RANK (X)10、其他函数Recode Group encodeX1 = recode gene (AGE, 24,28,32, ~) / / < < = 24 = 28Egen Group (x2 = x1)Strvar Gene ENCODE, 将字符变量转为数值变量 (numvar)Decode numvar, Gene (strvar)创建新的分类变量和定序变量假设有分类变量 (byte) type (1 - 3)Tab typeTab type, Gene (type) / / 产生type1 - 3三个哑变量2、将数值变量X1 = recode gene (AGE, 24,28,32, ~) / / 以 < < = 24 = 28～分组Egen Group (x2 = x1)变量下标Di x [4]Gene _ = X - X [N - 1] / / x与其前一个数值的差B gene _ = X - X [n + 1]从外部ascii文件导入数据以空格分隔, 字符串需带引号Str30 INFILE Place ulife tlife using aaa.raw / / 产生三个变量, place为30长度的字符变量COMPRESS / / 压缩place变量为最长的字符以tab或 "," 分隔Insheet Place ulife tlife using aaa.raw, comma (or tab).固定栏宽Infix Wood Year 1 - 4 5 - 8 9 - 10 aaa.raw using Water绘图Hist X, Bin (10) xlabel (0 (2) 10) ylabel (100 1000 xtick (100) (1) (2) 11) Norm fractionHist Start (50 x width (5) (FREQ by Group, total)Graph TwoWay Scatter and | | X Line and | | lfit X and X, mlabel (ID) msymbol (o / X)Graph TwoWay Scatter and x | | lfitci, STDFGraph Matrix X and ZGraph TwoWay line and year XGraph TwoWay line and yaxis (1 year) | yaxis | x Year (2)Graph TwoWay area and year XGraph box x and Z over (Group) yline (6.35).Graph pie x and Z, by (Group) foot (3, explode)Graph BAR (Mean) of X and Z, over (Group)Grapg DOT (median) x1 x2, over (Group) Marker (1, msymbol (OH) (2) Marker, msymbol (X))X Qnorm, GridPnorm X, Grid交叉表Tab B, SUM (X) meanTab B, All tabi B \ C D, All tab b] [FW = count, AllA B C 分布绘制abc的一维表 tab1A B C 建立所有可能的二维表 Tab2Sort by: a B C, Tab, All 以c的不同取值分别绘制a b的二维表Table Row col (col1, by 绘制多维表 row1)Sktest x swilk sfrancia正态性检验及数据变换Sktest x swilk sfrancia立方严重负偏态平方轻度负偏态平方根轻度正偏态对数正偏态平方根负倒数严重正偏态倒数非常严重正偏态平方倒数同上立方倒数同上X / / 产生以上8种变换后的正态性检验 LadderGladder X / / 针对ladder结果绘制直方图Bcskews newx = X / / 产生新变量newx, 是对x的变换方差齐性检验Sdtest X1 = x2Sdtest X1, by (Group)Robvar X, by levene检验, 返回值 (Group)W0: 均数 W50: 中位数 W10: 后的均数 trim10%方差分析单个样本TTEST (x = 10 signtest x = 10 二项分布ttest x1 = x2 signrank x1 = x2 wilcoxon符号检验ttest x city (group) ranksum x city (group) wilcoxon检验ttest x1 = x2, unpaired unequalbitest x = = pbitesti n c p单因素方差分析oneway x group, tabluate scheffe bonferroni sidak kwallis x city (group)多因素方差分析anova x a # # btest 1 (a = (test 2 (b = 3. bbonferonni: r (p) * c c: 比较次数, 组数x (组数 - 1) / 2scheffe: 1 - f (组数 - 1, 误差自由度, r (f) / (组数 - 1))regresspredict newvar 预测值predict newvar, stdp 预测值标准误anova x a b | aanova x a / id | a b a # banova x a b c.age相关分析 (town was:)cor x ypwcorr x y, bonferrior / sidakspearman x y, bonferrior / sidakpcorr y x1 - x3 去除其他x的影响后y与x的偏相关系数回归分析基本方法reg y x1 x2 x3, beta uncons预测值predict newvar, cooksd hat covratio dfits residuals rstudent rstandard stdp stdfhat > 2p / n 发现高杠杆值dfits > 2sqrt (p / n) 案例的自变量组合对回归直线的影响力cooksd > 4 / n 同上welsch > 3sqrt (p) 同上covratio: | r - 1 | > = 3p / nrvfplot, yline (0)假设检验reg x * ytest x1 x2 x1 和x2回归系数同时为0test x1 = x2虚拟变量loss region gene (reg) / / 产生reg1 - 4四个哑变量reg cmat reg2 / / reg2与其他3个地区的比较reg cmat reg1 reg2 reg3 reg4 = = xi: reg cmat i.region 此方法便于做交互分析char region [omit] 4 (与xi共同使用)xi: reg camt i.region逐步回归sw reg y x1 - x4, per (. 06) pe (0.05)sw reg y x1 x2 (x3, x4) lockterm1 per (. 06)面板数据iis regionten yearxtreg y x1 x2, rextmixed y 固定变量 | | school: 随机变量回归诊断estate ic 返回aic bic ll (null) ll (model) 值 (log likelihood 对数似然值)quietly reg y x1 - 85estimates of large fullquietly reg y x1 - x4lrtest fullovtest p < 0.05提示有二次、三次或四次方项目需要添加hettest p < 0.05提示方差不齐, 误差散点图不是随机分布的dwstat 一价自相关的durbin - watson检验kic 自变量共线性检查kic > 10 平均vif > 1 有问题宽容度 (vif的倒数) 表示该变量独立程度, 越大则越独立rvfplot 预测值与残差值的散点图rvpplot x 某一个自变量x与残差的散点图avplot x 去除其他变量影响后的x与y的线性关系, x轴上偏离的数值多为高杠杆值avplotsacprplot x, lowess 虚线在中间部分与直线不重和表示可能x与y 存在其他非线性关系,另外可以报告与x具有线性关系的其他自变量lvr2plot 注意拟合不好且具有较高杠杆作用的值可能是高杠杆值hat 较大值提示高杠杆值dfits cooksd covratio 提示对y影响较大的值logistic回归logit y x * logit y x *, orblogit n x * ylrocroctab y x, graphroccomp y x1 x2 比较y与x1的roc曲线和y与x2的是否相同rocgold y x x1 x2 比较y与x (金标准) 的roc曲线和y与x1的是否相同lsens, genprob (prob) gensens (sen) genspec (spec)lstatlfit, group (10) est gof, group (10)predict the phat, hat deviance ddeviance dbet dx2 dbetaclogit y x *, group (matchvar)ologit x * ymlogit y x *, b (1) mlogit y x *, rrrconstraint define 1 [3] x = 2 [2] xconstranit define 2 [4] x = 3 [2] xmlogit y x, c (1, 2) b (1)多元方差分析hotelling x *hotelling x *, city (group)manova x1 x2 x3 = g b g * b广义线性模型gaec y x1 x2家庭（高斯）链接（身份）* /正态分布线性回归GLM y x1 x2，家庭（二项式）链接（Logit）* /物流回归GLM y x1 x2，家庭（Poisson）链接（日志）lnoffset（暴露人年变量）泊松y x1 x2，曝光（暴露人年变量）poisgof [皮尔森]GLM y x1 x2，家庭（nbinomial）链接（日志）nbreg y x1 x2gnbreg y x1 x2，lnalpha（VAR）预测主成份分析PCA X点状图因子分析X因子*，PCF矿（0.5）主成份法X因子*，ML矿（0.5）最大似然法X因子*，IPF /迭代主因子法旋转方差极大旋转旋转，旋转斜交法生存分析认识时间，失败（结果）stsum，由（治疗）后缀树，由rmean（处理）STS列表，由（治疗）以损失为例的STS图STS图，通过gwood（治疗）STS测试组streg治疗组，诺尔公司（指数/ Weibull）预测new_var = = 1如果治疗，监测stcox治疗组，诺尔考克斯结果治疗组，死亡（结果）诺尔申银万国考克斯结果治疗组，死亡（结果）诺尔流行病队列研究IR案例的曝光时间（人年数）硝酸铵CS案例曝光[或数]CSI（A组）可使用结核病和精确（默认）计算RR可信区间，不能使用伍尔夫病例对照研究cc案例由（组）公开甲丙氨酯MCC的病例对照选择A B C D可使用精确（默认）、伍尔夫、麦田计算RR可信区间tabodds模型不如用物流、考克斯比例风险模型。
stata上机实验第八讲 似不相关回归(SUR)
use hsb2,clear sureg (read write math science) (socst write
math),corr isure test [read]math=[socst]math
例2：用三家公司的公司投资额对公司市值、 资本存量进行回归。（grunfeld2.dta）
Invest11011mvalue112kstock11 Invest22021mvalue222kstock22 Invest33031mvalue332kstock33
sureg (invest1 = mvalue1 kstock1) (invest2 = mvalue2 kstock2) (invest3 = mvalue3 kstock3)
这种迭代估计方法必须设定初始值和停止法 则。初始值的选择对于迅速找到最优解非常 重要。
例1：利用NLS方法估计非线性消费函数（数 据文件：usmacro）
csinc u
nl (realcons = {a} + {b}*realgdp^{gamma=1}) 如果不给定gamma的初始条件将无法达到收
分位数回归
传统回归模型着重考察解释变量x对被解释变量y的 条件期望的影响，实际上是“均值回归” 。但这种 方法容易产生如下问题：1。无法了解y的整体分布； 2。结果受极端值影响严重。
如果能够估计出条件分布的若干重要的“条件分位 数”（conditional quantiles），比如“中位数” （median）、“分位数”（lower quartile）、“分 位数”（upper quartile），就能对条件分布有更全 面的认识。
4。利用自助法重复100次计算[0.25，0.75]的 分位数回归。
Stata笔记 北京科技大学
li(st) if new2>=14 & new2<24 (按q可以退出，即quit）
replace new3=rep78 (输错了替换)
drop new new2 new3删除变量
list if new>10000
list make if new<10000|new>2000（竖线表示或者，回车上面那个）
reg因变量（因变量：被解释变量；自变量：解释变量）
set linesize 140加宽显示窗格
reg price rep78 headroom trunk weight (std. err.=standard error,T值,confidence interval置信区间)
P与T反向关系，标准误=T值
. save bs
file bs.dta saved
. clear
. import excel "D:\Stata\BS_AD8BD666A41_(1)0417\newIS_8CBBC6F6EBB_(1)0417.xls", sheet("IS") firstrow
. duplicates drop _Stkcd y,force删除重复的
data
而master data指主数据库，using data从数据库
（替换变量名Stked改成A_Stked）repA_Stked Stked
（替换原数据is3）save is3，replace
clear清除合并
use bs3
duplicates report A_Stkcd year报告重复的变量（两个）
stata上机实验第六讲
不同，其参数估计值并不直接可比。雨宫 （Amemiya）提出：同一个模型的logit和 probit模型大概具有如下关系： Blogit 约等于 1.6*BProbit 但利用mfx计算的两者的边际效应应该大致相 同。
一个综合例子
使用美国妇女就业数据集“womenwork.dta”，
估计决定就业的Probit与Logit模型。被解释 变量为work=1，就业；work=0，不就业。解 释变量为age（年龄）, married（婚否）, children（子女数）, education（教育年限）。
排序选择模型
根据GSS的调查数据，不同的家庭母亲与子 女之间的关系也不同。根据调查显示，有的 家庭母子（女）关系比较紧张，有的比较融 洽。变量包括：warm=关系融洽度（0、1、2、 3）；educ=子女接受教育的程度；age=子女 年龄；male=儿子；prst=职业威望；white= 白人；y89=89年调查结果。分析不同因素对 母子（女）关系的融洽程度有何影响。
sysuse auto,clear truncreg price weight length gear_ratio, ll(10000) reg price weight length gear_ratio if price>=10000
截取回归
对于线性模型yi =x’B+εi，如果满足yi ≥ c或者
nbreg daysabs langarts male 命令结果中将提供一个LR 检验， 原假设：不存在过度分散，应该使用泊松回 归。此时alpha=0。 备则假设： alpha<>0，不能使用泊松回归。
受限因变量模型
1。断尾回归
2。截取回归
stata统计分析与应用笔记汇总
第一章：Stata概述：help和search都是查找文件的命令但help用于查找精确的命令，search是模糊查找。
还可使用help|contents 来分类查找第二章：数据管理：2.1变量和变量的取值：1.变量的命名：不能以数字开头，区分大小写，不能命名为系统变量名2.变量的取值类型：（1)字符型：字符变量存储格式是str⋕，str表示格式⋕表示该变量的存储最多可容纳的字符数（2）数值型数据：存储格式：byte.int.long.float.double.Stata默认将数字存储为浮点数据，而将计算结果存为双浮点数据。
（3）缺失数据：一般仅用“.”表示3．变量的显示：（1）数值变量的显示格式：a.普通格式有%w.dg, %w.dgc（g表示普通，w表示整个显示所占的字符数，d表示显示的数字中小数点后的位数,c是要求Stata给出带逗号“,”数字显示格式如12345显示为12，345）b．固定格式有%w.df, %w.dfc（f表示固定）c.科学指数法格式：%w.de, (e表示科学计数)（2）字符变量的显示格式：仅有一种%⋕s，%是提示符，#表示显示字符数，s表示字符变量显示格式,默认右对齐，后加“-”可改为左对齐。
（3）使用format命令变量显示格式：format varlist %fmt 或者 format %fmt varlist 4.变量的标签（1）添加数据集的标签使用： label data [“lable”]（2）添加变量的标签使用：label variable varname [“lable”]（3）label为变量数值添加标签的语法有两部分，先定义数值标签：label define lblname#“lable” [#“lable”](lblname是标签名称) 然后将定义好的数值标签添加到变量上:label values varlist [lblnamel.]2.2创建一个新的数据集1.关于数据集操作的基本命令（1）browse 和edit 命令：browse 用于打开数据浏览器，edit命令用于打开数据编辑器Edit [varlist] [if] [in]browse [varlist] [if] [in](if和in 用于选择需要的子集)（2）rename:rename old_varname new_varname（3）save命令：save [filename] [,save_options]（[,save_options]可以指nolabel(不保存设定标签),replace（允许新文件覆盖原文件）,all主要用于编程（4）describe:用于产生一个对数据集的简明总结格式：describe [varlist] [,memory_options]（命令选项：simple,short,detail,fullnames）（5）list:用于显示变量的数值，其后可以跟需要显示的变量名称语法：list [varlist] [if] [in] [,options](命令选项包括：noobs(不显示观测值的数值)，clean,separator,sepby,nolabel)（6）codebook:用于详尽地描述变量的内容，包括变量名称、标签、赋值。
stata上机实验操作
第六章第二题：1. 建立完成的教育年数（ED ）对到最近大学的距离（Dist ）的回归：. reg ed dist, robust斜率估计值是：-0.0732. reg ed dist bytest female black hispanic incomehi ownhome dadcoll cue80 stwmfg80,robustDist 对ED 的效应估计是：-0.0323. 系数下降50%，存在很大差异，（1）中回归存在遗漏变量偏差4. di e(r2_a)（可看到调整后的R2）第一问中=0.0074 调整的2R =0.00718796_cons 13.95586 .0378112 369.09 0.000 13.88172 14.02999dist -.0733727 .0134334 -5.46 0.000 -.0997101 -.0470353ed Coef. Std. Err. t P>|t| [95% Conf. Interval]RobustRoot MSE = 1.8074R-squared = 0.0074Prob > F = 0.0000F( 1, 3794) = 29.83Linear regression Number of obs = 3796. reg ed dist , robust2R第二问中=0.2788 2R = 0.27693235可以得到第二问中的拟合效果要优于第一问。
第二问中相似的原因：因为n 很大。
5. Dadcoll 父亲有没有念过大学：系数为正（0.6961324）衡量父亲念过大学的学生接受的教育年数平均比其父亲没有年过大学的学生多。
-.0517777 1）原因：这些参数在一定程度上构成了上大学的机会成本。
2）它们的系数估计值的符号应该如此。
当Stwmfg80增加时，放弃的工资增加，所以大学入学率降低了；因而Stwmfg80的系数对应为负。
最新STATA实用学习笔记资料
最新STATA实⽤学习笔记资料北京科技⼤学STATA应⽤学习摘录第⼀章 STATA的基本操作⼀、设置内存容set mem 500m, perm⼀、显⽰输⼊内容Display 1Display “clive”⼆、显⽰数据集结构describeDescribe /d三、编辑editEdit四、重命名变量Rename var1 var2五、显⽰数据集内容list/browseList in 1List in 2/10六、数据导⼊:数据⽂件是⽂本类型（.csv）1、insheet: . insheet using “C:\Documents and Settings\Administrator\桌⾯\ST9007\dataset\Fees1.csv”, clear2、内存为空时才可以导⼊数据集，否则会出现（you must start with an empty dataset）（1）清空内存中的所有变量：.drop _all（2）导⼊语句后加⼊“clear”命令。
七、保存⽂件1、save “C:\Documents and Settings\Administrator\桌⾯\ST9007\dataset\Fees1.dta”2、save “C:\Documents and Settings\Administrator\桌⾯\ST9007\dataset\Fees1.dta”, replace ⼋、打开及退出已存⽂件use1、.Use ⽂件路径及⽂件名, clear2、. Drop _all/.exit九、记录命令和输出结果（log）1、开始建⽴记录⽂件：log using "J:\phd\output.log", replace2、暂停记录⽂件：log off3、重新打开记录⽂件：log on4、关闭记录⽂件：log close⼗⼀、创建和保存程序⽂件：（doedit, do）1、打开程序编辑窗⼝：doedit2、写⼊命令3、保存⽂件，.do.4、运⾏命令：.do 程序⽂件路径及⽂件名⼗⼆、多个数据集合并为⼀个数据集（变量和结构相同）纵向合并append insheet using "J:\phd\Fees1.csv", clearsave "J:\phd\Fees1.dta", replaceinsheet using "J:\phd\Fees2.csv", clearappend using "J:\phd\Fees1.dta"save "J:\phd\Fees1.dta", replace⼗三、横向合并，在原数据集基础上加上另外的变量merge1、insheet using "J:\phd\Fees1.csv", clearsort companyid yearendsave "J:\phd\Fees1.dta", replacedescribeinsheet using "J:\phd\Fees6.csv", clearsort companyid yearendmerge companyid yearend using "J:\phd\Fees1.dta"save "J:\phd\Fees1.dta", replacedescribe2、_merge==1 obs. From master data_merge==2 obs. From using data_merge==3 obs. From both master and using data⼗四、帮助⽂件：help1、. Help describe⼗五、描述性统计量1、summarize incorporationyear 单个summarize incorporationyear-big6 连续多个summarize _all or simply summarize 所有2、更详细的统计量summarize incorporationyear, detail3、centilecentile auditfees, centile(0(10)100)centile auditfees, centile(0(5)100)4、tabulate不同类型变量的频数和⽐例tabulate companytypetabulate companytype big6, column 按列计算百分⽐tabulate companytype big6, row 按⾏计算百分⽐tab companytype big6 if companytype<=3, row col 同时按⾏列和条件计算百分⽐5、计算满⾜条件观测的个数count if big6==1count if big6==0 | big6==16、按离散变量排序，对连续变量计算描述性统计量：（1）by companytype, sort: summarize auditfees, detail（2）sort companytypeBy companytype:summarize auditees⼗六、转换变量1、按公司类型将公开发⾏股票公司赋值为1，其他为0gen listed=0replace listed=1 if companytype==2replace listed=1 if companytype==3replace listed=1 if companytype==5replace listed=. if companytype==.⼗七、产⽣新变量genGenerate newvar=表达式⼗⼋、数据类型3、新建变量的过程中定义数据类型●gen str3 gender= "male"●list gender in 1/104、变量所占字节过长●drop gender●gen str30 gender= "male"●browse●describe gender●compress gender5、⽇期数据类型：%d dates, which is a count of the number of days elapsed since January 1, 1960。
stata笔记
1.一般检验假设系数为0，t比较大则拒绝假设，认为系数不为0.假设系数为0，P比较小则拒绝假设，认为系数不为0.假设方程不显著，F比较大则拒绝假设，认为方程显著。
2.小样本运用OLS进行估计的前提条件为：（1）线性假定。
即解释变量与被解释变量之间为线性关系。
这一前提可以通过将非线性转换为线性方程来解决。
（2）严格外生性。
即随机扰动项独立于所有解释变量：与解释变量之间所有时候都是正交关系，随机扰动项期望为0。
(工具变量法解决)（3）不存在严格的多重共线性。
一般在现实数据中不会出现，但是设置过多的虚拟变量时，可能会出现这种现象。
Stata可以自动剔除。
（4）扰动项为球型扰动项，即随即扰动项同方差，无自相关性。
3.大样本估计时，一般要求数据在30个以上就可以称为大样本了。
大样本的前提是（1）线性假定（2）渐进独立的平稳过程（3）前定解释变量，即解释变量与同期的扰动项正交。
（4）E（XiXit）为非退化矩阵。
（5）gt为鞅差分序列，且其协方差矩阵为非退化矩阵。
与小样本相比，其不需要严格的外生性和正太随机扰动项的要求。
4.命令稳健标准差回归：reg y x1 x2 x3, robust 回归系数与OLS一样，但标准差存在差异。
如果认为存在异方差，则使用稳健标准差。
使用稳健标准差可以对大样本进行检验。
只要样本容量足够大，在模型出现异方差的情况下，使用稳健标准差时参数估计、假设检验等均可正常进行，即可以很大程度上消除异方差带来的副作用对单个系数进行检验：test lnq=1线性检验：testnl _b[lnpl]=_b[lnq]^25.如果回归模型为非线性，不方便使用OLS,则可以采取最大似然估计法（MLE）,或者非线性最小二乘法（NLS）6.违背经典假设，即存在异方差的情况。
截面数据通常会出现异方差。
因此检验异方差可以：（1）看残差图，但只是直观，可能并不准确。
rvfplot (residual-versus-fitted plot) 与拟合值的散点图rvpplot varname (residual-versus-predictor plot) 与解释变量的散点图扰动项的方差随观测值而变动，表示可能存在异方差。
Stata上机实验笔记
Stata上机实验Stata 统计软件包是目前世界上最著名的统计软件之一，国外将Stata与SAS、SPSS 一起被并称为三大权威软件。
它同时具有数据管理软件、统计分析软件、绘图软件、矩阵计算软件和程序语言的特点，几乎可以完成全部复杂的统计分析工作。
Stata有什么优势？1。
Stata 的命令语句极为简洁明快，易学易记。
2。
强大的帮助信息。
本地帮助 Help 命令名在线帮助Findit 命令名3。
始终处于计量经济学和统计学的最前沿。
许多Stata 程序员会针对计量经济学发展编写一些最新的程序（ADO 文件）， Stata提供了严谨、简练而灵活的程序语句，用户可以编写自己的命令和函数，同时可随时到Stata 网站寻找并下载最新的升级文件。
下载后可以直接使用，也可以自行修改、添加功能。
（例如当前流行的面板单位根和面板门限数据，均可以安装下载使用）不同版本对样本容量、变量个数、矩阵阶数、宏的字符长度等有着不同的限制。
以SE版为例，其最大变量个数为32767，最大字符长度为244字节，最大矩阵阶数为11000（即11000 11000）。
Stata默认值为：最大变量个数为5000，最大矩阵阶数为400，最大内存为10兆。
如果用户需要更多的内存或者更多的变量，可以在命令栏输入如下命令进行扩展。
set maxvar 5000 <最大变量个数5000个。
>set memory 50m <占内存50兆。
>最重要的有三类文件1。
文件名.dta 数据文件2。
文件名.do 命令文件3。
文件名.ado 程序文件如果不加改变，安装时Stata会将系统程序安装到：C：\Program file\stata10 中。
将所用系统自带的一些系统数据、应用程序、帮助文件安装到C：\Program file\stata10\ado\base 中将所有升级程序安装到：C：\Program file\stata10\ado\update 中1。
stata上机实验第四讲
误差项存在自相关：非主对角线上的元素不
为0 。
2 . n 1
.
2
.. ... ... ...
n2
n2 . 2
n 1
考察英国政府如何根据长期利率（r20）的变
化来调整短期利率（rs），数据集为 ukrates.dta （1）做如下回归：rst r 20t 1 t ，其中：
5,000
国产拟合 整体拟合
1。仅截距发生变化。我们以国产还是进口作为虚
拟变量，为了熟悉虚拟变量的产生过程，这里不用 foreign，而是产生一个新的虚拟变量dummy（虚拟 变量的生成要灵活运用gen语句和replace语句）。 gen dummy = 0 replace dummy = 1 if foreign==1 reg price dummy weight 和前面做的回归比较，dummy反映了进口车和国产 车常数项即截距的差异，斜率基本没有变化。
理论分析：加入虚拟变量后，方程变为
price = b0+b1*dummy+b2*weight +u 国产车，即dummy=0时：
方程变为：price = b0+b2*weight+u 进口车，即dummy=1时： 方程变为：price = (b0+b1)+b2*weight+u 结论：截距变化但斜率未变。
rst rst rst 1 Nhomakorabea归方程为：
r 20t 1 r 20t 1 r 20t 2
use ukrates,clear tsset month reg D.rs LD.r20
Stata重点笔记
Stata重点笔记以下命令均采⽤⼩写字母Chapter 1 stata⼊门1.1打开数据use "D:\Stata9\auto.dta", clear ⽤use命令打开数据sysuse auto，clear auto 为系统数据sysuse为打开系统数据的命令1.2获取帮助Help summarize summarize为需要获取帮助对象可以改为其他的需要帮助的对象Findit summarize，net 寻找⽹络帮助summarize为需要获取帮助对象Search summarize ，net 寻找⽹络帮助summarize为需要获取帮助对象显⽰结果Display 5+91.3描述统计（summarize 可简写成sum）Use atuo，clearSummarize price 描述price的观察值个数、平均值、标准差、最⼩值、最⼤值Sum weight summarize可简写成sumSum weight price 同时完成上⾯两步1.4绘图Scatter price weight scatter 为绘制散点图命令Line price weight ，sort line 为绘制折线图命令，sort为排序，绘制折线图前需要先排序1.5⽣成新的数据（generate 可简写成gen）ClearSet obs 1000 设置观测值的组数Gen x=_n _n 为观察值得序号Gen y=x+1001.6控制结果输出显⽰List n1.7设置屏幕滚动Set more off 先设置此项则显⽰时，屏幕不停⽌Set more on 先设置此项则显⽰时，会使显⽰停⽌1.8清除内存中原有内容clear1.9设置⽂件存取路径（cd）Cd d:\stata d:\stata为路径1.10如果想知道当前路径下有哪些⽂件，可以⽤dir 命令来列⽰.dir1.11假设你想在D 盘的根⽬录下创建⼀个新的⽂件夹mydata 来存放数据⽂件，命令为mkdir。
stata上机练习总结
阶
—(1)
验异方差
rvfplot 自 量
predict e,residuals
朗
乘
异方差
验 验 方差
November 10, 2006 7 / 19
验 H0 :
)
White异方差 异方差 常
Econometrics Test
estat imtest,white 验
(
4、
法 predict r gen lagr=l.r
阶
OLS
立 系 系
regress y x1 x2 x3 F R2
方差分 分 判 依据p 最 显 符
验 进 经 验 一 显
显 验 除一 量
进
t
教
席尧
(
)
Econometrics Test
November 10, 2006
6 / 19
3、 立
阶
OLS
立 系 系 判
regress y x1 x2 x3 F R2
November 10, 2006
5 / 19
2、
方法 据 换
据
阶
要
gen lnq=log(q)
replace year=2001 in 23
量 summarize
教
席尧
(
)
Econometrics Test
November 10, 2006
5 / 19
2、
方法 据 换
据
阶
要
gen lnq=log(q)
量 summarize 相 系 矩 矩
pwcorr y x1 x2 x3,obs sig star(5)
线
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
例二：打开production reg lny lnl lnk testnl _b[lnl] * _b[lnk] = 0.25 testnl _b[lnl] * _b[lnk] = 0.5
大样本OLS
大样本OLS经常采用稳健标准差估计(robust) 稳健标准差是指其标准差对于模型中可能存在的异方差或自相关问 题不敏感，基于稳健标准差计算的稳健t统计量仍然渐进分布t分布。因 此，在Stata中利用robust选项可以得到异方差稳健估计量。
回归后预测值的获得 Predict 1。拟合值的获得： predict yhat, xb 或者 predict yhat 2。残差的获得 predict e , residuals 或者 predict e, res
回归的假设检验 Test命令 例一 sysuse auto, clear reg price mpg weight length 1。检验参数的联合显著性 2。分别检验各参数的显著性 3。三个参数对被解释变量的影响相同
如果不加改变，安装时Stata会将系统程序安装到： C：\Program file\stata10 中。 将所用系统自带的一些系统数据、应用程序、帮助文件安装到 C：\Program file\stata10\ado\base 中 将所有升级程序安装到： C：\Program file\stata10\ado\update 中
y = Xb +ε
假设2: X 是满秩的，i.e. rank(X) = k
假设3: 干扰项的条件期望为零（严格外生性）
* E[ε| X] = 0
通过上课的学at = X*b
/* 被解释变量的拟合值*/
e = y - y_hat = y - Xb /* 残差 */
建立回归方程 打开系统文件auto，建立如下方程：
例二： use wage2, clear reg lnwage educ tenure exper expersq 1。教育（educ）和工作时间（tenure）对工资的影响相同。
test educ=tenure 2。工龄（exper）对工资没有影响
test exper 或者 test exper =0 3。检验 educ和 tenure的联合显著性
test educ tenure 或者 test (educ=0) (tenure=0)
例三：生产函数production use production,clear reg lny lnl lnk
test lnl lnk test (lnl=0.8) (lnk=0.2) test lnk+lnl=1
本地帮助 Help 命令名 在线帮助Findit 命令名 3。始终处于计量经济学和统计学的最前沿。许多Stata 程序员会针对 计量经济学发展编写一些最新的程序（ADO 文件）， Stata提供了严 谨、简练而灵活的程序语句，用户可以编写自己的命令和函数，同时可 随时到Stata 网站寻找并下载最新的升级文件。下载后可以直接使用， 也可以自行修改、添加功能。（例如当前流行的面板单位根和面板门限 数据，均可以安装下载使用）
sysuse auto,clear regress price mpg weight foreign
Regress命令详解：
regress depvar [indepvars] [if] [in] [weight] [, options]
1。要求方程省略常数项(自己设置常数项) reg price mpg weight foreign, nocons(hascons) 2。稳健性估计（一般用于大样本OLS） reg price mpg weight foreign, vce(robust) 或者：reg price mpg weight foreign, r 3。设置置信区间（默认95%） reg price mpg weight foreign, level(99) 4。标准化系数 reg price mpg weight foreign, beta 5。部分数据回归 reg price mpg weight length foreign in 1/30 （为什么foreign被drop掉？） reg price mpg weight length if foreign==0
约束回归 定义约束条件 constraint define n 条件 约束回归语句 Cnsreg 被解释变量 解释变量, constraints(条件编号)
例一：打开production cons def 1 lnl+lnk = 1 cnsreg lny lnl lnk, c (1)
例二：sysuse auto,clear cons def 1 price = weight cons def 2 displacement = weight cons def 3 gear_ratio = -foreign cnsreg mpg price weight displacement gear_ratio foreign
回归结果解读 系数/标准误差= t值 P值 系数=0的概率为 p值 在5%的水准上显著不为0 否则和0的差异不显著 95%下限=估计值-t值*标准误差 95%下限=估计值+t值*标准误差 置信区间： 系数在95%的概率下会落在---之间 跨越0，则与0不显著
模型常用的其他形式： 对数 平方项 n次方 指数 交乘项 虽然对函数形式的选择有检验方法，但最好还是从“经济意义”角 度确定。
6。 in range 对数据进行范围筛选 7。 Weight 给数据赋一个权重 8。 , options 命令增加一些可选信息
Stata数据处理 1。Stata的数据格式为 文件名.dta 2。对于系统自带数据文件，任何情况下可以用 sysuse 命令打 开。 3。对于C:\data(或者D:\data)下的文件，可以使用use命令打开。 4。 放在其他目录的文件可以利用设置目录或者“打开” 菜单打 开。
Stata上机实验 Stata 统计软件包是目前世界上最著名的统计软件之一，国外将Stata 与SAS、SPSS 一起被并称为三大权威软件。它同时具有数据管理软件、 统计分析软件、绘图软件、矩阵计算软件和程序语言的特点，几乎可以 完成全部复杂的统计分析工作。 Stata有什么优势？ 1。Stata 的命令语句极为简洁明快，易学易记。 2。强大的帮助信息。
非线性检验：testnl
例一
.
sysuse auto
gen weight2 = weight^2
reg price mpg trunk length weight weight2 foreign
testnl _b[mpg] = 1/_b[weight]
testnl (_b[mpg] = 1/_b[weight]) (_b[trunk] = 1/_b[length])
1。所有的系统自带数据可以利用sysuse命令打开。 2。Use命令只能打开 C:\data 或者 D:\data
中的数据。 3。如果需要打开其他文件夹的数据，必须改变目录（例如，将自己的 数据放入D:\abc）
cd "D:\abc” 或者直接 file------open
特别注意：1。Stata的命令区分大小写。 2。大部分命令可以缩写。
Stata数据类型 数值型：用0、1、2…9 及+、–（正负号）与小数点和( )来表 示。 字符型：字符串变量由字母数字或一些特殊的符号组成（如地名、 住址，职业等等）。 注意：数字不表示大小信息，如电话、身份证号等。 字符型数据 必须要加 “”。 日期型：在Stata中，1960 年1 月1 日被认为是第0 天，因此1959 年12 月31 日为第-1天，表示形式为：jan/10/2001或者10jan2001。
命令： . label define marry 1 “married” 0 “unmarried" . label values married marry
Stata作图
STATA 提供各种曲线类型，包括点（scatter）、线（line）、面
（area），直方图（histogram）、
条形图（bar）、饼图（pie）、函数曲线（function）以及矩阵图
（matrix）等。
同时，对时间序列数据有以ts
开头的一系列特殊命令，如
tsline。还有一类是对双变量的回归拟合图（lfit、qfit 、lowess）
等。
小样本OLS
小样本OLS假设条件较为严格
假设1: 二者之间存在线性关系
y = a0 + a1*x1 + a2*x2 + ... + ak*xk +ε
例如：为变量marrid添加数值标签marry：1=married; 0=Unmarried
菜单：Data->Labels->Label values->Define or modify label values
Data->Labels->Label values->Assign label values to variable
举例：画出Y=X2的曲线图
drop _all set obs 100 gen x = _n gen y = x^2 scatter y x
(drop data from memory) (make 100 observations) (x = 1, 2, 3, .., 100) (y = 2, 4, 9, .., 10000) (make a graph)
use 打开数据文件，一般加clear选型清空内存中现有数据。 sysuse 打开系统数据文件。 describe 描述数据 edit 利用数据编辑器进行数据编辑 list 类似于edit，但只能显示不能修改数据。 display 显示计算结果。经常写为： di summarize 求某个变量的观察值个数、平均值、标准差、最小值和最 大值。经常写为：sum scatter 生成两个变量的散点图。 set obs 定义样本个数(使用前一定要用drop或者clear命令清空) generate 建立新变量并赋值。经常写为gen