



Marginal effect of dummy variable: effect on Y of going from 0 to 1
Marginal effect of continuous variable: effect on Y of a 1 unit change in X
Example: Consider the gender pay gap again.
Example: Consider the gender pay gap again.
Example: Consider the gender pay gap again.
Example: Consider the gender pay gap again.
^wagesi=β0+β1Genderi+β2Experiencei
Depending on the types of variables, there are 3 possible types of interaction effects
We will look at each in turn
Depending on the types of variables, there are 3 possible types of interaction effects
We will look at each in turn
Depending on the types of variables, there are 3 possible types of interaction effects
We will look at each in turn
Depending on the types of variables, there are 3 possible types of interaction effects
We will look at each in turn


Yi=β0+β1Xi+β2Di+β3(Xi×Di) where Di={0,1}
Yi=β0+β1Xi+β2Di+β3(Xi×Di) where Di={0,1}
Yi=β0+β1Xi+β2Di+β3(Xi×Di) where Di={0,1}
β3 estimates the interaction effect between Xi and Di on Yi
What do the different coefficients (β)’s tell us?
Yi=β0+β1Xi+β2Di+β3Xi×Di
Yi=β0+β1Xi+β2Di+β3Xi×Di
^Yi=^β0+^β1Xi+^β2(0)+^β3Xi×(0)^Yi=^β0+^β1Xi
Yi=β0+β1Xi+β2Di+β3Xi×Di
^Yi=^β0+^β1Xi+^β2(0)+^β3Xi×(0)^Yi=^β0+^β1Xi
^Yi=^β0+^β1Xi+^β2(1)+^β3Xi×(1)^Yi=(^β0+^β2)+(^β1+^β3)Xi
Yi=β0+β1Xi+β2Di+β3Xi×Di
^Yi=^β0+^β1Xi+^β2(0)+^β3Xi×(0)^Yi=^β0+^β1Xi
^Yi=^β0+^β1Xi+^β2(1)+^β3Xi×(1)^Yi=(^β0+^β2)+(^β1+^β3)Xi

