* https://bbs.pinggu.org/thread-437023-1-1.html
* 生成模拟数据 N=400,T=20
clear
*ssc install xtarsim, replace
xtarsim y x eta, n(400) t(20) gamma(0) beta(0.8) rho(0.2) ///
one(corr 1) sn(9) seed(1234)
tsset ivar tvar
xtdes
* 5年一平均
egen tt = seq(), from(1) to(4) block(5)
foreach v of varlist y x{
egen `v'_m = mean(`v'), by(tt)
}
keep *_m tt ivar
duplicates drop
tsset ivar tt
xtdes
*-------------results----------------
. * https://bbs.pinggu.org/thread-437023-1-1.html
. * 生成模拟数据 N=400,T=20
. clear
. ssc install xtarsim, replace
. xtarsim y x eta, n(400) t(20) gamma(0) beta(0.8) rho(0.2) ///
> one(corr 1) sn(9) seed(1234)
. tsset ivar tvar
panel variable: ivar (strongly balanced)
time variable: tvar, 1 to 20
delta: 1 unit
. xtdes
ivar: 1, 2, ..., 400 n = 400
tvar: 1, 2, ..., 20 T = 20
Delta(tvar) = 1 unit
Span(tvar) = 20 periods
(ivar*tvar uniquely identifies each observation)
Distribution of T_i: min 5% 25% 50% 75% 95% max
20 20 20 20 20 20 20
Freq. Percent Cum. | Pattern
---------------------------+----------------------
400 100.00 100.00 | 11111111111111111111
---------------------------+----------------------
400 100.00 | XXXXXXXXXXXXXXXXXXXX
. * 5年一平均
. egen tt = seq(), from(1) to(4) block(5)
. foreach v of varlist y x{
2. egen `v'_m = mean(`v'), by(tt)
3. }. keep *_m tt ivar
. duplicates drop
Duplicates in terms of all variables
(6400 observations deleted)
. tsset ivar tt
panel variable: ivar (strongly balanced)
time variable: tt, 1 to 4
delta: 1 unit
. xtdes
ivar: 1, 2, ..., 400 n = 400
tt: 1, 2, ..., 4 T = 4
Delta(tt) = 1 unit
Span(tt) = 4 periods
(ivar*tt uniquely identifies each observation)
Distribution of T_i: min 5% 25% 50% 75% 95% max
4 4 4 4 4 4 4
Freq. Percent Cum. | Pattern
---------------------------+---------
400 100.00 100.00 | 1111
---------------------------+---------
400 100.00 | XXXX

[em05]
[此贴子已经被作者于2009-3-31 14:46:31编辑过]