## Learn more about sjPlot with 'browseVignettes("sjPlot")'.
data1$cntry <- droplevels(data1$cntry)
countrylist <- unique(data1$cntry)
countrylist
##  [1] Belgium            Bulgaria           Switzerland       
##  [4] Cyprus             Czech Republic     Germany           
##  [7] Denmark            Estonia            Spain             
## [10] Finland            France             United Kingdom    
## [13] Greece             Croatia            Hungary           
## [16] Ireland            Israel             Lithuania         
## [19] Netherlands        Norway             Poland            
## [22] Portugal           Russian Federation Sweden            
## [25] Slovenia           Slovakia           Ukraine           
## 27 Levels: Belgium Bulgaria Switzerland Cyprus Czech Republic ... Ukraine

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -11.666  -2.129  -0.635   1.153  37.061 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 43.62629    2.88635  15.115  < 2e-16 ***
## agean        0.27266    0.04934   5.526 7.09e-08 ***
## gndrFemale  -0.24593    0.53325  -0.461  0.64499    
## eduyrsn      0.19987    0.07004   2.853  0.00462 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.629 on 302 degrees of freedom
## Multiple R-squared:  0.1103, Adjusted R-squared:  0.1015 
## F-statistic: 12.48 on 3 and 302 DF,  p-value: 1.028e-07

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -18.2801  -2.0815   0.1733   2.0047   9.7488 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 49.84008    1.83626  27.142  < 2e-16 ***
## agean        0.15355    0.03082   4.983 9.47e-07 ***
## gndrFemale  -1.93186    0.32792  -5.891 8.35e-09 ***
## eduyrsn      0.12844    0.05402   2.378   0.0179 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.178 on 386 degrees of freedom
## Multiple R-squared:  0.1473, Adjusted R-squared:  0.1407 
## F-statistic: 22.23 on 3 and 386 DF,  p-value: 2.678e-13

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -14.8263  -2.0437  -0.1128   1.6540  23.5018 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 44.12210    2.30961  19.104  < 2e-16 ***
## agean        0.32165    0.03947   8.149 6.49e-15 ***
## gndrFemale  -0.38646    0.42641  -0.906    0.365    
## eduyrsn      0.05505    0.05925   0.929    0.353    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.995 on 352 degrees of freedom
## Multiple R-squared:  0.1635, Adjusted R-squared:  0.1564 
## F-statistic: 22.94 on 3 and 352 DF,  p-value: 1.376e-13

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -8.5627 -1.7305 -0.0237  1.9216 13.4665 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 47.15011    2.60930  18.070  < 2e-16 ***
## agean        0.28313    0.04479   6.321 2.53e-09 ***
## gndrFemale  -1.61650    0.55848  -2.894  0.00434 ** 
## eduyrsn      0.02648    0.05189   0.510  0.61060    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.513 on 158 degrees of freedom
## Multiple R-squared:  0.2313, Adjusted R-squared:  0.2167 
## F-statistic: 15.85 on 3 and 158 DF,  p-value: 4.671e-09

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -14.5340  -1.4201  -0.0333   1.9842  10.9568 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 55.00295    2.18067  25.223  < 2e-16 ***
## agean        0.10061    0.03479   2.892  0.00406 ** 
## gndrFemale  -2.28749    0.33681  -6.792  4.4e-11 ***
## eduyrsn     -0.01767    0.07044  -0.251  0.80205    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.217 on 372 degrees of freedom
## Multiple R-squared:  0.1426, Adjusted R-squared:  0.1357 
## F-statistic: 20.63 on 3 and 372 DF,  p-value: 2.203e-12

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -26.200  -1.768  -0.245   1.625  36.850 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 49.41990    1.69320  29.187  < 2e-16 ***
## agean        0.17420    0.02800   6.221 8.25e-10 ***
## gndrFemale  -0.73741    0.32274  -2.285   0.0226 *  
## eduyrsn      0.21630    0.04845   4.464 9.29e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.349 on 739 degrees of freedom
## Multiple R-squared:  0.07815,    Adjusted R-squared:  0.07441 
## F-statistic: 20.88 on 3 and 739 DF,  p-value: 5.385e-13

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -14.7269  -2.5103  -0.5869   1.7089  17.0249 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 56.81711    2.06011  27.580  < 2e-16 ***
## agean        0.11255    0.03480   3.235  0.00133 ** 
## gndrFemale  -0.32972    0.41333  -0.798  0.42554    
## eduyrsn      0.06782    0.04477   1.515  0.13063    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.921 on 372 degrees of freedom
## Multiple R-squared:  0.03607,    Adjusted R-squared:  0.02829 
## F-statistic:  4.64 on 3 and 372 DF,  p-value: 0.00337

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -47.920  -2.495   0.898   2.850  19.203 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 51.40504    2.82933  18.169   <2e-16 ***
## agean        0.08640    0.04489   1.925   0.0550 .  
## gndrFemale  -1.51054    0.60105  -2.513   0.0124 *  
## eduyrsn      0.21939    0.09537   2.300   0.0220 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5.653 on 371 degrees of freedom
## Multiple R-squared:  0.03437,    Adjusted R-squared:  0.02656 
## F-statistic: 4.402 on 3 and 371 DF,  p-value: 0.004653

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -19.8755  -2.1463  -0.1528   3.0371  19.5481 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 47.48115    2.90089  16.368  < 2e-16 ***
## agean        0.23186    0.05044   4.596 6.39e-06 ***
## gndrFemale  -0.72402    0.56458  -1.282   0.2007    
## eduyrsn      0.11526    0.05209   2.213   0.0277 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.802 on 294 degrees of freedom
## Multiple R-squared:  0.0753, Adjusted R-squared:  0.06587 
## F-statistic:  7.98 on 3 and 294 DF,  p-value: 3.933e-05

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -13.791  -2.156  -0.515   1.362  56.320 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 37.55696    2.90323  12.936  < 2e-16 ***
## agean        0.39602    0.04705   8.417 6.55e-16 ***
## gndrFemale  -1.08271    0.52436  -2.065 0.039570 *  
## eduyrsn      0.25356    0.06513   3.893 0.000115 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5.312 on 410 degrees of freedom
## Multiple R-squared:  0.1624, Adjusted R-squared:  0.1563 
## F-statistic:  26.5 on 3 and 410 DF,  p-value: 1.097e-15

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -29.475  -2.068   0.257   1.867  19.178 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 43.72436    2.31173  18.914  < 2e-16 ***
## agean        0.26893    0.04059   6.625  1.3e-10 ***
## gndrFemale  -0.26483    0.40764  -0.650  0.51634    
## eduyrsn      0.13586    0.05216   2.605  0.00959 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.831 on 352 degrees of freedom
## Multiple R-squared:  0.1213, Adjusted R-squared:  0.1138 
## F-statistic:  16.2 on 3 and 352 DF,  p-value: 6.945e-10

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -18.968  -3.608  -0.259   2.260  38.597 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 36.57153    2.71356   13.48   <2e-16 ***
## agean        0.44967    0.04314   10.42   <2e-16 ***
## gndrFemale  -0.17672    0.57003   -0.31    0.757    
## eduyrsn      0.05493    0.07524    0.73    0.466    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 6.355 on 496 degrees of freedom
## Multiple R-squared:  0.1815, Adjusted R-squared:  0.1766 
## F-statistic: 36.67 on 3 and 496 DF,  p-value: < 2.2e-16

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -56.807  -2.331   0.333   2.309  13.337 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 48.31559    3.98054  12.138  < 2e-16 ***
## agean        0.20285    0.06623   3.063  0.00243 ** 
## gndrFemale  -1.11950    0.67968  -1.647  0.10076    
## eduyrsn      0.05885    0.08339   0.706  0.48100    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5.443 on 257 degrees of freedom
## Multiple R-squared:  0.04912,    Adjusted R-squared:  0.03802 
## F-statistic: 4.425 on 3 and 257 DF,  p-value: 0.004705

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -18.9513  -2.0608   0.5036   2.3889  11.3018 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 45.98874    3.48655  13.190  < 2e-16 ***
## agean        0.22904    0.06211   3.687 0.000308 ***
## gndrFemale  -1.69461    0.71655  -2.365 0.019202 *  
## eduyrsn      0.14295    0.08481   1.685 0.093808 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.523 on 164 degrees of freedom
## Multiple R-squared:  0.1333, Adjusted R-squared:  0.1175 
## F-statistic: 8.411 on 3 and 164 DF,  p-value: 3.115e-05

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -19.8287  -1.7261   0.6962   1.8535  10.5491 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 51.83133    2.51349  20.621  < 2e-16 ***
## agean        0.10495    0.04591   2.286   0.0230 *  
## gndrFemale  -1.96896    0.43188  -4.559 7.74e-06 ***
## eduyrsn      0.12596    0.05740   2.194   0.0291 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.592 on 275 degrees of freedom
## Multiple R-squared:  0.09829,    Adjusted R-squared:  0.08845 
## F-statistic: 9.992 on 3 and 275 DF,  p-value: 2.847e-06

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -21.4760  -3.0490  -0.4352   2.7808  25.8713 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 49.15151    3.38066  14.539  < 2e-16 ***
## agean        0.31879    0.04763   6.693 1.15e-10 ***
## gndrFemale  -2.67723    0.67625  -3.959 9.50e-05 ***
## eduyrsn     -0.15087    0.09492  -1.589    0.113    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5.657 on 287 degrees of freedom
## Multiple R-squared:  0.2305, Adjusted R-squared:  0.2225 
## F-statistic: 28.66 on 3 and 287 DF,  p-value: 3.07e-16

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -19.3792  -3.7500  -0.4648   2.6034  28.6717 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 37.78814    4.01429   9.413  < 2e-16 ***
## agean        0.41815    0.06608   6.328 9.06e-10 ***
## gndrFemale  -1.13871    0.77122  -1.476    0.141    
## eduyrsn      0.08444    0.10906   0.774    0.439    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 6.62 on 299 degrees of freedom
## Multiple R-squared:  0.1302, Adjusted R-squared:  0.1215 
## F-statistic: 14.92 on 3 and 299 DF,  p-value: 4.462e-09

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -13.1460  -2.6930   0.3751   2.1788  12.0719 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 56.33367    2.26092  24.916  < 2e-16 ***
## agean        0.04264    0.03474   1.228    0.221    
## gndrFemale  -2.03622    0.48988  -4.157 4.48e-05 ***
## eduyrsn      0.08588    0.07645   1.123    0.262    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.608 on 243 degrees of freedom
## Multiple R-squared:  0.07316,    Adjusted R-squared:  0.06171 
## F-statistic: 6.393 on 3 and 243 DF,  p-value: 0.0003471

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -15.0420  -2.1772  -0.2718   2.5519  22.6969 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 52.17638    2.40598  21.686  < 2e-16 ***
## agean        0.18018    0.03961   4.548 7.29e-06 ***
## gndrFemale  -0.83889    0.44499  -1.885   0.0602 .  
## eduyrsn      0.06986    0.05431   1.286   0.1991    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.277 on 378 degrees of freedom
## Multiple R-squared:  0.07088,    Adjusted R-squared:  0.0635 
## F-statistic: 9.612 on 3 and 378 DF,  p-value: 3.952e-06

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -13.8995  -2.1904  -0.2429   1.6722  23.0337 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 50.89164    1.99809  25.470  < 2e-16 ***
## agean        0.18755    0.03325   5.640 3.19e-08 ***
## gndrFemale  -0.03126    0.38078  -0.082    0.935    
## eduyrsn      0.22886    0.04889   4.681 3.89e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.823 on 407 degrees of freedom
## Multiple R-squared:  0.1126, Adjusted R-squared:  0.106 
## F-statistic: 17.21 on 3 and 407 DF,  p-value: 1.542e-10

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -13.0407  -2.1205  -0.0577   2.0468  10.5916 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 46.22117    2.19520  21.056  < 2e-16 ***
## agean        0.24841    0.03787   6.559 2.79e-10 ***
## gndrFemale  -2.62034    0.43560  -6.016 5.88e-09 ***
## eduyrsn      0.09596    0.06681   1.436    0.152    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.508 on 267 degrees of freedom
## Multiple R-squared:  0.2597, Adjusted R-squared:  0.2514 
## F-statistic: 31.22 on 3 and 267 DF,  p-value: < 2.2e-16

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -12.7744  -1.7040   0.0457   2.5784  11.4373 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 51.75203    2.26334  22.865  < 2e-16 ***
## agean        0.17543    0.03827   4.584 7.34e-06 ***
## gndrFemale  -0.36219    0.48542  -0.746   0.4563    
## eduyrsn     -0.09563    0.05522  -1.732   0.0846 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.753 on 241 degrees of freedom
## Multiple R-squared:  0.1026, Adjusted R-squared:  0.09148 
## F-statistic: 9.189 on 3 and 241 DF,  p-value: 8.85e-06

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -22.504  -3.026   1.483   2.663  16.998 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 53.83873    3.09135  17.416  < 2e-16 ***
## agean        0.08359    0.04980   1.679   0.0941 .  
## gndrFemale  -3.86308    0.52590  -7.346 1.31e-12 ***
## eduyrsn     -0.09783    0.09875  -0.991   0.3225    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5.046 on 372 degrees of freedom
## Multiple R-squared:  0.1404, Adjusted R-squared:  0.1335 
## F-statistic: 20.25 on 3 and 372 DF,  p-value: 3.531e-12

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -14.774  -2.140   0.015   2.047  32.410 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 46.02483    2.26902  20.284  < 2e-16 ***
## agean        0.27220    0.03564   7.638 1.82e-13 ***
## gndrFemale   0.50978    0.44407   1.148    0.252    
## eduyrsn      0.10008    0.06472   1.546    0.123    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.289 on 379 degrees of freedom
## Multiple R-squared:  0.1359, Adjusted R-squared:  0.1291 
## F-statistic: 19.87 on 3 and 379 DF,  p-value: 5.524e-12

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -8.214 -2.138  0.383  1.854 39.944 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 43.57818    3.56234  12.233  < 2e-16 ***
## agean        0.21508    0.06774   3.175  0.00173 ** 
## gndrFemale  -2.65578    0.58485  -4.541 9.64e-06 ***
## eduyrsn      0.32387    0.07741   4.184 4.29e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.112 on 201 degrees of freedom
## Multiple R-squared:  0.2176, Adjusted R-squared:  0.206 
## F-statistic: 18.64 on 3 and 201 DF,  p-value: 1.044e-10

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -12.237  -2.033   0.012   1.716  39.567 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 46.25992    2.20516  20.978  < 2e-16 ***
## agean        0.20594    0.03765   5.470 8.67e-08 ***
## gndrFemale  -2.44047    0.46178  -5.285 2.24e-07 ***
## eduyrsn      0.19266    0.08118   2.373   0.0182 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.295 on 344 degrees of freedom
## Multiple R-squared:  0.1618, Adjusted R-squared:  0.1545 
## F-statistic: 22.13 on 3 and 344 DF,  p-value: 3.982e-13

## 
## Call:
## lm(formula = agertrn ~ agean + gndr + eduyrsn, data = datasub)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -9.4062 -0.4445  0.7586  1.2636 11.6571 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 54.0154859  2.0780704  25.993   <2e-16 ***
## agean        0.0798759  0.0345060   2.315   0.0215 *  
## gndrFemale  -3.9358786  0.4045747  -9.728   <2e-16 ***
## eduyrsn      0.0008278  0.0758679   0.011   0.9913    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.979 on 233 degrees of freedom
## Multiple R-squared:  0.3047, Adjusted R-squared:  0.2957 
## F-statistic: 34.03 on 3 and 233 DF,  p-value: < 2.2e-16