R markdown_homework 2
2025-05-12
Part One - Likert Scale Frequency Table and Plot for Depression
Basic frequencies
## Item None or almost none of the time Some of the time Most of the time
## 1 fltdpr 64.915835 29.06631 4.557165
## 2 flteeff 48.395568 38.42383 9.814171
## 3 slprl 43.873854 39.87056 11.625059
## 4 wrhpp 4.003510 23.53973 48.886939
## 5 fltlnl 68.136458 24.27532 5.302253
## 6 enjlf 5.338783 24.82572 44.804153
## 7 fltsd 52.489933 41.07451 4.859808
## 8 cldgng 55.673484 36.10353 6.217928
## All or almost all of the time
## 1 1.460694
## 2 3.366431
## 3 4.630532
## 4 23.569817
## 5 2.285972
## 6 25.031346
## 7 1.575748
## 8 2.005056Basic plot 1
Append means and counts to table
## Item None or almost none of the time Some of the time Most of the time
## 1 fltdpr 64.915835 29.06631 4.557165
## 2 flteeff 48.395568 38.42383 9.814171
## 3 slprl 43.873854 39.87056 11.625059
## 4 wrhpp 4.003510 23.53973 48.886939
## 5 fltlnl 68.136458 24.27532 5.302253
## 6 enjlf 5.338783 24.82572 44.804153
## 7 fltsd 52.489933 41.07451 4.859808
## 8 cldgng 55.673484 36.10353 6.217928
## All or almost all of the time Mean Count
## 1 1.460694 1.425627 39981
## 2 3.366431 1.681515 39983
## 3 4.630532 1.770123 40017
## 4 23.569817 2.920231 39890
## 5 2.285972 1.417377 39983
## 6 25.031346 2.895281 39878
## 7 1.575748 1.555214 39981
## 8 2.005056 1.545546 39949Formatted table
| Item | None or almost none of the time | Some of the time | Most of the time | All or almost all of the time | Mean | Count |
|---|---|---|---|---|---|---|
| How much of the time during the past week you felt depressed? | 64.9 | 29.1 | 4.6 | 1.5 | 1.4 | 39981 |
| How much of the time during the past week you felt that everything you did was an effort? | 48.4 | 38.4 | 9.8 | 3.4 | 1.7 | 39983 |
| How much of the time during the past week your sleep was restless? | 43.9 | 39.9 | 11.6 | 4.6 | 1.8 | 40017 |
| How much of the time during the past week you were happy? | 4.0 | 23.5 | 48.9 | 23.6 | 2.9 | 39890 |
| How much of the time during the past week you felt lonely? | 68.1 | 24.3 | 5.3 | 2.3 | 1.4 | 39983 |
| How much of the time during the past week you enjoyed life? | 5.3 | 24.8 | 44.8 | 25.0 | 2.9 | 39878 |
| How much of the time during the past week you felt sad? | 52.5 | 41.1 | 4.9 | 1.6 | 1.6 | 39981 |
| How much of the time during the past week you could not get going?. | 55.7 | 36.1 | 6.2 | 2.0 | 1.5 | 39949 |
Basic plot 2
Part Two - Non-weighted vs weighted model
Non-weighted Model
##
## Call:
## lm(formula = depression ~ alcfreq + gndr + cgtsmok + gndr + agea_n +
## eisced + health, data = df_de)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.04472 -0.27157 -0.05028 0.20790 2.16524
##
## Coefficients:
## Estimate Std. Error
## (Intercept) 1.5204819 0.0786949
## alcfreqSeveral times a week -0.0221115 0.0399626
## alcfreqOnce a week -0.0123071 0.0406397
## alcfreq2-3 times a month 0.0137760 0.0435646
## alcfreqOnce a month 0.0075431 0.0464458
## alcfreqLess than once a month 0.0626197 0.0427665
## alcfreqNever 0.1115571 0.0431254
## gndrFemale 0.0806295 0.0177140
## cgtsmokI smoke daily, 9 or fewer cigarettes -0.0257815 0.0431473
## cgtsmokI smoke but not every day 0.0010662 0.0474679
## cgtsmokI don’t smoke now but I used to -0.0882779 0.0299806
## cgtsmokI have only smoked a few times -0.0688746 0.0352542
## cgtsmokI have never smoked -0.1048450 0.0281665
## agea_n -0.0023650 0.0005396
## eiscedES-ISCED II, lower secondary 0.1305194 0.0710685
## eiscedES-ISCED IIIb, lower tier upper secondary 0.0320965 0.0671653
## eiscedES-ISCED IIIa, upper tier upper secondary 0.0605257 0.0752024
## eiscedES-ISCED IV, advanced vocational, sub-degree 0.0121947 0.0678068
## eiscedES-ISCED V1, lower tertiary education, BA level 0.0108585 0.0706733
## eiscedES-ISCED V2, higher tertiary education, >= MA level 0.0051974 0.0690300
## healthGood 0.1644785 0.0242680
## healthFair 0.3517344 0.0270187
## healthBad 0.6969733 0.0366309
## healthVery bad 0.8613992 0.0850057
## t value Pr(>|t|)
## (Intercept) 19.321 < 2e-16 ***
## alcfreqSeveral times a week -0.553 0.580107
## alcfreqOnce a week -0.303 0.762043
## alcfreq2-3 times a month 0.316 0.751862
## alcfreqOnce a month 0.162 0.871000
## alcfreqLess than once a month 1.464 0.143267
## alcfreqNever 2.587 0.009746 **
## gndrFemale 4.552 5.59e-06 ***
## cgtsmokI smoke daily, 9 or fewer cigarettes -0.598 0.550216
## cgtsmokI smoke but not every day 0.022 0.982082
## cgtsmokI don’t smoke now but I used to -2.944 0.003266 **
## cgtsmokI have only smoked a few times -1.954 0.050861 .
## cgtsmokI have never smoked -3.722 0.000202 ***
## agea_n -4.382 1.23e-05 ***
## eiscedES-ISCED II, lower secondary 1.837 0.066406 .
## eiscedES-ISCED IIIb, lower tier upper secondary 0.478 0.632785
## eiscedES-ISCED IIIa, upper tier upper secondary 0.805 0.420995
## eiscedES-ISCED IV, advanced vocational, sub-degree 0.180 0.857290
## eiscedES-ISCED V1, lower tertiary education, BA level 0.154 0.877904
## eiscedES-ISCED V2, higher tertiary education, >= MA level 0.075 0.939989
## healthGood 6.778 1.54e-11 ***
## healthFair 13.018 < 2e-16 ***
## healthBad 19.027 < 2e-16 ***
## healthVery bad 10.133 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4132 on 2354 degrees of freedom
## (42 observations deleted due to missingness)
## Multiple R-squared: 0.2254, Adjusted R-squared: 0.2178
## F-statistic: 29.78 on 23 and 2354 DF, p-value: < 2.2e-16Weighted Model
##
## Call:
## lm(formula = depression ~ alcfreq + gndr + cgtsmok + gndr + agea_n +
## eisced + health, data = df_de, weights = pspwght)
##
## Weighted Residuals:
## Min 1Q Median 3Q Max
## -1.5726 -0.2070 -0.0391 0.1498 3.7044
##
## Coefficients:
## Estimate Std. Error
## (Intercept) 1.4930473 0.0635045
## alcfreqSeveral times a week -0.0183075 0.0411650
## alcfreqOnce a week -0.0039897 0.0413555
## alcfreq2-3 times a month -0.0385177 0.0439932
## alcfreqOnce a month 0.0058139 0.0470853
## alcfreqLess than once a month 0.0396259 0.0423457
## alcfreqNever 0.0988902 0.0426726
## gndrFemale 0.0889036 0.0180549
## cgtsmokI smoke daily, 9 or fewer cigarettes 0.0215546 0.0409516
## cgtsmokI smoke but not every day -0.0101073 0.0472732
## cgtsmokI don’t smoke now but I used to -0.0662052 0.0290395
## cgtsmokI have only smoked a few times -0.0891170 0.0342470
## cgtsmokI have never smoked -0.0966172 0.0266338
## agea_n -0.0028069 0.0005121
## eiscedES-ISCED II, lower secondary 0.1422677 0.0488486
## eiscedES-ISCED IIIb, lower tier upper secondary 0.0438513 0.0477990
## eiscedES-ISCED IIIa, upper tier upper secondary 0.1849768 0.0611179
## eiscedES-ISCED IV, advanced vocational, sub-degree 0.0124237 0.0496833
## eiscedES-ISCED V1, lower tertiary education, BA level 0.1033921 0.0563422
## eiscedES-ISCED V2, higher tertiary education, >= MA level 0.0010199 0.0523476
## healthGood 0.1864119 0.0260783
## healthFair 0.4087645 0.0288664
## healthBad 0.7640432 0.0367571
## healthVery bad 0.9240594 0.0832794
## t value Pr(>|t|)
## (Intercept) 23.511 < 2e-16 ***
## alcfreqSeveral times a week -0.445 0.656552
## alcfreqOnce a week -0.096 0.923154
## alcfreq2-3 times a month -0.876 0.381372
## alcfreqOnce a month 0.123 0.901741
## alcfreqLess than once a month 0.936 0.349487
## alcfreqNever 2.317 0.020567 *
## gndrFemale 4.924 9.06e-07 ***
## cgtsmokI smoke daily, 9 or fewer cigarettes 0.526 0.598700
## cgtsmokI smoke but not every day -0.214 0.830717
## cgtsmokI don’t smoke now but I used to -2.280 0.022707 *
## cgtsmokI have only smoked a few times -2.602 0.009321 **
## cgtsmokI have never smoked -3.628 0.000292 ***
## agea_n -5.481 4.67e-08 ***
## eiscedES-ISCED II, lower secondary 2.912 0.003620 **
## eiscedES-ISCED IIIb, lower tier upper secondary 0.917 0.359021
## eiscedES-ISCED IIIa, upper tier upper secondary 3.027 0.002500 **
## eiscedES-ISCED IV, advanced vocational, sub-degree 0.250 0.802565
## eiscedES-ISCED V1, lower tertiary education, BA level 1.835 0.066621 .
## eiscedES-ISCED V2, higher tertiary education, >= MA level 0.019 0.984457
## healthGood 7.148 1.17e-12 ***
## healthFair 14.161 < 2e-16 ***
## healthBad 20.786 < 2e-16 ***
## healthVery bad 11.096 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4153 on 2354 degrees of freedom
## (42 observations deleted due to missingness)
## Multiple R-squared: 0.2653, Adjusted R-squared: 0.2582
## F-statistic: 36.96 on 23 and 2354 DF, p-value: < 2.2e-16