El análisis se realizó comenzando por importar las librerias requeridas entre ellas readxl, nortest, ggplot2, dplyr, tidyverse, tableone, sandwich, lmtest, cowplot, gmodels, crosstable, haven, foreign y sandwich

Posteriormente se codificaron todas las variables correspondientes siendo edad, propofol, fentanilo, sangrado las variables cuantitativas numericas de tipo continuo y las variables de ASA, hipertensión, diabetes, hipotiroidismo, delirio y duración de la administración de la anestesia de tipo cualitativo.

Se aplicó la prueba de normalidad de Shapiro-Wilk a cada variable, debido a que la prueba de Kolmogorov-Smirnov en muestras grandes no es tan exacta y se opta por la prueba de lilliefors, en cambio para análisis estadisticos en R y Rstudio se recomienda emplear Shapiro-Wilk en muestras de 0-5mil observaciones. Se aprecia que las variables poseen una distribución libre.

## 
##  Shapiro-Wilk normality test
## 
## data:  edad
## W = 0.9442, p-value = 2.98e-05

## 
##  Shapiro-Wilk normality test
## 
## data:  propofol
## W = 0.87723, p-value = 3.471e-09

## 
##  Shapiro-Wilk normality test
## 
## data:  datos$fentanilo_mcg
## W = 0.9579, p-value = 0.0003631

## 
##  Shapiro-Wilk normality test
## 
## data:  datos$sangrado_ml
## W = 0.97186, p-value = 0.006732

Se realizó una tabla descriptiva que resume las caracteristicas generales de la muestra estudiada asi como de las variables, las variables cuantitativas se expresan en mediana (Q1-Q3) y las proporciones en frecuencias (%).

Tabla 1 Descripción general de las variables de estudio
level Overall
n 135
edad (median [IQR]) 68.00 [63.00, 73.00]
ASA (%) ASA II 82 (60.74)
ASA III 53 (39.26)
hipertension (%) sin hipertensión 61 (45.19)
Con hipertensión 74 (54.81)
diabetes_mellitus (%) Sin diabetes tipo 2 82 (60.74)
Con diabetes tipo 2 53 (39.26)
hipotiroidismo (%) Sin hipotiroidismo 127 (94.07)
Con hipotiroidismo 8 ( 5.93)
propofol_mg (median [IQR]) 110.00 [100.00, 120.00]
fentanilo_mcg (median [IQR]) 825.00 [622.50, 998.00]
sangrado_ml (median [IQR]) 450.00 [390.00, 515.00]
duracion_gt180 (%) <180 minutos 69 (51.11)
>180 minutos 66 (48.89)
delirio (%) sin delirio 110 (81.48)
con delirio 25 (18.52)

Se realizó la construcción de gráficos de cajas y brazos para cada una de las variables cuantitativas y permitir observar su distribución.

Se construyeron adicionalmente gráficos de barras para observar la distribución de las variables cualitativas.

Con motivo de efectuar una comparación entre grupos entre aquellos individuos que presentaron delirio y aquellos que no, se realizaron tablas cruzadas para conocer las frecuencias de cada variable cualitativa en cada grupo y contando con una frecuencia menor de 5 aplicar la prueba exacta de Fisher y para una frecuencia mayor de 5 la prueba de Chi-cuadrada de Pearson.

## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## |              Expected N |
## |-------------------------|
## 
##  
## Total Observations in Table:  135 
## 
##  
##              | delirio 
##          ASA | sin delirio | con delirio |   Row Total | 
## -------------|-------------|-------------|-------------|
##       ASA II |          64 |          18 |          82 | 
##              |      66.815 |      15.185 |             | 
## -------------|-------------|-------------|-------------|
##      ASA III |          46 |           7 |          53 | 
##              |      43.185 |       9.815 |             | 
## -------------|-------------|-------------|-------------|
## Column Total |         110 |          25 |         135 | 
## -------------|-------------|-------------|-------------|
## 
##  
## Statistics for All Table Factors
## 
## 
## Pearson's Chi-squared test 
## ------------------------------------------------------------
## Chi^2 =  1.631092     d.f. =  1     p =  0.201552 
## 
## Pearson's Chi-squared test with Yates' continuity correction 
## ------------------------------------------------------------
## Chi^2 =  1.103091     d.f. =  1     p =  0.2935887 
## 
## 
## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## |              Expected N |
## |-------------------------|
## 
##  
## Total Observations in Table:  135 
## 
##  
##                  | delirio 
##     hipertensión | sin delirio | con delirio |   Row Total | 
## -----------------|-------------|-------------|-------------|
## sin hipertensión |          55 |           6 |          61 | 
##                  |      49.704 |      11.296 |             | 
## -----------------|-------------|-------------|-------------|
## Con hipertensión |          55 |          19 |          74 | 
##                  |      60.296 |      13.704 |             | 
## -----------------|-------------|-------------|-------------|
##     Column Total |         110 |          25 |         135 | 
## -----------------|-------------|-------------|-------------|
## 
##  
## Statistics for All Table Factors
## 
## 
## Pearson's Chi-squared test 
## ------------------------------------------------------------
## Chi^2 =  5.559703     d.f. =  1     p =  0.01837853 
## 
## Pearson's Chi-squared test with Yates' continuity correction 
## ------------------------------------------------------------
## Chi^2 =  4.559519     d.f. =  1     p =  0.03273626 
## 
## 
## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## |              Expected N |
## |-------------------------|
## 
##  
## Total Observations in Table:  135 
## 
##  
##                     | delirio 
##            diabetes | sin delirio | con delirio |   Row Total | 
## --------------------|-------------|-------------|-------------|
## Sin diabetes tipo 2 |          67 |          15 |          82 | 
##                     |      66.815 |      15.185 |             | 
## --------------------|-------------|-------------|-------------|
## Con diabetes tipo 2 |          43 |          10 |          53 | 
##                     |      43.185 |       9.815 |             | 
## --------------------|-------------|-------------|-------------|
##        Column Total |         110 |          25 |         135 | 
## --------------------|-------------|-------------|-------------|
## 
##  
## Statistics for All Table Factors
## 
## 
## Pearson's Chi-squared test 
## ------------------------------------------------------------
## Chi^2 =  0.007059783     d.f. =  1     p =  0.9330385 
## 
## Pearson's Chi-squared test with Yates' continuity correction 
## ------------------------------------------------------------
## Chi^2 =  1.284923e-30     d.f. =  1     p =  1 
## 
## 
## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## |              Expected N |
## |-------------------------|
## 
##  
## Total Observations in Table:  135 
## 
##  
##                    | delirio 
##     hipotiroidismo | sin delirio | con delirio |   Row Total | 
## -------------------|-------------|-------------|-------------|
## Sin hipotiroidismo |         103 |          24 |         127 | 
##                    |     103.481 |      23.519 |             | 
## -------------------|-------------|-------------|-------------|
## Con hipotiroidismo |           7 |           1 |           8 | 
##                    |       6.519 |       1.481 |             | 
## -------------------|-------------|-------------|-------------|
##       Column Total |         110 |          25 |         135 | 
## -------------------|-------------|-------------|-------------|
## 
##  
## Statistics for All Table Factors
## 
## 
## Pearson's Chi-squared test 
## ------------------------------------------------------------
## Chi^2 =  0.2041428     d.f. =  1     p =  0.6513975 
## 
## Pearson's Chi-squared test with Yates' continuity correction 
## ------------------------------------------------------------
## Chi^2 =  2.541387e-31     d.f. =  1     p =  1 
## 
## 
## 
##  
##    Cell Contents
## |-------------------------|
## |                       N |
## |              Expected N |
## |-------------------------|
## 
##  
## Total Observations in Table:  135 
## 
##  
##              | delirio 
##     duracion | sin delirio | con delirio |   Row Total | 
## -------------|-------------|-------------|-------------|
## <180 minutos |          61 |           8 |          69 | 
##              |      56.222 |      12.778 |             | 
## -------------|-------------|-------------|-------------|
## >180 minutos |          49 |          17 |          66 | 
##              |      53.778 |      12.222 |             | 
## -------------|-------------|-------------|-------------|
## Column Total |         110 |          25 |         135 | 
## -------------|-------------|-------------|-------------|
## 
##  
## Statistics for All Table Factors
## 
## 
## Pearson's Chi-squared test 
## ------------------------------------------------------------
## Chi^2 =  4.484639     d.f. =  1     p =  0.03420077 
## 
## Pearson's Chi-squared test with Yates' continuity correction 
## ------------------------------------------------------------
## Chi^2 =  3.595109     d.f. =  1     p =  0.05794984 
## 
##

