las prevalencias de las caras oclusales
load("~/Dropbox/odontologia/maestria licet/octubre 2013/datos_licet_10112013.RData")
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
prevalencia_cara_oclusal <- read.csv("prevalencia_cara_oclusal.csv")
head(prevalencia_cara_oclusal)
## ID ERos17.O ERos16.O ERos15.O ERos14.O ERos13.I ERos12.I ERos11.I
## 1 1500 NA 1 NA NA NA 0 0
## 2 1501 NA 0 NA NA NA 0 0
## 3 1502 NA 0 NA NA NA 0 0
## 4 1503 NA 0 NA NA NA 0 0
## 5 1504 NA 0 NA NA NA 0 0
## 6 1505 NA 1 NA NA NA 0 0
## ERos21.I ERos22.I ERos23.I ERos24.O ERos25.O ERos26.O ERos27.O ERos37.O
## 1 0 0 NA NA NA 1 NA NA
## 2 0 0 NA NA NA 0 NA NA
## 3 0 0 NA NA NA 0 NA NA
## 4 0 0 NA NA NA 0 NA NA
## 5 0 0 NA NA NA 0 NA NA
## 6 0 0 NA NA NA 1 NA NA
## ERos36.O ERos35.O ERos34.O ERos33.I ERos32.I ERos31.I ERos41.I ERos42.I
## 1 0 NA NA NA 0 0 0 0
## 2 0 NA NA NA 0 0 0 0
## 3 0 NA NA NA 0 0 0 0
## 4 0 NA NA NA 0 0 0 0
## 5 1 NA NA NA 0 0 0 0
## 6 1 NA NA NA 0 0 0 0
## ERos43.I ERos44.O ERos45.O ERos46.O ERos47.O pOeros17 pOeros16 pOeros15
## 1 NA NA NA 0 NA NA 1 NA
## 2 NA NA NA 0 NA NA 0 NA
## 3 NA NA NA 0 NA NA 0 NA
## 4 NA NA NA 0 NA NA 0 NA
## 5 NA NA NA 1 NA NA 0 NA
## 6 NA NA NA 1 NA NA 1 NA
## pOeros14 pOeros13 pOeros12 pOeros11 pOeros21 pOeros22 pOeros23 pOeros24
## 1 NA NA 0 0 0 0 NA NA
## 2 NA NA 0 0 0 0 NA NA
## 3 NA NA 0 0 0 0 NA NA
## 4 NA NA 0 0 0 0 NA NA
## 5 NA NA 0 0 0 0 NA NA
## 6 NA NA 0 0 0 0 NA NA
## pOeros25 pOeros26 pOeros27 pOeros37 pOeros36 pOeros35 pOeros34 pOeros33
## 1 NA 1 NA NA 0 NA NA NA
## 2 NA 0 NA NA 0 NA NA NA
## 3 NA 0 NA NA 0 NA NA NA
## 4 NA 0 NA NA 0 NA NA NA
## 5 NA 0 NA NA 1 NA NA NA
## 6 NA 1 NA NA 1 NA NA NA
## pOeros32 pOeros31 pOeros41 pOeros42 pOeros43 pOeros44 pOeros45 pOeros46
## 1 0 0 0 0 NA NA NA 0
## 2 0 0 0 0 NA NA NA 0
## 3 0 0 0 0 NA NA NA 0
## 4 0 0 0 0 NA NA NA 0
## 5 0 0 0 0 NA NA NA 1
## 6 0 0 0 0 NA NA NA 1
## pOeros47
## 1 NA
## 2 NA
## 3 NA
## 4 NA
## 5 NA
## 6 NA
summary(diseniopost1$variables[, 572:599])
## pOeros17 pOeros16 pOeros15 pOeros14 pOeros13
## Min. :0 Min. :0,000 Min. :0,00 Min. :0,00 Min. :0
## 1st Qu.:0 1st Qu.:0,000 1st Qu.:0,00 1st Qu.:0,00 1st Qu.:0
## Median :0 Median :0,000 Median :0,00 Median :0,00 Median :0
## Mean :0 Mean :0,246 Mean :0,03 Mean :0,11 Mean :0
## 3rd Qu.:0 3rd Qu.:0,000 3rd Qu.:0,00 3rd Qu.:0,00 3rd Qu.:0
## Max. :1 Max. :1,000 Max. :1,00 Max. :1,00 Max. :0
## NA's :119 NA's :18 NA's :119 NA's :119 NA's :119
## pOeros12 pOeros11 pOeros21 pOeros22 pOeros23
## Min. :0 Min. :0 Min. :0 Min. :0,000 Min. :0
## 1st Qu.:0 1st Qu.:0 1st Qu.:0 1st Qu.:0,000 1st Qu.:0
## Median :0 Median :0 Median :0 Median :0,000 Median :0
## Mean :0 Mean :0 Mean :0 Mean :0,001 Mean :0
## 3rd Qu.:0 3rd Qu.:0 3rd Qu.:0 3rd Qu.:0,000 3rd Qu.:0
## Max. :0 Max. :0 Max. :0 Max. :1,000 Max. :0
## NA's :18 NA's :18 NA's :18 NA's :18 NA's :120
## pOeros24 pOeros25 pOeros26 pOeros27 pOeros37
## Min. :0,00 Min. :0,00 Min. :0,00 Min. :0 Min. :0
## 1st Qu.:0,00 1st Qu.:0,00 1st Qu.:0,00 1st Qu.:0 1st Qu.:0
## Median :0,00 Median :0,00 Median :0,00 Median :0 Median :0
## Mean :0,09 Mean :0,03 Mean :0,21 Mean :0 Mean :0
## 3rd Qu.:0,00 3rd Qu.:0,00 3rd Qu.:0,00 3rd Qu.:0 3rd Qu.:0
## Max. :1,00 Max. :1,00 Max. :1,00 Max. :1 Max. :1
## NA's :119 NA's :119 NA's :18 NA's :119 NA's :119
## pOeros36 pOeros35 pOeros34 pOeros33 pOeros32
## Min. :0,00 Min. :0,00 Min. :0,00 Min. :0 Min. :0
## 1st Qu.:0,00 1st Qu.:0,00 1st Qu.:0,00 1st Qu.:0 1st Qu.:0
## Median :0,00 Median :0,00 Median :0,00 Median :0 Median :0
## Mean :0,18 Mean :0,01 Mean :0,01 Mean :0 Mean :0
## 3rd Qu.:0,00 3rd Qu.:0,00 3rd Qu.:0,00 3rd Qu.:0 3rd Qu.:0
## Max. :1,00 Max. :1,00 Max. :1,00 Max. :0 Max. :0
## NA's :18 NA's :119 NA's :119 NA's :119 NA's :18
## pOeros31 pOeros41 pOeros42 pOeros43 pOeros44
## Min. :0 Min. :0 Min. :0,000 Min. :0 Min. :0,00
## 1st Qu.:0 1st Qu.:0 1st Qu.:0,000 1st Qu.:0 1st Qu.:0,00
## Median :0 Median :0 Median :0,000 Median :0 Median :0,00
## Mean :0 Mean :0 Mean :0,026 Mean :0 Mean :0,02
## 3rd Qu.:0 3rd Qu.:0 3rd Qu.:0,000 3rd Qu.:0 3rd Qu.:0,00
## Max. :0 Max. :0 Max. :1,000 Max. :0 Max. :1,00
## NA's :18 NA's :18 NA's :18 NA's :119 NA's :119
## pOeros45 pOeros46 pOeros47
## Min. :0,00 Min. :0,000 Min. :0,00
## 1st Qu.:0,00 1st Qu.:0,000 1st Qu.:0,00
## Median :0,00 Median :0,000 Median :0,00
## Mean :0,01 Mean :0,206 Mean :0,01
## 3rd Qu.:0,00 3rd Qu.:0,000 3rd Qu.:0,00
## Max. :1,00 Max. :1,000 Max. :1,00
## NA's :119 NA's :18 NA's :119
diseniopost1$variables <- cbind(diseniopost1$variables, prevalencia_cara_oclusal[,
30:57])
po <- names(diseniopost1$variables)
po[572:599]
## [1] "pOeros17" "pOeros16" "pOeros15" "pOeros14" "pOeros13" "pOeros12"
## [7] "pOeros11" "pOeros21" "pOeros22" "pOeros23" "pOeros24" "pOeros25"
## [13] "pOeros26" "pOeros27" "pOeros37" "pOeros36" "pOeros35" "pOeros34"
## [19] "pOeros33" "pOeros32" "pOeros31" "pOeros41" "pOeros42" "pOeros43"
## [25] "pOeros44" "pOeros45" "pOeros46" "pOeros47"
pOeros17 <- svymean(~factor(pOeros17), diseniopost1, na.rm = TRUE)
pOeros17 * 100
## mean SE
## factor(pOeros17)0 99,9077 0
## factor(pOeros17)1 0,0923 0
confint(pOeros17) * 100
## 2,5 % 97,5 %
## factor(pOeros17)0 99,78893 100,0264
## factor(pOeros17)1 -0,02643 0,2111
pOeros16 <- svymean(~factor(pOeros16), diseniopost1, na.rm = TRUE)
pOeros16 * 100
## mean SE
## factor(pOeros16)0 77,6 0,03
## factor(pOeros16)1 22,4 0,03
confint(pOeros16) * 100
## 2,5 % 97,5 %
## factor(pOeros16)0 72,56 82,65
## factor(pOeros16)1 17,35 27,44
pOeros15 <- svymean(~factor(pOeros15), diseniopost1, na.rm = TRUE)
pOeros15 * 100
## mean SE
## factor(pOeros15)0 96,25 0,01
## factor(pOeros15)1 3,75 0,01
confint(pOeros15) * 100
## 2,5 % 97,5 %
## factor(pOeros15)0 94,518 97,972
## factor(pOeros15)1 2,028 5,482
pOeros14 <- svymean(~factor(pOeros14), diseniopost1, na.rm = TRUE)
pOeros14 * 100
## mean SE
## factor(pOeros14)0 89,6 0,02
## factor(pOeros14)1 10,4 0,02
confint(pOeros14) * 100
## 2,5 % 97,5 %
## factor(pOeros14)0 86,643 92,56
## factor(pOeros14)1 7,437 13,36
pOeros13 <- svymean(~factor(pOeros13), diseniopost1, na.rm = TRUE)
## Error: contrasts can be applied only to factors with 2 or more levels
pOeros13 * 100
## Error: object 'pOeros13' not found
confint(pOeros13) * 100
## Error: object 'pOeros13' not found
pOeros12 <- svymean(~factor(pOeros12), diseniopost1, na.rm = TRUE)
## Error: contrasts can be applied only to factors with 2 or more levels
pOeros12 * 100
## Error: object 'pOeros12' not found
confint(pOeros12) * 100
## Error: object 'pOeros12' not found
pOeros11 <- svymean(~factor(pOeros11), diseniopost1, na.rm = TRUE)
## Error: contrasts can be applied only to factors with 2 or more levels
pOeros11 * 100
## Error: object 'pOeros11' not found
confint(pOeros1) * 100
## Error: object 'pOeros1' not found