mara.R
mauriciogutierrez
Sun Feb 12 12:24:12 2017
source("/Users/mauriciogutierrez/Library/Mobile Documents/com~apple~CloudDocs/estadisticas/general_functions_new.R")## Loading required package: XLConnectJars## XLConnect 0.2-12 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),
## Graph Builder [ctb, cph] (Curvesapi Java library)## http://www.mirai-solutions.com ,
## http://miraisolutions.wordpress.comdata <- readexcel("/Users/mauriciogutierrez/Library/Mobile Documents/com~apple~CloudDocs/estadisticas/Pilar_Serra/planilla_actualizada.xlsx","Hoja1")
labeldicotomico<-c("No","Si")
freq(data,data$sexo,using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## F 120 70.18 70.18 62.72
## M 51 29.82 100 23.08
## Total 171 100.00
## intervalo confianza inferior
## F 76.92
## M 37.28
## Totalmeans(data,data$edad)##
## edad
## n 171
## Media 61.78
## Mediana 62
## DS 11.31
## Min 29
## Max 89
## Rango 60
## Q1 54
## Q3 70.5
## RIQ 16.5
## Normal 0.048
## ICI 60.08
## ICS 63.49
## Resumen Mediana[Q1-Q3] -->62(54-70.5)
means(data,data$edad,data$sexo)##
## edad by sexo = F edad by sexo = M
## n 120 51
## Media 61.83 61.67
## Mediana 62 59
## DS 12.04 9.48
## Min 29 41
## Max 89 78
## Rango 60 37
## Q1 54.75 53.5
## Q3 70.25 70
## RIQ 15.5 16.5
## Normal 0.044 0.016
## ICI 59.66 59
## ICS 64.01 64.33
## Resumen2 62 (54.75-70.25) 59 (53.5-70)
## Resumen1 Mediana[Q1-Q3]
## grupos Homocedastica
## Normalidad No
## Diferencias No
## p-value 0.625
## Metodo Wilcoxon
freq(data,data$tipo,using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## 2 171 100 100 97.87
## Total 171 100.00
## intervalo confianza inferior
## 2 100
## Totalmeans(data,data$duracion)##
## duracion
## n 171
## Media 8.3
## Mediana 6
## DS 6.93
## Min 1
## Max 35
## Rango 34
## Q1 2
## Q3 12
## RIQ 10
## Normal 0
## ICI 7.25
## ICS 9.34
## Resumen Mediana[Q1-Q3] -->6(2-12)
freq(data,data$hd,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 79 46.2 46.2 38.56
## Si 92 53.8 100 46.03
## Total 171 100.00
## intervalo confianza inferior
## No 53.97
## Si 61.44
## Totalfreq(data,data$ejercicio,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 116 67.84 67.84 60.28
## Si 55 32.16 100 25.24
## Total 171 100.00
## intervalo confianza inferior
## No 74.76
## Si 39.72
## Totalfreq(data,data$basal,using_na=T,c("No","NPH","Glargina"))
##
## absoluta relativa acumulada intervalo confianza inferior
## No 133 77.78 77.78 70.8
## NPH 37 21.64 99.42 15.72
## Glargina 1 0.58 100 0.01
## Total 171 100.00
## intervalo confianza inferior
## No 83.77
## NPH 28.57
## Glargina 3.22
## Totalsolobasal=data[which(data$basal>0),]
freq(solobasal,solobasal$basal,using_na=T,c("NPH","Glargina"))
##
## absoluta relativa acumulada intervalo confianza inferior
## NPH 37 97.37 97.37 86.19
## Glargina 1 2.63 100 0.07
## Total 38 100.00
## intervalo confianza inferior
## NPH 99.93
## Glargina 13.81
## Totalfreq(data,data$prandial,using_na=T,c("No","Cristalina","Glulisina"))
##
## absoluta relativa acumulada intervalo confianza inferior
## No 160 93.57 93.57 88.78
## Cristalina 10 5.85 99.42 2.84
## Glulisina 1 0.58 100 0.01
## Total 171 100.00
## intervalo confianza inferior
## No 96.75
## Cristalina 10.49
## Glulisina 3.22
## Totalsoloprandial=data[which(data$prandial>0),]
freq(soloprandial,soloprandial$prandial,using_na=T,c("Cristalina","Glulisina"))
##
## absoluta relativa acumulada intervalo confianza inferior
## Cristalina 10 90.91 90.91 58.72
## Glulisina 1 9.09 100 0.23
## Total 11 100.00
## intervalo confianza inferior
## Cristalina 99.77
## Glulisina 41.28
## Totalfreq(data,data$metformina,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 36 21.05 21.05 15.2
## Si 135 78.95 100 72.07
## Total 171 100.00
## intervalo confianza inferior
## No 27.93
## Si 84.8
## Totalfreq(data,data$su,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 101 59.06 59.06 51.3
## Si 70 40.94 100 33.49
## Total 171 100.00
## intervalo confianza inferior
## No 66.51
## Si 48.7
## Totalfreq(data,data$metiglinida,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 170 99.42 99.42 96.78
## Si 1 0.58 100 0.01
## Total 171 100.00
## intervalo confianza inferior
## No 99.99
## Si 3.22
## Totalfreq(data,data$tiazolidinedionas,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 166 97.08 97.08 93.31
## Si 5 2.92 100 0.96
## Total 171 100.00
## intervalo confianza inferior
## No 99.04
## Si 6.69
## Totalfreq(data,data$idpp4,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 168 98.25 98.25 94.96
## Si 3 1.75 100 0.36
## Total 171 100.00
## intervalo confianza inferior
## No 99.64
## Si 5.04
## Totaldata$insulinabasal=ifelse(data$basal>0,1,0)
data$insulinaprandial=ifelse(data$prandial>0,1,0)
freq(data,data$insulinabasal,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 133 77.78 77.78 70.8
## Si 38 22.22 100 16.23
## Total 171 100.00
## intervalo confianza inferior
## No 83.77
## Si 29.2
## Totalfreq(data,data$insulinaprandial,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 160 93.57 93.57 88.78
## Si 11 6.43 100 3.25
## Total 171 100.00
## intervalo confianza inferior
## No 96.75
## Si 11.22
## Totaldata$tipoinsulinizacion=data$insulinabasal+data$insulinaprandial*2
freq(data,data$tipoinsulinizacion,using_na=T,c("Sin insulinizar","Solo Basal","Solo Prandial","Basal/Prandial"))
##
## absoluta relativa acumulada intervalo confianza inferior
## Sin insulinizar 132 77.19 77.19 70.17
## Solo Basal 28 16.37 93.57 11.16
## Solo Prandial 1 0.58 94.15 0.01
## Basal/Prandial 10 5.85 100 2.84
## Total 171 100.00
## intervalo confianza inferior
## Sin insulinizar 83.25
## Solo Basal 22.79
## Solo Prandial 3.22
## Basal/Prandial 10.49
## Totaldata$insulina=ifelse(data$tipoinsulinizacion==0,0,1)
freq(data,data$insulina,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 132 77.19 77.19 70.17
## Si 39 22.81 100 16.75
## Total 171 100.00
## intervalo confianza inferior
## No 83.25
## Si 29.83
## Totaldata$ados=data$metformina+data$su+data$metiglinida+data$tiazolidinedionas+data$idpp4
data$tratamiento=data$insulina*1+data$metformina*2+data$su*4+data$metiglinida*8+data$tiazolidinedionas*16+data$idpp4*32
data$recibeado=ifelse(data$ados==0,0,1)
freq(data,data$recibeado,labeldicotomico,using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 19 11.11 11.11 6.82
## Si 152 88.89 100 83.19
## Total 171 100.00
## intervalo confianza inferior
## No 16.81
## Si 93.18
## Totaldata$soloado=ifelse(data$tratamiento%%2==0&data$tratamiento>0,1,0)
freq(data,data$soloado,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 44 25.73 25.73 19.36
## Si 127 74.27 100 67.04
## Total 171 100.00
## intervalo confianza inferior
## No 32.96
## Si 80.64
## Totalsoloado=data[which(data$soloado==1),]
means(soloado,soloado$ados)##
## ados
## n 127
## Media 1.42
## Mediana 1
## DS 0.54
## Min 1
## Max 3
## Rango 2
## Q1 1
## Q3 2
## RIQ 1
## Normal 0
## ICI 1.32
## ICS 1.51
## Resumen Mediana[Q1-Q3] -->1(1-2)
labeltratamientos=c("Solo HD","Solo Insulina","Solo Metformina","Insulina+Metformina","Sulfonilureas","Insulina+Sulfonilureas","Metformina+Sulfonilureas","Insulina+Metformina+Sulfonilureas","Metiglinidas","Insulina+Metformina+tiazolidinedionas","Sulfonilureas+tiazolidinedionas","Metformina+Sulfonilureas+tiazolidinedionas","idpp4+Insulina","idpp4+sulfoniulureas+Metformina")
freq(data,data$tratamiento,using_na=T,labeltratamientos)
##
## absoluta relativa acumulada
## Solo HD 5 2.92 2.92
## Solo Insulina 14 8.19 11.11
## Solo Metformina 64 37.43 48.54
## Insulina+Metformina 14 8.19 56.73
## Sulfonilureas 12 7.02 63.74
## Insulina+Sulfonilureas 1 0.58 64.33
## Metformina+Sulfonilureas 45 26.32 90.64
## Insulina+Metformina+Sulfonilureas 7 4.09 94.74
## Metiglinidas 1 0.58 95.32
## Insulina+Metformina+tiazolidinedionas 2 1.17 96.49
## Sulfonilureas+tiazolidinedionas 2 1.17 97.66
## Metformina+Sulfonilureas+tiazolidinedionas 1 0.58 98.25
## idpp4+Insulina 1 0.58 98.83
## idpp4+sulfoniulureas+Metformina 2 1.17 100
## Total 171 100.00
## intervalo confianza inferior
## Solo HD 0.96
## Solo Insulina 4.55
## Solo Metformina 30.16
## Insulina+Metformina 4.55
## Sulfonilureas 3.68
## Insulina+Sulfonilureas 0.01
## Metformina+Sulfonilureas 19.89
## Insulina+Metformina+Sulfonilureas 1.66
## Metiglinidas 0.01
## Insulina+Metformina+tiazolidinedionas 0.14
## Sulfonilureas+tiazolidinedionas 0.14
## Metformina+Sulfonilureas+tiazolidinedionas 0.01
## idpp4+Insulina 0.01
## idpp4+sulfoniulureas+Metformina 0.14
## Total
## intervalo confianza inferior
## Solo HD 6.69
## Solo Insulina 13.36
## Solo Metformina 45.14
## Insulina+Metformina 13.36
## Sulfonilureas 11.94
## Insulina+Sulfonilureas 3.22
## Metformina+Sulfonilureas 33.58
## Insulina+Metformina+Sulfonilureas 8.25
## Metiglinidas 3.22
## Insulina+Metformina+tiazolidinedionas 4.16
## Sulfonilureas+tiazolidinedionas 4.16
## Metformina+Sulfonilureas+tiazolidinedionas 3.22
## idpp4+Insulina 3.22
## idpp4+sulfoniulureas+Metformina 4.16
## Totalfreq(data,data$ieca,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 67 39.18 39.18 31.82
## Si 104 60.82 100 53.07
## Total 171 100.00
## intervalo confianza inferior
## No 46.93
## Si 68.18
## Totalfreq(data,data$ara2,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 147 85.96 85.96 79.84
## Si 24 14.04 100 9.2
## Total 171 100.00
## intervalo confianza inferior
## No 90.8
## Si 20.16
## Totalfreq(data,data$bloqueoca,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 141 82.46 82.46 75.91
## Si 30 17.54 100 12.16
## Total 171 100.00
## intervalo confianza inferior
## No 87.84
## Si 24.09
## Totalfreq(data,data$bb,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 124 72.51 72.51 65.18
## Si 47 27.49 100 20.95
## Total 171 100.00
## intervalo confianza inferior
## No 79.05
## Si 34.82
## Totalfreq(data,data$alfabloqueante,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 170 99.42 99.42 96.78
## Si 1 0.58 100 0.01
## Total 171 100.00
## intervalo confianza inferior
## No 99.99
## Si 3.22
## Totalfreq(data,data$tiazidas,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 157 91.81 91.81 86.64
## Si 14 8.19 100 4.55
## Total 171 100.00
## intervalo confianza inferior
## No 95.45
## Si 13.36
## Totalfreq(data,data$asa,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 152 88.89 88.89 83.19
## Si 19 11.11 100 6.82
## Total 171 100.00
## intervalo confianza inferior
## No 93.18
## Si 16.81
## Totalfreq(data,data$ahorrak,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 169 98.83 98.83 95.84
## Si 2 1.17 100 0.14
## Total 171 100.00
## intervalo confianza inferior
## No 99.86
## Si 4.16
## Totaldata$tratamientohta=data$ieca*1+data$ara2*2+data$bloqueoca*4+data$bb*8+data$alfabloqueante*16+data$tiazidas*32+data$asa*64+data$ahorrak*128
freq(data,data$tratamientohta,using_na=T,c("Sin antihipertensivos",
"Solo+IECA",
"Solo+ARA2",
"IECA+ARA2",
"CA",
"IECA+CA",
"ARA2+CA",
"IECA+ARA2+CA",
"BB",
"BB+IECA",
"BB+ARA2",
"BB+CA",
"IECA+BB+CA",
"ARA2+BB+CA",
"IECA+Tiazidas",
"IECA+CA+Tiazidas",
"ARA2+CA+Tiazidas",
"BB+Tiazidas",
"IECA+BB+Tiazidas",
"ARA2+BB+Tiazidas",
"IECA+asa",
"ARA2+asa",
"asa+ieca+ara2",
"asa+CA+ieca",
"asa+CA+ara2",
"asa+BB",
"asa+BB+ieca",
"asa+AB",
"asa+AB+BB+CA",
"AK+asa+ieca",
"AK+asa+BB+ieca"))
##
## absoluta relativa acumulada
## Sin antihipertensivos 28 16.37 16.37
## Solo+IECA 50 29.24 45.61
## Solo+ARA2 9 5.26 50.88
## IECA+ARA2 1 0.58 51.46
## CA 3 1.75 53.22
## IECA+CA 8 4.68 57.89
## ARA2+CA 1 0.58 58.48
## IECA+ARA2+CA 1 0.58 59.06
## BB 8 4.68 63.74
## BB+IECA 17 9.94 73.68
## BB+ARA2 2 1.17 74.85
## BB+CA 3 1.75 76.61
## IECA+BB+CA 6 3.51 80.12
## ARA2+BB+CA 2 1.17 81.29
## IECA+Tiazidas 5 2.92 84.21
## IECA+CA+Tiazidas 1 0.58 84.8
## ARA2+CA+Tiazidas 2 1.17 85.96
## BB+Tiazidas 1 0.58 86.55
## IECA+BB+Tiazidas 2 1.17 87.72
## ARA2+BB+Tiazidas 2 1.17 88.89
## IECA+asa 8 4.68 93.57
## ARA2+asa 2 1.17 94.74
## asa+ieca+ara2 1 0.58 95.32
## asa+CA+ieca 1 0.58 95.91
## asa+CA+ara2 1 0.58 96.49
## asa+BB 1 0.58 97.08
## asa+BB+ieca 1 0.58 97.66
## asa+AB 1 0.58 98.25
## asa+AB+BB+CA 1 0.58 98.83
## AK+asa+ieca 1 0.58 99.42
## AK+asa+BB+ieca 1 0.58 100
## Total 171 100.00
## intervalo confianza inferior
## Sin antihipertensivos 11.16
## Solo+IECA 22.55
## Solo+ARA2 2.43
## IECA+ARA2 0.01
## CA 0.36
## IECA+CA 2.04
## ARA2+CA 0.01
## IECA+ARA2+CA 0.01
## BB 2.04
## BB+IECA 5.9
## BB+ARA2 0.14
## BB+CA 0.36
## IECA+BB+CA 1.3
## ARA2+BB+CA 0.14
## IECA+Tiazidas 0.96
## IECA+CA+Tiazidas 0.01
## ARA2+CA+Tiazidas 0.14
## BB+Tiazidas 0.01
## IECA+BB+Tiazidas 0.14
## ARA2+BB+Tiazidas 0.14
## IECA+asa 2.04
## ARA2+asa 0.14
## asa+ieca+ara2 0.01
## asa+CA+ieca 0.01
## asa+CA+ara2 0.01
## asa+BB 0.01
## asa+BB+ieca 0.01
## asa+AB 0.01
## asa+AB+BB+CA 0.01
## AK+asa+ieca 0.01
## AK+asa+BB+ieca 0.01
## Total
## intervalo confianza inferior
## Sin antihipertensivos 22.79
## Solo+IECA 36.67
## Solo+ARA2 9.76
## IECA+ARA2 3.22
## CA 5.04
## IECA+CA 9.01
## ARA2+CA 3.22
## IECA+ARA2+CA 3.22
## BB 9.01
## BB+IECA 15.44
## BB+ARA2 4.16
## BB+CA 5.04
## IECA+BB+CA 7.48
## ARA2+BB+CA 4.16
## IECA+Tiazidas 6.69
## IECA+CA+Tiazidas 3.22
## ARA2+CA+Tiazidas 4.16
## BB+Tiazidas 3.22
## IECA+BB+Tiazidas 4.16
## ARA2+BB+Tiazidas 4.16
## IECA+asa 9.01
## ARA2+asa 4.16
## asa+ieca+ara2 3.22
## asa+CA+ieca 3.22
## asa+CA+ara2 3.22
## asa+BB 3.22
## asa+BB+ieca 3.22
## asa+AB 3.22
## asa+AB+BB+CA 3.22
## AK+asa+ieca 3.22
## AK+asa+BB+ieca 3.22
## Totalfreq(data,data$estatinas,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 36 21.05 21.05 15.2
## Si 135 78.95 100 72.07
## Total 171 100.00
## intervalo confianza inferior
## No 27.93
## Si 84.8
## Totalfreq(data,data$fibratos,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 159 92.98 92.98 88.06
## Si 12 7.02 100 3.68
## Total 171 100.00
## intervalo confianza inferior
## No 96.32
## Si 11.94
## Totaldata$tratamientodislipemia=data$estatinas*1+data$fibratos*2
freq(data,data$tratamientodislipemia,using_na=T,c("Sin tratamiento dislipemia","Estatinas","Fibratos","Estatinas+Fibratos"))
##
## absoluta relativa acumulada
## Sin tratamiento dislipemia 31 18.13 18.13
## Estatinas 128 74.85 92.98
## Fibratos 5 2.92 95.91
## Estatinas+Fibratos 7 4.09 100
## Total 171 100.00
## intervalo confianza inferior
## Sin tratamiento dislipemia 12.66
## Estatinas 67.66
## Fibratos 0.96
## Estatinas+Fibratos 1.66
## Total
## intervalo confianza inferior
## Sin tratamiento dislipemia 24.73
## Estatinas 81.16
## Fibratos 6.69
## Estatinas+Fibratos 8.25
## Totalfreq(data,data$ci,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 127 74.27 74.27 67.04
## Si 44 25.73 100 19.36
## Total 171 100.00
## intervalo confianza inferior
## No 80.64
## Si 32.96
## Totalfreq(data,data$acv,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 167 97.66 97.66 94.12
## Si 4 2.34 100 0.64
## Total 171 100.00
## intervalo confianza inferior
## No 99.36
## Si 5.88
## Totalfreq(data,data$aoc,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 161 94.15 94.15 89.51
## Si 10 5.85 100 2.84
## Total 171 100.00
## intervalo confianza inferior
## No 97.16
## Si 10.49
## Totalfreq(data,data$hta,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 33 19.3 19.3 13.67
## Si 138 80.7 100 73.98
## Total 171 100.00
## intervalo confianza inferior
## No 26.02
## Si 86.33
## Totalfreq(data,data$ic,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 138 80.7 80.7 73.98
## Si 33 19.3 100 13.67
## Total 171 100.00
## intervalo confianza inferior
## No 86.33
## Si 26.02
## Totalfreq(data,data$erc,using_na=T,c("Sin ERC","Estadio 1","Estadio 3","Con ERC etapa desconocida"))
##
## absoluta relativa acumulada
## Sin ERC 157 91.81 91.81
## Estadio 1 3 1.75 93.57
## Estadio 3 1 0.58 94.15
## Con ERC etapa desconocida 10 5.85 100
## Total 171 100.00
## intervalo confianza inferior
## Sin ERC 86.64
## Estadio 1 0.36
## Estadio 3 0.01
## Con ERC etapa desconocida 2.84
## Total
## intervalo confianza inferior
## Sin ERC 95.45
## Estadio 1 5.04
## Estadio 3 3.22
## Con ERC etapa desconocida 10.49
## Totalfreq(data,data$retinopatia,using_na=T,c("Sin RD","RDNP","RDP","Con RD tipo desconocido"))
##
## absoluta relativa acumulada
## Sin RD 143 83.63 83.63
## RDNP 9 5.26 88.89
## RDP 8 4.68 93.57
## Con RD tipo desconocido 11 6.43 100
## Total 171 100.00
## intervalo confianza inferior
## Sin RD 77.21
## RDNP 2.43
## RDP 2.04
## Con RD tipo desconocido 3.25
## Total
## intervalo confianza inferior
## Sin RD 88.84
## RDNP 9.76
## RDP 9.01
## Con RD tipo desconocido 11.22
## Totalfreq(data,data$neuropatiaperiferica,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 107 62.57 62.57 54.86
## Si 64 37.43 100 30.16
## Total 171 100.00
## intervalo confianza inferior
## No 69.84
## Si 45.14
## Totalfreq(data,data$neuropatiaautonomica,using_na=T,labeldicotomico)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 152 88.89 88.89 83.19
## Si 19 11.11 100 6.82
## Total 171 100.00
## intervalo confianza inferior
## No 93.18
## Si 16.81
## Totaldata$imc=data$peso/((data$talla/100)^2)
means(data,data$imc)##
## imc
## n 171
## Media 31.4
## Mediana 30.48
## DS 6.07
## Min 19.78
## Max 49.77
## Rango 30
## Q1 26.8
## Q3 35.37
## RIQ 8.57
## Normal 0
## ICI 30.49
## ICS 32.32
## Resumen Mediana[Q1-Q3] -->30.48(26.8-35.37)
means(data,data$imc,data$sexo)##
## imc by sexo = F imc by sexo = M
## n 120 51
## Media 31.65 30.82
## Mediana 30.81 29.92
## DS 5.95 6.36
## Min 19.83 19.78
## Max 49.15 49.77
## Rango 29.33 30
## Q1 26.86 26.62
## Q3 35.72 32.99
## RIQ 8.86 6.37
## Normal 0.023 0.002
## ICI 30.58 29.03
## ICS 32.73 32.61
## Resumen2 30.81 (26.86-35.72) 29.92 (26.62-32.99)
## Resumen1 Mediana[Q1-Q3]
## grupos Homocedastica
## Normalidad No
## Diferencias No
## p-value 0.317
## Metodo Wilcoxon
means(data,data$cintura)##
## cintura
## n 171
## Media 104.34
## Mediana 103
## DS 12.95
## Min 80
## Max 143
## Rango 63
## Q1 96
## Q3 112
## RIQ 16
## Normal 0
## ICI 102.39
## ICS 106.3
## Resumen Mediana[Q1-Q3] -->103(96-112)
means(data,data$cintura,data$sexo)##
## cintura by sexo = F cintura by sexo = M
## n 120 51
## Media 103.29 106.82
## Mediana 101 105
## DS 12.47 13.83
## Min 81 80
## Max 141 143
## Rango 60 63
## Q1 94 100.5
## Q3 112 112
## RIQ 18 11.5
## Normal 0.005 0.011
## ICI 101.03 102.93
## ICS 105.54 110.71
## Resumen2 101 (94-112) 105 (100.5-112)
## Resumen1 Mediana[Q1-Q3]
## grupos Homocedastica
## Normalidad No
## Diferencias No
## p-value 0.113
## Metodo Wilcoxon
means(data,data$glucemia)##
## glucemia
## n 171
## Media 148.36
## Mediana 133
## DS 61.15
## Min 1.67
## Max 355
## Rango 353.33
## Q1 109
## Q3 176.5
## RIQ 67.5
## Normal 0
## ICI 139.13
## ICS 157.59
## Resumen Mediana[Q1-Q3] -->133(109-176.5)
means(data,data$HbA1c)##
## HbA1c
## n 171
## Media 7.99
## Mediana 7.3
## DS 2.16
## Min 5.1
## Max 16.8
## Rango 11.7
## Q1 6.55
## Q3 9.1
## RIQ 2.55
## Normal 0
## ICI 7.66
## ICS 8.31
## Resumen Mediana[Q1-Q3] -->7.3(6.55-9.1)
means(data,data$creatininemia)##
## creatininemia
## n 171
## Media 0.87
## Mediana 0.8
## DS 0.31
## Min 0.35
## Max 2.87
## Rango 2.52
## Q1 0.7
## Q3 0.96
## RIQ 0.26
## Normal 0
## ICI 0.82
## ICS 0.92
## Resumen Mediana[Q1-Q3] -->0.8(0.7-0.96)
means(data,data$ct)##
## ct
## n 171
## Media 191.29
## Mediana 186
## DS 51.31
## Min 11
## Max 366
## Rango 355
## Q1 160.5
## Q3 213.5
## RIQ 53
## Normal 0
## ICI 183.55
## ICS 199.04
## Resumen Mediana[Q1-Q3] -->186(160.5-213.5)
means(data,data$hdl,data$sexo)##
## hdl by sexo = F hdl by sexo = M
## n 120 51
## Media 48.45 46.78
## Mediana 47 43
## DS 12.69 18.51
## Min 14 19
## Max 86 116
## Rango 72 97
## Q1 39 34
## Q3 56 53.5
## RIQ 17 19.5
## Normal 0.062 0
## ICI 46.16 41.58
## ICS 50.74 51.99
## Resumen2 47 (39-56) 43 (34-53.5)
## Resumen1 Mediana[Q1-Q3]
## grupos Heteroscedastica
## Normalidad No
## Diferencias No
## p-value 0.115
## Metodo Wilcoxon
means(data,data$tg)##
## tg
## n 171
## Media 166.07
## Mediana 140
## DS 95.08
## Min 10
## Max 665
## Rango 655
## Q1 103
## Q3 201.5
## RIQ 98.5
## Normal 0
## ICI 151.72
## ICS 180.42
## Resumen Mediana[Q1-Q3] -->140(103-201.5)
data$ldl=data$ct-data$tg/5-data$hdl
means(data,data$ldl)##
## ldl
## n 171
## Media 110.13
## Mediana 105.8
## DS 45.31
## Min -50.2
## Max 261.6
## Rango 311.8
## Q1 83.5
## Q3 130.4
## RIQ 46.9
## Normal 0
## ICI 103.29
## ICS 116.96
## Resumen Mediana[Q1-Q3] -->105.8(83.5-130.4)
freq(data,data$anio,using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## 2010 38 22.22 22.22 16.23
## 2011 38 22.22 44.44 16.23
## 2012 45 26.32 70.76 19.89
## 2013 32 18.71 89.47 13.17
## 2014 18 10.53 100 6.36
## Total 171 100.00
## intervalo confianza inferior
## 2010 29.2
## 2011 29.2
## 2012 33.58
## 2013 25.38
## 2014 16.13
## Total