## Warning: NAs introduced by coercion
## Warning: NAs introduced by coercion

Analysis 1: Compare BMI among different social civil status after controlling for age and gender

## 
## Call:
## lm(formula = BMI ~ soc_civilstatus + VisitAge + Gender, data = data)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -10.7186  -2.7291  -0.5832   2.0990  27.9581 
## 
## Coefficients:
##                   Estimate Std. Error t value Pr(>|t|)    
## (Intercept)      24.609104   0.208235 118.179  < 2e-16 ***
## soc_civilstatus2  0.075695   0.082142   0.922  0.35679    
## soc_civilstatus3  0.426589   0.106593   4.002 6.30e-05 ***
## soc_civilstatus4  0.051651   0.096590   0.535  0.59283    
## soc_civilstatus5 -0.379644   0.134353  -2.826  0.00472 ** 
## soc_civilstatus6  0.133699   0.159130   0.840  0.40081    
## VisitAge          0.019125   0.003368   5.679 1.38e-08 ***
## GenderMale        1.065262   0.057387  18.563  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.952 on 20032 degrees of freedom
##   (503 observations deleted due to missingness)
## Multiple R-squared:  0.02093,    Adjusted R-squared:  0.02059 
## F-statistic: 61.17 on 7 and 20032 DF,  p-value: < 2.2e-16
## 
## Call:
## lm(formula = BMI ~ soc_civilstatus + VisitAge + Gender, data = data)
## 
## Standardized Coefficients::
##      (Intercept) soc_civilstatus2 soc_civilstatus3 soc_civilstatus4 
##      0.000000000      0.006684679      0.028666543      0.003866054 
## soc_civilstatus5 soc_civilstatus6         VisitAge       GenderMale 
##     -0.020062878      0.006036648      0.040725357      0.132197187
##                    Df Sum Sq Mean Sq F value   Pr(>F)    
## soc_civilstatus     5    426      85   5.454 5.10e-05 ***
## VisitAge            1    880     880  56.362 6.28e-14 ***
## Gender              1   5381    5381 344.576  < 2e-16 ***
## Residuals       20032 312812      16                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 503 observations deleted due to missingness

Analysis 2: Compare sle_needhrs_rested (Executed by number of hours of sleep per day) among different social civil status after controlling for age and gender. 4 hours or less [1] … 10 hours or more [7]

## 
## Call:
## lm(formula = sle_needhrs_rested ~ soc_civilstatus + VisitAge + 
##     Gender, data = data)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.6717 -0.5176 -0.0766  0.5730  2.9234 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)       5.3012814  0.0471916 112.335   <2e-16 ***
## soc_civilstatus2 -0.0101374  0.0186205  -0.544    0.586    
## soc_civilstatus3 -0.0066102  0.0242368  -0.273    0.785    
## soc_civilstatus4  0.0027383  0.0219323   0.125    0.901    
## soc_civilstatus5  0.0121980  0.0302610   0.403    0.687    
## soc_civilstatus6 -0.0042469  0.0361905  -0.117    0.907    
## VisitAge         -0.0123443  0.0007619 -16.203   <2e-16 ***
## GenderMale       -0.2988539  0.0130099 -22.971   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.892 on 19840 degrees of freedom
##   (695 observations deleted due to missingness)
## Multiple R-squared:  0.04288,    Adjusted R-squared:  0.04254 
## F-statistic:   127 on 7 and 19840 DF,  p-value: < 2.2e-16
## 
## Call:
## lm(formula = sle_needhrs_rested ~ soc_civilstatus + VisitAge + 
##     Gender, data = data)
## 
## Standardized Coefficients::
##      (Intercept) soc_civilstatus2 soc_civilstatus3 soc_civilstatus4 
##     0.0000000000    -0.0039234681    -0.0019404845     0.0008967544 
## soc_civilstatus5 soc_civilstatus6         VisitAge       GenderMale 
##     0.0028436614    -0.0008371297    -0.1153925252    -0.1625176747
##                    Df Sum Sq Mean Sq F value  Pr(>F)    
## soc_civilstatus     5     14     2.8   3.544 0.00333 ** 
## VisitAge            1    273   273.3 343.456 < 2e-16 ***
## Gender              1    420   419.9 527.682 < 2e-16 ***
## Residuals       19840  15786     0.8                    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 695 observations deleted due to missingness

