setwd("~/Google Drive/Research/Engage data")
library(corrplot)
library(plyr)
library(dplyr)
library(psych)
library(Hmisc)
library(rmarkdown)
library(jmv)
library(stats)
library(readr)
qol <- read_csv("CIW_ENS2021.csv")
qol$PC1 <- as.factor(qol$PC1)
qol$REGION <- as.factor(qol$REGION)
qol$SEX <- as.factor(qol$SEX)
qol$INCOME <- as.factor(qol$INCOME)
qol$MAINACT <- as.factor(qol$MAINACT)
qol$SEXORIENT <- as.factor(qol$SEXORIENT)
qol$MARSTAT <- as.factor(qol$MARSTAT)
qol$ETHNIC1 <- as.factor(qol$ETHNIC1)
qol$EDUCAT <- as.factor(qol$EDUCAT)
qol$BORN_CAN <- as.factor(qol$BORN_CAN)
qol$YEARS_CAN <- as.factor(qol$YEARS_CAN)
qol$DISABILITY <- as.factor(qol$DISABILITY)
qol$CITIZEN <- as.factor(qol$CITIZEN)
qol$PARTHAVE <- as.factor(qol$PARTHAVE)
qol$OLDER <- as.factor(qol$OLDER)
qol$NEW <- as.factor(qol$NEW)
qol$YOUNGER <- as.factor(qol$YOUNGER)
qol$SINGLE_PARENT <- as.factor(qol$SINGLE_PARENT)
qol$LOW_INCOME <- as.factor(qol$LOW_INCOME)
# screening predictors by running basic frequencies
table(qol$REGION, useNA="always")
##
## 1 2 3 4 5 6 7 8 9 10 <NA>
## 4409 1942 877 827 807 584 754 702 745 484 0
table(qol$PARRELIG, useNA="always")
##
## 1 2 <NA>
## 2920 8627 584
table(qol$RELATIVS, useNA="always")
##
## <NA>
## 0
table(qol$FRIENDS, useNA="always")
##
## 0 1 2 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13
## 714 1066 2171 1704 1 1547 1 1375 917 238 435 73 819 51 244 13
## 13.5 14 15 16 17 18 19 20 21 22 23 24 25 26 28 29
## 1 32 137 16 7 11 5 195 6 13 2 12 35 3 2 1
## 30 32 33 34 35 36 37 40 41 42 43 44 45 46 50 51
## 29 6 5 2 6 1 2 17 1 1 5 3 1 3 22 2
## 52 53 54 55 56 61 62 63 64 74 76 80 83 84 85 100
## 3 5 2 3 1 1 2 2 2 1 1 1 2 1 2 9
## 102 103 106 110 124 145 150 156 200 216 250 275 300 310 410 420
## 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1
## 512 530 566 699 <NA>
## 3 1 1 1 116
table(qol$NEIGHBRS, useNA="always")
##
## 0 1 2 2.5 3 4 5 6 7 8 9 10 11 12 13 14
## 1280 1246 2392 1 1391 1478 973 957 232 431 76 700 35 244 13 34
## 15 16 17 18 19 20 21 22 23 24 25 26 28 30 32 34
## 96 36 3 10 3 181 2 3 2 5 43 3 3 51 1 1
## 36 40 42 45 48 50 55 60 70 75 80 100 200 300 <NA>
## 1 17 1 2 1 26 1 4 1 2 1 13 4 1 130
table(qol$PHYSHLTH, useNA="always")
##
## 1 2 3 4 5 <NA>
## 464 1871 4491 4279 993 33
table(qol$MNTLHLTH, useNA="always")
##
## 1 2 3 4 5 <NA>
## 279 1200 3948 4991 1674 39
table(qol$HB_EXERCISE, useNA="always")
##
## 1 2 3 4 5 6 7 <NA>
## 527 567 1948 1749 3668 1918 1685 69
table(qol$POLICY, useNA="always")
##
## 1 2 3 4 5 6 7 8 9 <NA>
## 339 283 548 5287 2140 1616 561 1055 203 99
table(qol$WORKPAY, useNA="always")
##
## 1 2 <NA>
## 5670 6342 119
table(qol$EXP_BILL, useNA="always")
##
## 1 2 3 4 5 6 <NA>
## 9221 833 551 430 393 584 119
table(qol$EXP_HOME, useNA="always")
##
## 1 2 3 4 5 6 <NA>
## 9427 285 130 92 80 1981 136
table(qol$EXP_FOOD, useNA="always")
##
## 1 2 3 4 5 6 <NA>
## 9968 361 309 274 454 630 135
table(qol$EXP_NFOOD, useNA="always")
##
## 1 2 3 4 5 6 <NA>
## 9383 395 425 390 814 602 122
table(qol$EXP_FOODBK, useNA="always")
##
## 1 2 3 4 5 6 <NA>
## 10974 95 54 52 157 670 129
table(qol$EXP_TRANS, useNA="always")
##
## 1 2 3 4 5 6 <NA>
## 10337 288 171 138 226 836 135
table(qol$EXP_WANT, useNA="always")
##
## 1 2 3 4 5 6 <NA>
## 7220 860 1026 683 1559 657 126
table(qol$EXP_NEED, useNA="always")
##
## 1 2 3 4 5 6 <NA>
## 9135 661 490 440 693 568 144
table(qol$GENDER, useNA="always")
##
## 1 2 3 <NA>
## 5059 7024 29 19
table(qol$BORN_CAN, useNA="always")
##
## 1 2 <NA>
## 11117 879 135
table(qol$DISABILITY, useNA="always")
##
## 1 2 <NA>
## 3332 8612 187
table(qol$PARTHAVE, useNA="always")
##
## 1 2 <NA>
## 9204 2761 166
descriptives(qol, vars = vars("RELATVS", "FRIENDS", "NEIGHBRS", "PHYSHLTH", "MNTLHLTH"), missing=TRUE)
##
## DESCRIPTIVES
##
## Descriptives
## ───────────────────────────────────────────────────────────────────────
## RELATVS FRIENDS NEIGHBRS PHYSHLTH MNTLHLTH
## ───────────────────────────────────────────────────────────────────────
## N 12005 12015 12001 12098 12092
## Missing 126 116 130 33 39
## Mean 8.872928 5.799376 4.804641 3.286494 3.544244
## Median 4.000000 4.000000 3.000000 3.000000 4.000000
## Minimum 0.000000 0.000000 0.000000 1.000000 1.000000
## Maximum 13187.00 699.0000 300.0000 5.000000 5.000000
## ───────────────────────────────────────────────────────────────────────
descriptives(qol, vars = vars("HB_EXERCISE", "POLICY", "SLEEP", "VACADAYS", "AGE"), missing=TRUE, box=TRUE)
##
## DESCRIPTIVES
##
## Descriptives
## ──────────────────────────────────────────────────────────────────────────
## HB_EXERCISE POLICY SLEEP VACADAYS AGE
## ──────────────────────────────────────────────────────────────────────────
## N 12062 12032 12007 11853 12126
## Missing 69 99 124 278 5
## Mean 4.654618 4.844332 687.7746 19.58006 59.91062
## Median 5.000000 4.000000 7.000000 10.00000 62.00000
## Minimum 1.000000 1.000000 0.000000 0.000000 14.00000
## Maximum 7.000000 9.000000 8000000 365.0000 101.0000
## ──────────────────────────────────────────────────────────────────────────
descriptives(qol, vars = vars("SOCIAL_CLIMATE", "SOC_HELP", "SOC_NEEDS", "OVERALL_SOC", "SOCIAL_ISO", "TIME_ADEQ"), missing=TRUE, box=TRUE)
##
## DESCRIPTIVES
##
## Descriptives
## ────────────────────────────────────────────────────────────────────────────────────────────────
## SOCIAL_CLIMATE SOC_HELP SOC_NEEDS OVERALL_SOC SOCIAL_ISO TIME_ADEQ
## ────────────────────────────────────────────────────────────────────────────────────────────────
## N 11946 11965 11953 11978 11968 10759
## Missing 185 166 178 153 163 1372
## Mean 4.978926 4.915093 4.535960 4.809479 3.123886 7.694291
## Median 5.000000 5.000000 4.500000 4.833333 3.000000 8.333333
## Minimum 1.000000 1.000000 1.000000 1.000000 1.000000 1.000000
## Maximum 7.000000 7.000000 7.000000 7.000000 7.000000 10.00000
## ────────────────────────────────────────────────────────────────────────────────────────────────
qol.plot<- qol[,c("RELATVS", "FRIENDS", "NEIGHBRS", "PHYSHLTH", "MNTLHLTH", "HB_EXERCISE", "POLICY", "SLEEP",
"VACADAYS", "AGE", "SOCIAL_CLIMATE", "SOC_HELP", "SOC_NEEDS", "OVERALL_SOC", "SOCIAL_ISO", "TIME_ADEQ")]
#Calculating correlations and CIs
cor1 <- cor.mtest(qol.plot, use="pairwise.complete.obs", conf.level = 0.95)
cor1
## $p
## [,1] [,2] [,3] [,4] [,5]
## [1,] 0.000000e+00 3.215465e-02 5.999927e-02 1.934269e-01 5.952267e-02
## [2,] 3.215465e-02 0.000000e+00 3.911301e-62 5.313607e-08 1.252809e-11
## [3,] 5.999927e-02 3.911301e-62 0.000000e+00 1.423136e-29 4.086942e-44
## [4,] 1.934269e-01 5.313607e-08 1.423136e-29 0.000000e+00 0.000000e+00
## [5,] 5.952267e-02 1.252809e-11 4.086942e-44 0.000000e+00 0.000000e+00
## [6,] 5.826287e-02 5.172306e-07 1.275818e-17 0.000000e+00 8.936008e-184
## [7,] 2.837149e-02 2.343036e-01 2.548289e-09 8.440250e-44 1.365565e-29
## [8,] 9.754642e-01 7.601689e-01 8.741335e-01 1.759412e-01 6.190283e-01
## [9,] 9.874884e-01 1.238577e-01 8.907930e-04 3.705966e-21 1.329868e-32
## [10,] 6.347267e-01 5.388044e-06 6.950880e-06 6.526822e-16 9.135334e-119
## [11,] 4.286784e-11 1.429212e-40 1.255829e-232 2.566323e-154 0.000000e+00
## [12,] 1.720690e-07 3.440084e-25 9.123312e-176 2.377279e-124 3.141046e-191
## [13,] 5.654574e-05 1.121040e-20 5.579414e-72 7.420066e-151 3.698525e-229
## [14,] 1.394672e-09 3.459186e-37 3.094297e-199 4.315428e-195 0.000000e+00
## [15,] 2.836705e-10 6.641116e-26 1.305859e-105 1.056862e-150 0.000000e+00
## [16,] 9.775218e-02 8.965668e-09 1.463856e-17 2.329552e-70 2.148150e-293
## [,6] [,7] [,8] [,9] [,10]
## [1,] 5.826287e-02 2.837149e-02 0.9754642 9.874884e-01 6.347267e-01
## [2,] 5.172306e-07 2.343036e-01 0.7601689 1.238577e-01 5.388044e-06
## [3,] 1.275818e-17 2.548289e-09 0.8741335 8.907930e-04 6.950880e-06
## [4,] 0.000000e+00 8.440250e-44 0.1759412 3.705966e-21 6.526822e-16
## [5,] 8.936008e-184 1.365565e-29 0.6190283 1.329868e-32 9.135334e-119
## [6,] 0.000000e+00 8.466469e-26 0.2931558 2.663379e-18 5.392800e-01
## [7,] 8.466469e-26 0.000000e+00 0.6109871 2.386449e-03 3.984738e-03
## [8,] 2.931558e-01 6.109871e-01 0.0000000 6.502847e-01 7.614760e-01
## [9,] 2.663379e-18 2.386449e-03 0.6502847 0.000000e+00 2.182615e-40
## [10,] 5.392800e-01 3.984738e-03 0.7614760 2.182615e-40 0.000000e+00
## [11,] 2.059844e-114 5.139742e-91 0.2277202 3.035586e-18 2.733732e-51
## [12,] 2.254998e-64 4.515870e-100 0.2067218 1.645054e-09 1.312649e-24
## [13,] 4.780568e-116 1.362941e-208 0.3608242 1.049395e-20 6.485392e-86
## [14,] 1.404918e-132 6.402072e-178 0.1962249 3.658058e-21 2.996878e-70
## [15,] 5.517294e-92 1.431926e-55 0.4857339 1.059702e-22 3.878532e-22
## [16,] 1.179663e-86 2.945834e-52 0.4647070 4.555960e-75 0.000000e+00
## [,11] [,12] [,13] [,14] [,15]
## [1,] 4.286784e-11 1.720690e-07 5.654574e-05 1.394672e-09 2.836705e-10
## [2,] 1.429212e-40 3.440084e-25 1.121040e-20 3.459186e-37 6.641116e-26
## [3,] 1.255829e-232 9.123312e-176 5.579414e-72 3.094297e-199 1.305859e-105
## [4,] 2.566323e-154 2.377279e-124 7.420066e-151 4.315428e-195 1.056862e-150
## [5,] 0.000000e+00 3.141046e-191 3.698525e-229 0.000000e+00 0.000000e+00
## [6,] 2.059844e-114 2.254998e-64 4.780568e-116 1.404918e-132 5.517294e-92
## [7,] 5.139742e-91 4.515870e-100 1.362941e-208 6.402072e-178 1.431926e-55
## [8,] 2.277202e-01 2.067218e-01 3.608242e-01 1.962249e-01 4.857339e-01
## [9,] 3.035586e-18 1.645054e-09 1.049395e-20 3.658058e-21 1.059702e-22
## [10,] 2.733732e-51 1.312649e-24 6.485392e-86 2.996878e-70 3.878532e-22
## [11,] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [12,] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [13,] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [14,] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [15,] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [16,] 3.262273e-206 2.124869e-139 8.853780e-226 9.199175e-261 5.222997e-206
## [,16]
## [1,] 9.775218e-02
## [2,] 8.965668e-09
## [3,] 1.463856e-17
## [4,] 2.329552e-70
## [5,] 2.148150e-293
## [6,] 1.179663e-86
## [7,] 2.945834e-52
## [8,] 4.647070e-01
## [9,] 4.555960e-75
## [10,] 0.000000e+00
## [11,] 3.262273e-206
## [12,] 2.124869e-139
## [13,] 8.853780e-226
## [14,] 9.199175e-261
## [15,] 5.222997e-206
## [16,] 0.000000e+00
##
## $lowCI
## [,1] [,2] [,3] [,4] [,5]
## [1,] 1.0000000000 0.001669231 -0.0007246962 -0.006027305 -0.0006915401
## [2,] 0.0016692306 1.000000000 0.1337533526 0.031796310 0.0439817480
## [3,] -0.0007246962 0.133753353 1.0000000000 0.085200815 0.1091981739
## [4,] -0.0060273052 0.031796310 0.0852008148 1.000000000 0.4773586663
## [5,] -0.0006915401 0.043981748 0.1091981739 0.477358666 1.0000000000
## [6,] -0.0352587509 0.027995285 0.0602595073 0.434394502 0.2422886362
## [7,] 0.0021273565 -0.007057431 0.0366351213 0.108555977 0.0851380944
## [8,] -0.0182560724 -0.020763968 -0.0165258457 -0.030271088 -0.0133700213
## [9,] -0.0179407990 -0.003885254 0.0125839214 0.068747891 0.0912141959
## [10,] -0.0135559802 0.023629114 0.0231508380 -0.091080010 0.1913355413
## [11,] 0.0425684220 0.104260878 0.2761679654 0.221946348 0.3273151185
## [12,] 0.0299981060 0.077090586 0.2383610905 0.197423402 0.2485038914
## [13,] 0.0189883632 0.067569392 0.1461652989 0.219208671 0.2731987412
## [14,] 0.0375856639 0.098712387 0.2544732125 0.250930086 0.3321144178
## [15,] -0.0757716051 -0.114146163 -0.2156526628 -0.252837669 -0.3878170198
## [16,] -0.0029441342 0.036668148 0.0635898458 0.151495908 0.3259108342
## [,6] [,7] [,8] [,9] [,10]
## [1,] -0.03525875 0.002127357 -0.01825607 -0.017940799 -0.013555980
## [2,] 0.02799529 -0.007057431 -0.02076397 -0.003885254 0.023629114
## [3,] 0.06025951 0.036635121 -0.01652585 0.012583921 0.023150838
## [4,] 0.43439450 0.108555977 -0.03027109 0.068747891 -0.091080010
## [5,] 0.24228864 0.085138094 -0.01337002 0.091214196 0.191335541
## [6,] 1.00000000 0.078003722 -0.02754207 0.062248062 -0.012259739
## [7,] 0.07800372 1.000000000 -0.02261264 0.009931530 0.008385352
## [8,] -0.02754207 -0.022612639 1.00000000 -0.022229729 -0.015119888
## [9,] 0.06224806 0.009931530 -0.02222973 1.000000000 0.104007564
## [10,] -0.01225974 0.008385352 -0.01511989 0.104007564 1.000000000
## [11,] 0.18888348 0.166801621 -0.02910329 0.062421918 0.119586910
## [12,] 0.13667192 0.175690653 -0.02959747 0.037619621 0.075643354
## [13,] 0.19031482 0.260751253 -0.02640969 0.068128598 0.160875798
## [14,] 0.20439151 0.239643909 -0.02985545 0.069084984 0.143394628
## [15,] -0.20196873 -0.160939417 -0.01160055 -0.108405275 -0.106038685
## [16,] 0.17059514 0.127678962 -0.01188539 0.158329236 0.446828032
## [,11] [,12] [,13] [,14] [,15] [,16]
## [1,] 0.04256842 0.02999811 0.01898836 0.03758566 -0.07577161 -0.002944134
## [2,] 0.10426088 0.07709059 0.06756939 0.09871239 -0.11414616 0.036668148
## [3,] 0.27616797 0.23836109 0.14616530 0.25447321 -0.21565266 0.063589846
## [4,] 0.22194635 0.19742340 0.21920867 0.25093009 -0.25283767 0.151495908
## [5,] 0.32731512 0.24850389 0.27319874 0.33211442 -0.38781702 0.325910834
## [6,] 0.18888348 0.13667192 0.19031482 0.20439151 -0.20196873 0.170595142
## [7,] 0.16680162 0.17569065 0.26075125 0.23964391 -0.16093942 0.127678962
## [8,] -0.02910329 -0.02959747 -0.02640969 -0.02985545 -0.01160055 -0.011885388
## [9,] 0.06242192 0.03761962 0.06812860 0.06908498 -0.10840527 0.158329236
## [10,] 0.11958691 0.07564335 0.16087580 0.14339463 -0.10603869 0.446828032
## [11,] 1.00000000 0.68058153 0.60702780 0.88504036 -0.63356454 0.273462288
## [12,] 0.68058153 1.00000000 0.52706492 0.83942621 -0.52992973 0.222261645
## [13,] 0.60702780 0.52706492 1.00000000 0.84784458 -0.53461617 0.286533362
## [14,] 0.88504036 0.83942621 0.84784458 1.00000000 -0.65224897 0.308297541
## [15,] -0.63356454 -0.52992973 -0.53461617 -0.65224897 1.00000000 -0.307983543
## [16,] 0.27346229 0.22226164 0.28653336 0.30829754 -0.30798354 1.000000000
##
## $uppCI
## [,1] [,2] [,3] [,4] [,5]
## [1,] 1.000000000 0.03749681 0.03514803 0.029791004 0.03513018
## [2,] 0.037496814 1.00000000 0.16879475 0.067516406 0.07966220
## [3,] 0.035148026 0.16879475 1.00000000 0.120651958 0.14446131
## [4,] 0.029791004 0.06751641 0.12065196 1.000000000 0.50443222
## [5,] 0.035130180 0.07966220 0.14446131 0.504432222 1.00000000
## [6,] 0.000606317 0.06378057 0.09592397 0.462924065 0.27562285
## [7,] 0.038026339 0.02883370 0.07244473 0.143769284 0.12054993
## [8,] 0.017692105 0.01516882 0.01943133 0.005543812 0.02245710
## [9,] 0.018230132 0.03226301 0.04873470 0.104527100 0.12684595
## [10,] 0.022227884 0.05933718 0.05888111 -0.055625203 0.22544221
## [11,] 0.078449879 0.13972613 0.30910608 0.255809867 0.35900359
## [12,] 0.065902882 0.11273880 0.27201151 0.231656071 0.28187279
## [13,] 0.054946419 0.10329870 0.18121487 0.253108640 0.30609558
## [14,] 0.073444136 0.13417894 0.28780447 0.284226754 0.36364599
## [15,] -0.039907513 -0.07851067 -0.18107725 -0.218954560 -0.35690519
## [16,] 0.034990178 0.07448257 0.10129182 0.188238147 0.35931353
## [,6] [,7] [,8] [,9] [,10] [,11]
## [1,] 0.000606317 0.03802634 0.017692105 0.01823013 0.02222788 0.078449879
## [2,] 0.063780574 0.02883370 0.015168824 0.03226301 0.05933718 0.139726127
## [3,] 0.095923965 0.07244473 0.019431332 0.04873470 0.05888111 0.309106084
## [4,] 0.462924065 0.14376928 0.005543812 0.10452710 -0.05562520 0.255809867
## [5,] 0.275622851 0.12054993 0.022457096 0.12684595 0.22544221 0.359003588
## [6,] 1.000000000 0.11350741 0.008313940 0.09810407 0.02343899 0.223285637
## [7,] 0.113507413 1.00000000 0.013295524 0.04602959 0.04410507 0.201583964
## [8,] 0.008313940 0.01329552 1.000000000 0.01387913 0.02065990 0.006929977
## [9,] 0.098104068 0.04602959 0.013879134 1.00000000 0.13948695 0.098455509
## [10,] 0.023438993 0.04410507 0.020659902 0.13948695 1.00000000 0.154784410
## [11,] 0.223285637 0.20158396 0.006929977 0.09845551 0.15478441 1.000000000
## [12,] 0.171720919 0.21032550 0.006406480 0.07374456 0.11117483 0.699375850
## [13,] 0.224688699 0.29397365 0.009616339 0.10411552 0.19559848 0.629205052
## [14,] 0.238514560 0.27323445 0.006130060 0.10503069 0.17829239 0.892571045
## [15,] -0.167295904 -0.12570768 0.024403732 -0.07246749 -0.07047812 -0.611589632
## [16,] 0.207115326 0.16476113 0.026023043 0.19524990 0.47656745 0.308286398
## [,12] [,13] [,14] [,15] [,16]
## [1,] 0.06590288 0.054946419 0.07344414 -0.03990751 0.03499018
## [2,] 0.11273880 0.103298700 0.13417894 -0.07851067 0.07448257
## [3,] 0.27201151 0.181214873 0.28780447 -0.18107725 0.10129182
## [4,] 0.23165607 0.253108640 0.28422675 -0.21895456 0.18823815
## [5,] 0.28187279 0.306095576 0.36364599 -0.35690519 0.35931353
## [6,] 0.17172092 0.224688699 0.23851456 -0.16729590 0.20711533
## [7,] 0.21032550 0.293973653 0.27323445 -0.12570768 0.16476113
## [8,] 0.00640648 0.009616339 0.00613006 0.02440373 0.02602304
## [9,] 0.07374456 0.104115524 0.10503069 -0.07246749 0.19524990
## [10,] 0.11117483 0.195598478 0.17829239 -0.07047812 0.47656745
## [11,] 0.69937585 0.629205052 0.89257105 -0.61158963 0.30828640
## [12,] 1.00000000 0.552483596 0.84969879 -0.50364855 0.25809345
## [13,] 0.55248360 1.000000000 0.85762504 -0.50850281 0.32105694
## [14,] 0.84969879 0.857625041 1.00000000 -0.63116685 0.34228582
## [15,] -0.50364855 -0.508502808 -0.63116685 1.00000000 -0.27317254
## [16,] 0.25809345 0.321056941 0.34228582 -0.27317254 1.00000000
#Correlation Matrix
corrplot(cor(qol.plot, use="pairwise.complete.obs"), method="number", type="upper",
addCoef.col = "orchid4", tl.col="hotpink3", tl.srt=40, p.mat = cor1$p,
sig.level = 0.05, insig = "blank", diag=FALSE)
#correlations
qol.corr = rcorr(as.matrix(qol.plot,use="pairwise.complete.obs"), type = c("pearson"))
qol.corr
## RELATVS FRIENDS NEIGHBRS PHYSHLTH MNTLHLTH HB_EXERCISE POLICY
## RELATVS 1.00 0.02 0.02 0.01 0.02 -0.02 0.02
## FRIENDS 0.02 1.00 0.15 0.05 0.06 0.05 0.01
## NEIGHBRS 0.02 0.15 1.00 0.10 0.13 0.08 0.05
## PHYSHLTH 0.01 0.05 0.10 1.00 0.49 0.45 0.13
## MNTLHLTH 0.02 0.06 0.13 0.49 1.00 0.26 0.10
## HB_EXERCISE -0.02 0.05 0.08 0.45 0.26 1.00 0.10
## POLICY 0.02 0.01 0.05 0.13 0.10 0.10 1.00
## SLEEP 0.00 0.00 0.00 -0.01 0.00 -0.01 0.00
## VACADAYS 0.00 0.01 0.03 0.09 0.11 0.08 0.03
## AGE 0.00 0.04 0.04 -0.07 0.21 0.01 0.03
## SOCIAL_CLIMATE 0.06 0.12 0.29 0.24 0.34 0.21 0.18
## SOC_HELP 0.05 0.09 0.26 0.21 0.27 0.15 0.19
## SOC_NEEDS 0.04 0.09 0.16 0.24 0.29 0.21 0.28
## OVERALL_SOC 0.06 0.12 0.27 0.27 0.35 0.22 0.26
## SOCIAL_ISO -0.06 -0.10 -0.20 -0.24 -0.37 -0.18 -0.14
## TIME_ADEQ 0.02 0.06 0.08 0.17 0.34 0.19 0.15
## SLEEP VACADAYS AGE SOCIAL_CLIMATE SOC_HELP SOC_NEEDS
## RELATVS 0.00 0.00 0.00 0.06 0.05 0.04
## FRIENDS 0.00 0.01 0.04 0.12 0.09 0.09
## NEIGHBRS 0.00 0.03 0.04 0.29 0.26 0.16
## PHYSHLTH -0.01 0.09 -0.07 0.24 0.21 0.24
## MNTLHLTH 0.00 0.11 0.21 0.34 0.27 0.29
## HB_EXERCISE -0.01 0.08 0.01 0.21 0.15 0.21
## POLICY 0.00 0.03 0.03 0.18 0.19 0.28
## SLEEP 1.00 0.00 0.00 -0.01 -0.01 -0.01
## VACADAYS 0.00 1.00 0.12 0.08 0.06 0.09
## AGE 0.00 0.12 1.00 0.14 0.09 0.18
## SOCIAL_CLIMATE -0.01 0.08 0.14 1.00 0.69 0.62
## SOC_HELP -0.01 0.06 0.09 0.69 1.00 0.54
## SOC_NEEDS -0.01 0.09 0.18 0.62 0.54 1.00
## OVERALL_SOC -0.01 0.09 0.16 0.89 0.84 0.85
## SOCIAL_ISO 0.01 -0.09 -0.09 -0.62 -0.52 -0.52
## TIME_ADEQ 0.01 0.18 0.46 0.29 0.24 0.30
## OVERALL_SOC SOCIAL_ISO TIME_ADEQ
## RELATVS 0.06 -0.06 0.02
## FRIENDS 0.12 -0.10 0.06
## NEIGHBRS 0.27 -0.20 0.08
## PHYSHLTH 0.27 -0.24 0.17
## MNTLHLTH 0.35 -0.37 0.34
## HB_EXERCISE 0.22 -0.18 0.19
## POLICY 0.26 -0.14 0.15
## SLEEP -0.01 0.01 0.01
## VACADAYS 0.09 -0.09 0.18
## AGE 0.16 -0.09 0.46
## SOCIAL_CLIMATE 0.89 -0.62 0.29
## SOC_HELP 0.84 -0.52 0.24
## SOC_NEEDS 0.85 -0.52 0.30
## OVERALL_SOC 1.00 -0.64 0.33
## SOCIAL_ISO -0.64 1.00 -0.29
## TIME_ADEQ 0.33 -0.29 1.00
##
## n
## RELATVS FRIENDS NEIGHBRS PHYSHLTH MNTLHLTH HB_EXERCISE POLICY
## RELATVS 12005 11962 11934 11974 11968 11939 11914
## FRIENDS 11962 12015 11948 11984 11978 11949 11926
## NEIGHBRS 11934 11948 12001 11969 11963 11934 11912
## PHYSHLTH 11974 11984 11969 12098 12076 12042 12001
## MNTLHLTH 11968 11978 11963 12076 12092 12036 11996
## HB_EXERCISE 11939 11949 11934 12042 12036 12062 11968
## POLICY 11914 11926 11912 12001 11996 11968 12032
## SLEEP 11891 11901 11885 11976 11971 11950 11917
## VACADAYS 11745 11755 11736 11824 11817 11799 11774
## AGE 12000 12010 11996 12093 12087 12057 12027
## SOCIAL_CLIMATE 11848 11856 11841 11914 11910 11904 11854
## SOC_HELP 11865 11875 11860 11933 11927 11921 11873
## SOC_NEEDS 11852 11862 11847 11921 11916 11909 11862
## OVERALL_SOC 11877 11887 11872 11946 11940 11934 11886
## SOCIAL_ISO 11867 11877 11862 11936 11930 11925 11876
## TIME_ADEQ 10673 10680 10664 10735 10728 10714 10702
## SLEEP VACADAYS AGE SOCIAL_CLIMATE SOC_HELP SOC_NEEDS
## RELATVS 11891 11745 12000 11848 11865 11852
## FRIENDS 11901 11755 12010 11856 11875 11862
## NEIGHBRS 11885 11736 11996 11841 11860 11847
## PHYSHLTH 11976 11824 12093 11914 11933 11921
## MNTLHLTH 11971 11817 12087 11910 11927 11916
## HB_EXERCISE 11950 11799 12057 11904 11921 11909
## POLICY 11917 11774 12027 11854 11873 11862
## SLEEP 12007 11785 12003 11832 11851 11838
## VACADAYS 11785 11853 11848 11682 11702 11690
## AGE 12003 11848 12126 11941 11960 11948
## SOCIAL_CLIMATE 11832 11682 11941 11946 11938 11927
## SOC_HELP 11851 11702 11960 11938 11965 11941
## SOC_NEEDS 11838 11690 11948 11927 11941 11953
## OVERALL_SOC 11863 11713 11973 11946 11965 11953
## SOCIAL_ISO 11853 11704 11963 11931 11949 11940
## TIME_ADEQ 10692 10579 10755 10617 10627 10622
## OVERALL_SOC SOCIAL_ISO TIME_ADEQ
## RELATVS 11877 11867 10673
## FRIENDS 11887 11877 10680
## NEIGHBRS 11872 11862 10664
## PHYSHLTH 11946 11936 10735
## MNTLHLTH 11940 11930 10728
## HB_EXERCISE 11934 11925 10714
## POLICY 11886 11876 10702
## SLEEP 11863 11853 10692
## VACADAYS 11713 11704 10579
## AGE 11973 11963 10755
## SOCIAL_CLIMATE 11946 11931 10617
## SOC_HELP 11965 11949 10627
## SOC_NEEDS 11953 11940 10622
## OVERALL_SOC 11978 11959 10635
## SOCIAL_ISO 11959 11968 10629
## TIME_ADEQ 10635 10629 10759
##
## P
## RELATVS FRIENDS NEIGHBRS PHYSHLTH MNTLHLTH HB_EXERCISE POLICY
## RELATVS 0.0322 0.0600 0.1934 0.0595 0.0583 0.0284
## FRIENDS 0.0322 0.0000 0.0000 0.0000 0.0000 0.2343
## NEIGHBRS 0.0600 0.0000 0.0000 0.0000 0.0000 0.0000
## PHYSHLTH 0.1934 0.0000 0.0000 0.0000 0.0000 0.0000
## MNTLHLTH 0.0595 0.0000 0.0000 0.0000 0.0000 0.0000
## HB_EXERCISE 0.0583 0.0000 0.0000 0.0000 0.0000 0.0000
## POLICY 0.0284 0.2343 0.0000 0.0000 0.0000 0.0000
## SLEEP 0.9755 0.7602 0.8741 0.1759 0.6190 0.2932 0.6110
## VACADAYS 0.9875 0.1239 0.0009 0.0000 0.0000 0.0000 0.0024
## AGE 0.6347 0.0000 0.0000 0.0000 0.0000 0.5393 0.0040
## SOCIAL_CLIMATE 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
## SOC_HELP 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
## SOC_NEEDS 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
## OVERALL_SOC 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
## SOCIAL_ISO 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
## TIME_ADEQ 0.0978 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
## SLEEP VACADAYS AGE SOCIAL_CLIMATE SOC_HELP SOC_NEEDS
## RELATVS 0.9755 0.9875 0.6347 0.0000 0.0000 0.0000
## FRIENDS 0.7602 0.1239 0.0000 0.0000 0.0000 0.0000
## NEIGHBRS 0.8741 0.0009 0.0000 0.0000 0.0000 0.0000
## PHYSHLTH 0.1759 0.0000 0.0000 0.0000 0.0000 0.0000
## MNTLHLTH 0.6190 0.0000 0.0000 0.0000 0.0000 0.0000
## HB_EXERCISE 0.2932 0.0000 0.5393 0.0000 0.0000 0.0000
## POLICY 0.6110 0.0024 0.0040 0.0000 0.0000 0.0000
## SLEEP 0.6503 0.7615 0.2277 0.2067 0.3608
## VACADAYS 0.6503 0.0000 0.0000 0.0000 0.0000
## AGE 0.7615 0.0000 0.0000 0.0000 0.0000
## SOCIAL_CLIMATE 0.2277 0.0000 0.0000 0.0000 0.0000
## SOC_HELP 0.2067 0.0000 0.0000 0.0000 0.0000
## SOC_NEEDS 0.3608 0.0000 0.0000 0.0000 0.0000
## OVERALL_SOC 0.1962 0.0000 0.0000 0.0000 0.0000 0.0000
## SOCIAL_ISO 0.4857 0.0000 0.0000 0.0000 0.0000 0.0000
## TIME_ADEQ 0.4647 0.0000 0.0000 0.0000 0.0000 0.0000
## OVERALL_SOC SOCIAL_ISO TIME_ADEQ
## RELATVS 0.0000 0.0000 0.0978
## FRIENDS 0.0000 0.0000 0.0000
## NEIGHBRS 0.0000 0.0000 0.0000
## PHYSHLTH 0.0000 0.0000 0.0000
## MNTLHLTH 0.0000 0.0000 0.0000
## HB_EXERCISE 0.0000 0.0000 0.0000
## POLICY 0.0000 0.0000 0.0000
## SLEEP 0.1962 0.4857 0.4647
## VACADAYS 0.0000 0.0000 0.0000
## AGE 0.0000 0.0000 0.0000
## SOCIAL_CLIMATE 0.0000 0.0000 0.0000
## SOC_HELP 0.0000 0.0000 0.0000
## SOC_NEEDS 0.0000 0.0000 0.0000
## OVERALL_SOC 0.0000 0.0000
## SOCIAL_ISO 0.0000 0.0000
## TIME_ADEQ 0.0000 0.0000
print(qol.corr$P,3)
## RELATVS FRIENDS NEIGHBRS PHYSHLTH MNTLHLTH HB_EXERCISE
## RELATVS NA 3.22e-02 6.00e-02 1.93e-01 5.95e-02 5.83e-02
## FRIENDS 3.22e-02 NA 0.00e+00 5.31e-08 1.25e-11 5.17e-07
## NEIGHBRS 6.00e-02 0.00e+00 NA 0.00e+00 0.00e+00 0.00e+00
## PHYSHLTH 1.93e-01 5.31e-08 0.00e+00 NA 0.00e+00 0.00e+00
## MNTLHLTH 5.95e-02 1.25e-11 0.00e+00 0.00e+00 NA 0.00e+00
## HB_EXERCISE 5.83e-02 5.17e-07 0.00e+00 0.00e+00 0.00e+00 NA
## POLICY 2.84e-02 2.34e-01 2.55e-09 0.00e+00 0.00e+00 0.00e+00
## SLEEP 9.75e-01 7.60e-01 8.74e-01 1.76e-01 6.19e-01 2.93e-01
## VACADAYS 9.87e-01 1.24e-01 8.91e-04 0.00e+00 0.00e+00 0.00e+00
## AGE 6.35e-01 5.39e-06 6.95e-06 6.66e-16 0.00e+00 5.39e-01
## SOCIAL_CLIMATE 4.29e-11 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SOC_HELP 1.72e-07 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SOC_NEEDS 5.65e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## OVERALL_SOC 1.39e-09 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SOCIAL_ISO 2.84e-10 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## TIME_ADEQ 9.78e-02 8.97e-09 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## POLICY SLEEP VACADAYS AGE SOCIAL_CLIMATE SOC_HELP
## RELATVS 2.84e-02 0.975 9.87e-01 6.35e-01 4.29e-11 1.72e-07
## FRIENDS 2.34e-01 0.760 1.24e-01 5.39e-06 0.00e+00 0.00e+00
## NEIGHBRS 2.55e-09 0.874 8.91e-04 6.95e-06 0.00e+00 0.00e+00
## PHYSHLTH 0.00e+00 0.176 0.00e+00 6.66e-16 0.00e+00 0.00e+00
## MNTLHLTH 0.00e+00 0.619 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## HB_EXERCISE 0.00e+00 0.293 0.00e+00 5.39e-01 0.00e+00 0.00e+00
## POLICY NA 0.611 2.39e-03 3.98e-03 0.00e+00 0.00e+00
## SLEEP 6.11e-01 NA 6.50e-01 7.61e-01 2.28e-01 2.07e-01
## VACADAYS 2.39e-03 0.650 NA 0.00e+00 0.00e+00 1.65e-09
## AGE 3.98e-03 0.761 0.00e+00 NA 0.00e+00 0.00e+00
## SOCIAL_CLIMATE 0.00e+00 0.228 0.00e+00 0.00e+00 NA 0.00e+00
## SOC_HELP 0.00e+00 0.207 1.65e-09 0.00e+00 0.00e+00 NA
## SOC_NEEDS 0.00e+00 0.361 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## OVERALL_SOC 0.00e+00 0.196 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SOCIAL_ISO 0.00e+00 0.486 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## TIME_ADEQ 0.00e+00 0.465 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SOC_NEEDS OVERALL_SOC SOCIAL_ISO TIME_ADEQ
## RELATVS 5.65e-05 1.39e-09 2.84e-10 9.78e-02
## FRIENDS 0.00e+00 0.00e+00 0.00e+00 8.97e-09
## NEIGHBRS 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## PHYSHLTH 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## MNTLHLTH 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## HB_EXERCISE 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## POLICY 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SLEEP 3.61e-01 1.96e-01 4.86e-01 4.65e-01
## VACADAYS 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## AGE 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SOCIAL_CLIMATE 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SOC_HELP 0.00e+00 0.00e+00 0.00e+00 0.00e+00
## SOC_NEEDS NA 0.00e+00 0.00e+00 0.00e+00
## OVERALL_SOC 0.00e+00 NA 0.00e+00 0.00e+00
## SOCIAL_ISO 0.00e+00 0.00e+00 NA 0.00e+00
## TIME_ADEQ 0.00e+00 0.00e+00 0.00e+00 NA
#coefficients
qol.plot.coeff = qol.corr$r
#p-values
qol.plot.p = qol.corr$P
#plot
corrplot(qol.plot.coeff, type="upper",p.mat = qol.corr$P, sig.level = 0.05, insig = "blank")
#chart.Correlation(qol.plot, histogram=TRUE)