Yi=^β0+^β1Xi
Yi=(^β0+^β2)+(^β1+^β3)Xi
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi+ΔYi=β0+β1(Xi+ΔXi)β2Di+β3((Xi+ΔXi)Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi+ΔYi=β0+β1(Xi+ΔXi)β2Di+β3((Xi+ΔXi)Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi+ΔYi=β0+β1(Xi+ΔXi)β2Di+β3((Xi+ΔXi)Di)
ΔYi=β1ΔXi+β3DiΔXiΔYiΔXi=β1+β3Di
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi+ΔYi=β0+β1(Xi+ΔXi)β2Di+β3((Xi+ΔXi)Di)
ΔYi=β1ΔXi+β3DiΔXiΔYiΔXi=β1+β3Di
The effect of X→Y depends on the value of Di!
β3: increment to the effect of X→Y when Di=1 (vs. Di=0)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
^β0: E[Yi] for Xi=0 and Di=0
β1: Marginal effect of Xi→Yi for Di=0
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
^β0: E[Yi] for Xi=0 and Di=0
β1: Marginal effect of Xi→Yi for Di=0
β2: Marginal effect on Yi of difference between Di=0 and Di=1
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
^β0: E[Yi] for Xi=0 and Di=0
β1: Marginal effect of Xi→Yi for Di=0
β2: Marginal effect on Yi of difference between Di=0 and Di=1
β3: The difference of the marginal effect of Xi→Yi between Di=0 and Di=1
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
^β0: E[Yi] for Xi=0 and Di=0
β1: Marginal effect of Xi→Yi for Di=0
β2: Marginal effect on Yi of difference between Di=0 and Di=1
β3: The difference of the marginal effect of Xi→Yi between Di=0 and Di=1
This is a bit awkward, easier to think about the two regression lines:
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
For Di=0 Group: ^Yi=^β0+^β1Xi
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
For Di=0 Group: ^Yi=^β0+^β1Xi
For Di=1 Group: ^Yi=(^β0+^β2)+(^β1+^β3)Xi
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
For Di=0 Group: ^Yi=^β0+^β1Xi
For Di=1 Group: ^Yi=(^β0+^β2)+(^β1+^β3)Xi
^β2: difference in intercept between groups
^β3: difference in slope between groups
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
For Di=0 Group: ^Yi=^β0+^β1Xi
For Di=1 Group: ^Yi=(^β0+^β2)+(^β1+^β3)Xi
^β2: difference in intercept between groups
^β3: difference in slope between groups
Example: ^wagei=^β0+^β1experi+^β2femalei+^β3(experi×femalei)
Example: ^wagei=^β0+^β1experi+^β2femalei+^β3(experi×femalei)
Example: ^wagei=^β0+^β1experi+^β2femalei+^β3(experi×femalei)
For males (female=0): ^wagei=^β0+^β1exper
For females (female=1): ^wagei=(^β0+^β2)⏟intercept+(^β1+^β3)⏟slopeexper
interaction_plot <- ggplot(data = wages)+ aes(x = exper, y = wage, color = as.factor(Gender))+ # make factor geom_point(alpha = 0.5)+ scale_y_continuous(labels=scales::dollar)+ labs(x = "Experience (Years)", y = "Wage")+ scale_color_manual(values = c("Female" = "#e64173", "Male" = "#0047AB") )+ # setting custom colors guides(color=F)+ # hide legend theme_slidesinteraction_plot
color aesthetic uses a factor variableas.factor() in ggplot code
interaction_plot+ geom_smooth(method="lm")

interaction_plot+ geom_smooth(method="lm")+ facet_wrap(~Gender)

R: var1 * var2var1 * var2 (multiply)# both are identical in Rinteraction_reg <- lm(wage ~ exper * female, data = wages)interaction_reg <- lm(wage ~ exper + female + exper * female, data = wages)
| ABCDEFGHIJ0123456789 |
term <chr> | estimate <dbl> | std.error <dbl> | statistic <dbl> | p.value <dbl> |
|---|---|---|---|---|
| (Intercept) | 6.15827549 | 0.34167408 | 18.023830 | 7.998534e-57 |
| exper | 0.05360476 | 0.01543716 | 3.472450 | 5.585255e-04 |
| female | -1.54654677 | 0.48186030 | -3.209534 | 1.411253e-03 |
| exper:female | -0.05506989 | 0.02217496 | -2.483427 | 1.332533e-02 |
library(huxtable)huxreg(interaction_reg, coefs = c("Constant" = "(Intercept)", "Experience" = "exper", "Female" = "female", "Experience * Female" = "exper:female"), statistics = c("N" = "nobs", "R-Squared" = "r.squared", "SER" = "sigma"), number_format = 2)
| (1) | |
|---|---|
| Constant | 6.16 *** |
| (0.34) | |
| Experience | 0.05 *** |
| (0.02) | |
| Female | -1.55 ** |
| (0.48) | |
| Experience * Female | -0.06 * |
| (0.02) | |
| N | 526 |
| R-Squared | 0.14 |
| SER | 3.44 |
| *** p < 0.001; ** p < 0.01; * p < 0.05. | |
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1:
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^β2:
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^β2: Women with 0 years of experience earn $1.55 less than men
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^β2: Women with 0 years of experience earn $1.55 less than men
^β3:
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^β2: Women with 0 years of experience earn $1.55 less than men
^β3: Women earn $0.06 less than men for every additional year of experience
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for men (female=0) ^wagei=6.16+0.05Experiencei
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for men (female=0) ^wagei=6.16+0.05Experiencei
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for men (female=0) ^wagei=6.16+0.05Experiencei
Men with 0 years of experience earn $6.16 on average
For every additional year of experience, men earn $0.05 more on average
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for women (female=1) ^wagei=6.16+0.05Experiencei−1.55(1)−0.06Experiencei×(1)=(6.16−1.55)+(0.05−0.06)Experiencei=4.61−0.01Experiencei
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for women (female=1) ^wagei=6.16+0.05Experiencei−1.55(1)−0.06Experiencei×(1)=(6.16−1.55)+(0.05−0.06)Experiencei=4.61−0.01Experiencei
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for women (female=1) ^wagei=6.16+0.05Experiencei−1.55(1)−0.06Experiencei×(1)=(6.16−1.55)+(0.05−0.06)Experiencei=4.61−0.01Experiencei
Women with 0 years of experience earn $4.61 on average
For every additional year of experience, women earn $0.01 less on average
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) 6.16 0.342 18.0 8.00e-57## 2 exper 0.0536 0.0154 3.47 5.59e- 4## 3 female -1.55 0.482 -3.21 1.41e- 3## 4 exper:female -0.0551 0.0222 -2.48 1.33e- 2Are slopes & intercepts of the 2 regressions statistically significantly different?
Are intercepts different? H0:β2=0
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) 6.16 0.342 18.0 8.00e-57## 2 exper 0.0536 0.0154 3.47 5.59e- 4## 3 female -1.55 0.482 -3.21 1.41e- 3## 4 exper:female -0.0551 0.0222 -2.48 1.33e- 2Are slopes & intercepts of the 2 regressions statistically significantly different?
Are intercepts different? H0:β2=0
Are slopes different? H0:β3=0
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) 6.16 0.342 18.0 8.00e-57## 2 exper 0.0536 0.0154 3.47 5.59e- 4## 3 female -1.55 0.482 -3.21 1.41e- 3## 4 exper:female -0.0551 0.0222 -2.48 1.33e- 2

Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
D1i and D2i are dummy variables
^β1: effect on Y of going from D1i=0 to D1i=1 when D2i=0
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
D1i and D2i are dummy variables
^β1: effect on Y of going from D1i=0 to D1i=1 when D2i=0
^β2: effect on Y of going from D2i=0 to D2i=1 when D1i=0
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
D1i and D2i are dummy variables
^β1: effect on Y of going from D1i=0 to D1i=1 when D2i=0
^β2: effect on Y of going from D2i=0 to D2i=1 when D1i=0
^β3: effect on Y of going from D1i=0 to D1i=1 when D2i=1
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
D1i and D2i are dummy variables
^β1: effect on Y of going from D1i=0 to D1i=1 when D2i=0
^β2: effect on Y of going from D2i=0 to D2i=1 when D1i=0
^β3: effect on Y of going from D1i=0 to D1i=1 when D2i=1
As always, best to think logically about possibilities (when each dummy =0 or =1)
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
E(Yi|D1i=0,D2i=d2)=β0+β2d2E(Yi|D1i=1,D2i=d2)=β0+β1(1)+β2d2+β3(1)d2
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
E(Yi|D1i=0,D2i=d2)=β0+β2d2E(Yi|D1i=1,D2i=d2)=β0+β1(1)+β2d2+β3(1)d2
β1+β3d2
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
1) Unmarried men (femalei=0,marriedi=0) ^wagei=^β0
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
1) Unmarried men (femalei=0,marriedi=0) ^wagei=^β0
2) Married men (femalei=0,marriedi=1) ^wagei=^β0+^β2
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
1) Unmarried men (femalei=0,marriedi=0) ^wagei=^β0
2) Married men (femalei=0,marriedi=1) ^wagei=^β0+^β2
3) Unmarried women (femalei=1,marriedi=0) ^wagei=^β0+^β1
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
1) Unmarried men (femalei=0,marriedi=0) ^wagei=^β0
2) Married men (femalei=0,marriedi=1) ^wagei=^β0+^β2
3) Unmarried women (femalei=1,marriedi=0) ^wagei=^β0+^β1
4) Married women (femalei=1,marriedi=1) ^wagei=^β0+^β1+^β2+^β3
# get average wage for unmarried menwages %>% filter(female == 0, married == 0) %>% summarize(mean = mean(wage))
## mean## 1 5.168023# get average wage for married menwages %>% filter(female == 0, married == 1) %>% summarize(mean = mean(wage))
## mean## 1 7.983032# get average wage for unmarried womenwages %>% filter(female == 1, married == 0) %>% summarize(mean = mean(wage))
## mean## 1 4.611583# get average wage for married womenwages %>% filter(female == 1, married == 1) %>% summarize(mean = mean(wage))
## mean## 1 4.565909^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
reg_dummies <- lm(wage ~ female + married + female:married, data = wages)reg_dummies %>% tidy()
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) 5.17 0.361 14.3 2.26e-39## 2 female -0.556 0.474 -1.18 2.41e- 1## 3 married 2.82 0.436 6.45 2.53e-10## 4 female:married -2.86 0.608 -4.71 3.20e- 6library(huxtable)huxreg(reg_dummies, coefs = c("Constant" = "(Intercept)", "Female" = "female", "Married" = "married", "Female * Married" = "female:married"), statistics = c("N" = "nobs", "R-Squared" = "r.squared", "SER" = "sigma"), number_format = 2)
| (1) | |
|---|---|
| Constant | 5.17 *** |
| (0.36) | |
| Female | -0.56 |
| (0.47) | |
| Married | 2.82 *** |
| (0.44) | |
| Female * Married | -2.86 *** |
| (0.61) | |
| N | 526 |
| R-Squared | 0.18 |
| SER | 3.35 |
| *** p < 0.001; ** p < 0.01; * p < 0.05. | |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |


Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi+ΔYi=β0+β1(X1+ΔX1i)β2X2i+β3((X1i+ΔX1i)×X2i)
Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi+ΔYi=β0+β1(X1+ΔX1i)β2X2i+β3((X1i+ΔX1i)×X2i)
Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi+ΔYi=β0+β1(X1+ΔX1i)β2X2i+β3((X1i+ΔX1i)×X2i)
ΔYi=β1ΔX1i+β3X2iΔX1iΔYiΔX1i=β1+β3X2i
Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi+ΔYi=β0+β1(X1+ΔX1i)β2X2i+β3((X1i+ΔX1i)×X2i)
ΔYi=β1ΔX1i+β3X2iΔX1iΔYiΔX1i=β1+β3X2i
Example: Do education and experience interact in their determination of wages?
^wagei=^β0+^β1educi+^β2experi+^β3(educi×experi)
ΔwageΔeduc=^β1+β3experi
ΔwageΔexper=^β2+β3educi
reg_cont <- lm(wage ~ educ + exper + educ:exper, data = wages)reg_cont %>% tidy()
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) -2.86 1.18 -2.42 1.58e- 2## 2 educ 0.602 0.0899 6.69 5.64e-11## 3 exper 0.0458 0.0426 1.07 2.83e- 1## 4 educ:exper 0.00206 0.00349 0.591 5.55e- 1library(huxtable)huxreg(reg_cont, coefs = c("Constant" = "(Intercept)", "Education" = "educ", "Experience" = "exper", "Education * Experience" = "educ:exper"), statistics = c("N" = "nobs", "R-Squared" = "r.squared", "SER" = "sigma"), number_format = 3)
| (1) | |
|---|---|
| Constant | -2.860 * |
| (1.181) | |
| Education | 0.602 *** |
| (0.090) | |
| Experience | 0.046 |
| (0.043) | |
| Education * Experience | 0.002 |
| (0.003) | |
| N | 526 |
| R-Squared | 0.226 |
| SER | 3.259 |
| *** p < 0.001; ** p < 0.01; * p < 0.05. | |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Education on Wages by Years of Experience:
| Experience | ΔwageΔeduc=^β1+^β3exper |
|---|---|
| 5 years | 0.602+0.002(5)=0.612 |
| 10 years | 0.602+0.002(10)=0.622 |
| 15 years | 0.602+0.002(15)=0.632 |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Education on Wages by Years of Experience:
| Experience | ΔwageΔeduc=^β1+^β3exper |
|---|---|
| 5 years | 0.602+0.002(5)=0.612 |
| 10 years | 0.602+0.002(10)=0.622 |
| 15 years | 0.602+0.002(15)=0.632 |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Experience on Wages by Years of Education:
| Education | ΔwageΔexper=^β2+^β3educ |
|---|---|
| 5 years | 0.047+0.002(5)=0.057 |
| 10 years | 0.047+0.002(10)=0.067 |
| 15 years | 0.047+0.002(15)=0.077 |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Experience on Wages by Years of Education:
| Education | ΔwageΔexper=^β2+^β3educ |
|---|---|
| 5 years | 0.047+0.002(5)=0.057 |
| 10 years | 0.047+0.002(10)=0.067 |
| 15 years | 0.047+0.002(15)=0.077 |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Experience on Wages by Years of Education:
| Education | ΔwageΔexper=^β2+^β3educ |
|---|---|
| 5 years | 0.047+0.002(5)=0.057 |
| 10 years | 0.047+0.002(10)=0.067 |
| 15 years | 0.047+0.002(15)=0.077 |
Marginal effect of experience → wages increases with more education
If you want to estimate the marginal effects more precisely, and graph them, see the appendix in today’s class page
Keyboard shortcuts
| ↑, ←, Pg Up, k | Go to previous slide |
| ↓, →, Pg Dn, Space, j | Go to next slide |
| Home | Go to first slide |
| End | Go to last slide |
| Number + Return | Go to specific slide |
| b / m / f | Toggle blackout / mirrored / fullscreen mode |
| c | Clone slideshow |
| p | Toggle presenter mode |
| t | Restart the presentation timer |
| ?, h | Toggle this help |
| o | Tile View: Overview of Slides |
| Esc | Back to slideshow |