Analysis 3: Compare sle_sleephrs (Number of hours of sleep per day.) among different social civil status after controlling for age and gender. 4 hours or less [1] … 8 hours [5]

## 
## Call:
## lm(formula = sle_sleephrs ~ soc_civilstatus + VisitAge + Gender, 
##     data = data)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.0663 -0.8282  0.1011  0.9718  3.3382 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)       3.3424509  0.0531225  62.920  < 2e-16 ***
## soc_civilstatus2 -0.0064461  0.0209913  -0.307 0.758783    
## soc_civilstatus3 -0.1011991  0.0271642  -3.725 0.000195 ***
## soc_civilstatus4 -0.1288213  0.0246009  -5.236 1.65e-07 ***
## soc_civilstatus5 -0.0471216  0.0341955  -1.378 0.168216    
## soc_civilstatus6 -0.1483138  0.0403149  -3.679 0.000235 ***
## VisitAge          0.0095240  0.0008575  11.107  < 2e-16 ***
## GenderMale       -0.0270614  0.0146298  -1.850 0.064365 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.018 on 20466 degrees of freedom
##   (69 observations deleted due to missingness)
## Multiple R-squared:  0.00801,    Adjusted R-squared:  0.007671 
## F-statistic: 23.61 on 7 and 20466 DF,  p-value: < 2.2e-16
## 
## Call:
## lm(formula = sle_sleephrs ~ soc_civilstatus + VisitAge + Gender, 
##     data = data)
## 
## Standardized Coefficients::
##      (Intercept) soc_civilstatus2 soc_civilstatus3 soc_civilstatus4 
##      0.000000000     -0.002217868     -0.026574879     -0.037706129 
## soc_civilstatus5 soc_civilstatus6         VisitAge       GenderMale 
##     -0.009743334     -0.026324312      0.079309704     -0.013119225
##                    Df Sum Sq Mean Sq F value   Pr(>F)    
## soc_civilstatus     5     43    8.57   8.258 8.35e-08 ***
## VisitAge            1    125  125.04 120.544  < 2e-16 ***
## Gender              1      4    3.55   3.422   0.0644 .  
## Residuals       20466  21230    1.04                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 69 observations deleted due to missingness

Analysis 4: Compare sle_general_gen (How do you usually sleep?) among different social civil status after controlling for age and gender. very bad [1] … very good [5]

## 
## Call:
## lm(formula = sle_general_gen ~ soc_civilstatus + VisitAge + Gender, 
##     data = data)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -2.966 -0.708  0.062  1.042 67.037 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)       3.6438864  0.0590110  61.749  < 2e-16 ***
## soc_civilstatus2  0.0032260  0.0233180   0.138   0.8900    
## soc_civilstatus3 -0.1429499  0.0301545  -4.741 2.15e-06 ***
## soc_civilstatus4 -0.1246560  0.0273080  -4.565 5.03e-06 ***
## soc_civilstatus5 -0.0385761  0.0379512  -1.016   0.3094    
## soc_civilstatus6 -0.0762997  0.0447651  -1.704   0.0883 .  
## VisitAge          0.0009745  0.0009525   1.023   0.3063    
## GenderMale        0.2484590  0.0162519  15.288  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.133 on 20515 degrees of freedom
##   (20 observations deleted due to missingness)
## Multiple R-squared:  0.01522,    Adjusted R-squared:  0.01489 
## F-statistic: 45.31 on 7 and 20515 DF,  p-value: < 2.2e-16
## 
## Call:
## lm(formula = sle_general_gen ~ soc_civilstatus + VisitAge + Gender, 
##     data = data)
## 
## Standardized Coefficients::
##      (Intercept) soc_civilstatus2 soc_civilstatus3 soc_civilstatus4 
##     0.0000000000     0.0009944946    -0.0336546137    -0.0327138166 
## soc_civilstatus5 soc_civilstatus6         VisitAge       GenderMale 
##    -0.0071526703    -0.0121365453     0.0072699860     0.1079035511
##                    Df Sum Sq Mean Sq F value   Pr(>F)    
## soc_civilstatus     5     99   19.78  15.421 3.62e-15 ***
## VisitAge            1      8    8.10   6.317    0.012 *  
## Gender              1    300  299.81 233.723  < 2e-16 ***
## Residuals       20515  26316    1.28                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 20 observations deleted due to missingness

