encuesta.R
mauriciogutierrez
Sat Jun 3 11:28:07 2017
path="/Users/mauriciogutierrez/Library/Mobile Documents/com~apple~CloudDocs/estadisticas";
pathfunctions=paste(path,"/general_functions_new.R",sep="")
pathdatabase=paste(path,"/mercepin/datos.xlsx",sep="")
source(pathfunctions)## 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(pathdatabase,"Hoja1")
freq(data,data$sexo,c("Femenino","Masculino"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## Femenino 62 79.49 79.49 68.84
## Masculino 16 20.51 100 12.2
## Total 78 100.00
## intervalo confianza inferior
## Femenino 87.8
## Masculino 31.16
## Totalmeans(data,data$edad)##
## edad
## n 78
## Media 54.92
## Mediana 56
## DS 18.09
## Min 19
## Max 89
## Rango 70
## Q1 44
## Q3 68
## RIQ 24
## Normal 0.011
## ICI 50.84
## ICS 59
## Resumen Mediana[Q1-Q3] -->56(44-68)
freq(data,data$lugarresidencia,c("Montevideo","Interior"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## Montevideo 44 56.41 56.41 44.7
## Interior 34 43.59 100 32.39
## Total 78 100.00
## intervalo confianza inferior
## Montevideo 67.61
## Interior 55.3
## Totalfreq(data,data$urbanorural,c("Urbano","Rural"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## Urbano 66 84.62 84.62 74.67
## Rural 12 15.38 100 8.21
## Total 78 100.00
## intervalo confianza inferior
## Urbano 91.79
## Rural 25.33
## Totalfreq(data,data$estadocivil,c("Soltero","Casado","Divorciado","Viudo"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## Soltero 28 35.9 35.9 25.34
## Casado 30 38.46 74.36 27.66
## Divorciado 9 11.54 85.9 5.41
## Viudo 11 14.1 100 7.26
## Total 78 100.00
## intervalo confianza inferior
## Soltero 47.56
## Casado 50.17
## Divorciado 20.78
## Viudo 23.83
## Totalfreq(data,data$situacionlaboral,c("Empleado","Desempleado","Jubilado","Pensionista","Ama de casa","Estudiante"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## Empleado 17 21.79 21.79 13.24
## Desempleado 8 10.26 32.05 4.53
## Jubilado 22 28.21 60.26 18.59
## Pensionista 12 15.38 75.64 8.21
## Ama de casa 15 19.23 94.87 11.18
## Estudiante 4 5.13 100 1.41
## Total 78 100.00
## intervalo confianza inferior
## Empleado 32.59
## Desempleado 19.21
## Jubilado 39.53
## Pensionista 25.33
## Ama de casa 29.73
## Estudiante 12.61
## Totalfreq(data,data$escolaridad,c("Primaria incompleta","primaria completa","secundaria incompleta","secundaria completa","terciaria incompleta","terciaria completa"),using_na=T)
##
## absoluta relativa acumulada
## Primaria incompleta 14 17.95 17.95
## primaria completa 25 32.05 50
## secundaria incompleta 27 34.62 84.62
## secundaria completa 2 2.56 87.18
## terciaria incompleta 4 5.13 92.31
## terciaria completa 6 7.69 100
## Total 78 100.00
## intervalo confianza inferior
## Primaria incompleta 10.17
## primaria completa 21.93
## secundaria incompleta 24.2
## secundaria completa 0.31
## terciaria incompleta 1.41
## terciaria completa 2.88
## Total
## intervalo confianza inferior
## Primaria incompleta 28.28
## primaria completa 43.58
## secundaria incompleta 46.24
## secundaria completa 8.96
## terciaria incompleta 12.61
## terciaria completa 15.99
## Totalfreq(data,data$fuma,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 61 79.22 79.22 68.46
## Si 16 20.78 100 12.37
## Total 77 100.00
## intervalo confianza inferior
## No 87.63
## Si 31.54
## Totalfreq(data,data$consumodiariocigarros,using_na = T)
##
## absoluta relativa acumulada intervalo confianza inferior
## 2 1 5.88 5.88 0.15
## 3 2 11.76 17.65 1.46
## 5 2 11.76 29.41 1.46
## 10 5 29.41 58.82 10.31
## 20 6 35.29 94.12 14.21
## 40 1 5.88 100 0.15
## Total 17 100.00
## intervalo confianza inferior
## 2 28.69
## 3 36.44
## 5 36.44
## 10 55.96
## 20 61.67
## 40 28.69
## Totalfuman=data[which(data$fuma==2),]
freq(fuman,fuman$consumodiario,using_na = T)
##
## absoluta relativa acumulada intervalo confianza inferior
## 2 1 6.67 6.67 0.17
## 3 2 13.33 20 1.66
## 5 2 13.33 33.33 1.66
## 10 4 26.67 60 7.79
## 20 6 40 100 16.34
## Total 15 100.00
## intervalo confianza inferior
## 2 31.95
## 3 40.46
## 5 40.46
## 10 55.1
## 20 67.71
## Totalmeans(fuman,fuman$edadcomienzo)##
## edadcomienzo
## n 13
## Media 16.08
## Mediana 15
## DS 3.17
## Min 13
## Max 25
## Rango 12
## Q1 15
## Q3 16
## RIQ 1
## Normal 0.001
## ICI 14.16
## ICS 18
## Resumen Mediana[Q1-Q3] -->15(15-16)
freq(data,data$fumoalgunavez,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 40 65.57 65.57 52.31
## Si 21 34.43 100 22.73
## Total 61 100.00
## intervalo confianza inferior
## No 77.27
## Si 47.69
## Totalmeans(fuman,fuman$cuantosaniosfumo)##
## cuantosaniosfumo
## n 3
## Media 40.33
## Mediana 45
## DS 10.79
## Min 28
## Max 48
## Rango 20
## Q1 36.5
## Q3 46.5
## RIQ 10
## Normal 0.266
## ICI 13.54
## ICS 67.13
## Resumen Media[ICI-ICS] --> 40.33(13.54-67.13)
means(fuman,fuman$hacecuantodejo)##
## hacecuantodejo
## n 1
## Media 7
## Mediana 7
## DS <NA>
## Min 7
## Max 7
## Rango 0
## Q1 7
## Q3 7
## RIQ 0
## Normal -1
## ICI -1
## ICS -1
## Resumen Mediana[Q1-Q3] -->7(7-7)
freq(data,data$consumealcohol,using_na = T)
##
## absoluta relativa acumulada intervalo confianza inferior
## 1 63 81.82 81.82 71.38
## 2 14 18.18 100 10.31
## Total 77 100.00
## intervalo confianza inferior
## 1 89.69
## 2 28.62
## Totaltomanalcohol=data[which(data$consumealcohol==2),]
means(tomanalohol,tomanalcohol$frecuenciasemanal)##
## frecuenciasemanal
## n 14
## Media 2.5
## Mediana 1
## DS 2.47
## Min 1
## Max 7
## Rango 6
## Q1 1
## Q3 2
## RIQ 1
## Normal 0
## ICI 1.07
## ICS 3.93
## Resumen Mediana[Q1-Q3] -->1(1-2)
#means(data,data$diassemanacomefrutas)
#consumefrutas=data[which(data$diassemanacomefrutas==2),]
means(data,data$diassemanacomefrutas)##
## diassemanacomefrutas
## n 77
## Media 4.13
## Mediana 5
## DS 2.61
## Min 0
## Max 7
## Rango 7
## Q1 2
## Q3 7
## RIQ 5
## Normal 0
## ICI 3.54
## ICS 4.72
## Resumen Mediana[Q1-Q3] -->5(2-7)
means(data,data$diasemanacomeverduras)##
## diasemanacomeverduras
## n 77
## Media 3.99
## Mediana 4
## DS 2.58
## Min 0
## Max 7
## Rango 7
## Q1 2
## Q3 7
## RIQ 5
## Normal 0
## ICI 3.4
## ICS 4.57
## Resumen Mediana[Q1-Q3] -->4(2-7)
consumeverduras=data[which(data$diasemanacomeverduras==2),]
means(consumeverduras,consumeverduras$verdurasdia)##
## verdurasdia
## n 12
## Media 1.42
## Mediana 1
## DS 0.67
## Min 1
## Max 3
## Rango 2
## Q1 1
## Q3 2
## RIQ 1
## Normal 0
## ICI 0.99
## ICS 1.84
## Resumen Mediana[Q1-Q3] -->1(1-2)
freq(data,data$cocinacon,c("aceite de soja","aceite de girasol","aceite de oliva","Aceite de arroz","Manteca"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## aceite de soja 25 32.47 32.47 22.23
## aceite de girasol 38 49.35 81.82 37.76
## aceite de oliva 10 12.99 94.81 6.41
## Aceite de arroz 2 2.6 97.4 0.32
## Manteca 2 2.6 100 0.32
## Total 77 100.00
## intervalo confianza inferior
## aceite de soja 44.1
## aceite de girasol 61
## aceite de oliva 22.59
## Aceite de arroz 9.07
## Manteca 9.07
## Totalfreq(data,data$comidalistacomprados,using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## 0 48 64 64 52.09
## 1 23 30.67 94.67 20.53
## 2 2 2.67 97.33 0.32
## 6 1 1.33 98.67 0.03
## 7 1 1.33 100 0.03
## Total 75 100.00
## intervalo confianza inferior
## 0 74.77
## 1 42.38
## 2 9.3
## 6 7.21
## 7 7.21
## Totalfreq(data,data$trabajovigorosomayor10min,c("No"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 77 100 100 95.32
## Total 77 100.00
## intervalo confianza inferior
## No 100
## Totalfreq(data,data$trabajomoderadomayor10min,c("No","Si"),using_na = T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 67 87.01 87.01 77.41
## Si 10 12.99 100 6.41
## Total 77 100.00
## intervalo confianza inferior
## No 93.59
## Si 22.59
## Totaltrabajomoderadomayor10min=data[which(data$trabajomoderadomayor10min>=0),]
means(data,trabajomoderadomayor10min$minutos.dias.1)##
## minutos.dias.1
## n 77
## Media 19.35
## Mediana 0
## DS 55.23
## Min 0
## Max 180
## Rango 180
## Q1 0
## Q3 0
## RIQ 0
## Normal 0
## ICI 6.82
## ICS 31.89
## Resumen Mediana[Q1-Q3] -->0(0-0)
means(data,trabajomoderadomayor10min$dias.semana.1)##
## dias.semana.1
## n 77
## Media 0.56
## Mediana 0
## DS 1.58
## Min 0
## Max 7
## Rango 7
## Q1 0
## Q3 0
## RIQ 0
## Normal 0
## ICI 0.2
## ICS 0.92
## Resumen Mediana[Q1-Q3] -->0(0-0)
freq(data,data$desplazamientomayordiezminutobicicleta,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 50 64.94 64.94 53.22
## Si 27 35.06 100 24.53
## Total 77 100.00
## intervalo confianza inferior
## No 75.47
## Si 46.78
## Totaldesplazamientomayordiezminutobicicleta=data[which(data$desplazamientomayordiezminutobicicleta>=0),]
means(data,desplazamientomayordiezminutobicicleta$dias.semana.2)##
## dias.semana.2
## n 77
## Media 1.22
## Mediana 0
## DS 1.98
## Min 0
## Max 7
## Rango 7
## Q1 0
## Q3 3
## RIQ 3
## Normal 0
## ICI 0.77
## ICS 1.67
## Resumen Mediana[Q1-Q3] -->0(0-3)
means(data,desplazamientomayordiezminutobicicleta$minutos.semana.2)##
## minutos.semana.2
## n 30
## Media 76.33
## Mediana 60
## DS 69.54
## Min 0
## Max 300
## Rango 300
## Q1 30
## Q3 100
## RIQ 70
## Normal 0.001
## ICI 50.37
## ICS 102.3
## Resumen Mediana[Q1-Q3] -->60(30-100)
freq(data,data$deportevigoroso,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 74 96.1 96.1 89.03
## Si 3 3.9 100 0.81
## Total 77 100.00
## intervalo confianza inferior
## No 99.19
## Si 10.97
## Totaldeportevigoroso=data[which(data$deportevigoroso>=0),]
means(data,deportevigoroso$dias.semana.3)##
## dias.semana.3
## n 77
## Media 0.06
## Mediana 0
## DS 0.37
## Min 0
## Max 3
## Rango 3
## Q1 0
## Q3 0
## RIQ 0
## Normal 0
## ICI -0.02
## ICS 0.15
## Resumen Mediana[Q1-Q3] -->0(0-0)
means(data,deportevigoroso$min.semana.3)##
## min.semana.3
## n 77
## Media 3.51
## Mediana 0
## DS 21.75
## Min 0
## Max 180
## Rango 180
## Q1 0
## Q3 0
## RIQ 0
## Normal 0
## ICI -1.43
## ICS 8.44
## Resumen Mediana[Q1-Q3] -->0(0-0)
freq(data,data$deportemoderado,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 56 72.73 72.73 61.38
## Si 21 27.27 100 17.74
## Total 77 100.00
## intervalo confianza inferior
## No 82.26
## Si 38.62
## Totaldeportemoderado=data[which(data$deportemoderado>=2),]
means(data,deportemoderado$dias.semana.4)##
## dias.semana.4
## n 21
## Media 3.48
## Mediana 3
## DS 2.04
## Min 1
## Max 7
## Rango 6
## Q1 2
## Q3 5
## RIQ 3
## Normal 0.035
## ICI 2.55
## ICS 4.4
## Resumen Mediana[Q1-Q3] -->3(2-5)
means(data,deportemoderado$minutos.semana.4)##
## minutos.semana.4
## n 21
## Media 174.76
## Mediana 150
## DS 120.99
## Min 30
## Max 450
## Rango 420
## Q1 90
## Q3 210
## RIQ 120
## Normal 0.008
## ICI 119.69
## ICS 229.84
## Resumen Mediana[Q1-Q3] -->150(90-210)
data$totalactividadfisica=data$minutos.dias.1*data$dias.semana.1+data$minutos.semana.2+data$minutos.semana.4
means(data,data$totalactividadfisica)##
## totalactividadfisica
## n 30
## Media 268.83
## Mediana 150
## DS 292.95
## Min 20
## Max 1080
## Rango 1060
## Q1 60
## Q3 350
## RIQ 290
## Normal 0
## ICI 159.44
## ICS 378.22
## Resumen Mediana[Q1-Q3] -->150(60-350)
means(data,data$sedentario)##
## sedentario
## n 75
## Media 6.97
## Mediana 8
## DS 2.65
## Min 1
## Max 12
## Rango 11
## Q1 5
## Q3 8
## RIQ 3
## Normal 0.011
## ICI 6.36
## ICS 7.58
## Resumen Mediana[Q1-Q3] -->8(5-8)
freq(data,data$hta,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 42 55.26 55.26 43.41
## Si 34 44.74 100 33.31
## Total 76 100.00
## intervalo confianza inferior
## No 66.69
## Si 56.59
## Totalfreq(data,data$dieta.a,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 55 72.37 72.37 60.91
## Si 21 27.63 100 17.99
## Total 76 100.00
## intervalo confianza inferior
## No 82.01
## Si 39.09
## Totalfreq(data,data$farmacos.a,c("Sin farmacos","IECA/ARA II","beta-bloqueantes+IECA/ARA II","Beta-bloqueantes+IECA/ARA II+HCT","calcio antagonista","Calcioantagonista+IECA/ARA II","calcioantagonista+ieca/ARA II+HCT"),using_na=T)
##
## absoluta relativa acumulada
## Sin farmacos 50 65.79 65.79
## IECA/ARA II 20 26.32 92.11
## beta-bloqueantes+IECA/ARA II 2 2.63 94.74
## Beta-bloqueantes+IECA/ARA II+HCT 1 1.32 96.05
## calcio antagonista 1 1.32 97.37
## Calcioantagonista+IECA/ARA II 1 1.32 98.68
## calcioantagonista+ieca/ARA II+HCT 1 1.32 100
## Total 76 100.00
## intervalo confianza inferior
## Sin farmacos 54.01
## IECA/ARA II 16.87
## beta-bloqueantes+IECA/ARA II 0.32
## Beta-bloqueantes+IECA/ARA II+HCT 0.03
## calcio antagonista 0.03
## Calcioantagonista+IECA/ARA II 0.03
## calcioantagonista+ieca/ARA II+HCT 0.03
## Total
## intervalo confianza inferior
## Sin farmacos 76.29
## IECA/ARA II 37.68
## beta-bloqueantes+IECA/ARA II 9.18
## Beta-bloqueantes+IECA/ARA II+HCT 7.11
## calcio antagonista 7.11
## Calcioantagonista+IECA/ARA II 7.11
## calcioantagonista+ieca/ARA II+HCT 7.11
## Totalfreq(data,data$dm,c("No","Si"),using_na = T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 48 63.16 63.16 51.31
## Si 28 36.84 100 26.06
## Total 76 100.00
## intervalo confianza inferior
## No 73.94
## Si 48.69
## Totalfreq(data,data$dieta.b,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 63 82.89 82.89 72.53
## Si 13 17.11 100 9.43
## Total 76 100.00
## intervalo confianza inferior
## No 90.57
## Si 27.47
## Totalfreq(data,data$farmacos.b,c("Sin farmacos","insulinosensibilizadores","insulinosensibilizadores+insulinosecretores","Incretinas","Insulinas","Insulinas+insulinosensibilizadores"),using_na=T)
##
## absoluta relativa acumulada
## Sin farmacos 54 71.05 71.05
## insulinosensibilizadores 12 15.79 86.84
## insulinosensibilizadores+insulinosecretores 2 2.63 89.47
## Incretinas 1 1.32 90.79
## Insulinas 5 6.58 97.37
## Insulinas+insulinosensibilizadores 2 2.63 100
## Total 76 100.00
## intervalo confianza inferior
## Sin farmacos 59.51
## insulinosensibilizadores 8.43
## insulinosensibilizadores+insulinosecretores 0.32
## Incretinas 0.03
## Insulinas 2.17
## Insulinas+insulinosensibilizadores 0.32
## Total
## intervalo confianza inferior
## Sin farmacos 80.89
## insulinosensibilizadores 25.96
## insulinosensibilizadores+insulinosecretores 9.18
## Incretinas 7.11
## Insulinas 14.69
## Insulinas+insulinosensibilizadores 9.18
## Totalfreq(data,data$dislipemias,c("No","Si"),using_na = T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 44 57.89 57.89 46.02
## Si 32 42.11 100 30.86
## Total 76 100.00
## intervalo confianza inferior
## No 69.14
## Si 53.98
## Totalfreq(data,data$Dieta.c,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 65 85.53 85.53 75.58
## Si 11 14.47 100 7.45
## Total 76 100.00
## intervalo confianza inferior
## No 92.55
## Si 24.42
## Totalfreq(data,data$Farmacos.c,c("Sin farmacos","Estatinas","Omega 3"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## Sin farmacos 53 69.74 69.74 58.13
## Estatinas 22 28.95 98.68 19.11
## Omega 3 1 1.32 100 0.03
## Total 76 100.00
## intervalo confianza inferior
## Sin farmacos 79.75
## Estatinas 40.49
## Omega 3 7.11
## Totalfreq(data,data$cardiopatiaisquemica,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 69 92 92 83.4
## Si 6 8 100 2.99
## Total 75 100.00
## intervalo confianza inferior
## No 97.01
## Si 16.6
## Totalfreq(data,data$acv,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 72 96 96 88.75
## Si 3 4 100 0.83
## Total 75 100.00
## intervalo confianza inferior
## No 99.17
## Si 11.25
## Totalmeans(data,data$estatura)##
## estatura
## n 68
## Media 1.6
## Mediana 1.61
## DS 0.11
## Min 1.39
## Max 1.89
## Rango 0.5
## Q1 1.53
## Q3 1.65
## RIQ 0.12
## Normal 0.099
## ICI 1.58
## ICS 1.63
## Resumen Media[ICI-ICS] --> 1.6(1.58-1.63)
means(data,data$peso)##
## peso
## n 69
## Media 76.78
## Mediana 73
## DS 21
## Min 47
## Max 165
## Rango 118
## Q1 64
## Q3 86
## RIQ 22
## Normal 0
## ICI 71.74
## ICS 81.83
## Resumen Mediana[Q1-Q3] -->73(64-86)
means(data,data$imc)##
## imc
## n 68
## Media 30.06
## Mediana 28.61
## DS 8.77
## Min 19.72
## Max 64.97
## Rango 45.24
## Q1 24.36
## Q3 31.8
## RIQ 7.44
## Normal 0
## ICI 27.94
## ICS 32.18
## Resumen Mediana[Q1-Q3] -->28.61(24.36-31.8)
means(data,data$cinturaabdominal)##
## cinturaabdominal
## n 70
## Media 101.47
## Mediana 97.5
## DS 19.16
## Min 69
## Max 170
## Rango 101
## Q1 90.25
## Q3 111
## RIQ 20.75
## Normal 0
## ICI 96.9
## ICS 106.04
## Resumen Mediana[Q1-Q3] -->97.5(90.25-111)
means(data,data$cinturapelvica)##
## cinturapelvica
## n 68
## Media 106.76
## Mediana 105
## DS 12.84
## Min 87
## Max 156
## Rango 69
## Q1 98
## Q3 112
## RIQ 14
## Normal 0
## ICI 103.66
## ICS 109.87
## Resumen Mediana[Q1-Q3] -->105(98-112)
means(data,data$icc)##
## icc
## n 68
## Media 0.95
## Mediana 0.95
## DS 0.12
## Min 0.72
## Max 1.34
## Rango 0.62
## Q1 0.87
## Q3 0.99
## RIQ 0.12
## Normal 0.002
## ICI 0.92
## ICS 0.98
## Resumen Mediana[Q1-Q3] -->0.95(0.87-0.99)
means(data,data$frecuenciacardiaca)##
## frecuenciacardiaca
## n 45
## Media 77.07
## Mediana 80
## DS 7.81
## Min 60
## Max 90
## Rango 30
## Q1 70
## Q3 80
## RIQ 10
## Normal 0.019
## ICI 74.72
## ICS 79.41
## Resumen Mediana[Q1-Q3] -->80(70-80)
means(data,as.numeric(data$PAS))##
## PAS)
## n 63
## Media 119.52
## Mediana 120
## DS 23.17
## Min 90
## Max 200
## Rango 110
## Q1 100
## Q3 130
## RIQ 30
## Normal 0
## ICI 113.69
## ICS 125.36
## Resumen Mediana[Q1-Q3] -->120(100-130)
means(data,as.numeric(data$PAD))##
## PAD)
## n 63
## Media 72.38
## Mediana 70
## DS 10.88
## Min 50
## Max 100
## Rango 50
## Q1 60
## Q3 80
## RIQ 20
## Normal 0
## ICI 69.64
## ICS 75.12
## Resumen Mediana[Q1-Q3] -->70(60-80)
means(data,data$glicemia)##
## glicemia
## n 53
## Media 129.32
## Mediana 103
## DS 69.66
## Min 73
## Max 422
## Rango 349
## Q1 90
## Q3 135
## RIQ 45
## Normal 0
## ICI 110.12
## ICS 148.52
## Resumen Mediana[Q1-Q3] -->103(90-135)
means(data,data$hba1c)##
## hba1c
## n 22
## Media 7.69
## Mediana 7.05
## DS 2.22
## Min 5.3
## Max 13.1
## Rango 7.8
## Q1 5.82
## Q3 8.93
## RIQ 3.1
## Normal 0.013
## ICI 6.71
## ICS 8.67
## Resumen Mediana[Q1-Q3] -->7.05(5.82-8.93)
means(data,data$coltotal)##
## coltotal
## n 53
## Media 199.43
## Mediana 200
## DS 43.75
## Min 104
## Max 306
## Rango 202
## Q1 177
## Q3 221
## RIQ 44
## Normal 0.599
## ICI 187.37
## ICS 211.49
## Resumen Media[ICI-ICS] --> 199.43(187.37-211.49)
means(data,data$hdl)##
## hdl
## n 53
## Media 54.75
## Mediana 51
## DS 18.11
## Min 29
## Max 107
## Rango 78
## Q1 43
## Q3 62
## RIQ 19
## Normal 0
## ICI 49.76
## ICS 59.75
## Resumen Mediana[Q1-Q3] -->51(43-62)
means(data,data$ldl)##
## ldl
## n 52
## Media 114.54
## Mediana 117
## DS 39.62
## Min 37
## Max 226
## Rango 189
## Q1 88.25
## Q3 133.25
## RIQ 45
## Normal 0.297
## ICI 103.51
## ICS 125.57
## Resumen Media[ICI-ICS] --> 114.54(103.51-125.57)
means(data,data$colnohdl)##
## colnohdl
## n 45
## Media 116.4
## Mediana 125
## DS 66.64
## Min 0
## Max 248
## Rango 248
## Q1 86
## Q3 152
## RIQ 66
## Normal 0.005
## ICI 96.38
## ICS 136.42
## Resumen Mediana[Q1-Q3] -->125(86-152)
means(data,data$tg)##
## tg
## n 53
## Media 150
## Mediana 123
## DS 73.11
## Min 38
## Max 361
## Rango 323
## Q1 99
## Q3 181
## RIQ 82
## Normal 0.001
## ICI 129.85
## ICS 170.15
## Resumen Mediana[Q1-Q3] -->123(99-181)
diabeticos=data[which(data$dm==2),]
lineal(diabeticos$totalactividadfisica,diabeticos$hba1c,labely="Actividad Fisica",labelx="HbA1c")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 701.7 -43.3
##
## [1] "r2=0.14"
## [1] "p=0.535"
lineal(data$totalactividadfisica,data$coltotal,labely="Actividad Fisica",labelx="CT")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## -78.376 1.774
##
## [1] "r2=0.066"
## [1] "p=0.225"
lineal(data$totalactividadfisica,data$hdl,labely="Actividad Fisica",labelx="hdl")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 194.406 1.627
##
## [1] "r2=0.008"
## [1] "p=0.68"
lineal(data$totalactividadfisica,data$ldl,labely="Actividad Fisica",labelx="ldl")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 4.107 2.308
##
## [1] "r2=0.09"
## [1] "p=0.153"
lineal(data$totalactividadfisica,data$colnohdl,labely="Actividad Fisica",labelx="colnohdl")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## -9.020 1.686
##
## [1] "r2=0.133"
## [1] "p=0.124"
lineal(data$totalactividadfisica,data$tg,labely="Actividad Fisica",labelx="tg")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 493.155 -1.582
##
## [1] "r2=0.059"
## [1] "p=0.254"
means(data,data$totalactividadfisica,data$lugarresidencia,c("Montevideo","Interior"))##
## totalactividadfisica by lugarresidencia = Montevideo
## n 15
## Media 331
## Mediana 160
## DS 359.01
## Min 20
## Max 1080
## Rango 1060
## Q1 70
## Q3 435
## RIQ 365
## Normal 0.003
## ICI 132.19
## ICS 529.81
## Resumen2 160 (70-435)
## Resumen1 Mediana[Q1-Q3]
## grupos Heteroscedastica
## Normalidad No
## Diferencias No
## p-value 0.506
## Metodo Kruskal-Wallis
## totalactividadfisica by lugarresidencia = Interior
## n 15
## Media 206.67
## Mediana 100
## DS 201.49
## Min 20
## Max 705
## Rango 685
## Q1 60
## Q3 285
## RIQ 225
## Normal 0.01
## ICI 95.08
## ICS 318.25
## Resumen2 100 (60-285)
## Resumen1
## grupos
## Normalidad
## Diferencias
## p-value
## Metodo
means(data,data$totalactividadfisica,data$educacion,c("Primaria incompleta","primaria completa","secundaria incompleta","secundaria completa","terciaria incompleta","terciaria completa"))##
## totalactividadfisica
## n 30
## Media 268.83
## Mediana 150
## DS 292.95
## Min 20
## Max 1080
## Rango 1060
## Q1 60
## Q3 350
## RIQ 290
## Normal 0
## ICI 159.44
## ICS 378.22
## Resumen Mediana[Q1-Q3] -->150(60-350)
means(data,data$totalactividadfisica,data$escolaridad)##
## totalactividadfisica by escolaridad = 1
## n 5
## Media 385
## Mediana 320
## DS 362.97
## Min 60
## Max 930
## Rango 870
## Q1 75
## Q3 540
## RIQ 465
## Normal 0.436
## ICI -65.69
## ICS 835.69
## Resumen2 320 (75-540)
## Resumen1 Mediana[Q1-Q3]
## grupos Heteroscedastica
## Normalidad No
## Diferencias No
## p-value 0.337
## Metodo Kruskal-Wallis
## totalactividadfisica by escolaridad = 2
## n 6
## Media 148.33
## Mediana 135
## DS 134.67
## Min 20
## Max 360
## Rango 340
## Q1 37.5
## Q3 210
## RIQ 172.5
## Normal 0.266
## ICI 7
## ICS 289.66
## Resumen2 135 (37.5-210)
## Resumen1
## grupos
## Normalidad
## Diferencias
## p-value
## Metodo
## totalactividadfisica by escolaridad = 3
## n 13
## Media 220.77
## Mediana 100
## DS 273.47
## Min 20
## Max 900
## Rango 880
## Q1 60
## Q3 240
## RIQ 180
## Normal 0.001
## ICI 55.51
## ICS 386.03
## Resumen2 100 (60-240)
## Resumen1
## grupos
## Normalidad
## Diferencias
## p-value
## Metodo
## totalactividadfisica by escolaridad = 4
## n 1
## Media 360
## Mediana 360
## DS <NA>
## Min 360
## Max 360
## Rango 0
## Q1 360
## Q3 360
## RIQ 0
## Normal -1
## ICI -1
## ICS -1
## Resumen2 360 (360-360)
## Resumen1
## grupos
## Normalidad
## Diferencias
## p-value
## Metodo
## totalactividadfisica by escolaridad = 5
## n 1
## Media 80
## Mediana 80
## DS <NA>
## Min 80
## Max 80
## Rango 0
## Q1 80
## Q3 80
## RIQ 0
## Normal -1
## ICI -1
## ICS -1
## Resumen2 80 (80-80)
## Resumen1
## grupos
## Normalidad
## Diferencias
## p-value
## Metodo
## totalactividadfisica by escolaridad = 6
## n 4
## Media 485
## Mediana 380
## DS 431.32
## Min 100
## Max 1080
## Rango 980
## Q1 212.5
## Q3 652.5
## RIQ 440
## Normal 0.54
## ICI -201.32
## ICS 1171.32
## Resumen2 380 (212.5-652.5)
## Resumen1
## grupos
## Normalidad
## Diferencias
## p-value
## Metodo
lineal(data$totalactividadfisica,data$frutasdia,labely="Actividad Fisica",labelx="Frutas al dia")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 249.98 10.88
##
## [1] "r2=0.002"
## [1] "p=0.802"
lineal(data$totalactividadfisica,data$verdurasdia,labely="Actividad Fisica",labelx="Verduras al dia")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 372.88 -74.32
##
## [1] "r2=0.029"
## [1] "p=0.366"
means(data,data$totalactividadfisica,data$fuma,c("No","Si"))##
## totalactividadfisica by fuma = No
## n 24
## Media 310
## Mediana 210
## DS 312.86
## Min 20
## Max 1080
## Rango 1060
## Q1 71.25
## Q3 397.5
## RIQ 326.25
## Normal 0.001
## ICI 177.89
## ICS 442.11
## Resumen2 210 (71.25-397.5)
## Resumen1 Mediana[Q1-Q3]
## grupos Heteroscedastica
## Normalidad No
## Diferencias No
## p-value 0.107
## Metodo Kruskal-Wallis
## totalactividadfisica by fuma = Si
## n 6
## Media 104.17
## Mediana 80
## DS 82.61
## Min 20
## Max 240
## Rango 220
## Q1 48.75
## Q3 145
## RIQ 96.25
## Normal 0.512
## ICI 17.47
## ICS 190.86
## Resumen2 80 (48.75-145)
## Resumen1
## grupos
## Normalidad
## Diferencias
## p-value
## Metodo
# Nuevos solicitudes al 2/6
data$frutasyverdurasdia=data$frutasdia+data$verdurasdia
means(data,data$frutasyverdurasdia)##
## frutasyverdurasdia
## n 77
## Media 2.94
## Mediana 3
## DS 1.63
## Min 0
## Max 8
## Rango 8
## Q1 2
## Q3 4
## RIQ 2
## Normal 0
## ICI 2.56
## ICS 3.31
## Resumen Mediana[Q1-Q3] -->3(2-4)
data$frutasyverdurastotales=data$frutasdia*data$diassemanacomefrutas+data$verdurasdia*data$diasemanacomeverduras
means(data,data$frutasyverdurastotales)##
## frutasyverdurastotales
## n 77
## Media 14.26
## Mediana 10
## DS 12.08
## Min 0
## Max 56
## Rango 56
## Q1 5
## Q3 21
## RIQ 16
## Normal 0
## ICI 11.52
## ICS 17
## Resumen Mediana[Q1-Q3] -->10(5-21)
lineal(data$totalactividadfisica,data$frutasyverdurastotales,labely="Actividad Fisica",labelx="Frutas y Verduras totales")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 275.2529 -0.3867
##
## [1] "r2=0"
## [1] "p=0.923"
lineal(data$totalactividadfisica,data$frutasyverdurasdia,labely="Actividad Fisica",labelx="Frutas y Verduras dia")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 286.538 -5.651
##
## [1] "r2=0.001"
## [1] "p=0.865"
hist(data$totalactividadfisica,breaks = c(0,150,300,450,600,750,900,1050,1200),main=("Total Actividad Fisica"),xlab="actividad fisica")
hipertensos=data[which(data$hta==2),]
freq(hipertensos,hipertensos$farmacos.a,c("Sin farmacos","IECA/ARA II","beta-bloqueantes+IECA/ARA II","Beta-bloqueantes+IECA/ARA II+HCT","calcio antagonista","Calcioantagonista+IECA/ARA II","calcioantagonista+ieca/ARA II+HCT"),using_na=T)
##
## absoluta relativa acumulada
## Sin farmacos 8 23.53 23.53
## IECA/ARA II 20 58.82 82.35
## beta-bloqueantes+IECA/ARA II 2 5.88 88.24
## Beta-bloqueantes+IECA/ARA II+HCT 1 2.94 91.18
## calcio antagonista 1 2.94 94.12
## Calcioantagonista+IECA/ARA II 1 2.94 97.06
## calcioantagonista+ieca/ARA II+HCT 1 2.94 100
## Total 34 100.00
## intervalo confianza inferior
## Sin farmacos 10.75
## IECA/ARA II 40.7
## beta-bloqueantes+IECA/ARA II 0.72
## Beta-bloqueantes+IECA/ARA II+HCT 0.07
## calcio antagonista 0.07
## Calcioantagonista+IECA/ARA II 0.07
## calcioantagonista+ieca/ARA II+HCT 0.07
## Total
## intervalo confianza inferior
## Sin farmacos 41.17
## IECA/ARA II 75.35
## beta-bloqueantes+IECA/ARA II 19.68
## Beta-bloqueantes+IECA/ARA II+HCT 15.33
## calcio antagonista 15.33
## Calcioantagonista+IECA/ARA II 15.33
## calcioantagonista+ieca/ARA II+HCT 15.33
## Totalfreq(hipertensos,hipertensos$dieta.a,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 13 38.24 38.24 22.17
## Si 21 61.76 100 43.56
## Total 34 100.00
## intervalo confianza inferior
## No 56.44
## Si 77.83
## Totaldiabetico=data[which(data$dm==2),]
freq(diabetico,diabetico$dieta.b,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 15 53.57 53.57 33.87
## Si 13 46.43 100 27.51
## Total 28 100.00
## intervalo confianza inferior
## No 72.49
## Si 66.13
## Totalfreq(diabetico,diabetico$farmacos.b,c("Sin farmacos","insulinosensibilizadores","insulinosensibilizadores+insulinosecretores","Incretinas","Insulinas","Insulinas+insulinosensibilizadores"),using_na=T)
##
## absoluta relativa acumulada
## Sin farmacos 6 21.43 21.43
## insulinosensibilizadores 12 42.86 64.29
## insulinosensibilizadores+insulinosecretores 2 7.14 71.43
## Incretinas 1 3.57 75
## Insulinas 5 17.86 92.86
## Insulinas+insulinosensibilizadores 2 7.14 100
## Total 28 100.00
## intervalo confianza inferior
## Sin farmacos 8.3
## insulinosensibilizadores 24.46
## insulinosensibilizadores+insulinosecretores 0.88
## Incretinas 0.09
## Insulinas 6.06
## Insulinas+insulinosensibilizadores 0.88
## Total
## intervalo confianza inferior
## Sin farmacos 40.95
## insulinosensibilizadores 62.82
## insulinosensibilizadores+insulinosecretores 23.5
## Incretinas 18.35
## Insulinas 36.89
## Insulinas+insulinosensibilizadores 23.5
## Totaldislipemicos=data[which(data$dislipemias==2),]
freq(dislipemicos,dislipemicos$Dieta.c,c("No","Si"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## No 21 65.62 65.62 46.81
## Si 11 34.38 100 18.57
## Total 32 100.00
## intervalo confianza inferior
## No 81.43
## Si 53.19
## Totalfreq(dislipemicos,dislipemicos$Farmacos.c,c("Sin farmacos","Estatinas","Omega 3"),using_na=T)
##
## absoluta relativa acumulada intervalo confianza inferior
## Sin farmacos 11 34.38 34.38 18.57
## Estatinas 20 62.5 96.88 43.69
## Omega 3 1 3.12 100 0.08
## Total 32 100.00
## intervalo confianza inferior
## Sin farmacos 53.19
## Estatinas 78.9
## Omega 3 16.22
## Totalmeans(data,data$cinturaabdominal,data$sexo)##
## cinturaabdominal by sexo = 1 cinturaabdominal by sexo = 2
## n 54 16
## Media 100.87 103.5
## Mediana 97 109
## DS 20.52 13.94
## Min 69 77
## Max 170 120
## Rango 101 43
## Q1 90.25 91.5
## Q3 104.75 114.25
## RIQ 14.5 22.75
## Normal 0 0.091
## ICI 95.27 96.07
## ICS 106.47 110.93
## Resumen2 97 (90.25-104.75) 109 (91.5-114.25)
## Resumen1 Mediana[Q1-Q3]
## grupos Homocedastica
## Normalidad No
## Diferencias No
## p-value 0.231
## Metodo Wilcoxon
means(data,data$icc,data$sexo)##
## icc by sexo = 1 icc by sexo = 2
## n 52 16
## Media 0.93 1.01
## Mediana 0.91 0.99
## DS 0.12 0.12
## Min 0.72 0.82
## Max 1.34 1.28
## Rango 0.62 0.46
## Q1 0.86 0.96
## Q3 0.98 1.08
## RIQ 0.12 0.12
## Normal 0.001 0.646
## ICI 0.9 0.95
## ICS 0.96 1.08
## Resumen2 0.91 (0.86-0.98) 0.99 (0.96-1.08)
## Resumen1 Mediana[Q1-Q3]
## grupos Homocedastica
## Normalidad No
## Diferencias Si
## p-value 0.007
## Metodo Wilcoxon
means(data,data$hba1c)##
## hba1c
## n 22
## Media 7.69
## Mediana 7.05
## DS 2.22
## Min 5.3
## Max 13.1
## Rango 7.8
## Q1 5.82
## Q3 8.93
## RIQ 3.1
## Normal 0.013
## ICI 6.71
## ICS 8.67
## Resumen Mediana[Q1-Q3] -->7.05(5.82-8.93)
lineal(data$totalactividadfisica,data$hba1c,labely="Actividad Fisica",labelx="HbA1c")##
## Call:
## lm(formula = y ~ x)
##
## Coefficients:
## (Intercept) x
## 378.41 -13.85
##
## [1] "r2=0.014"
## [1] "p=0.777"