cps78_85
library(wooldridge)
## Warning: package 'wooldridge' was built under R version 4.1.3
library(rmarkdown)
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
data("cps78_85")
head(cps78_85)
## educ south nonwhite female married exper expersq union lwage age year y85
## 1 12 0 0 0 0 8 64 0 1.2150 25 78 0
## 2 12 0 0 1 1 30 900 1 1.6094 47 78 0
## 3 6 0 0 0 1 38 1444 1 2.1401 49 78 0
## 4 12 0 0 0 1 19 361 1 2.0732 36 78 0
## 5 12 0 0 0 1 11 121 0 1.6490 28 78 0
## 6 8 0 0 0 1 43 1849 0 1.7148 56 78 0
## y85fem y85educ y85union
## 1 0 0 0
## 2 0 0 0
## 3 0 0 0
## 4 0 0 0
## 5 0 0 0
## 6 0 0 0
tail(cps78_85)
## educ south nonwhite female married exper expersq union lwage age year y85
## 1079 17 0 1 1 1 25 625 1 3.1463 48 85 1
## 1080 12 0 1 1 1 18 324 0 2.1552 36 85 1
## 1081 14 0 0 1 0 13 169 0 2.9178 33 85 1
## 1082 18 0 0 0 1 8 64 0 3.1001 32 85 1
## 1083 12 0 1 0 0 14 196 1 2.7887 32 85 1
## 1084 12 0 0 0 0 9 81 0 2.9689 27 85 1
## y85fem y85educ y85union
## 1079 1 17 1
## 1080 1 12 0
## 1081 1 14 0
## 1082 0 18 0
## 1083 0 12 1
## 1084 0 12 0
summary(cps78_85)
## educ south nonwhite female
## Min. : 1.00 Min. :0.0000 Min. :0.0000 Min. :0.000
## 1st Qu.:12.00 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.:0.000
## Median :12.00 Median :0.0000 Median :0.0000 Median :0.000
## Mean :12.77 Mean :0.2943 Mean :0.1144 Mean :0.417
## 3rd Qu.:14.00 3rd Qu.:1.0000 3rd Qu.:0.0000 3rd Qu.:1.000
## Max. :18.00 Max. :1.0000 Max. :1.0000 Max. :1.000
## married exper expersq union
## Min. :0.0000 Min. : 0.00 Min. : 0.0 Min. :0.0000
## 1st Qu.:0.0000 1st Qu.: 8.00 1st Qu.: 64.0 1st Qu.:0.0000
## Median :1.0000 Median :15.00 Median : 225.0 Median :0.0000
## Mean :0.6541 Mean :18.28 Mean : 499.8 Mean :0.2435
## 3rd Qu.:1.0000 3rd Qu.:28.00 3rd Qu.: 784.0 3rd Qu.:0.0000
## Max. :1.0000 Max. :55.00 Max. :3025.0 Max. :1.0000
## lwage age year y85
## Min. :-0.470 Min. :18.00 Min. :78.00 Min. :0.0000
## 1st Qu.: 1.470 1st Qu.:27.00 1st Qu.:78.00 1st Qu.:0.0000
## Median : 1.833 Median :34.00 Median :78.00 Median :0.0000
## Mean : 1.867 Mean :36.54 Mean :81.45 Mean :0.4926
## 3rd Qu.: 2.225 3rd Qu.:46.00 3rd Qu.:85.00 3rd Qu.:1.0000
## Max. : 3.796 Max. :64.00 Max. :85.00 Max. :1.0000
## y85fem y85educ y85union
## Min. :0.000 Min. : 0.000 Min. :0.00000
## 1st Qu.:0.000 1st Qu.: 0.000 1st Qu.:0.00000
## Median :0.000 Median : 0.000 Median :0.00000
## Mean :0.226 Mean : 6.413 Mean :0.08856
## 3rd Qu.:0.000 3rd Qu.:12.000 3rd Qu.:0.00000
## Max. :1.000 Max. :18.000 Max. :1.00000
paged_table(cps78_85)
yaya <-cps78_85 %>%
group_by(year) %>%
summarise(yaya = n ())
summary(yaya)
## year yaya
## Min. :78.00 Min. :534
## 1st Qu.:79.75 1st Qu.:538
## Median :81.50 Median :542
## Mean :81.50 Mean :542
## 3rd Qu.:83.25 3rd Qu.:546
## Max. :85.00 Max. :550
paged_table(yaya)
summary(lm(lwage ~ y85+(educ + female) + exper + I(exper^2) + union, data = cps78_85))
##
## Call:
## lm(formula = lwage ~ y85 + (educ + female) + exper + I(exper^2) +
## union, data = cps78_85)
##
## Residuals:
## Min 1Q Median 3Q Max
## -2.5438 -0.2606 0.0085 0.2693 2.1522
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 3.322e-01 7.549e-02 4.400 1.19e-05 ***
## y85 3.900e-01 2.556e-02 15.260 < 2e-16 ***
## educ 8.333e-02 5.069e-03 16.438 < 2e-16 ***
## female -2.747e-01 2.583e-02 -10.638 < 2e-16 ***
## exper 2.968e-02 3.573e-03 8.306 2.94e-16 ***
## I(exper^2) -3.982e-04 7.771e-05 -5.124 3.54e-07 ***
## union 2.043e-01 3.032e-02 6.739 2.59e-11 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4136 on 1077 degrees of freedom
## Multiple R-squared: 0.4225, Adjusted R-squared: 0.4193
## F-statistic: 131.3 on 6 and 1077 DF, p-value: < 2.2e-16
summary(lm(lwage ~ y85+ y85educ+y85fem+(educ + female) + exper + I(exper^2) + union -1, data = cps78_85))
##
## Call:
## lm(formula = lwage ~ y85 + y85educ + y85fem + (educ + female) +
## exper + I(exper^2) + union - 1, data = cps78_85)
##
## Residuals:
## Min 1Q Median 3Q Max
## -2.56331 -0.24500 0.01702 0.27741 2.16580
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## y85 4.794e-01 1.006e-01 4.767 2.13e-06 ***
## y85educ -8.124e-03 7.710e-03 -1.054 0.292
## y85fem 5.859e-02 5.157e-02 1.136 0.256
## educ 1.046e-01 2.767e-03 37.816 < 2e-16 ***
## female -2.924e-01 3.667e-02 -7.973 3.92e-15 ***
## exper 3.355e-02 3.512e-03 9.552 < 2e-16 ***
## I(exper^2) -4.391e-04 7.794e-05 -5.634 2.24e-08 ***
## union 2.199e-01 3.040e-02 7.232 9.00e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4171 on 1076 degrees of freedom
## Multiple R-squared: 0.9543, Adjusted R-squared: 0.954
## F-statistic: 2810 on 8 and 1076 DF, p-value: < 2.2e-16
benimdatam <-cps78_85 %>%
mutate(benimsabit = 1) %>%
mutate(benimsabit = ifelse(married == 1,0,4) )
summary(benimdatam)
## educ south nonwhite female
## Min. : 1.00 Min. :0.0000 Min. :0.0000 Min. :0.000
## 1st Qu.:12.00 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.:0.000
## Median :12.00 Median :0.0000 Median :0.0000 Median :0.000
## Mean :12.77 Mean :0.2943 Mean :0.1144 Mean :0.417
## 3rd Qu.:14.00 3rd Qu.:1.0000 3rd Qu.:0.0000 3rd Qu.:1.000
## Max. :18.00 Max. :1.0000 Max. :1.0000 Max. :1.000
## married exper expersq union
## Min. :0.0000 Min. : 0.00 Min. : 0.0 Min. :0.0000
## 1st Qu.:0.0000 1st Qu.: 8.00 1st Qu.: 64.0 1st Qu.:0.0000
## Median :1.0000 Median :15.00 Median : 225.0 Median :0.0000
## Mean :0.6541 Mean :18.28 Mean : 499.8 Mean :0.2435
## 3rd Qu.:1.0000 3rd Qu.:28.00 3rd Qu.: 784.0 3rd Qu.:0.0000
## Max. :1.0000 Max. :55.00 Max. :3025.0 Max. :1.0000
## lwage age year y85
## Min. :-0.470 Min. :18.00 Min. :78.00 Min. :0.0000
## 1st Qu.: 1.470 1st Qu.:27.00 1st Qu.:78.00 1st Qu.:0.0000
## Median : 1.833 Median :34.00 Median :78.00 Median :0.0000
## Mean : 1.867 Mean :36.54 Mean :81.45 Mean :0.4926
## 3rd Qu.: 2.225 3rd Qu.:46.00 3rd Qu.:85.00 3rd Qu.:1.0000
## Max. : 3.796 Max. :64.00 Max. :85.00 Max. :1.0000
## y85fem y85educ y85union benimsabit
## Min. :0.000 Min. : 0.000 Min. :0.00000 Min. :0.000
## 1st Qu.:0.000 1st Qu.: 0.000 1st Qu.:0.00000 1st Qu.:0.000
## Median :0.000 Median : 0.000 Median :0.00000 Median :0.000
## Mean :0.226 Mean : 6.413 Mean :0.08856 Mean :1.384
## 3rd Qu.:0.000 3rd Qu.:12.000 3rd Qu.:0.00000 3rd Qu.:4.000
## Max. :1.000 Max. :18.000 Max. :1.00000 Max. :4.000
summary(lm(benimdatam))
##
## Call:
## lm(formula = benimdatam)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.0778 -0.0143 -0.0033 0.0077 3.9201
##
## Coefficients: (2 not defined because of singularities)
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 7.178e+00 4.216e-01 17.026 <2e-16 ***
## south -4.977e-03 8.312e-03 -0.599 0.5495
## nonwhite -6.889e-04 1.187e-02 -0.058 0.9537
## female 3.980e-03 1.126e-02 0.353 0.7239
## married -6.250e-03 8.408e-03 -0.743 0.4574
## exper -9.972e-01 2.353e-03 -423.803 <2e-16 ***
## expersq 4.442e-05 2.360e-05 1.882 0.0601 .
## union -3.296e-03 1.173e-02 -0.281 0.7787
## lwage 1.636e-02 9.002e-03 1.817 0.0695 .
## age 9.946e-01 2.071e-03 480.192 <2e-16 ***
## year -1.549e-01 5.115e-03 -30.281 <2e-16 ***
## y85 NA NA NA NA
## y85fem 1.663e-02 1.532e-02 1.086 0.2778
## y85educ 6.584e-03 2.751e-03 2.394 0.0169 *
## y85union -2.665e-03 1.791e-02 -0.149 0.8818
## benimsabit NA NA NA NA
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.1211 on 1070 degrees of freedom
## Multiple R-squared: 0.998, Adjusted R-squared: 0.998
## F-statistic: 4.153e+04 on 13 and 1070 DF, p-value: < 2.2e-16