Aca se analiza como se distribuye la extension
La tabla de datos importada es LIBROSINBORDESINCISALES.xlsx
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) extension[1150:1155,]
# extension<-extension[1:1154,]
# head(data.frame(diseniopost1$variables[,1:10],extension))
# colnames(extension) colnames(extension)[2]<-c('Erosinbord')
# colnames(extension)[3]<-c('Severas')
# diseniopost1$variables<-cbind(diseniopost1$variables,extension[,2:3])
# summary(diseniopost1$variables)
# diseniopost1$variables$Erosinbord<-recode(diseniopost1$variables$Erosinbord,'88=NA')
# diseniopost1$variables$Severas<-recode(diseniopost1$variables$Severas,'88=NA')
table(diseniopost1$variables$Severas)
##
## 0 1 2 3 4
## 1089 23 16 2 6
Severidad.1 <- svymean(~Severas, diseniopost1, na.rm = TRUE, deff = TRUE)
round(Severidad.1 * 1, 5)
## mean SE DEff
## Severas 0,0831 0,0238 3,34
confint(Severidad.1)
## 2,5 % 97,5 %
## Severas 0,03634 0,1298
svyboxplot(Severas ~ 1, diseniopost1, main = "Severidad en Extension sin bordes incisales",
col = "purple", horizontal = TRUE, cex.main = 1)
Severidad.6 <- svyby(~Severas, ~Sexo.rec, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
Severidad.6
## Sexo.rec Severas se DEff.Severas
## 1-F 1-F 0,03153 0,01131 1,698
## 2-M 2-M 0,13168 0,04348 3,096
Severidad.7 <- svyby(~Severas, ~Nivel.Socieconomico.rec, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
Severidad.7
## Nivel.Socieconomico.rec Severas se DEff.Severas
## 1-BAJO 1-BAJO 0,06513 0,02395 2,457
## 2-MEDIO 2-MEDIO 0,06959 0,03023 1,668
## 3-ALTO 3-ALTO 0,12703 0,06046 3,140
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_b, diseniopost1, svymean,
na.rm = TRUE, deff = TRUE)
round(Severidad16[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-No diario 0,07 0,03 4,46
## 2-Diario 0,11 0,05 1,25
round(confint(Severidad16) * 1, 2)
## 2,5 % 97,5 %
## 1-No diario 0,02 0,12
## 2-Diario 0,01 0,21
Severidad17 <- svyby(~Severas, ~Yogurt_b, diseniopost1, svymean, na.rm = TRUE,
deff = TRUE)
round(Severidad17[, 2:4] * 1, 2)
## Severas se DEff.Severas
## 1-Diario 0,08 0,04 1,32
## 2-No diario 0,08 0,03 3,15
round(confint(Severidad17) * 1, 2)
## 2,5 % 97,5 %
## 1-Diario 0,00 0,15
## 2-No diario 0,03 0,13
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