Analysis 5: Compare alc_drink_frq among different social civil status after controlling for age and gender.  4 times a week or more often [1]…  Never [6]

## 
## Call:
## lm(formula = alc_drink_frq ~ soc_civilstatus + VisitAge + Gender, 
##     data = data)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.9416 -1.0424 -0.0872  1.0140  3.2934 
## 
## Coefficients:
##                   Estimate Std. Error t value Pr(>|t|)    
## (Intercept)       3.791183   0.070450  53.814  < 2e-16 ***
## soc_civilstatus2  0.011284   0.027695   0.407    0.684    
## soc_civilstatus3  0.675669   0.036037  18.749  < 2e-16 ***
## soc_civilstatus4  0.413629   0.032589  12.692  < 2e-16 ***
## soc_civilstatus5 -0.012872   0.045023  -0.286    0.775    
## soc_civilstatus6  0.431926   0.053723   8.040 9.49e-16 ***
## VisitAge         -0.011176   0.001137  -9.830  < 2e-16 ***
## GenderMale       -0.235229   0.019367 -12.146  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.337 on 20106 degrees of freedom
##   (429 observations deleted due to missingness)
## Multiple R-squared:  0.04142,    Adjusted R-squared:  0.04108 
## F-statistic: 124.1 on 7 and 20106 DF,  p-value: < 2.2e-16
## 
## Call:
## lm(formula = alc_drink_frq ~ soc_civilstatus + VisitAge + Gender, 
##     data = data)
## 
## Standardized Coefficients::
##      (Intercept) soc_civilstatus2 soc_civilstatus3 soc_civilstatus4 
##      0.000000000      0.002920223      0.132678748      0.090663353 
## soc_civilstatus5 soc_civilstatus6         VisitAge       GenderMale 
##     -0.002005548      0.057030165     -0.069602508     -0.085448033
##                    Df Sum Sq Mean Sq F value Pr(>F)    
## soc_civilstatus     5   1070  214.07   119.8 <2e-16 ***
## VisitAge            1    218  218.15   122.1 <2e-16 ***
## Gender              1    264  263.54   147.5 <2e-16 ***
## Residuals       20106  35920    1.79                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 429 observations deleted due to missingness

Analysis 6: Compare alc_drink_summ among different social civil status after controlling for age and gender. 1 to 2 glasses [1] …10 glasses or more [5]

## 
## Call:
## lm(formula = alc_drink_summ ~ soc_civilstatus + VisitAge + Gender, 
##     data = data)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -0.8552 -0.4057 -0.2040  0.4308  3.6147 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)       1.8168924  0.0322414  56.353  < 2e-16 ***
## soc_civilstatus2  0.0352048  0.0125971   2.795   0.0052 ** 
## soc_civilstatus3  0.1265463  0.0168355   7.517 5.87e-14 ***
## soc_civilstatus4  0.0832261  0.0149707   5.559 2.74e-08 ***
## soc_civilstatus5  0.0919403  0.0205528   4.473 7.74e-06 ***
## soc_civilstatus6  0.0170238  0.0247574   0.688   0.4917    
## VisitAge         -0.0102155  0.0005209 -19.611  < 2e-16 ***
## GenderMale        0.3345495  0.0088669  37.730  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.6001 on 19337 degrees of freedom
##   (1198 observations deleted due to missingness)
## Multiple R-squared:  0.08425,    Adjusted R-squared:  0.08392 
## F-statistic: 254.1 on 7 and 19337 DF,  p-value: < 2.2e-16
## 
## Call:
## lm(formula = alc_drink_summ ~ soc_civilstatus + VisitAge + Gender, 
##     data = data)
## 
## Standardized Coefficients::
##      (Intercept) soc_civilstatus2 soc_civilstatus3 soc_civilstatus4 
##      0.000000000      0.019956630      0.052979985      0.039552527 
## soc_civilstatus5 soc_civilstatus6         VisitAge       GenderMale 
##      0.031274661      0.004860106     -0.138466982      0.264593990
##                    Df Sum Sq Mean Sq F value Pr(>F)    
## soc_civilstatus     5     37     7.3   20.38 <2e-16 ***
## VisitAge            1     91    91.3  253.60 <2e-16 ***
## Gender              1    513   512.6 1423.55 <2e-16 ***
## Residuals       19337   6963     0.4                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 1198 observations deleted due to missingness

