datos.R
mauriciogutierrez
Thu Nov 10 14:45:26 2016
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.com## Loading required package: reshape2## Loading required package: stringrdata <- readexcel("/Users/mauriciogutierrez/Library/Mobile Documents/com~apple~CloudDocs/estadisticas/bermudez/datos2.xls","Hoja1")
library(stringr)
labeldicotomico=c("no","yes")
# Datos sobre el total
means(data,data$edad)##
## edad
## n 63
## Mean 51.76
## Median 52
## Sd 10.94
## Min 29
## Max 74
## Range 45
## Q1 43
## Q3 58
## Iqr 15
## Normality 0.275
## ICI 49.01
## SCI 54.52
## Summary Mean[ICI-ICS] --> 51.76(49.01-54.52)
data$imc=as.numeric(data$IMC)## Warning: NAs introduced by coercionmeans(data,data$imc)##
## imc
## n 62
## Mean 32.37
## Median 31.68
## Sd 5.81
## Min 20.39
## Max 46.02
## Range 25.63
## Q1 28.32
## Q3 35.52
## Iqr 7.2
## Normality 0.212
## ICI 30.89
## SCI 33.84
## Summary Mean[ICI-ICS] --> 32.37(30.89-33.84)
means(data,data$imc,data$sexo)##
## sexo = F sexo = M
## n 25 37
## Mean 30.03 33.95
## Median 29.37 33.49
## Sd 5.47 5.56
## Min 20.39 24.84
## Max 45.23 46.02
## Range 24.84 21.18
## Q1 28.03 29.47
## Q3 32.04 37.57
## Iqr 4.01 8.1
## Normal 0.197 0.354
## IIC 27.77 32.09
## ICS 32.28 35.8
## Summary2 30.03 (27.77-32.28) 33.95 (32.09-35.8)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences Yes
## p-value 0.008
## Method t-Student(Pooled)
means(data,data$Cintura,data$sexo)##
## sexo = F sexo = M
## n 25 38
## Mean 97.32 114.82
## Median 98 115
## Sd 12.33 13.69
## Min 64 90
## Max 120 145
## Range 56 55
## Q1 92 105
## Q3 104 123.75
## Iqr 12 18.75
## Normal 0.767 0.529
## IIC 92.23 110.31
## ICS 102.41 119.32
## Summary2 97.32 (92.23-102.41) 114.82 (110.31-119.32)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences Yes
## p-value <0.001
## Method t-Student(Pooled)
freq(data,data$sexo,using_na=T)
##
## absolute relative cumulative ici ics
## F 25 39.68 39.68 27.57 52.8
## M 38 60.32 100 47.2 72.43
## Total 63 100.00freq(data,data$SM,using_na=T,labeldicotomico)
##
## absolute relative cumulative ici ics
## no 6 9.52 9.52 3.58 19.59
## yes 57 90.48 100 80.41 96.42
## Total 63 100.00freq(data,data$HTA,using_na=T,labeldicotomico)
##
## absolute relative cumulative ici ics
## no 20 31.75 31.75 20.58 44.69
## yes 43 68.25 100 55.31 79.42
## Total 63 100.00means(data,data$PAS)##
## PAS
## n 63
## Mean 125.71
## Median 120
## Sd 16.29
## Min 100
## Max 170
## Range 70
## Q1 120
## Q3 130
## Iqr 10
## Normality 0
## ICI 121.61
## SCI 129.82
## Summary Median[Q1-Q3] -->120(120-130)
means(data,data$PAD)##
## PAD
## n 63
## Mean 76.43
## Median 80
## Sd 11.16
## Min 50
## Max 100
## Range 50
## Q1 70
## Q3 80
## Iqr 10
## Normality 0.003
## ICI 73.62
## SCI 79.24
## Summary Median[Q1-Q3] -->80(70-80)
means(data,data$PP)##
## PP
## n 62
## Mean 49.27
## Median 50
## Sd 11.05
## Min 20
## Max 80
## Range 60
## Q1 40
## Q3 55
## Iqr 15
## Normality 0.002
## ICI 46.47
## SCI 52.08
## Summary Median[Q1-Q3] -->50(40-55)
means(data,data$PAS,data$SM,labeldicotomico)##
## SM = no SM = yes
## n 6 57
## Mean 115.83 126.75
## Median 120 125
## Sd 11.14 16.46
## Min 100 100
## Max 130 170
## Range 30 70
## Q1 108.75 120
## Q3 120 130
## Iqr 11.25 10
## Normal 0.316 0.001
## IIC 104.14 122.39
## ICS 127.53 131.12
## Summary2 120 (108.75-120) 125 (120-130)
## Summary Median[Q1-Q3]
## Groups Homocedastics
## Normality No
## Differences No
## p-value 0.109
## Method Wilcoxon
means(data,data$PAD,data$SM,labeldicotomico)##
## SM = no SM = yes
## n 6 57
## Mean 74.17 76.67
## Median 75 80
## Sd 9.17 11.39
## Min 60 50
## Max 85 100
## Range 25 50
## Q1 70 70
## Q3 80 80
## Iqr 10 10
## Normal 0.566 0.005
## IIC 64.54 73.64
## ICS 83.79 79.69
## Summary2 75 (70-80) 80 (70-80)
## Summary Median[Q1-Q3]
## Groups Homocedastics
## Normality No
## Differences No
## p-value 0.654
## Method Wilcoxon
means(data,data$PP,data$SM,labeldicotomico)##
## SM = no SM = yes
## n 6 56
## Mean 41.67 50.09
## Median 40 50
## Sd 6.83 11.14
## Min 35 20
## Max 50 80
## Range 15 60
## Q1 36.25 40
## Q3 47.5 60
## Iqr 11.25 20
## Normal 0.093 0.003
## IIC 34.5 47.11
## ICS 48.84 53.07
## Summary2 40 (36.25-47.5) 50 (40-60)
## Summary Median[Q1-Q3]
## Groups Homocedastics
## Normality No
## Differences Yes
## p-value 0.038
## Method Wilcoxon
freq(data,data$Glicemia,using_na=T,c("Normoglucemic","IFG","Diabetes"))
##
## absolute relative cumulative ici ics
## Normoglucemic 30 47.62 47.62 34.88 60.59
## IFG 25 39.68 87.3 27.57 52.8
## Diabetes 8 12.7 100 5.65 23.5
## Total 63 100.00freq(data,data$sedentarismo,using_na=T,labeldicotomico)
##
## absolute relative cumulative ici ics
## no 25 39.68 39.68 27.57 52.8
## yes 38 60.32 100 47.2 72.43
## Total 63 100.00freq(data,data$SAOS,using_na=T,labeldicotomico)
##
## absolute relative cumulative ici ics
## no 56 88.89 88.89 78.44 95.41
## yes 7 11.11 100 4.59 21.56
## Total 63 100.00freq(data,data$Diabetes,using_na=T,labeldicotomico)
##
## absolute relative cumulative ici ics
## no 56 88.89 88.89 78.44 95.41
## yes 7 11.11 100 4.59 21.56
## Total 63 100.00freq(data,data$Estatinas,using_na=T,labeldicotomico)
##
## absolute relative cumulative ici ics
## no 38 60.32 60.32 47.2 72.43
## yes 25 39.68 100 27.57 52.8
## Total 63 100.00data$score=as.numeric(str_replace(data$Score,"%",""))
means(data,data$score)##
## score
## n 63
## Mean 1.53
## Median 1
## Sd 1.98
## Min 0
## Max 12
## Range 12
## Q1 0.5
## Q3 2
## Iqr 1.5
## Normality 0
## ICI 1.03
## SCI 2.03
## Summary Median[Q1-Q3] -->1(0.5-2)
data$framminghan=as.numeric(str_replace(data$Framminghan,"%",""))
means(data,data$framminghan)##
## framminghan
## n 63
## Mean 7.38
## Median 5
## Sd 8.51
## Min 0.5
## Max 48
## Range 47.5
## Q1 1.5
## Q3 10
## Iqr 8.5
## Normality 0
## ICI 5.24
## SCI 9.53
## Summary Median[Q1-Q3] -->5(1.5-10)
data$capacidad_aerobica=as.numeric(data$capacidad_aerobica)
means(data,data$capacidad_aerobica)##
## capacidad_aerobica
## n 63
## Mean 72.32
## Median 71.9
## Sd 23.02
## Min 17.5
## Max 150.41
## Range 132.91
## Q1 54.12
## Q3 91.1
## Iqr 36.98
## Normality 0.102
## ICI 66.52
## SCI 78.12
## Summary Mean[ICI-ICS] --> 72.32(66.52-78.12)
data$METS=as.numeric(data$METS)## Warning: NAs introduced by coercionmeans(data,data$METS)##
## METS
## n 61
## Mean 5.24
## Median 5.05
## Sd 2.07
## Min 2.26
## Max 16.9
## Range 14.64
## Q1 3.81
## Q3 6.25
## Iqr 2.44
## Normality 0
## ICI 4.71
## SCI 5.77
## Summary Median[Q1-Q3] -->5.05(3.81-6.25)
means(data,data$HDL,data$sexo)##
## sexo = F sexo = M
## n 24 38
## Mean 51.12 45.95
## Median 53 45.5
## Sd 10.6 10.37
## Min 31 29
## Max 71 72
## Range 40 43
## Q1 45.5 39
## Q3 55.75 50
## Iqr 10.25 11
## Normal 0.396 0.048
## IIC 46.65 42.54
## ICS 55.6 49.36
## Summary2 53 (45.5-55.75) 45.5 (39-50)
## Summary Median[Q1-Q3]
## Groups Homocedastics
## Normality No
## Differences Yes
## p-value 0.036
## Method Kruskal-Wallis
means(data,data$LDL)##
## LDL
## n 62
## Mean 120.35
## Median 117.5
## Sd 42.48
## Min 40
## Max 236
## Range 196
## Q1 93.75
## Q3 148
## Iqr 54.25
## Normality 0.594
## ICI 109.57
## SCI 131.14
## Summary Mean[ICI-ICS] --> 120.35(109.57-131.14)
means(data,data$TG)##
## TG
## n 62
## Mean 157.94
## Median 137.5
## Sd 94.75
## Min 33
## Max 555
## Range 522
## Q1 93.75
## Q3 188.5
## Iqr 94.75
## Normality 0
## ICI 133.87
## SCI 182
## Summary Median[Q1-Q3] -->137.5(93.75-188.5)
# Comparaciones entre Normo e Hipertensos
means(data,data$edad,data$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 20 43
## Mean 46.3 54.3
## Median 46 54
## Sd 7.42 11.45
## Min 31 29
## Max 65 74
## Range 34 45
## Q1 42.75 45.5
## Q3 47.75 61.5
## Iqr 5 16
## Normal 0.409 0.345
## IIC 42.83 50.78
## ICS 49.77 57.83
## Summary2 46 (42.75-47.75) 54 (45.5-61.5)
## Summary Median[Q1-Q3]
## Groups Heteroscedastic
## Normality Yes
## Differences Yes
## p-value 0.002
## Method t-Student(Satterhwaite)
tables(data,data$HTA,data$sexo,labeldicotomico)##
## sexo
## HTA F M Sum
## no 11 9 20
## yes 14 29 43
## Sum 25 38 63
##
## by Row (%)
## F M
## no 55.00 45.00
## yes 32.55 67.44
## 39.68 60.31
##
## by Col (%)
## F M
## no 44 23.68 31.74
## yes 56 76.31 68.25
##
## Chi2 Pearson p = 0.09
means(data,data$imc,data$HTA)##
## HTA = 0 HTA = 1
## n 20 42
## Mean 30.73 33.15
## Median 29.29 33.02
## Sd 6.07 5.59
## Min 21.63 20.39
## Max 42.93 46.02
## Range 21.3 25.63
## Q1 26.93 29.37
## Q3 34.51 35.52
## Iqr 7.58 6.15
## Normal 0.249 0.174
## IIC 27.89 31.4
## ICS 33.57 34.89
## Summary2 30.73 (27.89-33.57) 33.15 (31.4-34.89)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences No
## p-value 0.126
## Method t-Student(Pooled)
tables(data,data$HTA,data$SM,labeldicotomico,labeldicotomico)##
## SM
## HTA no yes Sum
## no 3 17 20
## yes 3 40 43
## Sum 6 57 63
##
## by Row (%)
## no yes
## no 15.00 85.00
## yes 6.97 93.02
## 9.52 90.47
##
## by Col (%)
## no yes
## no 50 29.82 31.74
## yes 50 70.17 68.25
##
## Fischer Test p = 0.372
tables(data,data$HTA,data$Diabetes,labeldicotomico,labeldicotomico)##
## Diabetes
## HTA no yes Sum
## no 19 1 20
## yes 37 6 43
## Sum 56 7 63
##
## by Row (%)
## no yes
## no 95.00 5.00
## yes 86.04 13.95
## 88.88 11.11
##
## by Col (%)
## no yes
## no 33.92 14.28 31.74
## yes 66.07 85.71 68.25
##
## Fischer Test p = 0.415
tables(data,data$HTA,data$sedentarismo,labeldicotomico,labeldicotomico)##
## sedentarismo
## HTA no yes Sum
## no 11 9 20
## yes 14 29 43
## Sum 25 38 63
##
## by Row (%)
## no yes
## no 55.00 45.00
## yes 32.55 67.44
## 39.68 60.31
##
## by Col (%)
## no yes
## no 44 23.68 31.74
## yes 56 76.31 68.25
##
## Chi2 Pearson p = 0.09
tables(data,data$HTA,data$Estatinas,labeldicotomico,labeldicotomico)##
## Estatinas
## HTA no yes Sum
## no 14 6 20
## yes 24 19 43
## Sum 38 25 63
##
## by Row (%)
## no yes
## no 70.00 30.00
## yes 55.81 44.18
## 60.31 39.68
##
## by Col (%)
## no yes
## no 36.84 24 31.74
## yes 63.15 76 68.25
##
## Chi2 Pearson p = 0.284
means(data,data$score,data$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 20 43
## Mean 0.57 1.98
## Median 0.5 1
## Sd 0.44 2.25
## Min 0 0.5
## Max 2 12
## Range 2 11.5
## Q1 0.5 0.5
## Q3 0.5 2
## Iqr 0 1.5
## Normal 0 0
## IIC 0.37 1.28
## ICS 0.78 2.67
## Summary2 0.5 (0.5-0.5) 1 (0.5-2)
## Summary Median[Q1-Q3]
## Groups Heteroscedastic
## Normality No
## Differences No
## p-value <0.001
## Method Kruskal-Wallis
means(data,data$VDMF,data$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 20 43
## Mean 11.55 7.7
## Median 10 8
## Sd 4.77 3.97
## Min 4 0
## Max 23 16
## Range 19 16
## Q1 9 5
## Q3 13.25 10
## Iqr 4.25 5
## Normal 0.015 0.572
## IIC 9.32 6.48
## ICS 13.78 8.92
## Summary2 10 (9-13.25) 8 (5-10)
## Summary Median[Q1-Q3]
## Groups Homocedastics
## Normality No
## Differences Yes
## p-value 0.002
## Method Kruskal-Wallis
means(data,data$capacidad_aerobica,data$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 20 43
## Mean 76.89 70.2
## Median 79.7 63.2
## Sd 16.15 25.49
## Min 45.6 17.5
## Max 98 150.41
## Range 52.4 132.91
## Q1 71.68 53.6
## Q3 88.48 91.19
## Iqr 16.8 37.59
## Normal 0.037 0.091
## IIC 69.33 62.35
## ICS 84.44 78.04
## Summary2 79.7 (71.68-88.48) 63.2 (53.6-91.19)
## Summary Median[Q1-Q3]
## Groups Heteroscedastic
## Normality No
## Differences No
## p-value 0.243
## Method Kruskal-Wallis
means(data,data$LDL,data$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 19 43
## Mean 125.95 117.88
## Median 122 116
## Sd 48.72 39.79
## Min 40 44
## Max 236 221
## Range 196 177
## Q1 89 94.5
## Q3 153.5 146
## Iqr 64.5 51.5
## Normal 0.96 0.652
## IIC 102.47 105.64
## ICS 149.43 130.13
## Summary2 125.95 (102.47-149.43) 117.88 (105.64-130.13)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences No
## p-value 0.495
## Method t-Student(Pooled)
means(data,data$HDL,data$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 19 43
## Mean 47.89 47.98
## Median 53 46
## Sd 11 10.66
## Min 29 31
## Max 69 72
## Range 40 41
## Q1 37 39
## Q3 53.5 53.5
## Iqr 16.5 14.5
## Normal 0.254 0.061
## IIC 42.59 44.69
## ICS 53.2 51.26
## Summary2 47.89 (42.59-53.2) 47.98 (44.69-51.26)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences No
## p-value 0.978
## Method t-Student(Pooled)
means(data,data$TG,data$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 19 43
## Mean 140.21 165.77
## Median 134 138
## Sd 65.39 104.89
## Min 33 37
## Max 284 555
## Range 251 518
## Q1 89 98.5
## Q3 180 197.5
## Iqr 91 99
## Normal 0.918 0
## IIC 108.69 133.49
## ICS 171.73 198.05
## Summary2 134 (89-180) 138 (98.5-197.5)
## Summary Median[Q1-Q3]
## Groups Heteroscedastic
## Normality No
## Differences No
## p-value 0.598
## Method Kruskal-Wallis
# Comparacion de la VDMF segun SM
means(data,data$VDMF,data$SM,labeldicotomico)##
## SM = no SM = yes
## n 6 57
## Mean 12.67 8.53
## Median 10.5 8
## Sd 5.43 4.34
## Min 8 0
## Max 23 22
## Range 15 22
## Q1 10 6
## Q3 13.25 10
## Iqr 3.25 4
## Normal 0.059 0.263
## IIC 6.97 7.37
## ICS 18.36 9.68
## Summary2 12.67 (6.97-18.36) 8.53 (7.37-9.68)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences Yes
## p-value 0.034
## Method t-Student(Pooled)
# Comparaciones solo en grupo SM entre HTA y no HTA
soloSM=data[which(data$SM==1),]
means(soloSM,soloSM$edad,soloSM$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 17 40
## Mean 46.53 54.62
## Median 46 55.5
## Sd 8.01 11.79
## Min 31 29
## Max 65 74
## Range 34 45
## Q1 43 43.75
## Q3 50 62
## Iqr 7 18.25
## Normal 0.803 0.229
## IIC 42.41 50.86
## ICS 50.65 58.39
## Summary2 46.53 (42.41-50.65) 54.62 (50.86-58.39)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences Yes
## p-value 0.012
## Method t-Student(Pooled)
tables(soloSM,soloSM$HTA,soloSM$sexo,labeldicotomico)##
## sexo
## HTA F M Sum
## no 8 9 17
## yes 12 28 40
## Sum 20 37 57
##
## by Row (%)
## F M
## no 47.05 52.94
## yes 30.00 70.00
## 35.08 64.91
##
## by Col (%)
## F M
## no 40 24.32 29.82
## yes 60 75.67 70.17
##
## Chi2 Pearson p = 0.217
means(soloSM,soloSM$imc,soloSM$HTA)##
## HTA = 0 HTA = 1
## n 17 39
## Mean 31.7 33.48
## Median 30.34 33.04
## Sd 6.01 5.67
## Min 21.63 20.39
## Max 42.93 46.02
## Range 21.3 25.63
## Q1 28.07 29.59
## Q3 37.57 35.78
## Iqr 9.5 6.19
## Normal 0.521 0.367
## IIC 28.61 31.64
## ICS 34.79 35.31
## Summary2 31.7 (28.61-34.79) 33.48 (31.64-35.31)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences No
## p-value 0.294
## Method t-Student(Pooled)
tables(soloSM,soloSM$HTA,soloSM$Diabetes,labeldicotomico,labeldicotomico)##
## Diabetes
## HTA no yes Sum
## no 16 1 17
## yes 34 6 40
## Sum 50 7 57
##
## by Row (%)
## no yes
## no 94.11 5.88
## yes 85.00 15.00
## 87.71 12.28
##
## by Col (%)
## no yes
## no 32 14.28 29.82
## yes 68 85.71 70.17
##
## Fischer Test p = 0.662
tables(soloSM,soloSM$HTA,soloSM$sedentarismo,labeldicotomico,labeldicotomico)##
## sedentarismo
## HTA no yes Sum
## no 8 9 17
## yes 12 28 40
## Sum 20 37 57
##
## by Row (%)
## no yes
## no 47.05 52.94
## yes 30.00 70.00
## 35.08 64.91
##
## by Col (%)
## no yes
## no 40 24.32 29.82
## yes 60 75.67 70.17
##
## Chi2 Pearson p = 0.217
tables(soloSM,soloSM$HTA,soloSM$Estatinas,labeldicotomico,labeldicotomico)##
## Estatinas
## HTA no yes Sum
## no 11 6 17
## yes 22 18 40
## Sum 33 24 57
##
## by Row (%)
## no yes
## no 64.70 35.29
## yes 55.00 45.00
## 57.89 42.10
##
## by Col (%)
## no yes
## no 33.33 25 29.82
## yes 66.66 75 70.17
##
## Chi2 Pearson p = 0.497
means(soloSM,soloSM$score,soloSM$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 17 40
## Mean 0.59 2.08
## Median 0.5 1
## Sd 0.48 2.3
## Min 0 0.5
## Max 2 12
## Range 2 11.5
## Q1 0.5 0.5
## Q3 0.5 2.25
## Iqr 0 1.75
## Normal 0.001 0
## IIC 0.34 1.34
## ICS 0.83 2.81
## Summary2 0.5 (0.5-0.5) 1 (0.5-2.25)
## Summary Median[Q1-Q3]
## Groups Heteroscedastic
## Normality No
## Differences Yes
## p-value <0.001
## Method Kruskal-Wallis
means(soloSM,soloSM$VDMF,soloSM$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 17 40
## Mean 10.82 7.55
## Median 10 7.5
## Sd 4.22 4.06
## Min 4 0
## Max 22 16
## Range 18 16
## Q1 9 5
## Q3 13 9.25
## Iqr 4 4.25
## Normal 0.046 0.388
## IIC 8.66 6.25
## ICS 12.99 8.85
## Summary2 10 (9-13) 7.5 (5-9.25)
## Summary Median[Q1-Q3]
## Groups Homocedastics
## Normality No
## Differences Yes
## p-value 0.006
## Method Kruskal-Wallis
means(soloSM,soloSM$capacidad_aerobica,soloSM$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 17 40
## Mean 74.16 69.82
## Median 78.4 63.05
## Sd 15.89 25.94
## Min 45.6 17.5
## Max 96 150.41
## Range 50.4 132.91
## Q1 71 53.58
## Q3 81.6 90.95
## Iqr 10.6 37.37
## Normal 0.062 0.109
## IIC 66 61.52
## ICS 82.33 78.11
## Summary2 78.4 (71-81.6) 63.05 (53.58-90.95)
## Summary Median[Q1-Q3]
## Groups Heteroscedastic
## Normality Yes
## Differences No
## p-value 0.443
## Method t-Student(Satterhwaite)
means(soloSM,soloSM$LDL,soloSM$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 17 40
## Mean 127.41 118.92
## Median 122 116.5
## Sd 50 40.3
## Min 40 44
## Max 236 221
## Range 196 177
## Q1 99 97.5
## Q3 159 147.25
## Iqr 60 49.75
## Normal 0.989 0.598
## IIC 101.7 106.04
## ICS 153.12 131.81
## Summary2 127.41 (101.7-153.12) 118.92 (106.04-131.81)
## Summary Mean[ICI-ICS]
## Groups Homocedastics
## Normality Yes
## Differences No
## p-value 0.502
## Method t-Student(Pooled)
means(soloSM,soloSM$HDL,soloSM$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 17 40
## Mean 45.76 47.35
## Median 50 46
## Sd 9.44 10.61
## Min 29 31
## Max 58 72
## Range 29 41
## Q1 37 39
## Q3 53 53
## Iqr 16 14
## Normal 0.025 0.033
## IIC 40.91 43.96
## ICS 50.62 50.74
## Summary2 50 (37-53) 46 (39-53)
## Summary Median[Q1-Q3]
## Groups Homocedastics
## Normality No
## Differences No
## p-value 0.903
## Method Kruskal-Wallis
means(soloSM,soloSM$TG,soloSM$HTA,labeldicotomico)##
## HTA = no HTA = yes
## n 17 40
## Mean 147.24 169.55
## Median 140 141.5
## Sd 65.67 107.43
## Min 33 37
## Max 284 555
## Range 251 518
## Q1 116 99.75
## Q3 180 206.5
## Iqr 64 106.75
## Normal 0.955 0
## IIC 113.47 135.19
## ICS 181 203.91
## Summary2 140 (116-180) 141.5 (99.75-206.5)
## Summary Median[Q1-Q3]
## Groups Heteroscedastic
## Normality No
## Differences No
## p-value 0.807
## Method Kruskal-Wallis
# Para ver el peso de cada factor sobre la VDMF hacemos una regresion lineal multiple
reg=lm(data$VDMF ~ data$HDL + data$TG + data$HTA + data$Glicemia + data$Cintura + data$sexo + data$edad + data$LDL)
summary(reg)##
## Call:
## lm(formula = data$VDMF ~ data$HDL + data$TG + data$HTA + data$Glicemia +
## data$Cintura + data$sexo + data$edad + data$LDL)
##
## Residuals:
## Min 1Q Median 3Q Max
## -7.0918 -2.1430 -0.3988 2.5138 7.2188
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 25.391661 5.547312 4.577 2.88e-05 ***
## data$HDL 0.012601 0.051935 0.243 0.8092
## data$TG -0.007578 0.005813 -1.304 0.1980
## data$HTA -3.024003 1.180588 -2.561 0.0133 *
## data$Glicemia -0.648578 0.740941 -0.875 0.3853
## data$Cintura -0.091298 0.038773 -2.355 0.0223 *
## data$sexoM -1.567359 1.280170 -1.224 0.2262
## data$edad 0.002406 0.052030 0.046 0.9633
## data$LDL -0.022007 0.012277 -1.793 0.0787 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 3.842 on 53 degrees of freedom
## (1 observation deleted due to missingness)
## Multiple R-squared: 0.3962, Adjusted R-squared: 0.3051
## F-statistic: 4.347 on 8 and 53 DF, p-value: 0.0004435reg=lm(data$VDMF ~ data$PAS + data$PAD + data$PP)
summary(reg)##
## Call:
## lm(formula = data$VDMF ~ data$PAS + data$PAD + data$PP)
##
## Residuals:
## Min 1Q Median 3Q Max
## -9.7426 -1.7134 -0.1855 1.7646 11.4350
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 24.32165 4.37452 5.560 7.16e-07 ***
## data$PAS 0.07476 0.43411 0.172 0.864
## data$PAD -0.18826 0.43391 -0.434 0.666
## data$PP -0.21225 0.42592 -0.498 0.620
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 4.188 on 58 degrees of freedom
## (1 observation deleted due to missingness)
## Multiple R-squared: 0.2052, Adjusted R-squared: 0.1641
## F-statistic: 4.991 on 3 and 58 DF, p-value: 0.003791