source("/Users/mauriciogutierrez/Library/Mobile Documents/com~apple~CloudDocs/estadisticas/general_functions_new.R")
## Loading required package: XLConnectJars
## XLConnect 0.2-11 by Mirai Solutions GmbH [aut],
##   Martin Studer [cre],
##   The Apache Software Foundation [ctb, cph] (Apache POI, Apache Commons
##     Codec),
##   Stephen Colebourne [ctb, cph] (Joda-Time Java library)
## http://www.mirai-solutions.com ,
## http://miraisolutions.wordpress.com
## Loading required package: reshape2
## Loading required package: stringr
library("openxlsx")
## 
## Attaching package: 'openxlsx'
## 
## The following objects are masked from 'package:XLConnect':
## 
##     loadWorkbook, mergeCells, saveWorkbook
todos <- read.xlsx("/Users/mauriciogutierrez/Library/Mobile Documents/com~apple~CloudDocs/estadisticas/mercepin/dm2-2.xlsx", sheet = 1, startRow = 1, colNames = TRUE)

labelsexo=c("F","M")

labeldicotomico=c("Si","No")

labeltipodm=c("DM1","DM2","Farmacologica","Post")
# Demografico ---------------------------------------------------

freq(todos,todos$sexo,using_na=T,label=labelsexo)

## 
##       absolute relative cumulative   ici   ics
## F           23    44.23      44.23 30.47 58.67
## M           29    55.77        100 41.33 69.53
## Total       52   100.00
means(todos,todos$edad)
## 
##                                           edad
## n                                           52
## Mean                                     58.75
## Median                                      58
## Sd                                        14.6
## Min                                         24
## Max                                         88
## Range                                       64
## Q1                                          50
## Q3                                       69.25
## Iqr                                      19.25
## Normality                                0.851
## ICI                                      58.62
## SCI                                      58.88
## Summary   Mean[ICI-ICS] --> 58.75(58.62-58.88)

# Gestion Hospitalaria ---------------------------------------------------

