Aca se analiza como se distribuye la extension
La tabla de datos importada es LIBROSINBORDESINCISALES.xlsx Se toma en cuenta además la variable bewe donde se corrigio acumulando superificies de tipo 1 y 2 (en nivel de erosion)
library(survey)
##
## Attaching package: 'survey'
##
## The following object(s) are masked from 'package:graphics':
##
## dotchart
options(OutDec = ",")
library(car)
## Loading required package: MASS
## Loading required package: nnet
load("~/Dropbox/odontologia/maestria licet/octubre 2013/datos_licet_28102013.RData")
extension <- read.csv("sbi.csv")
str(extension)
## 'data.frame': 1155 obs. of 4 variables:
## $ ID : Factor w/ 1155 levels "","1500","1501",..: 2 3 4 5 6 7 8 9 10 11 ...
## $ bewes : int 14 0 6 0 15 14 10 8 2 0 ...
## $ ERODENTINASINBORDES: int 0 0 0 0 0 0 0 0 0 0 ...
## $ severas : int 0 0 0 0 0 0 0 0 0 0 ...
extension[1150:1155, ]
## ID bewes ERODENTINASINBORDES severas
## 1150 2747 0 0 0
## 1151 2748 0 0 0
## 1152 2749 0 0 0
## 1153 2750 0 0 0
## 1154 2751 0 0 0
## 1155 NA 47 NA
extension <- extension[1:1154, ]
diseniopost1$variables <- diseniopost1$variables[, 1:483]
head(data.frame(diseniopost1$variables[, 1:10], extension))
## Nro.sobre ID Escuela Tipo.Muestra Peso.cali1 wi.e1_fi1
## 1 1 1500 49 ma 8,295 9,556
## 2 1 1501 49 ma 7,162 9,556
## 3 1 1502 49 ma 8,295 9,556
## 4 1 1503 49 ma 7,162 9,556
## 5 1 1504 49 ma 8,295 9,556
## 6 1 1505 49 ma 8,295 9,556
## wi.e1_wi.e2_fi1_fi2 wi.e1_wi.e2_fi1_fi2_fi3 peso.cali2 Zona ID.1 bewes
## 1 5,063 4,701 18,16 1 1500 14
## 2 5,063 4,701 14,30 1 1501 0
## 3 5,063 4,701 18,16 1 1502 6
## 4 5,063 4,701 14,30 1 1503 0
## 5 5,063 4,701 18,16 1 1504 15
## 6 5,063 4,701 18,16 1 1505 14
## ERODENTINASINBORDES severas
## 1 0 0
## 2 0 0
## 3 0 0
## 4 0 0
## 5 0 0
## 6 0 0
colnames(extension)
## [1] "ID" "bewes" "ERODENTINASINBORDES"
## [4] "severas"
colnames(extension)[2] <- c("Bewes.ajust")
colnames(extension)[3] <- c("Erosinbord")
colnames(extension)[4] <- c("Severas")
diseniopost1$variables <- cbind(diseniopost1$variables, extension[, 2:4])
summary(diseniopost1$variables)[, 484:486]
## Bewes.ajust Erosinbord Severas
## "Min. : 0,00 " "Min. : 0,00 " "Min. : 0,00 "
## "1st Qu.: 0,00 " "1st Qu.: 0,00 " "1st Qu.: 0,00 "
## "Median : 4,00 " "Median : 0,00 " "Median : 0,00 "
## "Mean : 6,76 " "Mean : 1,41 " "Mean : 1,45 "
## "3rd Qu.:10,00 " "3rd Qu.: 0,00 " "3rd Qu.: 0,00 "
## "Max. :88,00 " "Max. :88,00 " "Max. :88,00 "
## NA NA NA
diseniopost1$variables$Bewes.ajust <- recode(diseniopost1$variables$Bewes.ajust,
"88=NA")
diseniopost1$variables$Erosinbord <- recode(diseniopost1$variables$Erosinbord,
"88=NA")
diseniopost1$variables$Severas <- recode(diseniopost1$variables$Severas, "88=NA")
table(diseniopost1$variables$Bewes.ajust)
##
## 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
## 539 7 7 8 36 12 50 38 88 37 50 48 46 32 37 26 23 12
## 18 19 20 21 22 23 24 26 27 29 31 32 35
## 6 6 5 4 7 2 3 1 2 1 1 1 1
table(diseniopost1$variables$Severas)
##
## 0 1 2 3 4
## 1089 23 16 2 6
Severidad.1 <- svymean(~Bewes.ajust, diseniopost1, na.rm = TRUE, deff = TRUE)
round(Severidad.1 * 1, 5)
## mean SE DEff
## Bewes.ajust 5,851 0,749 13,4
confint(Severidad.1)
## 2,5 % 97,5 %
## Bewes.ajust 4,382 7,319
Severidad.2 <- svymean(~Severas, diseniopost1, na.rm = TRUE, deff = TRUE)
round(Severidad.2 * 1, 5)
## mean SE DEff
## Severas 0,0831 0,0238 3,34
confint(Severidad.2)
## 2,5 % 97,5 %
## Severas 0,03634 0,1298
tabla.Bewes_ajust <- svycdf(~Bewes.ajust, diseniopost1)
tabla.Bewes_ajust$Bewes.ajust
## Step function
## Call: svycdf.default(formula = ~Bewes.ajust, design = diseniopost1)
## x[1:31] = 0, 1, 2, ..., 32, 35
## 32 plateau levels = 0, 0.48, 0.48, ..., 1, 1
plot(tabla.Bewes_ajust, cex.main = 1, 0, main = "Nivel de Bewes Ajustado")
abline(h = 0.4772, col = "red", lty = 2)
abline(v = 0, col = "red", lty = 2)
svyboxplot(Bewes.ajust ~ 1, diseniopost1, main = "Bewes ajustado en Extension sin bordes incisales",
col = "yellow", horizontal = TRUE, cex.main = 0.9, xlab = "Conteo de Bewes")
tabla.severidad <- svycdf(~Severas, diseniopost1)
tabla.severidad$Severas
## Step function
## Call: svycdf.default(formula = ~Severas, design = diseniopost1)
## x[1:5] = 0, 1, 2, 3, 4
## 6 plateau levels = 0, 0.96, 0.98, ..., 0.99, 1
plot(tabla.severidad, cex.main = 1, 0, main = "Nivel de Severidad")
abline(h = 0.95581, col = "red", lty = 2)
abline(v = 0, col = "red", lty = 2)
Bewes.ajust.6 <- svyby(~Bewes.ajust, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Bewes.ajust6[, 2:4] * 1, 2)
## Error: object 'Bewes.ajust6' not found
round(confint(Bewes.ajust6) * 1, 2)
## Error: object 'Bewes.ajust6' not found
Bewes.ajust.7 <- svyby(~Bewes.ajust, ~Nivel.Socieconomico.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Bewes.ajust7[, 2:4] * 1, 2)
## Error: object 'Bewes.ajust7' not found
round(confint(Bewes.ajust7) * 1, 2)
## Error: object 'Bewes.ajust7' not found
Bewes.ajust8 <- svyby(~Bewes.ajust, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Bewes.ajust8[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Incomplete Education 5,34 1,25 3,84
## 2-Basic 4,94 0,60 6,33
## 3-Medium 8,09 1,37 5,87
## 4-High 6,31 0,68 1,09
round(confint(Bewes.ajust8) * 1, 2)
## 2,5 % 97,5 %
## 1-Incomplete Education 2,88 7,79
## 2-Basic 3,76 6,12
## 3-Medium 5,41 10,76
## 4-High 4,97 7,65
Bewes.ajust8a <- svyby(~Bewes.ajust, ~Nive.Educativo.de.la.Madre2.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Bewes.ajust8a[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-ELEMENTARY SCHOOL 4,37 0,87 7,98
## 2-HIGH SCHOOL 6,10 0,75 6,85
## 3-COLLEGE-UNIVERSITY 7,07 1,12 3,79
round(confint(Bewes.ajust8a) * 1, 2)
## 2,5 % 97,5 %
## 1-ELEMENTARY SCHOOL 2,67 6,07
## 2-HIGH SCHOOL 4,62 7,58
## 3-COLLEGE-UNIVERSITY 4,89 9,26
Bewes.ajust9 <- svyby(~Bewes.ajust, ~Tipo.de.Escuela.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust9[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Public 6,49 0,56 1,70
## 2-Private 5,63 0,99 18,01
round(confint(Bewes.ajust9) * 1, 2)
## 2,5 % 97,5 %
## 1-Public 5,40 7,58
## 2-Private 3,69 7,57
Bewes.ajust10 <- svyby(~Bewes.ajust, ~FrCepDenti.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust10[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-1 vez al dia o menos 4,54 0,67 3,60
## 2- 2 veces al dia 5,71 0,59 3,29
## 3-3 o mas veces al dia 6,53 1,21 11,65
round(confint(Bewes.ajust10) * 1, 2)
## 2,5 % 97,5 %
## 1-1 vez al dia o menos 3,22 5,85
## 2- 2 veces al dia 4,56 6,86
## 3-3 o mas veces al dia 4,16 8,90
Bewes.ajust11 <- svyby(~Bewes.ajust, ~UsoDentifrico3.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust11[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Yes 5,92 0,75 13,07
## 2-No 1,73 0,77 1,23
round(confint(Bewes.ajust11) * 1, 2)
## 2,5 % 97,5 %
## 1-Yes 4,44 7,40
## 2-No 0,22 3,24
Bewes.ajust12 <- svyby(~Bewes.ajust, ~Consitencia_Cepillo.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Bewes.ajust12[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Soft 6,22 1,10 4,56
## 2-Medium 6,22 0,88 8,30
## 3-Hard 6,05 0,82 2,66
round(confint(Bewes.ajust12) * 1, 2)
## 2,5 % 97,5 %
## 1-Soft 4,07 8,38
## 2-Medium 4,49 7,95
## 3-Hard 4,44 7,65
Bewes.ajust13 <- svyby(~Bewes.ajust, ~IGS.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Bewes.ajust13[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Mas de 60 5,38 0,88 2,28
## 2-De 45 a 60 5,61 1,21 4,31
## 3 -Menos o igual a 45 5,93 0,79 11,57
round(confint(Bewes.ajust13) * 1, 2)
## 2,5 % 97,5 %
## 1-Mas de 60 3,66 7,10
## 2-De 45 a 60 3,24 7,97
## 3 -Menos o igual a 45 4,38 7,49
Bewes.ajust14 <- svyby(~Bewes.ajust, ~RefrCola.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Bewes.ajust14[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Nunca o raramente 6,29 0,82 2,62
## 2-Todos los dias 5,46 0,67 6,76
## 3-Mas de 3 veces al dia 6,74 1,42 9,84
round(confint(Bewes.ajust14) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 4,68 7,89
## 2-Todos los dias 4,15 6,76
## 3-Mas de 3 veces al dia 3,96 9,52
Bewes.ajust15 <- svyby(~Bewes.ajust, ~JugFrutas.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust15[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Nunca o raramente 5,44 0,99 1,48
## 2-Todos los dias 5,63 0,74 10,49
## 3-Mas de 3 veces al dia 7,17 1,02 3,40
round(confint(Bewes.ajust15) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 3,50 7,38
## 2-Todos los dias 4,18 7,08
## 3-Mas de 3 veces al dia 5,18 9,16
Bewes.ajust16 <- svyby(~Bewes.ajust, ~Bebidas_energizantes.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Bewes.ajust16[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Nunca o raramente 5,98 0,91 10,96
## 2-Todos los dias 5,77 0,72 4,87
## Mas de 3 veces al dia 2,58 1,52 1,10
round(confint(Bewes.ajust16) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 4,19 7,77
## 2-Todos los dias 4,37 7,18
## Mas de 3 veces al dia -0,40 5,56
Bewes.ajust17 <- svyby(~Bewes.ajust, ~Yogurt.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Bewes.ajust17[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Mas de 3 veces al dia 6,78 1,26 1,89
## 2-Todos los dias 5,81 0,91 13,00
## 3-Nunca o raramente 5,96 0,60 2,30
round(confint(Bewes.ajust17) * 1, 2)
## 2,5 % 97,5 %
## 1-Mas de 3 veces al dia 4,32 9,24
## 2-Todos los dias 4,02 7,60
## 3-Nunca o raramente 4,79 7,13
Bewes.ajust18 <- svyby(~Bewes.ajust, ~Bruxismo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Bewes.ajust18[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-No 5,84 0,81 12,52
## 2-Yes 5,65 0,73 2,27
round(confint(Bewes.ajust18) * 1, 2)
## 2,5 % 97,5 %
## 1-No 4,26 7,42
## 2-Yes 4,22 7,08
Bewes.ajust19 <- svyby(~Bewes.ajust, ~bedeportediario.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust19[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Water 6,01 0,71 7,98
## 2-Gatorade 6,51 1,78 7,25
## 3-other 6,63 1,14 1,35
round(confint(Bewes.ajust19) * 1, 2)
## 2,5 % 97,5 %
## 1-Water 4,62 7,41
## 2-Gatorade 3,03 9,99
## 3-other 4,40 8,85
Bewes.ajust21 <- svyby(~Bewes.ajust, ~BuchTragar.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust21[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-No 5,83 0,73 9,08
## 2-Yes 7,97 1,60 4,72
round(confint(Bewes.ajust21) * 1, 2)
## 2,5 % 97,5 %
## 1-No 4,40 7,25
## 2-Yes 4,83 11,12
Bewes.ajust22 <- svyby(~Bewes.ajust, ~FormBeber.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust22[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-Con sorbete 6,62 1,20 3,26
## 2-Por el pico 5,29 0,67 2,10
## 3-Con vaso 5,91 0,87 13,24
round(confint(Bewes.ajust22) * 1, 2)
## 2,5 % 97,5 %
## 1-Con sorbete 4,27 8,98
## 2-Por el pico 3,98 6,60
## 3-Con vaso 4,20 7,63
Bewes.ajust23 <- svyby(~Bewes.ajust, ~Natac2vec.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust23[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-No 6,05 0,81 10,63
## 2-Yes 6,22 0,94 2,69
round(confint(Bewes.ajust23) * 1, 2)
## 2,5 % 97,5 %
## 1-No 4,46 7,64
## 2-Yes 4,38 8,06
Bewes.ajust24 <- svyby(~Bewes.ajust, ~MedResp.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Bewes.ajust24[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-No 5,85 0,76 11,32
## 2-Yes 6,12 1,10 4,57
round(confint(Bewes.ajust24) * 1, 2)
## 2,5 % 97,5 %
## 1-No 4,36 7,35
## 2-Yes 3,95 8,28
Bewes.ajust25 <- svyby(~Bewes.ajust, ~AlterGastrica.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Bewes.ajust25[, 2:4] * 1, 2)
## Bewes.ajust se DEff.Bewes.ajust
## 1-No 5,95 0,80 13,73
## 2-Yes 4,70 0,87 1,61
round(confint(Bewes.ajust25) * 1, 2)
## 2,5 % 97,5 %
## 1-No 4,38 7,52
## 2-Yes 3,01 6,40
SEVERIDAD
Severidad.6 <- svyby(~Severas, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad.6[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-F 0,03 0,01 1,7
## 2-M 0,13 0,04 3,1
round(confint(Severidad.6) * 1, 2)
## 2,5 % 97,5 %
## 1-F 0,01 0,05
## 2-M 0,05 0,22
Severidad.7 <- svyby(~Severas, ~Nivel.Socieconomico.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad.7[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-BAJO 0,07 0,02 2,46
## 2-MEDIO 0,07 0,03 1,67
## 3-ALTO 0,13 0,06 3,14
round(confint(Severidad.7) * 1, 2)
## 2,5 % 97,5 %
## 1-BAJO 0,02 0,11
## 2-MEDIO 0,01 0,13
## 3-ALTO 0,01 0,25
Severidad8 <- svyby(~Severas, ~Nivel.Educativo.de.la.Madre1.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Severidad8[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Incomplete Education 0,26 0,12 2,46
## 2-Basic 0,07 0,03 2,92
## 3-Medium 0,06 0,02 0,78
## 4-High 0,11 0,05 1,41
round(confint(Severidad8) * 1, 2)
## 2,5 % 97,5 %
## 1-Incomplete Education 0,03 0,49
## 2-Basic 0,02 0,12
## 3-Medium 0,01 0,10
## 4-High 0,00 0,21
Severidad8a <- svyby(~Severas, ~Nive.Educativo.de.la.Madre2.rec, diseniopost1,
svymean, na.rm = TRUE, deff = TRUE)
round(Severidad8a[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-ELEMENTARY SCHOOL 0,11 0,04 2,10
## 2-HIGH SCHOOL 0,07 0,03 2,60
## 3-COLLEGE-UNIVERSITY 0,10 0,04 1,53
round(confint(Severidad8a) * 1, 2)
## 2,5 % 97,5 %
## 1-ELEMENTARY SCHOOL 0,03 0,18
## 2-HIGH SCHOOL 0,01 0,12
## 3-COLLEGE-UNIVERSITY 0,02 0,18
Severidad9 <- svyby(~Severas, ~Tipo.de.Escuela.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad9[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Public 0,06 0,04 1,61
## 2-Private 0,09 0,03 4,06
round(confint(Severidad9) * 1, 2)
## 2,5 % 97,5 %
## 1-Public -0,01 0,14
## 2-Private 0,03 0,15
Severidad10 <- svyby(~Severas, ~FrCepDenti.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad10[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-1 vez al dia o menos 0,07 0,03 1,44
## 2- 2 veces al dia 0,11 0,06 4,29
## 3-3 o mas veces al dia 0,08 0,02 1,42
round(confint(Severidad10) * 1, 2)
## 2,5 % 97,5 %
## 1-1 vez al dia o menos 0,01 0,13
## 2- 2 veces al dia -0,01 0,22
## 3-3 o mas veces al dia 0,04 0,11
Severidad11 <- svyby(~Severas, ~UsoDentifrico3.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad11[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Yes 0,08 0,02 3,50
## 2-No 0,18 0,15 0,96
round(confint(Severidad11) * 1, 2)
## 2,5 % 97,5 %
## 1-Yes 0,03 0,13
## 2-No -0,11 0,46
Severidad12 <- svyby(~Severas, ~Consitencia_Cepillo.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad12[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Soft 0,11 0,03 1,26
## 2-Medium 0,05 0,01 1,15
## 3-Hard 0,18 0,13 5,95
round(confint(Severidad12) * 1, 2)
## 2,5 % 97,5 %
## 1-Soft 0,05 0,17
## 2-Medium 0,02 0,07
## 3-Hard -0,07 0,42
Severidad13 <- svyby(~Severas, ~IGS.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad13[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Mas de 60 0,01 0,01 0,85
## 2-De 45 a 60 0,08 0,04 1,33
## 3 -Menos o igual a 45 0,09 0,03 3,51
round(confint(Severidad13) * 1, 2)
## 2,5 % 97,5 %
## 1-Mas de 60 -0,01 0,02
## 2-De 45 a 60 0,00 0,17
## 3 -Menos o igual a 45 0,03 0,15
Severidad14 <- svyby(~Severas, ~RefrCola.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad14[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Nunca o raramente 0,03 0,03 1,19
## 2-Todos los dias 0,08 0,03 2,58
## 3-Mas de 3 veces al dia 0,11 0,04 1,92
round(confint(Severidad14) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente -0,02 0,08
## 2-Todos los dias 0,03 0,14
## 3-Mas de 3 veces al dia 0,03 0,18
Severidad15 <- svyby(~Severas, ~JugFrutas.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad15[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Nunca o raramente 0,02 0,02 1,03
## 2-Todos los dias 0,08 0,03 3,30
## 3-Mas de 3 veces al dia 0,11 0,03 1,12
round(confint(Severidad15) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente -0,02 0,07
## 2-Todos los dias 0,03 0,13
## 3-Mas de 3 veces al dia 0,05 0,18
Severidad16 <- svyby(~Severas, ~Bebidas_energizantes.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad16[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Nunca o raramente 0,06 0,02 2,89
## 2-Todos los dias 0,09 0,03 2,14
## Mas de 3 veces al dia 0,51 0,28 0,94
round(confint(Severidad16) * 1, 2)
## 2,5 % 97,5 %
## 1-Nunca o raramente 0,01 0,10
## 2-Todos los dias 0,03 0,16
## Mas de 3 veces al dia -0,05 1,07
Severidad17 <- svyby(~Severas, ~Yogurt.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad17[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Mas de 3 veces al dia 0,17 0,09 0,94
## 2-Todos los dias 0,08 0,02 1,77
## 3-Nunca o raramente 0,08 0,04 2,16
round(confint(Severidad17) * 1, 2)
## 2,5 % 97,5 %
## 1-Mas de 3 veces al dia 0,00 0,34
## 2-Todos los dias 0,04 0,12
## 3-Nunca o raramente 0,01 0,15
Severidad18 <- svyby(~Severas, ~Bruxismo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad18[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 0,09 0,03 3,15
## 2-Yes 0,06 0,02 1,09
round(confint(Severidad18) * 1, 2)
## 2,5 % 97,5 %
## 1-No 0,04 0,14
## 2-Yes 0,01 0,11
Severidad19 <- svyby(~Severas, ~bedeportediario.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad19[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Water 0,08 0,03 4,28
## 2-Gatorade 0,12 0,06 2,69
## 3-other 0,02 0,02 0,54
round(confint(Severidad19) * 1, 2)
## 2,5 % 97,5 %
## 1-Water 0,02 0,15
## 2-Gatorade 0,00 0,24
## 3-other -0,02 0,07
Severidad21 <- svyby(~Severas, ~BuchTragar.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad21[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 0,08 0,03 4,09
## 2-Yes 0,11 0,03 0,82
round(confint(Severidad21) * 1, 2)
## 2,5 % 97,5 %
## 1-No 0,02 0,15
## 2-Yes 0,04 0,17
Severidad22 <- svyby(~Severas, ~FormBeber.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad22[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Con sorbete 0,06 0,04 0,87
## 2-Por el pico 0,04 0,02 1,10
## 3-Con vaso 0,09 0,03 3,33
round(confint(Severidad22) * 1, 2)
## 2,5 % 97,5 %
## 1-Con sorbete -0,01 0,14
## 2-Por el pico -0,01 0,09
## 3-Con vaso 0,04 0,15
Severidad23 <- svyby(~Severas, ~Natac2vec.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad23[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 0,09 0,03 3,84
## 2-Yes 0,07 0,02 1,15
round(confint(Severidad23) * 1, 2)
## 2,5 % 97,5 %
## 1-No 0,02 0,15
## 2-Yes 0,02 0,12
Severidad24 <- svyby(~Severas, ~MedResp.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad24[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 0,09 0,03 3,16
## 2-Yes 0,04 0,02 1,53
round(confint(Severidad24) * 1, 2)
## 2,5 % 97,5 %
## 1-No 0,04 0,14
## 2-Yes -0,01 0,09
Severidad25 <- svyby(~Severas, ~AlterGastrica.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad25[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No 0,09 0,03 3,24
## 2-Yes 0,00 0,00 NaN
round(confint(Severidad25) * 1, 2)
## 2,5 % 97,5 %
## 1-No 0,04 0,14
## 2-Yes 0,00 0,00