Los resultados de la prueba de normalidad indicaron una distribución libre por lo que las pruebas para comparación entre grupos se realizaron con la prueba U de Mann-Whitney exceptuando la prueba de homogeneidad de varianzas.

## 
##  Wilcoxon rank sum test with continuity correction
## 
## data:  edad by delirio
## W = 1539.5, p-value = 0.3521
## alternative hypothesis: true location shift is not equal to 0
## 
##  Wilcoxon rank sum test with continuity correction
## 
## data:  propofol_mg by delirio
## W = 1656, p-value = 0.1046
## alternative hypothesis: true location shift is not equal to 0
## 
##  Wilcoxon rank sum test with continuity correction
## 
## data:  fentanilo_mcg by delirio
## W = 1818, p-value = 0.01219
## alternative hypothesis: true location shift is not equal to 0
## 
##  Wilcoxon rank sum test with continuity correction
## 
## data:  sangrado_ml by delirio
## W = 1311, p-value = 0.7189
## alternative hypothesis: true location shift is not equal to 0

Habiéndose detectado únicamente detectado una diferencia estadistícamente significativa en la variable de fentanilo (microgramos). Los resultados de la comparación entre grupos con la aplicación de la prueba Chi-cuadrada de Pearson y Exacta de Fisher con corrección de Yates se observan a continuación.

## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  ASA and delirio
## X-squared = 1.1031, df = 1, p-value = 0.2936
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  hipertensión and delirio
## X-squared = 4.5595, df = 1, p-value = 0.03274
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  diabetes and delirio
## X-squared = 1.2849e-30, df = 1, p-value = 1
## 
##  Fisher's Exact Test for Count Data
## 
## data:  hipotiroidismo and delirio
## p-value = 1
## alternative hypothesis: true odds ratio is not equal to 1
## 95 percent confidence interval:
##  0.01307026 5.16796254
## sample estimates:
## odds ratio 
##  0.6150358
## 
##  Pearson's Chi-squared test with Yates' continuity correction
## 
## data:  duracion and delirio
## X-squared = 3.5951, df = 1, p-value = 0.05795

Para observar la distribución en ambos grupos de las variables cuantitativas y cualitativas se consturyeron gráficos de cajas comparandose por presencia y ausencia de delirio en la muestra estudiada, al igual que se construyeron gráficos de barras.