Analysis 7: Compare scr_sc23 (Lost interest> 2 weeks.) among different social civil status after controlling for age and gender. yes [1] no [0]

## 
## Call:
## glm(formula = scr_sc23 ~ soc_civilstatus + VisitAge + Gender, 
##     family = binomial, data = data)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -1.4873  -0.9566  -0.7485   1.2381   1.8634  
## 
## Coefficients:
##                   Estimate Std. Error z value Pr(>|z|)    
## (Intercept)       1.226357   0.110890  11.059  < 2e-16 ***
## soc_civilstatus2  0.318870   0.043344   7.357 1.89e-13 ***
## soc_civilstatus3  0.688281   0.055034  12.506  < 2e-16 ***
## soc_civilstatus4  0.727522   0.049784  14.613  < 2e-16 ***
## soc_civilstatus5  0.610356   0.069012   8.844  < 2e-16 ***
## soc_civilstatus6  0.835612   0.081184  10.293  < 2e-16 ***
## VisitAge         -0.029513   0.001807 -16.335  < 2e-16 ***
## GenderMale       -0.496119   0.031005 -16.001  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 26696  on 20231  degrees of freedom
## Residual deviance: 25580  on 20224  degrees of freedom
##   (311 observations deleted due to missingness)
## AIC: 25596
## 
## Number of Fisher Scoring iterations: 4
## 
## Call:
## glm(formula = scr_sc23 ~ soc_civilstatus + VisitAge + Gender, 
##     family = binomial, data = data)
## 
## Standardized Coefficients::
##      (Intercept) soc_civilstatus2 soc_civilstatus3 soc_civilstatus4 
##        0.0000000        0.2323992        0.3810884        0.4498753 
## soc_civilstatus5 soc_civilstatus6         VisitAge       GenderMale 
##        0.2673434        0.3132077       -0.5195312       -0.5088141
##                    Df Sum Sq Mean Sq F value Pr(>F)    
## soc_civilstatus     5    127   25.31   114.6 <2e-16 ***
## VisitAge            1     73   72.75   329.3 <2e-16 ***
## Gender              1     57   57.04   258.2 <2e-16 ***
## Residuals       20224   4468    0.22                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 311 observations deleted due to missingness

Analysis 8 : Compare sle_sleep_qual among different social civil status after controlling for age and gender. never or rarely[4] … 4 times or more a week [1]

