Asociaciones Univariadas
table1_all <- suppressWarnings(compareGroups(tamizaje_dep ~
p45_demo_sexo+
estamento+
p21_cond_trab_hr_trabajada+
p8_exp_c19_cerca_cont_rec,
method= c(
sexo=3,
estamento=3,
p21_cond_trab_hr_trabajada=2,
p11_exp_c19_miedo_pand_rec=3
),
data = bd_luciano,
include.miss = F, #baja el número de casos
p.corrected=T,
chisq.test.perm=T,
chisq.test.seed=1234,
var.equal=T,
byrow = T)
)
restab1_all <- createTable(table1_all, show.p.overall = T)
pvals1 <- getResults(table1_all)
#p.adjust(pvals, method = "BH")
export2md(restab1_all, size=12, first.strip=T, hide.no="no", position="center",
format="html",caption= "Tabla 1. Comparativo Asociaciones Bivariadas",
col.names=c("Variables", "Ausencia de Indicadores de Sintomatología Depresiva", "Presencia de Indicadores de Sintomatología Depresiva", "Sig.")
)%>%
kableExtra::add_footnote(c("Nota. N= Número de participantes con al menos una respuesta válida en las variables seleccionadas; Variables continuas se presentan como medianas y percentiles 25 y 75;",
"Variables categóricas son presentadas como el recuento y el porcentaje",
"Porcentajes por fila"), notation = "none")%>%
kableExtra::scroll_box(width = "100%", height = "800px")
Tabla 1. Comparativo Asociaciones Bivariadas
|
Variables
|
Ausencia de Indicadores de Sintomatología Depresiva
|
Presencia de Indicadores de Sintomatología Depresiva
|
Sig.
|
|
|
N=277
|
N=113
|
|
|
Sexo (Recodificación Binaria):
|
|
|
0.018
|
|
Hombre
|
68 (80.0%)
|
17 (20.0%)
|
|
|
Mujer
|
143 (65.3%)
|
76 (34.7%)
|
|
|
Estamento:
|
|
|
0.298
|
|
Administrativo
|
33 (70.2%)
|
14 (29.8%)
|
|
|
Enfermeros
|
23 (57.5%)
|
17 (42.5%)
|
|
|
Médicos
|
14 (77.8%)
|
4 (22.2%)
|
|
|
Otros Profesionales No Médicos
|
110 (68.8%)
|
50 (31.2%)
|
|
|
Técnicos y Auxiliares
|
41 (77.4%)
|
12 (22.6%)
|
|
|
¿Cuántas horas trabajó la semana pasada en su lugar de trabajo principal? (en números)
|
44.0 [30.0;46.0]
|
44.0 [24.0;45.0]
|
0.403
|
|
¿Qué tan cercano es el contacto que usted tiene con pacientes con COVID-19?(dicotomizada):
|
|
|
0.423
|
|
Contacto directo y estrecho con pacientes (menos de dos metros)
|
59 (67.0%)
|
29 (33.0%)
|
|
|
Menor o ningún contacto con pacientes
|
218 (72.2%)
|
84 (27.8%)
|
|
|
Nota. N= Número de participantes con al menos una respuesta válida en las variables seleccionadas; Variables continuas se presentan como medianas y percentiles 25 y 75;
|
|
Variables categóricas son presentadas como el recuento y el porcentaje
|
|
Porcentajes por fila
|
vec_y<-c("p45_demo_sexo", "estamento", "p21_cond_trab_hr_trabajada", "p8_exp_c19_cerca_cont_rec")
x<-"tamizaje_dep"
y<-vec_y[i]
paste0("#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:")
## [1] "#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:"
paste0("**### Exposición: Sexo;\n Var. de Resultado: Tamizaje Depresión**")
## [1] "**### Exposición: Sexo;\n Var. de Resultado: Tamizaje Depresión**"
chisq.test(table(bd_luciano$p45_demo_sexo,bd_luciano$tamizaje_dep)[,c(2,1)])
##
## Pearson's Chi-squared test with Yates' continuity correction
##
## data: table(bd_luciano$p45_demo_sexo, bd_luciano$tamizaje_dep)[, c(2, 1)]
## X-squared = 5.5612, df = 1, p-value = 0.01836
#epiR::epi.2by2(table(bd_luciano$p45_demo_sexo,bd_luciano$tamizaje_dep)[,c(2,1)], method = "cross.sectional", conf.level = 0.95, outcome = "as.columns")
paste0("#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:")
## [1] "#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:"
paste0("**### Exposición: Estamento;\n Var. de Resultado: Tamizaje Depresión**")
## [1] "**### Exposición: Estamento;\n Var. de Resultado: Tamizaje Depresión**"
table(bd_luciano$estamento,bd_luciano$tamizaje_dep)[,c(2,1)]
##
## Indicadores de Sintomatología Depresiva
## Administrativo 14
## Enfermeros 17
## Médicos 4
## Otros Profesionales No Médicos 50
## Técnicos y Auxiliares 12
##
## Ausencia de Indicadores de Sintomatología Depresiva
## Administrativo 33
## Enfermeros 23
## Médicos 14
## Otros Profesionales No Médicos 110
## Técnicos y Auxiliares 41
chisq.test(table(bd_luciano$estamento,bd_luciano$tamizaje_dep)[,c(2,1)])
##
## Pearson's Chi-squared test
##
## data: table(bd_luciano$estamento, bd_luciano$tamizaje_dep)[, c(2, 1)]
## X-squared = 4.8967, df = 4, p-value = 0.2981
paste0("#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:")
## [1] "#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:"
paste0("**### Exposición: Horas Trabajadas;\n Var. de Resultado: Tamizaje Depresión**")
## [1] "**### Exposición: Horas Trabajadas;\n Var. de Resultado: Tamizaje Depresión**"
#kruskal.test(p21_cond_trab_hr_trabajada ~ tamizaje_dep, data = bd_luciano)
coin::wilcox_test(p21_cond_trab_hr_trabajada ~ tamizaje_dep,
data = bd_luciano, paired=F)
##
## Asymptotic Wilcoxon-Mann-Whitney Test
##
## data: p21_cond_trab_hr_trabajada by
## tamizaje_dep (Ausencia de Indicadores de Sintomatología Depresiva, Indicadores de Sintomatología Depresiva)
## Z = 0.83545, p-value = 0.4035
## alternative hypothesis: true mu is not equal to 0
paste0("#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:")
## [1] "#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:#:"
paste0("**### Exposición: Cercanía Pacientes COVID-19;\n Var. de Resultado: Tamizaje Depresión**")
## [1] "**### Exposición: Cercanía Pacientes COVID-19;\n Var. de Resultado: Tamizaje Depresión**"
table(bd_luciano$p8_exp_c19_cerca_cont_rec,bd_luciano$tamizaje_dep)[,c(2,1)]
##
## Indicadores de Sintomatología Depresiva
## Contacto directo y estrecho con pacientes (menos de dos metros) 29
## Menor o ningún contacto con pacientes 84
##
## Ausencia de Indicadores de Sintomatología Depresiva
## Contacto directo y estrecho con pacientes (menos de dos metros) 59
## Menor o ningún contacto con pacientes 218
chisq.test(table(bd_luciano$p8_exp_c19_cerca_cont_rec,bd_luciano$tamizaje_dep)[,c(2,1)])
##
## Pearson's Chi-squared test with Yates' continuity correction
##
## data: table(bd_luciano$p8_exp_c19_cerca_cont_rec, bd_luciano$tamizaje_dep)[, c(2, 1)]
## X-squared = 0.64288, df = 1, p-value = 0.4227
#epiR::epi.2by2(table(bd_luciano$p8_exp_c19_cerca_cont_rec,bd_luciano$tamizaje_dep)[,c(2,1)], method = "cross.sectional", conf.level = 0.95, outcome = "as.columns")
Asociaciones Multivariadas
glm_60d <- glm(tamizaje_dep ~ p8_exp_c19_cerca_cont_rec*p45_demo_sexo ,#+estamento,
data=bd_luciano, family="binomial")
library(finalfit)
#explanatory = c("age.factor*sex.factor", "obstruct.factor", "perfor.factor")
bd_luciano %>%
dplyr::mutate(p45_demo_sexo= ff_label(p45_demo_sexo, "Sexo"),
p8_exp_c19_cerca_cont_rec= ff_label(p8_exp_c19_cerca_cont_rec, "Contacto Pac. COVID19"),
tamizaje_dep= ff_label(tamizaje_dep, "Presencia Sint. Depresiva"))%>%
ff_interaction(p45_demo_sexo, p8_exp_c19_cerca_cont_rec) %>%
finalfit("tamizaje_dep", c("p45_demo_sexo*p8_exp_c19_cerca_cont_rec"),dependent_label_prefix = "") -> t
## Note: dependent includes missing data. These are dropped.
## Waiting for profiling to be done...
## Waiting for profiling to be done...
knitr::kable(t,format = "html", format.args = list(decimal.mark = ".", big.mark = ","),
caption="Tabla 2. Asociación entre Sintomatología Depresiva y Sexo, Moderado por Exposición a COVID-19",
align =rep('c', 101)) %>%
kableExtra::kable_styling(bootstrap_options = c("striped", "hover"),font_size = 12) %>%
kableExtra::scroll_box(width = "100%", height = "350px")
Tabla 2. Asociación entre Sintomatología Depresiva y Sexo, Moderado por Exposición a COVID-19
|
|
Presencia Sint. Depresiva
|
|
Ausencia de Indicadores de Sintomatología Depresiva
|
Indicadores de Sintomatología Depresiva
|
OR (univariable)
|
OR (multivariable)
|
|
1
|
Sexo
|
Hombre
|
68 (80.0)
|
17 (20.0)
|
|
|
|
2
|
|
Mujer
|
143 (65.3)
|
76 (34.7)
|
3.59 (1.00-17.12, p=0.069)
|
3.59 (1.00-17.12, p=0.069)
|
|
4
|
Contacto Pac. COVID19
|
Contacto directo y estrecho con pacientes (menos de dos metros)
|
59 (67.0)
|
29 (33.0)
|
|
|
|
5
|
|
Menor o ningún contacto con pacientes
|
218 (72.2)
|
84 (27.8)
|
1.21 (0.34-5.76, p=0.786)
|
1.21 (0.34-5.76, p=0.786)
|
|
3
|
p45_demo_sexoMujer:p8_exp_c19_cerca_cont_recMenor o ningún contacto con pacientes
|
Interaction
|
|
|
0.51 (0.10-2.20, p=0.394)
|
0.51 (0.10-2.20, p=0.394)
|