Tabla 2: Comparación de variables según presencia de Delirio Postoperatorio
level sin delirio con delirio p test
n 110 25
edad (median [IQR]) 68.00 [63.00, 73.00] 66.00 [63.00, 71.00] 0.351 nonnorm
ASA (%) ASA II 64 (58.18) 18 (72.00) 0.294
ASA III 46 (41.82) 7 (28.00)
hipertension (%) sin hipertensión 55 (50.00) 6 (24.00) 0.033
Con hipertensión 55 (50.00) 19 (76.00)
diabetes_mellitus (%) Sin diabetes tipo 2 67 (60.91) 15 (60.00) 1.000
Con diabetes tipo 2 43 (39.09) 10 (40.00)
hipotiroidismo (%) Sin hipotiroidismo 103 (93.64) 24 (96.00) 1.000
Con hipotiroidismo 7 (6.36) 1 (4.00)
propofol_mg (median [IQR]) 110.00 [100.00, 120.00] 100.00 [100.00, 110.00] 0.104 nonnorm
fentanilo_mcg (median [IQR]) 853.50 [648.25, 1034.25] 669.00 [533.00, 856.00] 0.012 nonnorm
sangrado_ml (median [IQR]) 450.00 [390.00, 510.00] 470.00 [400.00, 520.00] 0.717 nonnorm
duracion_gt180 (%) <180 minutos 61 (55.45) 8 (32.00) 0.058
>180 minutos 49 (44.55) 17 (68.00)
delirio (%) sin delirio 110 (100.00) 0 (0.00) <0.001
con delirio 0 (0.00) 25 (100.00)

Con motivo de realizar un análisis de regresión de Poisson se recodifico la variable delirio para seleccionar como dependiente la variable con el valor 1 que fue de aquellos en los que se observó delirio.

#recodificación para modelos de regresión
datos$delirio_num = as.numeric(datos$delirio) - 1
delirionum = datos$delirio_num

Ya con la variable codificada nuevamente se procedió a realizar la construcción de un análisis univariado con cada una de las variables.

##        Variable    RR    LI    LS p_valor
##            Edad 0.967 0.898 1.040  0.3660
##        Propofol 0.979 0.952 1.007  0.1328
##       Fentanilo 0.998 0.996 1.000  0.0284
##        Sangrado 1.001 0.996 1.006  0.7646
##             ASA 0.602 0.251 1.441  0.2541
##    Hipertension 2.610 1.043 6.536  0.0405
##        Diabetes 1.031 0.463 2.296  0.9395
##  Hipotiroidismo 0.661 0.089 4.889  0.6855
##        Duracion 2.222 0.959 5.148  0.0626

Posteriormente se construyó un forest plot para observar el RR con su respectivo intervalo de confianza al 95% de cada variable incluida en el análisis de regresión, como se muestra acontinuación: Nuevamente se construyeron modelos de regresión de Poisson, para ello se debe mencionar que los resultados se expresan en RR y su respectivo intervalo de confianza al 95%.

##                                RRa 2.5 %   97.5 % p_valor
## (Intercept)                  7.381 0.036 1500.580  0.4610
## edad                         0.959 0.891    1.033  0.2743
## ASAASA III                   0.687 0.285    1.658  0.4039
## fentanilo_mcg                0.998 0.997    1.000  0.0371
## hipertensionCon hipertensión 2.314 0.914    5.857  0.0766

##                                RRa 2.5 % 97.5 % p_valor
## (Intercept)                  0.366 0.087  1.547  0.1717
## fentanilo_mcg                0.998 0.997  1.000  0.0426
## hipertensionCon hipertensión 2.444 0.973  6.141  0.0572

##                                    RRa        LI          LS     p_valor
## (Intercept)                  7.3814997 0.1003674 542.8708108 0.361975166
## edad                         0.9594736 0.9026934   1.0198254 0.183769002
## ASAASA III                   0.6873613 0.3142364   1.5035352 0.347844942
## fentanilo_mcg                0.9982458 0.9969926   0.9995006 0.006155303
## hipertensionCon hipertensión 2.3141215 0.9795739   5.4668239 0.055754792

En los modelos que se mostraron anteriormente, el primer modelo corresponde al modelo de regresión en el que se emplearon variables demográficas para su ajuste, posteriormente se tomaron unicamente las dos variables de hipertensión y fentanilo, en el último modelo se observa el modelo de regresión de Poisson ajustado por variables demográficas ajustado para dar un RR más robusto utilizando el enfoque sandwich.