Как построить график взаимодействия с моделью PLM? - PullRequest
0 голосов
/ 27 мая 2020

Я столкнулся с трудностями при планировании своих взаимодействий. Я использую модель plm с данными моей панели, но похоже, что большое количество пакетов построения графиков регрессии могут работать с lm или glm, но не с plm? Это выглядит как действительно хороший ресурс (использует ggplot для построения графиков взаимодействий), но когда я попытался использовать interact.plot , я получил эту ошибку: «Ошибка: ~ не похоже на односторонняя или двусторонняя формула ". Есть идеи, как это решить?

Вот мой режим:

FE_hours_interact <- plm(her_share_hours8 ~ her_share_hours4 + firstchild2 + gap4 + recession1 + recession2 + recession3 + firstchild2 * gap4 + firstchild2 * recession1 + firstchild2 * recession2 + firstchild2 * recession3 + firstchild2 * year_recode, data = Inflation_merge, index=c("statefip"), model="within")

Вот первые 5 строк моего df:

structure(list(X1 = c(1, 2, 3, 4, 5), X1_3 = c(1, 2, 3, 4, 5), 
X1_2 = c(1, 2, 3, 4, 5), X1_1 = c(1, 2, 3, 4, 5), cpsid = c(19821005481000, 
19820105935200, 19820904891200, 19820602523200, 19820704381600
), cpsidpf = c(19821005481002, 19820105935202, 19820904891202, 
19820602523202, 19820704381602), cpsidpm = c(19821005481001, 
19820105935201, 19820904891201, 19820602523201, 19820704381601
), totmonths = c(8, 8, 8, 8, 8), firstyear = c(1982, 1982, 
1982, 1982, 1982), firstmonth = c(10, 1, 9, 6, 7), lastyear = c(1984, 
1983, 1983, 1983, 1983), intv1976 = c(0, 0, 0, 0, 0), lastmonth = c(1, 
4, 12, 9, 10), numcplrecHH = c(1, 1, 1, 1, 1), statecensus = c(63, 
63, 63, 63, 63), statefip = c(1, 1, 1, 1, 1), match458 = c(1, 
1, 1, 1, 1), birthsample = c(1, 0, 0, 0, 0), firstchild = c(0, 
0, 0, 0, 0), secbsamp = c(0, 0, 0, 0, 0), secchild = c(0, 
0, 0, 0, 0), married = c(1, 1, 1, 1, 1), cohab = c(0, 0, 
0, 0, 0), cohabtomarr = c(0, 0, 0, 0, 0), f45 = c(1, 1, 1, 
1, 1), insampmo = c(1, 1, 1, 1, 1), inmis458 = c(1, 1, 1, 
1, 1), inmis1 = c(1, 1, 1, 1, 1), inmis2 = c(1, 1, 1, 1, 
1), inmis3 = c(1, 1, 1, 1, 1), inmis4 = c(1, 1, 1, 1, 1), 
inmis5 = c(1, 1, 1, 1, 1), inmis6 = c(1, 1, 1, 1, 1), inmis7 = c(1, 
1, 1, 1, 1), inmis8 = c(1, 1, 1, 1, 1), aff1 = c(0, 0, 0, 
0, 0), aff2 = c(0, 0, 0, 0, 0), aff3 = c(0, 0, 0, 0, 0), 
aff4 = c(0, 0, 0, 0, 0), aff5 = c(0, 0, 0, 0, 0), aff6 = c(0, 
0, 0, 0, 0), aff7 = c(0, 0, 0, 0, 0), aff8 = c(0, 0, 0, 0, 
0), afm1 = c(0, 0, 0, 0, 0), afm2 = c(0, 0, 0, 0, 0), afm3 = c(0, 
0, 0, 0, 0), afm4 = c(0, 0, 0, 0, 0), afm5 = c(0, 0, 0, 0, 
0), afm6 = c(0, 0, 0, 0, 0), afm7 = c(0, 0, 0, 0, 0), afm8 = c(0, 
0, 0, 0, 0), agef1 = c(31, 30, 38, 35, 27), agef2 = c(31, 
30, 38, 35, 27), agef3 = c(31, 30, 39, 35, 27), agef4 = c(31, 
30, 39, 35, 27), agef5 = c(32, 31, 39, 36, 28), agef6 = c(32, 
31, 39, 36, 28), agef7 = c(32, 31, 39, 36, 28), agef8 = c(32, 
31, 39, 36, 28), agem1 = c(33, 30, 46, 37, 28), agem2 = c(33, 
30, 46, 37, 28), agem3 = c(33, 30, 46, 37, 28), agem4 = c(33, 
30, 46, 38, 28), agem5 = c(34, 31, 47, 38, 29), agem6 = c(34, 
31, 47, 38, 29), agem7 = c(34, 31, 47, 38, 29), agem8 = c(34, 
31, 47, 39, 29), agematch_mis12 = c(1, 1, 1, 1, 1), agematch_mis14 = c(1, 
1, 1, 1, 1), agematch_mis15 = c(1, 1, 1, 1, 1), agematch_mis23 = c(1, 
1, 1, 1, 1), agematch_mis26 = c(1, 1, 1, 1, 1), agematch_mis34 = c(1, 
1, 1, 1, 1), agematch_mis37 = c(1, 1, 1, 1, 1), agematch_mis45 = c(1, 
1, 1, 1, 1), agematch_mis48 = c(1, 1, 1, 1, 1), agematch_mis56 = c(1, 
1, 1, 1, 1), agematch_mis58 = c(1, 1, 1, 1, 1), agematch_mis67 = c(1, 
1, 1, 1, 1), agematch_mis78 = c(1, 1, 1, 1, 1), agematchf_mis12 = c(1, 
1, 1, 1, 1), agematchf_mis14 = c(1, 1, 1, 1, 1), agematchf_mis15 = c(1, 
1, 1, 1, 1), agematchf_mis23 = c(1, 1, 1, 1, 1), agematchf_mis26 = c(1, 
1, 1, 1, 1), agematchf_mis34 = c(1, 1, 1, 1, 1), agematchf_mis37 = c(1, 
1, 1, 1, 1), agematchf_mis45 = c(1, 1, 1, 1, 1), agematchf_mis48 = c(1, 
1, 1, 1, 1), agematchf_mis56 = c(1, 1, 1, 1, 1), agematchf_mis58 = c(1, 
1, 1, 1, 1), agematchf_mis67 = c(1, 1, 1, 1, 1), agematchf_mis78 = c(1, 
1, 1, 1, 1), agematchm_mis12 = c(1, 1, 1, 1, 1), agematchm_mis14 = c(1, 
1, 1, 1, 1), agematchm_mis15 = c(1, 1, 1, 1, 1), agematchm_mis23 = c(1, 
1, 1, 1, 1), agematchm_mis26 = c(1, 1, 1, 1, 1), agematchm_mis34 = c(1, 
1, 1, 1, 1), agematchm_mis37 = c(1, 1, 1, 1, 1), agematchm_mis45 = c(1, 
1, 1, 1, 1), agematchm_mis48 = c(1, 1, 1, 1, 1), agematchm_mis56 = c(1, 
1, 1, 1, 1), agematchm_mis58 = c(1, 1, 1, 1, 1), agematchm_mis67 = c(1, 
1, 1, 1, 1), agematchm_mis78 = c(1, 1, 1, 1, 1), classwkrf1 = c(21, 
21, 0, 0, 0), classwkrf2 = c(21, 21, 0, 0, 0), classwkrf3 = c(21, 
21, 0, 0, 0), classwkrf4 = c(21, 21, 0, 24, 21), classwkrf5 = c(21, 
21, 0, 0, 0), classwkrf6 = c(21, 21, 0, 0, 0), classwkrf7 = c(21, 
21, 0, 0, 0), classwkrf8 = c(21, 21, 0, 25, 0), classwkrm1 = c(21, 
21, 21, 24, 21), classwkrm2 = c(21, 21, 21, 24, 21), classwkrm3 = c(21, 
21, 21, 24, 21), classwkrm4 = c(21, 21, 21, 24, 21), classwkrm5 = c(21, 
21, 21, 25, 21), classwkrm6 = c(21, 21, 21, 25, 21), classwkrm7 = c(21, 
21, 21, 25, 21), classwkrm8 = c(21, 21, 21, 25, 21), cohab1 = c(0, 
0, 0, 0, 0), cohab2 = c(0, 0, 0, 0, 0), cohab3 = c(0, 0, 
0, 0, 0), cohab4 = c(0, 0, 0, 0, 0), cohab5 = c(0, 0, 0, 
0, 0), cohab6 = c(0, 0, 0, 0, 0), cohab7 = c(0, 0, 0, 0, 
0), cohab8 = c(0, 0, 0, 0, 0), earnstatf4 = c(2, 2, 0, 0, 
0), earnstatf8 = c(2, 2, 0, 0, 0), earnstatm4 = c(2, 2, 2, 
2, 2), earnstatm8 = c(2, 2, 2, 2, 2), earnweek82onf4 = c(670, 
225, NA, NA, NA), earnweek82onf8 = c(275, 240, NA, NA, NA
), earnweek82onm4 = c(960, 550, 432, 450, 150), earnweek82onm8 = c(460, 
320, 680, 650, 150), earnweekf1 = c(9999.99, 9999.99, 9999.99, 
9999.99, 9999.99), earnweekf2 = c(9999.99, 9999.99, 9999.99, 
9999.99, 9999.99), earnweekf3 = c(9999.99, 9999.99, 9999.99, 
9999.99, 9999.99), earnweekf4 = c(670, 225, 9999.99, 9999.99, 
9999.99), earnweekf5 = c(9999.99, 9999.99, 9999.99, 9999.99, 
9999.99), earnweekf6 = c(9999.99, 9999.99, 9999.99, 9999.99, 
9999.99), earnweekf7 = c(9999.99, 9999.99, 9999.99, 9999.99, 
9999.99), earnweekf8 = c(275, 240, 9999.99, 9999.99, 9999.99
), earnweekm1 = c(9999.99, 9999.99, 9999.99, 9999.99, 9999.99
), earnweekm2 = c(9999.99, 9999.99, 9999.99, 9999.99, 9999.99
), earnweekm3 = c(9999.99, 9999.99, 9999.99, 9999.99, 9999.99
), earnweekm4 = c(960, 550, 432, 450, 150), earnweekm5 = c(9999.99, 
9999.99, 9999.99, 9999.99, 9999.99), earnweekm6 = c(9999.99, 
9999.99, 9999.99, 9999.99, 9999.99), earnweekm7 = c(9999.99, 
9999.99, 9999.99, 9999.99, 9999.99), earnweekm8 = c(460, 
320, 680, 650, 150), earnwtf4 = c(4752.14, 6119.42, 5005.91, 
6921.73, 6647.82), earnwtf8 = c(4530.97, 6340.22, 5291.25, 
6762.48, 6624.79), earnwtm4 = c(5179.23, 6346.1, 4577.67, 
6594.28, 6547.49), earnwtm8 = c(5086.8, 6607.46, 4889.55, 
6530.65, 6894.16), educf1 = c(90, 72, 73, 72, 72), educf2 = c(90, 
72, 73, 72, 72), educf3 = c(90, 72, 73, 72, 72), educf4 = c(90, 
72, 73, 72, 72), educf5 = c(90, 72, 80, 72, 72), educf6 = c(90, 
72, 80, 72, 72), educf7 = c(90, 72, 80, 72, 72), educf8 = c(90, 
72, 80, 72, 72), educm1 = c(122, 72, 80, 72, 72), educm2 = c(122, 
72, 80, 72, 72), educm3 = c(122, 72, 80, 72, 72), educm4 = c(122, 
72, 80, 72, 72), educm5 = c(122, 72, 80, 72, 72), educm6 = c(122, 
72, 80, 72, 72), educm7 = c(122, 72, 80, 72, 72), educm8 = c(122, 
72, 80, 72, 72), eldch1 = c(99, 8, 13, 13, 7), eldch2 = c(99, 
8, 13, 13, 7), eldch3 = c(99, 8, 13, 13, 7), eldch4 = c(99, 
8, 13, 11, 7), eldch5 = c(99, 9, 14, 14, 8), eldch6 = c(99, 
9, 14, 14, 8), eldch7 = c(99, 9, 14, 14, 8), eldch8 = c(99, 
9, 14, 14, 8), eligorgf4 = c(1, 1, 0, 0, 0), eligorgf8 = c(1, 
1, 0, 0, 0), eligorgm4 = c(1, 1, 1, 1, 1), eligorgm8 = c(1, 
1, 1, 1, 1), employedf1 = c(1, 1, 0, 0, 0), employedf2 = c(1, 
1, 0, 0, 0), employedf3 = c(1, 1, 0, 0, 0), employedf4 = c(1, 
1, 0, 0, 0), employedf5 = c(1, 1, 0, 0, 0), employedf6 = c(1, 
1, 0, 0, 0), employedf7 = c(1, 1, 0, 0, 0), employedf8 = c(1, 
1, 0, 0, 0), employedm1 = c(0, 1, 1, 1, 0), employedm2 = c(1, 
1, 1, 1, 1), employedm3 = c(1, 1, 0, 1, 1), employedm4 = c(1, 
1, 1, 1, 1), employedm5 = c(1, 1, 1, 1, 1), employedm6 = c(1, 
1, 1, 1, 1), employedm7 = c(1, 1, 1, 1, 1), employedm8 = c(1, 
1, 1, 1, 1), empstat3f1 = c(1, 1, 3, 3, 3), empstat3f2 = c(1, 
1, 3, 3, 3), empstat3f3 = c(1, 1, 3, 3, 3), empstat3f4 = c(1, 
1, 3, 3, 3), empstat3f5 = c(1, 1, 3, 3, 3), empstat3f6 = c(1, 
1, 3, 3, 3), empstat3f7 = c(1, 1, 3, 3, 3), empstat3f8 = c(1, 
1, 3, 3, 3), empstat3m1 = c(1, 1, 1, 1, 2), empstat3m2 = c(1, 
1, 1, 1, 1), empstat3m3 = c(1, 1, 1, 1, 1), empstat3m4 = c(1, 
1, 1, 1, 1), empstat3m5 = c(1, 1, 1, 1, 1), empstat3m6 = c(1, 
1, 1, 1, 1), empstat3m7 = c(1, 1, 1, 1, 1), empstat3m8 = c(1, 
1, 1, 1, 1), empstatf1 = c(10, 10, 33, 31, 31), empstatf2 = c(10, 
10, 33, 31, 31), empstatf3 = c(10, 10, 33, 31, 31), empstatf4 = c(10, 
10, 33, 31, 31), empstatf5 = c(10, 10, 31, 31, 31), empstatf6 = c(10, 
10, 31, 31, 31), empstatf7 = c(10, 10, 31, 31, 31), empstatf8 = c(10, 
10, 31, 31, 31), empstatm1 = c(12, 10, 10, 10, 20), empstatm2 = c(10, 
10, 10, 10, 10), empstatm3 = c(10, 10, 12, 10, 10), empstatm4 = c(10, 
10, 10, 10, 10), empstatm5 = c(10, 10, 10, 10, 10), empstatm6 = c(10, 
10, 10, 10, 10), empstatm7 = c(10, 10, 10, 10, 10), empstatm8 = c(10, 
10, 10, 10, 10), hispanf1 = c(0, 0, 0, 0, 0), hispanf2 = c(0, 
0, 0, 0, 0), hispanf3 = c(0, 0, 0, 0, 0), hispanf4 = c(0, 
0, 0, 0, 0), hispanf5 = c(0, 0, 0, 0, 0), hispanf6 = c(0, 
0, 0, 0, 0), hispanf7 = c(0, 0, 0, 0, 0), hispanf8 = c(0, 
0, 0, 0, 0), hispanm1 = c(0, 0, 901, 0, 0), hispanm2 = c(0, 
0, 901, 0, 0), hispanm3 = c(0, 0, 901, 0, 0), hispanm4 = c(0, 
0, 901, 0, 0), hispanm5 = c(0, 0, 901, 0, 0), hispanm6 = c(0, 
0, 901, 0, 0), hispanm7 = c(0, 0, 901, 0, 0), hispanm8 = c(0, 
0, 901, 0, 0), hourwagef4 = c(999.99, 999.99, 999.99, 999.99, 
999.99), hourwagef8 = c(999.99, 999.99, 999.99, 999.99, 999.99
), hourwagem4 = c(999.99, 999.99, 999.99, 999.99, 999.99), 
hourwagem8 = c(999.99, 8, 999.99, 999.99, 999.99), hwtfinlf1 = c(NA, 
NA, NA, NA, NA), hwtfinlf2 = c(NA, NA, NA, NA, NA), hwtfinlf3 = c(NA, 
NA, NA, NA, NA), hwtfinlf4 = c(NA, NA, NA, NA, NA), hwtfinlf5 = c(NA, 
NA, NA, NA, NA), hwtfinlf6 = c(NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_), hwtfinlf7 = c(NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_), hwtfinlf8 = c(NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_), hwtfinlm1 = c(NA, NA, NA, NA, NA), hwtfinlm2 = c(NA, 
NA, NA, NA, NA), hwtfinlm3 = c(NA, NA, NA, NA, NA), hwtfinlm4 = c(NA, 
NA, NA, NA, NA), hwtfinlm5 = c(NA, NA, NA, NA, NA), hwtfinlm6 = c(NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_), hwtfinlm7 = c(NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_), hwtfinlm8 = c(NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_), labforcef1 = c(2, 
2, 1, 1, 1), labforcef2 = c(2, 2, 1, 1, 1), labforcef3 = c(2, 
2, 1, 1, 1), labforcef4 = c(2, 2, 1, 1, 1), labforcef5 = c(2, 
2, 1, 1, 1), labforcef6 = c(2, 2, 1, 1, 1), labforcef7 = c(2, 
2, 1, 1, 1), labforcef8 = c(2, 2, 1, 1, 1), labforcem1 = c(2, 
2, 2, 2, 2), labforcem2 = c(2, 2, 2, 2, 2), labforcem3 = c(2, 
2, 2, 2, 2), labforcem4 = c(2, 2, 2, 2, 2), labforcem5 = c(2, 
2, 2, 2, 2), labforcem6 = c(2, 2, 2, 2, 2), labforcem7 = c(2, 
2, 2, 2, 2), labforcem8 = c(2, 2, 2, 2, 2), marstf1 = c(1, 
1, 1, 1, 1), marstf2 = c(1, 1, 1, 1, 1), marstf3 = c(1, 1, 
1, 1, 1), marstf4 = c(1, 1, 1, 1, 1), marstf5 = c(1, 1, 1, 
1, 1), marstf6 = c(1, 1, 1, 1, 1), marstf7 = c(1, 1, 1, 1, 
1), marstf8 = c(1, 1, 1, 1, 1), marstm1 = c(1, 1, 1, 1, 1
), marstm2 = c(1, 1, 1, 1, 1), marstm3 = c(1, 1, 1, 1, 1), 
marstm4 = c(1, 1, 1, 1, 1), marstm5 = c(1, 1, 1, 1, 1), marstm6 = c(1, 
1, 1, 1, 1), marstm7 = c(1, 1, 1, 1, 1), marstm8 = c(1, 1, 
1, 1, 1), merge2 = c(3, 3, 3, 3, 3), merge3 = c(3, 3, 3, 
3, 3), merge4 = c(3, 3, 3, 3, 3), merge5 = c(3, 3, 3, 3, 
3), merge6 = c(3, 3, 3, 3, 3), merge7 = c(3, 3, 3, 3, 3), 
merge8 = c(3, 3, 3, 3, 3), month1 = c(10, 1, 9, 6, 7), month2 = c(11, 
2, 10, 7, 8), month3 = c(12, 3, 11, 8, 9), month4 = c(1, 
4, 12, 9, 10), month5 = c(10, 1, 9, 6, 7), month6 = c(11, 
2, 10, 7, 8), month7 = c(12, 3, 11, 8, 9), month8 = c(1, 
4, 12, 9, 10), nativityf1 = c(NA, NA, NA, NA, NA), nativityf2 = c(NA, 
NA, NA, NA, NA), nativityf3 = c(NA, NA, NA, NA, NA), nativityf4 = c(NA, 
NA, NA, NA, NA), nativityf5 = c(NA, NA, NA, NA, NA), nativityf6 = c(NA, 
NA, NA, NA, NA), nativityf7 = c(NA, NA, NA, NA, NA), nativityf8 = c(NA, 
NA, NA, NA, NA), nativitym1 = c(NA, NA, NA, NA, NA), nativitym2 = c(NA, 
NA, NA, NA, NA), nativitym3 = c(NA, NA, NA, NA, NA), nativitym4 = c(NA, 
NA, NA, NA, NA), nativitym5 = c(NA, NA, NA, NA, NA), nativitym6 = c(NA, 
NA, NA, NA, NA), nativitym7 = c(NA, NA, NA, NA, NA), nativitym8 = c(NA, 
NA, NA, NA, NA), nchild1 = c(0, 2, 2, 3, 3), nchild2 = c(0, 
2, 2, 3, 3), nchild3 = c(0, 2, 2, 3, 3), nchild4 = c(0, 2, 
2, 3, 3), nchild5 = c(0, 2, 2, 3, 3), nchild6 = c(0, 2, 2, 
3, 3), nchild7 = c(0, 2, 2, 3, 3), nchild8 = c(0, 2, 2, 3, 
3), nchild_cohab01 = c(0, 0, 0, 0, 1), nchild_cohab02 = c(0, 
0, 0, 0, 0), nchild_cohab03 = c(0, 0, 0, 0, 0), nchild_cohab04 = c(0, 
0, 0, 0, 0), nchild_cohab05 = c(0, 0, 0, 0, 0), nchild_cohab06 = c(0, 
0, 0, 0, 0), nchild_cohab07 = c(0, 0, 0, 0, 0), nchild_cohab08 = c(0, 
0, 0, 0, 0), occ1990f1 = c(255, 458, 999, 999, 999), occ1990f2 = c(255, 
458, 999, 999, 999), occ1990f3 = c(255, 458, 999, 999, 999
), occ1990f4 = c(255, 458, 999, 335, 785), occ1990f5 = c(8, 
458, 999, 999, 999), occ1990f6 = c(25, 458, 999, 999, 999
), occ1990f7 = c(25, 458, 999, 999, 999), occ1990f8 = c(25, 
458, 999, 378, 999), occ1990m1 = c(22, 527, 213, 22, 628), 
occ1990m2 = c(22, 527, 549, 213, 274), occ1990m3 = c(64, 
527, 213, 213, 274), occ1990m4 = c(22, 527, 549, 779, 274
), occ1990m5 = c(22, 527, 549, 693, 274), occ1990m6 = c(22, 
527, 549, 693, 274), occ1990m7 = c(22, 527, 213, 213, 274
), occ1990m8 = c(22, 527, 213, 549, 274), occseif4 = c(55.9199981689453, 
26.0200004577637, NA, 26.0900001525879, 22.8600006103516), 
occseif5 = c(47.4300003051758, 26.0200004577637, NA, NA, 
NA), occseif8 = c(49.0900001525879, 26.0200004577637, NA, 
32.4700012207031, NA), occseim4 = c(47.5999984741211, 41.9099998474121, 
30.2600002288818, 23.6000003814697, 29.9599990844727), occseim5 = c(47.5999984741211, 
41.9099998474121, 30.2600002288818, 34.9599990844727, 29.9599990844727
), occseim8 = c(47.5999984741211, 41.9099998474121, 46.8400001525879, 
30.2600002288818, 29.9599990844727), oneracef1 = c(4, 4, 
4, 4, 4), oneracef2 = c(4, 4, 4, 4, 4), oneracef3 = c(4, 
4, 4, 4, 4), oneracef4 = c(4, 4, 4, 4, 4), oneracef5 = c(4, 
4, 4, 4, 4), oneracef6 = c(4, 4, 4, 4, 4), oneracef7 = c(4, 
4, 4, 4, 4), oneracef8 = c(4, 4, 4, 4, 4), oneracem1 = c(4, 
4, 4, 4, 4), oneracem2 = c(4, 4, 4, 4, 4), oneracem3 = c(4, 
4, 4, 4, 4), oneracem4 = c(4, 4, 4, 4, 4), oneracem5 = c(4, 
4, 4, 4, 4), oneracem6 = c(4, 4, 4, 4, 4), oneracem7 = c(4, 
4, 4, 4, 4), oneracem8 = c(4, 4, 4, 4, 4), orgeligf4 = c(1, 
1, 3, 3, 3), orgeligf8 = c(1, 1, 3, 3, 3), orgeligm4 = c(1, 
1, 1, 1, 1), orgeligm8 = c(1, 1, 1, 1, 1), paidhourf4 = c(1, 
1, 0, 0, 0), paidhourf8 = c(1, 1, 0, 0, 0), paidhourm4 = c(1, 
1, 1, 1, 1), paidhourm8 = c(1, 2, 1, 1, 1), racef1 = c(100, 
100, 100, 100, 100), racef2 = c(100, 100, 100, 100, 100), 
racef3 = c(100, 100, 100, 100, 100), racef4 = c(100, 100, 
100, 100, 100), racef5 = c(100, 100, 100, 100, 100), racef6 = c(100, 
100, 100, 100, 100), racef7 = c(100, 100, 100, 100, 100), 
racef8 = c(100, 100, 100, 100, 100), racem1 = c(100, 100, 
100, 100, 100), racem2 = c(100, 100, 100, 100, 100), racem3 = c(100, 
100, 100, 100, 100), racem4 = c(100, 100, 100, 100, 100), 
racem5 = c(100, 100, 100, 100, 100), racem6 = c(100, 100, 
100, 100, 100), racem7 = c(100, 100, 100, 100, 100), racem8 = c(100, 
100, 100, 100, 100), racematch_mis12 = c(1, 1, 1, 1, 1), 
racematch_mis14 = c(1, 1, 1, 1, 1), racematch_mis15 = c(1, 
1, 1, 1, 1), racematch_mis23 = c(1, 1, 1, 1, 1), racematch_mis26 = c(1, 
1, 1, 1, 1), racematch_mis34 = c(1, 1, 1, 1, 1), racematch_mis37 = c(1, 
1, 1, 1, 1), racematch_mis45 = c(1, 1, 1, 1, 1), racematch_mis48 = c(1, 
1, 1, 1, 1), racematch_mis56 = c(1, 1, 1, 1, 1), racematch_mis58 = c(1, 
1, 1, 1, 1), racematch_mis67 = c(1, 1, 1, 1, 1), racematch_mis78 = c(1, 
1, 1, 1, 1), racematchf_mis12 = c(1, 1, 1, 1, 1), racematchf_mis14 = c(1, 
1, 1, 1, 1), racematchf_mis15 = c(1, 1, 1, 1, 1), racematchf_mis23 = c(1, 
1, 1, 1, 1), racematchf_mis26 = c(1, 1, 1, 1, 1), racematchf_mis34 = c(1, 
1, 1, 1, 1), racematchf_mis37 = c(1, 1, 1, 1, 1), racematchf_mis45 = c(1, 
1, 1, 1, 1), racematchf_mis48 = c(1, 1, 1, 1, 1), racematchf_mis56 = c(1, 
1, 1, 1, 1), racematchf_mis58 = c(1, 1, 1, 1, 1), racematchf_mis67 = c(1, 
1, 1, 1, 1), racematchf_mis78 = c(1, 1, 1, 1, 1), racematchm_mis12 = c(1, 
1, 1, 1, 1), racematchm_mis14 = c(1, 1, 1, 1, 1), racematchm_mis15 = c(1, 
1, 1, 1, 1), racematchm_mis23 = c(1, 1, 1, 1, 1), racematchm_mis26 = c(1, 
1, 1, 1, 1), racematchm_mis34 = c(1, 1, 1, 1, 1), racematchm_mis37 = c(1, 
1, 1, 1, 1), racematchm_mis45 = c(1, 1, 1, 1, 1), racematchm_mis48 = c(1, 
1, 1, 1, 1), racematchm_mis56 = c(1, 1, 1, 1, 1), racematchm_mis58 = c(1, 
1, 1, 1, 1), racematchm_mis67 = c(1, 1, 1, 1, 1), racematchm_mis78 = c(1, 
1, 1, 1, 1), relatef1 = c(201, 201, 201, 201, 201), relatef2 = c(201, 
201, 201, 201, 201), relatef3 = c(201, 201, 201, 201, 201
), relatef4 = c(201, 201, 201, 201, 201), relatef5 = c(201, 
201, 201, 201, 201), relatef6 = c(201, 201, 201, 201, 201
), relatef7 = c(201, 201, 201, 201, 201), relatef8 = c(201, 
201, 201, 201, 201), relatem1 = c(101, 101, 101, 101, 101
), relatem2 = c(101, 101, 101, 101, 101), relatem3 = c(101, 
101, 101, 101, 101), relatem4 = c(101, 101, 101, 101, 101
), relatem5 = c(101, 101, 101, 101, 101), relatem6 = c(101, 
101, 101, 101, 101), relatem7 = c(101, 101, 101, 101, 101
), relatem8 = c(101, 101, 101, 101, 101), sef1 = c(0, 0, 
0, 0, 0), sef2 = c(0, 0, 0, 0, 0), sef3 = c(0, 0, 0, 0, 0
), sef4 = c(0, 0, 0, 0, 0), sef5 = c(0, 0, 0, 0, 0), sef6 = c(0, 
0, 0, 0, 0), sef7 = c(0, 0, 0, 0, 0), sef8 = c(0, 0, 0, 0, 
0), sem1 = c(0, 0, 0, 0, 0), sem2 = c(0, 0, 0, 0, 0), sem3 = c(0, 
0, 0, 0, 0), sem4 = c(0, 0, 0, 0, 0), sem5 = c(0, 0, 0, 0, 
0), sem6 = c(0, 0, 0, 0, 0), sem7 = c(0, 0, 0, 0, 0), sem8 = c(0, 
0, 0, 0, 0), statecensus1 = c(63, 63, 63, 63, 63), statecensus2 = c(63, 
63, 63, 63, 63), statecensus3 = c(63, 63, 63, 63, 63), statecensus4 = c(63, 
63, 63, 63, 63), statecensus5 = c(63, 63, 63, 63, 63), statecensus6 = c(63, 
63, 63, 63, 63), statecensus7 = c(63, 63, 63, 63, 63), statecensus8 = c(63, 
63, 63, 63, 63), statefip1 = c(1, 1, 1, 1, 1), statefip2 = c(1, 
1, 1, 1, 1), statefip3 = c(1, 1, 1, 1, 1), statefip4 = c(1, 
1, 1, 1, 1), statefip5 = c(1, 1, 1, 1, 1), statefip6 = c(1, 
1, 1, 1, 1), statefip7 = c(1, 1, 1, 1, 1), statefip8 = c(1, 
1, 1, 1, 1), wkstatf1 = c(10, 40, 99, 99, 99), wkstatf2 = c(10, 
40, 99, 99, 99), wkstatf3 = c(10, 40, 99, 99, 99), wkstatf4 = c(10, 
40, 99, 99, 99), wkstatf5 = c(10, 40, 99, 99, 99), wkstatf6 = c(10, 
40, 99, 99, 99), wkstatf7 = c(10, 40, 99, 99, 99), wkstatf8 = c(10, 
40, 99, 99, 99), wkstatm1 = c(10, 10, 10, 10, 50), wkstatm2 = c(10, 
10, 10, 10, 10), wkstatm3 = c(10, 10, 10, 10, 10), wkstatm4 = c(10, 
10, 10, 10, 10), wkstatm5 = c(10, 10, 10, 10, 10), wkstatm6 = c(10, 
10, 10, 10, 10), wkstatm7 = c(10, 10, 10, 10, 10), wkstatm8 = c(10, 
10, 10, 10, 10), workhr82onf4 = c(48, 31, NA, NA, NA), workhr82onf8 = c(50, 
30, NA, NA, NA), workhr82onm4 = c(55, 40, 50, 40, 40), workhr82onm8 = c(50, 
40, 50, 40, 40), workhrstatf4 = c(2, 2, NA, NA, NA), workhrstatf8 = c(2, 
2, NA, NA, NA), workhrstatm4 = c(2, 2, 2, 2, 2), workhrstatm8 = c(2, 
2, 2, 2, 2), wtfinlf1 = c(1248.33, 1645.3, 1203.23, 1865.61, 
1627.45), wtfinlf2 = c(1215.05, 1596.28, 1217.46, 1747.14, 
1650.36), wtfinlf3 = c(1172.82, 1529.08, 1234.23, 1697.64, 
1643.33), wtfinlf4 = c(1222.32, 1561.24, 1262.02, 1713.75, 
1629.05), wtfinlf5 = c(1246.07, 1634.79, 1362.73, 1711.81, 
1755.17), wtfinlf6 = c(1152.83, 1608.75, 1324.33, 1803.66, 
1714.75), wtfinlf7 = c(1136.77, 1637.46, 1370.03, 1800.02, 
1676.36), wtfinlf8 = c(1137.51, 1588.51, 1330.85, 1724.09, 
1672.49), wtfinlm1 = c(1261.05, 1619.03, 1137.29, 1809.66, 
1662.88), wtfinlm2 = c(1253.09, 1546.53, 1136.03, 1690.93, 
1623.57), wtfinlm3 = c(1240.68, 1521.73, 1132.68, 1685.34, 
1624.7), wtfinlm4 = c(1291.34, 1600.24, 1138.42, 1663.67, 
1650.74), wtfinlm5 = c(1373.51, 1656.8, 1301.41, 1681.27, 
1692.17), wtfinlm6 = c(1291.52, 1632.47, 1211.39, 1741.93, 
1692.3), wtfinlm7 = c(1274.28, 1646.33, 1260.04, 1725.2, 
1683.03), wtfinlm8 = c(1289.2, 1637.42, 1211.85, 1630.6, 
1722.4), year1 = c(1982, 1982, 1982, 1982, 1982), year2 = c(1982, 
1982, 1982, 1982, 1982), year3 = c(1982, 1982, 1982, 1982, 
1982), year4 = c(1983, 1982, 1982, 1982, 1982), year5 = c(1983, 
1983, 1983, 1983, 1983), year6 = c(1983, 1983, 1983, 1983, 
1983), year7 = c(1983, 1983, 1983, 1983, 1983), year8 = c(1984, 
1983, 1983, 1983, 1983), yngch1 = c(99, 7, 11, 8, 0), yngch2 = c(99, 
7, 11, 8, 1), yngch3 = c(99, 7, 11, 8, 1), yngch4 = c(99, 
7, 11, 2, 1), yngch5 = c(99, 8, 12, 9, 1), yngch6 = c(99, 
8, 12, 9, 2), yngch7 = c(99, 8, 12, 9, 2), yngch8 = c(99, 
8, 12, 10, 2), firstchild2 = c(0, 0, 0, 0, 0), earnweek4f_zero = c(670, 
225, 0, 0, 0), earnweek8f_zero = c(275, 240, 0, 0, 0), earnweek4m_zero = c(960, 
550, 432, 450, 150), earnweek8m_zero = c(460, 320, 680, 650, 
150), her_share_earnings4 = c(0.411042958498001, 0.290322571992874, 
0, 0, 0), her_share_earnings8 = c(0.374149650335312, 0.428571432828903, 
0, 0, 0), her_share_change = c(-0.0368933081626892, 0.138248860836029, 
0, 0, 0), his_share_earnings4 = c(0.588957071304321, 0.709677398204803, 
1, 1, 1), his_share_earnings8 = c(0.625850319862366, 0.571428596973419, 
1, 1, 1), his_share_change = c(0.0368932485580444, -0.138248801231384, 
0, 0, 0), has_empinfo = c(1, 1, 1, 1, 1), state = c("Alabama", 
"Alabama", "Alabama", "Alabama", "Alabama"), unemp4 = c(14.1000003814697, 
14.1000003814697, 14.1000003814697, 14.1000003814697, 14.1000003814697
), xyear = c(1983, 1983, 1983, 1983, 1983), unemp8 = c(13.8000001907349, 
13.8000001907349, 13.8000001907349, 13.8000001907349, 13.8000001907349
), fem_unemp4 = c(14.8999996185303, 14.8999996185303, 14.8999996185303, 
14.8999996185303, 14.8999996185303), fem_unemp8 = c(13.6000003814697, 
13.6000003814697, 13.6000003814697, 13.6000003814697, 13.6000003814697
), gap4 = c(0.899999976158142, 0.899999976158142, 0.899999976158142, 
0.899999976158142, 0.899999976158142), gap8 = c(-0.100000001490116, 
-0.100000001490116, -0.100000001490116, -0.100000001490116, 
-0.100000001490116), `_merge` = c(3, 3, 3, 3, 3), state2 = c(1, 
1, 1, 1, 1), recession = c(1, 1, 1, 1, 1), recession1 = c(0, 
0, 0, 0, 0), recession2 = c(0, 0, 0, 0, 0), recession3 = c(0, 
0, 0, 0, 0), firstyear_string = c(1982, 1982, 1982, 1982, 
1982), firstmonth_string = c(10, 1, 9, 6, 7), deflation_factor = c(2.58600818833163, 
2.58600818833163, 2.58600818833163, 2.58600818833163, 2.58600818833163
), wg_defl_f4 = c(1732.62548618219, 581.851842374616, NA, 
NA, NA), wg_diff = c(729.563766632549, 507.522620266121, 
1078.48556806551, 1030.90532241556, 237.901228249744), wg_defl_f8 = c(711.152251791198, 
620.641965199591, NA, NA, NA), wg_defl_m4 = c(2482.56786079836, 
1422.3045035824, 1117.15553735926, 1163.70368474923, 387.901228249744
), wg_defl_m8 = c(1189.56376663255, 827.522620266121, 1758.48556806551, 
1680.90532241556, 387.901228249744), wg_defl_f4_zero = c(1732.62548618219, 
581.851842374616, 0, 0, 0), wg_defl_f8_zero = c(711.152251791198, 
620.641965199591, 0, 0, 0), wg_defl_m4_zero = c(2482.56786079836, 
1422.3045035824, 1117.15553735926, 1163.70368474923, 387.901228249744
), wg_defl_m8_zero = c(1189.56376663255, 827.522620266121, 
1758.48556806551, 1680.90532241556, 387.901228249744), her_share_earnings4_adj = c(0.411042944785276, 
0.290322580645161, 0, 0, 0), her_share_earnings8_adj = c(0.374149659863946, 
0.428571428571429, 0, 0, 0), her_share_earnings_change_adj = c(-0.0368932849213305, 
0.138248847926267, 0, 0, 0), his_share_earnings4_adj = c(0.588957055214724, 
0.709677419354839, 1, 1, 1), his_share_earnings8_adj = c(0.625850340136054, 
0.571428571428571, 1, 1, 1), his_share_earnings_change_adj = c(0.0368932849213305, 
-0.138248847926267, 0, 0, 0), workhr4f_zero = c(48, 31, 0, 
0, 0), workhr8f_zero = c(50, 30, 0, 0, 0), workhr4m_zero = c(55, 
40, 50, 40, 40), workhr8m_zero = c(50, 40, 50, 40, 40), her_share_hours4 = c(0.466019417475728, 
0.436619718309859, 0, 0, 0), her_share_hours8 = c(0.5, 0.428571428571429, 
0, 0, 0), her_share_change_hours = c(0.0339805825242718, 
-0.00804828973843058, 0, 0, 0), his_share_hours4 = c(0.533980582524272, 
0.563380281690141, 1, 1, 1), his_share_hours8 = c(0.5, 0.571428571428571, 
1, 1, 1), his_share_change_hours = c(-0.0339805825242718, 
0.00804828973843053, 0, 0, 0), firstyear_factor = c("1982", 
"1982", "1982", "1982", "1982"), year_recode = c(0, 0, 0, 
0, 0), firstyear_character = c("1982", "1982", "1982", "1982", 
"1982")), row.names = c(NA, -5L), class = c("tbl_df", "tbl", "data.frame"))
...