## [1] "sle_sleep_qual_brea"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.9000  0.1396  0.2320  0.4101  2.1774 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.1700129  0.0438235  95.155  < 2e-16 ***
## data$soc_civilstatus2 -0.0095807  0.0171870  -0.557 0.577234    
## data$soc_civilstatus3 -0.0023276  0.0239823  -0.097 0.922683    
## data$soc_civilstatus4 -0.0784458  0.0213699  -3.671 0.000242 ***
## data$soc_civilstatus5  0.0006517  0.0280532   0.023 0.981466    
## data$soc_civilstatus6 -0.0549738  0.0358657  -1.533 0.125351    
## data$VisitAge         -0.0060009  0.0007101  -8.451  < 2e-16 ***
## data$GenderMale       -0.2440156  0.0121200 -20.133  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7878 on 17765 degrees of freedom
##   (2770 observations deleted due to missingness)
## Multiple R-squared:  0.02888,    Adjusted R-squared:  0.02849 
## F-statistic: 75.47 on 7 and 17765 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value Pr(>F)    
## data$soc_civilstatus     5      5    0.95   1.538  0.174    
## data$VisitAge            1     72   71.51 115.219 <2e-16 ***
## data$Gender              1    252  251.58 405.351 <2e-16 ***
## Residuals            17765  11026    0.62                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 2770 observations deleted due to missingness
## [1] "sle_sleep_qual_daya"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.6863 -0.5019  0.3909  0.4749  0.6570 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            3.8599792  0.0399154  96.704  < 2e-16 ***
## data$soc_civilstatus2 -0.0044647  0.0157615  -0.283  0.77697    
## data$soc_civilstatus3 -0.1063028  0.0204402  -5.201 2.00e-07 ***
## data$soc_civilstatus4 -0.1226080  0.0185064  -6.625 3.55e-11 ***
## data$soc_civilstatus5 -0.0783905  0.0256681  -3.054  0.00226 ** 
## data$soc_civilstatus6 -0.0517674  0.0303053  -1.708  0.08762 .  
## data$VisitAge         -0.0038598  0.0006445  -5.989 2.15e-09 ***
## data$GenderMale       -0.1010643  0.0109885  -9.197  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7641 on 20410 degrees of freedom
##   (125 observations deleted due to missingness)
## Multiple R-squared:  0.008917,   Adjusted R-squared:  0.008577 
## F-statistic: 26.23 on 7 and 20410 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5     30    5.98   10.23 8.18e-10 ***
## data$VisitAge            1     28   27.95   47.87 4.69e-12 ***
## data$Gender              1     49   49.39   84.59  < 2e-16 ***
## Residuals            20410  11916    0.58                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 125 observations deleted due to missingness
## [1] "sle_sleep_qual_dayt"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4411 -0.3745  0.5746  0.7476  1.0546 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            2.9981729  0.0487137  61.547  < 2e-16 ***
## data$soc_civilstatus2 -0.0431137  0.0192236  -2.243 0.024924 *  
## data$soc_civilstatus3 -0.2287600  0.0249847  -9.156  < 2e-16 ***
## data$soc_civilstatus4 -0.1727568  0.0225572  -7.659 1.96e-14 ***
## data$soc_civilstatus5 -0.1204336  0.0313350  -3.843 0.000122 ***
## data$soc_civilstatus6 -0.0969685  0.0369598  -2.624 0.008707 ** 
## data$VisitAge          0.0039119  0.0007865   4.974 6.61e-07 ***
## data$GenderMale        0.1455748  0.0134154  10.851  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9318 on 20367 degrees of freedom
##   (168 observations deleted due to missingness)
## Multiple R-squared:  0.01603,    Adjusted R-squared:  0.01569 
## F-statistic: 47.39 on 7 and 20367 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    154   30.81   35.49  < 2e-16 ***
## data$VisitAge            1     32   31.71   36.53 1.53e-09 ***
## data$Gender              1    102  102.23  117.75  < 2e-16 ***
## Residuals            20367  17682    0.87                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 168 observations deleted due to missingness
## [1] "sle_sleep_qual_earl"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.06016 -0.94657 -0.03466  0.95609  1.11531 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            3.0772253  0.0534032  57.622  < 2e-16 ***
## data$soc_civilstatus2  0.0017616  0.0210896   0.084  0.93343    
## data$soc_civilstatus3 -0.1456502  0.0273373  -5.328 1.00e-07 ***
## data$soc_civilstatus4 -0.1447739  0.0247749  -5.844 5.19e-09 ***
## data$soc_civilstatus5 -0.0984183  0.0343495  -2.865  0.00417 ** 
## data$soc_civilstatus6 -0.0840485  0.0407628  -2.062  0.03923 *  
## data$VisitAge         -0.0006169  0.0008627  -0.715  0.47457    
## data$GenderMale        0.0089360  0.0147141   0.607  0.54366    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.019 on 20263 degrees of freedom
##   (272 observations deleted due to missingness)
## Multiple R-squared:  0.003326,   Adjusted R-squared:  0.002982 
## F-statistic:  9.66 on 7 and 20263 DF,  p-value: 4.688e-12
## 
##                         Df Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5     69  13.883  13.363 4.95e-13 ***
## data$VisitAge            1      0   0.451   0.434    0.510    
## data$Gender              1      0   0.383   0.369    0.544    
## Residuals            20263  21051   1.039                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 272 observations deleted due to missingness
## [1] "sle_sleep_qual_hard"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.4870 -0.4148  0.0510  0.7622  1.2570 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            3.5675606  0.0482122  73.997  < 2e-16 ***
## data$soc_civilstatus2  0.0055868  0.0190423   0.293  0.76923    
## data$soc_civilstatus3 -0.1360880  0.0246728  -5.516 3.52e-08 ***
## data$soc_civilstatus4 -0.1387014  0.0223521  -6.205 5.56e-10 ***
## data$soc_civilstatus5 -0.0680497  0.0310075  -2.195  0.02820 *  
## data$soc_civilstatus6 -0.1150509  0.0366451  -3.140  0.00169 ** 
## data$VisitAge         -0.0090247  0.0007783 -11.595  < 2e-16 ***
## data$GenderMale        0.3199787  0.0132786  24.097  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9235 on 20428 degrees of freedom
##   (107 observations deleted due to missingness)
## Multiple R-squared:  0.03946,    Adjusted R-squared:  0.03913 
## F-statistic: 119.9 on 7 and 20428 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value Pr(>F)    
## data$soc_civilstatus     5    147    29.3   34.41 <2e-16 ***
## data$VisitAge            1     74    73.8   86.54 <2e-16 ***
## data$Gender              1    495   495.2  580.68 <2e-16 ***
## Residuals            20428  17421     0.9                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 107 observations deleted due to missingness
## [1] "sle_sleep_qual_legs"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.7501  0.2716  0.3130  0.4278  0.5065 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            3.6630591  0.0408212  89.734   <2e-16 ***
## data$soc_civilstatus2 -0.0281997  0.0161320  -1.748   0.0805 .  
## data$soc_civilstatus3  0.0231686  0.0208703   1.110   0.2670    
## data$soc_civilstatus4 -0.0581026  0.0189710  -3.063   0.0022 ** 
## data$soc_civilstatus5 -0.0384732  0.0263079  -1.462   0.1436    
## data$soc_civilstatus6  0.0023177  0.0310747   0.075   0.9405    
## data$VisitAge         -0.0014661  0.0006589  -2.225   0.0261 *  
## data$GenderMale        0.1312888  0.0112514  11.669   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7791 on 20262 degrees of freedom
##   (273 observations deleted due to missingness)
## Multiple R-squared:  0.008064,   Adjusted R-squared:  0.007721 
## F-statistic: 23.53 on 7 and 20262 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5     17    3.32   5.475 4.86e-05 ***
## data$VisitAge            1      1    0.72   1.191    0.275    
## data$Gender              1     83   82.66 136.158  < 2e-16 ***
## Residuals            20262  12300    0.61                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 273 observations deleted due to missingness
## [1] "sle_sleep_qual_notr"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.3135 -0.7397  0.1369  0.9099  1.7377 
## 
## Coefficients:
##                        Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            1.262924   0.052714  23.958  < 2e-16 ***
## data$soc_civilstatus2 -0.028702   0.020816  -1.379  0.16795    
## data$soc_civilstatus3 -0.118842   0.026971  -4.406 1.06e-05 ***
## data$soc_civilstatus4 -0.151991   0.024431  -6.221 5.03e-10 ***
## data$soc_civilstatus5 -0.075932   0.033998  -2.233  0.02553 *  
## data$soc_civilstatus6 -0.115439   0.040354  -2.861  0.00423 ** 
## data$VisitAge          0.025031   0.000852  29.380  < 2e-16 ***
## data$GenderMale        0.148300   0.014529  10.207  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.004 on 20168 degrees of freedom
##   (367 observations deleted due to missingness)
## Multiple R-squared:  0.05383,    Adjusted R-squared:  0.0535 
## F-statistic: 163.9 on 7 and 20168 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value Pr(>F)    
## data$soc_civilstatus     5    112    22.4   22.23 <2e-16 ***
## data$VisitAge            1    940   939.9  932.00 <2e-16 ***
## data$Gender              1    105   105.1  104.19 <2e-16 ***
## Residuals            20168  20339     1.0                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 367 observations deleted due to missingness
## [1] "sle_sleep_qual_seve"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -2.35001 -0.79067  0.07215  0.85758  1.26065 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            3.4582894  0.0494827  69.889   <2e-16 ***
## data$soc_civilstatus2  0.0313347  0.0195549   1.602    0.109    
## data$soc_civilstatus3  0.0155356  0.0252933   0.614    0.539    
## data$soc_civilstatus4  0.0174299  0.0229243   0.760    0.447    
## data$soc_civilstatus5  0.0040851  0.0318342   0.128    0.898    
## data$soc_civilstatus6 -0.0238793  0.0376895  -0.634    0.526    
## data$VisitAge         -0.0091455  0.0007988 -11.449   <2e-16 ***
## data$GenderMale        0.2877304  0.0136256  21.117   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9472 on 20403 degrees of freedom
##   (132 observations deleted due to missingness)
## Multiple R-squared:  0.02708,    Adjusted R-squared:  0.02674 
## F-statistic: 81.12 on 7 and 20403 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5     29     5.8   6.511 4.67e-06 ***
## data$VisitAge            1     80    80.1  89.331  < 2e-16 ***
## data$Gender              1    400   400.0 445.922  < 2e-16 ***
## Residuals            20403  18304     0.9                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 132 observations deleted due to missingness
## [1] "sle_sleep_qual_shor"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.3154 -0.6427  0.1062  0.7410  1.8689 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            1.0170207  0.0491508  20.692  < 2e-16 ***
## data$soc_civilstatus2  0.0011723  0.0193445   0.061 0.951679    
## data$soc_civilstatus3 -0.0703476  0.0252795  -2.783 0.005394 ** 
## data$soc_civilstatus4 -0.1472450  0.0227938  -6.460 1.07e-10 ***
## data$soc_civilstatus5 -0.1033559  0.0315377  -3.277 0.001050 ** 
## data$soc_civilstatus6 -0.1243226  0.0375051  -3.315 0.000919 ***
## data$VisitAge          0.0280284  0.0007947  35.270  < 2e-16 ***
## data$GenderMale        0.1670429  0.0135384  12.338  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.9316 on 19984 degrees of freedom
##   (551 observations deleted due to missingness)
## Multiple R-squared:  0.07468,    Adjusted R-squared:  0.07435 
## F-statistic: 230.4 on 7 and 19984 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value Pr(>F)    
## data$soc_civilstatus     5     99    19.9   22.91 <2e-16 ***
## data$VisitAge            1   1168  1168.1 1345.97 <2e-16 ***
## data$Gender              1    132   132.1  152.24 <2e-16 ***
## Residuals            19984  17343     0.9                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 551 observations deleted due to missingness
## [1] "sle_sleep_qual_snor"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.2946 -0.9518  0.2491  1.0017  1.5748 
## 
## Coefficients:
##                        Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            3.528202   0.061810  57.082  < 2e-16 ***
## data$soc_civilstatus2  0.016421   0.023974   0.685 0.493393    
## data$soc_civilstatus3  0.194049   0.034207   5.673 1.42e-08 ***
## data$soc_civilstatus4  0.030868   0.030270   1.020 0.307851    
## data$soc_civilstatus5  0.130472   0.039615   3.293 0.000991 ***
## data$soc_civilstatus6  0.038435   0.050819   0.756 0.449476    
## data$VisitAge         -0.009297   0.001002  -9.281  < 2e-16 ***
## data$GenderMale       -0.396408   0.017027 -23.281  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.141 on 18901 degrees of freedom
##   (1634 observations deleted due to missingness)
## Multiple R-squared:  0.03929,    Adjusted R-squared:  0.03893 
## F-statistic: 110.4 on 7 and 18901 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value Pr(>F)    
## data$soc_civilstatus     5    116    23.1   17.76 <2e-16 ***
## data$VisitAge            1    185   185.2  142.11 <2e-16 ***
## data$Gender              1    706   706.1  541.99 <2e-16 ***
## Residuals            18901  24626     1.3                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 1634 observations deleted due to missingness
## [1] "sle_sleep_qual_worr"
## 
## Call:
## lm(formula = data[, sle_qual_names[i]] ~ data$soc_civilstatus + 
##     data$VisitAge + data$Gender)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.1915 -0.8117  0.1594  0.9126  1.3095 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            2.4706946  0.0537858  45.936  < 2e-16 ***
## data$soc_civilstatus2 -0.0025420  0.0212439  -0.120    0.905    
## data$soc_civilstatus3  0.0013826  0.0276905   0.050    0.960    
## data$soc_civilstatus4 -0.0400920  0.0249721  -1.605    0.108    
## data$soc_civilstatus5 -0.0289146  0.0345822  -0.836    0.403    
## data$soc_civilstatus6  0.0120580  0.0411580   0.293    0.770    
## data$VisitAge          0.0057749  0.0008692   6.644 3.13e-11 ***
## data$GenderMale        0.2814180  0.0148385  18.965  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.024 on 20127 degrees of freedom
##   (408 observations deleted due to missingness)
## Multiple R-squared:  0.02205,    Adjusted R-squared:  0.02171 
## F-statistic: 64.84 on 7 and 20127 DF,  p-value: < 2.2e-16
## 
##                         Df Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5     22     4.4   4.186 0.000838 ***
## data$VisitAge            1     77    76.8  73.252  < 2e-16 ***
## data$Gender              1    377   377.1 359.687  < 2e-16 ***
## Residuals            20127  21102     1.0                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 408 observations deleted due to missingness