Marginal effect of dummy variable: effect on Y of going from 0 to 1
Marginal effect of continuous variable: effect on Y of a 1 unit change in X
Example: Consider the gender pay gap again.
Example: Consider the gender pay gap again.
Example: Consider the gender pay gap again.
Example: Consider the gender pay gap again.
^wagesi=β0+β1Genderi+β2Experiencei
Depending on the types of variables, there are 3 possible types of interaction effects
We will look at each in turn
Depending on the types of variables, there are 3 possible types of interaction effects
We will look at each in turn
Depending on the types of variables, there are 3 possible types of interaction effects
We will look at each in turn
Depending on the types of variables, there are 3 possible types of interaction effects
We will look at each in turn


Yi=β0+β1Xi+β2Di+β3(Xi×Di) where Di={0,1}
Yi=β0+β1Xi+β2Di+β3(Xi×Di) where Di={0,1}
Yi=β0+β1Xi+β2Di+β3(Xi×Di) where Di={0,1}
β3 estimates the interaction effect between Xi and Di on Yi
What do the different coefficients (β)’s tell us?
Yi=β0+β1Xi+β2Di+β3Xi×Di
Yi=β0+β1Xi+β2Di+β3Xi×Di
^Yi=^β0+^β1Xi+^β2(0)+^β3Xi×(0)^Yi=^β0+^β1Xi
Yi=β0+β1Xi+β2Di+β3Xi×Di
^Yi=^β0+^β1Xi+^β2(0)+^β3Xi×(0)^Yi=^β0+^β1Xi
^Yi=^β0+^β1Xi+^β2(1)+^β3Xi×(1)^Yi=(^β0+^β2)+(^β1+^β3)Xi
Yi=β0+β1Xi+β2Di+β3Xi×Di
^Yi=^β0+^β1Xi+^β2(0)+^β3Xi×(0)^Yi=^β0+^β1Xi
^Yi=^β0+^β1Xi+^β2(1)+^β3Xi×(1)^Yi=(^β0+^β2)+(^β1+^β3)Xi

