cd /Victor/stata
//二
//1a
sysuse "autodta",clear
foreach num of numlist 1/74{
export excel using "auto`num'xlsx" in `num'
}
//1b
foreach num of numlist 1/74{
import excel "auto`num'", clear
save "auto`num'dta", replace
}
//将excel转换成dta格式
use "auto1dta",clear
foreach num of numlist 2/70{
append using "auto`num'dta"
}
save "auto_newdta", replace
//1c
foreach num of numlist 1/74{
erase "auto`num'xlsx"
}
//2(因为我的stata软件无法识别跨行命令,只能一长串写在一行)
use "auto_newdta",clear
twoway(scatter weight mpg if foreign == 0, mcolor(navy blue) msymbol(x)) (scatter weight mpg if foreign == 1, mcolor(red)) , legend(on order(1 "foreign" 2 "domestic") ring(0)) title("The Relatonship between Car Weight and MPG") ylabel(1000 2000 3000(2000)5000) xlabel(10(5)35)
//3
sysuse "nlsw88dta",clear
tabstat wage hours ttl_exp tenure ,by(union) stat(n mean sd min median) c(s)
//4
sysuse "nlsw88dta",clear
order wage hours ttl_exp tenure collgrad
foreach v of varlist wage-collgrad{
ttest `v',by(union)
}
//5
sysuse "nlsw88dta",clear
reg wage age race industry occupation ttl_exp hours
ssc install logout
logout, save(result_2017302330037xlsx) excel dec(3) replace:reg wage age race industry occupation ttl_exp hours
//三
sysuse "firmdta",clear
sort ID date ipc patent_ID
//2
replace firm_name = subinstr(firm_name, "有限责任公司"," ",)
replace firm_name = subinstr(firm_name, "集团股份有限公司"," ",)
replace firm_name = subinstr(firm_name, "股份有限公司"," ",)
//3
//计算每个企业每年专利数量总和
use "firmdta",clear
encode patent_ID ,gen(patent_ID1)
bys ID date : egen num_patent=count(patent_ID1)
//计算不重复的领域数
duplicates drop ID date ipc,force
bys ID date: egen ipc_sum=count(ipc)
欢迎分享,转载请注明来源:品搜搜测评网