Analysis 9 : Compare die_brea_dummy among different social civil status after controlling for age and gender. Note: for Die_meal, there is no such variables in the dataset

## 
## Call:
## glm(formula = die_brea_dummy ~ soc_civilstatus + VisitAge + Gender, 
##     family = binomial, data = data)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -2.8076   0.2796   0.3537   0.4749   0.9057  
## 
## Coefficients:
##                   Estimate Std. Error z value Pr(>|z|)    
## (Intercept)      -1.800515   0.184861  -9.740  < 2e-16 ***
## soc_civilstatus2  0.079224   0.074230   1.067   0.2858    
## soc_civilstatus3 -0.029718   0.093154  -0.319   0.7497    
## soc_civilstatus4 -0.172371   0.083164  -2.073   0.0382 *  
## soc_civilstatus5 -0.105063   0.114679  -0.916   0.3596    
## soc_civilstatus6 -0.194919   0.161993  -1.203   0.2289    
## VisitAge          0.068879   0.003184  21.632  < 2e-16 ***
## GenderMale        0.408265   0.054645   7.471 7.95e-14 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for binomial family taken to be 1)
## 
##     Null deviance: 11877  on 20522  degrees of freedom
## Residual deviance: 11272  on 20515  degrees of freedom
##   (20 observations deleted due to missingness)
## AIC: 11288
## 
## Number of Fisher Scoring iterations: 5
## 
## Call:
## glm(formula = die_brea_dummy ~ soc_civilstatus + VisitAge + Gender, 
##     family = binomial, data = data)
## 
## Standardized Coefficients::
##      (Intercept) soc_civilstatus2 soc_civilstatus3 soc_civilstatus4 
##       0.00000000       0.10025419      -0.02872036      -0.18569251 
## soc_civilstatus5 soc_civilstatus6         VisitAge       GenderMale 
##      -0.07996718      -0.12727370       2.10935451       0.72784079
##                    Df Sum Sq Mean Sq F value   Pr(>F)    
## soc_civilstatus     5    1.4    0.29   3.821  0.00184 ** 
## VisitAge            1   40.9   40.86 544.503  < 2e-16 ***
## Gender              1    4.0    3.98  52.969 3.51e-13 ***
## Residuals       20515 1539.6    0.08                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 20 observations deleted due to missingness