Yi=^β0+^β1Xi
Yi=(^β0+^β2)+(^β1+^β3)Xi
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi+ΔYi=β0+β1(Xi+ΔXi)β2Di+β3((Xi+ΔXi)Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi+ΔYi=β0+β1(Xi+ΔXi)β2Di+β3((Xi+ΔXi)Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi+ΔYi=β0+β1(Xi+ΔXi)β2Di+β3((Xi+ΔXi)Di)
ΔYi=β1ΔXi+β3DiΔXiΔYiΔXi=β1+β3Di
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi+ΔYi=β0+β1(Xi+ΔXi)β2Di+β3((Xi+ΔXi)Di)
ΔYi=β1ΔXi+β3DiΔXiΔYiΔXi=β1+β3Di
The effect of X→Y depends on the value of Di!
β3: increment to the effect of X→Y when Di=1 (vs. Di=0)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
^β0: E[Yi] for Xi=0 and Di=0
β1: Marginal effect of Xi→Yi for Di=0
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
^β0: E[Yi] for Xi=0 and Di=0
β1: Marginal effect of Xi→Yi for Di=0
β2: Marginal effect on Yi of difference between Di=0 and Di=1
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
^β0: E[Yi] for Xi=0 and Di=0
β1: Marginal effect of Xi→Yi for Di=0
β2: Marginal effect on Yi of difference between Di=0 and Di=1
β3: The difference of the marginal effect of Xi→Yi between Di=0 and Di=1
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
^β0: E[Yi] for Xi=0 and Di=0
β1: Marginal effect of Xi→Yi for Di=0
β2: Marginal effect on Yi of difference between Di=0 and Di=1
β3: The difference of the marginal effect of Xi→Yi between Di=0 and Di=1
This is a bit awkward, easier to think about the two regression lines:
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
For Di=0 Group: ^Yi=^β0+^β1Xi
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
For Di=0 Group: ^Yi=^β0+^β1Xi
For Di=1 Group: ^Yi=(^β0+^β2)+(^β1+^β3)Xi
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
For Di=0 Group: ^Yi=^β0+^β1Xi
For Di=1 Group: ^Yi=(^β0+^β2)+(^β1+^β3)Xi
^β2: difference in intercept between groups
^β3: difference in slope between groups
Yi=β0+β1Xi+β2Di+β3(Xi×Di)
For Di=0 Group: ^Yi=^β0+^β1Xi
For Di=1 Group: ^Yi=(^β0+^β2)+(^β1+^β3)Xi
^β2: difference in intercept between groups
^β3: difference in slope between groups
Example: ^wagei=^β0+^β1experi+^β2femalei+^β3(experi×femalei)
Example: ^wagei=^β0+^β1experi+^β2femalei+^β3(experi×femalei)
Example: ^wagei=^β0+^β1experi+^β2femalei+^β3(experi×femalei)
For males (female=0): ^wagei=^β0+^β1exper
For females (female=1): ^wagei=(^β0+^β2)⏟intercept+(^β1+^β3)⏟slopeexper
interaction_plot <- ggplot(data = wages)+ aes(x = exper, y = wage, color = as.factor(Gender))+ # make factor geom_point(alpha = 0.5)+ scale_y_continuous(labels=scales::dollar)+ labs(x = "Experience (Years)", y = "Wage")+ scale_color_manual(values = c("Female" = "#e64173", "Male" = "#0047AB") )+ # setting custom colors guides(color=F)+ # hide legend theme_slidesinteraction_plot
color aesthetic uses a factor variableas.factor() in ggplot code
interaction_plot+ geom_smooth(method="lm")

interaction_plot+ geom_smooth(method="lm")+ facet_wrap(~Gender)

R: var1 * var2var1 * var2 (multiply)# both are identical in Rinteraction_reg <- lm(wage ~ exper * female, data = wages)interaction_reg <- lm(wage ~ exper + female + exper * female, data = wages)
| ABCDEFGHIJ0123456789 |
term <chr> | estimate <dbl> | std.error <dbl> | statistic <dbl> | p.value <dbl> |
|---|---|---|---|---|
| (Intercept) | 6.15827549 | 0.34167408 | 18.023830 | 7.998534e-57 |
| exper | 0.05360476 | 0.01543716 | 3.472450 | 5.585255e-04 |
| female | -1.54654677 | 0.48186030 | -3.209534 | 1.411253e-03 |
| exper:female | -0.05506989 | 0.02217496 | -2.483427 | 1.332533e-02 |
library(huxtable)huxreg(interaction_reg, coefs = c("Constant" = "(Intercept)", "Experience" = "exper", "Female" = "female", "Experience * Female" = "exper:female"), statistics = c("N" = "nobs", "R-Squared" = "r.squared", "SER" = "sigma"), number_format = 2)
| (1) | |
|---|---|
| Constant | 6.16 *** |
| (0.34) | |
| Experience | 0.05 *** |
| (0.02) | |
| Female | -1.55 ** |
| (0.48) | |
| Experience * Female | -0.06 * |
| (0.02) | |
| N | 526 |
| R-Squared | 0.14 |
| SER | 3.44 |
| *** p < 0.001; ** p < 0.01; * p < 0.05. | |
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1:
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^β2:
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^β2: Women with 0 years of experience earn $1.55 less than men
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^β2: Women with 0 years of experience earn $1.55 less than men
^β3:
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^β0: Men with 0 years of experience earn 6.16
^β1: For every additional year of experience, men earn $0.05
^β2: Women with 0 years of experience earn $1.55 less than men
^β3: Women earn $0.06 less than men for every additional year of experience
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for men (female=0) ^wagei=6.16+0.05Experiencei
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for men (female=0) ^wagei=6.16+0.05Experiencei
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for men (female=0) ^wagei=6.16+0.05Experiencei
Men with 0 years of experience earn $6.16 on average
For every additional year of experience, men earn $0.05 more on average
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for women (female=1) ^wagei=6.16+0.05Experiencei−1.55(1)−0.06Experiencei×(1)=(6.16−1.55)+(0.05−0.06)Experiencei=4.61−0.01Experiencei
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for women (female=1) ^wagei=6.16+0.05Experiencei−1.55(1)−0.06Experiencei×(1)=(6.16−1.55)+(0.05−0.06)Experiencei=4.61−0.01Experiencei
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
Regression for women (female=1) ^wagei=6.16+0.05Experiencei−1.55(1)−0.06Experiencei×(1)=(6.16−1.55)+(0.05−0.06)Experiencei=4.61−0.01Experiencei
Women with 0 years of experience earn $4.61 on average
For every additional year of experience, women earn $0.01 less on average
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) 6.16 0.342 18.0 8.00e-57## 2 exper 0.0536 0.0154 3.47 5.59e- 4## 3 female -1.55 0.482 -3.21 1.41e- 3## 4 exper:female -0.0551 0.0222 -2.48 1.33e- 2Are slopes & intercepts of the 2 regressions statistically significantly different?
Are intercepts different? H0:β2=0
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) 6.16 0.342 18.0 8.00e-57## 2 exper 0.0536 0.0154 3.47 5.59e- 4## 3 female -1.55 0.482 -3.21 1.41e- 3## 4 exper:female -0.0551 0.0222 -2.48 1.33e- 2Are slopes & intercepts of the 2 regressions statistically significantly different?
Are intercepts different? H0:β2=0
Are slopes different? H0:β3=0
^wagei=6.16+0.05Experiencei−1.55Femalei−0.06(Experiencei×Femalei)
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) 6.16 0.342 18.0 8.00e-57## 2 exper 0.0536 0.0154 3.47 5.59e- 4## 3 female -1.55 0.482 -3.21 1.41e- 3## 4 exper:female -0.0551 0.0222 -2.48 1.33e- 2

Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
D1i and D2i are dummy variables
^β1: effect on Y of going from D1i=0 to D1i=1 when D2i=0
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
D1i and D2i are dummy variables
^β1: effect on Y of going from D1i=0 to D1i=1 when D2i=0
^β2: effect on Y of going from D2i=0 to D2i=1 when D1i=0
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
D1i and D2i are dummy variables
^β1: effect on Y of going from D1i=0 to D1i=1 when D2i=0
^β2: effect on Y of going from D2i=0 to D2i=1 when D1i=0
^β3: effect on Y of going from D1i=0 to D1i=1 when D2i=1
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
D1i and D2i are dummy variables
^β1: effect on Y of going from D1i=0 to D1i=1 when D2i=0
^β2: effect on Y of going from D2i=0 to D2i=1 when D1i=0
^β3: effect on Y of going from D1i=0 to D1i=1 when D2i=1
As always, best to think logically about possibilities (when each dummy =0 or =1)
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
E(Yi|D1i=0,D2i=d2)=β0+β2d2E(Yi|D1i=1,D2i=d2)=β0+β1(1)+β2d2+β3(1)d2
Yi=β0+β1D1i+β2D2i+β3(D1i×D2i)
E(Yi|D1i=0,D2i=d2)=β0+β2d2E(Yi|D1i=1,D2i=d2)=β0+β1(1)+β2d2+β3(1)d2
β1+β3d2
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
1) Unmarried men (femalei=0,marriedi=0) ^wagei=^β0
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
1) Unmarried men (femalei=0,marriedi=0) ^wagei=^β0
2) Married men (femalei=0,marriedi=1) ^wagei=^β0+^β2
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
1) Unmarried men (femalei=0,marriedi=0) ^wagei=^β0
2) Married men (femalei=0,marriedi=1) ^wagei=^β0+^β2
3) Unmarried women (femalei=1,marriedi=0) ^wagei=^β0+^β1
Example: Does the gender pay gap change if a person is married vs. single?
^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
1) Unmarried men (femalei=0,marriedi=0) ^wagei=^β0
2) Married men (femalei=0,marriedi=1) ^wagei=^β0+^β2
3) Unmarried women (femalei=1,marriedi=0) ^wagei=^β0+^β1
4) Married women (femalei=1,marriedi=1) ^wagei=^β0+^β1+^β2+^β3
# get average wage for unmarried menwages %>% filter(female == 0, married == 0) %>% summarize(mean = mean(wage))
## mean## 1 5.168023# get average wage for married menwages %>% filter(female == 0, married == 1) %>% summarize(mean = mean(wage))
## mean## 1 7.983032# get average wage for unmarried womenwages %>% filter(female == 1, married == 0) %>% summarize(mean = mean(wage))
## mean## 1 4.611583# get average wage for married womenwages %>% filter(female == 1, married == 1) %>% summarize(mean = mean(wage))
## mean## 1 4.565909^wagei=^β0+^β1femalei+^β2marriedi+^β3(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
reg_dummies <- lm(wage ~ female + married + female:married, data = wages)reg_dummies %>% tidy()
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) 5.17 0.361 14.3 2.26e-39## 2 female -0.556 0.474 -1.18 2.41e- 1## 3 married 2.82 0.436 6.45 2.53e-10## 4 female:married -2.86 0.608 -4.71 3.20e- 6library(huxtable)huxreg(reg_dummies, coefs = c("Constant" = "(Intercept)", "Female" = "female", "Married" = "married", "Female * Married" = "female:married"), statistics = c("N" = "nobs", "R-Squared" = "r.squared", "SER" = "sigma"), number_format = 2)
| (1) | |
|---|---|
| Constant | 5.17 *** |
| (0.36) | |
| Female | -0.56 |
| (0.47) | |
| Married | 2.82 *** |
| (0.44) | |
| Female * Married | -2.86 *** |
| (0.61) | |
| N | 526 |
| R-Squared | 0.18 |
| SER | 3.35 |
| *** p < 0.001; ** p < 0.01; * p < 0.05. | |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |
^wagei=5.17−0.56femalei+2.82marriedi−2.86(femalei×marriedi)
| Men | Women | |
|---|---|---|
| Unmarried | $5.17 | $4.61 |
| Married | $7.98 | $4.57 |


Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi+ΔYi=β0+β1(X1+ΔX1i)β2X2i+β3((X1i+ΔX1i)×X2i)
Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi+ΔYi=β0+β1(X1+ΔX1i)β2X2i+β3((X1i+ΔX1i)×X2i)
Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi+ΔYi=β0+β1(X1+ΔX1i)β2X2i+β3((X1i+ΔX1i)×X2i)
ΔYi=β1ΔX1i+β3X2iΔX1iΔYiΔX1i=β1+β3X2i
Yi=β0+β1X1i+β2X2i+β3(X1i×X2i)
Yi+ΔYi=β0+β1(X1+ΔX1i)β2X2i+β3((X1i+ΔX1i)×X2i)
ΔYi=β1ΔX1i+β3X2iΔX1iΔYiΔX1i=β1+β3X2i
Example: Do education and experience interact in their determination of wages?
^wagei=^β0+^β1educi+^β2experi+^β3(educi×experi)
ΔwageΔeduc=^β1+β3experi
ΔwageΔexper=^β2+β3educi
reg_cont <- lm(wage ~ educ + exper + educ:exper, data = wages)reg_cont %>% tidy()
## # A tibble: 4 × 5## term estimate std.error statistic p.value## <chr> <dbl> <dbl> <dbl> <dbl>## 1 (Intercept) -2.86 1.18 -2.42 1.58e- 2## 2 educ 0.602 0.0899 6.69 5.64e-11## 3 exper 0.0458 0.0426 1.07 2.83e- 1## 4 educ:exper 0.00206 0.00349 0.591 5.55e- 1library(huxtable)huxreg(reg_cont, coefs = c("Constant" = "(Intercept)", "Education" = "educ", "Experience" = "exper", "Education * Experience" = "educ:exper"), statistics = c("N" = "nobs", "R-Squared" = "r.squared", "SER" = "sigma"), number_format = 3)
| (1) | |
|---|---|
| Constant | -2.860 * |
| (1.181) | |
| Education | 0.602 *** |
| (0.090) | |
| Experience | 0.046 |
| (0.043) | |
| Education * Experience | 0.002 |
| (0.003) | |
| N | 526 |
| R-Squared | 0.226 |
| SER | 3.259 |
| *** p < 0.001; ** p < 0.01; * p < 0.05. | |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Education on Wages by Years of Experience:
| Experience | ΔwageΔeduc=^β1+^β3exper |
|---|---|
| 5 years | 0.602+0.002(5)=0.612 |
| 10 years | 0.602+0.002(10)=0.622 |
| 15 years | 0.602+0.002(15)=0.632 |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Education on Wages by Years of Experience:
| Experience | ΔwageΔeduc=^β1+^β3exper |
|---|---|
| 5 years | 0.602+0.002(5)=0.612 |
| 10 years | 0.602+0.002(10)=0.622 |
| 15 years | 0.602+0.002(15)=0.632 |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Experience on Wages by Years of Education:
| Education | ΔwageΔexper=^β2+^β3educ |
|---|---|
| 5 years | 0.047+0.002(5)=0.057 |
| 10 years | 0.047+0.002(10)=0.067 |
| 15 years | 0.047+0.002(15)=0.077 |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Experience on Wages by Years of Education:
| Education | ΔwageΔexper=^β2+^β3educ |
|---|---|
| 5 years | 0.047+0.002(5)=0.057 |
| 10 years | 0.047+0.002(10)=0.067 |
| 15 years | 0.047+0.002(15)=0.077 |
^wagesi=−2.860+0.602educi+0.047experi+0.002(educi×experi)
Marginal Effect of Experience on Wages by Years of Education:
| Education | ΔwageΔexper=^β2+^β3educ |
|---|---|
| 5 years | 0.047+0.002(5)=0.057 |
| 10 years | 0.047+0.002(10)=0.067 |
| 15 years | 0.047+0.002(15)=0.077 |
Marginal effect of experience → wages increases with more education
If you want to estimate the marginal effects more precisely, and graph them, see the appendix in today’s class page