freq(todos,todos$sala,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## M           21    40.38      40.38 27.01  54.9
## Q           31    59.62        100  45.1 72.99
## Total       52   100.00
means(todos,todos$demorainterconsulta)
## 
##                demorainterconsulta
## n                               52
## Mean                          5.81
## Median                           4
## Sd                            6.38
## Min                              0
## Max                             37
## Range                           37
## Q1                               2
## Q3                               7
## Iqr                              5
## Normality                        0
## ICI                           5.75
## SCI                           5.86
## Summary   Median[Q1-Q3]  -->4(2-7)

means(todos,todos$demoravaloracion)
## 
##                      demoravaloracion
## n                                  52
## Mean                             1.79
## Median                              1
## Sd                               1.91
## Min                                 0
## Max                                 9
## Range                               9
## Q1                               0.75
## Q3                                  3
## Iqr                              2.25
## Normality                           0
## ICI                              1.77
## SCI                              1.81
## Summary   Median[Q1-Q3]  -->1(0.75-3)

means(todos,todos$diasinternacion)
## 
##                         diasinternacion
## n                                    52
## Mean                               26.5
## Median                               21
## Sd                                22.57
## Min                                   4
## Max                                 126
## Range                               122
## Q1                                   13
## Q3                                 29.5
## Iqr                                16.5
## Normality                             0
## ICI                                26.3
## SCI                                26.7
## Summary   Median[Q1-Q3]  -->21(13-29.5)

# Caracterizacion Diabetes  ---------------------------------------------------

freq(todos,todos$dmconocido,using_na=T,label=labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si          46    88.46      88.46 76.56 95.65
## No           6    11.54        100  4.35 23.44
## Total       52   100.00
freq(todos,todos$tipodm,using_na=T,label=labeltipodm)

## 
##               absolute relative cumulative   ici   ics
## DM1                  3     5.77       5.77  1.21 15.95
## DM2                 46    88.46      94.23 76.56 95.65
## Farmacologica        2     3.85      98.08  0.47 13.21
## Post                 1     1.92        100  0.05 10.26
## Total               52   100.00
diabeticosconocidos=todos[which(todos$dmconocido==1),]

diabeticosconocidos$tiempoevolucion=as.numeric(diabeticosconocidos$tiempoevolucion)
## Warning: NAs introduced by coercion
means(diabeticosconocidos,diabeticosconocidos$tiempoevolucion)
## 
##                     tiempoevolucion
## n                                45
## Mean                          11.84
## Median                            8
## Sd                            10.79
## Min                               1
## Max                              50
## Range                            49
## Q1                                4
## Q3                               20
## Iqr                              16
## Normality                         0
## ICI                           11.74
## SCI                           11.95
## Summary   Median[Q1-Q3]  -->8(4-20)

# Tratamiento previo ---------------------------------------------------

freq(todos,todos$adoprevio,c("Si","No"))

## 
##         absolute relative cumulative   ici   ics
## Si            26       50         50 35.81 64.19
## No            26       50        100 35.81 64.19
## Missing        0        0     100.00     0  6.85
## Total         52   100.00
freq(diabeticosconocidos,diabeticosconocidos$cumpledieta,using_na=T,labels=c("Si","No","NE"))

## 
##       absolute relative cumulative   ici   ics
## Si           4      8.7        8.7  2.42 20.79
## No          13    28.26      36.96 15.99 43.46
## NE          29    63.04        100 47.55 76.79
## Total       46   100.00
freq(diabeticosconocidos,diabeticosconocidos$cumpleejercicio,using_na=T,c("Si","No","NE"))

## 
##       absolute relative cumulative   ici   ics
## Si           1     2.17       2.17  0.06 11.53
## No          11    23.91      26.09 12.59 38.77
## NE          34    73.91        100 58.87 85.73
## Total       46   100.00
freq(diabeticosconocidos,diabeticosconocidos$metformina,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          24    52.17      52.17 36.95 67.11
## No          22    47.83        100 32.89 63.05
## Total       46   100.00
freq(diabeticosconocidos,diabeticosconocidos$sulfonilureas,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          10    21.74      21.74 10.95 36.36
## No          36    78.26        100 63.64 89.05
## Total       46   100.00
freq(diabeticosconocidos,diabeticosconocidos$glitazonas,using_na=T,c("No"))

## 
##       absolute relative cumulative   ici ics
## No          46      100        100 92.29 100
## Total       46   100.00
freq(diabeticosconocidos,diabeticosconocidos$inhibidoresdpp4,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si           2     4.35       4.35  0.53 14.84
## No          44    95.65        100 85.16 99.47
## Total       46   100.00
freq(diabeticosconocidos,diabeticosconocidos$insulina,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          24    52.17      52.17 36.95 67.11
## No          22    47.83        100 32.89 63.05
## Total       46   100.00
insulinizados=diabeticosconocidos[which(diabeticosconocidos$insulina==1),]

freq(insulinizados,insulinizados$insulinabasal,using_na=T,c("Si","N/E"))

## 
##       absolute relative cumulative   ici   ics
## Si          23    95.83      95.83 78.88 99.89
## N/E          1     4.17        100  0.11 21.12
## Total       24   100.00
freq(insulinizados,insulinizados$tipoinsulinabasal,using_na=T,c("NPH","Analogos","N/E"))

## 
##          absolute relative cumulative   ici   ics
## NPH            21     87.5       87.5 67.64 97.34
## Analogos        1     4.17      91.67  0.11 21.12
## N/E             2     8.33        100  1.03    27
## Total          24   100.00
freq(insulinizados,insulinizados$numerodosisinsulinabasal,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## 1            3     12.5       12.5  2.66 32.36
## 2            6       25       37.5  9.77 46.71
## 3           13    54.17      91.67 32.82 74.45
## 4            1     4.17      95.83  0.11 21.12
## NE           1     4.17        100  0.11 21.12
## Total       24   100.00
means(insulinizados,as.numeric(insulinizados$dosistotalinsulinabasal))
## Warning in descriptive_means(x, varname = namex): NAs introduced by
## coercion
## 
##                       dosistotalinsulinabasal)
## n                                           22
## Mean                                     58.14
## Median                                      54
## Sd                                       34.66
## Min                                          5
## Max                                        120
## Range                                      115
## Q1                                          28
## Q3                                          86
## Iqr                                         58
## Normality                                0.196
## ICI                                      57.67
## SCI                                      58.61
## Summary   Mean[ICI-ICS] --> 58.14(57.67-58.61)

freq(insulinizados,insulinizados$insulinaprandial,c("Si","No","NE"),using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          10    41.67      41.67 22.11 63.36
## No          12       50      91.67 29.12 70.88
## NE           2     8.33        100  1.03    27
## Total       24   100.00
coninsulinaprandial=diabeticosconocidos[which(diabeticosconocidos$insulinaprandial==1),]

freq(coninsulinaprandial,coninsulinaprandial$tipoinsulinaprandial,c("Cristalina","Analogos"),using_na=T)

## 
##            absolute relative cumulative   ici   ics
## Cristalina        7       70         70 34.75 93.33
## Analogos          3       30        100  6.67 65.25
## Total            10   100.00
freq(coninsulinaprandial,coninsulinaprandial$insulinaprandialpor,using_na=T,c("x valor","x Valor/Consumo"))

## 
##                 absolute relative cumulative   ici   ics
## x valor                6       60         60 26.24 87.84
## x Valor/Consumo        4       40        100 12.16 73.76
## Total                 10   100.00
freq(coninsulinaprandial,coninsulinaprandial$numerodosisinsulinaprandial)

## 
##         absolute relative cumulative   ici   ics
## 0              7       70         70 34.75 93.33
## 3              1       10         80  0.25  44.5
## 4              2       20        100  2.52 55.61
## Missing        0        0     100.00     0 30.85
## Total         10   100.00
diabeticosconocidos$tratamiento=""
diabeticosconocidos$tratamiento=paste(diabeticosconocidos$tratamiento,ifelse(diabeticosconocidos$metformina==1,"MET ",""),sep="")
diabeticosconocidos$tratamiento=paste(diabeticosconocidos$tratamiento,ifelse(diabeticosconocidos$sulfonilureas==1,"SUL ",""),sep="")
diabeticosconocidos$tratamiento=paste(diabeticosconocidos$tratamiento,ifelse(diabeticosconocidos$inhibidoresdpp4==1,"DPP ",""),sep="")
diabeticosconocidos$tratamiento=paste(diabeticosconocidos$tratamiento,ifelse(diabeticosconocidos$insulina==1,"INS ",""),sep="")
diabeticosconocidos$tratamiento=ifelse(diabeticosconocidos$tratamiento=="","No Trat. farmacologico",diabeticosconocidos$tratamiento)

freq(diabeticosconocidos,diabeticosconocidos$tratamiento,using_na=T)

## 
##                        absolute relative cumulative   ici   ics
## DPP INS                       1     2.17       2.17  0.06 11.53
## INS                          17    36.96      39.13 23.21 52.45
## MET                           8    17.39      56.52  7.82 31.42
## MET DPP                       1     2.17       58.7  0.06 11.53
## MET INS                       6    13.04      71.74  4.94 26.26
## MET SUL                       9    19.57       91.3  9.36 33.91
## No Trat. farmacologico        3     6.52      97.83  1.37  17.9
## SUL                           1     2.17        100  0.06 11.53
## Total                        46   100.00
# Control metabolico previo  ---------------------------------------------------

freq(todos,todos$tienehba1c36mesesprevios,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si           8    15.38      15.38  6.88 28.08
## No          44    84.62        100 71.92 93.12
## Total       52   100.00
means(todos,todos$valorhba1c36mesesprevios)
## 
##                      valorhba1c36mesesprevios
## n                                           9
## Mean                                    10.19
## Median                                    8.4
## Sd                                       4.03
## Min                                       6.3
## Max                                      18.3
## Range                                      12
## Q1                                        7.4
## Q3                                       11.8
## Iqr                                       4.4
## Normality                               0.134
## ICI                                      10.1
## SCI                                     10.28
## Summary   Mean[ICI-ICS] --> 10.19(10.1-10.28)

means(todos,todos$glucemia6mesesprevios)
## 
##                         glucemia6mesesprevios
## n                                          14
## Mean                                   247.93
## Median                                  178.5
## Sd                                      187.5
## Min                                        62
## Max                                       725
## Range                                     663
## Q1                                     130.75
## Q3                                        272
## Iqr                                    141.25
## Normality                               0.017
## ICI                                    244.73
## SCI                                    251.13
## Summary   Median[Q1-Q3]  -->178.5(130.75-272)

# Complicaciones Cronicas  ---------------------------------------------------

freq(todos,todos$repercusionmicro,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si          28    53.85      53.85 39.47 67.77
## No          24    46.15        100 32.23 60.53
## Total       52   100.00
freq(todos,todos$retinopatia,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          12    23.08      23.08 12.53 36.84
## No          40    76.92        100 63.16 87.47
## Total       52   100.00
freq(todos,todos$nefropatia,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          21    40.38      40.38 27.01  54.9
## No          31    59.62        100  45.1 72.99
## Total       52   100.00
freq(todos,todos$insuficienciarenal,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si          21    40.38      40.38 27.01  54.9
## No          31    59.62        100  45.1 72.99
## Total       52   100.00
freq(todos,todos$neuropatia,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si          23     45.1       45.1 31.13 59.66
## No          28     54.9        100 40.34 68.87
## Total       51   100.00
freq(todos,todos$neuropatiaautonomica,using_na=T,c("Si","No","NE"))

## 
##       absolute relative cumulative   ici   ics
## Si           1     1.92       1.92  0.05 10.26
## No          50    96.15      98.08 86.79 99.53
## NE           1     1.92        100  0.05 10.26
## Total       52   100.00
freq(todos,todos$cardiopatiaisquemica,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si          14    26.92      26.92 15.57 41.02
## No          38    73.08        100 58.98 84.43
## Total       52   100.00
freq(todos,todos$acvoait,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si           6    11.54      11.54  4.35 23.44
## No          46    88.46        100 76.56 95.65
## Total       52   100.00
freq(todos,todos$aocmmii,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si           6    11.54      11.54  4.35 23.44
## No          46    88.46        100 76.56 95.65
## Total       52   100.00
freq(todos,todos$lesionespie,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si          11    21.15      21.15 11.06  34.7
## No          41    78.85        100  65.3 88.94
## Total       52   100.00
freq(todos,todos$lesionespie,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si          11    21.15      21.15 11.06  34.7
## No          41    78.85        100  65.3 88.94
## Total       52   100.00
freq(todos,todos$amputacion,using_na=T,labeldicotomico)

## 
##       absolute relative cumulative   ici   ics
## Si           4     7.69       7.69  2.14 18.54
## No          48    92.31        100 81.46 97.86
## Total       52   100.00
freq(todos,todos$obesidad,using_na=T,c("Si","No","NE"))

## 
##       absolute relative cumulative   ici   ics
## Si          14    26.92      26.92 15.57 41.02
## No          11    21.15      48.08 11.06  34.7
## NE          27    51.92        100 37.63 65.99
## Total       52   100.00
obesos=todos[which(todos$obesidad==1),]
freq(obesos,obesos$gradoobesidad,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## 2            1     7.14       7.14  0.18 33.87
## 3            5    35.71      42.86 12.76 64.86
## 4            8    57.14        100 28.86 82.34
## Total       14   100.00
freq(todos,todos$hta,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          45    86.54      86.54 74.21 94.41
## No           7    13.46        100  5.59 25.79
## Total       52   100.00
freq(todos,todos$dislipidemia,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          22    42.31      42.31 28.73  56.8
## No          30    57.69        100  43.2 71.27
## Total       52   100.00
freq(todos,todos$tabaquista,using_na=T,c("Si","No","Ex"))

## 
##       absolute relative cumulative   ici   ics
## Si           9    17.31      17.31  8.23 30.33
## No          25    48.08      65.38 34.01 62.37
## Ex          18    34.62        100 21.97 49.09
## Total       52   100.00
freq(todos,todos$etlista,c("Si","No","Ex"),using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si           5     9.62       9.62   3.2 21.03
## No          46    88.46      98.08 76.56 95.65
## Ex           1     1.92        100  0.05 10.26
## Total       52   100.00
# Tratamiento en la internacion  ---------------------------------------------------

freq(todos,todos$adoduranteinternacion,using_na=T,c("Si","No"))

## 
##       absolute relative cumulative   ici   ics
## Si          19    36.54      36.54 23.62 51.04
## No          33    63.46        100 48.96 76.38
## Total       52   100.00
adoduranteinternacion=todos[which(todos$adoduranteinternacion==1),]
freq(adoduranteinternacion,adoduranteinternacion$metformina_int,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          18    94.74      94.74 73.97 99.87
## No           1     5.26        100  0.13 26.03
## Total       19   100.00
freq(adoduranteinternacion,adoduranteinternacion$sulfonilureas_int,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si           3    15.79      15.79  3.38 39.58
## No          16    84.21        100 60.42 96.62
## Total       19   100.00
freq(adoduranteinternacion,adoduranteinternacion$glitazonas_int,c("No"),using_na=T)

## 
##       absolute relative cumulative   ici ics
## No          19      100        100 82.35 100
## Total       19   100.00
freq(adoduranteinternacion,adoduranteinternacion$inhibidoresdpp4_int,labeldicotomico,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si           1     5.26       5.26  0.13 26.03
## No          18    94.74        100 73.97 99.87
## Total       19   100.00
insulinainternacion=todos[which(todos$insulinabasal_int==1  | todos$insulinaprandial_int==1),]

freq(insulinainternacion,insulinainternacion$insulinabasal_int,c("Si","No"),using_na=T)

## 
##       absolute relative cumulative   ici   ics
## Si          38    73.08      73.08 58.98 84.43
## No          14    26.92        100 15.57 41.02
## Total       52   100.00
freq(insulinainternacion,insulinainternacion$insulinaprandial_int,c("Si"),using_na=T)

## 
##       absolute relative cumulative   ici ics
## Si          52      100        100 93.15 100
## Total       52   100.00
freq(insulinainternacion,insulinainternacion$tipoinsulinaprandial_int,c("Cristalina","Analogos"),using_na=T)

## 
##            absolute relative cumulative   ici   ics
## Cristalina       49    94.23      94.23 84.05 98.79
## Analogos          3     5.77        100  1.21 15.95
## Total            52   100.00
freq(insulinainternacion,insulinainternacion$insulinaprandiapor_int,c("x valor","x valor/consumo"),using_na=T)

## 
##                 absolute relative cumulative   ici   ics
## x valor               41    78.85      78.85  65.3 88.94
## x valor/consumo       11    21.15        100 11.06  34.7
## Total                 52   100.00
freq(insulinainternacion,insulinainternacion$numerodosisinsulinaprandial_int,using_na=T)

## 
##       absolute relative cumulative  ici   ics
## 0           41    78.85      78.85 65.3 88.94
## 1            3     5.77      84.62 1.21 15.95
## 3            5     9.62      94.23  3.2 21.03
## 4            3     5.77        100 1.21 15.95
## Total       52   100.00
insulinabasalinternacion=todos[which(todos$insulinabasal_int==1),]
freq(insulinabasalinternacion,insulinabasalinternacion$tipoinsulinabasal_int,using_na=T,c("NPH","Analogos"))

## 
##          absolute relative cumulative   ici   ics
## NPH            37    97.37      97.37 86.19 99.93
## Analogos        1     2.63        100  0.07 13.81
## Total          38   100.00
freq(insulinabasalinternacion,insulinabasalinternacion$numeroinsulinabasal_int,using_na=T)

## 
##       absolute relative cumulative   ici   ics
## 1            8    21.05      21.05  9.55 37.32
## 2           11    28.95         50 15.42  45.9
## 3           19       50        100 33.38 66.62
## Total       38   100.00
means(insulinabasalinternacion,insulinabasalinternacion$dosisinsulinabasal_int)
## 
##                  dosisinsulinabasal_int
## n                                    38
## Mean                              37.61
## Median                               37
## Sd                                26.52
## Min                                   5
## Max                                 120
## Range                               115
## Q1                                 18.5
## Q3                                   48
## Iqr                                29.5
## Normality                         0.004
## ICI                               37.33
## SCI                               37.88
## Summary   Median[Q1-Q3]  -->37(18.5-48)

means(todos,as.numeric(todos$glucemiaalingreso))
## 
##                          glucemiaalingreso)
## n                                        38
## Mean                                 281.82
## Median                                  196
## Sd                                   230.45
## Min                                      46
## Max                                    1181
## Range                                  1135
## Q1                                      134
## Q3                                   331.25
## Iqr                                  197.25
## Normality                                 0
## ICI                                  279.46
## SCI                                  284.18
## Summary   Median[Q1-Q3]  -->196(134-331.25)

# Control metabolico en la internacion  ---------------------------------------------------

means(todos,todos$medianaglucemiainternacion)
## 
##                     medianaglucemiainternacion
## n                                           42
## Mean                                    217.13
## Median                                     191
## Sd                                      110.91
## Min                                         69
## Max                                        519
## Range                                      450
## Q1                                      133.25
## Q3                                      288.12
## Iqr                                     154.88
## Normality                                0.004
## ICI                                     216.05
## SCI                                     218.21
## Summary   Median[Q1-Q3]  -->191(133.25-288.12)

means(todos,todos$medianahgtinternacion)
## 
##                            medianahgtinternacion
## n                                             52
## Mean                                       178.6
## Median                                     172.5
## Sd                                         52.09
## Min                                           95
## Max                                          300
## Range                                        205
## Q1                                        145.62
## Q3                                        200.12
## Iqr                                         54.5
## Normality                                  0.011
## ICI                                       178.14
## SCI                                       179.05
## Summary   Median[Q1-Q3]  -->172.5(145.62-200.12)

means(todos,todos$medianahgtinternacion)
## 
##                            medianahgtinternacion
## n                                             52
## Mean                                       178.6
## Median                                     172.5
## Sd                                         52.09
## Min                                           95
## Max                                          300
## Range                                        205
## Q1                                        145.62
## Q3                                        200.12
## Iqr                                         54.5
## Normality                                  0.011
## ICI                                       178.14
## SCI                                       179.05
## Summary   Median[Q1-Q3]  -->172.5(145.62-200.12)

means(todos,todos$rihgtinternacion)
## 
##                               rihgtinternacion
## n                                           52
## Mean                                     81.28
## Median                                    78.5
## Sd                                       38.75
## Min                                         13
## Max                                        184
## Range                                      171
## Q1                                       51.19
## Q3                                      101.56
## Iqr                                      50.38
## Normality                                0.311
## ICI                                      80.95
## SCI                                      81.62
## Summary   Mean[ICI-ICS] --> 81.28(80.95-81.62)

means(todos,as.numeric(todos$porcentajecon3omashgtdia))
## 
##                    porcentajecon3omashgtdia)
## n                                         52
## Mean                                   70.56
## Median                                  76.6
## Sd                                     24.26
## Min                                        0
## Max                                      100
## Range                                    100
## Q1                                      58.3
## Q3                                     88.97
## Iqr                                    30.67
## Normality                                  0
## ICI                                    70.35
## SCI                                    70.77
## Summary   Median[Q1-Q3]  -->76.6(58.3-88.97)

means(todos,todos$porcentajehgtenobjetivo)
## 
##                       porcentajehgtenobjetivo
## n                                          52
## Mean                                    23.71
## Median                                   16.9
## Sd                                      20.31
## Min                                         0
## Max                                        90
## Range                                      90
## Q1                                      10.23
## Q3                                      34.95
## Iqr                                     24.73
## Normality                                   0
## ICI                                     23.53
## SCI                                     23.89
## Summary   Median[Q1-Q3]  -->16.9(10.23-34.95)

means(todos,todos$hba1c3mesespreviosointernacion)
## 
##           hba1c3mesespreviosointernacion
## n                                     11
## Mean                                 9.3
## Median                                 8
## Sd                                  3.43
## Min                                  6.3
## Max                                 18.3
## Range                                 12
## Q1                                   7.4
## Q3                                 10.45
## Iqr                                 3.05
## Normality                          0.003
## ICI                                 9.23
## SCI                                 9.37
## Summary   Median[Q1-Q3]  -->8(7.4-10.45)

means(todos,todos$porcentajedehipoglucemiasdurantelainternacion)
## 
##           porcentajedehipoglucemiasdurantelainternacion
## n                                                    52
## Mean                                               3.38
## Median                                             0.95
## Sd                                                  5.6
## Min                                                   0
## Max                                                23.6
## Range                                              23.6
## Q1                                                    0
## Q3                                                    5
## Iqr                                                   5
## Normality                                             0
## ICI                                                3.34
## SCI                                                3.43
## Summary                     Median[Q1-Q3]  -->0.95(0-5)

means(todos,todos$porcentajedehiperglucemiasduranteinternacion)
## 
##           porcentajedehiperglucemiasduranteinternacion
## n                                                   52
## Mean                                             42.51
## Median                                           44.05
## Sd                                               27.46
## Min                                                  0
## Max                                               98.5
## Range                                             98.5
## Q1                                               21.17
## Q3                                               60.12
## Iqr                                              38.95
## Normality                                        0.143
## ICI                                              42.27
## SCI                                              42.75
## Summary           Mean[ICI-ICS] --> 42.51(42.27-42.75)

# Ingreso / Complicaciones / Derivacion  ---------------------------------------------------

freq(todos,todos$resultadosinternacion)

## 
##         absolute relative cumulative   ici   ics
## 1             51    98.08      98.08 89.74 99.95
## 2              1     1.92        100  0.05 10.26
## Missing        0        0     100.00     0  6.85
## Total         52   100.00
freq(todos,todos$derivacion,c("CTI","No"))

## 
##         absolute relative cumulative   ici   ics
## CTI            1     1.92       1.92  0.05 10.26
## No            27    51.92      53.85 37.63 65.99
## Missing       24    46.15     100.00 32.23 60.53
## Total         28   100.00
freq(todos,todos$motivoingreso,using_na=T)

## 
##       absolute relative cumulative  ici   ics
## 1            6    11.54      11.54 4.35 23.44
## 2            6    11.54      23.08 4.35 23.44
## 3            1     1.92         25 0.05 10.26
## 4            4     7.69      32.69 2.14 18.54
## 5            3     5.77      38.46 1.21 15.95
## 6            9    17.31      55.77 8.23 30.33
## 7            2     3.85      59.62 0.47 13.21
## 8            9    17.31      76.92 8.23 30.33
## 9            1     1.92      78.85 0.05 10.26
## 11           4     7.69      86.54 2.14 18.54
## 12           3     5.77      92.31 1.21 15.95
## 13           2     3.85      96.15 0.47 13.21
## 14           2     3.85        100 0.47 13.21
## Total       52   100.00
freq(todos,todos$complicacionesinternacion)

## 
##              absolute relative cumulative   ici   ics
## -1                 29    55.77      55.77 41.33 69.53
## 1, 6, 7             1     1.92      57.69  0.05 10.26
## 11                  1     1.92      59.62  0.05 10.26
## 17                  1     1.92      61.54  0.05 10.26
## 18                  1     1.92      63.46  0.05 10.26
## 3                   1     1.92      65.38  0.05 10.26
## 3, 10               1     1.92      67.31  0.05 10.26
## 3, 6                1     1.92      69.23  0.05 10.26
## 3, 9                1     1.92      71.15  0.05 10.26
## 5, 6, 8, 12         1     1.92      73.08  0.05 10.26
## 6                   5     9.62      82.69   3.2 21.03
## 6, 7, 9             1     1.92      84.62  0.05 10.26
## 7                   2     3.85      88.46  0.47 13.21
## 7, 15               1     1.92      90.38  0.05 10.26
## 7, 16               1     1.92      92.31  0.05 10.26
## 7, 8                1     1.92      94.23  0.05 10.26
## 7, 8, 11            1     1.92      96.15  0.05 10.26
## 8                   2     3.85        100  0.47 13.21
## Missing             0        0     100.00     0  6.85
## Total              52   100.00
todos$complicacion=ifelse(todos$complicacionesinternacion==-1,0,1)

freq(todos,todos$complicacion,c("No","Si"),using_na=T)

## 
##       absolute relative cumulative   ici   ics
## No          29    55.77      55.77 41.33 69.53
## Si          23    44.23        100 30.47 58.67
## Total       52   100.00
# Vinculacion de complicaciones con confundentes --------------------------------------

means(todos,todos$edad,todos$complicacion,c("No","Si"))
## 
##               complicacion = No  complicacion = Si
## n                            29                 23
## Mean                      61.79              54.91
## Median                       61                 56
## Sd                        12.57               16.3
## Min                          40                 24
## Max                          88                 88
## Range                        48                 64
## Q1                           52                 43
## Q3                           72                 67
## Iqr                          20                 24
## Normal                    0.799              0.924
## IIC                       61.65               54.7
## ICS                       61.94              55.13
## Summary2    61.79 (61.65-61.94) 54.91 (54.7-55.13)
## Summary           Mean[ICI-ICS]                   
## Groups            Homocedastics                   
## Normality                   Yes                   
## Differences                  No                   
## p-value                   0.092                   
## Method        t-Student(Pooled)

tables(todos,todos$sexo,todos$complicacion,labelsexo,c("No","Si"))
## 
##      complicacion
## sexo  No Si Sum
##   F   13 10  23
##   M   16 13  29
##   Sum 29 23  52
## 
## by Row (%)
##      No    Si
## F 56.52 43.47
## M 55.17 44.82
##   55.76 44.23
## 
## by Col (%)
##      No    Si      
## F 44.82 43.47 44.23
## M 55.17 56.52 55.76
## 
## Chi2 Pearson  p = 0.922

tables(todos,todos$sala,todos$complicacion,c("M","Q"),c("No","Si"))
## 
##      complicacion
## sala  No Si Sum
##   M   10 11  21
##   Q   19 12  31
##   Sum 29 23  52
## 
## by Row (%)
##      No    Si
## M 47.61 52.38
## Q 61.29 38.70
##   55.76 44.23
## 
## by Col (%)
##      No    Si      
## M 34.48 47.82 40.38
## Q 65.51 52.17 59.61
## 
## Chi2 Pearson  p = 0.33

# Vinculacion de complicaciones tratamiento --------------------------------------

tables(todos,todos$insulina,todos$complicacion,c("Si","No"),c("No","Si"))
## 
##         complicacion
## insulina No Si Sum
##      Si  11 13  24
##      No  18 10  28
##      Sum 29 23  52
## 
## by Row (%)
##       No    Si
## Si 45.83 54.16
## No 64.28 35.71
##    55.76 44.23
## 
## by Col (%)
##       No    Si      
## Si 37.93 56.52 46.15
## No 62.06 43.47 53.84
## 
## Chi2 Pearson  p = 0.182

tables(todos,todos$adoduranteinternacion,todos$complicacion,c("Si","No"),c("No","Si"))
## 
##                      complicacion
## adoduranteinternacion No Si Sum
##                   Si  13  6  19
##                   No  16 17  33
##                   Sum 29 23  52
## 
## by Row (%)
##       No    Si
## Si 68.42 31.57
## No 48.48 51.51
##    55.76 44.23
## 
## by Col (%)
##       No    Si      
## Si 44.82 26.08 36.53
## No 55.17 73.91 63.46
## 
## Chi2 Pearson  p = 0.163

# Vinculacion de complicaciones con control metabolico --------------------------------------

means(todos,todos$hba1c3mesespreviosointernacion,todos$complicacion,c("No","Si"))
## 
##             complicacion = No complicacion = Si
## n                           3                 8
## Mean                     7.07             10.14
## Median                    6.8               8.8
## Sd                       0.93              3.69
## Min                       6.3               7.3
## Max                       8.1              18.3
## Range                     1.8                11
## Q1                       6.55              7.65
## Q3                       7.45             11.33
## Iqr                       0.9              3.68
## Normal                   0.52             0.016
## IIC                      7.03             10.05
## ICS                       7.1             10.22
## Summary2      6.8 (6.55-7.45)  8.8 (7.65-11.33)
## Summary         Median[Q1-Q3]                  
## Groups          Homocedastics                  
## Normality                  No                  
## Differences                No                  
## p-value                 0.133                  
## Method               Wilcoxon

means(todos,todos$porcentajehgtenobjetivo,todos$complicacion,c("No","Si"))
## 
##             complicacion = No complicacion = Si
## n                          29                23
## Mean                    24.64             22.54
## Median                   16.4              18.5
## Sd                      23.24             16.31
## Min                         0              2.86
## Max                        90              67.3
## Range                      90             64.44
## Q1                        9.7              14.1
## Q3                       36.3              28.9
## Iqr                      26.6              14.8
## Normal                  0.002              0.02
## IIC                     24.36             22.32
## ICS                     24.91             22.75
## Summary2      16.4 (9.7-36.3)  18.5 (14.1-28.9)
## Summary         Median[Q1-Q3]                  
## Groups          Homocedastics                  
## Normality                  No                  
## Differences                No                  
## p-value                 0.775                  
## Method         Kruskal-Wallis

means(todos,todos$porcentajedehipoglucemiasdurantelainternacion,todos$complicacion,c("No","Si"))
## 
##             complicacion = No complicacion = Si
## n                          29                23
## Mean                     3.42              3.34
## Median                      0               1.4
## Sd                       6.39              4.57
## Min                         0                 0
## Max                      23.6              16.7
## Range                    23.6              16.7
## Q1                          0                 0
## Q3                        5.3              4.25
## Iqr                       5.3              4.25
## Normal                      0                 0
## IIC                      3.35              3.28
## ICS                       3.5               3.4
## Summary2            0 (0-5.3)      1.4 (0-4.25)
## Summary         Median[Q1-Q3]                  
## Groups          Homocedastics                  
## Normality                  No                  
## Differences                No                  
## p-value                 0.211                  
## Method         Kruskal-Wallis

means(todos,todos$porcentajedehiperglucemiasduranteinternacion,todos$complicacion,c("No","Si"))
## 
##              complicacion = No  complicacion = Si
## n                           29                 23
## Mean                      40.7               44.8
## Median                    38.2                 45
## Sd                       31.33              22.11
## Min                          0                8.2
## Max                       98.5               85.7
## Range                     98.5               77.5
## Q1                        11.5               28.8
## Q3                        59.6              61.75
## Iqr                       48.1              32.95
## Normal                   0.072              0.608
## IIC                      40.33              44.51
## ICS                      41.07              45.09
## Summary2    40.7 (40.33-41.07) 44.8 (44.51-45.09)
## Summary          Mean[ICI-ICS]                   
## Groups           Homocedastics                   
## Normality                  Yes                   
## Differences                 No                   
## p-value                  0.598                   
## Method       t-Student(Pooled)

means(todos,todos$rihgtinternacion,todos$complicacion,c("No","Si"))
## 
##               complicacion = No  complicacion = Si
## n                            29                 23
## Mean                      67.77              98.33
## Median                    66.75                 91
## Sd                        32.79              39.59
## Min                          13                 30
## Max                         139                184
## Range                       126                154
## Q1                         39.5               74.5
## Q3                           81                121
## Iqr                        41.5               46.5
## Normal                    0.635              0.748
## IIC                       67.38               97.8
## ICS                       68.15              98.85
## Summary2    67.77 (67.38-68.15) 98.33 (97.8-98.85)
## Summary           Mean[ICI-ICS]                   
## Groups            Homocedastics                   
## Normality                   Yes                   
## Differences                 Yes                   
## p-value                   0.004                   
## Method        t-Student(Pooled)

means(todos,todos$medianahgtinternacion,todos$complicacion,c("No","Si"))
## 
##             complicacion = No  complicacion = Si
## n                          29                 23
## Mean                   179.57             177.37
## Median                    173                172
## Sd                      61.08              39.23
## Min                        95              122.5
## Max                       300                263
## Range                     205              140.5
## Q1                        130             147.75
## Q3                        211              193.5
## Iqr                        81              45.75
## Normal                  0.042              0.027
## IIC                    178.85             176.85
## ICS                    180.29             177.89
## Summary2        173 (130-211) 172 (147.75-193.5)
## Summary         Median[Q1-Q3]                   
## Groups        Heteroscedastic                   
## Normality                  No                   
## Differences                No                   
## p-value                 0.912                   
## Method         Kruskal-Wallis

# Vinculacion de complicaciones con gestion hospitalaria --------------------------------------

means(todos,todos$diasinternacion,todos$complicacion,c("No","Si"))
## 
##             complicacion = No complicacion = Si
## n                          29                23
## Mean                    16.97             38.52
## Median                     14                28
## Sd                      12.48             26.67
## Min                         4                 8
## Max                        67               126
## Range                      63               118
## Q1                          9                21
## Q3                         21              48.5
## Iqr                        12              27.5
## Normal                      0             0.001
## IIC                     16.82             38.17
## ICS                     17.11             38.87
## Summary2            14 (9-21)      28 (21-48.5)
## Summary         Median[Q1-Q3]                  
## Groups        Heteroscedastic                  
## Normality                  No                  
## Differences                No                  
## p-value                <0.001                  
## Method         Kruskal-Wallis

means(todos,todos$demoravaloracion,todos$complicacion,c("No","Si"))
## 
##             complicacion = No complicacion = Si
## n                          29                23
## Mean                     2.07              1.43
## Median                      1                 1
## Sd                       1.85              1.97
## Min                         0                 0
## Max                         9                 7
## Range                       9                 7
## Q1                          1                 0
## Q3                          3                 2
## Iqr                         2                 2
## Normal                      0                 0
## IIC                      2.05              1.41
## ICS                      2.09              1.46
## Summary2              1 (1-3)           1 (0-2)
## Summary         Median[Q1-Q3]                  
## Groups          Homocedastics                  
## Normality                  No                  
## Differences               Yes                  
## p-value                 0.047                  
## Method         Kruskal-Wallis

todos$tiempoevolucion=as.numeric(todos$tiempoevolucion)
## Warning: NAs introduced by coercion
means(todos,todos$tiempoevolucion,todos$complicacion,c("No","Si"))
## 
##             complicacion = No complicacion = Si
## n                          28                23
## Mean                    12.14              8.39
## Median                    7.5                 6
## Sd                      12.46              8.24
## Min                         0                 0
## Max                        50                26
## Range                      50                26
## Q1                          3               1.5
## Q3                      20.75              12.5
## Iqr                     17.75                11
## Normal                  0.001             0.007
## IIC                     11.99              8.28
## ICS                     12.29               8.5
## Summary2        7.5 (3-20.75)      6 (1.5-12.5)
## Summary         Median[Q1-Q3]                  
## Groups          Homocedastics                  
## Normality                  No                  
## Differences                No                  
## p-value                 0.301                  
## Method         Kruskal-Wallis

# Caracterizacion de grupos demografico / Gestion --------------------------------------

means(todos,todos$edad,todos$sexo,labelsexo)
## 
##                      sexo = F            sexo = M
## n                          23                  29
## Mean                     57.7               59.59
## Median                     56                  60
## Sd                      14.76               14.69
## Min                        33                  24
## Max                        85                  88
## Range                      52                  64
## Q1                         51                  50
## Q3                       70.5                  69
## Iqr                      19.5                  19
## Normal                  0.424               0.954
## IIC                      57.5               59.41
## ICS                     57.89               59.76
## Summary2    57.7 (57.5-57.89) 59.59 (59.41-59.76)
## Summary         Mean[ICI-ICS]                    
## Groups          Homocedastics                    
## Normality                 Yes                    
## Differences                No                    
## p-value                 0.647                    
## Method      t-Student(Pooled)

means(todos,todos$edad,todos$sala)
## 
##                      sala = M            sala = Q
## n                          21                  31
## Mean                       54               61.97
## Median                     52                  60
## Sd                       15.2                13.5
## Min                        24                  33
## Max                        77                  88
## Range                      53                  55
## Q1                         41                53.5
## Q3                         68                  71
## Iqr                        27                17.5
## Normal                  0.417               0.749
## IIC                     53.79               61.81
## ICS                     54.21               62.12
## Summary2     54 (53.79-54.21) 61.97 (61.81-62.12)
## Summary         Mean[ICI-ICS]                    
## Groups          Homocedastics                    
## Normality                 Yes                    
## Differences                No                    
## p-value                 0.053                    
## Method      t-Student(Pooled)

tables(todos,todos$sala,todos$sexo,c("M","Q"),labelsexo)
## 
##      sexo
## sala   F  M Sum
##   M   10 11  21
##   Q   13 18  31
##   Sum 23 29  52
## 
## by Row (%)
##       F     M
## M 47.61 52.38
## Q 41.93 58.06
##   44.23 55.76
## 
## by Col (%)
##       F     M      
## M 43.47 37.93 40.38
## Q 56.52 62.06 59.61
## 
## Chi2 Pearson  p = 0.686

means(todos,todos$demorainterconsulta,todos$sala)
## 
##                    sala = M    sala = Q
## n                        21          31
## Mean                   5.62        5.94
## Median                    6           4
## Sd                     4.09        7.62
## Min                       0           1
## Max                      14          37
## Range                    14          36
## Q1                        2         1.5
## Q3                        8         6.5
## Iqr                       6           5
## Normal                0.204           0
## IIC                    5.56        5.85
## ICS                    5.68        6.02
## Summary2            6 (2-8) 4 (1.5-6.5)
## Summary       Median[Q1-Q3]            
## Groups      Heteroscedastic            
## Normality                No            
## Differences              No            
## p-value               0.387            
## Method       Kruskal-Wallis

means(todos,todos$diasinternacion,todos$sala)
## 
##                    sala = M       sala = Q
## n                        21             31
## Mean                  23.67          28.42
## Median                   23             21
## Sd                    16.22           26.1
## Min                       5              4
## Max                      58            126
## Range                    53            122
## Q1                       12           13.5
## Q3                       28           32.5
## Iqr                      16             19
## Normal                0.011              0
## IIC                   23.44          28.12
## ICS                   23.89          28.72
## Summary2         23 (12-28) 21 (13.5-32.5)
## Summary       Median[Q1-Q3]               
## Groups      Heteroscedastic               
## Normality                No               
## Differences              No               
## p-value               0.758               
## Method       Kruskal-Wallis

means(todos,todos$diasinternacion,todos$sexo,labelsexo)
## 
##                    sexo = F   sexo = M
## n                        23         29
## Mean                  29.96      23.76
## Median                   25         19
## Sd                    28.21      16.87
## Min                       4          6
## Max                     126         88
## Range                   122         82
## Q1                      9.5         14
## Q3                       41         25
## Iqr                    31.5         11
## Normal                    0          0
## IIC                   29.58      23.56
## ICS                   30.33      23.96
## Summary2        25 (9.5-41) 19 (14-25)
## Summary       Median[Q1-Q3]           
## Groups      Heteroscedastic           
## Normality                No           
## Differences              No           
## p-value               0.685           
## Method       Kruskal-Wallis

# Caracterizacion demografica / Tratamientos. --------------------------------------

means(todos,todos$edad,todos$insulina)
## 
##                   insulina = 1        insulina = 2
## n                           24                  28
## Mean                     55.71               61.36
## Median                      54                  61
## Sd                        15.1               13.91
## Min                         33                  24
## Max                         85                  88
## Range                       52                  64
## Q1                        45.5               54.25
## Q3                          68                70.5
## Iqr                       22.5               16.25
## Normal                   0.319               0.678
## IIC                      55.51               61.19
## ICS                       55.9               61.52
## Summary2    55.71 (55.51-55.9) 61.36 (61.19-61.52)
## Summary          Mean[ICI-ICS]                    
## Groups           Homocedastics                    
## Normality                  Yes                    
## Differences                 No                    
## p-value                  0.167                    
## Method       t-Student(Pooled)

means(diabeticosconocidos,diabeticosconocidos$tiempoevolucion,diabeticosconocidos$sexo)
## 
##                   sexo = 1      sexo = 2
## n                       19            26
## Mean                 14.05         10.23
## Median                  15           7.5
## Sd                   10.93         10.61
## Min                      1             1
## Max                     34            50
## Range                   33            49
## Q1                       3          4.25
## Q3                    22.5            11
## Iqr                   19.5          6.75
## Normal               0.045             0
## IIC                  13.89          10.1
## ICS                  14.21         10.36
## Summary2       15 (3-22.5) 7.5 (4.25-11)
## Summary      Median[Q1-Q3]              
## Groups       Homocedastics              
## Normality               No              
## Differences             No              
## p-value              0.407              
## Method      Kruskal-Wallis

means(diabeticosconocidos,diabeticosconocidos$tiempoevolucion,diabeticosconocidos$insulina,labeldicotomico)
## 
##              insulina = Si insulina = No
## n                       24            21
## Mean                 13.67          9.76
## Median                12.5             5
## Sd                    9.58         11.92
## Min                      1             1
## Max                     30            50
## Range                   29            49
## Q1                    4.75             3
## Q3                      23            10
## Iqr                  18.25             7
## Normal               0.036             0
## IIC                  13.54           9.6
## ICS                  13.79          9.93
## Summary2    12.5 (4.75-23)      5 (3-10)
## Summary      Median[Q1-Q3]              
## Groups       Homocedastics              
## Normality               No              
## Differences             No              
## p-value              0.094              
## Method      Kruskal-Wallis

means(diabeticosconocidos,diabeticosconocidos$tiempoevolucion,diabeticosconocidos$metformina,labeldicotomico)
## 
##             metformina = Si metformina = No
## n                        23              22
## Mean                  12.13           11.55
## Median                    8               7
## Sd                     11.8            9.89
## Min                       1               1
## Max                      50              30
## Range                    49              29
## Q1                      4.5               3
## Q3                     13.5            21.5
## Iqr                       9            18.5
## Normal                    0           0.007
## IIC                   11.97           11.41
## ICS                   12.29           11.68
## Summary2       8 (4.5-13.5)      7 (3-21.5)
## Summary       Median[Q1-Q3]                
## Groups        Homocedastics                
## Normality                No                
## Differences              No                
## p-value               0.724                
## Method       Kruskal-Wallis