stata.docx

insheet using d:研究生statalianxicha3.csvrename v1 yearrename v2 ybrowse*查看数据browse if year20000 high, 15000 medium, 小于=小于等于10000*生成虚拟变量的快捷方法：gen p=price10000(price10000取1，否则取0)sort pricegen pp=_n*生成一个表示价格高低的序列help function*显示全部stata函数drop p pp ww www wdis round(_pi,0.00001)*显示pi四舍五入到0.00001dis normal(1.96)*返回标准正态分布t=1.96的累积概率分布dis invnormal(0.975)*返回标准正态分布累积分布临界概率对应的临界值dis Ftail(2,1500,3.6979)*返回自由度为2样本容量为1500的F分布的F=3.6979的累积概率dis invFtail(2,1500,0.025)*返回自由度为2样本容量为1500的F分布的累积概率为0.025的F值lookfor p*寻找含有p的变量sort pricelist make-mpg in 1/10,cleanlist make-mpg if price0histogram inc if inc0,normalladder inc if inc0boxcox inc if inc0bcskew0 incx=inc if inc0histogram incx,normalgladder inc if inc0*回归之后发现age对inc不显著，不符合常识，经过作图发现因变量不符合/*/正态分布假定，故而想办法将因变量转化为符合正态分布的变量，首先用/*/常见的幂指数进行试验，观察有无通过检验的形式，结果是没有（ladder）/*/然后用boxcox命令强制进行幂指数转换，计算出最佳幂指数为0.096/*/再用bcskew0命令计算出可以令偏度等于零的boxcox转换，并生成新的变量/*/incx，对incx进行画图正态比较发现结果较满意。gladder命令是用作图/*/显示调整结果ta edu,nolabelgen edug=recode(edu,2,3,4,6)ta edug,gen(sb)rename sb1 primaryrename sb2 juniorrename sb3 seniorrename sb4 collegegen incn=log(inc) if inc0reg incn age male junior senior college*一种生成多组别虚拟变量的方法,在ta命令后直接加入gen(*)选项，即可直接/*/根据ta分组创建等量的变量，每个变量都以*命名。adjust age=40 male=1,by(edug)*回归后的点预测！xi:reg incn age male g*最方便的设定虚拟变量用于拟合回归的方法，直接在回归命令前加xi:，然后在/*/需要转换为虚拟变量的类别变量前加i.，stata将自动生成虚拟变量，并且/*/自动选取n-1个虚拟变量（放弃一个）进入方程。use cgss2003,clearxi:reg incn i.ageg i.party*i.state i.male*eduyr*带交乘项的虚拟变量回归拟合（自动生成）。sw reg incn male age edu party state,pr(.05)sw reg incn male age edu party state,pe(.05)*逐步回归（stepwise）pr(probility of removing)向后逐步回归/*/pe（begin with empty model）向前逐步回归use hdr2004x,clearhist hdi,freq start(.2) width(.1) ylabel(0 (5) 35)/*/ xlabel(.2 (.1) 1.0) gap(8) normal addlabels /*/title(the distribution of human development index)*直方图绘制命令gen hdigroup=recode(hdi,.5,.7,.8,.9,1)ta hdigroup,gen(h)replace h1=. if h1=0replace h2=. if h2=0replace h3=. if h3=0replace h4=. if h4=0replace h5=. if h5=0graph pie h1-h5,legend(off) plabel(_all name)/*/ title(Human Development Index) /*/subtitle(by hdi group) /*/note(source:2004 UNDP Human Development Report)*饼图绘制命令use chnhf2000,cleargraph bar inc if occ1,over(edu) ylabel(0 (600) 2400) /*/ title(Monthly Average Income) /*/subtitle(by education) ytitle(Ave. Inc.) /*/note(Source: 2000 CHNHF Project Subsample)*简单条形图绘制命令graph bar inc if occ1,over(edu,label(alt) /*/over(male) ylabel(0 (600) 2400) /*/title(Monthly Average Income) /*/subtitle(by education and gender) /*/ytitle(Ave. Inc.) /*/note(Source: 2000 CHNHF Project Subsample)graph bar inc if occ1,over(male,gap(8) label(alt)/*/ over(edu) ylabel(0 (600) 2400) /*/title(Monthly Average Income) /*/subtitle(by education and gender) ytitle(Ave. Inc.) /*/note(Source: 2000 CHNHF Project Subsample)*绘制根据多个变量分组的单变量条形图use netview,clearreplace day=0 if day=7label define day 0 sunday 1 monday 2 tuesday,modifylabel val day daylist day view1-view3 if day3,cleangraph bar view1-view3 if day3|day=0,over(day) bargap(-25)/*/ legend(posi(1) rows(3) ring(0) lab(1 Department1) /*/lab(2 Department2) lab(3 Department3) /*/title(Internet Visits) subtitle(by development and day) /*/ytitle(Visits) note(Source: Undisclosed Sources)graph bar (sum) view1-view3 if day1,over(occ,sort(inc) descending) by(male) nofill /*/ylab(0 (2000) 6000) yline(943) title(Monthly Average Income) /*/ytitle(Ave. Inc.) note(Source: 2000 CHNHF Project Subsample)*横向条形图use bankchina3q,clearformat changex %3.1fgraph twoway bar changex day if month=8,yline(0) ylab(-5.0 (5) 10.0) /*/xlab(1 (1) 31,alter) ymtick(-5 (1) 10) title(Bank Of China) /*/subtitle(August 2007) ytitle(Change in %) xtitle(Date) /*/note(Source: /,extracted on 10/04/2007)*多方向条形图（股票走势图）use hdr2004x,clearscatter enrollf literacy,msymbol(oh)count if literacy95&contin=3scatter enrollf literacy if literacy95&contin=3,msymbol(oh)scatter enrollf literacy if literacy95&contin=3,msymbol(oh) jitter(5)*绘制散点图，处理重合点scatter hdi lifeexp if contin=1,mlabel(nation)gen pos=3replace pos=9 if rank=147replace pos=9 if rank=168replace pos=2 if rank=175replace pos=8 if rank=148replace pos=9 if rank=152scatter hdi lifeexp if contin=1,xlabel(25 (5) 75) mlabel(nation) mlabv(pos)*设置数据标签的位置gen gdpn=log(gdp)lab var gdpn gdp per capital in log formgraph matrix hdi gdpn literacy lifeexp,subtitle(Summary of HDR2004 data) /*/note(Source: UNDP Human Development Report 2004)*矩阵散点图use lifestats,clearscatter birth year ,connect(l) msymbol(o) subtitle(Birth Rate In China) /*/note(Source: from Statistical Yearbook 2006)*带连线的散点图line birth year,sort(year)line birth increase death year,sort(year) lpattern(solid dash_dot dash) /*/yaxis(1 2) xaxis(1 2) ylabel(0 10 (5) 25) ylabel(0 (2) 10,axis(2) /*/xlab(,axis(2) xtitle(,axis(2) legend(pos(2) ring(0) rows(3) size(small) /*/subtitle(Life Statistics in China,1980-2005) /*/note(Source: from Statistical Yearbook 2006)*直线图use hdr2004x,clearscatter hdi literacy|lfit hdi literacyscatter hdi literacy|lfitci hdi literacy*直线拟合图形gen gdpn=log(gdp)scatter internetk gdpn|qfit internetk gdpnscatter internetk gdpn|qfitci internetk gdpn*二次曲线拟合图形use bankchina3q,cleartwoway spike change datetm,yline(0)twoway dropline change datetm,yline(0)use lifestats,cleartwoway spike birth year,base(18) yline(18) ylabel(10 (2) 24)*长钉图use lifestats,cleartwoway area birth death year,ytitle(percent)use bankchina3q,cleartwoway area change datetm,xtitle(Date) ytitle(Price Change)*面积图use lifestats,cleartwoway rarea birth death year,ytitle(percent)*面积区间图use bankchina3q,cleartwoway rbar high low datetm|scatter close datetm,legend(pos(5) ring(0) row(2)*条形区间图use hdr2004x,cleartwoway lfitci hdi literacy|scatter hdi literacy|,legend(order(3 1 2) /*/lab(3 Observed) lab(2 Estd Line) rows(3) position(10) ring(0) ) /*/text(80 .2 I Love You Bei!)*图例的设计（legend）use bankchina3q,clearreplace volume=volume/1000000twoway rbar high low datetm|line close datetm|bar volume datetm, yaxis(2)|,/*/yscale(axis(1) range(4 6.5) yscale(axis(2) range(0 12) ylabel(4.5 (.5) 6.5) /*/ylabel(0 (2) 6,axis(2) legend(pos(3) ring(0) rows(3) xtitle(Date) /*/ytitle(Stock Volume in million axis(2) title(Bank of China Stock Trend) /*/subtitle(3rd Quarter,2007) /*/note(Source: /,exteracted on 10/04/2007)*双坐标轴的用法（axis）graph save bcstock,replacegraph dirgraph use bcstockgraph export bcstock.png,replace*图形的保存、调用与输出。use append1,clearlistappend using append2save data,replace*纵向合并数据（变量数目不变，样本数增加