Analysis 10: Compare qua_stress among different social civil status after controlling for age and gender. Never [5] … Very often [1]

## Warning: NAs introduced by coercion
## [1] "qua_stress_1"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_2"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_3"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_4"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_5"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_6"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_7"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_8"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_9"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness
## Warning: NAs introduced by coercion
## [1] "qua_stress_10"
## 
## Call:
## lm(formula = data$qua_stress_1 ~ data$soc_civilstatus + data$VisitAge + 
##     data$Gender)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -18.05  -7.08  -5.91  -4.98 990.47 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            4.28526    4.24881   1.009   0.3132    
## data$soc_civilstatus2  1.08700    1.67901   0.647   0.5174    
## data$soc_civilstatus3  9.52186    2.17091   4.386 1.16e-05 ***
## data$soc_civilstatus4  0.52700    1.96599   0.268   0.7887    
## data$soc_civilstatus5  4.00461    2.73222   1.466   0.1427    
## data$soc_civilstatus6  5.73127    3.22279   1.778   0.0754 .  
## data$VisitAge          0.06895    0.06858   1.005   0.3147    
## data$GenderMale        1.03908    1.17006   0.888   0.3745    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 81.54 on 20514 degrees of freedom
##   (21 observations deleted due to missingness)
## Multiple R-squared:  0.001168,   Adjusted R-squared:  0.0008269 
## F-statistic: 3.426 on 7 and 20514 DF,  p-value: 0.001152
## 
##                         Df    Sum Sq Mean Sq F value   Pr(>F)    
## data$soc_civilstatus     5    146216   29243   4.398 0.000528 ***
## data$VisitAge            1      7994    7994   1.202 0.272856    
## data$Gender              1      5243    5243   0.789 0.374519    
## Residuals            20514 136389466    6649                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 21 observations deleted due to missingness