常用的27個Stata命令
公司主營業(yè)務(wù):成都網(wǎng)站制作、成都做網(wǎng)站、移動網(wǎng)站開發(fā)等業(yè)務(wù)。幫助企業(yè)客戶真正實現(xiàn)互聯(lián)網(wǎng)宣傳,提高企業(yè)的競爭能力。創(chuàng)新互聯(lián)是一支青春激揚、勤奮敬業(yè)、活力青春激揚、勤奮敬業(yè)、活力澎湃、和諧高效的團(tuán)隊。公司秉承以“開放、自由、嚴(yán)謹(jǐn)、自律”為核心的企業(yè)文化,感謝他們對我們的高要求,感謝他們從不同領(lǐng)域給我們帶來的挑戰(zhàn),讓我們激情的團(tuán)隊有機(jī)會用頭腦與智慧不斷的給客戶帶來驚喜。創(chuàng)新互聯(lián)推出膠州免費做網(wǎng)站回饋大家。
【1】直接導(dǎo)入csv格式數(shù)據(jù)
insheet using name.csv, clear
【2】修改變量長度
format var %20.2g
【3】刪除重復(fù)值
sort var1 var2
duplicatesdrop var1 var2, force
【4】數(shù)據(jù)合并
use data1, clear
merge m:m var1var2 using data2
drop if _merge==2
drop if _merge==1
drop _merge
【5】生成一期滯后項
tsset stkcd accper
gen newvarname=L.varname
【6】將文字轉(zhuǎn)化為數(shù)字變量
genBigN=0
replaceBigN=1 if strmatch(dadtunit,"普華永道*")
【7】刪除有缺失值的記錄
egen mis=rowmiss(_all)
drop if mis
drop mis
【8】行業(yè)劃分
clonevarsic2=ind
order stkcd accper sic2
replace sic2=substr(sic2,1,1) if substr(sic2,1,1)!=”C”
replace sic2=substr(sic2,1,2) if substr(sic2,1,1)==”C”
tabulate sic2 accper
【9】日期只保留年份
drop if substr( reptdt ,6,2)!="12"
replace reptdt=substr(reptdt,1,4)
gen accper=real(reptdt)
【10】數(shù)據(jù)分列
split date ,parse(-) destring ignor("-")
【11】求兩個日期之間的間隔天數(shù)
g td=date(trading_date,"YMD")
g ed=date(eventdate,"YMD")
form td ed %td
g d=ed-td
【12】生成行業(yè)、年份啞變量
tab year, gen(year)
tab industry, gen(industry)
【13】對數(shù)據(jù)進(jìn)行Winsorize處理
findit winsor2
winsor2 varname, replace cut(1 99)
【14】描述性統(tǒng)計
tabstat var1var2, stat(n min mean median p25 p75 max sd),? if groupvar==0 or 1
logout, save(name) word replace: tabstat var, stat(n min mean p50 max sd) col(stat)f(%9.2g)
【15】兩變量列聯(lián)表
tabulate var1 var2, row chi2 taub gamma
【16】兩樣本間的均值T檢驗
ttest var, by(groupvar)
【17】兩樣本中位數(shù)Z檢驗
ranksum var, by(groupvar)
【18】Pearson/Spearman系數(shù)
spearmanx
n mata
x=st_data(.,"x")
c=correlation(x)
n=rows(c)
b=strofreal(lowertriangle(c)+uppertriangle(st_matrix("r(Rho)")),"%9.3f")
p=st_matrix("r(P)")
for (i=2; i<=n; i++) {
for (j=1; j<=i-1; j++) {
p[i,j]=2*ttail(rows(x)-2,abs(c[i,j]/sqrt((1-c[i,j]^2)/(rows(x)-2))))
b[i,j]=b[i,j]+(p[i,j]<0.01?"":(p[i,j]<0.05?"":(p[i,j]<0.1?"":"")))
b[j,i]=b[j,i]+(p[j,i]<0.01?"":(p[j,i]<0.05?"":(p[j,i]<0.1?"":"")))
}
}
c=editvalue(b, "2.000", "1")
c
end
直接導(dǎo)出結(jié)果
logout, save(pw) word replace:pwcorr_avars, star1(0.01) star5(0.05) star10(0.1)
【19】按年度按中位數(shù)分組
方法一
bysort year: egen g=xtile(var), n(2)
方法二
bys accper: cumul icindex, g(g) eq
levelsof accper, local(id)
display "`r(levels)'"
local cut1 = 1/2
foreach x of local id {
recode g (min/cut1'=0)(
cut1'/max=1) if accper==`x'
}
分三組
bys accper:cumul icindex, g(g) eq
levelsof accper, local(id)
display "`r(levels)'"
local cut1 = 1/3
local cut2 = 2/3
foreach x of local id {
recode g (min/cut1'=1)(
cut1'/cut2'=2)(
cut2'/max=3)if accper==`x'
}
【20】輸出回歸結(jié)果
安裝
ssc install estout, replace
單個回歸
reg
esttab using name.rtf, compress nogap r2 ar2 star( 0.1 0.05 0.01)
多個回歸一起
reg
est store m1
reg
est store m2
esttab m1 m2 using name.rtf, compress nogap r2 ar2 star( 0.1 0.05 0.01)
連續(xù)運行tobit模型結(jié)果導(dǎo)出:
esttab m1 m2, b(%9.4f) t scalars(N ll Fchi2 type), using name.rtf, compress nogap
連續(xù)運行OLS模型結(jié)果導(dǎo)出:
esttab m1m2, b(%9.4f) tscalars(N r2 F p), using name.rtf, compress nogap
【21】異方差檢驗及處理
檢驗:懷特檢驗
ssc install whitetst
reg
estat imtest, white
處理:“OLS+穩(wěn)健標(biāo)準(zhǔn)差”
reg y x1 x2 x3, robust
【22】DW檢驗(序列相關(guān)性一階)
gen id=_n
tsset id
estat dwatson
【23】多重共線性
reg y x1 x2 x3
vif
【24】是否遺漏高次項
例如,檢驗y對x的線性回歸有沒有遺漏高次項
reg y x
estat ovtest
或者estat ovtest, rhs
【25】逐步回歸
stepwise, pe(0.1): reg y x
【26】Maddala(1983)兩階段處理效應(yīng)模型
treatreg yx1-xn, tr(z=w1-wm)two
【27】Justified Jones Model
statsby _b, by(ind accper)saving(*.dta,replace):reg yx, noconstant
merge m:m indaccper using *.dta
gen yhat=y-_b*x
當(dāng)前名稱:27個常用stata命令(2)
文章網(wǎng)址:http://chinadenli.net/article20/pigijo.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供響應(yīng)式網(wǎng)站、用戶體驗、營銷型網(wǎng)站建設(shè)、App設(shè)計、外貿(mào)網(wǎng)站建設(shè)、定制網(wǎng)站
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請盡快告知,我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時需注明來源: 創(chuàng)新互聯(lián)