The following object(s) are masked _by_ '.GlobalEnv': one idhous99 idpers_new idpers sex birthm birthy age marstat intlang 1 251 1 25101 woman 12 1942 57 widower/widow french 2 261 1 26101 woman 11 1922 77 married german 3 301 1 30101 woman 4 1926 73 widower/widow french 4 371 1 37101 woman 8 1919 80 widower/widow french 5 491 1 49102 woman 9 1950 49 married french 6 571 1 57102 woman 6 1953 46 married french religion disrel_fr health work 1 roman catholic well not in labor force 2 roman catholic so, so (average) not in labor force 3 protestant or reformed church never well not in labor force 4 no denomination or religion never well not in labor force 5 roman catholic never so, so (average) active occupied 6 roman catholic well active occupied income cfriends reading inchild children 1 -2 4 at least once a week 0 0 2 -2 NA at least once a week 0 0 3 -2 8 at least once a week 0 2 4 45360 4 every day 0 0 5 66521 1 never 0 2 6 -2 12 every day 3 3 education fath15 moth15 stepp15 sib15 inear15 1 compulsory or less yes yes inapplicable yes both equally 2 compulsory or less yes yes inapplicable yes father 3 compulsory or less yes yes inapplicable yes both equally 4 high school yes yes inapplicable yes father 5 compulsory or less yes yes inapplicable yes father 6 general training/apprenticeship yes yes inapplicable yes father ocfath15 edufath mthwk15 1 skilled agricultural and fishery workers 1 yes 2 legislators, senior officials, managers 1 no 3 skilled agricultural and fishery workers 0 no 4 professionals 15 no 5 craft and related trades workers 3 no 6 technicians and associate professionals 1 yes ocmth15 edumth 1 elementary occupations compulsory or less 2 compulsory or less 3 skilled agricultural and fishery workers compulsory or less 4 compulsory or less 5 compulsory or less 6 elementary occupations general training/apprenticeship fathpol mothpol one 1 innaplicable innaplicable 1 2 don't know don't know 1 3 innaplicable innaplicable 1 4 right right 1 5 right centre 1 6 innaplicable innaplicable 1 [1] "idhous99" "idpers_new" "idpers" "sex" "birthm" [6] "birthy" "age" "marstat" "intlang" "religion" [11] "disrel_fr" "health" "work" "income" "cfriends" [16] "reading" "inchild" "children" "education" "fath15" [21] "moth15" "stepp15" "sib15" "inear15" "ocfath15" [26] "edufath" "mthwk15" "ocmth15" "edumth" "fathpol" [31] "mothpol" "one" idhous99 idpers_new idpers sex Min. : 251 Min. :1 Min. : 25101 other error : 0 1st Qu.: 37136 1st Qu.:1 1st Qu.: 3713601 filter error : 0 Median : 71891 Median :1 Median : 7189101 inapplicable : 0 Mean : 72703 Mean :1 Mean : 7270324 no answer : 0 3rd Qu.:108539 3rd Qu.:1 3rd Qu.:10853851 does not know: 0 Max. :146601 Max. :1 Max. :14660101 man : 0 woman :1878 birthm birthy age marstat Min. :-3.000 Min. :1908 Min. :45.00 married :1192 1st Qu.: 3.000 1st Qu.:1933 1st Qu.:50.00 widower/widow : 292 Median : 6.000 Median :1942 Median :57.00 divorced : 212 Mean : 6.381 Mean :1940 Mean :58.74 single, never married: 147 3rd Qu.: 9.000 3rd Qu.:1949 3rd Qu.:66.00 separated : 35 Max. :12.000 Max. :1954 Max. :91.00 other error : 0 (Other) : 0 intlang religion disrel_fr german :1269 protestant or reformed church:858 :1328 french : 515 roman catholic :763 never : 476 italian : 94 no denomination or religion :157 very often : 54 other error : 0 other :100 other error : 0 filter error: 0 other error : 0 filter error: 0 inapplicable: 0 filter error : 0 (Other) : 0 (Other) : 0 (Other) : 0 NA's : 20 health work income well :947 not in labor force:974 Min. : -8 very well :453 active occupied :904 1st Qu.: -2 so, so (average):417 other error : 0 Median : -2 not very well : 55 filter error : 0 Mean : 18717 not well at all : 6 inapplicable : 0 3rd Qu.: 26050 other error : 0 no answer : 0 Max. :1170000 (Other) : 0 (Other) : 0 cfriends reading inchild Min. : 1.000 every day :1600 Min. :0.0000 1st Qu.: 2.000 at least once a week : 187 1st Qu.:0.0000 Median : 4.000 at least once a month : 31 Median :0.0000 Mean : 6.102 never : 30 Mean :0.5351 3rd Qu.: 8.000 less than once a month: 26 3rd Qu.:1.0000 Max. :60.000 (Other) : 0 Max. :6.0000 NA's :304 NA's : 4 children education fath15 Min. : 0.000 compulsory or less :441 other error : 0 1st Qu.: 1.000 domestic science course, 1yr :197 filter error : 0 Median : 2.000 general training/apprenticeship:643 inapplicable : 0 Mean : 1.937 high school :483 no answer : 0 3rd Qu.: 3.000 university/postgrad :114 does not know: 0 Max. :11.000 yes :1645 no : 233 moth15 stepp15 sib15 other error : 0 other error : 0 other error : 0 filter error : 0 filter error : 0 filter error : 0 inapplicable : 0 inapplicable :1878 inapplicable : 0 no answer : 0 no answer : 0 no answer : 0 does not know: 0 does not know: 0 does not know: 0 yes :1731 yes : 0 yes :1616 no : 147 no : 0 no : 262 inear15 ocfath15 father :1353 craft and related trades workers :517 both equally: 276 skilled agricultural and fishery workers:368 mother : 156 elementary occupations :274 a guardian : 88 technicians and associate professionals :246 other error : 0 service workers, market sales workers :152 (Other) : 0 (Other) :248 NA's : 5 NA's : 73 edufath mthwk15 Min. : 0.000 no :1120 1st Qu.: 1.000 yes : 734 Median : 3.000 other error : 0 Mean : 3.451 filter error: 0 3rd Qu.: 3.000 inapplicable: 0 Max. :15.000 (Other) : 0 NA's :151 NA's : 24 ocmth15 service workers, market sales workers :387 craft and related trades workers :222 elementary occupations :215 technicians and associate professionals :167 skilled agricultural and fishery workers:156 (Other) : 76 NA's :655 edumth fathpol compulsory or less :967 other :156 general training/apprenticeship:548 don't know :158 high school :239 innaplicable :621 university/postgrad : 22 no particular tendency:144 NA's :102 left :178 centre :306 right :315 mothpol one other :158 Min. :1 don't know :153 1st Qu.:1 innaplicable :621 Median :1 no particular tendency:405 Mean :1 left :102 3rd Qu.:1 centre :236 Max. :1 right :203 idhous99 idpers_new idpers sex birthm birthy 7.270323e+04 1.000000e+00 7.270324e+06 NA 6.381257e+00 1.940257e+03 age marstat intlang religion disrel_fr health 5.874334e+01 NA NA NA NA NA work income cfriends reading inchild children NA 1.871670e+04 NA NA 5.351438e-01 1.937167e+00 education fath15 moth15 stepp15 sib15 inear15 NA NA NA NA NA NA ocfath15 edufath mthwk15 ocmth15 edumth fathpol NA NA NA NA NA NA mothpol one NA 1.000000e+00 idhous99 idpers_new idpers sex birthm birthy 4.217101e+04 0.000000e+00 4.217101e+06 NA 3.491143e+00 1.035302e+01 age marstat intlang religion disrel_fr health 1.035302e+01 NA NA NA NA NA work income cfriends reading inchild children NA 4.259607e+04 NA NA 9.293591e-01 1.452605e+00 education fath15 moth15 stepp15 sib15 inear15 NA NA NA NA NA NA ocfath15 edufath mthwk15 ocmth15 edumth fathpol NA NA NA NA NA NA mothpol one NA 0.000000e+00 children 0 1 2 3 4 5 6 7 8 9 10 11 379 262 684 353 128 35 16 8 10 1 1 1 children 0 1 2 3 4 5 0.2018104366 0.1395101171 0.3642172524 0.1879659212 0.0681576145 0.0186368477 6 7 8 9 10 11 0.0085197018 0.0042598509 0.0053248136 0.0005324814 0.0005324814 0.0005324814 Call: glm(formula = children ~ one, family = poisson(), data = data.ch06p3a) Deviance Residuals: Min 1Q Median 3Q Max -1.9683 -0.7429 0.0449 0.7062 4.4812 Coefficients: (1 not defined because of singularities) Estimate Std. Error z value Pr(>|z|) (Intercept) 0.66123 0.01658 39.88 <2e-16 *** one NA NA NA NA --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 2428.1 on 1877 degrees of freedom Residual deviance: 2428.1 on 1877 degrees of freedom AIC: 6478.7 Number of Fisher Scoring iterations: 5 Call: glm.nb(formula = children ~ one, data = data.ch06p3a, init.theta = 22.83989503, link = log) Deviance Residuals: Min 1Q Median 3Q Max -1.9284 -0.7183 0.0431 0.6736 4.1364 Coefficients: (1 not defined because of singularities) Estimate Std. Error z value Pr(>|z|) (Intercept) 0.66123 0.01727 38.29 <2e-16 *** one NA NA NA NA --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for Negative Binomial(22.8399) family taken to be 1) Null deviance: 2270.8 on 1877 degrees of freedom Residual deviance: 2270.8 on 1877 degrees of freedom AIC: 6474.2 Number of Fisher Scoring iterations: 1 Theta: 22.84 Std. Err.: 9.64 2 x log-likelihood: -6470.15 0 1 2 3 4 5 1.557684e-01 2.781575e-01 2.592283e-01 1.678140e-01 8.475713e-02 3.557168e-02 6 7 8 9 10 11 1.290440e-02 4.156721e-03 1.212203e-03 3.247605e-04 8.084480e-05 1.887029e-05 Call: glm(formula = children ~ age, family = poisson(), data = data.ch06p3a) Deviance Residuals: Min 1Q Median 3Q Max -2.0631 -0.7534 0.0568 0.6703 4.3733 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.488717 0.095143 5.137 2.8e-07 *** age 0.002929 0.001586 1.846 0.0648 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 2428.1 on 1877 degrees of freedom Residual deviance: 2424.7 on 1876 degrees of freedom AIC: 6477.3 Number of Fisher Scoring iterations: 5 Call: glm.nb(formula = children ~ age, data = data.ch06p3a, init.theta = 23.33720034, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.0184 -0.7289 0.0546 0.6395 4.0348 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.488324 0.099048 4.930 8.22e-07 *** age 0.002935 0.001652 1.777 0.0756 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for Negative Binomial(23.3372) family taken to be 1) Null deviance: 2273.9 on 1877 degrees of freedom Residual deviance: 2270.7 on 1876 degrees of freedom AIC: 6473 Number of Fisher Scoring iterations: 1 Theta: 23.3 Std. Err.: 10.0 2 x log-likelihood: -6467.02 Call: flexmix(formula = children ~ age, data = data.ch06p3a, k = 2, model = FLXMRglm(family = "poisson")) prior size post>0 ratio Comp.1 0.139 176 1657 0.106 Comp.2 0.861 1702 1878 0.906 'log Lik.' -3201.323 (df=5) AIC: 6412.645 BIC: 6440.335 Call: hurdle(formula = children ~ age, data = data.ch06p3a, dist = "poisson") Pearson residuals: Min 1Q Median 3Q Max -1.3456 -0.6181 0.0529 0.6234 5.8704 Count model coefficients (truncated poisson with log link): Estimate Std. Error z value Pr(>|z|) (Intercept) 0.201351 0.113481 1.774 0.076 . age 0.009501 0.001874 5.071 3.95e-07 *** Zero hurdle model coefficients (binomial with logit link): Estimate Std. Error z value Pr(>|z|) (Intercept) 2.436244 0.329243 7.400 1.37e-13 *** age -0.017892 0.005413 -3.305 0.000949 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Number of iterations in BFGS optimization: 7 Log-likelihood: -3185 on 4 Df Call: zeroinfl(formula = children ~ age, data = data.ch06p3a, dist = "poisson") Pearson residuals: Min 1Q Median 3Q Max -1.31295 -0.63277 0.03769 0.60155 5.87709 Count model coefficients (poisson with log link): Estimate Std. Error z value Pr(>|z|) (Intercept) 0.183790 0.108531 1.693 0.0904 . age 0.009777 0.001799 5.436 5.45e-08 *** Zero-inflation model coefficients (binomial with logit link): Estimate Std. Error z value Pr(>|z|) (Intercept) -6.31017 0.92629 -6.812 9.60e-12 *** age 0.06488 0.01307 4.965 6.89e-07 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Number of iterations in BFGS optimization: 14 Log-likelihood: -3183 on 4 Df Call: polr(formula = factor(children) ~ age, method = "probit") Coefficients: Value Std. Error t value age 0.002315 0.002341 0.9886 Intercepts: Value Std. Error t value 0|1 -0.6992 0.1414 -4.9455 1|2 -0.2737 0.1400 -1.9554 2|3 0.6758 0.1403 4.8153 3|4 1.3822 0.1437 9.6181 4|5 1.9084 0.1496 12.7601 5|6 2.1992 0.1561 14.0865 6|7 2.4245 0.1645 14.7421 7|8 2.6028 0.1743 14.9295 8|9 3.0921 0.2307 13.4058 9|10 3.2155 0.2569 12.5183 10|11 3.4174 0.3190 10.7118 Residual Deviance: 6162.133 AIC: 6186.133 logl.poiss.int logl.nb.int [1,] -3238.338 -3235.075 logl.poiss logl.nb logl.FMP2 logl.hurdle.poiss logl.zip.poiss [1,] -3236.645 -3233.51 -3201.323 -3184.836 -3182.777 logl.orderedprobit [1,] -3081.066 The following object(s) are masked _by_ '.GlobalEnv': one The following object(s) are masked from 'data.ch06p3a': age, children, one, sex pid sex dobm doby memorig plbornc race paju 1 10014578 female january 1937 original sample -8 white inapplicable 2 10016848 female august 1959 original sample -8 white inapplicable 3 10017933 female september 1942 original sample 11 white inapplicable 4 10020209 female august 1941 original sample -8 white inapplicable 5 10023526 female april 1953 original sample 37 chinese inapplicable 6 10029133 female november 1961 original sample -8 white inapplicable maju lprnt lnprnt ch1bm ch1by scend 1 mother not work not applicable -8 not applicable -8 15 2 inapplicable not applicable -8 not applicable -8 16 3 inapplicable not applicable -8 not applicable -8 17 4 mother not work not applicable -8 not applicable -8 15 5 mother not work not applicable -8 not applicable -8 17 6 mother not work not applicable -8 not applicable -8 19 bwtgm pagold pargsc 1 not applicable service class,higher grade managerial & technical occ 2 not applicable not applicable inapplicable 3 not applicable service class,higher grade professional occ 4 not applicable semi,unskilled manual workers partly skilled occ 5 not applicable service class,higher grade professional occ 6 not applicable service class,higher grade managerial & technical occ magold margsc j1gold j1rgsc 1 not applicable inapplicable inapplicable inapplicable 2 not applicable inapplicable inapplicable inapplicable 3 service class,lower grade managerial & technical occ inapplicable inapplicable 4 not applicable inapplicable inapplicable inapplicable 5 not applicable inapplicable inapplicable inapplicable 6 not applicable inapplicable inapplicable inapplicable fsource children ssex2 ssex3 twin1 twin2 twin3 age one 1 Wave 2 4 1 NA NA NA NA 72 1 2 Wave 2 1 NA NA NA NA NA 50 1 3 Wave 2 1 NA NA NA NA NA 67 1 4 Wave 2 2 NA NA NA NA NA 68 1 5 No children 0 NA NA NA NA NA 56 1 6 Wave 2 1 NA NA NA NA NA 48 1 children 0 1 2 3 4 5 6 7 8 9 10 11 1764 805 2160 1262 487 181 69 33 12 5 3 1 children 0 1 2 3 4 5 0.2601002654 0.1186965497 0.3184901209 0.1860808021 0.0718077263 0.0266882925 6 7 8 9 10 11 0.0101739900 0.0048658213 0.0017693896 0.0007372457 0.0004423474 0.0001474491 Call: glm(formula = children ~ one, family = poisson(), data = data.ch06p3a) Deviance Residuals: Min 1Q Median 3Q Max -1.9683 -0.7429 0.0449 0.7062 4.4812 Coefficients: (1 not defined because of singularities) Estimate Std. Error z value Pr(>|z|) (Intercept) 0.66123 0.01658 39.88 <2e-16 *** one NA NA NA NA --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 2428.1 on 1877 degrees of freedom Residual deviance: 2428.1 on 1877 degrees of freedom AIC: 6478.7 Number of Fisher Scoring iterations: 5 Call: glm.nb(formula = children ~ one, data = data.ch06p3a, init.theta = 22.83989503, link = log) Deviance Residuals: Min 1Q Median 3Q Max -1.9284 -0.7183 0.0431 0.6736 4.1364 Coefficients: (1 not defined because of singularities) Estimate Std. Error z value Pr(>|z|) (Intercept) 0.66123 0.01727 38.29 <2e-16 *** one NA NA NA NA --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for Negative Binomial(22.8399) family taken to be 1) Null deviance: 2270.8 on 1877 degrees of freedom Residual deviance: 2270.8 on 1877 degrees of freedom AIC: 6474.2 Number of Fisher Scoring iterations: 1 Theta: 22.84 Std. Err.: 9.64 2 x log-likelihood: -6470.15 0 1 2 3 4 5 1.557684e-01 2.781575e-01 2.592283e-01 1.678140e-01 8.475713e-02 3.557168e-02 6 7 8 9 10 11 1.290440e-02 4.156721e-03 1.212203e-03 3.247605e-04 8.084480e-05 1.887029e-05