ん，家族社会学会まであと9日しかないの？

Stata でイベントヒストリー分析．イベントヒストリー用のコマンドとかあるんだろうけど，よくわからないのでわりとアナログにやってます．

* JGSS-2009LCSの.csvを用意する

set more off

* 従属変数をつくる

* 問57で2「同居期間」と9無回答をドロップ
tab ptlvtgl
drop if ptlvtgl > 1
tab ptlvtgl

* 1回目の同居期間終了年月に2「現在まで」と9無回答をドロップ
tab plv01ong
drop if plv01ong > 1
tab plv01ong
* ここで2106人

* センサリングの処理
* 左センサー欠損
* 親との同居期間の終了年（俺的離家）に
* 具体的に答えてない9999をドロップ（4人）
tab plv01eny
drop if plv01eny == 9999
tab plv01eny
* ここで2102人

*右センサー欠損
* 親との2回目の同居期間の開始年（俺的帰家）に
* 具体的に答えてない9999をドロップ（13人）
tab plv02sty
drop if plv02sty == 9999
tab plv02sty
* ここで2089人

* 帰家イベント生起ダミーを作る
gen returning_home = plv02
recode returning_home 2 = 0
tab returning_home
* ここで464人が帰家経験者

* 生存時間を作る
* リスク終了年 – リスク開始年
tab plv02sty
gen lisk_end = plv02sty
* 2回目の同居開始8888非該当のひとたちは
* 調査時点においても「実家に帰ってきてない」人たちなので，
* この人たちを右センサーする
recode lisk_end 8888 = 2009
tab lisk_end
gen survival_time = lisk_end – plv01eny
tab survival_time
* １ずつ足す
replace survival_time = survival_time + 1
tab survival_time

* 記述的分析
stset survival_time, failure(returning_home == 1) id(id)
stsum
recode opffix15 9 = .
recode tp5loc15 9 = .
label define sex 1 “male” 2 “female”
label values sexa sex

sts graph, by(sex)
sts graph, by(edu)
sts graph, by(opffix15)
sts graph, by(tp5loc15)

* 離家理由をつくる
* 進学離家 高等教育在学開始年と離家年が一致していれば1
tab plv01eny
tab sch01sty
gen leaving_home_reason_school = 0
recode leaving_home_reason_school 0 = 1 if plv01eny == sch01sty
tab leaving_home_reason_school

* 就職離家 働きはじめた時期と離家年が一致していれば1
tab plv01eny
tab jb01sty
gen leaving_home_reason_job = 0
recode leaving_home_reason_job 0 = 1 if plv01eny == jb01sty | plv01eny == jb02sty | plv01eny == jb03sty
tab leaving_home_reason_job

tab leaving_home_reason_job plv01eny
* 闇を感じる

* 結婚離家
tab plv01eny
tab mrg01sty
gen leaving_home_reason_marriage = 0
recode leaving_home_reason_marriage 0 = 1 if plv01eny == mrg01sty
tab leaving_home_reason_marriage

gen lhr = 0
recode lhr 0 = 1 if leaving_home_reason_school == 1
recode lhr 0 = 2 if leaving_home_reason_job == 1
recode lhr 0 = 3 if leaving_home_reason_marriage == 1
sts graph, by(lhr)