-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path04_dmm_child.R
123 lines (83 loc) · 4.03 KB
/
04_dmm_child.R
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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
#===============================================================================
# 29/07/2024
# Diagonal Mobility Model on Children
# Educational Assortative Mating and Subjective Well-Being
# Yanwen Wang, [email protected]
#===============================================================================
######################################################################
# Loading packages
######################################################################
library(tidyverse)
library(haven)
library(gnm)
library(lmtest)
library(stargazer)
library(jtools)
######################################################################
# Loading data
######################################################################
load("to-be-filled")
# DMM All Children --------------------------------------------------------
# * Baseline --------------------------------------------------------------
set.seed(123)
mod.base <- gnm(RSATISFY ~ -1 + RGENDER + RAGE + RMSTAT + RPARTY + RHUKOU + RMINORITY + RRESIDE + RSINGLE +
Dref(PEDUA, REDUA),
data = child)
print(summary(mod.base), digits = 4)
print(DrefWeights(mod.base), digits = 4)
# * + Single-Child Weight Interaction --------------------------------------
mod.singleweight <- gnm(RSATISFY ~ -1 + RGENDER + RAGE + RMSTAT + RPARTY + RHUKOU + RMINORITY + RRESIDE + RSINGLE +
Dref(PEDUA, REDUA,
delta = ~ 1 + RSINGLE),
data = child)
print(summary(mod.singleweight), digits = 4)
print(DrefWeights(mod.singleweight), digits = 4)
lrtest(mod.singleweight, mod.base)
# * + Mobility Indicators -------------------------------------------------
# Main model with upward/downward mobility
mod.mob <- gnm(RSATISFY ~ -1 + upa + downa + RGENDER + RAGE + RMSTAT + RPARTY + RHUKOU + RMINORITY + RRESIDE + RSINGLE +
Dref(PEDUA, REDUA),
data = child)
print(summary(mod.mob), digits = 4)
print(DrefWeights(mod.mob), digits = 4)
lrtest(mod.mob, mod.base)
# Model with interaction with only-child status
mod.mob.single <- gnm(RSATISFY ~ -1 + RSINGLE*upa + RSINGLE*downa + RGENDER + RAGE + RMSTAT + RPARTY + RHUKOU + RMINORITY + RRESIDE +
Dref(PEDUA, REDUA),
data = child)
print(summary(mod.mob.single), digits = 4)
print(DrefWeights(mod.mob.single), digits = 4)
lrtest(mod.mob.single, mod.base)
# DMM Only Children -------------------------------------------------------
child <- child %>%
filter(RSINGLE == 1)
# * Baseline --------------------------------------------------------------
set.seed(123)
mod.base <- gnm(RSATISFY ~ -1 + RGENDER + RAGE + RMSTAT + RPARTY + RHUKOU + RMINORITY + RRESIDE +
Dref(PEDUA, REDUA),
data = child)
print(summary(mod.base), digits = 4)
print(DrefWeights(mod.base), digits = 4)
# * + Gender Weight Interaction -------------------------------------------
mod.genderweight <- gnm(RSATISFY ~ -1 + RGENDER + RAGE + RMSTAT + RPARTY + RHUKOU + RMINORITY + RRESIDE +
Dref(PEDUA, REDUA,
delta = ~ 1 + RGENDER),
data = child)
print(summary(mod.genderweight), digits = 4)
print(DrefWeights(mod.genderweight), digits = 4)
lrtest(mod.genderweight, mod.base)
# * + Mobility Indicators -------------------------------------------------
# Main model with upward/downward mobility
mod.mob <- gnm(RSATISFY ~ -1 + upa + downa + RGENDER + RAGE + RMSTAT + RPARTY + RHUKOU + RMINORITY + RRESIDE +
Dref(PEDUA, REDUA),
data = child)
print(summary(mod.mob), digits = 4)
print(DrefWeights(mod.mob), digits = 4)
lrtest(mod.mob, mod.base)
# Model with gender interaction
mod.mob.gender <- gnm(RSATISFY ~ -1 + RGENDER*upa + RGENDER*downa + RAGE + RMSTAT + RPARTY + RHUKOU + RMINORITY + RRESIDE +
Dref(PEDUA, REDUA),
data = child)
print(summary(mod.mob.gender), digits = 4)
print(DrefWeights(mod.mob.gender), digits = 4)
lrtest(mod.mob.gender, mod.base)