这个我一般是用循环做得,还没找到更好的办法。然后将结果存在一个矩阵中,最后再把这个矩阵转化为数据集。假设的你的时间变量是从1到20,那么命令是:
mat tteststat=J(20,5,.)
forvalues i=1/20 {
qui ttest RET if YEAR=`i',by(TREAT)
tteststat[`i',1]=r(mu_1)
tteststat[`i',2]=r(mu_2)
tteststat[`i',3]=r(se)
tteststat[`i',4]=r(t)
tteststat[`i',5]=r(p)
}
svmat tteststat,names(col)
brow c?