setwd("~/Google Drive/Research/Engage data")
library(readr)
library(GGally)
library(ggplot2)
library(PerformanceAnalytics)
library(corrplot)
library(plyr)
library(dplyr)
library(psych)
library(jmv)
library(Hmisc)
library(rmarkdown)
library(jmv)
qol <- read_csv("CIW-NovaScotiaQualityLife-SHAREDSurveyData-May2020.csv")
### NEW DATA FRAMES ###
### ALL TIME-USE AND WELL-BEING VARIABLES
qol.time<- qol[,c("RUSHED","SLEEP","VACADAYS", "WBI1MNTL", "WBI2PHYS","WBI3LEIS","WBI4PERS",
"WBE2EDUC", "WBI7BAL", "WBI8TIME", "WBE3ARTS", "WBE4REC", "WBE5NBHD","WBE6ENV",
"WBE7GOV", "WBE8DEM", "WBI5FIN", "WBI6WORK", "LIFEWRTH", "LIFESAT", "TIME_ADEQ",
"CV", "DE", "ED", "ENV", "HP", "LC", "LS", "TU", "CIW_WELLB", "MNTLHLTH", "PA_TEAMS",
"PA_INDIV", "PA_VIGOROUS", "PA_LIGHTEX", "SOC_FRND", "SOC_FILM", "SOC_CLUB",
"SOC_SPEC", "HM_READ", "HM_CARDS", "HM_PUZZL", "HM_MUSIC", "CULPAR_1", "CULPAR_2",
"CULPAR_3", "CULPAR_4", "CULPAR_5", "HM_SURF", "HM_CGAME", "HM_CSOC", "WL_ENHANCE",
"WK_INTERF","PL_INTERF", "WL_BALANCE", "BELONGING", "SOCIAL_CLIMATE", "SOC_HELP",
"SOC_NEEDS", "OVERALL_SOC", "SOCIAL_ISO", "TIMEREST", "TIMESELF", "TIMESOC","TIMEPHYS",
"TIMEMEAL","TIMECOMM", "TIMENURT", "TIMECHOR", "TIMEKIDS", "TIMESPOU", "TIMEFAML",
"TIMERELA", "TIMECREA")]
## TIME ADEQUACY AND WELL-BEING
qol.timead<- qol[,c("RUSHED","SLEEP","VACADAYS", "WBI1MNTL", "WBI2PHYS","WBI3LEIS","WBI4PERS",
"WBE2EDUC", "WBI7BAL", "WBI8TIME", "WBE3ARTS", "WBE4REC", "WBE5NBHD","WBE6ENV",
"WBE7GOV", "WBE8DEM", "WBI5FIN", "WBI6WORK", "LIFEWRTH", "LIFESAT", "TIME_ADEQ",
"TIMEREST", "TIMESELF", "TIMESOC","TIMEPHYS","TIMEMEAL","TIMECOMM","TIMENURT",
"TIMECHOR", "TIMEKIDS", "TIMESPOU", "TIMEFAML", "TIMERELA", "TIMECREA")]
## JUST MAIN VARIABLES FOR EXPLORATION
qol.timewb<- qol[,c("WBI1MNTL","WBI3LEIS","WBI4PERS","WBI7BAL", "WBI8TIME","LIFEWRTH", "LIFESAT",
"TIMEREST", "TIMESELF", "TIMESOC","TIMEPHYS","TIMEMEAL","TIMECOMM","TIMENURT",
"TIMECHOR", "TIMEKIDS", "TIMESPOU", "TIMEFAML", "TIMERELA", "TIMECREA",
"TIME_ADEQ")]
write.csv(file = "qol.timewb", qol.timewb, row.names = FALSE)
qol.plot<- qol[,c("WBI1MNTL","LIFESAT", "TIMEREST", "TIMESELF", "TIMESOC","TIMEPHYS",
"TIMEMEAL","TIMECOMM","TIMENURT","TIMECHOR", "TIMEKIDS", "TIMESPOU", "TIMEFAML",
"TIMERELA", "TIMECREA")]
#Calculating correlations and CIs
cor1 <- cor.mtest(qol.plot, use="pairwise.complete.obs", conf.level = 0.95)
cor1
## $p
## [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
## [1,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [2,] 0.000000e+00 0.000000e+00 0 0 0 0 0 1.193374e-301
## [3,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [4,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [5,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [6,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [7,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [8,] 0.000000e+00 1.193374e-301 0 0 0 0 0 0.000000e+00
## [9,] 0.000000e+00 2.436526e-266 0 0 0 0 0 0.000000e+00
## [10,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [11,] 8.667015e-206 2.365479e-144 0 0 0 0 0 0.000000e+00
## [12,] 1.969685e-299 3.693417e-184 0 0 0 0 0 0.000000e+00
## [13,] 0.000000e+00 1.758420e-296 0 0 0 0 0 0.000000e+00
## [14,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [15,] 0.000000e+00 0.000000e+00 0 0 0 0 0 0.000000e+00
## [,9] [,10] [,11] [,12] [,13] [,14] [,15]
## [1,] 0.000000e+00 0 8.667015e-206 1.969685e-299 0.00000e+00 0 0
## [2,] 2.436526e-266 0 2.365479e-144 3.693417e-184 1.75842e-296 0 0
## [3,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [4,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [5,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [6,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [7,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [8,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [9,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [10,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [11,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [12,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [13,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [14,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
## [15,] 0.000000e+00 0 0.000000e+00 0.000000e+00 0.00000e+00 0 0
##
## $lowCI
## [,1] [,2] [,3] [,4] [,5] [,6] [,7]
## [1,] 1.0000000 0.6503112 0.3896476 0.4587819 0.4306533 0.3833817 0.4179616
## [2,] 0.6503112 1.0000000 0.3352043 0.4111502 0.3948222 0.3379507 0.3746726
## [3,] 0.3896476 0.3352043 1.0000000 0.7007720 0.6611818 0.6283677 0.6384943
## [4,] 0.4587819 0.4111502 0.7007720 1.0000000 0.7738981 0.6855387 0.6628907
## [5,] 0.4306533 0.3948222 0.6611818 0.7738981 1.0000000 0.7653320 0.6925410
## [6,] 0.3833817 0.3379507 0.6283677 0.6855387 0.7653320 1.0000000 0.7303352
## [7,] 0.4179616 0.3746726 0.6384943 0.6628907 0.6925410 0.7303352 1.0000000
## [8,] 0.3626263 0.3116578 0.5779468 0.6232760 0.7307743 0.7249283 0.6878853
## [9,] 0.3502083 0.2925069 0.5117175 0.5714358 0.6113719 0.6106149 0.5953252
## [10,] 0.4118653 0.3601606 0.6422818 0.6830346 0.7155290 0.6955778 0.7265935
## [11,] 0.2573990 0.2133228 0.4203675 0.4386540 0.4473567 0.4475977 0.4790489
## [12,] 0.3111816 0.2420860 0.5152828 0.5418706 0.5664458 0.5403054 0.5661375
## [13,] 0.3516921 0.3091505 0.5108574 0.5565439 0.5906194 0.5602428 0.5722852
## [14,] 0.3695118 0.3411391 0.5042342 0.5609415 0.6155472 0.5653482 0.5819196
## [15,] 0.3931652 0.3380764 0.5695831 0.6390291 0.6955481 0.6642542 0.6269785
## [,8] [,9] [,10] [,11] [,12] [,13] [,14]
## [1,] 0.3626263 0.3502083 0.4118653 0.2573990 0.3111816 0.3516921 0.3695118
## [2,] 0.3116578 0.2925069 0.3601606 0.2133228 0.2420860 0.3091505 0.3411391
## [3,] 0.5779468 0.5117175 0.6422818 0.4203675 0.5152828 0.5108574 0.5042342
## [4,] 0.6232760 0.5714358 0.6830346 0.4386540 0.5418706 0.5565439 0.5609415
## [5,] 0.7307743 0.6113719 0.7155290 0.4473567 0.5664458 0.5906194 0.6155472
## [6,] 0.7249283 0.6106149 0.6955778 0.4475977 0.5403054 0.5602428 0.5653482
## [7,] 0.6878853 0.5953252 0.7265935 0.4790489 0.5661375 0.5722852 0.5819196
## [8,] 1.0000000 0.6791524 0.6857285 0.4558924 0.5224383 0.5550057 0.5869437
## [9,] 0.6791524 1.0000000 0.6111038 0.4417602 0.4917312 0.4923819 0.5329356
## [10,] 0.6857285 0.6111038 1.0000000 0.4776542 0.6006796 0.6198185 0.6010493
## [11,] 0.4558924 0.4417602 0.4776542 1.0000000 0.5670394 0.4594997 0.4706509
## [12,] 0.5224383 0.4917312 0.6006796 0.5670394 1.0000000 0.6265780 0.5507890
## [13,] 0.5550057 0.4923819 0.6198185 0.4594997 0.6265780 1.0000000 0.5826900
## [14,] 0.5869437 0.5329356 0.6010493 0.4706509 0.5507890 0.5826900 1.0000000
## [15,] 0.6782303 0.6354149 0.6870902 0.4500707 0.5667639 0.5976692 0.6692053
## [,15]
## [1,] 0.3931652
## [2,] 0.3380764
## [3,] 0.5695831
## [4,] 0.6390291
## [5,] 0.6955481
## [6,] 0.6642542
## [7,] 0.6269785
## [8,] 0.6782303
## [9,] 0.6354149
## [10,] 0.6870902
## [11,] 0.4500707
## [12,] 0.5667639
## [13,] 0.5976692
## [14,] 0.6692053
## [15,] 1.0000000
##
## $uppCI
## [,1] [,2] [,3] [,4] [,5] [,6] [,7]
## [1,] 1.0000000 0.6703178 0.4193499 0.4863979 0.4591540 0.4133050 0.4468555
## [2,] 0.6703178 1.0000000 0.3663404 0.4402559 0.4243866 0.3690657 0.4048074
## [3,] 0.4193499 0.3663404 1.0000000 0.7183652 0.6806556 0.6493486 0.6590077
## [4,] 0.4863979 0.4402559 0.7183652 1.0000000 0.7877346 0.7038966 0.6823083
## [5,] 0.4591540 0.4243866 0.6806556 0.7877346 1.0000000 0.7796303 0.7105358
## [6,] 0.4133050 0.3690657 0.6493486 0.7038966 0.7796303 1.0000000 0.7464648
## [7,] 0.4468555 0.4048074 0.6590077 0.6823083 0.7105358 0.7464648 1.0000000
## [8,] 0.3931286 0.3434188 0.6010890 0.6445140 0.7468834 0.7413500 0.7061300
## [9,] 0.3810325 0.3246926 0.5374245 0.5948678 0.6331136 0.6324111 0.6177553
## [10,] 0.4409851 0.3907319 0.6626587 0.7015312 0.7324200 0.7134700 0.7429280
## [11,] 0.2904876 0.2471968 0.4492869 0.4670370 0.4754406 0.4756898 0.5060785
## [12,] 0.3430262 0.2753583 0.5409190 0.5665318 0.5901185 0.5650270 0.5898274
## [13,] 0.3824981 0.3409911 0.5366156 0.5805921 0.6132625 0.5841498 0.5956893
## [14,] 0.3998940 0.3722568 0.5302694 0.5848476 0.6371556 0.5890794 0.6049633
## [15,] 0.4228492 0.3692300 0.5930873 0.6595883 0.7134381 0.6836548 0.6480656
## [,8] [,9] [,10] [,11] [,12] [,13] [,14]
## [1,] 0.3931286 0.3810325 0.4409851 0.2904876 0.3430262 0.3824981 0.3998940
## [2,] 0.3434188 0.3246926 0.3907319 0.2471968 0.2753583 0.3409911 0.3722568
## [3,] 0.6010890 0.5374245 0.6626587 0.4492869 0.5409190 0.5366156 0.5302694
## [4,] 0.6445140 0.5948678 0.7015312 0.4670370 0.5665318 0.5805921 0.5848476
## [5,] 0.7468834 0.6331136 0.7324200 0.4754406 0.5901185 0.6132625 0.6371556
## [6,] 0.7413500 0.6324111 0.7134700 0.4756898 0.5650270 0.5841498 0.5890794
## [7,] 0.7061300 0.6177553 0.7429280 0.5060785 0.5898274 0.5956893 0.6049633
## [8,] 1.0000000 0.6978304 0.7041030 0.4837222 0.5478470 0.5791265 0.6097907
## [9,] 0.6978304 1.0000000 0.6328903 0.4700463 0.5182599 0.5188447 0.5579378
## [10,] 0.7041030 0.6328903 1.0000000 0.5047551 0.6229636 0.6412435 0.6233130
## [11,] 0.4837222 0.4700463 0.5047551 1.0000000 0.5907861 0.4872088 0.4979984
## [12,] 0.5478470 0.5182599 0.6229636 0.5907861 1.0000000 0.6477288 0.5751289
## [13,] 0.5791265 0.5188447 0.6412435 0.4872088 0.6477288 1.0000000 0.6057148
## [14,] 0.6097907 0.5579378 0.6233130 0.4979984 0.5751289 0.6057148 1.0000000
## [15,] 0.6969748 0.6561435 0.7054166 0.4780972 0.5904572 0.6200488 0.6883918
## [,15]
## [1,] 0.4228492
## [2,] 0.3692300
## [3,] 0.5930873
## [4,] 0.6595883
## [5,] 0.7134381
## [6,] 0.6836548
## [7,] 0.6480656
## [8,] 0.6969748
## [9,] 0.6561435
## [10,] 0.7054166
## [11,] 0.4780972
## [12,] 0.5904572
## [13,] 0.6200488
## [14,] 0.6883918
## [15,] 1.0000000
#Correlation Matrix
corrplot(cor(qol.plot, use="pairwise.complete.obs"), method="number", type="upper",
addCoef.col = "orchid4", tl.col="hotpink3", tl.srt=45, p.mat = cor1$p,
sig.level = 0.05, insig = "blank", diag=FALSE)
ggpairs(qol.plot) + theme_classic()
