Theme NexT works best with JavaScript enabled
1
2
3
4
5
cd "F:\General Database of Social Science\Courses\Courses Materials\波士顿学院\纵贯调查分析\数据"
use hrs_hours,clear
reshape long r@workhours80 r@poorhealth r@married r@totalpar r@siblog h@childlg r@allparhelptw, i(hhid pn) j(wave)
tab wave
1
2
3
4
5
preserve
keep if wave < 3
xtset hhidpn wave
xtdes
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
egen miss = rowmiss(rworkhours80 rpoorhealth rmarried rtotalpar rsiblog hchildlg rallparhelptw)
tab miss
drop if miss == 7
miss | Freq. Percent Cum.
------------+-----------------------------------
0 | 11,327 85.93 85.93
1 | 1,017 7.72 93.64
2 | 115 0.87 94.52
3 | 90 0.68 95.20
4 | 7 0.05 95.25
5 | 3 0.02 95.27
6 | 364 2.76 98.04
7 | 259 1.96 100.00
------------+-----------------------------------
Total | 13,182 100.00
删除全部变量都缺失的个体
根据公式
计算两期相关系数
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
drop if miss == 7
xtset hhidpn wave
xtdes
Freq. Percent Cum. | Pattern
---------------------------+---------
6332 96.07 96.07 | 11
259 3.93 100.00 | 1.
---------------------------+---------
6591 100.00 | XX
xtsum rworkhours80
Variable | Mean Std. Dev. Min Max | Observations
-----------------+--------------------------------------------+----------------
rwork~80 overall | 29.53971 22.79859 0 80 | N = 12477
between | 21.33473 0 80 | n = 6580
within | 8.392351 -10.46029 69.53971 | T-bar = 1.8962
di 21.33473^2 / (21.33473^2 + 8.392351^2)
.86599838
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
xttab rmarried
Overall Between Within
rmarried | Freq. Percent Freq. Percent Percent
----------+-----------------------------------------------------
0 | 2662 21.20 1532 23.24 92.13
1 | 9895 78.80 5300 80.41 97.73
----------+-----------------------------------------------------
Total | 12557 100.00 6832 103.66 96.47
(n = 6591)
注意：
Between 中 80.41%两期内至少有一期结婚， 23.24至少有一期未婚。
within 中 97.73%表示对于两期任意一期结婚的人而言，两期平均的结婚概率
92.13%表示对于两期任意一期未婚的人而言，两期平均的未婚概率
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
xttrans rmarried, freq
| rmarried
rmarried | 0 1 | Total
-----------+----------------------+----------
0 | 1,130 87 | 1,217
| 92.85 7.15 | 100.00
-----------+----------------------+----------
1 | 154 4,595 | 4,749
| 3.24 96.76 | 100.00
-----------+----------------------+----------
Total | 1,284 4,682 | 5,966
| 21.52 78.48 | 100.00
di (1532*.9213+5300*.9773)/(1532+5300)
0.96474262
restore
1
2
3
4
5
6
egen miss=rowmiss( rworkhours80 rpoorhealth rmarried rtotalpar rsiblog hchildlg rallparhelptw)
drop if miss==7
(4,643 observations deleted)
xtset hhidpn wave
save hrs_hours_long.dta
呈现9期数据的模式类型
有5540个的个体9期数据均无缺失
有154个个案从第3期开始缺失
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
for var rworkhours80 rpoorhealth rmarried rtotalpar rsiblog hchildlg rallparhelptw:xtsum X
for var rpoorhealth rmarried : xttab X
xtdes
Freq. Percent Cum. | Pattern
---------------------------+-----------
5540 84.05 84.05 | 111111111
154 2.34 86.39 | 11.......
137 2.08 88.47 | 1........
84 1.27 89.74 | 1111.....
81 1.23 90.97 | 11111....
73 1.11 92.08 | 11111111.
69 1.05 93.13 | 111......
55 0.83 93.96 | 1111111..
49 0.74 94.70 | 111111...
349 5.30 100.00 | (other patterns)
---------------------------+-----------
6591 100.00 | XXXXXXXXX
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
use hrs_hours,clear
reshape long r@workhours80 r@poorhealth r@married r@totalpar r@siblog h@childlg r@allparhelptw, i(hhid pn) j(wave)
sqset rmarried hhidpn wave
sqtab,ranks(1/5)
Sequence-Pa |
ttern | Freq. Percent Cum.
------------+-----------------------------------
1:9 | 2,283 67.15 67.15
0:9 | 508 14.94 82.09
1:3 .:6 | 209 6.15 88.24
1 .:8 | 206 6.06 94.29
1:2 .:7 | 194 5.71 100.00
------------+-----------------------------------
Total | 3,400 100.00
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
preserve
keep if wave < 3
mean rworkhours80 rpoorhealth rmarried rtotalpar rsiblog hchildlg rallparhelptw,over(wave)
---------------------------------------------------------------
Over | Mean Std. Err. [95% Conf. Interval]
--------------+------------------------------------------------
rworkhours80 |
1 | 30.61686 .2907706 30.0469 31.18682
2 | 28.20681 .3175137 27.58442 28.82919
--------------+------------------------------------------------
rpoorhealth |
1 | .1954613 .0051042 .1854562 .2054665
2 | .1988658 .0054884 .1881076 .209624
--------------+------------------------------------------------
rmarried |
1 | .8189498 .0049563 .8092347 .8286649
2 | .8200378 .0052823 .8096836 .830392
--------------+------------------------------------------------
rtotalpar |
1 | 1.606427 .0098604 1.587099 1.625755
2 | 1.405104 .0111955 1.383159 1.427049
--------------+------------------------------------------------
rsiblog |
1 | 1.673055 .0081439 1.657091 1.689018
2 | 1.680897 .0086198 1.664001 1.697793
--------------+------------------------------------------------
hchildlg |
1 | 1.107862 .0070198 1.094102 1.121622
2 | 1.118057 .0074358 1.103481 1.132632
--------------+------------------------------------------------
rallparhelptw |
1 | .6292256 .0373714 .5559711 .70248
2 | 1.228155 .0464493 1.137107 1.319204
---------------------------------------------------------------
restore
1
2
3
4
egen mean_rwork = mean(rworkhours80), by(wave)
line mean_rwork wave, ytitle("Average hours of paid work") xtitle("Wave")
twoway scatter mean_rwork wave, ytitle("Average hours of paid work") xtitle("Wave")
计算两期相关系数