COVID-19 - Data Analysis
Luis Anunciação, Lucas Barrozo
Note: This R markdown relates to data analysis of the book chapter entitled “O custo psicológico e econômico das medidas de distanciamento social: prejuízo de um salário mínimo e possível impacto na depressão”. To achieve the same results, please download the R dataset file at https://osf.io/4rwkp/. Please have in mind that this dataset file is the processed one, but not the original version. If you want to check the data processing scripts, please see them at https://osf.io/kyfmd/
If you have any questions or queries, please reach us out at lucas.psic.rj@gmail.com or luisfca@puc-rio.br
last updated: 14 July, 2021
1 Data Analysis
load(url("https://osf.io/4rwkp/download"))1.1 Packages
pacman::p_load(knitr,
kableExtra,
multcomp,
tidyverse,
janitor,
summarytools,
DataExplorer,
readxl,
patchwork)1.2 resultados
1.3 correlação Depressão e Ansiedade/Idade e Ansiedade
1.4 Distribuição de depressão e ansiedade
ds %>% select(total_ces_t1, total_gad_t1) %>% tableone::CreateTableOne(data=.)
Overall
n 516
total_ces_t1 (mean (SD)) 23.19 (12.86)
total_gad_t1 (mean (SD)) 9.56 (5.75) 1.4.1 Depressão
ggplot(ds, aes(total_ces_t1)) +
geom_density(fill = "blue", alpha = 0.2) +
theme_bw() +
labs(x = "Depressão", y = "Densidade")
1.4.2 Ansiedade
ggplot(ds, aes(total_gad_t1)) +
geom_density(fill = "blue", alpha = 0.2) +
theme_bw() +
labs(x = "Ansiedade", y = "Densidade")
1.4.3 Formato da base e feature selection
ds %>% select(sexo , idade , escolaridade , estado_civil , quantos_filhos , plano_saude , como_vive , fica_em_casa , tem_saido_de_casa , falou_familia_pessoalmente , falou_familia_telefone , falou_amigos_telefone , testou_coronavirus , conhece_alguem_coronavirus , segue_recomendacoes , concordo_covid_preocupacao_saude , concordo_quarentena , area_profissao , preocupacao_pagar_conta , agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos , rio_sp) %>% strtibble [516 x 21] (S3: tbl_df/tbl/data.frame)
$ sexo : Factor w/ 2 levels "Feminino","Masculino": 1 1 1 2 1 1 1 1 1 1 ...
$ idade : num [1:516] 22 22 23 23 20 21 21 26 25 36 ...
$ escolaridade : Factor w/ 7 levels "Ensino fundamental incompleto",..: 5 5 4 5 5 5 4 6 5 6 ...
$ estado_civil : Factor w/ 5 levels "Casado(a)","Separado(a) ou divorciado(a)",..: 3 3 3 3 3 3 3 3 3 4 ...
$ quantos_filhos : Factor w/ 5 levels "Nenhum","1 filho (a)",..: 1 1 1 1 1 1 1 1 1 2 ...
$ plano_saude : Factor w/ 2 levels "Não","Sim": 2 2 1 1 2 NA 2 2 2 2 ...
$ como_vive : Factor w/ 4 levels "Vivo com a minha família (por exemplo, com companheiro, filhos)",..: 1 3 3 4 3 3 1 3 3 1 ...
$ fica_em_casa : Factor w/ 6 levels "Entre 1 a 4 horas por dia",..: 5 6 3 2 6 6 6 6 6 6 ...
$ tem_saido_de_casa : Factor w/ 5 levels "Nenhuma","1 vez por dia",..: 2 1 1 2 2 1 1 2 1 1 ...
$ falou_familia_pessoalmente : Factor w/ 4 levels "Nenhuma","1 pessoa",..: 1 4 4 4 3 4 4 3 4 3 ...
$ falou_familia_telefone : Factor w/ 4 levels "1 pessoa","2 pessoas",..: 1 4 4 4 4 4 4 1 4 4 ...
$ falou_amigos_telefone : Factor w/ 4 levels "1 pessoa","2 pessoas",..: 1 4 4 3 3 3 3 4 1 3 ...
$ testou_coronavirus : Factor w/ 2 levels "Não","Sim": 1 1 1 1 1 1 1 1 1 1 ...
$ conhece_alguem_coronavirus : Factor w/ 2 levels "Não","Sim": 1 2 1 1 1 1 1 1 1 1 ...
$ segue_recomendacoes : Factor w/ 2 levels "Não","Sim": 2 2 2 2 2 2 2 2 2 2 ...
$ concordo_covid_preocupacao_saude : Factor w/ 5 levels "Eu discordo fortemente",..: 5 5 4 5 5 5 5 5 5 5 ...
$ concordo_quarentena : Factor w/ 5 levels "Eu discordo fortemente",..: 5 5 4 5 5 5 5 5 5 5 ...
$ area_profissao : Factor w/ 8 levels "Comércio (vendedor, barbeiro, etc)",..: NA NA 1 5 3 7 NA 6 NA 7 ...
$ preocupacao_pagar_conta : Factor w/ 2 levels "Não","Sim": NA 1 2 2 2 2 1 1 1 1 ...
$ agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos: Factor w/ 4 levels "Não, nenhum",..: 2 1 1 1 2 4 4 1 1 3 ...
$ rio_sp : Factor w/ 3 levels "Outros","RJ",..: NA NA NA NA NA 2 2 2 2 2 ...1.4.4 Missing cases (within ds)
ds %>% select(sexo , idade , escolaridade , estado_civil , quantos_filhos , plano_saude , como_vive , fica_em_casa , tem_saido_de_casa , falou_familia_pessoalmente , falou_familia_telefone , falou_amigos_telefone , testou_coronavirus , conhece_alguem_coronavirus , segue_recomendacoes , concordo_covid_preocupacao_saude , concordo_quarentena , area_profissao , preocupacao_pagar_conta , agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos , rio_sp) %>%
profile_missing()1.5 Dados descritivos (Tabela 1)
ds %>% select(sexo , idade , escolaridade , estado_civil , quantos_filhos , plano_saude , como_vive , fica_em_casa , tem_saido_de_casa , falou_familia_pessoalmente , falou_familia_telefone , falou_amigos_telefone , testou_coronavirus , conhece_alguem_coronavirus , segue_recomendacoes , concordo_covid_preocupacao_saude , concordo_quarentena , area_profissao , preocupacao_pagar_conta , agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos , rio_sp, total_ces_t1, total_gad_t1) %>% arsenal::tableby(~., data = .) %>% summary()| Overall (N=516) | |
|---|---|
| sexo | |
| Feminino | 366 (70.9%) |
| Masculino | 150 (29.1%) |
| idade | |
| Mean (SD) | 37.531 (14.090) |
| Range | 13.000 - 77.000 |
| escolaridade | |
| Ensino fundamental incompleto | 6 (1.2%) |
| Ensino fundamental completo | 12 (2.3%) |
| Ensino médio incompleto | 17 (3.3%) |
| Ensino médio completo | 80 (15.5%) |
| Ensino superior incompleto | 111 (21.5%) |
| Ensino superior completo | 117 (22.7%) |
| Pós-graduação | 173 (33.5%) |
| estado_civil | |
| Casado(a) | 189 (36.6%) |
| Separado(a) ou divorciado(a) | 38 (7.4%) |
| Solteiro(a) | 239 (46.3%) |
| União estável (mesmo que sem reconhecimento legal) | 44 (8.5%) |
| Viuvo(a) | 6 (1.2%) |
| quantos_filhos | |
| Nenhum | 283 (54.8%) |
| 1 filho (a) | 95 (18.4%) |
| 2 filhos (as) | 86 (16.7%) |
| 3 filhos (as) | 45 (8.7%) |
| Mais que 3 filhos (as) | 7 (1.4%) |
| plano_saude | |
| N-Miss | 12 |
| Não | 192 (38.1%) |
| Sim | 312 (61.9%) |
| como_vive | |
| N-Miss | 5 |
| Vivo com a minha família (por exemplo, com companheiro, filhos) | 297 (58.1%) |
| Vivo com amigos/colegas de quarto | 10 (2.0%) |
| Vivo com os meus país | 159 (31.1%) |
| Vivo sozinho | 45 (8.8%) |
| fica_em_casa | |
| Entre 1 a 4 horas por dia | 3 (0.6%) |
| Entre 5 a 8 horas por dia | 21 (4.1%) |
| Entre 9 a 12 horas por dia | 32 (6.2%) |
| Entre 13 a 16 horas por dia | 35 (6.8%) |
| Entre 17 a 20 horas por dia | 58 (11.2%) |
| Mais de 20 horas por dia ou o dia todo | 367 (71.1%) |
| tem_saido_de_casa | |
| Nenhuma | 259 (50.2%) |
| 1 vez por dia | 210 (40.7%) |
| 2 vezes por dia | 25 (4.8%) |
| 3 vezes por dia | 7 (1.4%) |
| Mais de 3 vezes por dia | 15 (2.9%) |
| falou_familia_pessoalmente | |
| Nenhuma | 80 (15.5%) |
| 1 pessoa | 119 (23.1%) |
| 2 pessoas | 145 (28.1%) |
| Mais que 2 pessoas | 172 (33.3%) |
| falou_familia_telefone | |
| 1 pessoa | 126 (24.4%) |
| 2 pessoas | 82 (15.9%) |
| Mais que 2 pessoas | 118 (22.9%) |
| Nenhuma | 190 (36.8%) |
| falou_amigos_telefone | |
| 1 pessoa | 70 (13.6%) |
| 2 pessoas | 71 (13.8%) |
| Mais que 2 pessoas | 207 (40.1%) |
| Nenhuma | 168 (32.6%) |
| testou_coronavirus | |
| Não | 486 (94.2%) |
| Sim | 30 (5.8%) |
| conhece_alguem_coronavirus | |
| N-Miss | 1 |
| Não | 421 (81.7%) |
| Sim | 94 (18.3%) |
| segue_recomendacoes | |
| Não | 13 (2.5%) |
| Sim | 503 (97.5%) |
| concordo_covid_preocupacao_saude | |
| Eu discordo fortemente | 4 (0.8%) |
| Eu discordo | 14 (2.7%) |
| Estou neutro | 11 (2.1%) |
| Eu concordo | 106 (20.5%) |
| Eu concordo fortemente | 381 (73.8%) |
| concordo_quarentena | |
| Eu discordo fortemente | 9 (1.7%) |
| Eu discordo | 33 (6.4%) |
| Estou neutro | 29 (5.6%) |
| Eu concordo | 132 (25.6%) |
| Eu concordo fortemente | 313 (60.7%) |
| area_profissao | |
| N-Miss | 76 |
| Comércio (vendedor, barbeiro, etc) | 51 (11.6%) |
| Construção (pedreiro, etc) | 3 (0.7%) |
| Educação (professor, etc) | 66 (15.0%) |
| Indústria (operador de máquina, etc) | 8 (1.8%) |
| Outros | 172 (39.1%) |
| Saúde (médicos, fisioterapeutas, etc) | 98 (22.3%) |
| Serviço público | 38 (8.6%) |
| Turismo e hotelaria | 4 (0.9%) |
| preocupacao_pagar_conta | |
| N-Miss | 4 |
| Não | 228 (44.5%) |
| Sim | 284 (55.5%) |
| agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos | |
| Não, nenhum | 264 (51.2%) |
| Pouco (1 vez por semana) | 108 (20.9%) |
| Moderado (2 a 3 vezes por semana) | 100 (19.4%) |
| Muito (mais que 3 vezes por semana) | 44 (8.5%) |
| rio_sp | |
| N-Miss | 5 |
| Outros | 103 (20.2%) |
| RJ | 309 (60.5%) |
| SP | 99 (19.4%) |
| total_ces_t1 | |
| Mean (SD) | 23.186 (12.856) |
| Range | 0.000 - 57.000 |
| total_gad_t1 | |
| Mean (SD) | 9.558 (5.754) |
| Range | 0.000 - 21.000 |
NA1.6 modelos de regressão (Tabela 2 e 3)
1.6.1 Modelo irrestrito - depressão
reg_tudo_depre <- lm(total_ces_t1 ~ sexo + idade + escolaridade + estado_civil + quantos_filhos + plano_saude + como_vive + fica_em_casa + tem_saido_de_casa + falou_familia_pessoalmente + falou_familia_telefone + falou_amigos_telefone + testou_coronavirus + conhece_alguem_coronavirus + segue_recomendacoes + concordo_covid_preocupacao_saude + concordo_quarentena + area_profissao + preocupacao_pagar_conta + agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos + rio_sp, ds)1.6.1.1 VIF
car::vif(reg_tudo_depre) GVIF Df GVIF^(1/(2*Df))
sexo 1.402004 1 1.184063
idade 2.494342 1 1.579349
escolaridade 6.109971 6 1.162795
estado_civil 7.242210 4 1.280808
quantos_filhos 5.031859 4 1.223816
plano_saude 1.515407 1 1.231019
como_vive 5.312178 3 1.320927
fica_em_casa 4.168341 5 1.153443
tem_saido_de_casa 3.630577 4 1.174889
falou_familia_pessoalmente 2.341015 3 1.152305
falou_familia_telefone 2.273736 3 1.146718
falou_amigos_telefone 1.692486 3 1.091660
testou_coronavirus 1.383214 1 1.176101
conhece_alguem_coronavirus 1.475888 1 1.214861
segue_recomendacoes 1.235426 1 1.111497
concordo_covid_preocupacao_saude 5.174708 4 1.228106
concordo_quarentena 5.951119 4 1.249755
area_profissao 5.189227 7 1.124809
preocupacao_pagar_conta 1.299222 1 1.139834
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos 1.882504 3 1.111193
rio_sp 1.787295 2 1.1562431.6.1.2 Estimates (Betas)
olsrr::ols_regress(reg_tudo_depre) Model Summary
---------------------------------------------------------------
R 0.554 RMSE 11.251
R-Squared 0.307 Coef. Var 50.824
Adj. R-Squared 0.187 MSE 126.592
Pred R-Squared 0.035 MAE 8.182
---------------------------------------------------------------
RMSE: Root Mean Square Error
MSE: Mean Square Error
MAE: Mean Absolute Error
ANOVA
---------------------------------------------------------------------
Sum of
Squares DF Mean Square F Sig.
---------------------------------------------------------------------
Regression 20105.988 62 324.290 2.562 0.0000
Residual 45320.021 358 126.592
Total 65426.010 420
---------------------------------------------------------------------
Parameter Estimates
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
model Beta Std. Error Std. Beta t Sig lower upper
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
(Intercept) 39.284 14.405 2.727 0.007 10.955 67.614
sexoMasculino -4.465 1.431 -0.162 -3.119 0.002 -7.279 -1.650
idade -0.102 0.065 -0.109 -1.567 0.118 -0.230 0.026
escolaridadeEnsino fundamental completo 1.497 6.330 0.019 0.236 0.813 -10.952 13.946
escolaridadeEnsino médio incompleto -5.092 6.374 -0.065 -0.799 0.425 -17.627 7.443
escolaridadeEnsino médio completo 3.807 5.475 0.107 0.695 0.487 -6.960 14.574
escolaridadeEnsino superior incompleto 1.372 5.476 0.042 0.251 0.802 -9.397 12.140
escolaridadeEnsino superior completo 1.738 5.492 0.058 0.316 0.752 -9.062 12.538
escolaridadePós-graduação 0.912 5.545 0.036 0.164 0.869 -9.993 11.816
estado_civilSeparado(a) ou divorciado(a) -0.980 2.480 -0.022 -0.395 0.693 -5.858 3.898
estado_civilSolteiro(a) -0.617 2.170 -0.024 -0.284 0.776 -4.885 3.651
estado_civilUnião estável (mesmo que sem reconhecimento legal) 2.286 2.246 0.052 1.018 0.309 -2.131 6.703
estado_civilViuvo(a) 8.336 5.708 0.079 1.460 0.145 -2.890 19.561
quantos_filhos1 filho (a) -2.528 1.859 -0.081 -1.360 0.175 -6.184 1.128
quantos_filhos2 filhos (as) -2.030 2.238 -0.063 -0.907 0.365 -6.431 2.371
quantos_filhos3 filhos (as) -2.071 2.672 -0.048 -0.775 0.439 -7.325 3.184
quantos_filhosMais que 3 filhos (as) -1.977 5.064 -0.020 -0.390 0.696 -11.936 7.982
plano_saudeSim 0.185 1.402 0.007 0.132 0.895 -2.572 2.941
como_viveVivo com amigos/colegas de quarto -4.307 4.533 -0.047 -0.950 0.343 -13.221 4.607
como_viveVivo com os meus país 1.938 2.053 0.069 0.944 0.346 -2.099 5.975
como_viveVivo sozinho 0.578 2.794 0.013 0.207 0.836 -4.917 6.073
fica_em_casaEntre 5 a 8 horas por dia -4.966 7.835 -0.083 -0.634 0.527 -20.375 10.442
fica_em_casaEntre 9 a 12 horas por dia -0.687 7.728 -0.014 -0.089 0.929 -15.885 14.510
fica_em_casaEntre 13 a 16 horas por dia 1.351 7.709 0.030 0.175 0.861 -13.809 16.511
fica_em_casaEntre 17 a 20 horas por dia -2.224 7.515 -0.057 -0.296 0.767 -17.002 12.554
fica_em_casaMais de 20 horas por dia ou o dia todo -1.551 7.519 -0.058 -0.206 0.837 -16.337 13.235
tem_saido_de_casa1 vez por dia -2.099 1.507 -0.083 -1.393 0.164 -5.062 0.864
tem_saido_de_casa2 vezes por dia -8.069 2.951 -0.150 -2.734 0.007 -13.872 -2.265
tem_saido_de_casa3 vezes por dia -6.477 5.264 -0.062 -1.230 0.219 -16.829 3.876
tem_saido_de_casaMais de 3 vezes por dia -10.025 3.812 -0.144 -2.630 0.009 -17.522 -2.528
falou_familia_pessoalmente1 pessoa -0.486 2.124 -0.016 -0.229 0.819 -4.662 3.691
falou_familia_pessoalmente2 pessoas 1.171 2.071 0.043 0.565 0.572 -2.902 5.244
falou_familia_pessoalmenteMais que 2 pessoas -1.732 2.021 -0.065 -0.857 0.392 -5.706 2.243
falou_familia_telefone2 pessoas -0.621 1.961 -0.018 -0.317 0.752 -4.478 3.235
falou_familia_telefoneMais que 2 pessoas -2.313 1.866 -0.080 -1.239 0.216 -5.984 1.357
falou_familia_telefoneNenhuma 0.646 1.643 0.025 0.393 0.694 -2.585 3.877
falou_amigos_telefone2 pessoas -3.069 2.179 -0.088 -1.409 0.160 -7.354 1.216
falou_amigos_telefoneMais que 2 pessoas -1.803 1.870 -0.071 -0.964 0.336 -5.480 1.875
falou_amigos_telefoneNenhuma -2.379 1.931 -0.088 -1.232 0.219 -6.177 1.419
testou_coronavirusSim -4.788 2.679 -0.092 -1.787 0.075 -10.057 0.481
conhece_alguem_coronavirusSim 5.327 1.674 0.170 3.181 0.002 2.034 8.620
segue_recomendacoesSim -0.177 3.821 -0.002 -0.046 0.963 -7.691 7.337
concordo_covid_preocupacao_saudeEu discordo -13.065 8.801 -0.167 -1.485 0.139 -30.373 4.243
concordo_covid_preocupacao_saudeEstou neutro -11.228 8.459 -0.130 -1.327 0.185 -27.864 5.407
concordo_covid_preocupacao_saudeEu concordo -15.220 7.512 -0.510 -2.026 0.043 -29.993 -0.448
concordo_covid_preocupacao_saudeEu concordo fortemente -10.491 7.389 -0.378 -1.420 0.157 -25.022 4.040
concordo_quarentenaEu discordo 5.275 5.545 0.112 0.951 0.342 -5.630 16.180
concordo_quarentenaEstou neutro 9.005 5.828 0.174 1.545 0.123 -2.457 20.467
concordo_quarentenaEu concordo 6.622 5.500 0.233 1.204 0.229 -4.195 17.438
concordo_quarentenaEu concordo fortemente 5.943 5.515 0.235 1.078 0.282 -4.903 16.788
area_profissaoConstrução (pedreiro, etc) -12.442 7.255 -0.084 -1.715 0.087 -26.709 1.825
area_profissaoEducação (professor, etc) -2.683 2.482 -0.077 -1.081 0.280 -7.564 2.198
area_profissaoIndústria (operador de máquina, etc) -7.030 4.630 -0.077 -1.519 0.130 -16.135 2.074
area_profissaoOutros -2.762 2.066 -0.108 -1.337 0.182 -6.825 1.301
area_profissaoSaúde (médicos, fisioterapeutas, etc) -6.101 2.385 -0.205 -2.558 0.011 -10.791 -1.411
area_profissaoServiço público -2.060 2.861 -0.046 -0.720 0.472 -7.686 3.567
area_profissaoTurismo e hotelaria 5.321 7.076 0.041 0.752 0.453 -8.594 19.237
preocupacao_pagar_contaSim 4.389 1.259 0.175 3.485 0.001 1.912 6.866
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosPouco (1 vez por semana) -4.543 1.549 -0.151 -2.934 0.004 -7.588 -1.498
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosModerado (2 a 3 vezes por semana) -2.095 1.589 -0.069 -1.318 0.188 -5.221 1.031
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosMuito (mais que 3 vezes por semana) -4.377 2.424 -0.093 -1.805 0.072 -9.144 0.391
rio_spRJ 0.360 1.625 0.014 0.221 0.825 -2.836 3.555
rio_spSP 1.185 2.032 0.038 0.583 0.560 -2.810 5.181
------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------1.6.1.3 Estimates (SSQ)
apaTables::apa.aov.table(reg_tudo_depre)
ANOVA results using total_ces_t1 as the dependent variable
Note: Values in square brackets indicate the bounds of the 90% confidence interval for partial eta-squared 1.6.2 Reg teste com tudo - Ansiedade
reg_tudo_ansi <- lm(total_gad_t1 ~ sexo + idade + escolaridade + estado_civil + quantos_filhos + plano_saude + como_vive + fica_em_casa + tem_saido_de_casa + falou_familia_pessoalmente + falou_familia_telefone + falou_amigos_telefone + testou_coronavirus + conhece_alguem_coronavirus + segue_recomendacoes + concordo_covid_preocupacao_saude + concordo_quarentena + area_profissao + preocupacao_pagar_conta + agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos + rio_sp, ds )
1.6.3 VIF
car::vif(reg_tudo_ansi) GVIF Df GVIF^(1/(2*Df))
sexo 1.402004 1 1.184063
idade 2.494342 1 1.579349
escolaridade 6.109971 6 1.162795
estado_civil 7.242210 4 1.280808
quantos_filhos 5.031859 4 1.223816
plano_saude 1.515407 1 1.231019
como_vive 5.312178 3 1.320927
fica_em_casa 4.168341 5 1.153443
tem_saido_de_casa 3.630577 4 1.174889
falou_familia_pessoalmente 2.341015 3 1.152305
falou_familia_telefone 2.273736 3 1.146718
falou_amigos_telefone 1.692486 3 1.091660
testou_coronavirus 1.383214 1 1.176101
conhece_alguem_coronavirus 1.475888 1 1.214861
segue_recomendacoes 1.235426 1 1.111497
concordo_covid_preocupacao_saude 5.174708 4 1.228106
concordo_quarentena 5.951119 4 1.249755
area_profissao 5.189227 7 1.124809
preocupacao_pagar_conta 1.299222 1 1.139834
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos 1.882504 3 1.111193
rio_sp 1.787295 2 1.1562431.6.3.1 Estimates (beta)
olsrr::ols_regress(reg_tudo_ansi) Model Summary
--------------------------------------------------------------
R 0.541 RMSE 5.262
R-Squared 0.293 Coef. Var 56.559
Adj. R-Squared 0.170 MSE 27.691
Pred R-Squared 0.012 MAE 3.948
--------------------------------------------------------------
RMSE: Root Mean Square Error
MSE: Mean Square Error
MAE: Mean Absolute Error
ANOVA
---------------------------------------------------------------------
Sum of
Squares DF Mean Square F Sig.
---------------------------------------------------------------------
Regression 4101.580 62 66.155 2.389 0.0000
Residual 9913.503 358 27.691
Total 14015.083 420
---------------------------------------------------------------------
Parameter Estimates
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
model Beta Std. Error Std. Beta t Sig lower upper
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
(Intercept) 16.776 6.737 2.490 0.013 3.527 30.026
sexoMasculino -1.849 0.669 -0.145 -2.763 0.006 -3.166 -0.533
idade -0.095 0.030 -0.219 -3.126 0.002 -0.155 -0.035
escolaridadeEnsino fundamental completo 1.254 2.961 0.035 0.423 0.672 -4.569 7.076
escolaridadeEnsino médio incompleto -0.968 2.981 -0.027 -0.325 0.745 -6.831 4.894
escolaridadeEnsino médio completo 3.476 2.561 0.212 1.358 0.175 -1.559 8.512
escolaridadeEnsino superior incompleto 2.994 2.561 0.197 1.169 0.243 -2.043 8.030
escolaridadeEnsino superior completo 2.025 2.568 0.147 0.788 0.431 -3.026 7.076
escolaridadePós-graduação 2.032 2.593 0.172 0.783 0.434 -3.068 7.131
estado_civilSeparado(a) ou divorciado(a) -0.667 1.160 -0.032 -0.575 0.565 -2.949 1.614
estado_civilSolteiro(a) -0.671 1.015 -0.057 -0.661 0.509 -2.667 1.325
estado_civilUnião estável (mesmo que sem reconhecimento legal) 0.690 1.050 0.034 0.657 0.512 -1.376 2.756
estado_civilViuvo(a) -0.277 2.670 -0.006 -0.104 0.917 -5.528 4.973
quantos_filhos1 filho (a) -0.172 0.869 -0.012 -0.198 0.843 -1.882 1.538
quantos_filhos2 filhos (as) -0.191 1.047 -0.013 -0.182 0.856 -2.249 1.868
quantos_filhos3 filhos (as) 0.378 1.250 0.019 0.303 0.762 -2.079 2.836
quantos_filhosMais que 3 filhos (as) 1.810 2.368 0.040 0.764 0.445 -2.848 6.468
plano_saudeSim 0.556 0.655 0.046 0.849 0.397 -0.733 1.845
como_viveVivo com amigos/colegas de quarto -0.771 2.120 -0.018 -0.364 0.716 -4.940 3.398
como_viveVivo com os meus país 0.335 0.960 0.026 0.349 0.727 -1.553 2.223
como_viveVivo sozinho -0.622 1.307 -0.030 -0.476 0.634 -3.192 1.948
fica_em_casaEntre 5 a 8 horas por dia -5.165 3.664 -0.186 -1.409 0.160 -12.371 2.042
fica_em_casaEntre 9 a 12 horas por dia -1.376 3.614 -0.061 -0.381 0.704 -8.483 5.732
fica_em_casaEntre 13 a 16 horas por dia -1.408 3.605 -0.066 -0.391 0.696 -8.499 5.682
fica_em_casaEntre 17 a 20 horas por dia -2.222 3.515 -0.124 -0.632 0.528 -9.134 4.690
fica_em_casaMais de 20 horas por dia ou o dia todo -2.324 3.516 -0.188 -0.661 0.509 -9.240 4.591
tem_saido_de_casa1 vez por dia -0.777 0.705 -0.066 -1.103 0.271 -2.163 0.608
tem_saido_de_casa2 vezes por dia -2.601 1.380 -0.105 -1.884 0.060 -5.315 0.113
tem_saido_de_casa3 vezes por dia -2.202 2.462 -0.045 -0.895 0.372 -7.044 2.639
tem_saido_de_casaMais de 3 vezes por dia -2.152 1.783 -0.067 -1.207 0.228 -5.659 1.354
falou_familia_pessoalmente1 pessoa 0.313 0.993 0.022 0.316 0.753 -1.640 2.267
falou_familia_pessoalmente2 pessoas -0.078 0.969 -0.006 -0.080 0.936 -1.983 1.827
falou_familia_pessoalmenteMais que 2 pessoas -0.226 0.945 -0.018 -0.239 0.811 -2.085 1.633
falou_familia_telefone2 pessoas 0.385 0.917 0.025 0.419 0.675 -1.419 2.188
falou_familia_telefoneMais que 2 pessoas 0.167 0.873 0.013 0.191 0.848 -1.550 1.884
falou_familia_telefoneNenhuma 0.645 0.768 0.053 0.839 0.402 -0.866 2.156
falou_amigos_telefone2 pessoas -0.991 1.019 -0.062 -0.972 0.332 -2.995 1.013
falou_amigos_telefoneMais que 2 pessoas -0.523 0.875 -0.044 -0.598 0.551 -2.242 1.197
falou_amigos_telefoneNenhuma -0.746 0.903 -0.060 -0.826 0.410 -2.522 1.031
testou_coronavirusSim 0.036 1.253 0.001 0.029 0.977 -2.428 2.500
conhece_alguem_coronavirusSim 1.734 0.783 0.120 2.214 0.027 0.194 3.274
segue_recomendacoesSim -0.467 1.787 -0.013 -0.261 0.794 -3.981 3.048
concordo_covid_preocupacao_saudeEu discordo -4.076 4.116 -0.113 -0.990 0.323 -12.170 4.019
concordo_covid_preocupacao_saudeEstou neutro -1.561 3.956 -0.039 -0.394 0.693 -9.341 6.220
concordo_covid_preocupacao_saudeEu concordo -3.578 3.513 -0.259 -1.018 0.309 -10.487 3.331
concordo_covid_preocupacao_saudeEu concordo fortemente -1.870 3.456 -0.146 -0.541 0.589 -8.667 4.926
concordo_quarentenaEu discordo 0.182 2.593 0.008 0.070 0.944 -4.919 5.282
concordo_quarentenaEstou neutro 2.767 2.726 0.115 1.015 0.311 -2.594 8.128
concordo_quarentenaEu concordo 0.813 2.572 0.062 0.316 0.752 -4.246 5.872
concordo_quarentenaEu concordo fortemente 1.689 2.579 0.144 0.655 0.513 -3.383 6.761
area_profissaoConstrução (pedreiro, etc) -5.084 3.393 -0.074 -1.498 0.135 -11.756 1.589
area_profissaoEducação (professor, etc) -1.024 1.161 -0.064 -0.882 0.378 -3.308 1.259
area_profissaoIndústria (operador de máquina, etc) -3.481 2.165 -0.082 -1.608 0.109 -7.739 0.777
area_profissaoOutros -2.094 0.966 -0.177 -2.167 0.031 -3.994 -0.194
area_profissaoSaúde (médicos, fisioterapeutas, etc) -2.329 1.115 -0.169 -2.088 0.037 -4.522 -0.136
area_profissaoServiço público -1.814 1.338 -0.088 -1.356 0.176 -4.446 0.818
area_profissaoTurismo e hotelaria 6.367 3.309 0.107 1.924 0.055 -0.141 12.875
preocupacao_pagar_contaSim 1.695 0.589 0.146 2.878 0.004 0.537 2.854
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosPouco (1 vez por semana) -1.437 0.724 -0.103 -1.984 0.048 -2.861 -0.013
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosModerado (2 a 3 vezes por semana) -0.497 0.743 -0.035 -0.669 0.504 -1.959 0.965
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosMuito (mais que 3 vezes por semana) -3.410 1.134 -0.157 -3.008 0.003 -5.640 -1.180
rio_spRJ -0.302 0.760 -0.026 -0.397 0.691 -1.797 1.192
rio_spSP 0.807 0.950 0.056 0.849 0.396 -1.062 2.676
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------1.6.3.2 Estimates (SSQ)
apaTables::apa.reg.table(reg_tudo_ansi)
Regression results using total_gad_t1 as the criterion
Note. A significant b-weight indicates the semi-partial correlation is also significant.
b represents unstandardized regression weights.
sr2 represents the semi-partial correlation squared.
Square brackets are used to enclose the lower and upper limits of a confidence interval.
* indicates p < .05. ** indicates p < .01.
summary(reg_tudo_ansi)
Call:
lm(formula = total_gad_t1 ~ sexo + idade + escolaridade + estado_civil +
quantos_filhos + plano_saude + como_vive + fica_em_casa +
tem_saido_de_casa + falou_familia_pessoalmente + falou_familia_telefone +
falou_amigos_telefone + testou_coronavirus + conhece_alguem_coronavirus +
segue_recomendacoes + concordo_covid_preocupacao_saude +
concordo_quarentena + area_profissao + preocupacao_pagar_conta +
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos +
rio_sp, data = ds)
Residuals:
Min 1Q Median 3Q Max
-10.3174 -3.7066 -0.2752 3.3496 14.2821
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 16.77627 6.73734 2.490 0.01323 *
sexoMasculino -1.84942 0.66938 -2.763 0.00603 **
idade -0.09516 0.03045 -3.126 0.00192 **
escolaridadeEnsino fundamental completo 1.25380 2.96062 0.423 0.67219
escolaridadeEnsino médio incompleto -0.96837 2.98112 -0.325 0.74550
escolaridadeEnsino médio completo 3.47623 2.56057 1.358 0.17545
escolaridadeEnsino superior incompleto 2.99355 2.56099 1.169 0.24322
escolaridadeEnsino superior completo 2.02510 2.56848 0.788 0.43096
escolaridadePós-graduação 2.03154 2.59327 0.783 0.43392
estado_civilSeparado(a) ou divorciado(a) -0.66744 1.16013 -0.575 0.56544
estado_civilSolteiro(a) -0.67063 1.01501 -0.661 0.50922
estado_civilUnião estável (mesmo que sem reconhecimento legal) 0.68992 1.05047 0.657 0.51175
estado_civilViuvo(a) -0.27733 2.66967 -0.104 0.91732
quantos_filhos1 filho (a) -0.17219 0.86949 -0.198 0.84313
quantos_filhos2 filhos (as) -0.19067 1.04668 -0.182 0.85556
quantos_filhos3 filhos (as) 0.37812 1.24955 0.303 0.76237
quantos_filhosMais que 3 filhos (as) 1.80999 2.36837 0.764 0.44523
plano_saudeSim 0.55619 0.65548 0.849 0.39672
como_viveVivo com amigos/colegas de quarto -0.77105 2.11990 -0.364 0.71628
como_viveVivo com os meus país 0.33484 0.96004 0.349 0.72746
como_viveVivo sozinho -0.62239 1.30681 -0.476 0.63417
fica_em_casaEntre 5 a 8 horas por dia -5.16486 3.66444 -1.409 0.15957
fica_em_casaEntre 9 a 12 horas por dia -1.37567 3.61420 -0.381 0.70371
fica_em_casaEntre 13 a 16 horas por dia -1.40812 3.60540 -0.391 0.69636
fica_em_casaEntre 17 a 20 horas por dia -2.22211 3.51456 -0.632 0.52762
fica_em_casaMais de 20 horas por dia ou o dia todo -2.32447 3.51644 -0.661 0.50902
tem_saido_de_casa1 vez por dia -0.77722 0.70461 -1.103 0.27075
tem_saido_de_casa2 vezes por dia -2.60102 1.38027 -1.884 0.06032 .
tem_saido_de_casa3 vezes por dia -2.20239 2.46194 -0.895 0.37162
tem_saido_de_casaMais de 3 vezes por dia -2.15227 1.78298 -1.207 0.22818
falou_familia_pessoalmente1 pessoa 0.31340 0.99325 0.316 0.75254
falou_familia_pessoalmente2 pessoas -0.07767 0.96873 -0.080 0.93614
falou_familia_pessoalmenteMais que 2 pessoas -0.22587 0.94512 -0.239 0.81125
falou_familia_telefone2 pessoas 0.38463 0.91720 0.419 0.67521
falou_familia_telefoneMais que 2 pessoas 0.16689 0.87294 0.191 0.84849
falou_familia_telefoneNenhuma 0.64507 0.76840 0.839 0.40176
falou_amigos_telefone2 pessoas -0.99086 1.01902 -0.972 0.33152
falou_amigos_telefoneMais que 2 pessoas -0.52256 0.87457 -0.598 0.55055
falou_amigos_telefoneNenhuma -0.74568 0.90330 -0.826 0.40963
testou_coronavirusSim 0.03583 1.25306 0.029 0.97720
conhece_alguem_coronavirusSim 1.73383 0.78315 2.214 0.02746 *
segue_recomendacoesSim -0.46665 1.78704 -0.261 0.79414
concordo_covid_preocupacao_saudeEu discordo -4.07556 4.11618 -0.990 0.32278
concordo_covid_preocupacao_saudeEstou neutro -1.56073 3.95627 -0.394 0.69345
concordo_covid_preocupacao_saudeEu concordo -3.57782 3.51327 -1.018 0.30919
concordo_covid_preocupacao_saudeEu concordo fortemente -1.87043 3.45575 -0.541 0.58867
concordo_quarentenaEu discordo 0.18153 2.59342 0.070 0.94424
concordo_quarentenaEstou neutro 2.76716 2.72586 1.015 0.31072
concordo_quarentenaEu concordo 0.81266 2.57238 0.316 0.75225
concordo_quarentenaEu concordo fortemente 1.68902 2.57930 0.655 0.51299
area_profissaoConstrução (pedreiro, etc) -5.08374 3.39302 -1.498 0.13494
area_profissaoEducação (professor, etc) -1.02446 1.16090 -0.882 0.37812
area_profissaoIndústria (operador de máquina, etc) -3.48080 2.16525 -1.608 0.10881
area_profissaoOutros -2.09404 0.96626 -2.167 0.03088 *
area_profissaoSaúde (médicos, fisioterapeutas, etc) -2.32891 1.11530 -2.088 0.03749 *
area_profissaoServiço público -1.81397 1.33816 -1.356 0.17609
area_profissaoTurismo e hotelaria 6.36724 3.30930 1.924 0.05514 .
preocupacao_pagar_contaSim 1.69539 0.58900 2.878 0.00424 **
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosPouco (1 vez por semana) -1.43715 0.72425 -1.984 0.04798 *
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosModerado (2 a 3 vezes por semana) -0.49716 0.74340 -0.669 0.50407
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicosMuito (mais que 3 vezes por semana) -3.41018 1.13380 -3.008 0.00282 **
rio_spRJ -0.30203 0.75995 -0.397 0.69128
rio_spSP 0.80679 0.95028 0.849 0.39645
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 5.262 on 358 degrees of freedom
(95 observations deleted due to missingness)
Multiple R-squared: 0.2927, Adjusted R-squared: 0.1702
F-statistic: 2.389 on 62 and 358 DF, p-value: 3.331e-071.6.3.3 Post hocs
glht(reg_tudo_depre, linfct = mcp(sexo = 'Tukey')) %>% summary()
Simultaneous Tests for General Linear Hypotheses
Multiple Comparisons of Means: Tukey Contrasts
Fit: lm(formula = total_ces_t1 ~ sexo + idade + escolaridade + estado_civil +
quantos_filhos + plano_saude + como_vive + fica_em_casa +
tem_saido_de_casa + falou_familia_pessoalmente + falou_familia_telefone +
falou_amigos_telefone + testou_coronavirus + conhece_alguem_coronavirus +
segue_recomendacoes + concordo_covid_preocupacao_saude +
concordo_quarentena + area_profissao + preocupacao_pagar_conta +
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos +
rio_sp, data = ds)
Linear Hypotheses:
Estimate Std. Error t value Pr(>|t|)
Masculino - Feminino == 0 -4.465 1.431 -3.119 0.00196 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Adjusted p values reported -- single-step method)TukeyHSD(aov(total_ces_t1 ~ sexo + escolaridade + estado_civil + quantos_filhos + plano_saude + como_vive + fica_em_casa + tem_saido_de_casa + falou_familia_pessoalmente + falou_familia_telefone + falou_amigos_telefone + testou_coronavirus + conhece_alguem_coronavirus + segue_recomendacoes + concordo_covid_preocupacao_saude + concordo_quarentena + area_profissao + preocupacao_pagar_conta + agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos + rio_sp, ds)) Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = total_ces_t1 ~ sexo + escolaridade + estado_civil + quantos_filhos + plano_saude + como_vive + fica_em_casa + tem_saido_de_casa + falou_familia_pessoalmente + falou_familia_telefone + falou_amigos_telefone + testou_coronavirus + conhece_alguem_coronavirus + segue_recomendacoes + concordo_covid_preocupacao_saude + concordo_quarentena + area_profissao + preocupacao_pagar_conta + agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos + rio_sp, data = ds)
$sexo
diff lwr upr p adj
Masculino-Feminino -4.591178 -6.973076 -2.209281 0.0001762
$escolaridade
diff lwr upr p adj
Ensino fundamental completo-Ensino fundamental incompleto -0.5364973 -17.502576 16.4295813 0.9999999
Ensino médio incompleto-Ensino fundamental incompleto -5.6274064 -22.593485 11.3386722 0.9573517
Ensino médio completo-Ensino fundamental incompleto 1.7078355 -12.595122 16.0107928 0.9998425
Ensino superior incompleto-Ensino fundamental incompleto 1.5410179 -12.648985 15.7310211 0.9999096
Ensino superior completo-Ensino fundamental incompleto -2.2625850 -16.334469 11.8092988 0.9991207
Pós-graduação-Ensino fundamental incompleto -3.8148591 -17.711278 10.0815601 0.9834633
Ensino médio incompleto-Ensino fundamental completo -5.0909091 -19.345266 9.1634481 0.9394510
Ensino médio completo-Ensino fundamental completo 2.2443328 -8.706168 13.1948333 0.9965558
Ensino superior incompleto-Ensino fundamental completo 2.0775152 -8.725034 12.8800642 0.9975818
Ensino superior completo-Ensino fundamental completo -1.7260877 -12.373002 8.9208268 0.9990783
Pós-graduação-Ensino fundamental completo -3.2783618 -13.692263 7.1355396 0.9669439
Ensino médio completo-Ensino médio incompleto 7.3352419 -3.615259 18.2857424 0.4252449
Ensino superior incompleto-Ensino médio incompleto 7.1684243 -3.634125 17.9709733 0.4372161
Ensino superior completo-Ensino médio incompleto 3.3648213 -7.282093 14.0117359 0.9663021
Pós-graduação-Ensino médio incompleto 1.8125473 -8.601354 12.2264487 0.9986203
Ensino superior incompleto-Ensino médio completo -0.1668176 -5.947979 5.6143438 1.0000000
Ensino superior completo-Ensino médio completo -3.9704205 -9.455268 1.5144270 0.3278359
Pós-graduação-Ensino médio completo -5.5226946 -10.540274 -0.5051150 0.0204378
Ensino superior completo-Ensino superior incompleto -3.8036030 -8.986764 1.3795577 0.3112430
Pós-graduação-Ensino superior incompleto -5.3558770 -10.041782 -0.6699719 0.0136153
Pós-graduação-Ensino superior completo -1.5522741 -5.867298 2.7627501 0.9373861
$estado_civil
diff lwr upr p adj
Separado(a) ou divorciado(a)-Casado(a) -3.1958553 -8.9391493 2.547439 0.5464691
Solteiro(a)-Casado(a) 2.3575138 -0.9811685 5.696196 0.3001972
União estável (mesmo que sem reconhecimento legal)-Casado(a) 2.2418554 -3.3715075 7.855218 0.8091018
Viuvo(a)-Casado(a) 3.0117761 -9.8306198 15.854172 0.9678629
Solteiro(a)-Separado(a) ou divorciado(a) 5.5533691 -0.1700862 11.276824 0.0620709
União estável (mesmo que sem reconhecimento legal)-Separado(a) ou divorciado(a) 5.4377107 -1.8507037 12.726125 0.2465538
Viuvo(a)-Separado(a) ou divorciado(a) 6.2076314 -7.4502701 19.865533 0.7242463
União estável (mesmo que sem reconhecimento legal)-Solteiro(a) -0.1156584 -5.7087217 5.477405 0.9999976
Viuvo(a)-Solteiro(a) 0.6542623 -12.1792737 13.487798 0.9999138
Viuvo(a)-União estável (mesmo que sem reconhecimento legal) 0.7699207 -12.8338541 14.373695 0.9998694
$quantos_filhos
diff lwr upr p adj
1 filho (a)-Nenhum -1.2752837 -5.252493 2.701925 0.9044761
2 filhos (as)-Nenhum -1.9659930 -6.071138 2.139152 0.6832983
3 filhos (as)-Nenhum -2.0956883 -7.535087 3.343711 0.8286432
Mais que 3 filhos (as)-Nenhum -5.1930344 -17.064653 6.678584 0.7518072
2 filhos (as)-1 filho (a) -0.6907093 -5.567450 4.186031 0.9951737
3 filhos (as)-1 filho (a) -0.8204047 -6.863371 5.222562 0.9959023
Mais que 3 filhos (as)-1 filho (a) -3.9177507 -16.077755 8.242254 0.9029532
3 filhos (as)-2 filhos (as) -0.1296954 -6.257621 5.998230 0.9999974
Mais que 3 filhos (as)-2 filhos (as) -3.2270414 -15.429489 8.975406 0.9506203
Mais que 3 filhos (as)-3 filhos (as) -3.0973461 -15.810905 9.616213 0.9631273
$plano_saude
diff lwr upr p adj
Sim-Não -0.07774821 -2.32123 2.165734 0.9457021
$como_vive
diff lwr upr p adj
Vivo com amigos/colegas de quarto-Vivo com a minha família (por exemplo, com companheiro, filhos) -3.1299185 -13.574698 7.314861 0.8664204
Vivo com os meus país-Vivo com a minha família (por exemplo, com companheiro, filhos) 0.7120391 -2.535167 3.959245 0.9420959
Vivo sozinho-Vivo com a minha família (por exemplo, com companheiro, filhos) -1.3609325 -6.534559 3.812695 0.9050476
Vivo com os meus país-Vivo com amigos/colegas de quarto 3.8419576 -6.795191 14.479107 0.7875475
Vivo sozinho-Vivo com amigos/colegas de quarto 1.7689860 -9.605152 13.143124 0.9780841
Vivo sozinho-Vivo com os meus país -2.0729716 -7.624722 3.478779 0.7701241
$fica_em_casa
diff lwr upr p adj
Entre 5 a 8 horas por dia-Entre 1 a 4 horas por dia -3.2189513 -23.287166 16.849264 0.9974222
Entre 9 a 12 horas por dia-Entre 1 a 4 horas por dia 0.9075703 -18.652506 20.467647 0.9999942
Entre 13 a 16 horas por dia-Entre 1 a 4 horas por dia 1.5545251 -17.900671 21.009721 0.9999136
Entre 17 a 20 horas por dia-Entre 1 a 4 horas por dia 0.7976389 -18.414594 20.009872 0.9999967
Mais de 20 horas por dia ou o dia todo-Entre 1 a 4 horas por dia 2.3059904 -16.441369 21.053350 0.9992825
Entre 9 a 12 horas por dia-Entre 5 a 8 horas por dia 4.1265216 -5.344470 13.597513 0.8125564
Entre 13 a 16 horas por dia-Entre 5 a 8 horas por dia 4.7734763 -4.478969 14.025921 0.6784622
Entre 17 a 20 horas por dia-Entre 5 a 8 horas por dia 4.0165902 -4.713414 12.746594 0.7748332
Mais de 20 horas por dia ou o dia todo-Entre 5 a 8 horas por dia 5.5249417 -2.127939 13.177822 0.3062677
Entre 13 a 16 horas por dia-Entre 9 a 12 horas por dia 0.6469548 -7.444465 8.738375 0.9999133
Entre 17 a 20 horas por dia-Entre 9 a 12 horas por dia -0.1099314 -7.598343 7.378480 1.0000000
Mais de 20 horas por dia ou o dia todo-Entre 9 a 12 horas por dia 1.3984201 -4.800767 7.597607 0.9873418
Entre 17 a 20 horas por dia-Entre 13 a 16 horas por dia -0.7568862 -7.966905 6.453132 0.9996689
Mais de 20 horas por dia ou o dia todo-Entre 13 a 16 horas por dia 0.7514653 -5.108396 6.611327 0.9991215
Mais de 20 horas por dia ou o dia todo-Entre 17 a 20 horas por dia 1.5083515 -3.485966 6.502669 0.9544516
$tem_saido_de_casa
diff lwr upr p adj
1 vez por dia-Nenhuma -0.7770095 -3.967945 2.413926 0.9631926
2 vezes por dia-Nenhuma -5.2656778 -11.941262 1.409907 0.1964259
3 vezes por dia-Nenhuma -7.6797748 -20.485766 5.126217 0.4701512
Mais de 3 vezes por dia-Nenhuma -7.1265451 -15.670499 1.417408 0.1514985
2 vezes por dia-1 vez por dia -4.4886683 -11.214635 2.237299 0.3580466
3 vezes por dia-1 vez por dia -6.9027654 -19.735093 5.929562 0.5795523
Mais de 3 vezes por dia-1 vez por dia -6.3495357 -14.932912 2.233840 0.2546213
3 vezes por dia-2 vezes por dia -2.4140971 -16.522605 11.694410 0.9900448
Mais de 3 vezes por dia-2 vezes por dia -1.8608674 -12.255852 8.534117 0.9881878
Mais de 3 vezes por dia-3 vezes por dia 0.5532297 -14.529399 15.635859 0.9999768
$falou_familia_pessoalmente
diff lwr upr p adj
1 pessoa-Nenhuma -0.6525857 -5.337038 4.031867 0.9840616
2 pessoas-Nenhuma 0.5369240 -3.875293 4.949141 0.9892537
Mais que 2 pessoas-Nenhuma -1.3237074 -5.651607 3.004192 0.8592362
2 pessoas-1 pessoa 1.1895096 -2.827252 5.206271 0.8704313
Mais que 2 pessoas-1 pessoa -0.6711217 -4.595078 3.252835 0.9712183
Mais que 2 pessoas-2 pessoas -1.8606313 -5.455210 1.733948 0.5405626
$falou_familia_telefone
diff lwr upr p adj
2 pessoas-1 pessoa -1.6007400 -6.154889 2.953409 0.8010034
Mais que 2 pessoas-1 pessoa -1.8933820 -5.959413 2.172649 0.6260233
Nenhuma-1 pessoa 0.2447022 -3.527354 4.016759 0.9983281
Mais que 2 pessoas-2 pessoas -0.2926420 -4.802302 4.217018 0.9983266
Nenhuma-2 pessoas 1.8454422 -2.401067 6.091951 0.6763822
Nenhuma-Mais que 2 pessoas 2.1380842 -1.580138 5.856306 0.4481068
$falou_amigos_telefone
diff lwr upr p adj
2 pessoas-1 pessoa -2.3593562 -7.634865 2.916152 0.6560295
Mais que 2 pessoas-1 pessoa -0.8380072 -5.266378 3.590364 0.9616424
Nenhuma-1 pessoa -0.7067320 -5.301691 3.888227 0.9787756
Mais que 2 pessoas-2 pessoas 1.5213490 -2.749685 5.792383 0.7945065
Nenhuma-2 pessoas 1.6526242 -2.790900 6.096149 0.7722536
Nenhuma-Mais que 2 pessoas 0.1312752 -3.263484 3.526034 0.9996435
$testou_coronavirus
diff lwr upr p adj
Sim-Não -0.8927708 -5.38181 3.596268 0.6959464
$conhece_alguem_coronavirus
diff lwr upr p adj
Sim-Não 4.956827 2.240755 7.672898 0.000378
$segue_recomendacoes
diff lwr upr p adj
Sim-Não 0.8363692 -5.937711 7.610449 0.8082931
$concordo_covid_preocupacao_saude
diff lwr upr p adj
Eu discordo-Eu discordo fortemente -18.6355830 -38.7685945 1.4974285 0.0846881
Estou neutro-Eu discordo fortemente -13.7393718 -34.3461660 6.8674225 0.3590211
Eu concordo-Eu discordo fortemente -18.4339648 -36.5595611 -0.3083686 0.0440291
Eu concordo fortemente-Eu discordo fortemente -14.9303592 -32.8644955 3.0037771 0.1529376
Estou neutro-Eu discordo 4.8962112 -8.9968867 18.7893090 0.8701461
Eu concordo-Eu discordo 0.2016181 -9.6429457 10.0461820 0.9999977
Eu concordo fortemente-Eu discordo 3.7052238 -5.7822122 13.1926598 0.8214675
Eu concordo-Estou neutro -4.6945931 -15.4750038 6.0858176 0.7549015
Eu concordo fortemente-Estou neutro -1.1909874 -11.6462856 9.2643107 0.9979297
Eu concordo fortemente-Eu concordo 3.5036057 -0.1310236 7.1382349 0.0649342
$concordo_quarentena
diff lwr upr p adj
Eu discordo-Eu discordo fortemente 4.3291101 -9.422162 18.080382 0.9101333
Estou neutro-Eu discordo fortemente 6.6371452 -7.362415 20.636705 0.6915008
Eu concordo-Eu discordo fortemente 3.8574717 -9.101147 16.816091 0.9255873
Eu concordo fortemente-Eu discordo fortemente 4.1768702 -8.594511 16.948251 0.8980351
Estou neutro-Eu discordo 2.3080351 -5.853165 10.469235 0.9375871
Eu concordo-Eu discordo -0.4716384 -6.679965 5.736688 0.9995799
Eu concordo fortemente-Eu discordo -0.1522398 -5.959621 5.655141 0.9999939
Eu concordo-Estou neutro -2.7796735 -9.520112 3.960765 0.7902233
Eu concordo fortemente-Estou neutro -2.4602750 -8.833333 3.912783 0.8276102
Eu concordo fortemente-Eu concordo 0.3193985 -3.223763 3.862560 0.9991741
$area_profissao
diff lwr upr p adj
Construção (pedreiro, etc)-Comércio (vendedor, barbeiro, etc) -11.8794911 -32.348735 8.5897526 0.6412125
Educação (professor, etc)-Comércio (vendedor, barbeiro, etc) -2.1816348 -8.784751 4.4214812 0.9731825
Indústria (operador de máquina, etc)-Comércio (vendedor, barbeiro, etc) -6.5227232 -19.669333 6.6238866 0.8003288
Outros-Comércio (vendedor, barbeiro, etc) -2.0597547 -7.750816 3.6313066 0.9558281
Saúde (médicos, fisioterapeutas, etc)-Comércio (vendedor, barbeiro, etc) -5.2165584 -11.335972 0.9028554 0.1596253
Serviço público-Comércio (vendedor, barbeiro, etc) -2.5750768 -10.188241 5.0380878 0.9694391
Turismo e hotelaria-Comércio (vendedor, barbeiro, etc) 2.1427408 -15.760536 20.0460171 0.9999588
Educação (professor, etc)-Construção (pedreiro, etc) 9.6978563 -10.607618 30.0033306 0.8299900
Indústria (operador de máquina, etc)-Construção (pedreiro, etc) 5.3567678 -17.914341 28.6278764 0.9969139
Outros-Construção (pedreiro, etc) 9.8197364 -10.207728 29.8472006 0.8099206
Saúde (médicos, fisioterapeutas, etc)-Construção (pedreiro, etc) 6.6629327 -13.490439 26.8163039 0.9730839
Serviço público-Construção (pedreiro, etc) 9.3044143 -11.351612 29.9604405 0.8688621
Turismo e hotelaria-Construção (pedreiro, etc) 14.0222319 -12.231119 40.2755830 0.7328000
Indústria (operador de máquina, etc)-Educação (professor, etc) -4.3410885 -17.231228 8.5490511 0.9701651
Outros-Educação (professor, etc) 0.1218801 -4.948675 5.1924350 1.0000000
Saúde (médicos, fisioterapeutas, etc)-Educação (professor, etc) -3.0349237 -8.581957 2.5121092 0.7078696
Serviço público-Educação (professor, etc) -0.3934420 -7.554631 6.7677466 0.9999998
Turismo e hotelaria-Educação (professor, etc) 4.3243755 -13.391427 22.0401778 0.9955466
Outros-Indústria (operador de máquina, etc) 4.4629685 -7.984631 16.9105681 0.9579985
Saúde (médicos, fisioterapeutas, etc)-Indústria (operador de máquina, etc) 1.3061648 -11.343016 13.9553459 0.9999850
Serviço público-Indústria (operador de máquina, etc) 3.9476464 -9.487934 17.3832273 0.9863259
Turismo e hotelaria-Indústria (operador de máquina, etc) 8.6654640 -12.384046 29.7149737 0.9144170
Saúde (médicos, fisioterapeutas, etc)-Outros -3.1568037 -7.579090 1.2654822 0.3685593
Serviço público-Outros -0.5153221 -6.845382 5.8147381 0.9999971
Turismo e hotelaria-Outros 4.2024955 -13.193961 21.5989517 0.9958254
Serviço público-Saúde (médicos, fisioterapeutas, etc) 2.6414816 -4.076309 9.3592720 0.9319518
Turismo e hotelaria-Saúde (médicos, fisioterapeutas, etc) 7.3592992 -10.181959 24.9005573 0.9061314
Turismo e hotelaria-Serviço público 4.7178176 -13.398716 22.8343510 0.9933665
$preocupacao_pagar_conta
diff lwr upr p adj
Sim-Não 3.703518 1.526326 5.88071 0.0009088
$agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos
diff lwr upr p adj
Pouco (1 vez por semana)-Não, nenhum -3.594220 -7.237759 0.04931934 0.0547484
Moderado (2 a 3 vezes por semana)-Não, nenhum -1.503646 -5.189515 2.18222355 0.7183042
Muito (mais que 3 vezes por semana)-Não, nenhum -3.232326 -8.758026 2.29337492 0.4326238
Moderado (2 a 3 vezes por semana)-Pouco (1 vez por semana) 2.090574 -2.235947 6.41709438 0.5972181
Muito (mais que 3 vezes por semana)-Pouco (1 vez por semana) 0.361894 -5.610251 6.33403887 0.9986355
Muito (mais que 3 vezes por semana)-Moderado (2 a 3 vezes por semana) -1.728680 -7.726744 4.26938397 0.8792286
$rio_sp
diff lwr upr p adj
RJ-Outros 0.1899082 -3.166408 3.546225 0.9902716
SP-Outros 0.7857919 -3.308688 4.880271 0.8936964
SP-RJ 0.5958837 -2.730568 3.922335 0.9067138# regressão de depressão
glht(reg_tudo_depre, linfct = mcp(agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos = 'Tukey')) %>% summary()
Simultaneous Tests for General Linear Hypotheses
Multiple Comparisons of Means: Tukey Contrasts
Fit: lm(formula = total_ces_t1 ~ sexo + idade + escolaridade + estado_civil +
quantos_filhos + plano_saude + como_vive + fica_em_casa +
tem_saido_de_casa + falou_familia_pessoalmente + falou_familia_telefone +
falou_amigos_telefone + testou_coronavirus + conhece_alguem_coronavirus +
segue_recomendacoes + concordo_covid_preocupacao_saude +
concordo_quarentena + area_profissao + preocupacao_pagar_conta +
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos +
rio_sp, data = ds)
Linear Hypotheses:
Estimate Std. Error t value Pr(>|t|)
Pouco (1 vez por semana) - Não, nenhum == 0 -4.5431 1.5485 -2.934 0.0178 *
Moderado (2 a 3 vezes por semana) - Não, nenhum == 0 -2.0948 1.5895 -1.318 0.5416
Muito (mais que 3 vezes por semana) - Não, nenhum == 0 -4.3768 2.4242 -1.805 0.2631
Moderado (2 a 3 vezes por semana) - Pouco (1 vez por semana) == 0 2.4483 1.8282 1.339 0.5280
Muito (mais que 3 vezes por semana) - Pouco (1 vez por semana) == 0 0.1663 2.5960 0.064 0.9999
Muito (mais que 3 vezes por semana) - Moderado (2 a 3 vezes por semana) == 0 -2.2820 2.5590 -0.892 0.8028
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Adjusted p values reported -- single-step method)# regressão de ansiedade
glht(reg_tudo_ansi, linfct = mcp(agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos = 'Tukey')) %>% summary()
Simultaneous Tests for General Linear Hypotheses
Multiple Comparisons of Means: Tukey Contrasts
Fit: lm(formula = total_gad_t1 ~ sexo + idade + escolaridade + estado_civil +
quantos_filhos + plano_saude + como_vive + fica_em_casa +
tem_saido_de_casa + falou_familia_pessoalmente + falou_familia_telefone +
falou_amigos_telefone + testou_coronavirus + conhece_alguem_coronavirus +
segue_recomendacoes + concordo_covid_preocupacao_saude +
concordo_quarentena + area_profissao + preocupacao_pagar_conta +
agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos +
rio_sp, data = ds)
Linear Hypotheses:
Estimate Std. Error t value Pr(>|t|)
Pouco (1 vez por semana) - Não, nenhum == 0 -1.4372 0.7242 -1.984 0.1881
Moderado (2 a 3 vezes por semana) - Não, nenhum == 0 -0.4972 0.7434 -0.669 0.9054
Muito (mais que 3 vezes por semana) - Não, nenhum == 0 -3.4102 1.1338 -3.008 0.0142 *
Moderado (2 a 3 vezes por semana) - Pouco (1 vez por semana) == 0 0.9400 0.8551 1.099 0.6813
Muito (mais que 3 vezes por semana) - Pouco (1 vez por semana) == 0 -1.9730 1.2142 -1.625 0.3553
Muito (mais que 3 vezes por semana) - Moderado (2 a 3 vezes por semana) == 0 -2.9130 1.1968 -2.434 0.0689 .
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Adjusted p values reported -- single-step method)1.6.4 Plots
1.6.4.1 Depressão / ansiedade - Sexo
pairwise.t.test(x = ds$total_ces_t1, g = ds$agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos, p.adjust.method = "bonferroni", ... = ds)
Pairwise comparisons using t tests with pooled SD
data: ds$total_ces_t1 and ds$agora_durante_a_quarentena_coronavirus_voce_faz_exercicios_fisicos
Não, nenhum Pouco (1 vez por semana) Moderado (2 a 3 vezes por semana)
Pouco (1 vez por semana) 0.0119 - -
Moderado (2 a 3 vezes por semana) 0.0032 1.0000 -
Muito (mais que 3 vezes por semana) 0.2211 1.0000 1.0000
P value adjustment method: bonferroni ggplot(ds, aes(x = sexo, y = total_ces_t1)) +
geom_boxplot() +
theme_bw()+ labs(x = "Sexo", y = "Depressão")
ggplot(ds, aes(x = sexo, y = total_gad_t1)) +
geom_boxplot() +
theme_bw()+ labs(x = "Sexo", y = "Ansiedade")1.6.4.2 Depressão / ansiedade e idade
ggplot(ds,aes(x = idade , y = total_ces_t1 )) +1.6.4.3 Depressão / ansiedade - Tem saído de casa
ds %>% count(tem_saido_de_casa)
ggplot(ds, aes(x = tem_saido_de_casa, y = total_ces_t1)) +
geom_boxplot() +
theme_bw()+ labs(x = "Tem saído de casa", y = "Depressão")
ggplot(ds, aes(x = tem_saido_de_casa, y = total_gad_t1)) +
geom_boxplot() +
theme_bw()+ labs(x = "Tem saído de casa", y = "Ansiedade")
1.6.4.4 Unir gráficos em uma única imagem
theme_set(theme_minimal())
plot_ces_gad <- ggplot(ds,aes(x = total_ces_t1, y = total_gad_t1 )) +
geom_jitter() +
geom_smooth(method = "lm") +
theme_bw() + labs(x = "Depressão", y = "Ansiedade")
plot_gad_idad <- ggplot(ds,aes(x = idade, y = total_gad_t1)) +
geom_jitter() +
geom_smooth(method = "lm") +
theme_bw() + labs(x = "Idade", y = "Ansiedade")
(plot_ces_gad + plot_gad_idad ) +theme(legend.justification = "left") 
1.7 Cronbach’s Alfa
library(psych)package 㤼㸱psych㤼㸲 was built under R version 4.0.5
Attaching package: 㤼㸱psych㤼㸲
The following objects are masked from 㤼㸱package:ggplot2㤼㸲:
%+%, alpha# CES - D
ds %>%
select(ces_t1_1:ces_t1_20) %>% psych::alpha(.) # 0,92
Reliability analysis
Call: psych::alpha(x = .)
lower alpha upper 95% confidence boundaries
0.91 0.92 0.93
Reliability if an item is dropped:
Item statistics
Non missing response frequency for each item
0 1 2 3 miss
ces_t1_1 0.39 0.28 0.24 0.10 0
ces_t1_2 0.55 0.22 0.18 0.05 0
ces_t1_3 0.58 0.17 0.14 0.12 0
ces_t1_4 0.23 0.27 0.25 0.24 0
ces_t1_5 0.28 0.25 0.26 0.21 0
ces_t1_6 0.35 0.27 0.21 0.16 0
ces_t1_7 0.41 0.29 0.19 0.11 0
ces_t1_8 0.22 0.28 0.30 0.21 0
ces_t1_9 0.54 0.23 0.11 0.12 0
ces_t1_10 0.24 0.24 0.29 0.22 0
ces_t1_11 0.22 0.25 0.31 0.22 0
ces_t1_12 0.16 0.33 0.32 0.19 0
ces_t1_13 0.34 0.29 0.22 0.15 0
ces_t1_14 0.40 0.24 0.19 0.18 0
ces_t1_15 0.59 0.24 0.13 0.05 0
ces_t1_16 0.13 0.29 0.31 0.26 0
ces_t1_17 0.55 0.24 0.13 0.08 0
ces_t1_18 0.31 0.29 0.25 0.16 0
ces_t1_19 0.62 0.22 0.09 0.07 0
ces_t1_20 0.30 0.24 0.25 0.21 0# GAD 7
ds %>% select(gad_t1_1:gad_t1_7) %>% psych::alpha(.) # 0,91
Reliability analysis
Call: psych::alpha(x = .)
lower alpha upper 95% confidence boundaries
0.9 0.91 0.93
Reliability if an item is dropped:
Item statistics
Non missing response frequency for each item
0 1 2 3 miss
gad_t1_1 0.13 0.41 0.24 0.22 0
gad_t1_2 0.21 0.39 0.21 0.19 0
gad_t1_3 0.11 0.33 0.32 0.24 0
gad_t1_4 0.16 0.31 0.31 0.22 0
gad_t1_5 0.45 0.30 0.16 0.09 0
gad_t1_6 0.29 0.35 0.20 0.16 0
gad_t1_7 0.34 0.29 0.19 0.18 01.8 Done
LS0tDQp0aXRsZTogIkNPVklELTE5IC0gRGF0YSBBbmFseXNpcyINCmF1dGhvcjogIkx1aXMgQW51bmNpYcOnw6NvLCBMdWNhcyBCYXJyb3pvIg0Kb3V0cHV0Og0KICBodG1sX25vdGVib29rOg0KICAgIHRvYzogeWVzDQogICAgdG9jX2Zsb2F0OiB5ZXMNCiAgICBudW1iZXJfc2VjdGlvbnM6IHllcw0KICAgIHRoZW1lOiB1bml0ZWQNCiAgICBoaWdobGlnaHQ6IHRleHRtYXRlDQplZGl0b3Jfb3B0aW9uczogDQogIGNodW5rX291dHB1dF90eXBlOiBpbmxpbmUNCi0tLQ0KDQo8ZGl2IGNsYXNzPSJhbGVydCBhbGVydC1pbmZvIiByb2xlPSJhbGVydCI+DQoqKk5vdGUqKjogVGhpcyBSIG1hcmtkb3duIHJlbGF0ZXMgdG8gZGF0YSBhbmFseXNpcyBvZiB0aGUgYm9vayBjaGFwdGVyIGVudGl0bGVkICJPIGN1c3RvIHBzaWNvbMOzZ2ljbyBlIGVjb27DtG1pY28gZGFzIG1lZGlkYXMgZGUgZGlzdGFuY2lhbWVudG8gc29jaWFsOiBwcmVqdcOtem8gZGUgdW0gc2Fsw6FyaW8gbcOtbmltbyBlIHBvc3PDrXZlbCBpbXBhY3RvIG5hIGRlcHJlc3PDo28iLiBUbyBhY2hpZXZlIHRoZSBzYW1lIHJlc3VsdHMsIHBsZWFzZSBkb3dubG9hZCB0aGUgUiBkYXRhc2V0IGZpbGUgYXQgaHR0cHM6Ly9vc2YuaW8vNHJ3a3AvLiBQbGVhc2UgaGF2ZSBpbiBtaW5kIHRoYXQgdGhpcyBkYXRhc2V0IGZpbGUgaXMgdGhlIHByb2Nlc3NlZCBvbmUsIGJ1dCBub3QgdGhlIG9yaWdpbmFsIHZlcnNpb24uIElmIHlvdSB3YW50IHRvIGNoZWNrIHRoZSBkYXRhIHByb2Nlc3Npbmcgc2NyaXB0cywgcGxlYXNlIHNlZSB0aGVtIGF0IGh0dHBzOi8vb3NmLmlvL2t5Zm1kLw0KDQpJZiB5b3UgaGF2ZSBhbnkgcXVlc3Rpb25zIG9yIHF1ZXJpZXMsIHBsZWFzZSByZWFjaCB1cyBvdXQgYXQgDQpsdWNhcy5wc2ljLnJqQGdtYWlsLmNvbSBvciBsdWlzZmNhQHB1Yy1yaW8uYnINCg0KbGFzdCB1cGRhdGVkOiBgciBmb3JtYXQoU3lzLnRpbWUoKSwgJyVkICVCLCAlWScpYA0KPC9kaXY+DQoNCg0KIyBEYXRhIEFuYWx5c2lzDQoNCmBgYHtyfQ0KbG9hZCh1cmwoImh0dHBzOi8vb3NmLmlvLzRyd2twL2Rvd25sb2FkIikpDQpgYGANCg0KDQoNCiMjIFBhY2thZ2VzDQoNCg0KYGBge3J9DQpwYWNtYW46OnBfbG9hZChrbml0ciwNCiAgICAgICAgICAgICAgIGthYmxlRXh0cmEsDQogICAgICAgICAgICAgICBtdWx0Y29tcCwNCiAgICAgICAgICAgICAgIHRpZHl2ZXJzZSwNCiAgICAgICAgICAgICAgIGphbml0b3IsDQogICAgICAgICAgICAgICBzdW1tYXJ5dG9vbHMsDQogICAgICAgICAgICAgICBEYXRhRXhwbG9yZXIsDQogICAgICAgICAgICAgICByZWFkeGwsDQogICAgICAgICAgICAgICBwYXRjaHdvcmspDQpgYGANCg0KDQojIyByZXN1bHRhZG9zDQoNCg0KIyMgY29ycmVsYcOnw6NvIERlcHJlc3PDo28gZSBBbnNpZWRhZGUvSWRhZGUgZSBBbnNpZWRhZGUNCg0KDQojIyBEaXN0cmlidWnDp8OjbyBkZSBkZXByZXNzw6NvIGUgYW5zaWVkYWRlDQoNCg0KYGBge3J9DQpkcyAlPiUgc2VsZWN0KHRvdGFsX2Nlc190MSwgdG90YWxfZ2FkX3QxKSAlPiUgdGFibGVvbmU6OkNyZWF0ZVRhYmxlT25lKGRhdGE9LikNCg0KYGBgDQoNCiMjIyBEZXByZXNzw6NvDQpgYGB7cn0NCmdncGxvdChkcywgYWVzKHRvdGFsX2Nlc190MSkpICsNCiAgICAgICAgIGdlb21fZGVuc2l0eShmaWxsID0gImJsdWUiLCBhbHBoYSA9IDAuMikgKw0KICB0aGVtZV9idygpICsNCiAgbGFicyh4ID0gIkRlcHJlc3PDo28iLCB5ID0gIkRlbnNpZGFkZSIpDQpgYGANCg0KIyMjIEFuc2llZGFkZQ0KDQpgYGB7cn0NCmdncGxvdChkcywgYWVzKHRvdGFsX2dhZF90MSkpICsNCiAgICAgICAgIGdlb21fZGVuc2l0eShmaWxsID0gImJsdWUiLCBhbHBoYSA9IDAuMikgKw0KICB0aGVtZV9idygpICsNCiAgbGFicyh4ID0gIkFuc2llZGFkZSIsIHkgPSAiRGVuc2lkYWRlIikNCmBgYA0KDQojIyMgRm9ybWF0byBkYSBiYXNlIGUgZmVhdHVyZSBzZWxlY3Rpb24NCg0KYGBge3J9DQpkcyAlPiUgc2VsZWN0KHNleG8gLCBpZGFkZSAsIGVzY29sYXJpZGFkZSAsIGVzdGFkb19jaXZpbCAsIHF1YW50b3NfZmlsaG9zICwgcGxhbm9fc2F1ZGUgLCBjb21vX3ZpdmUgLCBmaWNhX2VtX2Nhc2EgLCB0ZW1fc2FpZG9fZGVfY2FzYSAsIGZhbG91X2ZhbWlsaWFfcGVzc29hbG1lbnRlICwgZmFsb3VfZmFtaWxpYV90ZWxlZm9uZSAsIGZhbG91X2FtaWdvc190ZWxlZm9uZSAsIHRlc3RvdV9jb3JvbmF2aXJ1cyAsIGNvbmhlY2VfYWxndWVtX2Nvcm9uYXZpcnVzICwgc2VndWVfcmVjb21lbmRhY29lcyAsIGNvbmNvcmRvX2NvdmlkX3ByZW9jdXBhY2FvX3NhdWRlICwgY29uY29yZG9fcXVhcmVudGVuYSAsIGFyZWFfcHJvZmlzc2FvICwgcHJlb2N1cGFjYW9fcGFnYXJfY29udGEgLCBhZ29yYV9kdXJhbnRlX2FfcXVhcmVudGVuYV9jb3JvbmF2aXJ1c192b2NlX2Zhel9leGVyY2ljaW9zX2Zpc2ljb3MgLCByaW9fc3ApICU+JSBzdHINCmBgYA0KDQojIyMgTWlzc2luZyBjYXNlcyAod2l0aGluIGRzKQ0KDQpgYGB7cn0NCmRzICU+JSBzZWxlY3Qoc2V4byAsIGlkYWRlICwgZXNjb2xhcmlkYWRlICwgZXN0YWRvX2NpdmlsICwgcXVhbnRvc19maWxob3MgLCBwbGFub19zYXVkZSAsIGNvbW9fdml2ZSAsIGZpY2FfZW1fY2FzYSAsIHRlbV9zYWlkb19kZV9jYXNhICwgZmFsb3VfZmFtaWxpYV9wZXNzb2FsbWVudGUgLCBmYWxvdV9mYW1pbGlhX3RlbGVmb25lICwgZmFsb3VfYW1pZ29zX3RlbGVmb25lICwgdGVzdG91X2Nvcm9uYXZpcnVzICwgY29uaGVjZV9hbGd1ZW1fY29yb25hdmlydXMgLCBzZWd1ZV9yZWNvbWVuZGFjb2VzICwgY29uY29yZG9fY292aWRfcHJlb2N1cGFjYW9fc2F1ZGUgLCBjb25jb3Jkb19xdWFyZW50ZW5hICwgYXJlYV9wcm9maXNzYW8gLCBwcmVvY3VwYWNhb19wYWdhcl9jb250YSAsIGFnb3JhX2R1cmFudGVfYV9xdWFyZW50ZW5hX2Nvcm9uYXZpcnVzX3ZvY2VfZmF6X2V4ZXJjaWNpb3NfZmlzaWNvcyAsIHJpb19zcCkgJT4lIA0KICBwcm9maWxlX21pc3NpbmcoKQ0KYGBgDQoNCiMjIERhZG9zIGRlc2NyaXRpdm9zIChUYWJlbGEgMSkNCg0KDQpgYGB7cn0NCmRzICU+JSBzZWxlY3Qoc2V4byAsIGlkYWRlICwgZXNjb2xhcmlkYWRlICwgZXN0YWRvX2NpdmlsICwgcXVhbnRvc19maWxob3MgLCBwbGFub19zYXVkZSAsIGNvbW9fdml2ZSAsIGZpY2FfZW1fY2FzYSAsIHRlbV9zYWlkb19kZV9jYXNhICwgZmFsb3VfZmFtaWxpYV9wZXNzb2FsbWVudGUgLCBmYWxvdV9mYW1pbGlhX3RlbGVmb25lICwgZmFsb3VfYW1pZ29zX3RlbGVmb25lICwgdGVzdG91X2Nvcm9uYXZpcnVzICwgY29uaGVjZV9hbGd1ZW1fY29yb25hdmlydXMgLCBzZWd1ZV9yZWNvbWVuZGFjb2VzICwgY29uY29yZG9fY292aWRfcHJlb2N1cGFjYW9fc2F1ZGUgLCBjb25jb3Jkb19xdWFyZW50ZW5hICwgYXJlYV9wcm9maXNzYW8gLCBwcmVvY3VwYWNhb19wYWdhcl9jb250YSAsIGFnb3JhX2R1cmFudGVfYV9xdWFyZW50ZW5hX2Nvcm9uYXZpcnVzX3ZvY2VfZmF6X2V4ZXJjaWNpb3NfZmlzaWNvcyAsIHJpb19zcCwgdG90YWxfY2VzX3QxLCB0b3RhbF9nYWRfdDEpICU+JSBhcnNlbmFsOjp0YWJsZWJ5KH4uLCBkYXRhID0gLikgJT4lIHN1bW1hcnkoKQ0KYGBgDQoNCg0KIyMgbW9kZWxvcyBkZSByZWdyZXNzw6NvIChUYWJlbGEgMiBlIDMpIA0KDQojIyMgTW9kZWxvIGlycmVzdHJpdG8gLSBkZXByZXNzw6NvDQpgYGB7cn0NCnJlZ190dWRvX2RlcHJlIDwtIGxtKHRvdGFsX2Nlc190MSB+IHNleG8gKyBpZGFkZSArIGVzY29sYXJpZGFkZSArIGVzdGFkb19jaXZpbCArIHF1YW50b3NfZmlsaG9zICsgcGxhbm9fc2F1ZGUgKyBjb21vX3ZpdmUgKyBmaWNhX2VtX2Nhc2EgKyB0ZW1fc2FpZG9fZGVfY2FzYSArIGZhbG91X2ZhbWlsaWFfcGVzc29hbG1lbnRlICsgZmFsb3VfZmFtaWxpYV90ZWxlZm9uZSArIGZhbG91X2FtaWdvc190ZWxlZm9uZSArIHRlc3RvdV9jb3JvbmF2aXJ1cyArIGNvbmhlY2VfYWxndWVtX2Nvcm9uYXZpcnVzICsgc2VndWVfcmVjb21lbmRhY29lcyArIGNvbmNvcmRvX2NvdmlkX3ByZW9jdXBhY2FvX3NhdWRlICsgY29uY29yZG9fcXVhcmVudGVuYSArIGFyZWFfcHJvZmlzc2FvICsgcHJlb2N1cGFjYW9fcGFnYXJfY29udGEgKyBhZ29yYV9kdXJhbnRlX2FfcXVhcmVudGVuYV9jb3JvbmF2aXJ1c192b2NlX2Zhel9leGVyY2ljaW9zX2Zpc2ljb3MgKyByaW9fc3AsIGRzKQ0KYGBgDQoNCiMjIyMgVklGDQoNCmBgYHtyfQ0KY2FyOjp2aWYocmVnX3R1ZG9fZGVwcmUpDQpgYGANCg0KIyMjIyBFc3RpbWF0ZXMgKEJldGFzKQ0KDQpgYGB7cn0NCm9sc3JyOjpvbHNfcmVncmVzcyhyZWdfdHVkb19kZXByZSkgDQpgYGANCg0KIyMjIyBFc3RpbWF0ZXMgKFNTUSkNCg0KYGBge3J9DQphcGFUYWJsZXM6OmFwYS5hb3YudGFibGUocmVnX3R1ZG9fZGVwcmUpDQpgYGANCg0KIyMjIFJlZyB0ZXN0ZSBjb20gdHVkbyAtIEFuc2llZGFkZQ0KDQpgYGB7cn0NCnJlZ190dWRvX2Fuc2kgPC0gbG0odG90YWxfZ2FkX3QxIH4gc2V4byArIGlkYWRlICsgZXNjb2xhcmlkYWRlICsgZXN0YWRvX2NpdmlsICsgcXVhbnRvc19maWxob3MgKyBwbGFub19zYXVkZSArIGNvbW9fdml2ZSArIGZpY2FfZW1fY2FzYSArIHRlbV9zYWlkb19kZV9jYXNhICsgZmFsb3VfZmFtaWxpYV9wZXNzb2FsbWVudGUgKyBmYWxvdV9mYW1pbGlhX3RlbGVmb25lICsgZmFsb3VfYW1pZ29zX3RlbGVmb25lICsgdGVzdG91X2Nvcm9uYXZpcnVzICsgY29uaGVjZV9hbGd1ZW1fY29yb25hdmlydXMgKyBzZWd1ZV9yZWNvbWVuZGFjb2VzICsgY29uY29yZG9fY292aWRfcHJlb2N1cGFjYW9fc2F1ZGUgKyBjb25jb3Jkb19xdWFyZW50ZW5hICsgYXJlYV9wcm9maXNzYW8gKyBwcmVvY3VwYWNhb19wYWdhcl9jb250YSArIGFnb3JhX2R1cmFudGVfYV9xdWFyZW50ZW5hX2Nvcm9uYXZpcnVzX3ZvY2VfZmF6X2V4ZXJjaWNpb3NfZmlzaWNvcyArIHJpb19zcCwgZHMgKQ0KDQpgYGANCg0KIyMjIFZJRg0KYGBge3J9DQpjYXI6OnZpZihyZWdfdHVkb19hbnNpKQ0KDQpgYGANCg0KIyMjIyBFc3RpbWF0ZXMgKGJldGEpDQoNCmBgYHtyfQ0Kb2xzcnI6Om9sc19yZWdyZXNzKHJlZ190dWRvX2Fuc2kpDQpgYGANCg0KDQojIyMjIEVzdGltYXRlcyAoU1NRKQ0KDQpgYGB7cn0NCmFwYVRhYmxlczo6YXBhLnJlZy50YWJsZShyZWdfdHVkb19hbnNpKQ0KDQoNCnN1bW1hcnkocmVnX3R1ZG9fYW5zaSkNCmBgYA0KDQoNCg0KIyMjIyBQb3N0IGhvY3MNCg0KYGBge3J9DQpnbGh0KHJlZ190dWRvX2RlcHJlLCBsaW5mY3QgPSBtY3Aoc2V4byA9ICdUdWtleScpKSAlPiUgc3VtbWFyeSgpDQpgYGANCg0KYGBge3J9DQpUdWtleUhTRChhb3YodG90YWxfY2VzX3QxIH4gc2V4byArIGVzY29sYXJpZGFkZSArIGVzdGFkb19jaXZpbCArIHF1YW50b3NfZmlsaG9zICsgcGxhbm9fc2F1ZGUgKyBjb21vX3ZpdmUgKyBmaWNhX2VtX2Nhc2EgKyB0ZW1fc2FpZG9fZGVfY2FzYSArIGZhbG91X2ZhbWlsaWFfcGVzc29hbG1lbnRlICsgZmFsb3VfZmFtaWxpYV90ZWxlZm9uZSArIGZhbG91X2FtaWdvc190ZWxlZm9uZSArIHRlc3RvdV9jb3JvbmF2aXJ1cyArIGNvbmhlY2VfYWxndWVtX2Nvcm9uYXZpcnVzICsgc2VndWVfcmVjb21lbmRhY29lcyArIGNvbmNvcmRvX2NvdmlkX3ByZW9jdXBhY2FvX3NhdWRlICsgY29uY29yZG9fcXVhcmVudGVuYSArIGFyZWFfcHJvZmlzc2FvICsgcHJlb2N1cGFjYW9fcGFnYXJfY29udGEgKyBhZ29yYV9kdXJhbnRlX2FfcXVhcmVudGVuYV9jb3JvbmF2aXJ1c192b2NlX2Zhel9leGVyY2ljaW9zX2Zpc2ljb3MgKyByaW9fc3AsIGRzKSkNCmBgYA0KDQoNCmBgYHtyfQ0KIyByZWdyZXNzw6NvIGRlIGRlcHJlc3PDo28NCmdsaHQocmVnX3R1ZG9fZGVwcmUsIGxpbmZjdCA9IG1jcChhZ29yYV9kdXJhbnRlX2FfcXVhcmVudGVuYV9jb3JvbmF2aXJ1c192b2NlX2Zhel9leGVyY2ljaW9zX2Zpc2ljb3MgPSAnVHVrZXknKSkgJT4lIHN1bW1hcnkoKQ0KDQojIHJlZ3Jlc3PDo28gZGUgYW5zaWVkYWRlDQpnbGh0KHJlZ190dWRvX2Fuc2ksIGxpbmZjdCA9IG1jcChhZ29yYV9kdXJhbnRlX2FfcXVhcmVudGVuYV9jb3JvbmF2aXJ1c192b2NlX2Zhel9leGVyY2ljaW9zX2Zpc2ljb3MgPSAnVHVrZXknKSkgJT4lIHN1bW1hcnkoKQ0KDQpgYGANCg0KDQoNCg0KDQojIyMgUGxvdHMNCg0KIyMjIyBEZXByZXNzw6NvIC8gYW5zaWVkYWRlIC0gU2V4byANCmBgYHtyfQ0KcGFpcndpc2UudC50ZXN0KHggPSBkcyR0b3RhbF9jZXNfdDEsIGcgPSBkcyRhZ29yYV9kdXJhbnRlX2FfcXVhcmVudGVuYV9jb3JvbmF2aXJ1c192b2NlX2Zhel9leGVyY2ljaW9zX2Zpc2ljb3MsIHAuYWRqdXN0Lm1ldGhvZCAgPSAiYm9uZmVycm9uaSIsIC4uLiA9IGRzKSANCg0KZ2dwbG90KGRzLCBhZXMoeCA9IHNleG8sIHkgPSB0b3RhbF9jZXNfdDEpKSArDQogIGdlb21fYm94cGxvdCgpICsNCiAgdGhlbWVfYncoKSsgbGFicyh4ID0gIlNleG8iLCB5ID0gIkRlcHJlc3PDo28iKQ0KDQpnZ3Bsb3QoZHMsIGFlcyh4ID0gc2V4bywgeSA9IHRvdGFsX2dhZF90MSkpICsNCiAgZ2VvbV9ib3hwbG90KCkgKw0KICB0aGVtZV9idygpKyBsYWJzKHggPSAiU2V4byIsIHkgPSAiQW5zaWVkYWRlIikNCg0KYGBgDQoNCiMjIyMgRGVwcmVzc8OjbyAvIGFuc2llZGFkZSBlIGlkYWRlIA0KYGBge3J9DQpnZ3Bsb3QoZHMsYWVzKHggPSBpZGFkZSAsIHkgPSB0b3RhbF9jZXNfdDEgKSkgKw0KICBnZW9tX2ppdHRlcigpICsNCiAgZ2VvbV9zbW9vdGgobWV0aG9kID0gImxtIikgKw0KICB0aGVtZV9idygpICsgbGFicyh4ID0gIklkYWRlIiwgeSA9ICJEZXByZXNzw6NvIikNCg0KZ2dwbG90KGRzLGFlcyh4ID0gaWRhZGUsIHkgPSB0b3RhbF9nYWRfdDEpKSArDQogIGdlb21faml0dGVyKCkgKw0KICBnZW9tX3Ntb290aChtZXRob2QgPSAibG0iKSArDQogIHRoZW1lX2J3KCkgKyBsYWJzKHggPSAiSWRhZGUiLCB5ID0gIkFuc2llZGFkZSIpDQoNCmBgYA0KDQoNCiMjIyMgRGVwcmVzc8OjbyAvIGFuc2llZGFkZSAtIFRlbSBzYcOtZG8gZGUgY2FzYSANCmBgYHtyfQ0KDQpkcyAlPiUgY291bnQodGVtX3NhaWRvX2RlX2Nhc2EpDQoNCmdncGxvdChkcywgYWVzKHggPSB0ZW1fc2FpZG9fZGVfY2FzYSwgeSA9IHRvdGFsX2Nlc190MSkpICsNCiAgZ2VvbV9ib3hwbG90KCkgKw0KICB0aGVtZV9idygpKyBsYWJzKHggPSAiVGVtIHNhw61kbyBkZSBjYXNhIiwgeSA9ICJEZXByZXNzw6NvIikNCg0KZ2dwbG90KGRzLCBhZXMoeCA9IHRlbV9zYWlkb19kZV9jYXNhLCB5ID0gdG90YWxfZ2FkX3QxKSkgKw0KICBnZW9tX2JveHBsb3QoKSArDQogIHRoZW1lX2J3KCkrIGxhYnMoeCA9ICJUZW0gc2HDrWRvIGRlIGNhc2EiLCB5ID0gIkFuc2llZGFkZSIpDQoNCmBgYA0KDQoNCmBgYHtyfQ0KZHMgJT4lIA0KICBmaWx0ZXIoIWlzLm5hKGNvbmhlY2VfYWxndWVtX2Nvcm9uYXZpcnVzKSkgJT4lIA0KICBnZ3Bsb3QoYWVzKHggPSBjb25oZWNlX2FsZ3VlbV9jb3JvbmF2aXJ1cywgeSA9IHRvdGFsX2Nlc190MSkpICsNCiAgZ2VvbV9ib3hwbG90KCkgKw0KICB0aGVtZV9idygpKyBsYWJzKHggPSAiQ29uaGVjZSBBbGd1w6ltIGNvbSBDb3ZpZCIsIHkgPSAiRGVwcmVzc8OjbyIpDQpgYGANCg0KIyMjIyBVbmlyIGdyw6FmaWNvcyBlbSB1bWEgw7puaWNhIGltYWdlbQ0KYGBge3J9DQp0aGVtZV9zZXQodGhlbWVfbWluaW1hbCgpKQ0KDQpwbG90X2Nlc19nYWQgPC0gIGdncGxvdChkcyxhZXMoeCA9IHRvdGFsX2Nlc190MSwgeSA9IHRvdGFsX2dhZF90MSApKSArDQogIGdlb21faml0dGVyKCkgKw0KICBnZW9tX3Ntb290aChtZXRob2QgPSAibG0iKSArDQogIHRoZW1lX2J3KCkgKyBsYWJzKHggPSAiRGVwcmVzc8OjbyIsIHkgPSAiQW5zaWVkYWRlIikNCg0KDQpwbG90X2dhZF9pZGFkIDwtICBnZ3Bsb3QoZHMsYWVzKHggPSBpZGFkZSwgeSA9IHRvdGFsX2dhZF90MSkpICsNCiAgZ2VvbV9qaXR0ZXIoKSArDQogIGdlb21fc21vb3RoKG1ldGhvZCA9ICJsbSIpICsNCiAgdGhlbWVfYncoKSArIGxhYnMoeCA9ICJJZGFkZSIsIHkgPSAiQW5zaWVkYWRlIikNCg0KDQoNCihwbG90X2Nlc19nYWQgKyBwbG90X2dhZF9pZGFkICkgK3RoZW1lKGxlZ2VuZC5qdXN0aWZpY2F0aW9uID0gImxlZnQiKSAgICAgICAgICAgICAgDQpgYGANCg0KIyMgQ3JvbmJhY2gncyBBbGZhDQoNCmBgYHtyfQ0KbGlicmFyeShwc3ljaCkNCiMgQ0VTIC0gRA0KZHMgJT4lIA0KICBzZWxlY3QoY2VzX3QxXzE6Y2VzX3QxXzIwKSAlPiUgcHN5Y2g6OmFscGhhKC4pICMgMCw5Mg0KDQojIEdBRCA3IA0KDQpkcyAlPiUgc2VsZWN0KGdhZF90MV8xOmdhZF90MV83KSAlPiUgcHN5Y2g6OmFscGhhKC4pICMgMCw5MQ0KDQpgYGANCg0KDQojIyBEb25lDQo=